CaloProfile/inc/CaloProfile.h

/**
 * \file CaloProfile.h
 * \author Emiliano Mocchiutti
 */
#ifndef caloprofile_h
#define caloprofile_h

#define EMPTY -1
#define MAX 50

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <ctype.h>
#include <TTree.h>
#include <TFriendElement.h>
#include <TChain.h>
#include <TCanvas.h>
#include <TGraph.h>
#include <TH1F.h>
#include <TH2F.h>
#include <TFile.h>
#include <TPolyLine.h>
#include <TMath.h>
#include <TStyle.h>
#include <TList.h>
#include <TKey.h>
#include <TSystemFile.h>
#include <TSystemDirectory.h>
#include <TSQLServer.h>
#include <TF1.h>
#include <TGraphErrors.h>
#include <TMinuit.h>

#include <PamLevel2.h>

using namespace std;

struct stack{
    char data[MAX];
    int top;
};

/**
 *
 */
class CaloLat : public TObject {

 private:
    //
    PamLevel2 *L2;
    Bool_t debug;
    //
    // needed to avoid reprocessing the same event over and over to obtain the variables
    //
    UInt_t OBT;
    UInt_t PKT;
    UInt_t atime;
    //
    Float_t estrip[2][22][96];
    TString suf;
    Bool_t usepl18x;
    //

 public:
    //
    //
    void Draw();
    void Draw(Int_t,Int_t);
    //
    CaloLat();
    CaloLat(PamLevel2 *L2);
    ~CaloLat(){ Delete(); };
    //
    void SetDebug(Bool_t d){ debug=d; };
    //
    void Clear();
    void Clear(Option_t *option){Clear();};
    void Delete();
    void Delete(Option_t *option){Delete();};
    //
    void Process(); ///< Process data 
    void Print();
    void Print(Option_t *option){Print();};
    //
    
    void UsePlane18X(Bool_t use){usepl18x = use;};
    //
    //
    void SetSuffix(TString suffix){ suf = suffix;};
    //
    ClassDef(CaloLat,3);
};

/**
 *
 */
class CaloLong : public TObject {

 private:
    //
    PamLevel2 *L2;
    Bool_t debug;
    //
    // needed to avoid reprocessing the same event over and over to obtain the variables
    //
    UInt_t OBT;
    UInt_t PKT;
    UInt_t atime;
    UInt_t fOBT;
    UInt_t fPKT;
    UInt_t fatime;
    //
    Int_t N;
    Int_t NC;
    Bool_t sel;
    Bool_t cont;
    Int_t mask18b;
    Bool_t usepl18x;
    //
    Float_t chi2;
    Float_t ndf;
    Float_t E0;
    Float_t a;
    Float_t b;
    Float_t errE0;
    Float_t erra;
    Float_t errb;
    Float_t etmax;
    Float_t asymm;
    Float_t X0pl;
    Float_t defE0;
    Float_t umax;
    Float_t lmax;
    TString sumax;
    TString slmax;
    Int_t fitresult;
    //
    Bool_t no18x;
    Bool_t maskXE;
    Bool_t maskYE;
    Bool_t maskXO;
    Bool_t maskYO;
    //
    Bool_t xyaverage;
    //
    Bool_t heavytail;
    Float_t letmax;
    Float_t lmipth;
    //
    Float_t eplane[2][22];
    //
    CaloLevel2 *clp;
    //
    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"...
    //
    TString suf;
    TF1 *lfit;

