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 class AliEMCALRecParam;
  26 #include "AliEMCALUnfolding.h"
  27
  28 class AliEMCALClusterizer : public TObject {
  29
  30 public:
  31
  32   AliEMCALClusterizer();
  33   virtual ~AliEMCALClusterizer();
  34   AliEMCALClusterizer(AliEMCALGeometry *geometry);
  35   AliEMCALClusterizer(AliEMCALGeometry *geometry, AliEMCALCalibData *calib, AliCaloCalibPedestal *pedestal);
  36
  37   virtual void    Digits2Clusters(Option_t *option) = 0;
  38
  39   virtual Float_t Calibrate(const Float_t amp, const Float_t time, const Int_t cellId);
  40   virtual void    Init();
  41   virtual void    InitParameters();
  42   virtual void    InitParameters(const AliEMCALRecParam* recParam);
  43
  44   //Get/Set reconstruction parameters
  45   virtual void  GetCalibrationParameters(void);
  46   virtual void  GetCaloCalibPedestal(void);
  47   virtual void  SetCalibrationParameters(AliEMCALCalibData *calib)   { fCalibData = calib;   }
  48   virtual void  SetCaloCalibPedestal(AliCaloCalibPedestal  *caloped) { fCaloPed   = caloped; }
  49
  50   virtual Float_t GetTimeMin()           const { return fTimeMin; }
  51   virtual Float_t GetTimeMax()           const { return fTimeMax; }
  52   virtual Float_t GetTimeCut()           const { return fTimeCut; }
  53   virtual Float_t GetECAClusteringThreshold()   const { return fECAClusteringThreshold; }
  54   virtual Float_t GetECALocalMaxCut()           const { return fECALocMaxCut;           }
  55   virtual Float_t GetECALogWeight()             const { return fECAW0;                  }
  56   virtual Float_t GetMinECut()                  const { return fMinECut;                }
  57
  58   virtual void SetTimeMin(Float_t t)                  { fTimeMin = t; }
  59   virtual void SetTimeMax(Float_t t)                  { fTimeMax = t; }
  60   virtual void SetTimeCut(Float_t t)                  { fTimeCut = t; }
  61   virtual void SetECAClusteringThreshold(Float_t th)  { fECAClusteringThreshold = th; }
  62   virtual void SetMinECut(Float_t mine)               { fMinECut = mine;     }
  63   virtual void SetECALocalMaxCut(Float_t cut)         { fECALocMaxCut = cut; }
  64   virtual void SetECALogWeight(Float_t w)             { fECAW0 = w;          }
  65
  66   //Unfolding
  67   virtual void SetUnfolding(Bool_t toUnfold = kTRUE ) { fToUnfold = toUnfold; }
  68   virtual void SetSSPars   (Int_t ipar, Double_t par) { fSSPars[ipar] = par;  }
  69   virtual void SetPar5     (Int_t ipar, Double_t par) { fPar5  [ipar] = par;  }
  70   virtual void SetPar6     (Int_t ipar, Double_t par) { fPar6  [ipar] = par;  }
  71   virtual void InitClusterUnfolding() {
  72     fClusterUnfolding=new AliEMCALUnfolding(fGeom,fECALocMaxCut,fSSPars,fPar5,fPar6);
  73   }
  74
  75   virtual void SetInput(TTree *digitsTree);
  76   virtual void SetOutput(TTree *clustersTree);
  77
  78   virtual void Print(Option_t *option) const ;
  79   virtual void PrintRecPoints(Option_t *option);
  80   virtual void PrintRecoInfo();
  81
  82   void         SetInputCalibrated(Bool_t val);
  83   void         SetJustClusters(Bool_t val);
  84
  85   virtual const char *Version() const { Warning("Version", "Not Defined"); return 0; }
  86
  87   // add for clusterizing task
  88   virtual void              SetDigitsArr(TClonesArray *arr) { fDigitsArr = arr;  }
  89   virtual const TObjArray  *GetRecPoints() const            { return fRecPoints; }
  90
  91 protected:
  92
  93   virtual void MakeClusters() = 0;
  94
  95   Bool_t fIsInputCalibrated; // to enable reclusterization from ESD cells
  96   Bool_t fJustClusters;      // false for standard reco
  97   TClonesArray *fDigitsArr;  // array with EMCAL digits
  98   TTree        *fTreeR;      // tree with output clusters
  99   TObjArray    *fRecPoints;  // array with EMCAL clusters
 100
 101   AliEMCALGeometry     *fGeom;      //!pointer to geometry for utilities
 102   AliEMCALCalibData    *fCalibData; //!calibration database if aval
 103   AliCaloCalibPedestal *fCaloPed;   //!tower status map if aval
 104
 105   Float_t fADCchannelECA;           // width of one ADC channel for EC section (GeV)
 106   Float_t fADCpedestalECA;          // pedestal of ADC for EC section (GeV)
 107
 108   Float_t fTimeMin;                 // minimum time of physical signal in a cell/digit
 109   Float_t fTimeMax;                 // maximum time of physical signal in a cell/digit
 110   Float_t fTimeCut;                 // maximum time difference between the digits inside EMC cluster
 111
 112   Bool_t  fDefaultInit;             //!says if the task was created by defaut ctor (only parameters are initialized)
 113   Bool_t  fToUnfold;                // says if unfolding should be performed
 114   Int_t   fNumberOfECAClusters;     // number of clusters found in EC section
 115
 116   Float_t fECAClusteringThreshold;  // minimum energy to seed a EC digit in a cluster
 117   Float_t fECALocMaxCut;            // minimum energy difference to distinguish local maxima in a cluster
 118   Float_t fECAW0;                   // logarithmic weight for the cluster center of gravity calculation
 119   Float_t fMinECut;                 // minimum energy for a digit to be a member of a cluster
 120
 121   AliEMCALUnfolding *fClusterUnfolding; //!pointer to unfolding object
 122   Double_t fSSPars[8];                  // shower shape parameters
 123   Double_t fPar5[3];                    // shower shape parameter 5
 124   Double_t fPar6[3];                    // shower shape parameter 6
 125
 126  private:
 127   AliEMCALClusterizer(const AliEMCALClusterizer &);
 128   AliEMCALClusterizer & operator = (const AliEMCALClusterizer &);
 129
 130   ClassDef(AliEMCALClusterizer,6)  // Clusterization algorithm class
 131 };
 132 #endif // AliEMCALCLUSTERIZER_H