]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - EMCAL/AliEMCALClusterizer.h
Allow running without tree by using
[u/mrichter/AliRoot.git] / EMCAL / AliEMCALClusterizer.h
index 2e2e52be1a64d95b88718c5e7dff0e166bd61455..77866d78695d28b8f9d0b52be89c019cd21c2978 100644 (file)
@@ -19,54 +19,55 @@ class TTree;
 // --- Standard library ---
 
 // --- AliRoot header files ---
-class AliEMCALGeometry ;
-class AliEMCALCalibData ;
-class AliCaloCalibPedestal ;
+class AliEMCALGeometry;
+class AliEMCALCalibData;
+class AliCaloCalibPedestal;
+class AliEMCALRecParam;
 #include "AliEMCALUnfolding.h"
 
 class AliEMCALClusterizer : public TObject {
 
 public:
 
-  AliEMCALClusterizer() ;        // default ctor
-  virtual ~AliEMCALClusterizer() ; // dtorEM
-  AliEMCALClusterizer(AliEMCALGeometrygeometry);
-  AliEMCALClusterizer(AliEMCALGeometry* geometry, AliEMCALCalibData * calib, AliCaloCalibPedestal * pedestal);
+  AliEMCALClusterizer();
+  virtual ~AliEMCALClusterizer();
+  AliEMCALClusterizer(AliEMCALGeometry *geometry);
+  AliEMCALClusterizer(AliEMCALGeometry *geometry, AliEMCALCalibData *calib, AliCaloCalibPedestal *pedestal);
 
   virtual void    Digits2Clusters(Option_t *option) = 0;
 
-  virtual Float_t Calibrate(const Float_t amp, const Float_t time, const Int_t cellId) ;  // Tranforms Amp to energy 
-  virtual void    Init() ;
-  virtual void    InitParameters() ; //{ AliInfo("Overload this method."); }
+  virtual Float_t Calibrate(const Float_t amp, const Float_t time, const Int_t cellId);
+  virtual void    Init();
+  virtual void    InitParameters();
+  virtual void    InitParameters(const AliEMCALRecParam* recParam);
 
   //Get/Set reconstruction parameters
-  virtual void  GetCalibrationParameters(void) ;
-  virtual void  GetCaloCalibPedestal(void) ;
-  virtual void  SetCalibrationParameters(AliEMCALCalibData * calib)   { fCalibData = calib ; }
-  virtual void  SetCaloCalibPedestal(AliCaloCalibPedestal  * caloped) { fCaloPed   = caloped ; }
+  virtual void  GetCalibrationParameters(void);
+  virtual void  GetCaloCalibPedestal(void);
+  virtual void  SetCalibrationParameters(AliEMCALCalibData *calib)   { fCalibData = calib;   }
+  virtual void  SetCaloCalibPedestal(AliCaloCalibPedestal  *caloped) { fCaloPed   = caloped; }
   
-  virtual Float_t GetTimeMin()           const { return fTimeMin ; }
-  virtual Float_t GetTimeMax()           const { return fTimeMax ; }
-  virtual Float_t GetTimeCut()           const { return fTimeCut ; }
-  //virtual void    GetNumberOfClustersFound(int numb )const { numb = fNumberOfECAClusters ;} 
-  virtual Float_t GetECAClusteringThreshold()        const { return fECAClusteringThreshold;}  
-  virtual Float_t GetECALocalMaxCut()                const { return fECALocMaxCut;} 
-  virtual Float_t GetECALogWeight()                  const { return fECAW0;}
-  virtual Float_t GetMinECut()                       const { return fMinECut;}
-
-  virtual void SetTimeMin(Float_t t)                                           { fTimeMin = t ;}
-  virtual void SetTimeMax(Float_t t)                                           { fTimeMax = t ;}
-  virtual void SetTimeCut(Float_t t)                                           { fTimeCut = t ;}
-  virtual void SetECAClusteringThreshold(Float_t th)  { fECAClusteringThreshold = th ; }
-  virtual void SetMinECut(Float_t mine)               { fMinECut = mine; }
-  virtual void SetECALocalMaxCut(Float_t cut)         { fECALocMaxCut = cut ; }
-  virtual void SetECALogWeight(Float_t w)             { fECAW0 = w ; }
+  virtual Float_t GetTimeMin()           const { return fTimeMin; }
+  virtual Float_t GetTimeMax()           const { return fTimeMax; }
+  virtual Float_t GetTimeCut()           const { return fTimeCut; }
+  virtual Float_t GetECAClusteringThreshold()   const { return fECAClusteringThreshold; }  
+  virtual Float_t GetECALocalMaxCut()           const { return fECALocMaxCut;           } 
+  virtual Float_t GetECALogWeight()             const { return fECAW0;                  }
+  virtual Float_t GetMinECut()                  const { return fMinECut;                }
+
+  virtual void SetTimeMin(Float_t t)                 { fTimeMin = t; }
+  virtual void SetTimeMax(Float_t t)                 { fTimeMax = t; }
+  virtual void SetTimeCut(Float_t t)                 { fTimeCut = t; }
+  virtual void SetECAClusteringThreshold(Float_t th)  { fECAClusteringThreshold = th; }
+  virtual void SetMinECut(Float_t mine)               { fMinECut = mine;     }
+  virtual void SetECALocalMaxCut(Float_t cut)         { fECALocMaxCut = cut; }
+  virtual void SetECALogWeight(Float_t w)             { fECAW0 = w;          }
   
   //Unfolding
-  virtual void SetUnfolding(Bool_t toUnfold = kTRUE ) { fToUnfold = toUnfold ;}  
-  virtual void SetSSPars   (Int_t ipar, Double_t par) { fSSPars[ipar] = par  ;}
-  virtual void SetPar5     (Int_t ipar, Double_t par) { fPar5  [ipar] = par  ;}
-  virtual void SetPar6     (Int_t ipar, Double_t par) { fPar6  [ipar] = par  ;}
+  virtual void SetUnfolding(Bool_t toUnfold = kTRUE ) { fToUnfold = toUnfold}  
+  virtual void SetSSPars   (Int_t ipar, Double_t par) { fSSPars[ipar] = par;  }
+  virtual void SetPar5     (Int_t ipar, Double_t par) { fPar5  [ipar] = par;  }
+  virtual void SetPar6     (Int_t ipar, Double_t par) { fPar6  [ipar] = par;  }
   virtual void InitClusterUnfolding() {
     fClusterUnfolding=new AliEMCALUnfolding(fGeom,fECALocMaxCut,fSSPars,fPar5,fPar6);
   }
@@ -74,56 +75,57 @@ public:
   virtual void SetInput(TTree *digitsTree);
   virtual void SetOutput(TTree *clustersTree);
   
-  virtual void Print(Option_t * option)const ;
-  virtual void PrintRecPoints(Option_t * option);
-  virtual void PrintRecoInfo();                        //*MENU*
+  virtual void Print(Option_t *option) const ;
+  virtual void PrintRecPoints(Option_t *option);
+  virtual void PrintRecoInfo();
 
-  static void     SetInputCalibrated(Bool_t val);
+  void    SetInputCalibrated(Bool_t val);
   
-  virtual const char * Version() const {Warning("Version", "Not Defined") ; return 0 ; } 
+  virtual const char *Version() const { Warning("Version", "Not Defined"); return 0; } 
+
+  // add for clusterizing task
+  virtual void              SetDigitsArr(TClonesArray *arr) { fDigitsArr = arr;  }
+  virtual const TObjArray  *GetRecPoints() const            { return fRecPoints; }
 
 protected:
 
   virtual void MakeClusters() = 0;
   
-  static Bool_t fgkIsInputCalibrated; // to enable reclusterization from ESD cells
+  Bool_t fIsInputCalibrated; // to enable reclusterization from ESD cells
   
-  TClonesArray *fDigitsArr; // Array with EMCAL digits
-  TTree *fTreeR;            // Tree with output clusters
-  TObjArray    *fRecPoints; // Array with EMCAL clusters
+  TClonesArray *fDigitsArr;  // array with EMCAL digits
+  TTree        *fTreeR;      // tree with output clusters
+  TObjArray    *fRecPoints;  // array with EMCAL clusters
   
-  AliEMCALGeometry     * fGeom;       //! pointer to geometry for utilities
-  AliEMCALCalibData    * fCalibData ; //! Calibration database if aval
-  AliCaloCalibPedestal * fCaloPed   ; //! Tower status map if aval
+  AliEMCALGeometry     *fGeom;      //!pointer to geometry for utilities
+  AliEMCALCalibData    *fCalibData; //!calibration database if aval
+  AliCaloCalibPedestal *fCaloPed;   //!tower status map if aval
   
-  Float_t fADCchannelECA ;           // width of one ADC channel for EC section (GeV)
-  Float_t fADCpedestalECA ;          // pedestal of ADC for EC section (GeV) 
+  Float_t fADCchannelECA;           // width of one ADC channel for EC section (GeV)
+  Float_t fADCpedestalECA;          // pedestal of ADC for EC section (GeV) 
 
-  Float_t fTimeMin ;                 // Minimum time of physical signal in a cell/digit
-  Float_t fTimeMax ;                 // Maximum time of physical signal in a cell/digit
-  Float_t fTimeCut ;                 // Maximum time difference between the digits inside EMC cluster
+  Float_t fTimeMin;                 // minimum time of physical signal in a cell/digit
+  Float_t fTimeMax;                 // maximum time of physical signal in a cell/digit
+  Float_t fTimeCut;                 // maximum time difference between the digits inside EMC cluster
 
-  Bool_t  fDefaultInit;              //! Says if the task was created by defaut ctor (only parameters are initialized)
-  Bool_t  fToUnfold ;                // To perform unfolding 
-  Int_t   fNumberOfECAClusters ;     // number of clusters found in EC section
+  Bool_t  fDefaultInit;             //!says if the task was created by defaut ctor (only parameters are initialized)
+  Bool_t  fToUnfold;                // says if unfolding should be performed 
+  Int_t   fNumberOfECAClusters;     // number of clusters found in EC section
   
-  Float_t fECAClusteringThreshold ;  // minimum energy to seed a EC digit in a cluster
-  Float_t fECALocMaxCut ;            // minimum energy difference to distinguish local maxima in a cluster
-  Float_t fECAW0 ;                   // logarithmic weight for the cluster center of gravity calculation
-  Float_t fMinECut;                  // Minimum energy for a digit to be a member of a cluster
+  Float_t fECAClusteringThreshold;  // minimum energy to seed a EC digit in a cluster
+  Float_t fECALocMaxCut;            // minimum energy difference to distinguish local maxima in a cluster
+  Float_t fECAW0;                   // logarithmic weight for the cluster center of gravity calculation
+  Float_t fMinECut;                 // minimum energy for a digit to be a member of a cluster
   
-  AliEMCALUnfolding * fClusterUnfolding ; //! pointer to unfolding object
-  Double_t fSSPars[8];// Shower shape parameters 
-  Double_t fPar5[3];  // Shower shape parameter 5
-  Double_t fPar6[3];  // Shower shape parameter 6
+  AliEMCALUnfolding *fClusterUnfolding; //!pointer to unfolding object
+  Double_t fSSPars[8];                  // shower shape parameters 
+  Double_t fPar5[3];                    // shower shape parameter 5
+  Double_t fPar6[3];                    // shower shape parameter 6
 
-private:
-  AliEMCALClusterizer(const AliEMCALClusterizer &); //copy ctor
+ private:
+  AliEMCALClusterizer(const AliEMCALClusterizer &);
   AliEMCALClusterizer & operator = (const AliEMCALClusterizer &);
   
-  
-  
-  ClassDef(AliEMCALClusterizer,4)  // Clusterization algorithm class 
-} ;
-
+  ClassDef(AliEMCALClusterizer,5)  // Clusterization algorithm class 
+};
 #endif // AliEMCALCLUSTERIZER_H