5 |
#ifndef caloprofile_h |
#ifndef caloprofile_h |
6 |
#define caloprofile_h |
#define caloprofile_h |
7 |
|
|
8 |
|
#define EMPTY -1 |
9 |
|
#define MAX 50 |
10 |
|
|
11 |
#include <iostream> |
#include <iostream> |
12 |
|
#include <stdio.h> |
13 |
|
#include <string.h> |
14 |
|
#include <ctype.h> |
15 |
#include <TTree.h> |
#include <TTree.h> |
16 |
#include <TFriendElement.h> |
#include <TFriendElement.h> |
17 |
#include <TChain.h> |
#include <TChain.h> |
29 |
#include <TSystemDirectory.h> |
#include <TSystemDirectory.h> |
30 |
#include <TSQLServer.h> |
#include <TSQLServer.h> |
31 |
#include <TF1.h> |
#include <TF1.h> |
32 |
|
#include <TGraphErrors.h> |
33 |
|
#include <TMinuit.h> |
34 |
|
|
35 |
#include <PamLevel2.h> |
#include <PamLevel2.h> |
36 |
|
|
37 |
using namespace std; |
using namespace std; |
38 |
|
|
39 |
|
struct stack{ |
40 |
|
char data[MAX]; |
41 |
|
int top; |
42 |
|
}; |
43 |
|
|
44 |
/** |
/** |
45 |
* |
* |
46 |
*/ |
*/ |
58 |
UInt_t atime; |
UInt_t atime; |
59 |
// |
// |
60 |
Float_t estrip[2][22][96]; |
Float_t estrip[2][22][96]; |
61 |
|
TString suf; |
62 |
|
Bool_t usepl18x; |
63 |
// |
// |
64 |
|
|
65 |
public: |
public: |
83 |
void Print(); |
void Print(); |
84 |
void Print(Option_t *option){Print();}; |
void Print(Option_t *option){Print();}; |
85 |
// |
// |
86 |
ClassDef(CaloLat,1); |
|
87 |
|
void UsePlane18X(Bool_t use){usepl18x = use;}; |
88 |
|
// |
89 |
|
// |
90 |
|
void SetSuffix(TString suffix){ suf = suffix;}; |
91 |
|
// |
92 |
|
ClassDef(CaloLat,2); |
93 |
}; |
}; |
94 |
|
|
95 |
/** |
/** |
116 |
Bool_t sel; |
Bool_t sel; |
117 |
Bool_t cont; |
Bool_t cont; |
118 |
Int_t mask18b; |
Int_t mask18b; |
119 |
|
Bool_t usepl18x; |
120 |
// |
// |
121 |
Float_t chi2; |
Float_t chi2; |
122 |
Float_t ndf; |
Float_t ndf; |
129 |
Float_t etmax; |
Float_t etmax; |
130 |
Float_t asymm; |
Float_t asymm; |
131 |
Float_t X0pl; |
Float_t X0pl; |
132 |
|
Float_t defE0; |
133 |
|
Float_t umax; |
134 |
|
Float_t lmax; |
135 |
|
TString sumax; |
136 |
|
TString slmax; |
137 |
Int_t fitresult; |
Int_t fitresult; |
138 |
// |
// |
139 |
Bool_t no18x; |
Bool_t no18x; |
142 |
Bool_t maskXO; |
Bool_t maskXO; |
143 |
Bool_t maskYO; |
Bool_t maskYO; |
144 |
// |
// |
145 |
|
Bool_t xyaverage; |
146 |
|
// |
147 |
|
Bool_t heavytail; |
148 |
|
Float_t letmax; |
149 |
|
Float_t lmipth; |
150 |
|
// |
151 |
Float_t eplane[2][22]; |
Float_t eplane[2][22]; |
152 |
// |
// |
153 |
CaloLevel2 *clp; |
CaloLevel2 *clp; |
154 |
// |
// |
155 |
|
Float_t Evaluate(TString s, Float_t tmax, Float_t X0pl ); // expression must be of the form "tmax+2.*X0pl", "5*tmax"."tmax+10","tmax-(4*tmax)+3.*X0pl"... |
156 |
|
// |
157 |
|
TString suf; |
158 |
|
|
159 |
public: |
public: |
160 |
// |
// |
167 |
void SetCaloLevel2Pointer(CaloLevel2 *cp){ clp = cp;}; |
void SetCaloLevel2Pointer(CaloLevel2 *cp){ clp = cp;}; |
168 |
// |
// |
169 |
Float_t Get_E0(){this->Fit(); return E0;}; |
Float_t Get_E0(){this->Fit(); return E0;}; |
170 |
|
Float_t Get_defE0(){this->Fit(); return defE0;}; |
171 |
Float_t Get_a(){this->Fit(); return a;}; |
Float_t Get_a(){this->Fit(); return a;}; |
172 |
Float_t Get_b(){this->Fit(); return b;}; |
Float_t Get_b(){this->Fit(); return b;}; |
173 |
Float_t Get_errE0(){this->Fit(); return errE0;}; |
Float_t Get_errE0(){this->Fit(); return errE0;}; |
180 |
Float_t Get_asymm(){this->Fit(); return asymm;}; |
Float_t Get_asymm(){this->Fit(); return asymm;}; |
181 |
Float_t Get_exptmax(){this->Process(); return etmax;}; |
Float_t Get_exptmax(){this->Process(); return etmax;}; |
182 |
Float_t Get_X0pl(){this->Process(); return X0pl;}; |
Float_t Get_X0pl(){this->Process(); return X0pl;}; |
183 |
Int_t Get_fitresult(){this->Fit(); return fitresult;}; |
// |
184 |
|
Float_t Get_letmax(){ return letmax;}; |
185 |
|
Float_t Get_lmipth(){ return lmipth;}; |
186 |
|
Int_t Get_fitresult(){this->Fit(); return fitresult;}; |
187 |
// |
// |
188 |
void ForceNextFit(){atime=0;fatime=0;}; |
void ForceNextFit(){atime=0;fatime=0;}; |
189 |
void Draw(); |
void Draw(); |
194 |
~CaloLong(){ Delete(); }; |
~CaloLong(){ Delete(); }; |
195 |
// |
// |
196 |
void SetDebug(Bool_t d){ debug=d; }; |
void SetDebug(Bool_t d){ debug=d; }; |
197 |
void UsePlane18X(){ no18x=false; }; |
// void UsePlane18X(){ no18x=false; }; |
198 |
|
|
199 |
|
void UsePlane18X(Bool_t use="false"){usepl18x = use; no18x=!use;}; |
200 |
|
// |
201 |
|
// |
202 |
|
void UseAverage(){ xyaverage = true;}; |
203 |
|
void UseAllMeas(){ xyaverage = false;}; |
204 |
// |
// |
205 |
void MaskSection(TString); |
void MaskSection(TString); |
206 |
void UnMaskSection(TString); |
void UnMaskSection(TString); |
214 |
Int_t GetNoWcalo(){return NC;}; ///< Get the number of W planes used as calorimeter. |
Int_t GetNoWcalo(){return NC;}; ///< Get the number of W planes used as calorimeter. |
215 |
void SetEnergies(Float_t myene[][22]); |
void SetEnergies(Float_t myene[][22]); |
216 |
// |
// |
217 |
|
void SetLowerLimit(Float_t l){ lmax = l; }; |
218 |
|
void SetUpperLimit(Float_t u){ umax = u; }; |
219 |
|
void SetLowerLimit(TString sl){ slmax = sl; };// expression must be of the form "5*t"."t+10","t-(4*t)"... where t will be replaced by the fitted maximum (X0) |
220 |
|
void SetUpperLimit(TString su){ sumax = su; };// expression must be of the form "5*t"."t+10","t-(4*t)"... where t will be replaced by the fitted maximum (X0) |
221 |
|
// |
222 |
|
void Setletmax(Float_t l){ letmax = l;}; |
223 |
|
void Setlmipth(Float_t l){ lmipth = l;}; |
224 |
|
void HeavyTail(Bool_t b){ heavytail=b;}; |
225 |
|
// |
226 |
|
Float_t GetLowerLimit(){ return lmax;}; |
227 |
|
Float_t GetUpperLimit(){ return umax;}; |
228 |
|
// |
229 |
|
void SetSuffix(TString suffix){ suf = suffix;}; |
230 |
|
// |
231 |
void Clear(); |
void Clear(); |
232 |
void Clear(Option_t *option){Clear();}; |
void Clear(Option_t *option){Clear();}; |
233 |
void Delete(); |
void Delete(); |
237 |
void Print(); |
void Print(); |
238 |
void Print(Option_t *option){Print();}; |
void Print(Option_t *option){Print();}; |
239 |
// |
// |
240 |
ClassDef(CaloLong,2); |
ClassDef(CaloLong,3); |
241 |
}; |
}; |
242 |
|
|
243 |
/** |
/** |
259 |
Float_t estrip[23][96][96]; |
Float_t estrip[23][96][96]; |
260 |
Int_t smax[23]; |
Int_t smax[23]; |
261 |
Int_t smay[23]; |
Int_t smay[23]; |
262 |
|
Bool_t usepl18x; |
263 |
// |
// |
264 |
|
TString suf; |
265 |
|
|
266 |
public: |
public: |
267 |
// |
// |
280 |
void Delete(); |
void Delete(); |
281 |
void Delete(Option_t *option){Delete();}; |
void Delete(Option_t *option){Delete();}; |
282 |
// |
// |
283 |
|
void SetSuffix(TString suffix){ suf = suffix;}; |
284 |
|
|
285 |
|
void UsePlane18X(Bool_t use){usepl18x = use;}; |
286 |
|
// |
287 |
|
// |
288 |
void Process(); ///< Process data |
void Process(); ///< Process data |
289 |
void Print(); |
void Print(); |
290 |
void Print(Option_t *option){Print();}; |
void Print(Option_t *option){Print();}; |
291 |
// |
// |
292 |
ClassDef(Calo2D,1); |
ClassDef(Calo2D,2); |
293 |
}; |
}; |
294 |
|
|
295 |
#endif |
#endif |