update of Jenn and Marco
[u/mrichter/AliRoot.git] / EMCAL / AliEMCALRecPoint.h
CommitLineData
ab48128d 1#ifndef ALIEMCALRECPOINT_H
2#define ALIEMCALRECPOINT_H
3/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
5//_________________________________________________________________________
6// Base Class for EMCAL Reconstructed Points
7// A recpoint being equivalent to a cluster in encal terminology
d64c959b 8//*-- Author: Yves Schutz (SUBATECH)
70a93198 9//*-- Author: Dmitri Peressounko (RRC KI & SUBATECH)
10//*-- Author: Heather Gray (LBL): merged AliEMCALRecPoint and AliEMCALTowerRecPoint 02/04
ab48128d 11
12// --- ROOT system ---
70a93198 13class TVector3 ;
ab48128d 14
ab48128d 15// --- Standard library ---
16
17// --- AliRoot header files ---
18
19#include "AliRecPoint.h"
20#include "AliEMCALDigit.h"
21
22class AliEMCALRecPoint : public AliRecPoint {
23
24 public:
25
26 typedef TObjArray RecPointsList ;
27
28 AliEMCALRecPoint() ; // ctor
29 AliEMCALRecPoint(const char * opt) ; // ctor
d64c959b 30 AliEMCALRecPoint(const AliEMCALRecPoint & rp):AliRecPoint(rp) { Fatal("cpy ctor", "not implemented") ; }
ab48128d 31
70a93198 32 virtual ~AliEMCALRecPoint();
33 virtual void AddDigit(AliDigitNew &){ Fatal("AddDigit", "use AddDigit(AliEMCALDigit & digit, Float_t Energy )") ; }
34 virtual void AddDigit(AliEMCALDigit & digit, Float_t Energy);
35 virtual Int_t Compare(const TObject * obj) const;
ab48128d 36 virtual Int_t DistancetoPrimitive(Int_t px, Int_t py);
37 virtual void Draw(Option_t * option="") ;
9e5d2067 38 virtual void ExecuteEvent(Int_t event, Int_t, Int_t) ;
70a93198 39
40 virtual void EvalAll(Float_t logWeight, TClonesArray * digits);
41 virtual void EvalLocalPosition(Float_t logWeight, TClonesArray * digits) ;
42 virtual void EvalPrimaries(TClonesArray * digits) ;
87cdc3be 43 virtual void EvalParents(TClonesArray * digits) ;
70a93198 44
45 // virtual void GetGlobalPosition(TVector3 & gpos, TMatrix & /*gmat*/) const; // return global position in ALICE
46 virtual void GetGlobalPosition(TVector3 & gpos) const; // return global position (x, y, z) in ALICE
47 virtual void GetLocalPosition(TVector3 & lpos) const; // return local position (eta, phi, r) in EMCAL
ab48128d 48 virtual Int_t * GetPrimaries(Int_t & number) const {number = fMulTrack ;
49 return fTracksList ; }
87cdc3be 50 virtual Int_t * GetParents(Int_t & number) const {number = fMulParent ;
51 return fParentsList ; }
70a93198 52 Float_t GetCoreEnergy()const {return fCoreEnergy ;}
53 virtual Float_t GetDispersion()const {return fDispersion ;}
54 virtual void GetElipsAxis(Float_t * lambda)const {lambda[0] = fLambda[0]; lambda[1] = fLambda[1];};
55
56 Float_t * GetEnergiesList() const {return fEnergyList ;} // gets the list of energies making this recpoint
57 Float_t GetMaximalEnergy(void) const ; // get the highest energy in the cluster
58 Int_t GetMaximumMultiplicity() const {return fMaxDigit ;} // gets the maximum number of digits allowed
59 Int_t GetMultiplicity(void) const { return fMulDigit ; } // gets the number of digits making this recpoint
60 Int_t GetMultiplicityAtLevel(Float_t level) const ; // computes multiplicity of digits with
61 // energy above relative level
62 virtual Int_t GetNumberOfLocalMax(AliEMCALDigit ** maxAt, Float_t * maxAtEnergy,
63 Float_t locMaxCut,TClonesArray * digits ) const ;
64 // searches for the local maxima
65 Float_t GetTime(void) const{return fTime ; }
66
7ee5c5be 67 virtual Bool_t IsEmc(void)const { return kTRUE ; }
ab48128d 68 virtual Bool_t IsSortable() const {
69 // tells that this is a sortable object
70 return kTRUE ;
71 }
72 virtual void Paint(Option_t * option="");
5dee926e 73 virtual void Print(Option_t * option="") const ;
7ee5c5be 74
ab48128d 75 AliEMCALRecPoint & operator = (const AliEMCALRecPoint & ) {
d64c959b 76 Fatal("operator =", "not implemented") ;
ab48128d 77 return *this ;
78 }
79
80protected:
70a93198 81 void EvalCoreEnergy(Float_t logWeight,TClonesArray * digits) ;
82 virtual void EvalDispersion(Float_t logWeight,TClonesArray * digits) ; // computes the dispersion of the shower
83 virtual void EvalElipsAxis(Float_t logWeight, TClonesArray * digits ); // computes the axis of shower ellipsoide
84 void EvalTime( TClonesArray * digits );
85 virtual Bool_t AreNeighbours(AliEMCALDigit * digit1, AliEMCALDigit * digit2 ) const;
86 Float_t ThetaToEta(Float_t arg) const; //Converts Theta (Radians) to Eta(Radians)
87 Float_t EtaToTheta(Float_t arg) const; //Converts Eta (Radians) to Theta(Radians)
88
89 Float_t fCoreEnergy ; // energy in a shower core
90 Float_t fLambda[2] ; // shower ellipse axes
91 Float_t fDispersion ; // shower dispersion
92 Float_t *fEnergyList ; //[fMulDigit] energy of digits
93 Float_t fTime ; // Time of the digit with maximal energy deposition
94 Float_t fCoreRadius; // The radius in which the core energy is evaluated
87cdc3be 95 Int_t fMulParent; // Multiplicity of the parents
96 Int_t fMaxParent; // Maximum number of parents allowed
73f2fed2 97 Int_t * fParentsList; // [fMulParent] list of the parents of the digits
ab48128d 98
87cdc3be 99 ClassDef(AliEMCALRecPoint,6) // RecPoint for EMCAL (Base Class)
ab48128d 100
101};
102
103#endif // AliEMCALRECPOINT_H