Reconstruction part now handle all geometry options
[u/mrichter/AliRoot.git] / PHOS / AliPHOSClusterizerv1.h
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 #include "AliPHOSGeometry.h" 
22
23
24
25 class AliPHOSClusterizerv1 : public AliPHOSClusterizer {
26   
27 public:
28   
29   AliPHOSClusterizerv1() ;             // ctor            
30   virtual ~AliPHOSClusterizerv1(){} ;  // dtor
31   
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) ; 
36  
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 GetLocalMaxCutCPV()          { return fLocMaxCutCPV;} 
43   virtual Float_t GetLogWeightCutCPV()         { return fW0CPV;}  
44   virtual Float_t GetPpsdClusteringThreshold() { return fPpsdClusteringThreshold;  } 
45   virtual Float_t GetPpsdEnergyThreshold()     { return fPpsdEnergyThreshold;  }
46   virtual Float_t GetCpvClusteringThreshold()  { return fCpvClusteringThreshold;  } 
47   virtual Float_t GetCpvEnergyThreshold()      { return fCpvEnergyThreshold;  }
48
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
52   virtual void MakeClusters(const DigitsList * dl, 
53                             AliPHOSRecPoint::RecPointsList * emcl, 
54                             AliPHOSRecPoint::RecPointsList * ppsdl) ; // does the job 
55   virtual void PrintParameters() ;
56   virtual void SetCalibrationParameters(Float_t A,Float_t B){ fA = A ; fB = B;} 
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 ; }
63   virtual void SetPpsdClusteringThreshold(Float_t cluth) { fPpsdClusteringThreshold = cluth ; }
64   virtual void SetPpsdEnergyThreshold(Float_t enth)      { fPpsdEnergyThreshold = enth ; } 
65   virtual void SetCpvClusteringThreshold(Float_t cluth)  { fCpvClusteringThreshold = cluth ; }
66   virtual void SetCpvEnergyThreshold(Float_t enth)       { fCpvEnergyThreshold = enth ; } 
67   
68 private:
69   
70   Float_t fA ;                       // offset of the energy calibration
71   Float_t fB ;                       // gain of the energy calibration
72   AliPHOSGeometry * fGeom ;          // pointer to geometry
73   Int_t   fNumberOfEmcClusters ;     // number of EMC clusters found 
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
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
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
82   Float_t fLocMaxCut ;               // minimum energy difference to distinguish local maxima in a cluster
83   Float_t fW0 ;                      // logarithmic weight for the cluster center of gravity calculation
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
86     
87   ClassDef(AliPHOSClusterizerv1,1)   // Clusterizer implementation version 1
88
89 };
90
91 #endif // AliPHOSCLUSTERIZERV1_H