35 |
PamLevel2 *L2; |
PamLevel2 *L2; |
36 |
Bool_t debug; |
Bool_t debug; |
37 |
Bool_t usetrack; |
Bool_t usetrack; |
38 |
|
Bool_t startZero; |
39 |
|
|
40 |
// |
// |
41 |
// needed to avoid reprocessing the same event over and over to obtain the variables |
// needed to avoid reprocessing the same event over and over to obtain the variables |
42 |
// |
// |
61 |
Float_t MIP;//Mev g/cm2 energia al minimo nel silicio per 0.38 mm |
Float_t MIP;//Mev g/cm2 energia al minimo nel silicio per 0.38 mm |
62 |
|
|
63 |
Float_t emin; |
Float_t emin; |
64 |
|
Float_t fzeta; |
65 |
|
Float_t fene; |
66 |
|
|
67 |
// |
// |
68 |
Bool_t usepl18x; |
Bool_t usepl18x; |
69 |
|
Bool_t newchi2; |
70 |
|
Bool_t usenewBB; |
71 |
|
Int_t ndf; |
72 |
|
|
73 |
|
|
74 |
//parametri bethe-bloch |
//parametri bethe-bloch |
75 |
Float_t pigr; |
Float_t pigr; |
76 |
Float_t Na; |
Float_t Na; |
77 |
Float_t ZA; /*Z/A per Si*/ |
Float_t ZA; /*Z/A per Si*/ |
78 |
Float_t ISi; /*MeV*/ |
Float_t ISi; /*MeV*/ |
79 |
|
Float_t IW; /*MeV*/ |
80 |
Float_t Me; /* MeV*/ |
Float_t Me; /* MeV*/ |
81 |
Float_t MassP;/*MeV*/ |
Float_t MassP;/*MeV*/ |
82 |
Float_t r2; /*ro*ro in cm */ |
Float_t r2; /*ro*ro in cm */ |
83 |
|
|
84 |
//variabili globali caricabili |
//variabili globali caricabili |
85 |
Float_t Integrale; |
Float_t Integrale; |
86 |
Float_t bestchi2[4]; |
Float_t bestchi2[5]; |
87 |
/* Float_t qtchi2; */ |
/* Float_t qtchi2; */ |
88 |
/* Float_t qtz; */ |
/* Float_t qtz; */ |
89 |
/* Float_t qtetot; */ |
/* Float_t qtetot; */ |
90 |
/* Float_t qtpskip; */ |
/* Float_t qtpskip; */ |
91 |
Float_t lpchi2; |
Float_t lpchi2; |
92 |
Float_t lpz; |
Float_t lpz; |
93 |
|
Float_t lpisotope; |
94 |
Float_t lpetot; |
Float_t lpetot; |
95 |
Float_t lppskip; |
Float_t lppskip; |
96 |
Float_t calorimetro[44][2]; |
Float_t calorimetro[44][2]; |
97 |
Float_t spessore[3]; |
Float_t spessore[4]; |
98 |
Float_t estremi[2][2];//[0][0]primo piano [0][1]energia primo |
Float_t estremi[2][2];//[0][0]primo piano [0][1]energia primo |
99 |
//[1][0]ultim piano [1][1]energia ultim |
//[1][0]ultim piano [1][1]energia ultim |
100 |
//Float_t shift[2][22];//shift dei piani |
//Float_t shift[2][22];//shift dei piani |
101 |
//Float_t piano[22]; //posizione z dei piani |
//Float_t piano[22]; //posizione z dei piani |
102 |
Float_t elem[32];//A degli elementi |
Float_t elem[32][7];//A degli elementi compresi isotopi |
103 |
// |
// |
104 |
|
|
105 |
|
Float_t dEplan[2*22];//energia rilasciata calcolata |
106 |
|
|
107 |
/* void SWAP(Float_t*,Float_t*); */ |
/* void SWAP(Float_t*,Float_t*); */ |
108 |
void LoadParam(); |
void LoadParam(); |
109 |
void conversione(); |
void conversione(); |
110 |
void BetheBloch(Float_t *,Float_t *,Float_t *,Float_t *,Float_t *,Float_t *); |
void BetheBloch(Float_t *,Float_t *,Float_t *,Float_t *,Float_t *,Float_t *, Float_t); |
111 |
void ELOSS(Float_t *,Int_t *, Float_t *, Float_t *); |
void ELOSS(Float_t *,Int_t *,Int_t *, Float_t *, Float_t *, Float_t); |
112 |
void Enetrack(Int_t *, Float_t *, Float_t *,Float_t *,Float_t * ); |
void Enetrack(Int_t *,Int_t *, Float_t *, Float_t *,Float_t *,Float_t * ); |
113 |
void chiquadro(Float_t *, Float_t *); |
void chiquadro(Float_t *, Float_t *); |
114 |
void loopze(Float_t ,Float_t ,Float_t ,Float_t, Int_t ); |
void loopze(Float_t ,Float_t ,Float_t ,Float_t, Int_t ); |
115 |
void mediatroncata(); |
void mediatroncata(); |
118 |
|
|
119 |
public: |
public: |
120 |
// |
// |
121 |
|
void SetStartZero(Bool_t yes=false){ startZero = yes;}; |
122 |
// |
// |
123 |
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_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) |
Float_t Get_lastplane(){ Process(); return estremi[1][0];}; ///< 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) */ |
/* Float_t Get_qtpskip(){ Process(); return qtpskip;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane) */ |
129 |
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_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) |
Float_t Get_lpz(){ Process(); return lpz;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane) |
131 |
|
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 |
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_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) |
Float_t Get_lppskip(){ Process(); return lppskip;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane) |
134 |
|
|
135 |
|
Float_t Integral(); |
136 |
|
Float_t LastIntegral(); |
137 |
|
|
138 |
|
void ForceZeta(Float_t zeta){ fzeta = zeta;}; |
139 |
|
void ForceEne(Float_t ener ){ fene = ener;}; |
140 |
|
void ForceProcessing(){ sntr = -1000.;}; |
141 |
|
|
142 |
void UsePlane18X(Bool_t use){usepl18x = use;}; |
void UsePlane18X(Bool_t use){usepl18x = use;}; |
143 |
|
|
144 |
|
void UseNewChi2(Bool_t use){newchi2 = use;}; |
145 |
|
void UseNewBB(Bool_t use){usenewBB = use;}; |
146 |
// |
// |
147 |
|
|
148 |
// |
// |
151 |
~CaloBragg(){ Delete(); }; |
~CaloBragg(){ Delete(); }; |
152 |
// |
// |
153 |
void Draw(); |
void Draw(); |
154 |
|
void Draw(Int_t,Int_t, Float_t); |
155 |
void SetDebug(Bool_t d){ debug=d; }; |
void SetDebug(Bool_t d){ debug=d; }; |
156 |
void UseTrack(Bool_t d){ usetrack=d; }; |
void UseTrack(Bool_t d){ usetrack=d; }; |
157 |
// |
// |
165 |
void Print(); |
void Print(); |
166 |
void Print(Option_t *option){Print();}; |
void Print(Option_t *option){Print();}; |
167 |
// |
// |
168 |
ClassDef(CaloBragg,4); |
ClassDef(CaloBragg,6); |
169 |
}; |
}; |
170 |
|
|
171 |
#endif |
#endif |