1 #ifndef ALIRDHFCUTSLCTOPKPI_H
2 #define ALIRDHFCUTSLCTOPKPI_H
3 /* Copyright(c) 1998-2010, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
8 //***********************************************************
9 // Class AliRDHFCutsLctopKpi
10 // class for cuts on AOD reconstructed Lc->pKpi
11 // Author: A.Dainese, andrea.dainese@pd.infn.it
12 //***********************************************************
14 #include "AliRDHFCuts.h"
15 #include "AliAODPidHF.h"
16 #include "AliAODRecoDecayHF3Prong.h"
18 class AliRDHFCutsLctopKpi : public AliRDHFCuts
34 AliRDHFCutsLctopKpi(const char* name="CutsLctopKpi");
36 virtual ~AliRDHFCutsLctopKpi();
38 AliRDHFCutsLctopKpi(const AliRDHFCutsLctopKpi& source);
39 AliRDHFCutsLctopKpi& operator=(const AliRDHFCutsLctopKpi& source);
41 using AliRDHFCuts::GetCutVarsForOpt;
42 virtual void GetCutVarsForOpt(AliAODRecoDecayHF *d,Float_t *vars,Int_t nvars,Int_t *pdgdaughters){
43 return GetCutVarsForOpt(d,vars,nvars,pdgdaughters,0x0);
45 virtual void GetCutVarsForOpt(AliAODRecoDecayHF *d,Float_t *vars,Int_t nvars,Int_t *pdgdaughters,AliAODEvent *aod);
47 void SetPidpion(AliAODPidHF* pidPion) {
48 if(fPidObjpion) delete fPidObjpion;
49 fPidObjpion=new AliAODPidHF(*pidPion);
51 void SetPidprot(AliAODPidHF* pidProt) {
52 if(fPidObjprot) delete fPidObjprot;
53 fPidObjprot=new AliAODPidHF(*pidProt);
56 virtual void SetStandardCutsPP2010();
57 virtual void SetStandardCutsPbPb2010();
58 virtual void SetStandardCutsPbPb2011();
61 AliAODPidHF* GetPidpion() const {return fPidObjpion;}
62 AliAODPidHF* GetPidprot() const {return fPidObjprot;}
63 void SetPIDStrategy(EPIDStrategy pidStrategy) {
64 fPIDStrategy=pidStrategy;
66 EPIDStrategy GetPIDStrategy() const {
69 void SetCutsStrategy(ECutsStrategy cutsStrategy) {
70 fCutsStrategy=cutsStrategy;
72 ECutsStrategy GetCutsStrategy() const {
75 void SetPIDThreshold(AliPID::EParticleType species,Double_t threshold) {
76 fPIDThreshold[static_cast<Int_t>(species)]=threshold;
78 Double_t GetPIDThreshold(AliPID::EParticleType species) const {
79 return fPIDThreshold[static_cast<Int_t>(species)];
81 Bool_t GetUseSpecialCut(){return fUseSpecialCut;}
82 void SetUseSpecialCut(Bool_t useSpecialCut=kTRUE){fUseSpecialCut=useSpecialCut;}
84 using AliRDHFCuts::IsSelected;
85 virtual Int_t IsSelected(TObject* obj,Int_t selectionLevel)
86 {return IsSelected(obj,selectionLevel,0);}
87 virtual Int_t IsSelected(TObject* obj,Int_t selectionLevel,AliAODEvent *aod);
88 using AliRDHFCuts::IsSelectedPID;
89 virtual Int_t IsSelectedPID(AliAODRecoDecayHF* obj);
90 Int_t IsSelectedCombinedPID(AliAODRecoDecayHF* obj);
91 Int_t IsSelectedCombinedPIDSoft(AliAODRecoDecayHF* obj);
92 Int_t IsSelectedPIDStrong(AliAODRecoDecayHF* obj);
93 Int_t CombinePIDCuts (Int_t returnvalue, Int_t returnvaluePID) const;
95 virtual Bool_t IsInFiducialAcceptance(Double_t pt,Double_t y) const;
97 Float_t GetMassCut(Int_t iPtBin=0) const { return (GetCuts() ? fCutsRD[GetGlobalIndex(0,iPtBin)] : 1.e6);}
98 Float_t GetDCACut(Int_t iPtBin=0) const { return (GetCuts() ? fCutsRD[GetGlobalIndex(11,iPtBin)] : 1.e6);}
100 void SetUseImpParProdCorrCut(Bool_t use){
101 fUseImpParProdCorrCut=use;
103 Bool_t GetUseImpParProdCorrCut() const {
104 return fUseImpParProdCorrCut;
107 AliKFParticle* ReconstructKF(AliAODRecoDecayHF3Prong *d,Int_t *pdgs,Double_t field,Bool_t constraint) const;
109 AliAODPidHF *fPidObjprot;
110 AliAODPidHF *fPidObjpion;
111 Bool_t fUseImpParProdCorrCut; //switch for cut on d0p*d0K vs. d0K*d0pi
114 EPIDStrategy fPIDStrategy; // PIS strategy (nsigma, combined)
115 Double_t fPIDThreshold[AliPID::kSPECIES]; // PID threshold for each species
116 ECutsStrategy fCutsStrategy; // cut strategy (standard or KF)
117 Bool_t fUseSpecialCut;
119 ClassDef(AliRDHFCutsLctopKpi,6); // class for cuts on AOD reconstructed Lc->pKpi