--- calo/flight/CaloBragg/inc/CaloBragg.h 2008/06/18 15:22:33 1.2 +++ calo/flight/CaloBragg/inc/CaloBragg.h 2011/06/15 12:28:44 1.8 @@ -35,6 +35,8 @@ 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 // @@ -59,12 +61,21 @@ Float_t MIP;//Mev g/cm2 energia al minimo nel silicio per 0.38 mm Float_t emin; - + Float_t fzeta; + + // + 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 */ @@ -72,46 +83,64 @@ //variabili globali caricabili Float_t Integrale; Float_t bestchi2[4]; - Float_t qtchi2; - Float_t qtz; - Float_t qtetot; - Float_t qtpskip; +/* Float_t qtchi2; */ +/* Float_t qtz; */ +/* Float_t qtetot; */ +/* Float_t qtpskip; */ Float_t lpchi2; Float_t lpz; Float_t lpetot; Float_t lppskip; Float_t calorimetro[44][2]; - Float_t spessore[3]; + 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];//A degli elementi // + + 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 *); - void ELOSS(Float_t *,Int_t *, Float_t *, Float_t *); + void BetheBloch(Float_t *,Float_t *,Float_t *,Float_t *,Float_t *,Float_t *, Float_t); + void ELOSS(Float_t *,Int_t *, Float_t *, Float_t *, Float_t); void Enetrack(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 ); + void loopze(Float_t ,Float_t ,Float_t ,Float_t, Int_t ); void mediatroncata(); void Zdaloop(); - + void CleanPlanes(Float_t epiano[22][2]); public: // + void SetStartZero(Bool_t yes=false){ startZero = yes;}; // - 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_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_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 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); @@ -131,7 +160,7 @@ void Print(); void Print(Option_t *option){Print();}; // - ClassDef(CaloBragg,2); + ClassDef(CaloBragg,6); }; #endif