| 1 | #ifndef __MMcEvt__
|
|---|
| 2 | #define __MMcEvt__
|
|---|
| 3 |
|
|---|
| 4 | #ifndef MARS_MParContainer
|
|---|
| 5 | #include "MParContainer.h"
|
|---|
| 6 | #endif
|
|---|
| 7 |
|
|---|
| 8 | class MMcEvt : public MParContainer
|
|---|
| 9 | {
|
|---|
| 10 | private:
|
|---|
| 11 | UInt_t fEvtNumber;
|
|---|
| 12 | UShort_t fPartId; // Type of particle
|
|---|
| 13 | Float_t fEnergy; // [GeV] Energy
|
|---|
| 14 | Float_t fThick0; // [g/cm2]
|
|---|
| 15 | Float_t fFirstTarget; // []
|
|---|
| 16 | Float_t fZFirstInteraction; // [cm]
|
|---|
| 17 |
|
|---|
| 18 | Float_t fTheta; // [rad] Theta angle of event
|
|---|
| 19 | Float_t fPhi; // [rad] Phi angle of event
|
|---|
| 20 |
|
|---|
| 21 | Float_t fCoreD; // [cm] Core d pos
|
|---|
| 22 | Float_t fCoreX; // [cm] Core x pos
|
|---|
| 23 | Float_t fCoreY; // [cm] Core y pos
|
|---|
| 24 | Float_t fImpact; // [cm] impact parameter
|
|---|
| 25 |
|
|---|
| 26 | // Up to here, the info from the CORSIKA event header.
|
|---|
| 27 |
|
|---|
| 28 | // Telescope orientation:
|
|---|
| 29 | Float_t fTelescopePhi; // [rad]
|
|---|
| 30 | Float_t fTelescopeTheta; // [rad]
|
|---|
| 31 |
|
|---|
| 32 | // Time of first and last photon:
|
|---|
| 33 | Float_t fTimeFirst; // [ns]
|
|---|
| 34 | Float_t fTimeLast; // [ns]
|
|---|
| 35 |
|
|---|
| 36 | // 6 parameters and chi2 of the NKG fit to the longitudinal
|
|---|
| 37 | // particle distribution. See CORSIKA manual for explanation,
|
|---|
| 38 | // section 4.42 "Longitudinal shower development":
|
|---|
| 39 | //
|
|---|
| 40 | Float_t fLongiNmax; // [particles]
|
|---|
| 41 | Float_t fLongit0; // [g/cm2]
|
|---|
| 42 | Float_t fLongitmax; // [g/cm2]
|
|---|
| 43 | Float_t fLongia; // [g/cm2]
|
|---|
| 44 | Float_t fLongib; // []
|
|---|
| 45 | Float_t fLongic; // [cm2/g]
|
|---|
| 46 | Float_t fLongichi2;
|
|---|
| 47 |
|
|---|
| 48 | UInt_t fPhotIni; // [ph] Initial number of photons
|
|---|
| 49 | UInt_t fPassPhotAtm; // [ph] Passed atmosphere
|
|---|
| 50 | UInt_t fPassPhotRef; // [ph] Passed reflector(reflectivity + effective area)
|
|---|
| 51 | UInt_t fPassPhotCone; // [ph] Passed glas
|
|---|
| 52 | UInt_t fPhotElfromShower; // [phe] Passed qe coming from the shower
|
|---|
| 53 | UInt_t fPhotElinCamera; // [phe] usPhotElfromShower + mean of phe
|
|---|
| 54 | // from NSB
|
|---|
| 55 |
|
|---|
| 56 | public:
|
|---|
| 57 | MMcEvt() ;
|
|---|
| 58 |
|
|---|
| 59 | MMcEvt( UInt_t, UShort_t, Float_t, Float_t, Float_t,
|
|---|
| 60 | Float_t, Float_t, Float_t, Float_t, Float_t, Float_t,
|
|---|
| 61 | Float_t, Float_t, Float_t, Float_t, Float_t, Float_t,
|
|---|
| 62 | Float_t, Float_t, Float_t, Float_t, Float_t, Float_t,
|
|---|
| 63 | UInt_t, UInt_t, UInt_t, UInt_t, UInt_t, UInt_t ) ;
|
|---|
| 64 |
|
|---|
| 65 | ~MMcEvt();
|
|---|
| 66 |
|
|---|
| 67 | void Clear(Option_t *opt=NULL);
|
|---|
| 68 |
|
|---|
| 69 | void Fill( UInt_t, UShort_t, Float_t, Float_t, Float_t,
|
|---|
| 70 | Float_t, Float_t, Float_t, Float_t, Float_t, Float_t,
|
|---|
| 71 | Float_t, Float_t, Float_t, Float_t, Float_t, Float_t,
|
|---|
| 72 | Float_t, Float_t, Float_t, Float_t, Float_t, Float_t,
|
|---|
| 73 | UInt_t, UInt_t, UInt_t, UInt_t, UInt_t, UInt_t ) ;
|
|---|
| 74 |
|
|---|
| 75 | //virtual void AsciiWrite(ofstream &fout) const;
|
|---|
| 76 |
|
|---|
| 77 | void Print(Option_t *opt=NULL) const;
|
|---|
| 78 |
|
|---|
| 79 | Short_t GetPartId() const { return fPartId; } //Get Type of particle
|
|---|
| 80 | Float_t GetEnergy() const { return fEnergy; } //Get Energy
|
|---|
| 81 |
|
|---|
| 82 | Float_t GetTheta() const { return fTheta; } //Get Theta angle
|
|---|
| 83 | Float_t GetPhi() const { return fPhi ; } //Get Phi angle
|
|---|
| 84 |
|
|---|
| 85 | /* Float_t GetCoreD() { return fCoreD; } //Get Core d pos */
|
|---|
| 86 | /* Float_t GetCoreX() { return fCoreX; } //Get Core x pos */
|
|---|
| 87 | /* Float_t GetCoreY() { return fCoreY; } //Get Core y pos */
|
|---|
| 88 | Float_t GetImpact() const { return fImpact;} //Get impact parameter
|
|---|
| 89 |
|
|---|
| 90 | /* UInt_t GetPhotIni() { return fPhotIni; } //Get Initial photons */
|
|---|
| 91 | /* UInt_t GetPassPhotAtm() { return fPassPhotAtm;} //Get Passed atmosphere */
|
|---|
| 92 | /* UInt_t GetPassPhotRef() { return fPassPhotRef; } //Get Passed reflector */
|
|---|
| 93 | /* UInt_t GetPassPhotCone() { return fPassPhotCone; } //Get Passed glas */
|
|---|
| 94 | UInt_t GetPhotElfromShower() { return fPhotElfromShower; } //Get Passed qe from shower
|
|---|
| 95 | /* UInt_t GetPhotElinCamera() { return fPhotElinCamera; } //Get Passed qe total*/
|
|---|
| 96 |
|
|---|
| 97 |
|
|---|
| 98 | void SetPartId(Short_t PartId)
|
|---|
| 99 | {fPartId=PartId;} //Set Type of particle
|
|---|
| 100 |
|
|---|
| 101 | void SetEnergy(Float_t Energy)
|
|---|
| 102 | { fEnergy=Energy; } //Set Energy
|
|---|
| 103 |
|
|---|
| 104 | void SetTheta(Float_t Theta)
|
|---|
| 105 | { fTheta=Theta; } //Set Theta angle
|
|---|
| 106 |
|
|---|
| 107 | void SetPhi(Float_t Phi)
|
|---|
| 108 | { fPhi=Phi; } //Set Phi angle
|
|---|
| 109 |
|
|---|
| 110 | void SetCoreD(Float_t CoreD)
|
|---|
| 111 | { fCoreD=CoreD; } //Set Core d pos
|
|---|
| 112 |
|
|---|
| 113 | void SetCoreX(Float_t CoreX)
|
|---|
| 114 | { fCoreX=CoreX; } //Set Core x pos
|
|---|
| 115 |
|
|---|
| 116 | void SetCoreY(Float_t CoreY )
|
|---|
| 117 | { fCoreY=CoreY; } //Set Core y pos
|
|---|
| 118 |
|
|---|
| 119 | void SetImpact(Float_t Impact)
|
|---|
| 120 | { fImpact=Impact;} //Set impact parameter
|
|---|
| 121 |
|
|---|
| 122 | /* void SetPhotIni(Short_t PhotIni) */
|
|---|
| 123 | /* { fPhotIni=PhotIni; } //Set Initial photons */
|
|---|
| 124 | /* void SetPassPhotAtm(Short_t PassPhotAtm) */
|
|---|
| 125 | /* { fPassPhotAtm=PassPhotAtm;} //Set Passed atmosphere */
|
|---|
| 126 | /* void SetPassPhotRef(Short_t PassPhotRef) */
|
|---|
| 127 | /* { fPassPhotRef=PassPhotRef ; } //Set Passed reflector */
|
|---|
| 128 | /* void SetPassPhotCone(Short_t PhotCon) */
|
|---|
| 129 | /* { fPassPhotCone=PhotCon; } //Set Passed glas */
|
|---|
| 130 |
|
|---|
| 131 |
|
|---|
| 132 | ClassDef(MMcEvt, 2) //Stores Montecarlo Information of one event (eg. the energy)
|
|---|
| 133 |
|
|---|
| 134 | };
|
|---|
| 135 |
|
|---|
| 136 | #endif
|
|---|