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

Last change on this file since 1809 was 1809, 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
112 void SetPartId(Short_t PartId)
113 {fPartId=PartId;} //Set Type of particle
114
115 void SetEnergy(Float_t Energy)
116 { fEnergy=Energy; } //Set Energy
117
118 void SetTheta(Float_t Theta)
119 { fTheta=Theta; } //Set Theta angle
120
121 void SetPhi(Float_t Phi)
122 { fPhi=Phi; } //Set Phi angle
123
124 void SetCoreD(Float_t CoreD)
125 { fCoreD=CoreD; } //Set Core d pos
126
127 void SetCoreX(Float_t CoreX)
128 { fCoreX=CoreX; } //Set Core x pos
129
130 void SetCoreY(Float_t CoreY )
131 { fCoreY=CoreY; } //Set Core y pos
132
133 void SetImpact(Float_t Impact)
134 { fImpact=Impact;} //Set impact parameter
135
136/* void SetPhotIni(Short_t PhotIni) */
137/* { fPhotIni=PhotIni; } //Set Initial photons */
138/* void SetPassPhotAtm(Short_t PassPhotAtm) */
139/* { fPassPhotAtm=PassPhotAtm;} //Set Passed atmosphere */
140/* void SetPassPhotRef(Short_t PassPhotRef) */
141/* { fPassPhotRef=PassPhotRef ; } //Set Passed reflector */
142/* void SetPassPhotCone(Short_t PhotCon) */
143/* { fPassPhotCone=PhotCon; } //Set Passed glas */
144
145
146 ClassDef(MMcEvt, 3) //Stores Montecarlo Information of one event (eg. the energy)
147
148};
149
150#endif
151
152
153
154
155
Note: See TracBrowser for help on using the repository browser.