| 1 | #ifndef MARS_MHMcRate | 
|---|
| 2 | #define MARS_MHMcRate | 
|---|
| 3 |  | 
|---|
| 4 | #ifndef MARS_MParContainer | 
|---|
| 5 | #include "MParContainer.h" | 
|---|
| 6 | #endif | 
|---|
| 7 |  | 
|---|
| 8 | class MHMcRate : public MParContainer | 
|---|
| 9 | { | 
|---|
| 10 |  | 
|---|
| 11 | private: | 
|---|
| 12 | UShort_t fPartId;           // Type of particle | 
|---|
| 13 |  | 
|---|
| 14 | Float_t fEnergyMax;         // Maximum Energy [TeV] | 
|---|
| 15 | Float_t fEnergyMin;         // Minimum Energy [TeV] | 
|---|
| 16 |  | 
|---|
| 17 | Float_t fThetaMax;          // Maximum theta angle of run [?] | 
|---|
| 18 | Float_t fThetaMin;          // Minimum theta angle of run [?] | 
|---|
| 19 | Float_t fPhiMax;            // Maximum phi angle of run [?] | 
|---|
| 20 | Float_t fPhiMin;            // Minimum phi angle of run [?] | 
|---|
| 21 |  | 
|---|
| 22 | Float_t fSolidAngle;        // Solid angle within which incident directions are distributed [sr] | 
|---|
| 23 |  | 
|---|
| 24 | Float_t fImpactMax;         //[cm] Maximum impact parameter [cm] | 
|---|
| 25 | Float_t fImpactMin;         //[cm] Minimum impact parameter [cm] | 
|---|
| 26 |  | 
|---|
| 27 | Float_t fBackTrig;          // Number of triggers from background | 
|---|
| 28 | Float_t fBackSim;           // Number of simulated showers for the background | 
|---|
| 29 |  | 
|---|
| 30 | Float_t fSpecIndex;         // dn/dE = k * e^{- fSpecIndex} | 
|---|
| 31 | Float_t fFlux0;             // dn/dE = fFlux0 * E^{-a} | 
|---|
| 32 |  | 
|---|
| 33 | Float_t fShowerRate;        // Showers rate in Hz | 
|---|
| 34 | Float_t fShowerRateError;   // Estimated error of shower rate in Hz | 
|---|
| 35 |  | 
|---|
| 36 | Float_t fTriggerRate;       // Trigger rate in Hz | 
|---|
| 37 | Float_t fTriggerRateError;  // Estimated error for the trigger rate in Hz | 
|---|
| 38 |  | 
|---|
| 39 | Float_t fMeanThreshold;     // Mean discriminator threshold of trigger pixels [mV] | 
|---|
| 40 |  | 
|---|
| 41 | Short_t fMultiplicity;      // L1 trigger multiplicity. | 
|---|
| 42 |  | 
|---|
| 43 | Short_t fTriggerCondNum;    // Trigger condition number, for the case of running over camra files containing several. | 
|---|
| 44 |  | 
|---|
| 45 | void Init(const char *name, const char *title); | 
|---|
| 46 |  | 
|---|
| 47 | public: | 
|---|
| 48 |  | 
|---|
| 49 | MHMcRate(const char *name=NULL, const char *title=NULL); | 
|---|
| 50 | MHMcRate(Float_t showrate, | 
|---|
| 51 | const char *name=NULL, const char *title=NULL); | 
|---|
| 52 | MHMcRate(Float_t specindex, Float_t flux0, | 
|---|
| 53 | const char *name=NULL, const char *title=NULL); | 
|---|
| 54 |  | 
|---|
| 55 | void SetParticle(UShort_t part); | 
|---|
| 56 | void SetBackground(Float_t showers, Float_t triggers); | 
|---|
| 57 | void SetFlux(Float_t flux0, Float_t specindx); | 
|---|
| 58 | void SetIncidentRate(Float_t showerrate); | 
|---|
| 59 |  | 
|---|
| 60 | void SetImpactMax(Float_t Impact) {fImpactMax=Impact;} | 
|---|
| 61 | void SetImpactMin(Float_t Impact) {fImpactMin=Impact;} | 
|---|
| 62 |  | 
|---|
| 63 | void SetThetaMax(Float_t Theta) {fThetaMax=Theta;} | 
|---|
| 64 | void SetThetaMin(Float_t Theta) {fThetaMin=Theta;} | 
|---|
| 65 | void SetPhiMax(Float_t Phi) {fPhiMax=Phi;} | 
|---|
| 66 | void SetPhiMin(Float_t Phi) {fPhiMin=Phi;} | 
|---|
| 67 |  | 
|---|
| 68 | void SetSolidAngle(Float_t Solid) {fSolidAngle=Solid;} | 
|---|
| 69 | void SetEnergyMax(Float_t Energy) {fEnergyMax=Energy;} | 
|---|
| 70 | void SetEnergyMin(Float_t Energy) {fEnergyMin=Energy;} | 
|---|
| 71 |  | 
|---|
| 72 | void SetMultiplicity(Short_t nMul) {fMultiplicity = nMul;} | 
|---|
| 73 | void SetMeanThreshold(Float_t thresh) {fMeanThreshold = thresh;} | 
|---|
| 74 |  | 
|---|
| 75 | void SetTriggerCondNum(Short_t num) {fTriggerCondNum = num;} | 
|---|
| 76 |  | 
|---|
| 77 | void UpdateBoundaries(Float_t energy, Float_t theta, Float_t phi, Float_t impact); | 
|---|
| 78 |  | 
|---|
| 79 | Float_t GetTriggerRate() {return fTriggerRate;} | 
|---|
| 80 | Float_t GetTriggerRateError() {return fTriggerRateError;} | 
|---|
| 81 |  | 
|---|
| 82 | Short_t GetMultiplicity() {return fMultiplicity;} | 
|---|
| 83 | Float_t GetMeanThreshold() {return fMeanThreshold;} | 
|---|
| 84 | Short_t GetTriggerCondNum() {return fTriggerCondNum;} | 
|---|
| 85 |  | 
|---|
| 86 |  | 
|---|
| 87 | void CalcRate(Float_t trig, Float_t anal, Float_t simu); | 
|---|
| 88 |  | 
|---|
| 89 | void Print(Option_t *o=NULL) const; | 
|---|
| 90 |  | 
|---|
| 91 | void Draw(Option_t *o=NULL); | 
|---|
| 92 |  | 
|---|
| 93 | ClassDef(MHMcRate, 1)  // Data Container to calculate trigger rate | 
|---|
| 94 | }; | 
|---|
| 95 |  | 
|---|
| 96 | #endif | 
|---|
| 97 |  | 
|---|
| 98 |  | 
|---|
| 99 |  | 
|---|