source: trunk/MagicSoft/include-Classes/MMcFormat/MMcEvt.hxx@ 1950

Last change on this file since 1950 was 1868, checked in by wittek, 22 years ago
*** empty log message ***
File size: 5.4 KB
Line 
1#ifndef __MMcEvt__
2#define __MMcEvt__
3
4#ifndef MARS_MParContainer
5#include "MParContainer.h"
6#endif
7
8class 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 // Now follow the fraction of photons reaching the camera produced by
57 // electrons, muons and other particles respectively:
58
59 Float_t fElecCphFraction;
60 Float_t fMuonCphFraction;
61 Float_t fOtherCphFraction;
62
63
64
65 public:
66 MMcEvt() ;
67
68 MMcEvt( UInt_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 Float_t, Float_t, Float_t) ;
74
75 ~MMcEvt();
76
77 void Clear(Option_t *opt=NULL);
78
79 void Fill( UInt_t, UShort_t, Float_t, Float_t, Float_t,
80 Float_t, Float_t, Float_t, Float_t, Float_t, Float_t,
81 Float_t, Float_t, Float_t, Float_t, Float_t, Float_t,
82 Float_t, Float_t, Float_t, Float_t, Float_t, Float_t,
83 UInt_t, UInt_t, UInt_t, UInt_t, UInt_t, UInt_t,
84 Float_t, Float_t, Float_t) ;
85
86 //virtual void AsciiWrite(ofstream &fout) const;
87
88 void Print(Option_t *opt=NULL) const;
89
90 UInt_t GetEvtNumber() const { return fEvtNumber; } //Get Event Number
91 Short_t GetPartId() const { return fPartId; } //Get Type of particle
92 Float_t GetEnergy() const { return fEnergy; } //Get Energy
93
94 Float_t GetTheta() const { return fTheta; } //Get Theta angle
95 Float_t GetPhi() const { return fPhi ; } //Get Phi angle
96
97/* Float_t GetCoreD() { return fCoreD; } //Get Core d pos */
98/* Float_t GetCoreX() { return fCoreX; } //Get Core x pos */
99/* Float_t GetCoreY() { return fCoreY; } //Get Core y pos */
100 Float_t GetImpact() const { return fImpact;} //Get impact parameter
101
102/* UInt_t GetPhotIni() { return fPhotIni; } //Get Initial photons */
103/* UInt_t GetPassPhotAtm() { return fPassPhotAtm;} //Get Passed atmosphere */
104/* UInt_t GetPassPhotRef() { return fPassPhotRef; } //Get Passed reflector */
105/* UInt_t GetPassPhotCone() { return fPassPhotCone; } //Get Passed glas */
106 UInt_t GetPhotElfromShower() { return fPhotElfromShower; } //Get Passed qe from shower
107/* UInt_t GetPhotElinCamera() { return fPhotElinCamera; } //Get Passed qe total*/
108
109 Float_t GetTelescopePhi() const { return fTelescopePhi; }
110 Float_t GetTelescopeTheta() const { return fTelescopeTheta; }
111 Float_t GetOtherCphFraction() const { return fOtherCphFraction; }
112
113 void SetPartId(Short_t PartId)
114 {fPartId=PartId;} //Set Type of particle
115
116 void SetEnergy(Float_t Energy)
117 { fEnergy=Energy; } //Set Energy
118
119 void SetTheta(Float_t Theta)
120 { fTheta=Theta; } //Set Theta angle
121
122 void SetPhi(Float_t Phi)
123 { fPhi=Phi; } //Set Phi angle
124
125 void SetCoreD(Float_t CoreD)
126 { fCoreD=CoreD; } //Set Core d pos
127
128 void SetCoreX(Float_t CoreX)
129 { fCoreX=CoreX; } //Set Core x pos
130
131 void SetCoreY(Float_t CoreY )
132 { fCoreY=CoreY; } //Set Core y pos
133
134 void SetImpact(Float_t Impact)
135 { fImpact=Impact;} //Set impact parameter
136
137/* void SetPhotIni(Short_t PhotIni) */
138/* { fPhotIni=PhotIni; } //Set Initial photons */
139/* void SetPassPhotAtm(Short_t PassPhotAtm) */
140/* { fPassPhotAtm=PassPhotAtm;} //Set Passed atmosphere */
141/* void SetPassPhotRef(Short_t PassPhotRef) */
142/* { fPassPhotRef=PassPhotRef ; } //Set Passed reflector */
143/* void SetPassPhotCone(Short_t PhotCon) */
144/* { fPassPhotCone=PhotCon; } //Set Passed glas */
145
146
147 ClassDef(MMcEvt, 3) //Stores Montecarlo Information of one event (eg. the energy)
148
149};
150
151#endif
152
153
154
155
156
Note: See TracBrowser for help on using the repository browser.