 public:
    //
    //
    void Fit();
    void Fit(Bool_t draw);
    TF1 *GetFit();
    //
    //    Double_t ccurve(Double_t *t, Double_t *par); 
    //
    void SetCaloLevel2Pointer(CaloLevel2 *cp){ clp = cp;};
    //
    Float_t Get_E0(){this->Fit(); return E0;};
    Float_t Get_defE0(){this->Fit(); return defE0;};
    Float_t Get_a(){this->Fit(); return a;};
    Float_t Get_b(){this->Fit(); return b;};
    Float_t Get_errE0(){this->Fit(); return errE0;};
    Float_t Get_erra(){this->Fit(); return erra;};
    Float_t Get_errb(){this->Fit(); return errb;};
    Float_t Get_chi2(){this->Fit(); return chi2;};
    Float_t Get_ndf(){this->Fit(); return ndf;};
    Float_t Get_nchi2(){this->Fit(); if ( ndf > 0 ) return (chi2/ndf); return 0;};
    Float_t Get_tmax(){this->Fit(); if ( b != 0 ) return ((a-1.)/b); return 0;};
    Float_t Get_asymm(){this->Fit(); return asymm;};
    Float_t Get_exptmax(){this->Process(); return etmax;};
    Float_t Get_X0pl(){this->Process(); return X0pl;};
    //
    Float_t Get_letmax(){ return letmax;};
    Float_t Get_lmipth(){ return lmipth;};
    Int_t Get_fitresult(){this->Fit(); return fitresult;};    
    //
    void ForceNextFit(){atime=0;fatime=0;};
    void Draw();
    void Draw(Int_t);
    //
    CaloLong();
    CaloLong(PamLevel2 *L2);
    ~CaloLong(){ Delete(); };
    //
    void SetDebug(Bool_t d){ debug=d; };
    //    void UsePlane18X(){ no18x=false; };
    
    void UsePlane18X(Bool_t use="false"){usepl18x = use; no18x=!use;};
    //
    //
    void UseAverage(){ xyaverage = true;};
    void UseAllMeas(){ xyaverage = false;};
    //
    void MaskSection(TString);
    void UnMaskSection(TString);
    void UnMaskSections();
    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 Contamination(){sel = false; cont = true;}; ///< Set contamination mode: planes from N to 22 are used.
    void SetNoWpreSampler(Int_t n);
    void SetNoWcalo(Int_t n);
    void SplitInto(Int_t NoWpreSampler, Int_t NoWcalo);
    Int_t GetNoWpreSampler(){return N;}; ///< Get the number of W planes used as presampler. 
    Int_t GetNoWcalo(){return NC;}; ///< Get the number of W planes used as calorimeter.
    void SetEnergies(Float_t myene[][22]);
    //
    void SetLowerLimit(Float_t l){ lmax = l; };
    void SetUpperLimit(Float_t u){ umax = u; };
    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)
    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)
    //
    void Setletmax(Float_t l){ letmax = l;};
    void Setlmipth(Float_t l){ lmipth = l;};
    void HeavyTail(Bool_t b){ heavytail=b;};
    //
    Float_t GetLowerLimit(){ return lmax;};
    Float_t GetUpperLimit(){ return umax;};
    //
    void SetSuffix(TString suffix){ suf = suffix;};
    //
    void Clear();
    void Clear(Option_t *option){Clear();};
    void Delete();
    void Delete(Option_t *option){Delete();};
    //
    void Process(); ///< Process data 
    void Print();
    void Print(Option_t *option){Print();};
    //
    ClassDef(CaloLong,4);
};

/**
 *
 */
class Calo2D : public TObject {

 private:
    //
    PamLevel2 *L2;
    Bool_t debug;
    //
    // needed to avoid reprocessing the same event over and over to obtain the variables
    //
    UInt_t OBT;
    UInt_t PKT;
    UInt_t atime;
    //
    Float_t estrip[23][96][96];
    Int_t smax[23];
    Int_t smay[23];
    Bool_t usepl18x;
    //
    TString suf;

 public:
    //
    //
    void Draw();
    void Draw(Int_t);
    //
    Calo2D();
    Calo2D(PamLevel2 *L2);
    ~Calo2D(){ Delete(); };
    //
    void SetDebug(Bool_t d){ debug=d; };
    //
    void Clear();
    void Clear(Option_t *option){Clear();};
    void Delete();
    void Delete(Option_t *option){Delete();};
    //
    void SetSuffix(TString suffix){ suf = suffix;};
    
