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 |
|
|
34 |
#include <PamLevel2.h> |
#include <PamLevel2.h> |
35 |
|
|
36 |
using namespace std; |
using namespace std; |
37 |
|
|
38 |
|
struct stack{ |
39 |
|
char data[MAX]; |
40 |
|
int top; |
41 |
|
}; |
42 |
|
|
43 |
/** |
/** |
44 |
* |
* |
45 |
*/ |
*/ |
119 |
Float_t etmax; |
Float_t etmax; |
120 |
Float_t asymm; |
Float_t asymm; |
121 |
Float_t X0pl; |
Float_t X0pl; |
122 |
|
Float_t defE0; |
123 |
|
Float_t umax; |
124 |
|
Float_t lmax; |
125 |
|
TString sumax; |
126 |
|
TString slmax; |
127 |
Int_t fitresult; |
Int_t fitresult; |
128 |
// |
// |
129 |
Bool_t no18x; |
Bool_t no18x; |
132 |
Bool_t maskXO; |
Bool_t maskXO; |
133 |
Bool_t maskYO; |
Bool_t maskYO; |
134 |
// |
// |
135 |
|
Bool_t xyaverage; |
136 |
|
// |
137 |
Float_t eplane[2][22]; |
Float_t eplane[2][22]; |
138 |
// |
// |
139 |
|
CaloLevel2 *clp; |
140 |
|
// |
141 |
|
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"... |
142 |
|
// |
143 |
|
|
144 |
public: |
public: |
145 |
// |
// |
149 |
// |
// |
150 |
// Double_t ccurve(Double_t *t, Double_t *par); |
// Double_t ccurve(Double_t *t, Double_t *par); |
151 |
// |
// |
152 |
|
void SetCaloLevel2Pointer(CaloLevel2 *cp){ clp = cp;}; |
153 |
|
// |
154 |
Float_t Get_E0(){this->Fit(); return E0;}; |
Float_t Get_E0(){this->Fit(); return E0;}; |
155 |
|
Float_t Get_defE0(){this->Fit(); return defE0;}; |
156 |
Float_t Get_a(){this->Fit(); return a;}; |
Float_t Get_a(){this->Fit(); return a;}; |
157 |
Float_t Get_b(){this->Fit(); return b;}; |
Float_t Get_b(){this->Fit(); return b;}; |
158 |
Float_t Get_errE0(){this->Fit(); return errE0;}; |
Float_t Get_errE0(){this->Fit(); return errE0;}; |
167 |
Float_t Get_X0pl(){this->Process(); return X0pl;}; |
Float_t Get_X0pl(){this->Process(); return X0pl;}; |
168 |
Int_t Get_fitresult(){this->Fit(); return fitresult;}; |
Int_t Get_fitresult(){this->Fit(); return fitresult;}; |
169 |
// |
// |
170 |
|
void ForceNextFit(){atime=0;fatime=0;}; |
171 |
void Draw(); |
void Draw(); |
172 |
void Draw(Int_t); |
void Draw(Int_t); |
173 |
// |
// |
178 |
void SetDebug(Bool_t d){ debug=d; }; |
void SetDebug(Bool_t d){ debug=d; }; |
179 |
void UsePlane18X(){ no18x=false; }; |
void UsePlane18X(){ no18x=false; }; |
180 |
// |
// |
181 |
|
void UseAverage(){ xyaverage = true;}; |
182 |
|
void UseAllMeas(){ xyaverage = false;}; |
183 |
|
// |
184 |
void MaskSection(TString); |
void MaskSection(TString); |
185 |
|
void UnMaskSection(TString); |
186 |
|
void UnMaskSections(); |
187 |
void Selection(){sel = true; cont = false;}; ///< Set selection mode: planes from 1 to 22-N are used, plane 18 - N is masked if "emulate18" variable is true (DEFAULT); |
void Selection(){sel = true; cont = false;}; ///< Set selection mode: planes from 1 to 22-N are used, plane 18 - N is masked if "emulate18" variable is true (DEFAULT); |
188 |
void Contamination(){sel = false; cont = true;}; ///< Set contamination mode: planes from N to 22 are used. |
void Contamination(){sel = false; cont = true;}; ///< Set contamination mode: planes from N to 22 are used. |
189 |
void SetNoWpreSampler(Int_t n); |
void SetNoWpreSampler(Int_t n); |
191 |
void SplitInto(Int_t NoWpreSampler, Int_t NoWcalo); |
void SplitInto(Int_t NoWpreSampler, Int_t NoWcalo); |
192 |
Int_t GetNoWpreSampler(){return N;}; ///< Get the number of W planes used as presampler. |
Int_t GetNoWpreSampler(){return N;}; ///< Get the number of W planes used as presampler. |
193 |
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. |
194 |
|
void SetEnergies(Float_t myene[][22]); |
195 |
|
// |
196 |
|
void SetLowerLimit(Float_t l){ lmax = l; }; |
197 |
|
void SetUpperLimit(Float_t u){ umax = u; }; |
198 |
|
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) |
199 |
|
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) |
200 |
|
// |
201 |
|
Float_t GetLowerLimit(){ return lmax;}; |
202 |
|
Float_t GetUpperLimit(){ return umax;}; |
203 |
// |
// |
204 |
void Clear(); |
void Clear(); |
205 |
void Clear(Option_t *option){Clear();}; |
void Clear(Option_t *option){Clear();}; |