Update in LambdaC cut class (Rossella)

diff --git a/PWG3/vertexingHF/AliRDHFCutsLctopKpi.cxx b/PWG3/vertexingHF/AliRDHFCutsLctopKpi.cxx
index 75323252bcfa7aa56a831e019632920c24e2bec9..be90913e034a2a1c0444cbbca9402041218d644d 100644 (file)
--- a/PWG3/vertexingHF/AliRDHFCutsLctopKpi.cxx
+++ b/PWG3/vertexingHF/AliRDHFCutsLctopKpi.cxx
@@ -40,7 +40,8 @@ AliRDHFCutsLctopKpi::AliRDHFCutsLctopKpi(const char* name) :
 AliRDHFCuts(name),
 fPidObjprot(0),
 fPidObjpion(0),
-fRecoKF(kFALSE)
+fRecoKF(kFALSE),
+fUseImpParProdCorrCut(kFALSE)
 {
   //
   // Default Constructor
@@ -97,7 +98,8 @@ AliRDHFCutsLctopKpi::AliRDHFCutsLctopKpi(const AliRDHFCutsLctopKpi &source) :
   AliRDHFCuts(source),
   fPidObjprot(0),
   fPidObjpion(0),
-  fRecoKF(kFALSE)
+  fRecoKF(kFALSE),
+  fUseImpParProdCorrCut(source.fUseImpParProdCorrCut)
 {
   //
   // Copy constructor
@@ -307,8 +309,10 @@ Int_t AliRDHFCutsLctopKpi::IsSelected(TObject* obj,Int_t selectionLevel,AliAODEv
     if(d->GetDist12toPrim()<fCutsRD[GetGlobalIndex(5,ptbin)]|| d->GetDist23toPrim()<fCutsRD[GetGlobalIndex(5,ptbin)]) return 0;
     if(d->GetDist12toPrim()>1.) return 0;
     if(d->GetDist23toPrim()>1.) return 0;
-    if(d->Getd0Prong(0)*d->Getd0Prong(1)<0. && d->Getd0Prong(2)*d->Getd0Prong(1)<0.) return 0;
-    
+    if(fUseImpParProdCorrCut){
+      if(d->Getd0Prong(0)*d->Getd0Prong(1)<0. && d->Getd0Prong(2)*d->Getd0Prong(1)<0.) return 0;
+    }
+
     //sec vert
     if(d->GetSigmaVert()>fCutsRD[GetGlobalIndex(6,ptbin)]) return 0;
 

diff --git a/PWG3/vertexingHF/AliRDHFCutsLctopKpi.h b/PWG3/vertexingHF/AliRDHFCutsLctopKpi.h
index 8d8d1313a8bf489874008dc4c5de2066c47db7a2..44144abaa1b6df696093de046c16d889a006d1a7 100644 (file)
--- a/PWG3/vertexingHF/AliRDHFCutsLctopKpi.h
+++ b/PWG3/vertexingHF/AliRDHFCutsLctopKpi.h
@@ -60,13 +60,21 @@ class AliRDHFCutsLctopKpi : public AliRDHFCuts
   Float_t GetMassCut(Int_t iPtBin=0) const { return (GetCuts() ? fCutsRD[GetGlobalIndex(0,iPtBin)] : 1.e6);}
   Float_t GetDCACut(Int_t iPtBin=0) const { return (GetCuts() ? fCutsRD[GetGlobalIndex(11,iPtBin)] : 1.e6);}
 
+  void SetUseImpParProdCorrCut(Bool_t use){
+    fUseImpParProdCorrCut=use;
+  }
+  Bool_t GetUseImpParProdCorrCut() const {
+    return fUseImpParProdCorrCut;
+  }
+
   Bool_t ReconstructKF(AliAODRecoDecayHF3Prong *d,Int_t *pdgs,Double_t field) const;
  protected:
   AliAODPidHF *fPidObjprot;
   AliAODPidHF *fPidObjpion;
   Bool_t fRecoKF;
+  Bool_t fUseImpParProdCorrCut; //switch for cut on d0p*d0K vs. d0K*d0pi 
 
-  ClassDef(AliRDHFCutsLctopKpi,3);  // class for cuts on AOD reconstructed Lc->pKpi
+  ClassDef(AliRDHFCutsLctopKpi,4);  // class for cuts on AOD reconstructed Lc->pKpi
 };
 
 #endif