9 |
#include <TClonesArray.h> |
#include <TClonesArray.h> |
10 |
#include <TArrayI.h> |
#include <TArrayI.h> |
11 |
// |
// |
12 |
|
#include <math.h> |
13 |
|
// |
14 |
#include <CaloStruct.h> |
#include <CaloStruct.h> |
15 |
// |
// |
16 |
|
|
40 |
Int_t ntr; ///< the same as "ncyl" but with radius 4 strips |
Int_t ntr; ///< the same as "ncyl" but with radius 4 strips |
41 |
Int_t planetot; ///< number of planes used to calculate the energy truncated mean "qmean" |
Int_t planetot; ///< number of planes used to calculate the energy truncated mean "qmean" |
42 |
Int_t nlow; ///< the same as "nstrip" but only after the calculated electromagnetic shower maximum |
Int_t nlow; ///< the same as "nstrip" but only after the calculated electromagnetic shower maximum |
43 |
Int_t tibar[22][2]; ///< strip traversed by the trajectory as measured by the tracker |
Int_t tibar[22][2]; ///< strip traversed by the trajectory as measured by the tracker or by the selftrigger when trkseqno = -1 |
44 |
Float_t tbar[22][2]; ///< position in cm as measured by the tracker |
Float_t tbar[22][2]; ///< position in cm as measured by the tracker or by the selftrigger when trkseqno = -1 |
45 |
Float_t qcore; ///< SUM(j=1,2)SUM(i=1,PLmax) Qhit(i,j)*i , where Qhit(i,j) is the energy released (MIP) in a cylinder of radius 2 Rm (Moliere radius) around the track in the i-th plane (where the top plane is number 1 and the sum runs up to plane number PLmax, closest to the calculated electromagnetic shower maximum of the j-th view). |
Float_t qcore; ///< SUM(j=1,2)SUM(i=1,PLmax) Qhit(i,j)*i , where Qhit(i,j) is the energy released (MIP) in a cylinder of radius 2 Rm (Moliere radius) around the track in the i-th plane (where the top plane is number 1 and the sum runs up to plane number PLmax, closest to the calculated electromagnetic shower maximum of the j-th view). |
46 |
Float_t qcyl; ///< the measured energy deposited in a cylinder of radius 8 strips around the shower axis |
Float_t qcyl; ///< the measured energy deposited in a cylinder of radius 8 strips around the shower axis |
47 |
Float_t qlast; ///< the same as "qcyl" but only for the last four planes and radius 4 strips. |
Float_t qlast; ///< the same as "qcyl" but only for the last four planes and radius 4 strips. |
63 |
*/ |
*/ |
64 |
CaloTrkVar(const CaloTrkVar &trkvar); ///< copy values from trkvar to this |
CaloTrkVar(const CaloTrkVar &trkvar); ///< copy values from trkvar to this |
65 |
// |
// |
66 |
void Clear(); ///< clear variables |
void Clear(Option_t *t=""); ///< clear variables |
67 |
CaloTrkVar* GetCaloTrkVar(){return this;}; ///< returns pointer to this object |
CaloTrkVar* GetCaloTrkVar(){return this;}; ///< returns pointer to this object |
68 |
// |
// |
69 |
ClassDef(CaloTrkVar,2); |
ClassDef(CaloTrkVar,2); |
88 |
Int_t perr[4]; ///< processing errors (one for each calorimeter section) |
Int_t perr[4]; ///< processing errors (one for each calorimeter section) |
89 |
Int_t swerr[4];///< DSP status word |
Int_t swerr[4];///< DSP status word |
90 |
Int_t crc[4]; ///< CRC errors on data |
Int_t crc[4]; ///< CRC errors on data |
91 |
Int_t selftrigger;///< self-trigger flag (1 selftrigger event, 0 normal event) |
Int_t selftrigger;///< self-trigger flag: 0 tof trigger, 2 selftrigger event, 3 selftrigger + tof trigger, 102 selftrigger event not stored in the selfdelay array NB: selftrigger == S4 trigger || no trigger || calo trigger in the triggerconf array |
92 |
// |
// |
93 |
// common variables (not related to tracks) |
// common variables (not related to tracks) |
94 |
// |
// |
95 |
Int_t nstrip; ///< total number of strip hit |
Int_t nstrip; ///< total number of strip hit |
96 |
|
Int_t nsatstrip; ///< total number of strip hit with saturated signal |
97 |
Int_t nx22; ///< number of strip hit in the last silicon plane of the calorimeter (x view number 22) |
Int_t nx22; ///< number of strip hit in the last silicon plane of the calorimeter (x view number 22) |
98 |
Int_t planemax[2]; ///< plane of maximum energy release (x and y) |
Int_t planemax[2]; ///< plane of maximum energy release (x and y) |
99 |
|
Int_t selfdelay[4][7]; ///< Delay of the selftrigger planes + coincidence for every section |
100 |
Float_t qtot; ///< total energy detected (MIP) |
Float_t qtot; ///< total energy detected (MIP) |
101 |
Float_t qx22; ///< energy detected in the last silicon plane of the calorimeter (x view number 22) |
Float_t qx22; ///< energy detected in the last silicon plane of the calorimeter (x view number 22) |
102 |
Float_t qmax; ///< the maximum energy detected in a strip |
Float_t qmax; ///< the maximum energy detected in a strip |
104 |
// |
// |
105 |
// Fit variables |
// Fit variables |
106 |
// |
// |
107 |
Int_t npcfit[2]; ///< number of point used to perform the fit for the two views |
Int_t npcfit[4]; ///< number of point used to perform the fit for the two views (0,1 calo fit, 2,3 selftrigger fit if any) |
108 |
Int_t cibar[22][2]; ///< strip traversed by the trajectory as measured by the calorimeter |
Float_t varcfit[4]; ///< variance of the calorimeter fit for the two views (0,1 calo fit, 2,3 selftrigger fit if any) |
109 |
Float_t cbar[22][2]; ///< position in cm as measured by the calorimeter |
Float_t tanx[2]; ///< the tangent of the angle in the x direction as determined by the track fitted in the calorimeter (0 calo fit, 1 selftrigger fit) |
110 |
Float_t impx; ///< the x impact position on the first plane as determined by the track fitted in the calorimeter |
Float_t tany[2]; ///< the tangent of the angle in the x direction as determined by the track fitted in the calorimeter (0 calo fit, 1 selftrigger fit) |
111 |
Float_t impy; ///< the y impact position on the first plane as determined by the track fitted in the calorimeter |
Int_t fitmode[2]; ///< for x and y is 0 if the fit was performed with the "electron" algorithm, is 1 if the fit was performed with the "nuclei" algorithm |
112 |
Float_t tanx; ///< the tangent of the angle in the x direction as determined by the track fitted in the calorimeter |
Int_t cibar[22][2]; ///< strip traversed by the trajectory as measured by the calorimeter (calo fit) |
113 |
Float_t tany; ///< the tangent of the angle in the x direction as determined by the track fitted in the calorimeter |
Float_t cbar[22][2]; ///< position in cm as measured by the calorimeter (calo fit) |
|
Float_t varcfit[2]; ///< variance of the calorimeter fit for the two views |
|
114 |
// |
// |
115 |
// Energy variables |
// Energy variables |
116 |
// |
// |
119 |
// |
// |
120 |
// track related variables: inline methods |
// track related variables: inline methods |
121 |
// |
// |
122 |
Int_t ntrk() {return CaloTrk->GetEntries();}; ///< number of saved blocks of track-related variables |
Int_t ntrk(){return CaloTrk->GetEntries();}; ///< number of saved blocks of track-related variables |
123 |
// |
// |
124 |
// METHODS |
// METHODS |
125 |
// |
// |
126 |
|
Float_t impx(Int_t tr); ///< the x impact position on the first plane as determined by the track fitted in the calorimeter ( tr = 0 calo fit, tr = 1 selftrigger fit) |
127 |
|
Float_t impy(Int_t tr); ///< the y impact position on the first plane as determined by the track fitted in the calorimeter ( tr = 0 calo fit, tr = 1 selftrigger fit) |
128 |
|
|
129 |
|
// |
130 |
void GetElectronEnergy(Float_t &energy, Float_t &sigma); ///< returns energy and sigma using qtot and assuming the particle being an electron |
void GetElectronEnergy(Float_t &energy, Float_t &sigma); ///< returns energy and sigma using qtot and assuming the particle being an electron |
131 |
// |
// |
132 |
CaloTrkVar *GetCaloTrkVar(Int_t notrack); ///< returns a pointer to the CaloTrkVar class containing track related variables for track number notrack |
CaloTrkVar *GetCaloTrkVar(Int_t notrack); ///< returns a pointer to the CaloTrkVar class containing track related variables for track number notrack |
133 |
|
CaloTrkVar* GetCaloStoredTrack(Int_t seqno); ///< returns pointer to the track set related to the seqno number |
134 |
// |
// |
135 |
TClonesArray *GetTrackArray(){return CaloTrk;}; ///< returns a pointer to the track related variables array |
TClonesArray *GetTrackArray(){return CaloTrk;}; ///< returns a pointer to the track related variables array |
136 |
CaloLevel2* GetCaloLevel2(){return this;}; ///< returns pointer to this object |
CaloLevel2* GetCaloLevel2(){return this;}; ///< returns pointer to this object |
137 |
// |
// |
138 |
void GetLevel2Struct(cCaloLevel2 *l2) const; |
void GetLevel2Struct(cCaloLevel2 *l2) const; |
139 |
// |
// |
140 |
void Clear(); |
void Clear(Option_t *t=""); |
141 |
|
void Delete(Option_t *t=""); //ELENA |
142 |
|
void Set(); //ELENA |
143 |
|
Bool_t IsGood(Bool_t strict = false); // method to check if the event is good or not |
144 |
// |
// |
145 |
// constructor |
// constructor |
146 |
// |
// |
147 |
CaloLevel2(); ///< Constructor. |
CaloLevel2(); ///< Constructor. |
148 |
|
~CaloLevel2(){Delete();}; //ELENA |
149 |
// |
// |
150 |
friend class CaloProcessing; |
friend class CaloLevel0; |
151 |
// |
// |
152 |
ClassDef(CaloLevel2,3); |
ClassDef(CaloLevel2,6); |
153 |
}; |
}; |
154 |
|
|
155 |
#endif |
#endif |