--- calo/flight/CaloBragg/inc/CaloBragg.h	2008/10/14 13:08:28	1.3
+++ calo/flight/CaloBragg/inc/CaloBragg.h	2011/06/17 11:56:41	1.9
@@ -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 */
@@ -81,27 +92,31 @@
     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_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) 
@@ -114,12 +129,25 @@
     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);
     ~CaloBragg(){ Delete(); };
     //
     void Draw();
+    void Draw(Int_t, Float_t);
     void SetDebug(Bool_t d){ debug=d; };
     void UseTrack(Bool_t d){ usetrack=d; };
     //
@@ -133,7 +161,7 @@
     void Print();
     void Print(Option_t *option){Print();};
     //
-    ClassDef(CaloBragg,2);
+    ClassDef(CaloBragg,6);
 };
 
 #endif