source: trunk/MagicSoft/include-Classes/MMcFormat/MMcEvtBasic.h@ 9312

Last change on this file since 9312 was 9312, checked in by tbretz, 16 years ago
*** empty log message ***
File size: 3.7 KB
Line 
1#ifndef MARS_MMcEvtBasic
2#define MARS_MMcEvtBasic
3
4#ifndef MARS_MParContainer
5#include "MParContainer.h"
6#endif
7
8class MMcEvtBasic : public MParContainer
9{
10public:
11 enum ParticleId_t
12 {
13 kUNDEFINED = -1,
14 kGAMMA = 1,
15 kPOSITRON = 2,
16 kELECTRON = 3,
17 kANTIMUON = 5,
18 kMUON = 6,
19 kPI0 = 7,
20 kNEUTRON = 13,
21 kPROTON = 14,
22 kHELIUM = 402,
23 kOXYGEN = 1608,
24 kIRON = 5626,
25 kArtificial = 9998,
26 kNightSky = 9999
27 };
28
29protected:
30 ParticleId_t fPartId; // Type of particle
31 Float_t fEnergy; // [GeV] Energy
32 Float_t fImpact; // [cm] impact parameter
33
34 // Telescope orientation (see TDAS 02-11 regarding the
35 // precise meaning of these angles):
36 Float_t fTelescopePhi; // [rad]
37 Float_t fTelescopeTheta; // [rad]
38
39public:
40 MMcEvtBasic();
41 MMcEvtBasic(ParticleId_t, Float_t, Float_t, Float_t, Float_t);
42 void operator=(const MMcEvtBasic &evt);
43
44 // Getter
45 ParticleId_t GetPartId() const { return fPartId; }
46
47 Float_t GetEnergy() const { return fEnergy; }
48 Float_t GetImpact() const { return fImpact; }
49
50 Float_t GetTelescopePhi() const { return fTelescopePhi; }
51 Float_t GetTelescopeTheta() const { return fTelescopeTheta; }
52
53 static TString GetParticleName(Int_t id)
54 {
55 switch (id)
56 {
57 case kUNDEFINED: return "Undefined";
58 case kGAMMA: return "Gamma";
59 case kPOSITRON: return "Positron";
60 case kELECTRON: return "Electron";
61 case kANTIMUON: return "Anti-Muon";
62 case kMUON: return "Muon";
63 case kPI0: return "Pi-0";
64 case kNEUTRON: return "Neutron";
65 case kPROTON: return "Proton";
66 case kHELIUM: return "Helium";
67 case kOXYGEN: return "Oxygen";
68 case kIRON: return "Iron";
69 case kArtificial: return "Artificial";
70 case kNightSky: return "NightSky";
71 }
72
73 return Form("Id:%d", id);
74 }
75
76 static TString GetParticleSymbol(Int_t id)
77 {
78 switch (id)
79 {
80 case kUNDEFINED:return "N/A";
81 case kGAMMA: return "\\gamma";
82 case kPOSITRON: return "e^{+}";
83 case kELECTRON: return "e^{-}";
84 case kANTIMUON: return "\\mu^{+}";
85 case kMUON: return "\\mu^{-}";
86 case kPI0: return "\\pi^{0}";
87 case kNEUTRON: return "n";
88 case kPROTON: return "p";
89 case kHELIUM: return "He";
90 case kOXYGEN: return "O";
91 case kIRON: return "Fe";
92 case kNightSky: return "\\gamma_{NSB}";
93 }
94
95 return Form("Id:%d", id);
96 }
97
98 static TString GetEnergyStr(Float_t e)
99 {
100 if (e>=1000)
101 return Form("%.1fTeV", e/1000);
102
103 if (e>=10)
104 return Form("%dGeV", (Int_t)(e+.5));
105
106 if (e>=1)
107 return Form("%.1fGeV", e);
108
109 return Form("%dMeV", (Int_t)(e*1000+.5));
110 }
111
112 TString GetParticleSymbol() const
113 {
114 return GetParticleSymbol(fPartId);
115 }
116
117 TString GetParticleName() const
118 {
119 return GetParticleName(fPartId);
120 }
121
122 TString GetEnergyStr() const
123 {
124 return GetEnergyStr(fEnergy);
125 }
126
127 // Setter
128 void SetPartId(ParticleId_t id) { fPartId = id; }
129 void SetEnergy(Float_t Energy) { fEnergy=Energy; } //Set Energy
130 void SetImpact(Float_t Impact) { fImpact=Impact;} //Set impact parameter
131
132 void SetTelescopeTheta(Float_t Theta) { fTelescopeTheta=Theta; }
133 void SetTelescopePhi (Float_t Phi) { fTelescopePhi=Phi; }
134
135 void Fill(ParticleId_t, Float_t, Float_t, Float_t, Float_t);
136
137 // TObject
138 void Clear(Option_t *opt=NULL);
139 void Print(Option_t *opt=NULL) const;
140
141 ClassDef(MMcEvtBasic, 2) //Stores Basic Montecarlo Information of one event
142
143};
144
145#endif
Note: See TracBrowser for help on using the repository browser.