MdShowerRegenerator.h

Go to the documentation of this file.

#ifndef _MdShowerRegeneratorAG_MdShowerRegenerator_h_
#define _MdShowerRegeneratorAG_MdShowerRegenerator_h_

#include <fwk/VModule.h>
#include <utl/Particle.h>
#include <utl/RandomEngine.h>
#include <utl/ShowerParticleIterator.h>
#include <utl/TimeDistribution.h>
#include <utl/AugerUnits.h>


namespace utl {
  class ShowerParticleIterator;
  class TimeInterval;
  class TabulatedFunction;
}

namespace mevt {
  class Counter;
}

namespace mdet {
  class Counter;
}


namespace MdShowerRegeneratorAG {

  struct InStats;

  class MdShowerRegenerator : public fwk::VModule {
   /* Regenerates particle for ground propagation 
    * the idea is similar to what is done for the tank 
    * but regeneration algorithm must consider that in this
    * case the detector volume is not the tank itself but the soil
    * beneath the tank till the muon counter (which is the active 
    * component of the detector. In some sense the equivalence is:
    * water  -> soil
    * pmts   -> scintillators
    */
  public:
    MdShowerRegenerator();

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

  private:

    bool IsParticleEnergyLow(const int, const double) const;
    bool IsMuon( const int );
    bool IsIn( double phi, double r, InStats & );
    void SetRadRange(double r1,  double r2);
    void SetPhiRange(double phi, double dphi);
    bool IsPhiIn(double phi);
    bool IsRIn(double r);

    // Parameters taken from xml file
    double fInnerRadiusCut;
    double fOuterRadiusCut;

    double fElectronEnergyCut;
    double fMuonEnergyCut;
    double fPhotonEnergyCut;
    double fHadronEnergyCut;
    double fMesonEnergyCut;

    double fDeltaROverR;
    double fDeltaPhi;
    double fHorizontalParticleCut;

    double fPhi1;
    double fPhi2;
    double fR1;
    double fR2;

    utl::RandomEngine::RandomEngineType* fRandomEngine;

    bool fOnlyMuons;
    int fNMuons;
    int fNInjected; 
    int fNHorizontalRejected;
    double fMuonWeightScale;

    REGISTER_MODULE("ShowerRegeneratorAG", MdShowerRegenerator);
  };

} 

#endif

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