--- calo/flight/CaloNuclei/inc/CaloNuclei.h	2007/04/03 15:55:16	1.1
+++ calo/flight/CaloNuclei/inc/CaloNuclei.h	2007/04/04 11:04:55	1.2
@@ -41,6 +41,7 @@
     //
     Int_t interplane; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane) 
     Int_t N; ///< Number of dE/dx measurements to be used to calculate qpremeanN, default N = 5
+    Int_t R; ///< Number of strip to be used around the trajectory to calculate qpremeanN, default R = 3
     Float_t preq; ///< Energy release (MIP) up to the interaction plane (included)
     Float_t postq; ///< Energy release (MIP) from the interaction plane (excluded) up to the last plane
     Float_t dedx1; ///< Energy release (MIP) along the track on the first Silicon detector (Y EVEN).
@@ -57,16 +58,18 @@
     //
     Int_t GetInterplane(){ Process(); return interplane;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane) 
     Int_t GetN(){ return N;}; ///< Number of available dE/dx measurements before interaction or exit from the calo (interaction plane) 
+    Int_t GetR(){ return R;}; ///< Number of strip to be used around the trajectory to calculate qpremeanN, default R = 3 
     Float_t GetPreq(){ Process(); return preq;}; ///< Energy release (MIP) up to the interaction plane (included)
     Float_t GetPostq(){ Process(); return postq;}; ///< Energy release (MIP) from the interaction plane (excluded) up to the last plane
     Float_t GetDedx1(){ Process(); return dedx1;}; ///< Energy release (MIP) along the track on the first Silicon detector (Y EVEN).
     Float_t GetDedx3(){ Process(); return dedx3;}; ///< Energy release (MIP) along the track on the first three Silicon detectors (Y EVEN, X EVEN, Y ODD).
-    Float_t GetQpremean(){ Process(); return qpremean;}; ///< Truncated mean (MIP) along the track up to the interaction plane preq.
-    Float_t GetQpremeanN(){ Process(); return qpremeanN;}; ///< Truncated mean (MIP) along the track up to the interaction plane preq.
+    Float_t GetQpremean(){ Process(); return qpremean;}; ///< Truncated mean (MIP) along the track up to the interaction plane preq using three points
+    Float_t GetQpremeanN(){ Process(); return qpremeanN;}; ///< Truncated mean (MIP) along the track up to the interaction plane preq using N points
     Float_t GetEthr(){ Process(); return ethr;}; ///< Threshold (MIP) needed to find the interaction plane
     Bool_t IsMulthit(){ Process(); return multhit;}; ///< True if the interaction plane has been determined by multiple hit counting
     //
     void SetN(Int_t n){ N=n;};
+    void SetR(Int_t r){ R=r;};
     //
     CaloNuclei();
     CaloNuclei(PamLevel2 *L2);