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 **************************************************************************/
15 //-----------------------------------------------------------------
16 // Implementation of the ESD track class
17 // ESD = Event Summary Data
18 // This is the class to deal with during the phisical analysis of data
19 // Origin: Iouri Belikov, CERN
20 // e-mail: Jouri.Belikov@cern.ch
21 //-----------------------------------------------------------------
25 #include "AliESDtrack.h"
26 #include "AliKalmanTrack.h"
30 //_______________________________________________________________________
31 AliESDtrack::AliESDtrack() :
52 fTPCClusterMap(159),//number of padrows
66 // The default ESD constructor
68 for (Int_t i=0; i<kSPECIES; i++) {
78 for (Int_t i=0; i<kSPECIESN; i++) {
84 fPHOSpos[0]=fPHOSpos[1]=fPHOSpos[2]=0.;
85 fEMCALpos[0]=fEMCALpos[1]=fEMCALpos[2]=0.;
87 for (i=0; i<5; i++) { fRp[i]=0.; fCp[i]=0.; fIp[i]=0.; fOp[i]=0.; fTp[i]=0.;}
88 for (i=0; i<15; i++) { fRc[i]=0.; fCc[i]=0.; fIc[i]=0.; fOc[i]=0.; fTc[i]=0.; }
89 for (i=0; i<6; i++) { fITSindex[i]=0; }
90 for (i=0; i<180; i++){ fTPCindex[i]=0; }
91 for (i=0; i<130; i++) { fTRDindex[i]=0; }
99 //_______________________________________________________________________
101 AliESDtrack::AliESDtrack(const AliESDtrack& track):TObject(track){
105 fFlags = track.fFlags;
106 fLabel =track.fLabel;
107 fTrackLength =track.fTrackLength;
108 for (Int_t i=0;i<kSPECIES;i++) fTrackTime[i] =track.fTrackTime[i];
109 for (Int_t i=0;i<kSPECIES;i++) fR[i] =track.fR[i];
110 fStopVertex =track.fStopVertex;
112 fRalpha =track.fRalpha;
114 for (Int_t i=0;i<5;i++) fRp[i] =track.fRp[i];
115 for (Int_t i=0;i<15;i++) fRc[i] =track.fRc[i];
117 fCalpha =track.fCalpha;
119 for (Int_t i=0;i<5;i++) fCp[i] =track.fCp[i];
120 for (Int_t i=0;i<15;i++) fCc[i] =track.fCc[i];
121 fCchi2 =track.fCchi2;
123 fIalpha =track.fIalpha;
125 for (Int_t i=0;i<5;i++) fIp[i] =track.fIp[i];
126 for (Int_t i=0;i<15;i++) fIc[i] =track.fIc[i];
128 fTalpha =track.fTalpha;
130 for (Int_t i=0;i<5;i++) fTp[i] =track.fTp[i];
131 for (Int_t i=0;i<15;i++) fTc[i] =track.fTc[i];
133 fOalpha =track.fOalpha;
135 for (Int_t i=0;i<5;i++) fOp[i] =track.fOp[i];
136 for (Int_t i=0;i<15;i++) fOc[i] =track.fOc[i];
138 fXalpha =track.fXalpha;
140 for (Int_t i=0;i<5;i++) fXp[i] =track.fXp[i];
141 for (Int_t i=0;i<15;i++) fXc[i] =track.fXc[i];
143 fITSchi2 =track.fITSchi2;
144 for (Int_t i=0;i<12;i++) fITSchi2MIP[i] =track.fITSchi2MIP[i];
145 fITSncls =track.fITSncls;
146 for (Int_t i=0;i<6;i++) fITSindex[i]=track.fITSindex[i];
147 fITSsignal =track.fITSsignal;
148 for (Int_t i=0;i<kSPECIES;i++) fITSr[i]=track.fITSr[i];
149 fITSLabel =track.fITSLabel;
150 fITSFakeRatio =track.fITSFakeRatio;
151 fITStrack =0; //coping separatelly - in user code
153 fTPCchi2 =track.fTPCchi2;
154 fTPCncls =track.fTPCncls;
155 for (Int_t i=0;i<180;i++) fTPCindex[i]=track.fTPCindex[i];
156 fTPCClusterMap=track.fTPCClusterMap;
157 fTPCsignal=track.fTPCsignal;
158 for (Int_t i=0;i<kSPECIES;i++) fTPCr[i]=track.fTPCr[i];
159 fTPCLabel=track.fTPCLabel;
161 fTRDchi2=track.fTRDchi2;
162 fTRDncls=track.fTRDncls;
163 fTRDncls0=track.fTRDncls0;
164 for (Int_t i=0;i<130;i++) fTRDindex[i]=track.fTRDindex[i];
165 fTRDsignal=track.fTRDsignal;
166 for (Int_t i=0;i<kSPECIES;i++) fTRDr[i]=track.fTRDr[i];
167 fTRDLabel=track.fTRDLabel;
170 fTOFchi2=track.fTOFchi2;
171 fTOFindex=track.fTOFindex;
172 fTOFsignal=track.fTOFsignal;
173 for (Int_t i=0;i<kSPECIES;i++) fTOFr[i]=track.fTOFr[i];
175 for (Int_t i=0;i<3;i++) fPHOSpos[i]=track.fPHOSpos[i];
176 fPHOSsignal=track.fPHOSsignal;
177 for (Int_t i=0;i<kSPECIESN;i++) fPHOSr[i]=track.fPHOSr[i];
179 for (Int_t i=0;i<3;i++) fEMCALpos[i]=track.fEMCALpos[i];
180 fEMCALsignal=track.fEMCALsignal;
181 for (Int_t i=0;i<kSPECIESN;i++) fEMCALr[i]=track.fEMCALr[i];
183 fRICHsignal=track.fRICHsignal;
184 for (Int_t i=0;i<kSPECIES;i++) fRICHr[i]=track.fRICHr[i];
188 //_______________________________________________________________________
189 AliESDtrack::~AliESDtrack(){
191 // This is destructor according Coding Conventrions
193 //printf("Delete track\n");
198 //_______________________________________________________________________
199 Double_t AliESDtrack::GetMass() const {
200 // Returns the mass of the most probable particle type
203 for (Int_t i=0; i<kSPECIES; i++) {
204 if (fR[i]>max) {k=i; max=fR[i];}
206 if (k==0) { // dE/dx "crossing points" in the TPC
208 if ((p>0.38)&&(p<0.48))
209 if (fR[0]<fR[3]*10.) return 0.49368;
210 if ((p>0.75)&&(p<0.85))
211 if (fR[0]<fR[4]*10.) return 0.93827;
214 if (k==1) return 0.10566;
215 if (k==2||k==-1) return 0.13957;
216 if (k==3) return 0.49368;
217 if (k==4) return 0.93827;
218 Warning("GetMass()","Undefined mass !");
222 //_______________________________________________________________________
223 Bool_t AliESDtrack::UpdateTrackParams(AliKalmanTrack *t, ULong_t flags) {
225 // This function updates track's running parameters
228 fLabel=t->GetLabel();
230 if (t->IsStartedTimeIntegral()) {
232 Double_t times[10];t->GetIntegratedTimes(times); SetIntegratedTimes(times);
233 SetIntegratedLength(t->GetIntegratedLength());
236 fRalpha=t->GetAlpha();
237 t->GetExternalParameters(fRx,fRp);
238 t->GetExternalCovariance(fRc);
242 case kITSin: case kITSout: case kITSrefit:
243 fITSncls=t->GetNumberOfClusters();
244 fITSchi2=t->GetChi2();
245 for (Int_t i=0;i<fITSncls;i++) fITSindex[i]=t->GetClusterIndex(i);
246 fITSsignal=t->GetPIDsignal();
247 fITSLabel = t->GetLabel();
248 fITSFakeRatio = t->GetFakeRatio();
251 case kTPCin: case kTPCrefit:
252 fTPCLabel = t->GetLabel();
257 for (i=0; i<5; i++) fIp[i]=fRp[i];
258 for (i=0; i<15;i++) fIc[i]=fRc[i];
262 fTPCncls=t->GetNumberOfClusters();
263 fTPCchi2=t->GetChi2();
265 {//prevrow must be declared in separate namespace, otherwise compiler cries:
266 //"jump to case label crosses initialization of `Int_t prevrow'"
268 // for (Int_t i=0;i<fTPCncls;i++)
269 for (Int_t i=0;i<160;i++)
271 fTPCindex[i]=t->GetClusterIndex(i);
273 // Piotr's Cluster Map for HBT
274 // ### please change accordingly if cluster array is changing
275 // to "New TPC Tracking" style (with gaps in array)
276 Int_t idx = fTPCindex[i];
277 Int_t sect = (idx&0xff000000)>>24;
278 Int_t row = (idx&0x00ff0000)>>16;
279 if (sect > 18) row +=63; //if it is outer sector, add number of inner sectors
281 fTPCClusterMap.SetBitNumber(row,kTRUE);
283 //Fill the gap between previous row and this row with 0 bits
284 //In case ### pleas change it as well - just set bit 0 in case there
285 //is no associated clusters for current "i"
288 prevrow = row;//if previous bit was not assigned yet == this is the first one
291 { //we don't know the order (inner to outer or reverse)
292 //just to be save in case it is going to change
305 for (Int_t j = n+1; j < m; j++)
307 fTPCClusterMap.SetBitNumber(j,kFALSE);
311 // End Of Piotr's Cluster Map for HBT
314 fTPCsignal=t->GetPIDsignal();
315 {Double_t mass=t->GetMass(); // preliminary mass setting
316 if (mass>0.5) fR[4]=1.; // used by
317 else if (mass<0.4) fR[2]=1.; // the ITS reconstruction
323 { //requested by the PHOS/EMCAL ("temporary solution")
325 Double_t x=460.,buf[15];
326 if (t->PropagateTo(x,30.,0.)) {
327 fOalpha=t->GetAlpha();
328 t->GetExternalParameters(x,buf); fOx=x;
329 for (i=0; i<5; i++) fOp[i]=buf[i];
330 t->GetExternalCovariance(buf);
331 for (i=0; i<15; i++) fOc[i]=buf[i];
334 if (t->PropagateTo(x,30.,0.)) {
335 fXalpha=t->GetAlpha();
336 t->GetExternalParameters(x,buf); fXx=x;
337 for (i=0; i<5; i++) fXp[i]=buf[i];
338 t->GetExternalCovariance(buf);
339 for (i=0; i<15; i++) fXc[i]=buf[i];
342 case kTRDin: case kTRDrefit:
343 fTRDLabel = t->GetLabel();
345 fTRDncls=t->GetNumberOfClusters();
346 fTRDchi2=t->GetChi2();
347 for (Int_t i=0;i<fTRDncls;i++) fTRDindex[i]=t->GetClusterIndex(i);
348 fTRDsignal=t->GetPIDsignal();
351 t->GetExternalParameters(fTx,fTp);
352 t->GetExternalCovariance(fTc);
353 fTRDncls0 = t->GetNumberOfClusters();
362 Error("UpdateTrackParams()","Wrong flag !\n");
369 //_______________________________________________________________________
371 AliESDtrack::SetConstrainedTrackParams(AliKalmanTrack *t, Double_t chi2) {
373 // This function sets the constrained track parameters
377 fCalpha=t->GetAlpha();
378 t->GetExternalParameters(x,buf); fCx=x;
379 for (i=0; i<5; i++) fCp[i]=buf[i];
380 t->GetExternalCovariance(buf);
381 for (i=0; i<15; i++) fCc[i]=buf[i];
386 //_______________________________________________________________________
387 void AliESDtrack::GetExternalParameters(Double_t &x, Double_t p[5]) const {
388 //---------------------------------------------------------------------
389 // This function returns external representation of the track parameters
390 //---------------------------------------------------------------------
392 for (Int_t i=0; i<5; i++) p[i]=fRp[i];
394 //_______________________________________________________________________
395 void AliESDtrack::GetExternalCovariance(Double_t cov[15]) const {
396 //---------------------------------------------------------------------
397 // This function returns external representation of the cov. matrix
398 //---------------------------------------------------------------------
399 for (Int_t i=0; i<15; i++) cov[i]=fRc[i];
403 //_______________________________________________________________________
405 AliESDtrack::GetConstrainedExternalParameters(Double_t &x, Double_t p[5])const{
406 //---------------------------------------------------------------------
407 // This function returns the constrained external track parameters
408 //---------------------------------------------------------------------
410 for (Int_t i=0; i<5; i++) p[i]=fCp[i];
412 //_______________________________________________________________________
414 AliESDtrack::GetConstrainedExternalCovariance(Double_t c[15]) const {
415 //---------------------------------------------------------------------
416 // This function returns the constrained external cov. matrix
417 //---------------------------------------------------------------------
418 for (Int_t i=0; i<15; i++) c[i]=fCc[i];
422 Double_t AliESDtrack::GetP() const {
423 //---------------------------------------------------------------------
424 // This function returns the track momentum
425 //---------------------------------------------------------------------
426 if (TMath::Abs(fRp[4])<=0) return 0;
427 Double_t pt=1./TMath::Abs(fRp[4]);
428 return pt*TMath::Sqrt(1.+ fRp[3]*fRp[3]);
431 void AliESDtrack::GetConstrainedPxPyPz(Double_t *p) const {
432 //---------------------------------------------------------------------
433 // This function returns the constrained global track momentum components
434 //---------------------------------------------------------------------
435 if (TMath::Abs(fCp[4])<=0) {
439 Double_t phi=TMath::ASin(fCp[2]) + fCalpha;
440 Double_t pt=1./TMath::Abs(fCp[4]);
441 p[0]=pt*TMath::Cos(phi); p[1]=pt*TMath::Sin(phi); p[2]=pt*fCp[3];
443 void AliESDtrack::GetConstrainedXYZ(Double_t *xyz) const {
444 //---------------------------------------------------------------------
445 // This function returns the global track position
446 //---------------------------------------------------------------------
447 Double_t phi=TMath::ATan2(fCp[0],fCx) + fCalpha;
448 Double_t r=TMath::Sqrt(fCx*fCx + fCp[0]*fCp[0]);
449 xyz[0]=r*TMath::Cos(phi); xyz[1]=r*TMath::Sin(phi); xyz[2]=fCp[1];
452 void AliESDtrack::GetPxPyPz(Double_t *p) const {
453 //---------------------------------------------------------------------
454 // This function returns the global track momentum components
455 //---------------------------------------------------------------------
456 if (TMath::Abs(fRp[4])<=0) {
460 Double_t phi=TMath::ASin(fRp[2]) + fRalpha;
461 Double_t pt=1./TMath::Abs(fRp[4]);
462 p[0]=pt*TMath::Cos(phi); p[1]=pt*TMath::Sin(phi); p[2]=pt*fRp[3];
464 void AliESDtrack::GetXYZ(Double_t *xyz) const {
465 //---------------------------------------------------------------------
466 // This function returns the global track position
467 //---------------------------------------------------------------------
468 Double_t phi=TMath::ATan2(fRp[0],fRx) + fRalpha;
469 Double_t r=TMath::Sqrt(fRx*fRx + fRp[0]*fRp[0]);
470 xyz[0]=r*TMath::Cos(phi); xyz[1]=r*TMath::Sin(phi); xyz[2]=fRp[1];
474 void AliESDtrack::GetInnerPxPyPz(Double_t *p) const {
475 //---------------------------------------------------------------------
476 // This function returns the global track momentum components
477 // af the entrance of the TPC
478 //---------------------------------------------------------------------
479 if (fIx==0) {p[0]=p[1]=p[2]=0.; return;}
480 Double_t phi=TMath::ASin(fIp[2]) + fIalpha;
481 Double_t pt=1./TMath::Abs(fIp[4]);
482 p[0]=pt*TMath::Cos(phi); p[1]=pt*TMath::Sin(phi); p[2]=pt*fIp[3];
485 void AliESDtrack::GetInnerXYZ(Double_t *xyz) const {
486 //---------------------------------------------------------------------
487 // This function returns the global track position
488 // af the entrance of the TPC
489 //---------------------------------------------------------------------
490 if (fIx==0) {xyz[0]=xyz[1]=xyz[2]=0.; return;}
491 Double_t phi=TMath::ATan2(fIp[0],fIx) + fIalpha;
492 Double_t r=TMath::Sqrt(fIx*fIx + fIp[0]*fIp[0]);
493 xyz[0]=r*TMath::Cos(phi); xyz[1]=r*TMath::Sin(phi); xyz[2]=fIp[1];
496 void AliESDtrack::GetInnerExternalParameters(Double_t &x, Double_t p[5]) const
499 //---------------------------------------------------------------------
500 // This function returns external representation of the track parameters at Inner Layer of TPC
501 //---------------------------------------------------------------------
503 for (Int_t i=0; i<5; i++) p[i]=fIp[i];
505 void AliESDtrack::GetInnerExternalCovariance(Double_t cov[15]) const
508 //---------------------------------------------------------------------
509 // This function returns external representation of the cov. matrix at Inner Layer of TPC
510 //---------------------------------------------------------------------
511 for (Int_t i=0; i<15; i++) cov[i]=fIc[i];
515 void AliESDtrack::GetTRDExternalParameters(Double_t &x, Double_t p[5], Double_t cov[15]) const
518 //this function returns TRD parameters
521 for (Int_t i=0; i<5; i++) p[i]=fTp[i];
522 for (Int_t i=0; i<15; i++) cov[i]=fTc[i];
525 void AliESDtrack::GetOuterPxPyPzPHOS(Double_t *p) const {
526 //---------------------------------------------------------------------
527 // This function returns the global track momentum components
528 // af the radius of the PHOS
529 //---------------------------------------------------------------------
533 Double_t phi=TMath::ASin(fOp[2]) + fOalpha;
534 Double_t pt=1./TMath::Abs(fOp[4]);
535 p[0]=pt*TMath::Cos(phi);
536 p[1]=pt*TMath::Sin(phi);
539 void AliESDtrack::GetOuterPxPyPzEMCAL(Double_t *p) const {
540 //---------------------------------------------------------------------
541 // This function returns the global track momentum components
542 // af the radius of the EMCAL
543 //---------------------------------------------------------------------
546 Double_t phi=TMath::ASin(fXp[2]) + fXalpha;
547 Double_t pt=1./TMath::Abs(fXp[4]);
548 p[0]=pt*TMath::Cos(phi);
549 p[1]=pt*TMath::Sin(phi);
553 void AliESDtrack::GetOuterXYZPHOS(Double_t *xyz) const {
554 //---------------------------------------------------------------------
555 // This function returns the global track position
556 // af the radius of the PHOS
557 //---------------------------------------------------------------------
558 xyz[0]=xyz[1]=xyz[2]=0.;
561 Double_t phi=TMath::ATan2(fOp[0],fOx) + fOalpha;
562 Double_t r=TMath::Sqrt(fOx*fOx + fOp[0]*fOp[0]);
563 xyz[0]=r*TMath::Cos(phi); xyz[1]=r*TMath::Sin(phi); xyz[2]=fOp[1];
565 void AliESDtrack::GetOuterXYZEMCAL(Double_t *xyz) const {
566 //---------------------------------------------------------------------
567 // This function returns the global track position
568 // af the radius of the EMCAL
569 //---------------------------------------------------------------------
572 Double_t phi=TMath::ATan2(fXp[0],fOx) + fXalpha;
573 Double_t r=TMath::Sqrt(fXx*fXx + fXp[0]*fXp[0]);
574 xyz[0]=r*TMath::Cos(phi);
575 xyz[1]=r*TMath::Sin(phi);
579 //_______________________________________________________________________
580 void AliESDtrack::GetIntegratedTimes(Double_t *times) const {
581 // Returns the array with integrated times for each particle hypothesis
582 for (Int_t i=0; i<kSPECIES; i++) times[i]=fTrackTime[i];
585 //_______________________________________________________________________
586 void AliESDtrack::SetIntegratedTimes(const Double_t *times) {
587 // Sets the array with integrated times for each particle hypotesis
588 for (Int_t i=0; i<kSPECIES; i++) fTrackTime[i]=times[i];
591 //_______________________________________________________________________
592 void AliESDtrack::SetITSpid(const Double_t *p) {
593 // Sets values for the probability of each particle type (in ITS)
594 for (Int_t i=0; i<kSPECIES; i++) fITSr[i]=p[i];
595 SetStatus(AliESDtrack::kITSpid);
598 void AliESDtrack::SetITSChi2MIP(const Float_t *chi2mip){
599 for (Int_t i=0; i<12; i++) fITSchi2MIP[i]=chi2mip[i];
601 //_______________________________________________________________________
602 void AliESDtrack::GetITSpid(Double_t *p) const {
603 // Gets the probability of each particle type (in ITS)
604 for (Int_t i=0; i<kSPECIES; i++) p[i]=fITSr[i];
607 //_______________________________________________________________________
608 Int_t AliESDtrack::GetITSclusters(UInt_t *idx) const {
609 //---------------------------------------------------------------------
610 // This function returns indices of the assgined ITS clusters
611 //---------------------------------------------------------------------
612 for (Int_t i=0; i<fITSncls; i++) idx[i]=fITSindex[i];
616 //_______________________________________________________________________
617 Int_t AliESDtrack::GetTPCclusters(Int_t *idx) const {
618 //---------------------------------------------------------------------
619 // This function returns indices of the assgined ITS clusters
620 //---------------------------------------------------------------------
622 for (Int_t i=0; i<180; i++) idx[i]=fTPCindex[i]; // MI I prefer some constant
626 //_______________________________________________________________________
627 void AliESDtrack::SetTPCpid(const Double_t *p) {
628 // Sets values for the probability of each particle type (in TPC)
629 for (Int_t i=0; i<kSPECIES; i++) fTPCr[i]=p[i];
630 SetStatus(AliESDtrack::kTPCpid);
633 //_______________________________________________________________________
634 void AliESDtrack::GetTPCpid(Double_t *p) const {
635 // Gets the probability of each particle type (in TPC)
636 for (Int_t i=0; i<kSPECIES; i++) p[i]=fTPCr[i];
639 //_______________________________________________________________________
640 Int_t AliESDtrack::GetTRDclusters(UInt_t *idx) const {
641 //---------------------------------------------------------------------
642 // This function returns indices of the assgined TRD clusters
643 //---------------------------------------------------------------------
645 for (Int_t i=0; i<130; i++) idx[i]=fTRDindex[i]; // MI I prefer some constant
649 //_______________________________________________________________________
650 void AliESDtrack::SetTRDpid(const Double_t *p) {
651 // Sets values for the probability of each particle type (in TRD)
652 for (Int_t i=0; i<kSPECIES; i++) fTRDr[i]=p[i];
653 SetStatus(AliESDtrack::kTRDpid);
656 //_______________________________________________________________________
657 void AliESDtrack::GetTRDpid(Double_t *p) const {
658 // Gets the probability of each particle type (in TRD)
659 for (Int_t i=0; i<kSPECIES; i++) p[i]=fTRDr[i];
662 //_______________________________________________________________________
663 void AliESDtrack::SetTRDpid(Int_t iSpecies, Float_t p)
665 // Sets the probability of particle type iSpecies to p (in TRD)
669 Float_t AliESDtrack::GetTRDpid(Int_t iSpecies) const
671 // Returns the probability of particle type iSpecies (in TRD)
672 return fTRDr[iSpecies];
675 //_______________________________________________________________________
676 void AliESDtrack::SetTOFpid(const Double_t *p) {
677 // Sets the probability of each particle type (in TOF)
678 for (Int_t i=0; i<kSPECIES; i++) fTOFr[i]=p[i];
679 SetStatus(AliESDtrack::kTOFpid);
682 //_______________________________________________________________________
683 void AliESDtrack::GetTOFpid(Double_t *p) const {
684 // Gets probabilities of each particle type (in TOF)
685 for (Int_t i=0; i<kSPECIES; i++) p[i]=fTOFr[i];
690 //_______________________________________________________________________
691 void AliESDtrack::SetPHOSpid(const Double_t *p) {
692 // Sets the probability of each particle type (in PHOS)
693 for (Int_t i=0; i<kSPECIESN; i++) fPHOSr[i]=p[i];
694 SetStatus(AliESDtrack::kPHOSpid);
697 //_______________________________________________________________________
698 void AliESDtrack::GetPHOSpid(Double_t *p) const {
699 // Gets probabilities of each particle type (in PHOS)
700 for (Int_t i=0; i<kSPECIESN; i++) p[i]=fPHOSr[i];
703 //_______________________________________________________________________
704 void AliESDtrack::SetEMCALpid(const Double_t *p) {
705 // Sets the probability of each particle type (in EMCAL)
706 for (Int_t i=0; i<kSPECIESN; i++) fEMCALr[i]=p[i];
707 SetStatus(AliESDtrack::kEMCALpid);
710 //_______________________________________________________________________
711 void AliESDtrack::GetEMCALpid(Double_t *p) const {
712 // Gets probabilities of each particle type (in EMCAL)
713 for (Int_t i=0; i<kSPECIESN; i++) p[i]=fEMCALr[i];
716 //_______________________________________________________________________
717 void AliESDtrack::SetRICHpid(const Double_t *p) {
718 // Sets the probability of each particle type (in RICH)
719 for (Int_t i=0; i<kSPECIES; i++) fRICHr[i]=p[i];
720 SetStatus(AliESDtrack::kRICHpid);
723 //_______________________________________________________________________
724 void AliESDtrack::GetRICHpid(Double_t *p) const {
725 // Gets probabilities of each particle type (in RICH)
726 for (Int_t i=0; i<kSPECIES; i++) p[i]=fRICHr[i];
731 //_______________________________________________________________________
732 void AliESDtrack::SetESDpid(const Double_t *p) {
733 // Sets the probability of each particle type for the ESD track
734 for (Int_t i=0; i<kSPECIES; i++) fR[i]=p[i];
735 SetStatus(AliESDtrack::kESDpid);
738 //_______________________________________________________________________
739 void AliESDtrack::GetESDpid(Double_t *p) const {
740 // Gets probability of each particle type for the ESD track
741 for (Int_t i=0; i<kSPECIES; i++) p[i]=fR[i];
744 //_______________________________________________________________________
745 void AliESDtrack::Print(Option_t *) const {
746 // Prints info on the track
748 Info("Print","Track info") ;
749 Double_t p[kSPECIESN] ;
752 printf("From ITS: ") ;
754 for(index = 0 ; index < kSPECIES; index++)
755 printf("%f, ", p[index]) ;
756 printf("\n signal = %f\n", GetITSsignal()) ;
759 printf("From TPC: ") ;
761 for(index = 0 ; index < kSPECIES; index++)
762 printf("%f, ", p[index]) ;
763 printf("\n signal = %f\n", GetTPCsignal()) ;
766 printf("From TRD: ") ;
768 for(index = 0 ; index < kSPECIES; index++)
769 printf("%f, ", p[index]) ;
770 printf("\n signal = %f\n", GetTRDsignal()) ;
773 printf("From TOF: ") ;
775 for(index = 0 ; index < kSPECIES; index++)
776 printf("%f, ", p[index]) ;
777 printf("\n signal = %f\n", GetTOFsignal()) ;
779 if( IsOn(kRICHpid) ){
780 printf("From TOF: ") ;
782 for(index = 0 ; index < kSPECIES; index++)
783 printf("%f, ", p[index]) ;
784 printf("\n signal = %f\n", GetRICHsignal()) ;
786 if( IsOn(kPHOSpid) ){
787 printf("From PHOS: ") ;
789 for(index = 0 ; index < kSPECIESN; index++)
790 printf("%f, ", p[index]) ;
791 printf("\n signal = %f\n", GetPHOSsignal()) ;
793 if( IsOn(kEMCALpid) ){
794 printf("From EMCAL: ") ;
796 for(index = 0 ; index < kSPECIESN; index++)
797 printf("%f, ", p[index]) ;
798 printf("\n signal = %f\n", GetEMCALsignal()) ;