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 **************************************************************************/
17 // Implements the abstract base class AliHFEpidBase
18 // IsInitialized() does the PID decision
21 // Markus Fasel <M.Fasel@gsi.de>
22 // Matus Kalisky <matus.kalisky@cern.ch> (contact)
27 #include "AliAODTrack.h"
28 #include "AliESDtrack.h"
30 #include "AliPIDResponse.h"
32 #include "AliHFEdetPIDqa.h"
33 #include "AliHFEpidTOF.h"
34 #include "AliHFEpidQAmanager.h"
37 ClassImp(AliHFEpidTOF)
39 //___________________________________________________________________
40 AliHFEpidTOF::AliHFEpidTOF():
49 //___________________________________________________________________
50 AliHFEpidTOF::AliHFEpidTOF(const Char_t *name):
59 //___________________________________________________________________
60 AliHFEpidTOF::AliHFEpidTOF(const AliHFEpidTOF &c):
70 //___________________________________________________________________
71 AliHFEpidTOF &AliHFEpidTOF::operator=(const AliHFEpidTOF &ref){
73 // Assignment operator
82 //___________________________________________________________________
83 AliHFEpidTOF::~AliHFEpidTOF(){
88 //___________________________________________________________________
89 void AliHFEpidTOF::Copy(TObject &ref) const {
91 // Performs the copying of the object
93 AliHFEpidTOF &target = dynamic_cast<AliHFEpidTOF &>(ref);
95 target.fNsigmaTOF = fNsigmaTOF;
97 AliHFEpidBase::Copy(ref);
99 //___________________________________________________________________
100 Bool_t AliHFEpidTOF::InitializePID(Int_t /*run*/){
102 // InitializePID: TOF experts have to implement code here
107 //___________________________________________________________________
108 Int_t AliHFEpidTOF::IsSelected(const AliHFEpidObject *track, AliHFEpidQAmanager *pidqa) const
111 // TOF PID based on n-Sigma cut
112 // Selects Protons and Kaons via n-sigma cut up to 3 GeV/c
113 // In addition histos for n-sigma before (all species) and after (only closest species) are filled
115 if(!fkPIDResponse) return 0;
116 AliDebug(2, "PID object available");
118 const AliVTrack *vtrack = dynamic_cast<const AliVTrack *>(track->GetRecTrack());
119 if(!(vtrack && (vtrack->GetStatus() & AliESDtrack::kTOFpid))) return 0;
120 AliDebug(2, "Track Has TOF PID");
122 if(pidqa) pidqa->ProcessTrack(track, AliHFEpid::kTOFpid, AliHFEdetPIDqa::kBeforePID);
124 // Fill before selection
125 Double_t sigEle = fkPIDResponse->NumberOfSigmasTOF(track->GetRecTrack(), AliPID::kElectron);
126 AliDebug(2, Form("Number of sigmas in TOF: %f", sigEle));
128 if(TMath::Abs(sigEle) < fNsigmaTOF){
130 if(pidqa) pidqa->ProcessTrack(track, AliHFEpid::kTOFpid, AliHFEdetPIDqa::kAfterPID);