X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=EMCAL%2FAliEMCALClusterizer.h;h=77866d78695d28b8f9d0b52be89c019cd21c2978;hb=ac084be7bc414c5a0f7dbf45412e7f8efa590c93;hp=2e2e52be1a64d95b88718c5e7dff0e166bd61455;hpb=0d0d6b983a5bd02ff0618a84fc99cde8020a6f24;p=u%2Fmrichter%2FAliRoot.git diff --git a/EMCAL/AliEMCALClusterizer.h b/EMCAL/AliEMCALClusterizer.h index 2e2e52be1a6..77866d78695 100644 --- a/EMCAL/AliEMCALClusterizer.h +++ b/EMCAL/AliEMCALClusterizer.h @@ -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(AliEMCALGeometry* geometry); - 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