]>
Commit | Line | Data |
---|---|---|
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 | |
8 | //*-- Author: Yves Schutz (SUBATECH) | |
9 | //*-- Author: Dmitri Peressounko (RRC KI & SUBATECH) | |
10 | //*-- Author: Heather Gray (LBL): merged AliEMCALRecPoint and AliEMCALTowerRecPoint 02/04 | |
11 | ||
12 | // --- ROOT system --- | |
13 | class TVector3 ; | |
14 | ||
15 | // --- Standard library --- | |
16 | ||
17 | // --- AliRoot header files --- | |
18 | ||
19 | #include "AliRecPoint.h" | |
20 | #include "AliEMCALDigit.h" | |
21 | ||
22 | class AliEMCALRecPoint : public AliRecPoint { | |
23 | ||
24 | public: | |
25 | ||
26 | typedef TObjArray RecPointsList ; | |
27 | ||
28 | AliEMCALRecPoint() ; // ctor | |
29 | AliEMCALRecPoint(const char * opt) ; // ctor | |
30 | AliEMCALRecPoint(const AliEMCALRecPoint & rp):AliRecPoint(rp) { Fatal("cpy ctor", "not implemented") ; } | |
31 | ||
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; | |
36 | virtual Int_t DistancetoPrimitive(Int_t px, Int_t py); | |
37 | virtual void Draw(Option_t * option="") ; | |
38 | virtual void ExecuteEvent(Int_t event, Int_t, Int_t) ; | |
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) ; | |
43 | virtual void EvalParents(TClonesArray * digits) ; | |
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 | |
48 | virtual Int_t * GetPrimaries(Int_t & number) const {number = fMulTrack ; | |
49 | return fTracksList ; } | |
50 | virtual Int_t * GetParents(Int_t & number) const {number = fMulParent ; | |
51 | return fParentsList ; } | |
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 | ||
67 | virtual Bool_t IsEmc(void)const { return kTRUE ; } | |
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=""); | |
73 | virtual void Print(Option_t * ) const ; | |
74 | ||
75 | AliEMCALRecPoint & operator = (const AliEMCALRecPoint & ) { | |
76 | Fatal("operator =", "not implemented") ; | |
77 | return *this ; | |
78 | } | |
79 | ||
80 | protected: | |
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 | |
95 | Int_t fMulParent; // Multiplicity of the parents | |
96 | Int_t fMaxParent; // Maximum number of parents allowed | |
97 | Int_t * fParentsList; // [fMaxParent] list of the parents of the digits | |
98 | ||
99 | ClassDef(AliEMCALRecPoint,6) // RecPoint for EMCAL (Base Class) | |
100 | ||
101 | }; | |
102 | ||
103 | #endif // AliEMCALRECPOINT_H |