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
32 AliRDHFCutsLctopKpi(const char* name="CutsLctopKpi");
34 virtual ~AliRDHFCutsLctopKpi();
36 AliRDHFCutsLctopKpi(const AliRDHFCutsLctopKpi& source);
37 AliRDHFCutsLctopKpi& operator=(const AliRDHFCutsLctopKpi& source);
39 using AliRDHFCuts::GetCutVarsForOpt;
40 virtual void GetCutVarsForOpt(AliAODRecoDecayHF *d,Float_t *vars,Int_t nvars,Int_t *pdgdaughters){
41 return GetCutVarsForOpt(d,vars,nvars,pdgdaughters,0x0);
43 virtual void GetCutVarsForOpt(AliAODRecoDecayHF *d,Float_t *vars,Int_t nvars,Int_t *pdgdaughters,AliAODEvent *aod);
45 void SetPidpion(AliAODPidHF* pidPion) {
46 if(fPidObjpion) delete fPidObjpion;
47 fPidObjpion=new AliAODPidHF(*pidPion);
49 void SetPidprot(AliAODPidHF* pidProt) {
50 if(fPidObjprot) delete fPidObjprot;
51 fPidObjprot=new AliAODPidHF(*pidProt);
54 virtual void SetStandardCutsPP2010();
55 virtual void SetStandardCutsPbPb2010();
58 AliAODPidHF* GetPidpion() const {return fPidObjpion;}
59 AliAODPidHF* GetPidprot() const {return fPidObjprot;}
60 void SetPIDStrategy(EPIDStrategy pidStrategy) {
61 fPIDStrategy=pidStrategy;
63 EPIDStrategy GetPIDStrategy() const {
66 void SetCutsStrategy(ECutsStrategy cutsStrategy) {
67 fCutsStrategy=cutsStrategy;
69 ECutsStrategy GetCutsStrategy() const {
72 void SetPIDThreshold(AliPID::EParticleType species,Double_t threshold) {
73 fPIDThreshold[static_cast<Int_t>(species)]=threshold;
75 Double_t GetPIDThreshold(AliPID::EParticleType species) const {
76 return fPIDThreshold[static_cast<Int_t>(species)];
80 using AliRDHFCuts::IsSelected;
81 virtual Int_t IsSelected(TObject* obj,Int_t selectionLevel)
82 {return IsSelected(obj,selectionLevel,0);}
83 virtual Int_t IsSelected(TObject* obj,Int_t selectionLevel,AliAODEvent *aod);
84 using AliRDHFCuts::IsSelectedPID;
85 virtual Int_t IsSelectedPID(AliAODRecoDecayHF* obj);
86 Int_t IsSelectedCombinedPID(AliAODRecoDecayHF* obj);
87 Int_t CombinePIDCuts (Int_t returnvalue, Int_t returnvaluePID) const;
89 virtual Bool_t IsInFiducialAcceptance(Double_t pt,Double_t y) const;
91 Float_t GetMassCut(Int_t iPtBin=0) const { return (GetCuts() ? fCutsRD[GetGlobalIndex(0,iPtBin)] : 1.e6);}
92 Float_t GetDCACut(Int_t iPtBin=0) const { return (GetCuts() ? fCutsRD[GetGlobalIndex(11,iPtBin)] : 1.e6);}
94 void SetUseImpParProdCorrCut(Bool_t use){
95 fUseImpParProdCorrCut=use;
97 Bool_t GetUseImpParProdCorrCut() const {
98 return fUseImpParProdCorrCut;
101 AliKFParticle* ReconstructKF(AliAODRecoDecayHF3Prong *d,Int_t *pdgs,Double_t field,Bool_t constraint) const;
103 AliAODPidHF *fPidObjprot;
104 AliAODPidHF *fPidObjpion;
105 Bool_t fUseImpParProdCorrCut; //switch for cut on d0p*d0K vs. d0K*d0pi
108 EPIDStrategy fPIDStrategy; // PIS strategy (nsigma, combined)
109 Double_t fPIDThreshold[AliPID::kSPECIES]; // PID threshold for each species
110 ECutsStrategy fCutsStrategy; // cut strategy (standard or KF)
112 ClassDef(AliRDHFCutsLctopKpi,5); // class for cuts on AOD reconstructed Lc->pKpi