1 #ifndef ALIRDHFCUTSD0TOKPI_H
2 #define ALIRDHFCUTSD0TOKPI_H
3 /* Copyright(c) 1998-2010, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
8 //***********************************************************
9 // Class AliRDHFCutsD0toKpi
10 // class for cuts on AOD reconstructed D0->Kpi
11 // Author: A.Dainese, andrea.dainese@pd.infn.it
12 //***********************************************************
14 #include "AliRDHFCuts.h"
17 class AliAODRecoDecayHF;
18 class AliAODRecoDecayHF2Prong;
20 class AliRDHFCutsD0toKpi : public AliRDHFCuts
24 AliRDHFCutsD0toKpi(const char* name="CutsD0toKpi");
26 virtual ~AliRDHFCutsD0toKpi(){}
28 AliRDHFCutsD0toKpi(const AliRDHFCutsD0toKpi& source);
29 AliRDHFCutsD0toKpi& operator=(const AliRDHFCutsD0toKpi& source);
31 using AliRDHFCuts::GetCutVarsForOpt;
32 virtual void GetCutVarsForOpt(AliAODRecoDecayHF *d,Float_t *vars,Int_t nvars,Int_t *pdgdaughters){
33 return GetCutVarsForOpt(d,vars,nvars,pdgdaughters,0x0);
35 virtual void GetCutVarsForOpt(AliAODRecoDecayHF *d,Float_t *vars,Int_t nvars,Int_t *pdgdaughters,AliAODEvent *aod);
37 using AliRDHFCuts::IsSelected;
38 virtual Int_t IsSelected(TObject* obj,Int_t selectionLevel)
39 {return IsSelected(obj,selectionLevel,0);}
40 virtual Int_t IsSelected(TObject* obj,Int_t selectionLevel,AliAODEvent* aod);
42 Float_t GetMassCut(Int_t iPtBin=0) const { return (GetCuts() ? fCutsRD[GetGlobalIndex(0,iPtBin)] : 1.e6);}
43 Float_t GetDCACut(Int_t iPtBin=0) const { return (GetCuts() ? fCutsRD[GetGlobalIndex(1,iPtBin)] : 1.e6);}
44 Int_t CombineSelectionLevels(Int_t selectionvalTrack,Int_t selectionvalCand,Int_t selectionvalPID)const;
45 virtual Bool_t IsInFiducialAcceptance(Double_t pt,Double_t y) const;
46 virtual void SetStandardCutsPP2010();
47 virtual void SetStandardCutsPP2011_276TeV();
48 virtual void SetStandardCutsPbPb2010();
49 virtual void SetStandardCutsPbPb2011();
50 void SetStandardCutsPbPb2010Peripherals();
51 virtual Int_t IsSelectedPID(AliAODRecoDecayHF *rd);
52 Int_t IsSelectedPIDdefault(AliAODRecoDecayHF *rd);
53 Int_t IsSelectedSpecialCuts(AliAODRecoDecayHF *d) const;
54 void SetUseSpecialCuts(Bool_t useSpecialCuts) {fUseSpecialCuts=useSpecialCuts;}
55 void SetMaximumPtSpecialCuts(Double_t pt) { fPtMaxSpecialCuts=pt; }
56 void SetMaximumPforPID(Double_t p){fmaxPtrackForPID=p;}
57 Double_t GetMaximumPforPID(){return fmaxPtrackForPID;}
58 Double_t GetMaximumPtSpecialCuts() const { return fPtMaxSpecialCuts; }
59 void SetLowPt(Bool_t lowpt,Double_t ptlow=2.) {fLowPt=lowpt;fPtLowPID=ptlow;}
60 Bool_t GetUseSpecialCuts() const {return fUseSpecialCuts;}
61 void SetUseDefaultPID(Bool_t defPID){fDefaultPID=defPID;}
62 Bool_t GetIsUsedDefPID(){return fDefaultPID;}
63 Double_t GetPtForPIDtight()const {return fPtLowPID;}
64 void SetUseKF(Bool_t useKF);
65 Bool_t GetIsUsedKF() const {return fUseKF;}
69 Int_t IsSelectedKF(AliAODRecoDecayHF2Prong* d,AliAODEvent* aod) const;
71 Bool_t fUseSpecialCuts; // flag to switch on/off special cuts
72 Bool_t fLowPt; // flag to switch on/off different pid for low pt D0
73 Bool_t fDefaultPID; // flag to switch on/off the default pid
74 Bool_t fUseKF; // flag to switch on/off D0 selection via KF
75 Double_t fPtLowPID; // transverse momentum below which the strong PID is applied
76 Double_t fPtMaxSpecialCuts; // transverse momentum below which the special cuts are applied
77 // if set to zero, used for all pt
78 Double_t fmaxPtrackForPID; // max momentum for applying PID
80 ClassDef(AliRDHFCutsD0toKpi,8); // class for cuts on AOD reconstructed D0->Kpi