| d15a28e7 |
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 | |
| 6ad0bfa0 |
6 | /* $Id$ */ |
| 7 | |
| b2a60966 |
8 | //_________________________________________________________________________ |
| 9 | // Implementation version 1 of the clusterization algorithm |
| 10 | // |
| 11 | //*-- Author: Yves Schutz (SUBATECH) |
| d15a28e7 |
12 | |
| 13 | // --- ROOT system --- |
| 14 | |
| 15 | // --- Standard library --- |
| 16 | |
| 17 | // --- AliRoot header files --- |
| 18 | |
| 19 | #include "AliPHOSClusterizer.h" |
| 20 | #include "AliPHOSDigit.h" |
| 6b9b768e |
21 | #include "AliPHOSGeometry.h" |
| d15a28e7 |
22 | |
| 23 | |
| 24 | |
| 25 | class AliPHOSClusterizerv1 : public AliPHOSClusterizer { |
| 26 | |
| 27 | public: |
| 28 | |
| 29 | AliPHOSClusterizerv1() ; // ctor |
| 30 | virtual ~AliPHOSClusterizerv1(){} ; // dtor |
| 31 | |
| 9f616d61 |
32 | Int_t AreNeighbours(AliPHOSDigit * d1, AliPHOSDigit * d2) ; // Checks if digits are in neighbour cells |
| 33 | Float_t Calibrate(Int_t Amp){ return (fA + fB * Amp) ;} // Tranforms Amp to energy |
| 34 | void FillandSort(const DigitsList * dl, TObjArray * tl) ; // Sorts the list according to increasing id |
| 4a646085 |
35 | virtual void GetNumberOfClustersFound(Int_t * numb) ; |
| 36 | |
| 37 | virtual void GetCalibrationParameters(Float_t & A, Float_t &B) { A = fA; B = fB; } |
| d72dfbc3 |
38 | virtual Float_t GetEmcClusteringThreshold() { return fEmcClusteringThreshold;} |
| 39 | virtual Float_t GetEmcEnergyThreshold() { return fEmcEnergyThreshold; } |
| 40 | virtual Float_t GetLocalMaxCut() { return fLocMaxCut;} |
| 41 | virtual Float_t GetLogWeightCut() { return fW0;} |
| 42 | virtual Float_t GetLocalMaxCutCPV() { return fLocMaxCutCPV;} |
| 43 | virtual Float_t GetLogWeightCutCPV() { return fW0CPV;} |
| 416692f4 |
44 | virtual Float_t GetPpsdClusteringThreshold() { return fPpsdClusteringThreshold; } |
| d72dfbc3 |
45 | virtual Float_t GetPpsdEnergyThreshold() { return fPpsdEnergyThreshold; } |
| ed4205d8 |
46 | virtual Float_t GetCpvClusteringThreshold() { return fCpvClusteringThreshold; } |
| 47 | virtual Float_t GetCpvEnergyThreshold() { return fCpvEnergyThreshold; } |
| 4a646085 |
48 | |
| ed4205d8 |
49 | virtual Bool_t IsInEmc (AliPHOSDigit * digit) ; // Tells if id digit is in EMC |
| 50 | virtual Bool_t IsInPpsd(AliPHOSDigit * digit) ; // Tells if id digit is in PPSD |
| 51 | virtual Bool_t IsInCpv (AliPHOSDigit * digit) ; // Tells if id digit is in CPV |
| 88714635 |
52 | virtual void MakeClusters(const DigitsList * dl, |
| 53 | AliPHOSRecPoint::RecPointsList * emcl, |
| 54 | AliPHOSRecPoint::RecPointsList * ppsdl) ; // does the job |
| 6ad0bfa0 |
55 | virtual void PrintParameters() ; |
| 56 | virtual void SetCalibrationParameters(Float_t A,Float_t B){ fA = A ; fB = B;} |
| d72dfbc3 |
57 | virtual void SetEmcClusteringThreshold(Float_t cluth) { fEmcClusteringThreshold = cluth ; } |
| 58 | virtual void SetEmcEnergyThreshold(Float_t enth) { fEmcEnergyThreshold = enth ; } |
| 59 | virtual void SetLocalMaxCut(Float_t cut) { fLocMaxCut = cut ; } |
| 60 | virtual void SetLogWeightCut(Float_t w) { fW0 = w ; } |
| 61 | virtual void SetLocalMaxCutCPV(Float_t cut) { fLocMaxCutCPV = cut ; } |
| 62 | virtual void SetLogWeightCutCPV(Float_t w) { fW0CPV = w ; } |
| 6ad0bfa0 |
63 | virtual void SetPpsdClusteringThreshold(Float_t cluth) { fPpsdClusteringThreshold = cluth ; } |
| d72dfbc3 |
64 | virtual void SetPpsdEnergyThreshold(Float_t enth) { fPpsdEnergyThreshold = enth ; } |
| ed4205d8 |
65 | virtual void SetCpvClusteringThreshold(Float_t cluth) { fCpvClusteringThreshold = cluth ; } |
| 66 | virtual void SetCpvEnergyThreshold(Float_t enth) { fCpvEnergyThreshold = enth ; } |
| d15a28e7 |
67 | |
| 68 | private: |
| 69 | |
| b2a60966 |
70 | Float_t fA ; // offset of the energy calibration |
| 71 | Float_t fB ; // gain of the energy calibration |
| 6b9b768e |
72 | AliPHOSGeometry * fGeom ; // pointer to geometry |
| ed4205d8 |
73 | Int_t fNumberOfEmcClusters ; // number of EMC clusters found |
| b2a60966 |
74 | Float_t fEmcClusteringThreshold ; // minimum energy to include a EMC digit in a cluster |
| 75 | Float_t fEmcEnergyThreshold ; // minimum energy of EMC digit to be considered |
| b2a60966 |
76 | Int_t fNumberOfPpsdClusters ; // number of PPSD clusters found |
| 77 | Float_t fPpsdClusteringThreshold ; // minimum energy to include a PPSD digit in a cluster |
| 78 | Float_t fPpsdEnergyThreshold ; // minimum energy of PPSD digit to be considered |
| ed4205d8 |
79 | Int_t fNumberOfCpvClusters ; // number of CPV clusters found |
| 80 | Float_t fCpvClusteringThreshold ; // minimum energy to include a CPV digit in a cluster |
| 81 | Float_t fCpvEnergyThreshold ; // minimum energy of CPV digit to be considered |
| d72dfbc3 |
82 | Float_t fLocMaxCut ; // minimum energy difference to distinguish local maxima in a cluster |
| b2a60966 |
83 | Float_t fW0 ; // logarithmic weight for the cluster center of gravity calculation |
| d72dfbc3 |
84 | Float_t fLocMaxCutCPV ; // minimum energy difference to distinguish local maxima in a CPV cluster |
| 85 | Float_t fW0CPV ; // logarithmic weight for the CPV cluster center of gravity calculation |
| 92862013 |
86 | |
| ed4205d8 |
87 | ClassDef(AliPHOSClusterizerv1,1) // Clusterizer implementation version 1 |
| d15a28e7 |
88 | |
| 89 | }; |
| 90 | |
| 91 | #endif // AliPHOSCLUSTERIZERV1_H |
git://git.uio.no / u / mrichter / AliRoot.git / blame