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 */
8 //_________________________________________________________________________
9 // Implementation version 1 of the clusterization algorithm
11 //*-- Author: Yves Schutz (SUBATECH)
13 // --- ROOT system ---
15 // --- Standard library ---
17 // --- AliRoot header files ---
19 #include "AliPHOSClusterizer.h"
20 #include "AliPHOSDigit.h"
21 #include "AliPHOSGeometry.h"
25 class AliPHOSClusterizerv1 : public AliPHOSClusterizer {
29 AliPHOSClusterizerv1() ; // ctor
30 virtual ~AliPHOSClusterizerv1(){} ; // dtor
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
35 virtual void GetNumberOfClustersFound(Int_t * numb) ;
37 virtual void GetCalibrationParameters(Float_t & A, Float_t &B) { A = fA; B = fB; }
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 GetPpsdClusteringThreshold() { return fPpsdClusteringThreshold; }
43 virtual Float_t GetPpsdEnergyThreshold() { return fPpsdEnergyThreshold; }
45 virtual Bool_t IsInEmc(AliPHOSDigit * digit) ; // Tells if id digit is in EMC
46 virtual void MakeClusters(const DigitsList * dl,
47 AliPHOSRecPoint::RecPointsList * emcl,
48 AliPHOSRecPoint::RecPointsList * ppsdl) ; // does the job
49 virtual void PrintParameters() ;
50 virtual void SetCalibrationParameters(Float_t A,Float_t B){ fA = A ; fB = B;}
51 virtual void SetEmcClusteringThreshold(Float_t cluth) { fEmcClusteringThreshold = cluth ; }
52 virtual void SetEmcEnergyThreshold(Float_t enth) { fEmcEnergyThreshold = enth ; }
53 virtual void SetLocalMaxCut(Float_t cut) { fLocMaxCut = cut ; }
54 virtual void SetLogWeightCut(Float_t w) { fW0 = w ; }
55 virtual void SetPpsdClusteringThreshold(Float_t cluth) { fPpsdClusteringThreshold = cluth ; }
56 virtual void SetPpsdEnergyThreshold(Float_t enth) { fPpsdEnergyThreshold = enth ; }
60 Float_t fA ; // offset of the energy calibration
61 Float_t fB ; // gain of the energy calibration
62 AliPHOSGeometry * fGeom ; // pointer to geometry
63 Float_t fEmcClusteringThreshold ; // minimum energy to include a EMC digit in a cluster
64 Float_t fEmcEnergyThreshold ; // minimum energy of EMC digit to be considered
65 Float_t fLocMaxCut ; // minimum energy difference to distinguish local maxima in a cluster
66 Int_t fNumberOfEmcClusters ; // number of EMC clusters found
67 Int_t fNumberOfPpsdClusters ; // number of PPSD clusters found
68 Float_t fPpsdClusteringThreshold ; // minimum energy to include a PPSD digit in a cluster
69 Float_t fPpsdEnergyThreshold ; // minimum energy of PPSD digit to be considered
70 Float_t fW0 ; // logarithmic weight for the cluster center of gravity calculation
72 ClassDef(AliPHOSClusterizerv1,1) // Clusterizer implementation version 1
76 #endif // AliPHOSCLUSTERIZERV1_H