1 #ifndef ALIRDHFCUTSLCTOV0_H
2 #define ALIRDHFCUTSLCTOV0_H
3 /* Copyright(c) 1998-2010, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
8 //***********************************************************
9 // Class AliRDHFCutsLctoV0
10 // class for cuts on AOD reconstructed Lc-> V0 + bachelor
11 //***********************************************************
13 #include "AliAODPidHF.h"
14 #include "AliRDHFCuts.h"
16 class AliRDHFCutsLctoV0 : public AliRDHFCuts
26 enum ELctoV0pidStrategy {
31 AliRDHFCutsLctoV0(const char* name="CutsLctoV0", Short_t v0channel=0);
33 virtual ~AliRDHFCutsLctoV0();
35 AliRDHFCutsLctoV0(const AliRDHFCutsLctoV0& source);
36 AliRDHFCutsLctoV0& operator=(const AliRDHFCutsLctoV0& source);
38 using AliRDHFCuts::GetCutVarsForOpt;
39 virtual void GetCutVarsForOpt(AliAODRecoDecayHF *d,Float_t *vars,Int_t nvars,Int_t *pdgdaughters);
41 using AliRDHFCuts::IsSelected;
42 virtual Int_t IsSelected(TObject* obj,Int_t selectionLevel);
44 using AliRDHFCuts::IsSelectedPID;
45 virtual Int_t IsSelectedPID(AliAODRecoDecayHF* obj);
47 Int_t IsSelectedSingleCut(TObject* obj, Int_t selectionLevel, Int_t cutIndex);
49 Int_t CombinePIDCuts (Int_t returnvalue, Int_t returnvaluePID) const;
51 Float_t GetMassCut(Int_t iPtBin=0) const { return (GetCuts() ? fCutsRD[GetGlobalIndex(0,iPtBin)] : 1.e6);}
52 Float_t GetDCACut(Int_t iPtBin=0) const { return (GetCuts() ? fCutsRD[GetGlobalIndex(7,iPtBin)] : 1.e6);}
54 void SetPidSelectionFlag(Int_t a) {fPidSelectionFlag=a;}
55 Int_t GetPidSelectionFlag() {return fPidSelectionFlag;}
57 virtual void SetStandardCutsPP2010();
58 virtual void SetStandardCutsPbPb2010();
59 virtual void SetStandardCutsPbPb2011();
61 void SetPidV0pos(AliAODPidHF* pidV0pos) {
62 if (fPidHFV0pos) delete fPidHFV0pos;
63 fPidHFV0pos = new AliAODPidHF(*pidV0pos);
65 void SetPidV0neg(AliAODPidHF* pidV0neg) {
66 if (fPidHFV0neg) delete fPidHFV0neg;
67 fPidHFV0neg = new AliAODPidHF(*pidV0neg);
70 AliAODPidHF * GetPidV0pos() { return fPidHFV0pos; }
71 AliAODPidHF * GetPidV0neg() { return fPidHFV0neg; }
73 void AddTrackCutsV0daughters(AliESDtrackCuts* v0daug) {
74 fV0daughtersCuts = new AliESDtrackCuts(*v0daug);
76 virtual AliESDtrackCuts *GetTrackCutsV0daughters() const {return fV0daughtersCuts;}
80 void CheckPID(AliAODTrack *bachelor, AliAODTrack *v0Neg, AliAODTrack *v0Pos,
81 Bool_t &isBachelorID1, Bool_t &isV0NegID2, Bool_t &isV0PosID4);
85 Int_t fPidSelectionFlag;
86 AliAODPidHF *fPidHFV0pos;
87 AliAODPidHF *fPidHFV0neg;
88 AliESDtrackCuts *fV0daughtersCuts; // cuts for v0 daughters (AOD converted to ESD on the flight!)
89 //UShort_t fV0channel;
91 ClassDef(AliRDHFCutsLctoV0,3); // class for cuts on AOD reconstructed Lc->V0+bachelor