1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
19 //-----------------------------------------------------------------
20 // Implementation of the ESD Calorimeter cluster class
21 // ESD = Event Summary Data
22 // This is the class to deal with during the phisics analysis of data
25 //-----------------------------------------------------------------
27 #include <TLorentzVector.h>
28 #include "AliESDCaloCluster.h"
30 ClassImp(AliESDCaloCluster)
32 //_______________________________________________________________________
33 AliESDCaloCluster::AliESDCaloCluster() :
39 fCellsAmpFraction(0x0),
40 fDigitAmplitude(0x0),//not in use
41 fDigitTime(0x0),//not in use
42 fDigitIndex(0x0),//not in use
48 fEmcCpvDistance(1024),
49 fDistToBadChannel(1024),
52 fClusterType(kUndef), fTOF(0.)
55 // The default ESD constructor
57 fGlobalPos[0] = fGlobalPos[1] = fGlobalPos[2] = 0.;
58 for(Int_t i=0; i<AliPID::kSPECIESN; i++) fPID[i] = 0.;
61 //_______________________________________________________________________
62 AliESDCaloCluster::AliESDCaloCluster(const AliESDCaloCluster& clus) :
64 fTracksMatched(clus.fTracksMatched?new TArrayI(*clus.fTracksMatched):0x0),
65 fLabels(clus.fLabels?new TArrayI(*clus.fLabels):0x0),
66 fNCells(clus.fNCells),
69 fDigitAmplitude(clus.fDigitAmplitude?new TArrayS(*clus.fDigitAmplitude):0x0),//not in use
70 fDigitTime(clus.fDigitTime?new TArrayS(*clus.fDigitTime):0x0),//not in use
71 fDigitIndex(clus.fDigitIndex?new TArrayS(*clus.fDigitIndex):0x0),//not in use
72 fEnergy(clus.fEnergy),
73 fDispersion(clus.fDispersion),
77 fEmcCpvDistance(clus.fEmcCpvDistance),
78 fDistToBadChannel(clus.fDistToBadChannel),
80 fNExMax(clus.fNExMax),
81 fClusterType(clus.fClusterType),
85 // The copy constructor
87 fGlobalPos[0] = clus.fGlobalPos[0];
88 fGlobalPos[1] = clus.fGlobalPos[1];
89 fGlobalPos[2] = clus.fGlobalPos[2];
91 for(Int_t i=0; i<AliPID::kSPECIESN; i++) fPID[i] = clus.fPID[i];
93 if (clus.fNCells > 0) {
96 fCellsAbsId = new UShort_t[clus.fNCells];
97 for (Int_t i=0; i<clus.fNCells; i++)
98 fCellsAbsId[i]=clus.fCellsAbsId[i];
101 if(clus.fCellsAmpFraction){
102 fCellsAmpFraction = new Double32_t[clus.fNCells];
103 for (Int_t i=0; i<clus.fNCells; i++)
104 fCellsAmpFraction[i]=clus.fCellsAmpFraction[i];
111 //_______________________________________________________________________
112 AliESDCaloCluster &AliESDCaloCluster::operator=(const AliESDCaloCluster& source)
114 // assignment operator
116 if(&source == this) return *this;
117 TObject::operator=(source);
118 fGlobalPos[0] = source.fGlobalPos[0];
119 fGlobalPos[1] = source.fGlobalPos[1];
120 fGlobalPos[2] = source.fGlobalPos[2];
122 fEnergy = source.fEnergy;
123 fDispersion = source.fDispersion;
124 fChi2 = source.fChi2;
127 fEmcCpvDistance = source.fEmcCpvDistance;
128 fDistToBadChannel = source.fDistToBadChannel ;
129 for(Int_t i=0; i<AliPID::kSPECIESN; i++) fPID[i] = source.fPID[i];
132 fNCells= source.fNCells;
134 if (source.fNCells > 0) {
135 if(source.fCellsAbsId){
136 if(fNCells != source.fNCells||!fCellsAbsId){
137 if(fCellsAbsId)delete [] fCellsAbsId;
138 fCellsAbsId = new UShort_t[source.fNCells];
140 for (Int_t i=0; i<source.fNCells; i++){
141 fCellsAbsId[i]=source.fCellsAbsId[i];
145 if(source.fCellsAmpFraction){
146 if(fNCells != source.fNCells||!fCellsAmpFraction){
147 if(fCellsAmpFraction) delete [] fCellsAmpFraction;
148 fCellsAmpFraction = new Double32_t[source.fNCells];
150 for (Int_t i=0; i<source.fNCells; i++)
151 fCellsAmpFraction[i]=source.fCellsAmpFraction[i];
155 fNExMax = source.fNExMax;
156 fClusterType = source.fClusterType;
160 if(source.fTracksMatched){
161 // assign or copy construct
163 *fTracksMatched = *source.fTracksMatched;
165 else fTracksMatched = new TArrayI(*source.fTracksMatched);
168 if(fTracksMatched)delete fTracksMatched;
173 // assign or copy construct
175 *fLabels = *source.fLabels;
177 else fLabels = new TArrayI(*source.fLabels);
180 if(fLabels)delete fLabels;
185 if(source.fDigitAmplitude){
186 // assign or copy construct
187 if(fDigitAmplitude) *fDigitAmplitude = *source.fDigitAmplitude;
188 else fDigitAmplitude = new TArrayS(*source.fDigitAmplitude);
191 delete fDigitAmplitude;
197 if(source.fDigitTime){
198 // assign or copy construct
199 if(fDigitTime) *fDigitTime = *source.fDigitTime;
200 else fDigitTime = new TArrayS(*source.fDigitTime);
209 if(source.fDigitIndex){
210 // assign or copy construct
211 if(fDigitIndex) *fDigitIndex = *source.fDigitIndex;
212 else fDigitIndex = new TArrayS(*source.fDigitIndex);
223 void AliESDCaloCluster::Copy(TObject &obj) const {
225 // this overwrites the virtual TOBject::Copy()
226 // to allow run time copying without casting
229 if(this==&obj)return;
230 AliESDCaloCluster *robj = dynamic_cast<AliESDCaloCluster*>(&obj);
231 if(!robj)return; // not an AliESDCluster
236 //_______________________________________________________________________
237 AliESDCaloCluster::~AliESDCaloCluster(){
239 // This is destructor according Coding Conventions
241 if(fTracksMatched)delete fTracksMatched;fTracksMatched = 0;
242 if(fLabels) delete fLabels; fLabels = 0;
243 delete fDigitAmplitude; //not in use
244 delete fDigitTime; //not in use
245 delete fDigitIndex; //not in use
246 if(fCellsAmpFraction){ delete[] fCellsAmpFraction; fCellsAmpFraction=0;}
247 if(fCellsAbsId){ delete[] fCellsAbsId; fCellsAbsId = 0;}
250 //_______________________________________________________________________
251 void AliESDCaloCluster::SetPid(const Float_t *p) {
252 // Sets the probability of each particle type
253 // Copied from AliESDtrack SetPIDValues
254 // This function copies "n" PID weights from "scr" to "dest"
255 // and normalizes their sum to 1 thus producing conditional
257 // The negative weights are set to 0.
258 // In case all the weights are non-positive they are replaced by
259 // uniform probabilities
261 Int_t n = AliPID::kSPECIESN;
263 Float_t uniform = 1./(Float_t)n;
266 for (Int_t i=0; i<n; i++)
276 for (Int_t i=0; i<n; i++) fPID[i] /= sum;
278 for (Int_t i=0; i<n; i++) fPID[i] = uniform;
282 //_______________________________________________________________________
283 void AliESDCaloCluster::GetMomentum(TLorentzVector& p, Double_t *vertex ) {
284 // Returns TLorentzVector with momentum of the cluster. Only valid for clusters
285 // identified as photons or pi0 (overlapped gamma) produced on the vertex
286 //Vertex can be recovered with esd pointer doing:
287 //" Double_t vertex[3] ; esd->GetVertex()->GetXYZ(vertex) ; "
289 if(vertex){//calculate direction from vertex
290 fGlobalPos[0]-=vertex[0];
291 fGlobalPos[1]-=vertex[1];
292 fGlobalPos[2]-=vertex[2];
295 Double_t r = TMath::Sqrt(fGlobalPos[0]*fGlobalPos[0]+
296 fGlobalPos[1]*fGlobalPos[1]+
297 fGlobalPos[2]*fGlobalPos[2] ) ;
299 p.SetPxPyPzE( fEnergy*fGlobalPos[0]/r, fEnergy*fGlobalPos[1]/r, fEnergy*fGlobalPos[2]/r, fEnergy) ;
303 //_______________________________________________________________________
304 void AliESDCaloCluster::SetCellsAbsId(UShort_t *array)
306 // Set the array of cell absId numbers
308 fCellsAbsId = new UShort_t[fNCells];
309 for (Int_t i = 0; i < fNCells; i++) fCellsAbsId[i] = array[i];
313 //_______________________________________________________________________
314 void AliESDCaloCluster::SetCellsAmplitudeFraction(Double32_t *array)
316 // Set the array of cell amplitude fraction
318 fCellsAmpFraction = new Double32_t[fNCells];
319 for (Int_t i = 0; i < fNCells; i++) fCellsAmpFraction[i] = array[i];