EMCAL/AliEMCALClusterizer.h

   1 #ifndef ALIEMCALCLUSTERIZER_H
   2 #define ALIEMCALCLUSTERIZER_H
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice                               */
   5
   6 /* $Id$ */
   7
   8 //_________________________________________________________________________
   9 //  Base class for the clusterization algorithm (pure abstract)
  10 //*-- Author: Yves Schutz (SUBATECH) & Dmitri Peressounko (SUBATECH & Kurchatov Institute)
  11 // Modif:
  12 //  August 2002 Yves Schutz: clone PHOS as closely as possible and intoduction
  13 //                           of new  IO (à la PHOS)
  14 // --- ROOT system ---
  15 #include "AliLog.h"
  16 #include "TObject.h"
  17 class TTree;
  18
  19 // --- Standard library ---
  20
  21 // --- AliRoot header files ---
  22 class AliEMCALGeometry ;
  23 class AliEMCALCalibData ;
  24 class AliCaloCalibPedestal ;
  25
  26 class AliEMCALClusterizer : public TObject {
  27
  28 public:
  29
  30   AliEMCALClusterizer() ;        // default ctor
  31   virtual ~AliEMCALClusterizer() ; // dtorEM
  32   AliEMCALClusterizer(AliEMCALGeometry* geometry);
  33   AliEMCALClusterizer(AliEMCALGeometry* geometry, AliEMCALCalibData * calib, AliCaloCalibPedestal * pedestal);
  34
  35   virtual void    Digits2Clusters(Option_t *option) = 0;
  36
  37   virtual Float_t Calibrate(const Float_t amp, const Float_t time, const Int_t cellId) ;  // Tranforms Amp to energy
  38   virtual void    Init() ;
  39   virtual void    InitParameters() ; //{ AliInfo("Overload this method."); }
  40
  41   //Get/Set reconstruction parameters
  42   virtual void  GetCalibrationParameters(void) ;
  43   virtual void  GetCaloCalibPedestal(void) ;
  44   virtual void  SetCalibrationParameters(AliEMCALCalibData * calib)   { fCalibData = calib ; }
  45   virtual void  SetCaloCalibPedestal(AliCaloCalibPedestal  * caloped) { fCaloPed   = caloped ; }
  46
  47   virtual Float_t GetTimeMin()           const { return fTimeMin ; }
  48   virtual Float_t GetTimeMax()           const { return fTimeMax ; }
  49   virtual Float_t GetTimeCut()           const { return fTimeCut ; }
  50   virtual void    GetNumberOfClustersFound(int numb )const { numb = fNumberOfECAClusters ;}
  51   virtual Float_t GetECAClusteringThreshold()        const { return fECAClusteringThreshold;}
  52   virtual Float_t GetECALocalMaxCut()                const { return fECALocMaxCut;}
  53   virtual Float_t GetECALogWeight()                  const { return fECAW0;}
  54   virtual Float_t GetMinECut()                       const { return fMinECut;}
  55
  56   virtual void SetTimeMin(Float_t t)                                            { fTimeMin = t ;}
  57   virtual void SetTimeMax(Float_t t)                                            { fTimeMax = t ;}
  58   virtual void SetTimeCut(Float_t t)                                            { fTimeCut = t ;}
  59   virtual void SetECAClusteringThreshold(Float_t th)  { fECAClusteringThreshold = th ; }
  60   virtual void SetMinECut(Float_t mine)               { fMinECut = mine; }
  61   virtual void SetECALocalMaxCut(Float_t cut)         { fECALocMaxCut = cut ; }
  62   virtual void SetECALogWeight(Float_t w)             { fECAW0 = w ; }
  63   virtual void SetUnfolding(Bool_t toUnfold = kTRUE ) {fToUnfold = toUnfold ;}
  64
  65   virtual void SetInput(TTree *digitsTree);
  66   virtual void SetOutput(TTree *clustersTree);
  67
  68   virtual void Print(Option_t * option)const ;
  69   virtual void PrintRecPoints(Option_t * option);
  70   virtual void PrintRecoInfo();                        //*MENU*
  71
  72   virtual const char * Version() const {Warning("Version", "Not Defined") ; return 0 ; }
  73
  74 protected:
  75
  76   virtual void MakeClusters() = 0;
  77
  78   TClonesArray *fDigitsArr; // Array with EMCAL digits
  79   TTree *fTreeR;            // Tree with output clusters
  80   TObjArray    *fRecPoints; // Array with EMCAL clusters
  81
  82   AliEMCALGeometry     * fGeom;       //! pointer to geometry for utilities
  83   AliEMCALCalibData    * fCalibData ; //! Calibration database if aval
  84   AliCaloCalibPedestal * fCaloPed   ; //! Tower status map if aval
  85
  86   Float_t fADCchannelECA ;           // width of one ADC channel for EC section (GeV)
  87   Float_t fADCpedestalECA ;          // pedestal of ADC for EC section (GeV)
  88
  89   Float_t fTimeMin ;                 // Minimum time of physical signal in a cell/digit
  90   Float_t fTimeMax ;                 // Maximum time of physical signal in a cell/digit
  91   Float_t fTimeCut ;                 // Maximum time difference between the digits inside EMC cluster
  92
  93   Bool_t  fDefaultInit;              //! Says if the task was created by defaut ctor (only parameters are initialized)
  94   Bool_t  fToUnfold ;                // To perform unfolding
  95   Int_t   fNumberOfECAClusters ;     // number of clusters found in EC section
  96
  97   Float_t fECAClusteringThreshold ;  // minimum energy to seed a EC digit in a cluster
  98   Float_t fECALocMaxCut ;            // minimum energy difference to distinguish local maxima in a cluster
  99   Float_t fECAW0 ;                   // logarithmic weight for the cluster center of gravity calculation
 100   Float_t fMinECut;                  // Minimum energy for a digit to be a member of a cluster
 101
 102 private:
 103   AliEMCALClusterizer(const AliEMCALClusterizer &); //copy ctor
 104   AliEMCALClusterizer & operator = (const AliEMCALClusterizer &);
 105
 106
 107
 108   ClassDef(AliEMCALClusterizer,2)  // Clusterization algorithm class
 109 } ;
 110
 111 #endif // AliEMCALCLUSTERIZER_H