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 //-----------------------------------------------------------------
19 // Implementation of the combined PID class
20 // For the Event Summary Data Class
21 // produced by the reconstruction process
22 // and containing information on the particle identification
23 // Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch
24 //-----------------------------------------------------------------
26 #include "AliESDpid.h"
27 #include "AliESDEvent.h"
28 #include "AliESDtrack.h"
32 //_________________________________________________________________________
33 Int_t AliESDpid::MakePID(AliESDEvent *event)
36 // Combine the information of various detectors
37 // to determine the Particle Identification
39 Int_t ntrk=event->GetNumberOfTracks();
40 for (Int_t i=0; i<ntrk; i++) {
41 Int_t ns=AliPID::kSPECIES;
42 Double_t p[10]={1.,1.,1.,1.,1.,1.,1.,1.,1.,1.};
44 AliESDtrack *t=event->GetTrack(i);
46 if (t->IsOn(AliESDtrack::kITSpid)) {
49 for (Int_t j=0; j<ns; j++) p[j]*=d[j];
52 if (t->IsOn(AliESDtrack::kTPCpid)) {
55 for (Int_t j=0; j<ns; j++) p[j]*=d[j];
58 if (t->IsOn(AliESDtrack::kTRDpid)) {
61 for (Int_t j=0; j<ns; j++) p[j]*=d[j];
64 if (t->IsOn(AliESDtrack::kTOFpid)) {
67 for (Int_t j=0; j<ns; j++) p[j]*=d[j];
70 if (t->IsOn(AliESDtrack::kHMPIDpid)) {
73 for (Int_t j=0; j<ns; j++) p[j]*=d[j];