Class AliHMPIDRecoParamV1 introduced (29-11-2009) as a copy of AliHMPIDRecoParam...

[u/mrichter/AliRoot.git] / HMPID / AliHMPIDCluster.h

diff --git a/HMPID/AliHMPIDCluster.h b/HMPID/AliHMPIDCluster.h
index 9d1f644ca4c1db4815d8ef800747a9e2f9d358b5..e5ee8951160b20628f18986c3ecb37ae617bc7f5 100644 (file)
--- a/HMPID/AliHMPIDCluster.h
+++ b/HMPID/AliHMPIDCluster.h
@@ -19,13 +19,11 @@ class TClonesArray;        //Solve()
 class AliHMPIDCluster :public AliCluster3D
 {
 public:
-  enum EClusterStatus {kFrm,kCoG,kLo1,kUnf,kMax,kNot,kEdg,kSi1,kNoLoc,kAbn,kEmp=-1};      //status flags    
+  enum EClusterStatus {kFrm,kCoG,kLo1,kUnf,kMax,kNot,kEdg,kSi1,kNoLoc,kAbn,kBig,kEmp=-1};      //status flags    
       AliHMPIDCluster():AliCluster3D(),
                           fCh(-1),fSi(-1),fSt(kEmp),fBox(-1),fNlocMax(-1),fMaxQpad(-1),fMaxQ(-1),fQRaw(0),
-                          fQ(0),fErrQ(-1),fXX(0),fErrX(-1),fYY(0),fErrY(-1),fChi2(-1),fDigs(0) 
+                          fQ(0),fErrQ(-1),fXX(0),fErrX(-1),fYY(0),fErrY(-1),fChi2(-1),fDigs(0),fParam(AliHMPIDParam::InstanceNoGeo())
       {
-        AliHMPIDParam *pParam=AliHMPIDParam::Instance();
-        fParam = pParam;
       }//ctor
   
 
@@ -33,14 +31,14 @@ public:
                         fCh(c.fCh),fSi(c.fSi),fSt(c.fSt),fBox(c.fBox),fNlocMax(c.fNlocMax),fMaxQpad(c.fMaxQpad),fMaxQ(c.fMaxQ),fQRaw(c.fQRaw),
                         fQ (c.fQ ),fErrQ(c.fErrQ),
                         fXX (c.fXX ),fErrX(c.fErrX),
-                        fYY (c.fYY ),fErrY(c.fErrY),fChi2(c.fChi2),fDigs(0)                  {}//copy ctor
+                        fYY (c.fYY ),fErrY(c.fErrY),fChi2(c.fChi2),fDigs(0),fParam(c.fParam)                  {}//copy ctor
    virtual ~AliHMPIDCluster();//dtor   {if(fDigs) delete fDigs; fDigs=0;}
 //framework part                   
          void           Draw   (Option_t *opt=""                                  );                       //overloaded TObject::Print() to draw cluster in current canvas
          void           Print  (Option_t *opt=""                                  )const;                  //overloaded TObject::Print() to print cluster info
   static void           FitFunc(Int_t &iNpars, Double_t* deriv, Double_t &chi2, Double_t *par, Int_t iflag);//fit function to be used by MINUIT
 //private part  
-         Int_t          Box      (                                         )const{return fBox;                                  }  //Dimension of the cluster
+         Int_t          Box      (                                         )const{return fBox;                                   } //Dimension of the cluster
          void           CoG      (                                         );                                                      //calculates center of gravity
          void           CorrSin  (                                         );                                                      //sinoidal correction   
          Int_t          Ch       (                                         )const{return fCh;                                    } //chamber number
@@ -50,19 +48,27 @@ public:
   inline void           Reset    (                                         );                                                      //cleans the cluster
   void           SetClusterParams(Double_t xL,Double_t yL,Int_t iCh  );                                                            //Set AliCluster3D part
          Int_t          Size     (                                         )const{return fSi;                                    } //returns number of pads in formed cluster 
-         Int_t          Solve    (TClonesArray *pCluLst,Bool_t isUnfold    );                                                      //solve cluster: MINUIT fit or CoG
+         Int_t          Solve    (TClonesArray *pCluLst,Int_t *pSigmaCut, Bool_t isUnfold);                                        //solve cluster: MINUIT fit or CoG
          Int_t          Status   (                                         ) const{return fSt;}                                    //Status of cluster                                  
          Double_t       QRaw     (                                         )const{return fQRaw;                                  } //raw cluster charge in QDC channels 
          Double_t       Q        (                                         )const{return fQ;                                     } //given cluster charge in QDC channels 
          Double_t       Qe       (                                         )const{return fErrQ;                                  } //Error in cluster charge in QDC channels 
-         Double_t       X        (                                         )const{return fXX;                                     } //cluster x position in LRS
+         Double_t       X        (                                         )const{return fXX;                                    } //cluster x position in LRS
          Double_t       Xe       (                                         )const{return fErrX;                                  } //cluster charge in QDC channels 
-         Double_t       Y        (                                         )const{return fYY;                                     } //cluster y position in LRS 
+         Double_t       Y        (                                         )const{return fYY;                                    } //cluster y position in LRS 
          Double_t       Ye       (                                         )const{return fErrY;                                  } //cluster charge in QDC channels 
          Double_t       Chi2     (                                         )const{return fChi2;                                  } //chi2 of the fit
-         void           DoCorrSin(Bool_t doCorrSin                         ){fgDoCorrSin=doCorrSin;}                                // Set sinoidal correction
-         void           SetX     (Double_t x                               ){fXX=x;}                                                // Setter
-         void           SetY     (Double_t y                               ){fYY=y;}                                                // Setter
+         void           DoCorrSin(Bool_t doCorrSin                         ){fgDoCorrSin=doCorrSin;}                               // Set sinoidal correction
+         void           SetX     (Double_t x                               ){fXX=x;}                                               // Setter
+         void           SetY     (Double_t y                               ){fYY=y;}                                               // Setter
+         void           SetQ     (Double_t q                               ){fQ=q;if(fQ>4095)fQ=4095;}                             // Setter         
+         void           SetQRaw  (Double_t qRaw                            ){fQRaw=qRaw;if(fQRaw>4095)fQRaw=4095;}                 // Setter         
+         void           SetSize  (Int_t size                               ){fSi=size;}                                            // Setter
+         void           SetCh    (Int_t chamber                            ){fCh=chamber;}                                         // Setter
+         void           SetChi2  (Double_t chi2                            ){fChi2=chi2;}                                          // Setter
+         void           SetStatus(Int_t status                             ){fSt=status;}                                              // Setter
+         void           FindClusterSize(Int_t i,Int_t *pSigmaCut);                                                                 //Find the clusterSize of deconvoluted clusters 
+ virtual void          Clear(const Option_t*) { delete [] fDigs; fDigs=0; delete [] fParam; fParam=0; }
          
 protected:
   Int_t         fCh;          //chamber number
@@ -82,7 +88,7 @@ protected:
   Double_t      fChi2;        //some estimator of the fit quality
   TObjArray    *fDigs;        //! list of digits forming this cluster
   static  Bool_t fgDoCorrSin; //flag to switch on/off correction for Sinusoidal to cluster reco
-  AliHMPIDParam *fParam;      // Pointer to AliHMPIDParam
+  AliHMPIDParam *fParam;      //!Pointer to AliHMPIDParam
   
 private:
 /*
@@ -96,7 +102,7 @@ private:
   AliHMPIDCluster &operator=(const AliHMPIDCluster& c);   //dummy assignment operator
       
   
-  ClassDef(AliHMPIDCluster,8) //HMPID cluster class
+  ClassDef(AliHMPIDCluster,9) //HMPID cluster class
 };
 
 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++