source: trunk/Mars/mmc/MMcEvt.hxx@ 19992

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