CaloBragg/inc/CaloBragg.h

/**
 * \file CaloBragg.h
 * \author Emiliano Mocchiutti & Valeria Malvezzi & Daniele Tavagnacco
 */
#ifndef calobragg_h
#define calobragg_h
#define SQ(a) ( a * a )
#include <PamLevel2.h>

#include <TTree.h>
#include <TFriendElement.h>
#include <TChain.h>
#include <TFile.h>
#include <TList.h>
#include <TKey.h>
#include <TSystemFile.h>
#include <TStyle.h>
#include <TSystemDirectory.h>
#include <TSQLServer.h>
#include <TCanvas.h>
#include <TH2F.h>
#include <iostream>

using namespace std;


/**
 *
 * Class to store and calculate variables useful for nuclei analysis
 */
class CaloBragg : public TObject {

 private:
    //
    PamLevel2 *L2;
    Bool_t debug;
    Bool_t usetrack;
    Bool_t startZero;

    //
    // needed to avoid reprocessing the same event over and over to obtain the variables
    //
    UInt_t OBT;
    UInt_t PKT;
    UInt_t atime;
    Int_t tr;
    Int_t sntr;
    //parametri calorimetro
    Int_t NPLA;
    Int_t NCHA;
    Int_t nView;
    
    Float_t AA;//mm larghezza strip
    Float_t ADIST;//mm distanza tra pad
    Float_t PIANO;//mm distanza 
    
    Float_t ySi;//mm spessore silicio
    Float_t yW;//mm spessore tungsteno
    Float_t rhoSi;//g/cm3 densita' silicio
    Float_t rhoW;//g/cm3 densita' tugsteno
    Float_t MIP;//Mev g/cm2 energia al minimo nel silicio per 0.38 mm
    
    Float_t emin;
    Float_t fzeta;
    Float_t fene;
        Int_t   fiso;
    //
    Bool_t usepl18x;
    Bool_t newchi2;
    Bool_t usenewBB;
    Int_t ndf;


    //parametri bethe-bloch
    Float_t pigr;
    Float_t Na;
    Float_t ZA; /*Z/A per Si*/
    Float_t ISi; /*MeV*/
    Float_t IW; /*MeV*/
    Float_t Me; /* MeV*/
    Float_t MassP;/*MeV*/
    Float_t r2; /*ro*ro in cm */
    
    //variabili globali caricabili
    Float_t Integrale;
    Float_t bestchi2[5];
/*     Float_t qtchi2; */
/*     Float_t qtz; */
/*     Float_t qtetot; */
/*     Float_t qtpskip; */
    Float_t lpchi2;
    Float_t lpz;
    Float_t lpisotope;
    Float_t lpetot;
    Float_t lppskip;
    Float_t calorimetro[44][2];
    Float_t spessore[4];
    Float_t estremi[2][2];//[0][0]primo piano [0][1]energia primo
                        //[1][0]ultim piano [1][1]energia ultim
    //Float_t shift[2][22];//shift dei piani
    //Float_t piano[22]; //posizione z dei piani
    Float_t elem[32][7];//A degli elementi compresi isotopi
    //

    Float_t dEplan[2*22];//energia rilasciata calcolata

/*     void SWAP(Float_t*,Float_t*); */
    void LoadParam();
    void conversione();
    void BetheBloch(Float_t *,Float_t *,Float_t *,Float_t *,Float_t *,Float_t *, Float_t);
    void ELOSS(Float_t *,Int_t *,Int_t *, Float_t *, Float_t *, Float_t);
    void Enetrack(Int_t *,Int_t *, Float_t *, Float_t *,Float_t *,Float_t * );
    void chiquadro(Float_t *, Float_t *);
    void loopze(Float_t ,Float_t ,Float_t ,Float_t, Int_t  );
    void mediatroncata();
    void Zdaloop();
    void CleanPlanes(Float_t epiano[22][2], Bool_t zpiano[22][2]);

 public:
    //
    void SetStartZero(Bool_t yes=false){ startZero = yes;};
    //
    Float_t Get_firstplane(){ Process(); return estremi[0][0];}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane) 
    Float_t Get_lastplane(){ Process(); return estremi[1][0];}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane) 
