1 #ifndef ALIAODPIDUTIL_H
2 #define ALIAODPIDUTIL_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
6 /* $Id: AliAODpidUtil.h 38493 2010-01-26 16:33:03Z hristov $ */
8 //-------------------------------------------------------
11 // Origin: Rosa Romita, GSI, r.romita@gsi.de
12 // Modified: Jens Wiechula, Uni Tuebingen, jens.wiechula@cern.ch
13 //-------------------------------------------------------
16 #include "AliAODTrack.h" // Needed for inline functions
17 #include "AliAODPid.h" // Needed for inline functions
18 //#include "HMPID/AliHMPID.h"
20 #include "AliPIDResponse.h"
25 class AliAODpidUtil : public AliPIDResponse {
28 AliAODpidUtil(Bool_t isMC = kFALSE): AliPIDResponse(isMC), fRange(5.) {;}
29 virtual ~AliAODpidUtil() {;}
31 Int_t MakePID(AliAODTrack *track,Double_t *p) const;
32 void MakeTPCPID(const AliAODTrack *track,Double_t *p) const;
33 void MakeITSPID(const AliAODTrack *track,Double_t *p) const;
34 void MakeTOFPID(const AliAODTrack *track,Double_t *p) const;
35 // void MakeHMPIDPID(AliESDtrack *track);
36 void MakeTRDPID(const AliAODTrack *track,Double_t *p) const;
38 virtual Float_t NumberOfSigmasTOF(const AliVParticle *vtrack, AliPID::EParticleType type) const;
41 Float_t fRange; // nSigma max in likelihood
43 ClassDef(AliAODpidUtil,2) // PID calculation class
46 inline Float_t AliAODpidUtil::NumberOfSigmasTOF(const AliVParticle *vtrack, AliPID::EParticleType type) const {
47 AliAODTrack *track=(AliAODTrack*)vtrack;
48 Double_t times[AliPID::kSPECIES];
49 Double_t sigmaTOFPid[AliPID::kSPECIES];
50 AliAODPid *pidObj = track->GetDetPid();
51 if (!pidObj) return -999.;
52 pidObj->GetIntegratedTimes(times);
53 pidObj->GetTOFpidResolution(sigmaTOFPid);
54 if (sigmaTOFPid[type]>0) return (pidObj->GetTOFsignal() - times[type])/sigmaTOFPid[type];
55 else return (pidObj->GetTOFsignal() - times[type])/fTOFResponse.GetExpectedSigma(track->P(),times[type],AliPID::ParticleMass(type));