source: tags/Mars-V2.0/mmc/MMcEvt.hxx

Last change on this file was 8438, checked in by tbretz, 17 years ago
*** empty log message ***
File size: 5.2 KB
Line 
1#ifndef MARS_MMcEvt
2#define MARS_MMcEvt
3
4#ifndef MARS_MMcEvtBasic
5#include "MMcEvtBasic.h"
6#endif
7
8
9class MMcEvt : public MMcEvtBasic
10{
11private:
12 UInt_t fEvtNumber;
13 Float_t fThick0; // [g/cm2]
14 Float_t fFirstTarget; // []
15 Float_t fZFirstInteraction; // [cm]
16
17 Float_t fTheta; // [rad] Theta angle of event
18 Float_t fPhi; // [rad] Phi angle of event (see class description)
19
20 Float_t fCoreD; // [cm] Core d pos
21 Float_t fCoreX; // [cm] Core x pos
22 Float_t fCoreY; // [cm] Core y pos
23
24 // Up to here, the info from the CORSIKA event header.
25
26 // Time of first and last photon:
27 Float_t fTimeFirst; // [ns]
28 Float_t fTimeLast; // [ns]
29
30 // 6 parameters and chi2 of the NKG fit to the longitudinal
31 // particle distribution. See CORSIKA manual for explanation,
32 // section 4.42 "Longitudinal shower development":
33 //
34 Float_t fLongiNmax; // [particles]
35 Float_t fLongit0; // [g/cm2]
36 Float_t fLongitmax; // [g/cm2]
37 Float_t fLongia; // [g/cm2]
38 Float_t fLongib; // []
39 Float_t fLongic; // [cm2/g]
40 Float_t fLongichi2;
41
42 UInt_t fPhotIni; // [ph] Initial number of photons
43 UInt_t fPassPhotAtm; // [ph] Passed atmosphere
44 UInt_t fPassPhotRef; // [ph] Passed reflector(reflectivity + effective area)
45 UInt_t fPassPhotCone; // [ph] Within any valid pixel, before plexiglas
46 UInt_t fPhotElfromShower; // [phe] Passed qe, coming from the shower
47 UInt_t fPhotElinCamera; // [phe] usPhotElfromShower + mean of phe from NSB
48
49 // Now follow the fraction of photons reaching the camera produced by
50 // electrons, muons and other particles respectively:
51
52 Float_t fElecCphFraction;
53 Float_t fMuonCphFraction;
54 Float_t fOtherCphFraction;
55
56 Float_t fFadcTimeJitter;
57
58 Int_t fEventReuse; // Number running from 0 to N-1, being N the number
59 // of times a Corsika event has been reused, by
60 // orienting the telescope in different ways or by
61 // setting it at a different location on the ground.
62
63public:
64 MMcEvt();
65 MMcEvt(UInt_t, ParticleId_t, Float_t, Float_t, Float_t,
66 Float_t, Float_t, Float_t, Float_t, Float_t, Float_t,
67 Float_t, Float_t, Float_t, Float_t, Float_t, Float_t,
68 Float_t, Float_t, Float_t, Float_t, Float_t, Float_t,
69 UInt_t, UInt_t, UInt_t, UInt_t, UInt_t, UInt_t,
70 Float_t, Float_t, Float_t, Float_t, Int_t ireuse=0) ;
71
72 // Getter
73 UInt_t GetEvtNumber() const { return fEvtNumber; } //Get Event Number
74 Float_t GetTheta() const { return fTheta; } //Get Theta angle
75 Float_t GetPhi() const { return fPhi ; } //Get Phi angle
76
77 Float_t GetCoreX() const { return fCoreX; } //Get Core x pos
78 Float_t GetCoreY() const { return fCoreY; } //Get Core y pos
79
80 UInt_t GetPhotIni() const { return fPhotIni; } //Get Initial photons
81 UInt_t GetPassPhotAtm() const { return fPassPhotAtm;} //Get Passed atmosphere
82 UInt_t GetPassPhotRef() const { return fPassPhotRef; } //Get Passed reflector
83 UInt_t GetPassPhotCone() const { return fPassPhotCone; } //Get Passed glas
84 UInt_t GetPhotElfromShower() const { return fPhotElfromShower; } //Get Passed qe from shower
85 UInt_t GetPhotElinCamera() const { return fPhotElinCamera; } //Get Passed qe total
86 Float_t GetZFirstInteraction() const { return fZFirstInteraction; }
87
88 Float_t GetOtherCphFraction() const { return fOtherCphFraction; }
89
90 Float_t GetLongiNmax() const { return fLongiNmax; }
91 Float_t GetLongia() const { return fLongia; }
92 Float_t GetLongib() const { return fLongib; }
93 Float_t GetLongic() const { return fLongic; }
94 Float_t GetLongichi2() const { return fLongichi2; }
95 Float_t GetLongit0() const { return fLongit0; }
96 Float_t GetLongitmax() const { return fLongitmax; }
97
98 Float_t GetFadcTimeJitter() const { return fFadcTimeJitter; }
99
100 Float_t GetMuonCphFraction() const { return fMuonCphFraction; }
101
102 TString GetDescription(const TString &s="") const;
103
104 // Setter
105 void SetTheta(Float_t Theta) { fTheta=Theta; } //Set Theta angle
106 void SetPhi(Float_t Phi) { fPhi=Phi; } //Set Phi angle
107 void SetCoreD(Float_t CoreD) { fCoreD=CoreD; } //Set Core d pos
108 void SetCoreX(Float_t CoreX) { fCoreX=CoreX; } //Set Core x pos
109 void SetCoreY(Float_t CoreY) { fCoreY=CoreY; } //Set Core y pos
110
111 void Fill( UInt_t, ParticleId_t, Float_t, Float_t, Float_t,
112 Float_t, Float_t, Float_t, Float_t, Float_t, Float_t,
113 Float_t, Float_t, Float_t, Float_t, Float_t, Float_t,
114 Float_t, Float_t, Float_t, Float_t, Float_t, Float_t,
115 UInt_t, UInt_t, UInt_t, UInt_t, UInt_t, UInt_t,
116 Float_t, Float_t, Float_t, Float_t, Int_t ireuse=0);
117
118 // TObject
119 void Print(Option_t *opt=NULL) const;
120 void Clear(Option_t *opt=NULL);
121
122 ClassDef(MMcEvt, 6) //Stores Montecarlo Information of one event (eg. the energy)
123};
124
125#endif
Note: See TracBrowser for help on using the repository browser.