1 #ifndef ALIESDCALOCLUSTER_H
2 #define ALIESDCALOCLUSTER_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
8 //-------------------------------------------------------------------------
9 // Class AliESDCaloCluster
10 // This is the class to deal with during the physics analysis of data
12 // New container for calorimeter clusters, which are the effective
13 // "tracks" for calorimeter detectors. Can be used by PHOS and EMCAL
16 //-------------------------------------------------------------------------
25 class AliESDCaloCluster : public TObject
31 AliESDCaloCluster(const AliESDCaloCluster& clus);
32 AliESDCaloCluster & operator=(const AliESDCaloCluster& source);
33 virtual ~AliESDCaloCluster();
34 virtual void Copy(TObject &) const;
36 void SetID(Int_t id) {fID = id;}
37 Int_t GetID() const {return fID;}
39 //similar to AliAODCluster but offset by one for
40 // backward comp. -1 was undefined, which only applied
41 // for PHOS clusters before
42 enum ESDClu_t {kUndef = -2,
44 kEMCALPseudoCluster, //Not any more in use, keep for backward comp.
47 void SetClusterType(Int_t type) { fClusterType = type; }
48 Char_t GetClusterType() const {return fClusterType; }
50 Bool_t IsEMCAL() const {return (fClusterType == kEMCALClusterv1);}
51 Bool_t IsPHOS() const {return (fClusterType == kPHOSCluster);}
53 void SetPosition(const Float_t *pos) {
54 fGlobalPos[0] = pos[0]; fGlobalPos[1] = pos[1]; fGlobalPos[2] = pos[2];
56 void GetPosition(Float_t *pos) const {
57 pos[0] = fGlobalPos[0]; pos[1] = fGlobalPos[1]; pos[2] = fGlobalPos[2];
60 void SetE(Float_t ene) { fEnergy = ene;}
61 Double_t E() const { return fEnergy;}
63 void SetClusterDisp(Float_t disp) { fDispersion = disp; }
64 Double_t GetClusterDisp() const { return fDispersion; }
66 void SetClusterChi2(Float_t chi2) { fChi2 = chi2; }
67 Double_t GetClusterChi2() const { return fChi2; }
69 void SetPid(const Float_t *p);
70 Double_t *GetPid() {return fPID;}
72 void SetM20(Float_t m20) { fM20 = m20; }
73 Double_t GetM20() const { return fM20; }
75 void SetM02(Float_t m02) { fM02 = m02; }
76 Double_t GetM02() const { return fM02; }
78 void SetNExMax(UChar_t nExMax) { fNExMax = nExMax; }
79 UChar_t GetNExMax() const { return fNExMax; }
81 void SetEmcCpvDistance(Float_t dEmcCpv) { fEmcCpvDistance = dEmcCpv; }
82 Double_t GetEmcCpvDistance() const { return fEmcCpvDistance; }
83 void SetTrackDistance(Double_t dx, Double_t dz){fTrackDx=dx; fTrackDz=dz;}
84 Double_t GetTrackDx(void)const {return fTrackDx;}
85 Double_t GetTrackDz(void)const {return fTrackDz;}
87 void SetDistanceToBadChannel(Float_t dist) {fDistToBadChannel=dist;}
88 Double_t GetDistanceToBadChannel() const {return fDistToBadChannel;}
90 void SetTOF(Double_t tof) { fTOF = tof; }
91 Double_t GetTOF() const { return fTOF; }
93 void AddTracksMatched(TArrayI & array) {
94 if(!fTracksMatched)fTracksMatched = new TArrayI(array);
95 else *fTracksMatched = array;
97 void AddLabels(TArrayI & array) {
98 if(!fLabels)fLabels = new TArrayI(array) ;
99 else *fLabels = array;
102 TArrayI * GetTracksMatched() const {return fTracksMatched;}
103 TArrayI * GetLabels() const {return fLabels;}
105 Int_t GetTrackMatched() const
106 {if( fTracksMatched && fTracksMatched->GetSize() >0) return fTracksMatched->At(0);
107 else return -1;} //Most likely the track associated to the cluster
108 Int_t GetLabel() const
109 {if( fLabels && fLabels->GetSize() >0) return fLabels->At(0);
110 else return -1;} //Most likely the track associated to the cluster
112 Int_t GetNTracksMatched() const {if (fTracksMatched) return fTracksMatched->GetSize();
114 Int_t GetNLabels() const { if (fLabels) return fLabels->GetSize();
117 void GetMomentum(TLorentzVector& p, Double_t * vertexPosition );
120 void SetNCells(Int_t n) { fNCells = n;}
121 Int_t GetNCells() const { return fNCells;}
123 void SetCellsAbsId(UShort_t *array) ;
124 UShort_t *GetCellsAbsId() {return fCellsAbsId;}
126 void SetCellsAmplitudeFraction(Double32_t *array) ;
127 Double32_t *GetCellsAmplitudeFraction() {return fCellsAmpFraction;}
129 Int_t GetCellAbsId(Int_t i) const {
130 if (fCellsAbsId && i >=0 && i < fNCells ) return fCellsAbsId[i];
133 Double_t GetCellAmplitudeFraction(Int_t i) const {
134 if (fCellsAmpFraction && i >=0 && i < fNCells ) return fCellsAmpFraction[i];
137 //_____________________________________________________
138 //Not used anymore, kept to avoid backward incompatibility
139 void AddDigitIndex(TArrayS & array) { fDigitIndex = new TArrayS(array) ; Warning("AddDigitAmplitude","This method is no more in use") ;}
140 void AddDigitAmplitude(TArrayS & array) { fDigitAmplitude = new TArrayS(array) ; Warning("AddDigitAmplitude","This method is no more in use") ;}
141 void AddDigitTime(TArrayS & array) { fDigitTime = new TArrayS(array) ;Warning("AddDigitTime","This method is no more in use") ;}
142 TArrayS * GetDigitAmplitude() const {return fDigitAmplitude;}
143 TArrayS * GetDigitTime() const {return fDigitTime;}
144 TArrayS * GetDigitIndex() const {return fDigitIndex;}
145 Int_t GetNumberOfDigits() const { return -1;}
146 //_____________________________________________________
150 TArrayI * fTracksMatched; //Index of tracks close to cluster. First entry is the most likely match.
151 TArrayI * fLabels; //list of primaries that generated the cluster, ordered in deposited energy.
155 UShort_t *fCellsAbsId; //[fNCells] array of cell absId numbers
156 Double32_t *fCellsAmpFraction; //[fNCells][0.,1.,16] array with cell amplitudes fraction.
158 //__________________________________________________________
160 TArrayS * fDigitAmplitude; //digit energy (integer units)
161 TArrayS * fDigitTime; //time of this digit (integer units)
162 TArrayS * fDigitIndex; //calorimeter digit index
163 //_________________________________________________________
165 Double32_t fGlobalPos[3]; // position in global coordinate systemD
166 Double32_t fEnergy; // energy measured by calorimeter
167 Double32_t fDispersion; // cluster dispersion, for shape analysis
168 Double32_t fChi2; // chi2 of cluster fi
169 Double32_t fM20; // 2-nd moment along the main eigen axis
170 Double32_t fM02; // 2-nd moment along the second eigen axis
172 Double32_t fEmcCpvDistance; // the distance from PHOS EMC rec.point to the closest CPV rec.point
173 Double32_t fTrackDx ; // Distance to closest track in phi
174 Double32_t fTrackDz ; // Distance to closest track in z
176 Double32_t fDistToBadChannel; // Distance to nearest bad channel
177 Double32_t fPID[AliPID::kSPECIESN]; //[0,1,8]"detector response probabilities" (for the PID)
178 Int_t fID; // Unique Id of the cluster
179 UChar_t fNExMax ; // number of (Ex-)maxima before unfolding
180 Char_t fClusterType; // Flag for different cluster type/versions
181 Double32_t fTOF; //[0,0,12] time-of-flight
183 ClassDef(AliESDCaloCluster,9) //ESDCaloCluster