#ifndef ALIEMCALCLUSTERIZER_H
#define ALIEMCALCLUSTERIZER_H
#include <TObject.h>
#include <TClonesArray.h>
class TTree;
#include "AliLog.h"
class AliEMCALGeometry;
class AliEMCALCalibData;
class AliCaloCalibPedestal;
class AliEMCALRecParam;
#include "AliEMCALUnfolding.h"
class AliEMCALClusterizer : public TObject {
public:
AliEMCALClusterizer();
AliEMCALClusterizer(AliEMCALGeometry *geometry);
AliEMCALClusterizer(AliEMCALGeometry *geometry, AliEMCALCalibData *calib, AliCaloCalibPedestal *pedestal);
virtual ~AliEMCALClusterizer();
virtual void DeleteDigits();
virtual void DeleteRecPoints();
virtual void Digits2Clusters(Option_t *option) = 0;
virtual void Calibrate(Float_t & amp, Float_t & time, const Int_t cellId);
virtual void Init();
virtual void InitParameters();
virtual void InitParameters(const AliEMCALRecParam* recParam);
virtual void Print (Option_t *option) const ;
virtual void PrintRecPoints(Option_t *option);
virtual void PrintRecoInfo();
virtual const char *Version() const { Warning("Version", "Not Defined");
return 0 ; }
virtual void SetInput (TTree *digitsTree );
virtual void SetOutput(TTree *clustersTree);
virtual void GetCalibrationParameters(void);
virtual void GetCaloCalibPedestal(void);
virtual void SetCalibrationParameters(AliEMCALCalibData *calib) { fCalibData = calib; }
virtual void SetCaloCalibPedestal(AliCaloCalibPedestal *caped) { fCaloPed = caped; }
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 Bool_t GetRejectBelowThreshold() const { return fRejectBelowThreshold; }
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 void SetRejectBelowThreshold(Bool_t reject) { fRejectBelowThreshold = reject; }
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);
fClusterUnfolding->SetThreshold(fMinECut);
fClusterUnfolding->SetRejectBelowThreshold(fRejectBelowThreshold); }
virtual void SetNRowDiff(Int_t ) { ; }
virtual void SetNColDiff(Int_t ) { ; }
virtual void SetEnergyGrad(Bool_t ) { ; }
virtual Int_t GetNRowDiff() const { return -1 ; }
virtual Int_t GetNColDiff() const { return -1 ; }
virtual Bool_t GetEnergyGrad() const { return -1 ; }
virtual void SetDigitsArr(TClonesArray *arr) { fDigitsArr = arr ; }
virtual TClonesArray *GetDigits() { if (!fDigitsArr)
fDigitsArr = new TClonesArray("AliEMCALDigit",12000);
return fDigitsArr ; }
virtual const TObjArray *GetRecPoints() const { return fRecPoints ; }
void SetInputCalibrated(Bool_t val);
void SetJustClusters (Bool_t val);
protected:
virtual void MakeClusters() = 0;
Bool_t fIsInputCalibrated;
Bool_t fJustClusters;
TClonesArray *fDigitsArr;
TTree *fTreeR;
TObjArray *fRecPoints;
AliEMCALGeometry *fGeom;
AliEMCALCalibData *fCalibData;
AliCaloCalibPedestal *fCaloPed;
Float_t fADCchannelECA;
Float_t fADCpedestalECA;
Float_t fTimeECA;
Float_t fTimeMin;
Float_t fTimeMax;
Float_t fTimeCut;
Bool_t fDefaultInit;
Bool_t fToUnfold;
Int_t fNumberOfECAClusters;
Float_t fECAClusteringThreshold;
Float_t fECALocMaxCut;
Float_t fECAW0;
Float_t fMinECut;
Bool_t fRejectBelowThreshold;
AliEMCALUnfolding *fClusterUnfolding;
Double_t fSSPars[8];
Double_t fPar5[3];
Double_t fPar6[3];
private:
AliEMCALClusterizer( const AliEMCALClusterizer &);
AliEMCALClusterizer & operator = (const AliEMCALClusterizer &);
ClassDef(AliEMCALClusterizer,8)
};
#endif // AliEMCALCLUSTERIZER_H
AliEMCALClusterizer.h:100 AliEMCALClusterizer.h:101 AliEMCALClusterizer.h:102 AliEMCALClusterizer.h:103 AliEMCALClusterizer.h:104 AliEMCALClusterizer.h:105 AliEMCALClusterizer.h:106 AliEMCALClusterizer.h:107 AliEMCALClusterizer.h:108 AliEMCALClusterizer.h:109 AliEMCALClusterizer.h:110 AliEMCALClusterizer.h:111 AliEMCALClusterizer.h:112 AliEMCALClusterizer.h:113 AliEMCALClusterizer.h:114 AliEMCALClusterizer.h:115 AliEMCALClusterizer.h:116 AliEMCALClusterizer.h:117 AliEMCALClusterizer.h:118 AliEMCALClusterizer.h:119 AliEMCALClusterizer.h:120 AliEMCALClusterizer.h:121 AliEMCALClusterizer.h:122 AliEMCALClusterizer.h:123 AliEMCALClusterizer.h:124 AliEMCALClusterizer.h:125 AliEMCALClusterizer.h:126 AliEMCALClusterizer.h:127 AliEMCALClusterizer.h:128 AliEMCALClusterizer.h:129 AliEMCALClusterizer.h:130 AliEMCALClusterizer.h:131 AliEMCALClusterizer.h:132 AliEMCALClusterizer.h:133 AliEMCALClusterizer.h:134 AliEMCALClusterizer.h:135 AliEMCALClusterizer.h:136 AliEMCALClusterizer.h:137 AliEMCALClusterizer.h:138 AliEMCALClusterizer.h:139 AliEMCALClusterizer.h:140 AliEMCALClusterizer.h:141 AliEMCALClusterizer.h:142 AliEMCALClusterizer.h:143 AliEMCALClusterizer.h:144 AliEMCALClusterizer.h:145 AliEMCALClusterizer.h:146 AliEMCALClusterizer.h:147 AliEMCALClusterizer.h:148 AliEMCALClusterizer.h:149 AliEMCALClusterizer.h:150 AliEMCALClusterizer.h:151 AliEMCALClusterizer.h:152 AliEMCALClusterizer.h:153 AliEMCALClusterizer.h:154 AliEMCALClusterizer.h:155 AliEMCALClusterizer.h:156 AliEMCALClusterizer.h:157 AliEMCALClusterizer.h:158 AliEMCALClusterizer.h:159 AliEMCALClusterizer.h:160 AliEMCALClusterizer.h:161