    void UsePlane18X(Bool_t use){usepl18x = use;};
    //
    //
    void Process(); ///< Process data 
    void Print();
    void Print(Option_t *option){Print();};
    //
    ClassDef(Calo2D,3);
};

#endif

1	mocchiut	1.1	/**
2			* \file CaloProfile.h
3			* \author Emiliano Mocchiutti
4			*/
5			#ifndef caloprofile_h
6			#define caloprofile_h
7
8	mocchiut	1.7	#define EMPTY -1
9			#define MAX 50
10
11	mocchiut	1.3	#include <iostream>
12	mocchiut	1.7	#include <stdio.h>
13			#include <string.h>
14			#include <ctype.h>
15	mocchiut	1.1	#include <TTree.h>
16			#include <TFriendElement.h>
17			#include <TChain.h>
18			#include <TCanvas.h>
19			#include <TGraph.h>
20			#include <TH1F.h>
21			#include <TH2F.h>
22			#include <TFile.h>
23			#include <TPolyLine.h>
24			#include <TMath.h>
25			#include <TStyle.h>
26			#include <TList.h>
27			#include <TKey.h>
28			#include <TSystemFile.h>
29			#include <TSystemDirectory.h>
30			#include <TSQLServer.h>
31			#include <TF1.h>
32	mocchiut	1.9	#include <TGraphErrors.h>
33	mocchiut	1.13	#include <TMinuit.h>
34	mocchiut	1.3
35			#include <PamLevel2.h>
36	mocchiut	1.1
37			using namespace std;
38
39	mocchiut	1.7	struct stack{
40			char data[MAX];
41			int top;
42			};
43
44	mocchiut	1.1	/**
45			*
46			*/
47			class CaloLat : public TObject {
48
49			private:
50			//
51			PamLevel2 *L2;
52			Bool_t debug;
53			//
54			// needed to avoid reprocessing the same event over and over to obtain the variables
55			//
56			UInt_t OBT;
57			UInt_t PKT;
58			UInt_t atime;
59			//
60			Float_t estrip[2][22][96];
61	mocchiut	1.11	TString suf;
62	mocchiut	1.14	Bool_t usepl18x;
63	mocchiut	1.1	//
64
65			public:
66			//
67			//
68			void Draw();
69			void Draw(Int_t,Int_t);
70			//
71			CaloLat();
72			CaloLat(PamLevel2 *L2);
73			~CaloLat(){ Delete(); };
74			//
75			void SetDebug(Bool_t d){ debug=d; };
76			//
77			void Clear();
78			void Clear(Option_t *option){Clear();};
79			void Delete();
80	mocchiut	1.2	void Delete(Option_t *option){Delete();};
81	mocchiut	1.1	//
82			void Process(); ///< Process data
83			void Print();
84	mocchiut	1.2	void Print(Option_t *option){Print();};
85	mocchiut	1.1	//
86	mocchiut	1.14
87			void UsePlane18X(Bool_t use){usepl18x = use;};
88			//
89			//
90	mocchiut	1.11	void SetSuffix(TString suffix){ suf = suffix;};
91			//
92	mocchiut	1.15	ClassDef(CaloLat,3);
93	mocchiut	1.1	};
94
95			/**
96			*
97			*/
98			class CaloLong : public TObject {
99
100			private:
101			//
102			PamLevel2 *L2;
103			Bool_t debug;
104			//
105			// needed to avoid reprocessing the same event over and over to obtain the variables
106			//
107			UInt_t OBT;
108			UInt_t PKT;
109			UInt_t atime;
110	mocchiut	1.3	UInt_t fOBT;
111			UInt_t fPKT;
112			UInt_t fatime;
113			//
114			Int_t N;
115			Int_t NC;
116			Bool_t sel;
117			Bool_t cont;
118			Int_t mask18b;
119	mocchiut	1.14	Bool_t usepl18x;
120	mocchiut	1.3	//
121			Float_t chi2;
122			Float_t ndf;
123			Float_t E0;
124			Float_t a;
125			Float_t b;
126			Float_t errE0;
127			Float_t erra;
128			Float_t errb;
129			Float_t etmax;
130			Float_t asymm;
131			Float_t X0pl;
132	mocchiut	1.7	Float_t defE0;
133			Float_t umax;
134			Float_t lmax;
135			TString sumax;
136			TString slmax;
137	mocchiut	1.3	Int_t fitresult;
138			//
139			Bool_t no18x;
140	mocchiut	1.5	Bool_t maskXE;
141			Bool_t maskYE;
142			Bool_t maskXO;
143			Bool_t maskYO;
144	mocchiut	1.1	//
145	mocchiut	1.9	Bool_t xyaverage;
146			//
147	mocchiut	1.12	Bool_t heavytail;
148			Float_t letmax;
149			Float_t lmipth;
150			//
151	mocchiut	1.1	Float_t eplane[2][22];
152			//
153	mocchiut	1.6	CaloLevel2 *clp;
154			//
155	mocchiut	1.10	Float_t Evaluate(TString s, Float_t tmax, Float_t X0pl ); // expression must be of the form "tmax+2.X0pl", "5tmax"."tmax+10","tmax-(4tmax)+3.X0pl"...
156	mocchiut	1.7	//
157	mocchiut	1.11	TString suf;
158	mocchiut	1.17	TF1 *lfit;
159	mocchiut	1.1
160			public:
161			//
162			//
163	mocchiut	1.3	void Fit();
164			void Fit(Bool_t draw);
165	mocchiut	1.16	TF1 *GetFit();
166	mocchiut	1.3	//
167			// Double_t ccurve(Double_t t, Double_t par);
168			//
169	mocchiut	1.6	void SetCaloLevel2Pointer(CaloLevel2 *cp){ clp = cp;};
170			//
171	mocchiut	1.3	Float_t Get_E0(){this->Fit(); return E0;};
172	mocchiut	1.7	Float_t Get_defE0(){this->Fit(); return defE0;};
173	mocchiut	1.3	Float_t Get_a(){this->Fit(); return a;};
174			Float_t Get_b(){this->Fit(); return b;};
175			Float_t Get_errE0(){this->Fit(); return errE0;};
176			Float_t Get_erra(){this->Fit(); return erra;};
177			Float_t Get_errb(){this->Fit(); return errb;};
178			Float_t Get_chi2(){this->Fit(); return chi2;};
179			Float_t Get_ndf(){this->Fit(); return ndf;};
180			Float_t Get_nchi2(){this->Fit(); if ( ndf > 0 ) return (chi2/ndf); return 0;};
181			Float_t Get_tmax(){this->Fit(); if ( b != 0 ) return ((a-1.)/b); return 0;};
182			Float_t Get_asymm(){this->Fit(); return asymm;};
183			Float_t Get_exptmax(){this->Process(); return etmax;};
184			Float_t Get_X0pl(){this->Process(); return X0pl;};
185	mocchiut	1.12	//
186			Float_t Get_letmax(){ return letmax;};
187			Float_t Get_lmipth(){ return lmipth;};
188			Int_t Get_fitresult(){this->Fit(); return fitresult;};
189	mocchiut	1.3	//
190	mocchiut	1.6	void ForceNextFit(){atime=0;fatime=0;};
191	mocchiut	1.1	void Draw();
192			void Draw(Int_t);
193			//
194			CaloLong();
195			CaloLong(PamLevel2 *L2);
196			~CaloLong(){ Delete(); };
197			//
198			void SetDebug(Bool_t d){ debug=d; };
199	mocchiut	1.14	// void UsePlane18X(){ no18x=false; };
200
201			void UsePlane18X(Bool_t use="false"){usepl18x = use; no18x=!use;};
202			//
203	mocchiut	1.3	//
204	mocchiut	1.9	void UseAverage(){ xyaverage = true;};
205			void UseAllMeas(){ xyaverage = false;};
206			//
207	mocchiut	1.5	void MaskSection(TString);
208	mocchiut	1.6	void UnMaskSection(TString);
209			void UnMaskSections();
210	mocchiut	1.3	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);
211			void Contamination(){sel = false; cont = true;}; ///< Set contamination mode: planes from N to 22 are used.
212			void SetNoWpreSampler(Int_t n);
213			void SetNoWcalo(Int_t n);
214	mocchiut	1.4	void SplitInto(Int_t NoWpreSampler, Int_t NoWcalo);
215	mocchiut	1.3	Int_t GetNoWpreSampler(){return N;}; ///< Get the number of W planes used as presampler.
216			Int_t GetNoWcalo(){return NC;}; ///< Get the number of W planes used as calorimeter.
217	mocchiut	1.6	void SetEnergies(Float_t myene[][22]);
218	mocchiut	1.3	//
219	mocchiut	1.7	void SetLowerLimit(Float_t l){ lmax = l; };
220			void SetUpperLimit(Float_t u){ umax = u; };
221			void SetLowerLimit(TString sl){ slmax = sl; };// expression must be of the form "5t"."t+10","t-(4t)"... where t will be replaced by the fitted maximum (X0)
222			void SetUpperLimit(TString su){ sumax = su; };// expression must be of the form "5t"."t+10","t-(4t)"... where t will be replaced by the fitted maximum (X0)
223			//
224	mocchiut	1.12	void Setletmax(Float_t l){ letmax = l;};
225			void Setlmipth(Float_t l){ lmipth = l;};
226			void HeavyTail(Bool_t b){ heavytail=b;};
227			//
228	mocchiut	1.8	Float_t GetLowerLimit(){ return lmax;};
229			Float_t GetUpperLimit(){ return umax;};
230			//
231	mocchiut	1.11	void SetSuffix(TString suffix){ suf = suffix;};
232			//
233	mocchiut	1.3	void Clear();
234			void Clear(Option_t *option){Clear();};
235			void Delete();
236			void Delete(Option_t *option){Delete();};
237			//
238			void Process(); ///< Process data
239			void Print();
240			void Print(Option_t *option){Print();};
241			//
242	mocchiut	1.15	ClassDef(CaloLong,4);
243	mocchiut	1.3	};
244
245			/**
246			*
247			*/
248			class Calo2D : public TObject {
249
250			private:
251			//
252			PamLevel2 *L2;
253			Bool_t debug;
254			//
255			// needed to avoid reprocessing the same event over and over to obtain the variables
256			//
257			UInt_t OBT;
258			UInt_t PKT;
259			UInt_t atime;
260			//
261			Float_t estrip[23][96][96];
262			Int_t smax[23];
263			Int_t smay[23];
264	mocchiut	1.14	Bool_t usepl18x;
265	mocchiut	1.3	//
266	mocchiut	1.11	TString suf;
267	mocchiut	1.3
268			public:
269			//
270			//
271			void Draw();
272			void Draw(Int_t);
273			//
274			Calo2D();
275			Calo2D(PamLevel2 *L2);
276			~Calo2D(){ Delete(); };
277			//
278			void SetDebug(Bool_t d){ debug=d; };
279	mocchiut	1.1	//
280			void Clear();
281			void Clear(Option_t *option){Clear();};
282			void Delete();
283	mocchiut	1.2	void Delete(Option_t *option){Delete();};
284	mocchiut	1.1	//
285	mocchiut	1.11	void SetSuffix(TString suffix){ suf = suffix;};
286	mocchiut	1.14
287			void UsePlane18X(Bool_t use){usepl18x = use;};
288			//
289	mocchiut	1.11	//
290	mocchiut	1.1	void Process(); ///< Process data
291			void Print();
292	mocchiut	1.2	void Print(Option_t *option){Print();};
293	mocchiut	1.1	//
294	mocchiut	1.15	ClassDef(Calo2D,3);
295	mocchiut	1.1	};
296
297			#endif
298