X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=EMCAL%2FAliEMCALClusterizer.h;h=4aa0956328ae4a78ac1d2ba122615f2551011e6c;hb=3a282863bcc86e188d6f4d7697c46efef2facbbd;hp=a547ffa337f7acfa1579f5fa43fdf6d294f1fd93;hpb=d64c959b2712e5f0872a5c693ef5d9f0e0af07ee;p=u%2Fmrichter%2FAliRoot.git diff --git a/EMCAL/AliEMCALClusterizer.h b/EMCAL/AliEMCALClusterizer.h index a547ffa337f..4aa0956328a 100644 --- a/EMCAL/AliEMCALClusterizer.h +++ b/EMCAL/AliEMCALClusterizer.h @@ -12,53 +12,110 @@ // August 2002 Yves Schutz: clone PHOS as closely as possible and intoduction // of new IO (à la PHOS) // --- ROOT system --- - -#include "TTask.h" +#include "AliLog.h" +#include "TObject.h" +class TTree; // --- Standard library --- // --- AliRoot header files --- +class AliEMCALGeometry ; +class AliEMCALCalibData ; +class AliCaloCalibPedestal ; +class AliEMCALUnfolding ; -class AliEMCALClusterizer : public TTask { +class AliEMCALClusterizer : public TObject { public: AliEMCALClusterizer() ; // default ctor - AliEMCALClusterizer(const TString alirunFileName, const TString eventFolderName = AliConfig::fgkDefaultEventFolderName) ; - virtual ~AliEMCALClusterizer() ; // dtor - - virtual Float_t GetTowerClusteringThreshold()const {Warning("GetTowerClusteringThreshold", "Not Defined") ; return 0. ; } - virtual Float_t GetTowerLocalMaxCut()const {Warning("GetTowerLocalMaxCut", "Not Defined") ; return 0. ; } - virtual Float_t GetTowerLogWeight()const {Warning("GetTowerLogWeight", "Not Defined") ; return 0. ; } - virtual Float_t GetTimeGate() const {Warning("GetTimeGate", "Not Defined") ; return 0. ; } - virtual Float_t GetPreShoClusteringThreshold()const {Warning("GetPreShoClusteringThreshold", "Not Defined") ; return 0. ; } - virtual Float_t GetPreShoLocalMaxCut()const {Warning("GetPreShoLocalMaxCut", "Not Defined") ; return 0. ; } - virtual Float_t GetPreShoLogWeight()const {Warning("GetPreShoLogWeight", "Not Defined") ; return 0. ; } - virtual const char * GetRecPointsBranch() const {Warning("GetRecPointsBranch", "Not Defined") ; return 0 ; } - virtual const Int_t GetRecPointsInRun() const {Warning("GetRecPointsInRun", "Not Defined") ; return 0 ; } - virtual const char * GetDigitsBranch() const {Warning("GetDigitsBranch", "Not Defined") ; return 0 ; } - - virtual void MakeClusters() {Warning("MakeClusters", "Not Defined") ; } - virtual void Print(Option_t * /*option*/)const {Warning("Print", "Not Defined") ; } - - virtual void SetECAClusteringThreshold(Float_t) = 0; - virtual void SetECALocalMaxCut(Float_t) = 0; - virtual void SetECALogWeight(Float_t) = 0; - virtual void SetHCAClusteringThreshold(Float_t) = 0; - virtual void SetHCALocalMaxCut(Float_t) = 0; - virtual void SetHCALogWeight(Float_t) = 0; - virtual void SetTimeGate(Float_t) = 0; - virtual void SetPREClusteringThreshold(Float_t) = 0; - virtual void SetPRELocalMaxCut(Float_t) = 0; - virtual void SetPRELogWeight(Float_t) = 0; - virtual void SetUnfolding(Bool_t) = 0; + virtual ~AliEMCALClusterizer() ; // dtorEM + 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."); } + + //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 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 void SetUnfolding(Bool_t toUnfold = kTRUE ) {fToUnfold = toUnfold ;} + + 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* + + static void SetInputCalibrated(Bool_t val); + virtual const char * Version() const {Warning("Version", "Not Defined") ; return 0 ; } protected: - TString fEventFolderName ; // event folder name - ClassDef(AliEMCALClusterizer,3) // Clusterization algorithm class + virtual void MakeClusters() = 0; + + static Bool_t fgkIsInputCalibrated; // to enable reclusterization from ESD cells + + 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 + + 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 + + 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 + + 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 +private: + AliEMCALClusterizer(const AliEMCALClusterizer &); //copy ctor + AliEMCALClusterizer & operator = (const AliEMCALClusterizer &); + + + + ClassDef(AliEMCALClusterizer,4) // Clusterization algorithm class } ; #endif // AliEMCALCLUSTERIZER_H