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 **************************************************************************/
18 /////////////////////////////////////////////////////////////////////////////
20 // AliTOFtrack class //
22 // Authors: Bologna-CERN-ITEP-Salerno Group //
24 // Description: class for handling ESD extracted tracks for TOF matching. //
26 /////////////////////////////////////////////////////////////////////////////
28 #include "AliESDtrack.h"
29 #include "AliTracker.h"
31 #include "AliTOFGeometry.h"
32 #include "AliTOFtrack.h"
36 //_____________________________________________________________________________
37 AliTOFtrack::AliTOFtrack() :
43 // Default constructor.
47 //_____________________________________________________________________________
48 AliTOFtrack::AliTOFtrack(const AliTOFtrack& t) :
58 //_____________________________________________________________________________
59 AliTOFtrack::AliTOFtrack(const AliESDtrack& t) :
65 // Constructor from AliESDtrack
67 SetLabel(t.GetLabel());
71 Set(t.GetX(),t.GetAlpha(),t.GetParameter(),t.GetCovariance());
73 if ((t.GetStatus()&AliESDtrack::kTIME) == 0) return;
75 Double_t times[10]; t.GetIntegratedTimes(times); SetIntegratedTimes(times);
76 SetIntegratedLength(t.GetIntegratedLength());
80 //____________________________________________________________________________
81 AliTOFtrack& AliTOFtrack::operator=(const AliTOFtrack &/*source*/)
89 //_____________________________________________________________________________
90 Bool_t AliTOFtrack::PropagateTo(Double_t xk,Double_t /*x0*/,Double_t /*rho*/)
93 // Propagates a track of particle with mass=pm to a reference plane
94 // defined by x=xk through media of density=rho and radiationLength=x0
97 if (xk == GetX()) return kTRUE;
99 Double_t oldX=GetX();//, oldY=GetY(), oldZ=GetZ();
100 Double_t start[3], end[3], mparam[7];
102 /* get start position */
105 /* propagate the track */
106 Double_t b[3];GetBxByBz(b);
107 if (!AliExternalTrackParam::PropagateToBxByBz(xk,b)) return kFALSE;
109 //Double_t bz=GetBz();
110 //if (!AliExternalTrackParam::PropagateTo(xk,bz)) return kFALSE;
112 /* get end position */
115 /* add time step to integral */
117 Double_t d = TMath::Sqrt((GetX()-oldX)*(GetX()-oldX) +
118 (GetY()-oldY)*(GetY()-oldY) +
119 (GetZ()-oldZ)*(GetZ()-oldZ));
120 if (IsStartedTimeIntegral() && GetX()>oldX) AddTimeStep(d);
122 Double_t d = TMath::Sqrt((end[0]-start[0])*(end[0]-start[0]) +
123 (end[1]-start[1])*(end[1]-start[1]) +
124 (end[2]-start[2])*(end[2]-start[2]));
125 if (IsStartedTimeIntegral() && GetX()>oldX) AddTimeStep(d);
127 /* get material budget from tracker */
128 AliTracker::MeanMaterialBudget(start, end, mparam);
129 Double_t xTimesRho = mparam[4]*mparam[0];
130 Double_t xOverX0 = mparam[1];
132 /* correct for mean material */
133 if (!AliExternalTrackParam::CorrectForMeanMaterial(xOverX0,xTimesRho,GetMass())) return kFALSE;
137 if (!AliExternalTrackParam::CorrectForMaterial(d*rho/x0,x0,GetMass()))
142 //Energy losses************************
143 if((5940*beta2/(1-beta2+1e-10) - beta2) < 0){return 0;}
145 Double_t dE=0.153e-3/beta2*(log(5940*beta2/(1-beta2+1e-10)) - beta2)*d*rho;
147 // suspicious part - think about it ?
148 Double_t kinE = TMath::Sqrt(p2);
149 if (dE>0.8*kinE) dE = 0.8*kinE; //
150 if (dE<0) dE = 0.0; // not valid region for Bethe bloch
156 //_____________________________________________________________________________
157 Bool_t AliTOFtrack::PropagateToInnerTOF()
159 // Propagates a track of particle with mass=pm to a reference plane
160 // defined by x=xk through media of density=rho and radiationLength=x0
163 Double_t ymax=AliTOFGeometry::RinTOF()*TMath::Tan(0.5*AliTOFGeometry::GetAlpha());
164 Bool_t skip = kFALSE;
165 Double_t y=GetYat(AliTOFGeometry::RinTOF(),skip);
170 if (!Rotate(AliTOFGeometry::GetAlpha())) {
173 } else if (y <-ymax) {
174 if (!Rotate(-AliTOFGeometry::GetAlpha())) {
180 Int_t nsteps=Int_t((AliTOFGeometry::Rmin()-x)/0.5); // 0.5 cm Steps
181 for (Int_t istep=0;istep<nsteps;istep++){
182 Float_t xp = x+istep*0.5;
183 // GetPropagationParameters(param);
184 PropagateTo(xp,0.,0.); /* passing 0.,0. as arguments since now
185 this method queries TGeo for material budget
190 if(!PropagateTo(AliTOFGeometry::RinTOF()))return 0;
196 //_________________________________________________________________________
197 Double_t AliTOFtrack::GetPredictedChi2(const AliCluster3D *c) const {
199 // Estimate the chi2 of the space point "c" with its cov. matrix elements
202 Double_t p[3]={c->GetX(), c->GetY(), c->GetZ()};
203 Double_t covyz[3]={c->GetSigmaY2(), c->GetSigmaYZ(), c->GetSigmaZ2()};
204 Double_t covxyz[3]={c->GetSigmaX2(), c->GetSigmaXY(), c->GetSigmaXZ()};
205 return AliExternalTrackParam::GetPredictedChi2(p, covyz, covxyz);
207 //_________________________________________________________________________
208 Bool_t AliTOFtrack::PropagateTo(const AliCluster3D *c) {
210 // Propagates a track to the plane containing cluster 'c'
212 Double_t oldX=GetX(), oldY=GetY(), oldZ=GetZ();
213 Double_t p[3]={c->GetX(), c->GetY(), c->GetZ()};
214 Double_t covyz[3]={c->GetSigmaY2(), c->GetSigmaYZ(), c->GetSigmaZ2()};
215 Double_t covxyz[3]={c->GetSigmaX2(), c->GetSigmaXY(), c->GetSigmaXZ()};
217 if (!AliExternalTrackParam::PropagateTo(p, covyz, covxyz, bz)) return kFALSE;
218 if (IsStartedTimeIntegral()) {
219 Double_t d = TMath::Sqrt((GetX()-oldX)*(GetX()-oldX) +
220 (GetY()-oldY)*(GetY()-oldY) +
221 (GetZ()-oldZ)*(GetZ()-oldZ));
222 if (GetX()<oldX) d=-d;
227 //_________________________________________________________________________
228 Double_t AliTOFtrack::GetYat(Double_t xk, Bool_t & skip) const {
229 //-----------------------------------------------------------------
230 // This function calculates the Y-coordinate of a track at the plane x=xk.
231 // Needed for matching with the TOF (I.Belikov)
232 //-----------------------------------------------------------------
234 skip=(!GetYAt(xk,GetBz(),y));
238 //_____________________________________________________________________________
239 Int_t AliTOFtrack::Compare(const TObject *o) const {
240 //-----------------------------------------------------------------
241 // This function compares tracks according to the their curvature
242 //-----------------------------------------------------------------
243 AliTOFtrack *t=(AliTOFtrack*)o;
244 Double_t co=t->GetSigmaY2()*t->GetSigmaZ2();
245 Double_t c =GetSigmaY2()*GetSigmaZ2();
247 else if (c<co) return -1;
252 //_____________________________________________________________________________
253 void AliTOFtrack::GetPropagationParameters(Double_t *param) {
255 //Get average medium density, x0 while propagating the track
257 //For TRD holes description
259 Double_t thetamin = (90.-31.1) * TMath::Pi()/180.;
260 Double_t thetamax = (90.+31.1) * TMath::Pi()/180.;
262 Double_t zmin = -55.;
266 // Detector inner/outer radii
267 Double_t rTPC = 261.53;
268 Double_t rTPCTRD = 294.5;
269 Double_t rTRD = 369.1;
272 Double_t x0TPC = 40.;
273 Double_t rhoTPC =0.06124;
275 Double_t x0Air = 36.66;
276 Double_t rhoAir =1.2931e-3;
278 Double_t x0TRD = 171.7;
279 Double_t rhoTRD =0.33;
281 // Int_t isec = GetSector();
282 Double_t r[3]; GetXYZ(r);
283 // Float_t thetatr = TMath::ATan2(TMath::Sqrt(r[0]*r[0]+r[1]*r[1]),r[2]);
287 if (isec == 0 || isec == 1 || isec == 2 ) {
288 if( thetatr>=thetamin && thetatr<=thetamax){
293 if (isec == 11 || isec == 12 || isec == 13 || isec == 14 || isec == 15 ) {
294 if( r[2]>=zmin && r[2]<=zmax){
302 {param[0]=x0TPC;param[1]=rhoTPC;}
303 else if(GetX() > rTPC && GetX() < rTPCTRD)
304 {param[0]=x0Air;param[1]=rhoAir;}
305 else if(GetX() >= rTPCTRD && GetX() < rTRD)
306 {param[0]=x0TRD;param[1]=rhoTRD;}
307 else if(GetX() >= rTRD )
308 {param[0]=x0Air;param[1]=rhoAir;}