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 GetPpsdClusteringThreshold() { return fPpsdClusteringThreshold;  }
  43   virtual Float_t GetPpsdEnergyThreshold() { return  fPpsdEnergyThreshold;  }
  44
  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 ; }
  57
  58 private:
  59
  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
  71
  72   ClassDef(AliPHOSClusterizerv1,1)  // Clusterizer implementation version 1
  73
  74 };
  75
  76 #endif // AliPHOSCLUSTERIZERV1_H