/*     Float_t Get_qtchi2(){ Process(); return qtchi2;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane)  */
/*     Float_t Get_qtz(){ Process(); return qtz;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane)  */
/*     Float_t Get_qtetot(){ Process(); return qtetot;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane)  */
/*     Float_t Get_qtpskip(){ Process(); return qtpskip;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane)  */
    Float_t Get_lpchi2(){ Process(); return lpchi2;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane) 
    Float_t Get_lpz(){ Process(); return lpz;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane) 
    Float_t Get_mass(){ Process(); return elem[(Int_t)lpz - 1][(Int_t)lpisotope];}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane) 
    Float_t Get_lpetot(){ Process(); return lpetot;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane) 
    Float_t Get_lppskip(){ Process(); return lppskip;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane) 

        Float_t Integral();
        Float_t LastIntegral();

    void ForceZeta(Float_t zeta){ fzeta = zeta;};
    void ForceEne(Float_t ener ){ fene = ener;};
    void ForceIso(Int_t isor ){ fiso = isor;};
    void ForceProcessing(){ sntr = -1000.;};

    void UsePlane18X(Bool_t use){usepl18x = use;};

    void UseNewChi2(Bool_t use){newchi2 = use;};
    void UseNewBB(Bool_t use){usenewBB = use;};
    //

    //
    CaloBragg();
    CaloBragg(PamLevel2 *L2);
    ~CaloBragg(){ Delete(); };
    //
    void Draw();
    void Draw(Int_t,Int_t, Float_t);
    void SetDebug(Bool_t d){ debug=d; };
    void UseTrack(Bool_t d){ usetrack=d; };
    //
    void Clear();
    void Clear(Option_t *option){Clear();};
    void Delete();
    void Delete(Option_t *option){Delete();};
    //
    void Process(); ///< Process data for track number 0
    void Process(Int_t ntr); ///< Process data for track number ntr
    void Print();
    void Print(Option_t *option){Print();};
    //
    ClassDef(CaloBragg,6);
};

#endif

