FdSimEventChecker.h

Go to the documentation of this file.

#ifndef _FdSimEventChecker_h
#define _FdSimEventChecker_h

#include <fwk/VModule.h>
#include <TH2.h>

namespace utl {
  class RandomEngine;
}

namespace fdet{
  class Telescope;
  class Eye;
}

namespace evt {
  class ShowerSimData;
}

namespace FdSimEventCheckerOG {

  //#warning Rename module, to something like FdSimEventBuilder  

  class FdSimEventChecker : public fwk::VModule{
    
  public:
    FdSimEventChecker() { }
   ~FdSimEventChecker() { }
    
    fwk::VModule::ResultFlag Init();
    fwk::VModule::ResultFlag Run(evt::Event& event);
    fwk::VModule::ResultFlag Finish();

  private :
    double Rcutoff(const double lgE) const;
    bool PCGFTriggerProbabilityCut(const double lgE, const double dist, const double angle) const;
    bool PCGFTriggerProbabilityCutHECO(const double lgE, const double dist, const double angle) const;
    void DumpDB() const;

    double InterpolatePCGFAcceptanceMap(const double lgE, const double dist) const;
    
    double CalculateMinViewingAngle(const fdet::Telescope& tel, const evt::ShowerSimData& shower) const;
    double CalculateXmaxViewingAngle(const fdet::Eye& eye, const evt::ShowerSimData& shower) const;   
      
    // options
    int fVerbosity;
    
    bool fUpTimeCheck;
    bool fGeometryCheck; // for normal Fd events, results from Sergio
    bool fGeometryCheckCherenkoHEAT; // for HEAT Cherenkov, PCGF
    bool fGeometryCheckCherenkoHECO; // for HECO Cherenkov, PCGF
    TH2D* hMaxVA;
    bool fFastSkip;
    bool fFastRepeat;
    bool fUseCDAS;

    double fMinViewingAngle;
    
    utl::RandomEngine* fRandomEngine;
    
    REGISTER_MODULE("FdSimEventCheckerOG", FdSimEventChecker);
  };

}

#endif