SdSimpleSim.h

Go to the documentation of this file.

#ifndef _SdSimpleSimKG_SdSimpleSim_h_
#define _SdSimpleSimKG_SdSimpleSim_h_


namespace utl {
  class RandomEngine;
  class TimeStamp;
}



#include <fwk/VModule.h>
#include <utl/Vector.h>
#include <utl/TabulatedFunction.h>
#include <evt/ShowerSimData.h>

#include <map>

class TTree;
class TProfile;

namespace sevt {
  class SEvent;
  class Station;
  class StationRecData;
}

namespace atm {
  class ProfileResult;
}

namespace SdSimpleSimKG {

  class SdSimpleSim : public fwk::VModule {

    struct FitParam {
      FitParam() : p0(0), p1(0) { }
      FitParam(double a, double b) : p0(a), p1(b) { }
      double p0;
      double p1;
    };

  public:
    SdSimpleSim();

    fwk::VModule::ResultFlag Init();
    fwk::VModule::ResultFlag Run(evt::Event& event);
    fwk::VModule::ResultFlag Finish();

  private:
    utl::TabulatedFunction GetLogZdist(const utl::TabulatedFunction& prof,
                                       const utl::Point& core,
                                       const utl::Vector& axis);

    utl::TabulatedFunction CalculateLogZdist (const utl::TabulatedFunction& MuonProfile,
                                              //const utl::Point& core,
                                              //const utl::Vector& axis);
                                              double cosTheta,
                                              double XobsVert);

    double SampleEnergy (double Emin,
                         double Emax,
                         double Age);

    bool Noise (sevt::SEvent& sevent,
                const utl::TimeStamp& T0,
                int TankID,
                double TimeOffset, double TimeWindow);

    bool IsNoise (double TimeWindow);

    double GenerateNoiseStation (sevt::Station& tank,
                                 const utl::TimeStamp& T0,
                                 double TimeOffset,
                                 double TimeWindow);

    double NKG (double N, double Rm, double R, double s);

    double CalculateTankSignal (double CoreDistance,
                                double TankRadius, double TankHeight, double Zenith,
                                double Age, double Rm,
                                double Ne, double Ng, double Nmu,
                                bool fluctuations=true);

    double TankIntersection (double r,
                             double phi,
                             double z,
                             double theta,
                             double TankRadius);

    double LTP (double r,
                double theta,
                double energy);

    double T1TriggerProbability (double signal, double S1000, double theta);

    double GetCDASTriggerTimeWindow(const utl::TimeStamp& time);

    double NerlingF1 (double E,
                      double a1);

    double NerlingF2 (double E,
                      double a2,
                      double s);

    // parameters for nerling model
    static const double fNerlingA1a;
    static const double fNerlingA1b;
    static const double fNerlingA2a;
    static const double fNerlingA2b;

    FitParam FitDecay (const utl::TabulatedFunction& prof, bool dEdX=false);
    FitParam FitAtmosphere (const atm::ProfileResult& prof);
    std::map <evt::ShowerSimData::ProfileType, FitParam> fProfileExtrapolation;
    FitParam fSlantDepthExtrapolation;

  private:
    utl::RandomEngine& fRandomEngine;

    int fVerbosity;
    bool fUseT2LifeTime;

    int fSdId;

    double fMuonRScale;
    double fEMAgeFactor;
    double fLateralCut;

    bool fNoise;
    bool fMuonProductionHeightFromProfile;

    int fNSample;

    // number of particles
    int fNMuons;
    int fNGammas;
    int fNElectrons;

    // VEM of particles
    double fSMuons;
    double fSGammas;
    double fSElectrons;

    double fNoiseRatePerStation;

    bool fForceClosestTankToTrigger;
    bool fOnlyClosestTank;
    bool fPlaneFront;

    // for debug root tree
    TTree *fStatTree;
    TProfile *fLorenzo_vs_X;
    TProfile *fdNdX;
    TProfile *fCharged;
    TProfile *fMuon;
    TProfile *fElectron;
    TProfile *fGamma;
    TProfile *fdEdX;

    TProfile *fLorenzo;
    TProfile *fdNdz;
    TProfile *fCharged_vs_z;
    TProfile *fMuon_vs_z;
    TProfile *fElectron_vs_z;
    TProfile *fGamma_vs_z;
    TProfile *fdEdX_vs_z;

    REGISTER_MODULE ("SdSimpleSimKG", SdSimpleSim);
  };


} // SdSimpleSimKG


#endif // _SdSimpleSimKG_SdSimpleSim_h_

// Configure (x)emacs for this file ...
// Local Variables:
// mode:c++
// compile-command: "make -C ..
// End: