Modification for new Eve (Marian)
[u/mrichter/AliRoot.git] / EMCAL / AliEMCALRecPoint.h
CommitLineData
ab48128d 1#ifndef ALIEMCALRECPOINT_H
2#define ALIEMCALRECPOINT_H
801f6c04 3/* Copyright(c) 1998-2007, ALICE Experiment at CERN, All rights reserved. *
ab48128d 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 ---
9aa6a5f6 13#include <TVector3.h>
37890aaf 14class TGeoManager;
15class TGeoPhysicalNode;
16class TPad;
17class TPaveText;
18class TGraph;
19class Riostream;
ab48128d 20
ab48128d 21// --- Standard library ---
22
23// --- AliRoot header files ---
24
9aa6a5f6 25#include "AliCluster.h"
1d59832c 26class AliEMCALDigit;
9aa6a5f6 27class AliDigitNew;
a5c60732 28class AliEMCALGeometry;
37890aaf 29class AliEMCALHit;
a5c60732 30
9aa6a5f6 31class AliEMCALRecPoint : public AliCluster {
ab48128d 32
33 public:
34
35 typedef TObjArray RecPointsList ;
36
37 AliEMCALRecPoint() ; // ctor
38 AliEMCALRecPoint(const char * opt) ; // ctor
0a4cb131 39 AliEMCALRecPoint(const AliEMCALRecPoint & rp);
9aa6a5f6 40
41 AliEMCALRecPoint& operator= (const AliEMCALRecPoint &rp);
42
70a93198 43 virtual ~AliEMCALRecPoint();
9aa6a5f6 44
37890aaf 45 virtual void AddDigit(AliDigitNew &) const {
9aa6a5f6 46 Fatal("AddDigit", "use AddDigit(AliEMCALDigit & digit, Float_t Energy )") ;
47 }
6e49fcaa 48 virtual void AddDigit(AliEMCALDigit & digit, Float_t Energy);
70a93198 49 virtual Int_t Compare(const TObject * obj) const;
ab48128d 50 virtual Int_t DistancetoPrimitive(Int_t px, Int_t py);
51 virtual void Draw(Option_t * option="") ;
9e5d2067 52 virtual void ExecuteEvent(Int_t event, Int_t, Int_t) ;
70a93198 53
85c60a8e 54 virtual void SetClusterType(Int_t ver) { fClusterType = ver; }
55 virtual Int_t GetClusterType() const { return fClusterType; }
56
70a93198 57 virtual void EvalAll(Float_t logWeight, TClonesArray * digits);
1ae500a2 58 virtual void EvalLocalPosition(Float_t logWeight, TClonesArray * digits);
70a93198 59 virtual void EvalPrimaries(TClonesArray * digits) ;
87cdc3be 60 virtual void EvalParents(TClonesArray * digits) ;
70a93198 61
9aa6a5f6 62 virtual int * GetDigitsList(void) const { return fDigitsList ; }
63 virtual Float_t GetEnergy() const {return fAmp; }
64
65 void EvalLocal2TrackingCSTransform();
1ae500a2 66 void EvalLocalPositionFit(Double_t deff, Double_t w0, Double_t phiSlope,TClonesArray * digits);
67 Bool_t EvalLocalPosition2(TClonesArray *digits, TArrayD &ed);
68
69 static Bool_t EvalLocalPositionFromDigits(const Double_t esum, const Double_t deff, const Double_t w0,
70 TClonesArray *digits, TArrayD &ed, TVector3 &locPos);
71 static Bool_t EvalLocalPositionFromDigits(TClonesArray *digits, TArrayD &ed, TVector3 &locPos);
72 static void GetDeffW0(const Double_t esum, Double_t &deff, Double_t &w0);
73
9aa6a5f6 74 virtual void GetGlobalPosition(TVector3 & gpos, TMatrixF & gmat) const; // return global position (x, y, z) in ALICE
70a93198 75 virtual void GetGlobalPosition(TVector3 & gpos) const; // return global position (x, y, z) in ALICE
e52475ed 76 virtual void GetLocalPosition(TVector3 & lpos) const; // return local position (x, y, z) in EMCAL SM
ab48128d 77 virtual Int_t * GetPrimaries(Int_t & number) const {number = fMulTrack ;
78 return fTracksList ; }
e52475ed 79 virtual Int_t * GetParents(Int_t & number) const {number = fMulParent ;
87cdc3be 80 return fParentsList ; }
9aa6a5f6 81
82 virtual Int_t GetDigitsMultiplicity(void) const { return fMulDigit ; }
83 Int_t GetIndexInList() const { return fIndexInList ; }
84
70a93198 85 Float_t GetCoreEnergy()const {return fCoreEnergy ;}
86 virtual Float_t GetDispersion()const {return fDispersion ;}
87 virtual void GetElipsAxis(Float_t * lambda)const {lambda[0] = fLambda[0]; lambda[1] = fLambda[1];};
88
89 Float_t * GetEnergiesList() const {return fEnergyList ;} // gets the list of energies making this recpoint
1d46d1f6 90 Double_t GetPointEnergy() const; // gets point energy (sum of energy list)
85c60a8e 91 Float_t * GetTimeList() const {return fTimeList ;} // gets the list of digit times in this recpoint
70a93198 92 Float_t GetMaximalEnergy(void) const ; // get the highest energy in the cluster
93 Int_t GetMaximumMultiplicity() const {return fMaxDigit ;} // gets the maximum number of digits allowed
94 Int_t GetMultiplicity(void) const { return fMulDigit ; } // gets the number of digits making this recpoint
95 Int_t GetMultiplicityAtLevel(Float_t level) const ; // computes multiplicity of digits with
e52475ed 96 Int_t * GetAbsId() const {return fAbsIdList;}
97 Int_t GetAbsId(int i) const {if(i>=0 && i<fMulDigit)return fAbsIdList[i]; else return -1;}
37890aaf 98 Int_t GetAbsIdMaxDigit() const {return GetAbsId(fDigitIndMax);}
99 Int_t GetIndMaxDigit() const {return fDigitIndMax;}
1ae500a2 100 void SetIndMaxDigit(const Int_t ind) {fDigitIndMax = ind;}
9aa6a5f6 101 void SetIndexInList(Int_t val) { fIndexInList = val ; }
102
103 virtual Int_t GetSuperModuleNumber(void) const { return fSuperModuleNumber;}
104
e52475ed 105 // energy above relative level
70a93198 106 virtual Int_t GetNumberOfLocalMax(AliEMCALDigit ** maxAt, Float_t * maxAtEnergy,
107 Float_t locMaxCut,TClonesArray * digits ) const ;
108 // searches for the local maxima
4635df1f 109
110 Int_t GetPrimaryIndex() const ;
70a93198 111 Float_t GetTime(void) const{return fTime ; }
112
7ee5c5be 113 virtual Bool_t IsEmc(void)const { return kTRUE ; }
ab48128d 114 virtual Bool_t IsSortable() const {
115 // tells that this is a sortable object
116 return kTRUE ;
117 }
118 virtual void Paint(Option_t * option="");
5dee926e 119 virtual void Print(Option_t * option="") const ;
7ee5c5be 120
1ae500a2 121 static Double_t TmaxInCm(const Double_t e=0.0, const Int_t key=0);
ab48128d 122
123protected:
70a93198 124 void EvalCoreEnergy(Float_t logWeight,TClonesArray * digits) ;
125 virtual void EvalDispersion(Float_t logWeight,TClonesArray * digits) ; // computes the dispersion of the shower
126 virtual void EvalElipsAxis(Float_t logWeight, TClonesArray * digits ); // computes the axis of shower ellipsoide
127 void EvalTime( TClonesArray * digits );
128 virtual Bool_t AreNeighbours(AliEMCALDigit * digit1, AliEMCALDigit * digit2 ) const;
129 Float_t ThetaToEta(Float_t arg) const; //Converts Theta (Radians) to Eta(Radians)
130 Float_t EtaToTheta(Float_t arg) const; //Converts Eta (Radians) to Theta(Radians)
131
85c60a8e 132private:
9aa6a5f6 133
134 //JLK do we need this?
14ce0a6e 135 AliEMCALGeometry* fGeomPtr; //! Pointer to geometry for utilities
a5c60732 136
9aa6a5f6 137 Float_t fAmp ; // summed amplitude of digits
138 Int_t fIndexInList ; // the index of this RecPoint in the
139 // list stored in TreeR (to be set by analysis)
140 TVector3 fLocPos ; // local position in the sub-detector coordinate
141 TMatrixF * fLocPosM ; // covariance matrix ;
142 Int_t fMaxDigit ; //! max initial size of digits array (not saved)
143 Int_t fMulDigit ; // total multiplicity of digits
144 Int_t fMaxTrack ; //! max initial size of tracks array (not saved)
145 Int_t fMulTrack ; // total multiplicity of tracks
146 Int_t * fDigitsList ; //[fMulDigit] list of digit's indexes from which the point was reconstructed
147 Int_t * fTracksList ; //[fMulTrack] list of tracks to which the point was assigned
148
149 Int_t fClusterType; // type of cluster stored: pseudocluster or v1
70a93198 150 Float_t fCoreEnergy ; // energy in a shower core
151 Float_t fLambda[2] ; // shower ellipse axes
1ae500a2 152 Float_t fDispersion ; // shower dispersion
70a93198 153 Float_t *fEnergyList ; //[fMulDigit] energy of digits
85c60a8e 154 Float_t *fTimeList ; //[fMulDigit] time of digits
e52475ed 155 Int_t *fAbsIdList; //[fMulDigit] absId of digits
70a93198 156 Float_t fTime ; // Time of the digit with maximal energy deposition
157 Float_t fCoreRadius; // The radius in which the core energy is evaluated
af5bdd85 158 Float_t *fDETracksList ; //[fMulTrack] list of tracks to which the point was assigned
159 Int_t fMulParent; // Multiplicity of the parents
87cdc3be 160 Int_t fMaxParent; // Maximum number of parents allowed
73f2fed2 161 Int_t * fParentsList; // [fMulParent] list of the parents of the digits
af5bdd85 162 Float_t * fDEParentsList; // [fMulParent] list of the parents of the digits
14ce0a6e 163 Int_t fSuperModuleNumber; // number identifying supermodule containing recpoint
1ae500a2 164 // Aug 16, 2007
165 Int_t fDigitIndMax; // Index of digit with max energy in array fAbsIdList
9aa6a5f6 166
167 ClassDef(AliEMCALRecPoint,10) // RecPoint for EMCAL (Base Class)
ab48128d 168
169};
170
171#endif // AliEMCALRECPOINT_H