| 29c81279 |
1 | #ifndef ALIPHOSCLUSTERIZERV1_H |
| 2 | #define ALIPHOSCLUSTERIZERV1_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 | // Implementation version 1 of the clusterization algorithm |
| 10 | // |
| 11 | //*-- Author: Yves Schutz (SUBATECH) |
| 12 | |
| 13 | // --- ROOT system --- |
| 14 | |
| 15 | // --- Standard library --- |
| 16 | |
| 17 | // --- AliRoot header files --- |
| 18 | |
| 19 | #include "AliPHOSClusterizer.h" |
| 20 | #include "AliPHOSDigit.h" |
| 21 | |
| 22 | |
| 23 | |
| 24 | class AliPHOSClusterizerv1 : public AliPHOSClusterizer { |
| 25 | |
| 26 | public: |
| 27 | |
| 28 | AliPHOSClusterizerv1() ; // ctor |
| 29 | virtual ~AliPHOSClusterizerv1(){} ; // dtor |
| 30 | |
| 31 | Int_t AreNeighbours(AliPHOSDigit * d1, AliPHOSDigit * d2) ; // Checks if digits are in neighbour cells |
| 32 | Float_t Calibrate(Int_t Amp){ return (fA + fB * Amp) ;} // Tranforms Amp to energy |
| 33 | void FillandSort(const DigitsList * dl, TObjArray * tl) ; // Sorts the list according to increasing id |
| 34 | virtual void GetNumberOfClustersFound(Int_t * numb) ; |
| 35 | |
| 36 | virtual void GetCalibrationParameters(Float_t & A, Float_t &B) { A = fA; B = fB; } |
| 37 | virtual Float_t GetEmcClusteringThreshold() { return fEmcClusteringThreshold;} |
| 38 | virtual Float_t GetEmcEnergyThreshold() { return fEmcEnergyThreshold; } |
| 39 | virtual Float_t GetLocalMaxCut() { return fLocMaxCut;} |
| 40 | virtual Float_t GetLogWeightCut() { return fW0;} |
| 41 | virtual Float_t GetPpsdClusteringThreshold() { return fPpsdClusteringThreshold; } |
| 42 | virtual Float_t GetPpsdEnergyThreshold() { return fPpsdEnergyThreshold; } |
| 43 | |
| 44 | virtual Bool_t IsInEmc(AliPHOSDigit * digit) ; // Tells if id digit is in EMC |
| 45 | virtual void MakeClusters(const DigitsList * dl, RecPointsList * emcl, RecPointsList * ppsdl) ; // does the job |
| 46 | virtual void PrintParameters() ; |
| 47 | virtual void SetCalibrationParameters(Float_t A,Float_t B){ fA = A ; fB = B;} |
| 48 | virtual void SetEmcClusteringThreshold(Float_t cluth) { fEmcClusteringThreshold = cluth ; } |
| 49 | virtual void SetEmcEnergyThreshold(Float_t enth) { fEmcEnergyThreshold = enth ; } |
| 50 | virtual void SetLocalMaxCut(Float_t cut) { fLocMaxCut = cut ; } |
| 51 | virtual void SetLogWeightCut(Float_t w) { fW0 = w ; } |
| 52 | virtual void SetPpsdClusteringThreshold(Float_t cluth) { fPpsdClusteringThreshold = cluth ; } |
| 53 | virtual void SetPpsdEnergyThreshold(Float_t enth) { fPpsdEnergyThreshold = enth ; } |
| 54 | |
| 55 | private: |
| 56 | |
| 57 | Float_t fA ; // offset of the energy calibration |
| 58 | Float_t fB ; // gain of the energy calibration |
| 59 | Float_t fEmcClusteringThreshold ; // minimum energy to include a EMC digit in a cluster |
| 60 | Float_t fEmcEnergyThreshold ; // minimum energy of EMC digit to be considered |
| 61 | Float_t fLocMaxCut ; // minimum energy difference to distinguish local maxima in a cluster |
| 62 | Int_t fNumberOfEmcClusters ; // number of EMC clusters found |
| 63 | Int_t fNumberOfPpsdClusters ; // number of PPSD clusters found |
| 64 | Float_t fPpsdClusteringThreshold ; // minimum energy to include a PPSD digit in a cluster |
| 65 | Float_t fPpsdEnergyThreshold ; // minimum energy of PPSD digit to be considered |
| 66 | Float_t fW0 ; // logarithmic weight for the cluster center of gravity calculation |
| 67 | |
| 68 | ClassDef(AliPHOSClusterizerv1,1) // Clusterizer implementation version 1 |
| 69 | |
| 70 | }; |
| 71 | |
| 72 | #endif // AliPHOSCLUSTERIZERV1_H |
git://git.uio.no / u / mrichter / AliRoot.git / blame