1 /**************************************************************************
2 * Copyright(c) 1998-2007, 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 **************************************************************************/
18 //-------------------------------------------------------------------------
19 // AOD Pid class to store additional pid information
20 // Author: Annalisa Mastroserio
21 //-------------------------------------------------------------------------
23 #include "AliAODPid.h"
24 #include "AliESDtrack.h"
26 #include "AliTPCdEdxInfo.h"
31 //______________________________________________________________________________
32 AliAODPid::AliAODPid():
46 // default constructor
47 for(Int_t i=0; i<AliPID::kSPECIES; i++) fIntTime[i] = 0;
48 for(Int_t i=0; i<5; i++) fTOFpidResolution[i] = 0.;
49 for(Int_t i=0; i<6; i++) {
53 for(Int_t i=0; i<4; i++) fITSdEdxSamples[i] = 0.;
56 //______________________________________________________________________________
57 AliAODPid::~AliAODPid()
66 //______________________________________________________________________________
67 AliAODPid::AliAODPid(const AliAODPid& pid) :
69 fITSsignal(pid.fITSsignal),
70 fTPCsignal(pid.fTPCsignal),
71 fTPCsignalN(pid.fTPCsignalN),
72 fTPCmomentum(pid.fTPCmomentum),
74 fTRDnSlices(pid.fTRDnSlices),
75 fTRDntls(pid.fTRDntls),
77 fTRDsignal(pid.fTRDsignal),
78 fTRDChi2(pid.fTRDChi2),
79 fTOFesdsignal(pid.fTOFesdsignal),
83 SetTRDslices(fTRDnSlices, pid.fTRDslices);
84 for(Int_t i=0; i<AliPID::kSPECIES; i++) fIntTime[i]=pid.fIntTime[i];
86 for(Int_t i=0; i<6; i++){
87 fTRDmomentum[i]=pid.fTRDmomentum[i];
91 for(Int_t i=0; i<5; i++) fTOFpidResolution[i]=pid.fTOFpidResolution[i];
93 for(Int_t i=0; i<4; i++) fITSdEdxSamples[i]=pid.fITSdEdxSamples[i];
95 if (pid.fTPCdEdxInfo) fTPCdEdxInfo=new AliTPCdEdxInfo(*pid.fTPCdEdxInfo);
98 //______________________________________________________________________________
99 AliAODPid& AliAODPid::operator=(const AliAODPid& pid)
101 // Assignment operator
104 TObject::operator=(pid);
106 fITSsignal = pid.fITSsignal;
107 for (Int_t i = 0; i < 4; i++) fITSdEdxSamples[i]=pid.fITSdEdxSamples[i];
108 fTPCsignal = pid.fTPCsignal;
109 fTPCsignalN = pid.fTPCsignalN;
110 fTPCmomentum = pid.fTPCmomentum;
111 fTPCTgl = pid.fTPCTgl;
113 fTRDsignal = pid.fTRDsignal;
114 if(fTRDnSlices != pid.fTRDnSlices) {
115 // only delete if number changed or is 0
116 delete [] fTRDslices;
118 fTRDnSlices = pid.fTRDnSlices;
119 if(pid.fTRDnSlices > 0) fTRDslices = new Double32_t[fTRDnSlices];
122 if (fTRDslices && pid.fTRDslices)
123 memcpy(fTRDslices, pid.fTRDslices, fTRDnSlices*sizeof(Double32_t));
125 fTRDntls = pid.fTRDntls;
126 for(Int_t i = 0; i < 6; i++){
127 fTRDmomentum[i] = pid.fTRDmomentum[i];
128 fTRDncls[i] = pid.fTRDncls[i];
131 fTRDChi2 = pid.fTRDChi2;
133 fTOFesdsignal=pid.fTOFesdsignal;
134 for (Int_t i = 0; i < 5; i++) fTOFpidResolution[i]=pid.fTOFpidResolution[i];
135 for (Int_t i = 0; i < 5; i++) fIntTime[i]=pid.fIntTime[i];
137 SetTPCdEdxInfo(pid.fTPCdEdxInfo);
142 //_______________________________________________________________________________
143 void AliAODPid::GetIntegratedTimes(Double_t *timeint, Int_t nspec) const
145 // Returns the array with integrated times for each particle hypothesis
147 for(Int_t i=0; i<AliPID::kSPECIES; i++) timeint[i]=fIntTime[i];
148 //Note: at the moment only kSPECIES entries are available
149 if (nspec>AliPID::kSPECIES) for (int i=AliPID::kSPECIES;i<AliPID::kSPECIESC;i++) timeint[i]=0;
152 //_______________________________________________________________________________
153 void AliAODPid::SetIntegratedTimes(Double_t timeint[AliPID::kSPECIES])
155 // Returns the array with integrated times for each particle hypothesis
156 for(Int_t i=0; i<AliPID::kSPECIES; i++) fIntTime[i]=timeint[i];
158 //______________________________________________________________________________
159 void AliAODPid::SetTOFpidResolution(Double_t tofPIDres[5])
161 for (Int_t i=0; i<5; i++) fTOFpidResolution[i]=tofPIDres[i];
164 //______________________________________________________________________________
165 void AliAODPid::GetTOFpidResolution(Double_t tofRes[5]) const
167 for (Int_t i=0; i<5; i++) tofRes[i]=fTOFpidResolution[i];
170 //______________________________________________________________________________
171 void AliAODPid::SetITSdEdxSamples(const Double_t s[4])
174 // Set the 4 values of dE/dx from individual ITS layers that are read from ESD
176 for (Int_t i=0; i<4; i++) fITSdEdxSamples[i]=s[i];
178 //______________________________________________________________________________
179 void AliAODPid::GetITSdEdxSamples(Double_t s[4]) const
182 // Get the 4 values of dE/dx from individual ITS layers that are read from ESD
184 for (Int_t i=0; i<4; i++) s[i]=fITSdEdxSamples[i];
186 //______________________________________________________________________________
187 void AliAODPid::SetTPCdEdxInfo(AliTPCdEdxInfo * dEdxInfo)
197 if (!fTPCdEdxInfo) fTPCdEdxInfo=new AliTPCdEdxInfo;
198 (*fTPCdEdxInfo)=(*dEdxInfo);