PHOS/AliPHOSEmcRecPoint.h

   1 #ifndef ALIPHOSEMCRECPOINT_H
   2 #define ALIPHOSEMCRECPOINT_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 //  RecPoint implementation for PHOS-EMC
  10 //  An EmcRecPoint is a cluster of digits
  11 //
  12 //*-- Author: Dmitri Peressounko (RRC KI & SUBATECH)
  13
  14 // --- ROOT system ---
  15
  16 #include "TObject.h"
  17 #include "TArrayI.h"
  18
  19 // --- Standard library ---
  20
  21 // --- AliRoot header files ---
  22
  23 #include "AliPHOSDigit.h"
  24 #include "AliPHOSRecPoint.h"
  25
  26 class AliPHOSEmcRecPoint : public AliPHOSRecPoint  {
  27
  28 public:
  29
  30   AliPHOSEmcRecPoint(){
  31    // default ctor
  32     fEnergyList = 0;
  33   } ;
  34   AliPHOSEmcRecPoint(Float_t W0, Float_t LocMaxCut) ;
  35   AliPHOSEmcRecPoint(const AliPHOSEmcRecPoint & rp) {
  36     // cpy ctor requested by Coding Convention
  37     // but not yet needed
  38     assert(0==1) ;
  39   }
  40
  41   virtual ~AliPHOSEmcRecPoint() ;
  42
  43   virtual void  AddDigit(AliPHOSDigit & digit, Float_t Energy) ;  // add a digit to the digits list
  44   Int_t       Compare(TObject * obj) ;                         // method for sorting
  45   Float_t CoreEnergy() ;
  46
  47   Float_t     GetDelta (){
  48     // gets the fDelta data member
  49     return fDelta ; }
  50   Float_t     GetDispersion() ;                               // computes the dispersion of the shower
  51   void        GetElipsAxis(Float_t * lambda) ;                // computes the axis of shower ellipsoide
  52   Float_t *   GetEnergiesList(){
  53     // gets the list of energies makink this recpoint
  54     return fEnergyList ;}
  55   virtual void ExecuteEvent(Int_t event, Int_t px, Int_t py) ;
  56   Float_t     GetLocMaxCut () {
  57     // gets the cut of the local maximum search
  58     return fLocMaxCut ; }
  59   Float_t     GetLogWeightCut (){
  60     // gets the logarythmic weight for the center of gravity calculation
  61     return fW0 ; }
  62   Float_t     GetMaximalEnergy(void) ;                        // get the highest energy in the cluster
  63   Int_t       GetMaximumMultiplicity() {
  64     // gets the maximum number of digits allowed
  65     return   fMaxDigit ; }
  66   Int_t       GetMultiplicity(void) const {
  67     // gets the number of digits making this recpoint
  68     return fMulDigit ; }
  69   Int_t       GetMultiplicityAtLevel(const Float_t level) ;   // computes multiplicity of digits with energy above relative level
  70   Int_t       GetNumberOfLocalMax(Int_t *  maxAt, Float_t * maxAtEnergy) ; // searches for the local maxima
  71
  72   Float_t     GetTotalEnergy(void) const {
  73     // gets the total amplitude of this recpoint (in EMC RecPoint Amp = Energy)
  74     return fAmp ; }
  75   void        GetLocalPosition(TVector3 &Lpos) ;  // computes the position in the PHOS module
  76   Bool_t      IsEmc(void) {
  77     // true if the recpoint is in EMC
  78     return kTRUE ; }
  79   Bool_t      IsSortable() const {
  80     // says that emcrecpoints are sortable objects
  81     return kTRUE ; }
  82   void        Print(Option_t * opt = "void") ;
  83
  84   AliPHOSEmcRecPoint & operator = (const AliPHOSEmcRecPoint & rvalue)  {
  85     // assignement operator requested by coding convention
  86     // but not needed
  87     assert(0==1) ;
  88     return *this ;
  89   }
  90
  91  private:
  92
  93   Bool_t AreNeighbours(AliPHOSDigit * digit1, AliPHOSDigit * digit2 ) ;
  94
  95   Float_t  fDelta ;          // parameter used to sort the clusters
  96   Float_t  *fEnergyList ;    //[fMulDigit] energy of digits
  97   Float_t  fLocMaxCut ;      // minimum energy difference to distinguish two maxima
  98   Float_t  fW0 ;             // logarithmic weight factor for center of gravity calculation
  99
 100   ClassDef(AliPHOSEmcRecPoint,1)  // EMC RecPoint (cluster)
 101
 102 };
 103
 104 #endif // AliPHOSEMCRECPOINT_H