1	pamelats	1.1	/**
2			* \file CaloBragg.h
3			* \author Emiliano Mocchiutti & Valeria Malvezzi & Daniele Tavagnacco
4			*/
5			#ifndef calobragg_h
6			#define calobragg_h
7			#define SQ(a) ( a * a )
8			#include <PamLevel2.h>
9
10			#include <TTree.h>
11			#include <TFriendElement.h>
12			#include <TChain.h>
13			#include <TFile.h>
14			#include <TList.h>
15			#include <TKey.h>
16			#include <TSystemFile.h>
17			#include <TStyle.h>
18			#include <TSystemDirectory.h>
19			#include <TSQLServer.h>
20			#include <TCanvas.h>
21			#include <TH2F.h>
22			#include <iostream>
23
24			using namespace std;
25
26
27			/**
28			*
29			* Class to store and calculate variables useful for nuclei analysis
30			*/
31			class CaloBragg : public TObject {
32
33			private:
34			//
35			PamLevel2 *L2;
36			Bool_t debug;
37			Bool_t usetrack;
38	mocchiut	1.7	Bool_t startZero;
39
40	pamelats	1.1	//
41			// needed to avoid reprocessing the same event over and over to obtain the variables
42			//
43			UInt_t OBT;
44			UInt_t PKT;
45			UInt_t atime;
46			Int_t tr;
47			Int_t sntr;
48			//parametri calorimetro
49			Int_t NPLA;
50			Int_t NCHA;
51			Int_t nView;
52
53			Float_t AA;//mm larghezza strip
54			Float_t ADIST;//mm distanza tra pad
55			Float_t PIANO;//mm distanza
56
57			Float_t ySi;//mm spessore silicio
58			Float_t yW;//mm spessore tungsteno
59			Float_t rhoSi;//g/cm3 densita' silicio
60			Float_t rhoW;//g/cm3 densita' tugsteno
61			Float_t MIP;//Mev g/cm2 energia al minimo nel silicio per 0.38 mm
62
63			Float_t emin;
64	mocchiut	1.7	Float_t fzeta;
65	mocchiut	1.10	Float_t fene;
66	biancoa	1.12	Int_t fiso;
67	mocchiut	1.4	//
68			Bool_t usepl18x;
69	mocchiut	1.7	Bool_t newchi2;
70			Bool_t usenewBB;
71			Int_t ndf;
72
73	mocchiut	1.4
74	pamelats	1.1	//parametri bethe-bloch
75			Float_t pigr;
76			Float_t Na;
77			Float_t ZA; /Z/A per Si/
78			Float_t ISi; /MeV/
79	mocchiut	1.7	Float_t IW; /MeV/
80	pamelats	1.1	Float_t Me; /* MeV*/
81			Float_t MassP;/MeV/
82			Float_t r2; /roro in cm */
83
84			//variabili globali caricabili
85			Float_t Integrale;
86	mocchiut	1.10	Float_t bestchi2[5];
87	pamelats	1.3	/* Float_t qtchi2; */
88			/* Float_t qtz; */
89			/* Float_t qtetot; */
90			/* Float_t qtpskip; */
91	pamelats	1.1	Float_t lpchi2;
92			Float_t lpz;
93	mocchiut	1.10	Float_t lpisotope;
94	pamelats	1.1	Float_t lpetot;
95			Float_t lppskip;
96			Float_t calorimetro[44][2];
97	mocchiut	1.8	Float_t spessore[4];
98	pamelats	1.1	Float_t estremi[2][2];//[0][0]primo piano [0][1]energia primo
99			//[1][0]ultim piano [1][1]energia ultim
100			//Float_t shift[2][22];//shift dei piani
101			//Float_t piano[22]; //posizione z dei piani
102	mocchiut	1.10	Float_t elem[32][7];//A degli elementi compresi isotopi
103	pamelats	1.1	//
104	mocchiut	1.7
105			Float_t dEplan[2*22];//energia rilasciata calcolata
106
107	pamelats	1.2	/* void SWAP(Float_t,Float_t); */
108	pamelats	1.1	void LoadParam();
109			void conversione();
110	mocchiut	1.7	void BetheBloch(Float_t ,Float_t ,Float_t ,Float_t ,Float_t ,Float_t , Float_t);
111	mocchiut	1.10	void ELOSS(Float_t ,Int_t ,Int_t , Float_t , Float_t *, Float_t);
112			void Enetrack(Int_t ,Int_t , Float_t , Float_t ,Float_t ,Float_t );
113	pamelats	1.1	void chiquadro(Float_t , Float_t );
114	mocchiut	1.6	void loopze(Float_t ,Float_t ,Float_t ,Float_t, Int_t );
115	pamelats	1.1	void mediatroncata();
116			void Zdaloop();
117	mocchiut	1.11	void CleanPlanes(Float_t epiano[22][2], Bool_t zpiano[22][2]);
118	pamelats	1.1
119			public:
120			//
121	mocchiut	1.7	void SetStartZero(Bool_t yes=false){ startZero = yes;};
122	pamelats	1.1	//
123	pamelats	1.3	Float_t Get_firstplane(){ Process(); return estremi[0][0];}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane)
124			Float_t Get_lastplane(){ Process(); return estremi[1][0];}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane)
125			/* Float_t Get_qtchi2(){ Process(); return qtchi2;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane) */
126			/* Float_t Get_qtz(){ Process(); return qtz;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane) */
127			/* Float_t Get_qtetot(){ Process(); return qtetot;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane) */
128			/* Float_t Get_qtpskip(){ Process(); return qtpskip;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane) */
129	pamelats	1.1	Float_t Get_lpchi2(){ Process(); return lpchi2;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane)
130			Float_t Get_lpz(){ Process(); return lpz;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane)
131	mocchiut	1.10	Float_t Get_mass(){ Process(); return elem[(Int_t)lpz - 1][(Int_t)lpisotope];}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane)
132	pamelats	1.1	Float_t Get_lpetot(){ Process(); return lpetot;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane)
133			Float_t Get_lppskip(){ Process(); return lppskip;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane)
134
135	mocchiut	1.7	Float_t Integral();
136			Float_t LastIntegral();
137
138			void ForceZeta(Float_t zeta){ fzeta = zeta;};
139	mocchiut	1.10	void ForceEne(Float_t ener ){ fene = ener;};
140	biancoa	1.12	void ForceIso(Int_t isor ){ fiso = isor;};
141	mocchiut	1.7	void ForceProcessing(){ sntr = -1000.;};
142
143	mocchiut	1.4	void UsePlane18X(Bool_t use){usepl18x = use;};
144	mocchiut	1.7
145			void UseNewChi2(Bool_t use){newchi2 = use;};
146			void UseNewBB(Bool_t use){usenewBB = use;};
147	mocchiut	1.4	//
148
149	pamelats	1.1	//
150			CaloBragg();
151			CaloBragg(PamLevel2 *L2);
152			~CaloBragg(){ Delete(); };
153			//
154			void Draw();
155	mocchiut	1.10	void Draw(Int_t,Int_t, Float_t);
156	pamelats	1.1	void SetDebug(Bool_t d){ debug=d; };
157			void UseTrack(Bool_t d){ usetrack=d; };
158			//
159			void Clear();
160			void Clear(Option_t *option){Clear();};
161			void Delete();
162			void Delete(Option_t *option){Delete();};
163			//
164			void Process(); ///< Process data for track number 0
165			void Process(Int_t ntr); ///< Process data for track number ntr
166			void Print();
167			void Print(Option_t *option){Print();};
168			//
169	mocchiut	1.7	ClassDef(CaloBragg,6);
170	pamelats	1.1	};
171
172			#endif
173