UnivParamNSTester.cc

Go to the documentation of this file.

#include "UnivParamNSTester.h"

#include <utl/ErrorLogger.h>

#include <tls/UnivParam.h>

#include <iostream>

using namespace fwk;
using namespace utl;

using namespace std;


namespace UnivParamNSTester {


  VModule::ResultFlag
  UnivParamNSTester::Init()
  {
    return eSuccess;
  }


  VModule::ResultFlag
  UnivParamNSTester::Run(evt::Event& /*event*/)
  {
    UnivParamNS::UnivParam fUnivParam(0);  // WCD(0) (see UnivParamNS/UnivParam.h)

    double logTotalEnergy = 19.;
    double r = 1000 * meter;
    double Xmax_EM = 800 * gram / cm2;
    double Xmax_MU = 500 * gram / cm2;
    double theta = 35 * degree;
    double psi = 90 * degree;
    double Nmu = 1.0;
    int monthlyModel = 4; // Malargue April
    double hGround = 1452 * meter;
    double rhoGround = 1.04e-3 * gram / cm3;

    int icomp;
    double componentSignal;

    cout << endl;

    cout << "log10 E/eV = " << logTotalEnergy << endl;
    cout << "r/m = " << r / meter << endl;
    cout << "Xmax_EM/gcm^-2 = " << Xmax_EM / (gram / cm2) << endl;
    cout << "Xmax_MU/gcm^-2 = " << Xmax_MU / (gram / cm2) << endl;
    cout << "theta/deg = " << theta / degree << endl;
    cout << "psi/deg = " << psi / degree << endl;
    cout << "height/m = " << hGround / meter << endl;
    cout << "density/gcm^-3) = " << rhoGround / (gram / cm3) << endl;
    cout << "Nmu = " << Nmu << endl;

    cout << endl;

    cout << "component signals / VEM" << endl;

    icomp = 0; // pure muonic
    componentSignal = fUnivParam.GetSignal(r / cm, Xmax_EM / (gram / cm2), logTotalEnergy,
                                           theta, psi, rhoGround / (gram / cm3), hGround / cm,
                                           Nmu, icomp, monthlyModel);
    cout << "pure muonic " << componentSignal << endl;

    icomp = 1; // pure electromagnetic
    componentSignal = fUnivParam.GetSignal(r / cm, Xmax_EM / (gram / cm2), logTotalEnergy,
                                           theta, psi, rhoGround / (gram / cm3), hGround / cm,
                                           Nmu, icomp, monthlyModel);
    cout << "pure electromagnetic " << componentSignal << endl;

    icomp = 2; // muon halo
    componentSignal = fUnivParam.GetSignal(r / cm, Xmax_EM / (gram / cm2), logTotalEnergy,
                                           theta, psi, rhoGround / (gram / cm3), hGround / cm,
                                           Nmu, icomp, monthlyModel);
    cout << "muon halo " << componentSignal << endl;

    icomp = 3; // hadron jets
    componentSignal = fUnivParam.GetSignal(r / cm, Xmax_EM / (gram / cm2), logTotalEnergy,
                                           theta, psi, rhoGround / (gram / cm3), hGround / cm,
                                           Nmu, icomp, monthlyModel);
    cout << "hadron jets " << componentSignal << endl;

    cout << endl;

    return eSuccess;
  }


  VModule::ResultFlag
  UnivParamNSTester::Finish()
  {
    return eSuccess;
  }

}