X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=EMCAL%2FAliEMCALClusterizerv1.h;h=ef68008651dda54246176319fcd621c1a20b7338;hb=00ee79995be2555766b76dcac7641455e991a4dd;hp=31e841a121312beeef0afcefb65f82dc64ac3926;hpb=9e5d206744cddf451f4a32a4928872ad2f52d2cc;p=u%2Fmrichter%2FAliRoot.git diff --git a/EMCAL/AliEMCALClusterizerv1.h b/EMCAL/AliEMCALClusterizerv1.h index 31e841a1213..ef68008651d 100644 --- a/EMCAL/AliEMCALClusterizerv1.h +++ b/EMCAL/AliEMCALClusterizerv1.h @@ -24,120 +24,110 @@ // --- AliRoot header files --- #include "AliEMCALClusterizer.h" -class AliEMCALTowerRecPoint ; +class TH1F; +class AliEMCALRecPoint ; class AliEMCALDigit ; class AliEMCALDigitizer ; class AliEMCALGeometry ; - +class AliEMCALCalibData ; class AliEMCALClusterizerv1 : public AliEMCALClusterizer { public: AliEMCALClusterizerv1() ; - AliEMCALClusterizerv1(const TString alirunFileNameFile, const TString eventFolderName = AliConfig::fgkDefaultEventFolderName); + virtual ~AliEMCALClusterizerv1() ; - + virtual Int_t AreNeighbours(AliEMCALDigit * d1, AliEMCALDigit * d2)const ; // Checks if digits are in neighbour cells - virtual Float_t Calibrate(Int_t amp, Int_t where)const ; // Tranforms Amp to energy - - virtual void GetNumberOfClustersFound(int * numb )const{ numb[0] = fNumberOfPREClusters ; - numb[1] = fNumberOfECAClusters ; - numb[2] = fNumberOfHCAClusters ; } + virtual Float_t Calibrate(Int_t amp, Int_t cellId) ; // Tranforms Amp to energy - virtual Float_t GetPREClusteringThreshold()const{ return fPREClusteringThreshold; } - virtual Float_t GetECAClusteringThreshold()const{ return fECAClusteringThreshold;} - virtual Float_t GetHCAClusteringThreshold()const{ return fHCAClusteringThreshold;} - - virtual Float_t GetPRELocalMaxCut()const { return fPRELocMaxCut;} - virtual Float_t GetPREShoLogWeight()const { return fPREW0;} + 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 GetHCALocalMaxCut()const { return fHCALocMaxCut;} - virtual Float_t GetHCALogWeight()const { return fHCAW0;} + virtual Float_t GetECALogWeight()const { return fECAW0;} + virtual Float_t GetMinECut()const { return fMinECut;} - virtual Float_t GetTimeGate() const { return fTimeGate ; } - virtual const char * GetRecPointsBranch() const{ return GetName() ;} - virtual const Int_t GetRecPointsInRun() const {return fRecPointsInRun ;} + virtual Float_t GetTimeCut() const { return fTimeCut ; } - void Exec(Option_t *option); // Does the job + virtual void Digits2Clusters(Option_t *option); // Does the job virtual void Print(Option_t * option)const ; virtual void SetECAClusteringThreshold(Float_t cluth) { fECAClusteringThreshold = cluth ; } + 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 SetHCAClusteringThreshold(Float_t cluth) { fHCAClusteringThreshold = cluth ; } - virtual void SetHCALocalMaxCut(Float_t cut) { fHCALocMaxCut = cut ; } - virtual void SetHCALogWeight(Float_t w) { fHCAW0 = w ; } - virtual void SetTimeGate(Float_t gate) { fTimeGate = gate ;} - virtual void SetPREClusteringThreshold(Float_t cluth) { fPREClusteringThreshold = cluth ; } - virtual void SetPRELocalMaxCut(Float_t cut) { fPRELocMaxCut = cut ; } - virtual void SetPRELogWeight(Float_t w) { fPREW0 = w ; } + virtual void SetTimeCut(Float_t gate) { fTimeCut = gate ;} virtual void SetUnfolding(Bool_t toUnfold = kTRUE ) {fToUnfold = toUnfold ;} static Double_t ShowerShape(Double_t r) ; // Shape of EM shower used in unfolding; //class member function (not object member function) static void UnfoldingChiSquare(Int_t & nPar, Double_t * Grad, Double_t & fret, Double_t * x, Int_t iflag) ; // Chi^2 of the fit. Should be static to be passes to MINUIT - void Unload() ; virtual const char * Version() const { return "clu-v1" ; } - + + TList* BookHists(); + void SaveHists(const char *fn="reco.root"); //*MENU* + void PrintRecoInfo(); //*MENU* + void DrawLambdasHists(); //*MENU* protected: - void WriteRecPoints() ; - virtual void MakeClusters( ) ; - + virtual void MakeClusters(); + +///////////////////// + TList *fHists; //! + TH1F* fPointE; //histogram of point energy + TH1F* fPointL1; //histogram of point L1 + TH1F* fPointL2; //histogram of point L2 + TH1F* fPointDis; //histogram of point dispersion + TH1F* fPointMult; //histogram of point multiplicity + TH1F* fDigitAmp; //histogram of digit ADC Amplitude + TH1F* fMaxE; //histogram of maximum point energy + TH1F* fMaxL1; //histogram of largest (first) of eigenvalue of covariance matrix + TH1F* fMaxL2; //histogram of smalest (second) of eigenvalue of covariace matrix + TH1F* fMaxDis; //histogram of point dispersion +/////////////////////// + + private: + AliEMCALClusterizerv1(const AliEMCALClusterizerv1 &); //copy ctor + AliEMCALClusterizerv1 & operator = (const AliEMCALClusterizerv1 &); - const TString BranchName() const ; void GetCalibrationParameters(void) ; - Bool_t FindFit(AliEMCALTowerRecPoint * emcRP, AliEMCALDigit ** MaxAt, Float_t * maxAtEnergy, + Bool_t FindFit(AliEMCALRecPoint * emcRP, AliEMCALDigit ** MaxAt, Float_t * maxAtEnergy, Int_t NPar, Float_t * FitParametres) const; //Used in UnfoldClusters, calls TMinuit void Init() ; void InitParameters() ; - virtual void MakeUnfolding() ; - void UnfoldCluster(AliEMCALTowerRecPoint * /*iniEmc*/, Int_t /*Nmax*/, + virtual void MakeUnfolding() const; + void UnfoldCluster(AliEMCALRecPoint * /*iniEmc*/, Int_t /*Nmax*/, AliEMCALDigit ** /*maxAt*/, - Float_t * /*maxAtEnergy*/ ) ; //Unfolds cluster using TMinuit package - void PrintRecPoints(Option_t * /*option*/) ; + Float_t * /*maxAtEnergy*/ ) const; //Unfolds cluster using TMinuit package + void PrintRecPoints(Option_t * option) ; private: + AliEMCALGeometry* fGeom; //! pointer to geometry for utilities Bool_t fDefaultInit; //! Says if the task was created by defaut ctor (only parameters are initialized) - - Int_t fNTowers ; // number of Towers in EMCAL - Bool_t fToUnfold ; // To perform unfolding - - Int_t fNumberOfPREClusters ; // number of clusters found in PRE section Int_t fNumberOfECAClusters ; // number of clusters found in EC section - Int_t fNumberOfHCAClusters ; // number of clusters found in HC section - + //Calibration parameters... to be replaced by database - Float_t fADCchannelPRE ; // width of one ADC channel for PRE section (GeV) - Float_t fADCpedestalPRE ; // pedestal of ADC for PRE section (GeV) + + AliEMCALCalibData * fCalibData ; //! Calibration database 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 fADCchannelHCA ; // width of one ADC channel for HC section (GeV) - Float_t fADCpedestalHCA ; // pedestal of ADC for HC section (GeV) - Float_t fECAClusteringThreshold ; // minimum energy to include a EC digit in a cluster - Float_t fHCAClusteringThreshold ; // minimum energy to include a HC digit in a cluster - Float_t fPREClusteringThreshold ; // minimum energy to include a PRE digit in 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 fHCALocMaxCut ; // minimum energy difference to distinguish local maxima in a cluster - Float_t fHCAW0 ; // logarithmic weight for the cluster center of gravity calculation - Float_t fPRELocMaxCut ; // minimum energy difference to distinguish local maxima in a CPV cluster - Float_t fPREW0 ; // logarithmic weight for the CPV cluster center of gravity calculation - Int_t fRecPointsInRun ; //! Total number of recpoints in one run - Float_t fTimeGate ; // Maximum time difference between the digits in ont EMC cluster - - ClassDef(AliEMCALClusterizerv1,2) // Clusterizer implementation version 1 + Float_t fTimeCut ; // Maximum time difference between the digits in ont EMC cluster + Float_t fMinECut; // Minimum energy for a digit to be a member of a cluster + + ClassDef(AliEMCALClusterizerv1,6) // Clusterizer implementation version 1 };