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 | Float_t fEvtNumber;
|
---|
12 | UShort_t fPartId; // Type of particle
|
---|
13 | Float_t fEnergy; // [GeV] Energy
|
---|
14 | Float_t fThick0;
|
---|
15 | Float_t fFirstTarget;
|
---|
16 | Float_t fzFirstInt;
|
---|
17 |
|
---|
18 | Float_t fTheta; // [rad] Theta angle of event
|
---|
19 | Float_t fPhi; // [rad] Phi angle of event
|
---|
20 |
|
---|
21 | Float_t fCoreD; // Core d pos
|
---|
22 | Float_t fCoreX; // Core x pos
|
---|
23 | Float_t fCoreY; // Core y pos
|
---|
24 | Float_t fImpact; // [mm] 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;
|
---|
34 | Float_t fTimeLast;
|
---|
35 |
|
---|
36 | /* 6 parameters and chi2 of the NKG fit to the longitudinal
|
---|
37 | * particle distribution (see CORSIKA manual for explanation):
|
---|
38 | */
|
---|
39 | Float_t flongi_Nmax;
|
---|
40 | Float_t flongi_t0;
|
---|
41 | Float_t flongi_tmax;
|
---|
42 | Float_t flongi_a;
|
---|
43 | Float_t flongi_b;
|
---|
44 | Float_t flongi_c;
|
---|
45 | Float_t flongi_chi2;
|
---|
46 |
|
---|
47 | UInt_t fPhotIni; // Initial number of photons
|
---|
48 | UInt_t fPassPhotAtm; // Passed atmosphere
|
---|
49 | UInt_t fPassPhotRef; // Passed reflector(reflectivity + effective area)
|
---|
50 | UInt_t fPassPhotCone; // Passed glas
|
---|
51 | UInt_t fPhotElfromShower; // Passed qe coming from the shower
|
---|
52 | UInt_t fPhotElinCamera; // usPhotElfromShower + mean of phe
|
---|
53 | // from NSB
|
---|
54 |
|
---|
55 | public:
|
---|
56 | MMcEvt() ;
|
---|
57 |
|
---|
58 | MMcEvt( Float_t, UShort_t, Float_t, Float_t, Float_t,
|
---|
59 | Float_t, Float_t, Float_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 | UInt_t, UInt_t, UInt_t, UInt_t, UInt_t, UInt_t ) ;
|
---|
63 |
|
---|
64 | ~MMcEvt();
|
---|
65 |
|
---|
66 | void Clear(Option_t *opt=NULL);
|
---|
67 |
|
---|
68 | void Fill( Float_t, UShort_t, Float_t, Float_t, Float_t,
|
---|
69 | Float_t, Float_t, Float_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 | UInt_t, UInt_t, UInt_t, UInt_t, UInt_t, UInt_t ) ;
|
---|
73 |
|
---|
74 | //virtual void AsciiWrite(ofstream &fout) const;
|
---|
75 |
|
---|
76 | void Print(Option_t *opt=NULL) const;
|
---|
77 |
|
---|
78 | Short_t GetPartId() const { return fPartId; } //Get Type of particle
|
---|
79 | Float_t GetEnergy() const { return fEnergy; } //Get Energy
|
---|
80 |
|
---|
81 | Float_t GetTheta() const { return fTheta; } //Get Theta angle
|
---|
82 | Float_t GetPhi() const { return fPhi ; } //Get Phi angle
|
---|
83 |
|
---|
84 | /* Float_t GetCoreD() { return fCoreD; } //Get Core d pos */
|
---|
85 | /* Float_t GetCoreX() { return fCoreX; } //Get Core x pos */
|
---|
86 | /* Float_t GetCoreY() { return fCoreY; } //Get Core y pos */
|
---|
87 | Float_t GetImpact() const { return fImpact;} //Get impact parameter
|
---|
88 |
|
---|
89 | /* UInt_t GetPhotIni() { return fPhotIni; } //Get Initial photons */
|
---|
90 | /* UInt_t GetPassPhotAtm() { return fPassPhotAtm;} //Get Passed atmosphere */
|
---|
91 | /* UInt_t GetPassPhotRef() { return fPassPhotRef; } //Get Passed reflector */
|
---|
92 | /* UInt_t GetPassPhotCone() { return fPassPhotCone; } //Get Passed glas */
|
---|
93 | UInt_t GetPhotElfromShower() { return fPhotElfromShower; } //Get Passed qe from shower
|
---|
94 | /* UInt_t GetPhotElinCamera() { return fPhotElinCamera; } //Get Passed qe total*/
|
---|
95 |
|
---|
96 |
|
---|
97 | void SetPartId(Short_t PartId)
|
---|
98 | {fPartId=PartId;} //Set Type of particle
|
---|
99 |
|
---|
100 | void SetEnergy(Float_t Energy)
|
---|
101 | { fEnergy=Energy; } //Set Energy
|
---|
102 |
|
---|
103 | void SetTheta(Float_t Theta)
|
---|
104 | { fTheta=Theta; } //Set Theta angle
|
---|
105 |
|
---|
106 | void SetPhi(Float_t Phi)
|
---|
107 | { fPhi=Phi; } //Set Phi angle
|
---|
108 |
|
---|
109 | void SetCoreD(Float_t CoreD)
|
---|
110 | { fCoreD=CoreD; } //Set Core d pos
|
---|
111 |
|
---|
112 | void SetCoreX(Float_t CoreX)
|
---|
113 | { fCoreX=CoreX; } //Set Core x pos
|
---|
114 |
|
---|
115 | void SetCoreY(Float_t CoreY )
|
---|
116 | { fCoreY=CoreY; } //Set Core y pos
|
---|
117 |
|
---|
118 | void SetImpact(Float_t Impact)
|
---|
119 | { fImpact=Impact;} //Set impact parameter
|
---|
120 |
|
---|
121 | /* void SetPhotIni(Short_t PhotIni) */
|
---|
122 | /* { fPhotIni=PhotIni; } //Set Initial photons */
|
---|
123 | /* void SetPassPhotAtm(Short_t PassPhotAtm) */
|
---|
124 | /* { fPassPhotAtm=PassPhotAtm;} //Set Passed atmosphere */
|
---|
125 | /* void SetPassPhotRef(Short_t PassPhotRef) */
|
---|
126 | /* { fPassPhotRef=PassPhotRef ; } //Set Passed reflector */
|
---|
127 | /* void SetPassPhotCone(Short_t PhotCon) */
|
---|
128 | /* { fPassPhotCone=PhotCon; } //Set Passed glas */
|
---|
129 |
|
---|
130 |
|
---|
131 | ClassDef(MMcEvt, 2) //Stores Montecarlo Information of one event (eg. the energy)
|
---|
132 |
|
---|
133 | };
|
---|
134 |
|
---|
135 | #endif
|
---|