Added IsInFiducualAcceptance() (Alessandro)

diff --git a/PWG3/vertexingHF/AliRDHFCutsDStartoKpipi.cxx b/PWG3/vertexingHF/AliRDHFCutsDStartoKpipi.cxx
index 6887747da31cd8a54f9a1330c65427c361953a7c..fd538dcc056cc8879bea6ea0b7230e4b12fc524d 100644 (file)
--- a/PWG3/vertexingHF/AliRDHFCutsDStartoKpipi.cxx
+++ b/PWG3/vertexingHF/AliRDHFCutsDStartoKpipi.cxx
@@ -376,5 +376,28 @@ Int_t AliRDHFCutsDStartoKpipi::IsD0FromDStarSelected(Double_t pt, TObject* obj,I
   //if((Charge()==+1 && !okD0) || (Charge()==-1 && !okD0bar)) return kFALSE; 
   return returnvalue;
 }
+//----------------------------------------------------------------------------------
+Bool_t AliRDHFCutsDStartoKpipi::IsInFiducialAcceptance(Double_t pt, Double_t y) const
+{
+  //
+  // D* fiducial acceptance region 
+  //
 
+  if(pt > 5.) {
+    // applying cut for pt > 5 GeV
+    AliDebug(4,Form("pt of D* = %f (> 5), cutting at |y| < 0.8\n",pt)); 
+    if (TMath::Abs(y) > 0.8){
+      return kFALSE;
+    }
+  } else {
+    // appliying smooth cut for pt < 5 GeV
+    Double_t maxFiducialY = -0.151/15*pt*pt+1.9/15*pt+0.4; 
+    Double_t minFiducialY = 0.151/15*pt*pt-1.9/15*pt-0.4;              
+    AliDebug(4,Form("pt of D* = %f (< 5), cutting  according to the fiducial zone [%f, %f]\n",pt,minFiducialY,maxFiducialY)); 
+    if (y < minFiducialY || y > maxFiducialY){
+      return kFALSE;
+    }
+  }
 
+  return kTRUE;
+}

diff --git a/PWG3/vertexingHF/AliRDHFCutsDStartoKpipi.h b/PWG3/vertexingHF/AliRDHFCutsDStartoKpipi.h
index dc091aea1ccf65f242c3c551ec85246e673bc615..9290992c213ccd23f38be696bd8c19a463dcd7bf 100644 (file)
--- a/PWG3/vertexingHF/AliRDHFCutsDStartoKpipi.h
+++ b/PWG3/vertexingHF/AliRDHFCutsDStartoKpipi.h
@@ -27,7 +27,7 @@ class AliRDHFCutsDStartoKpipi : public AliRDHFCuts
   using AliRDHFCuts::IsSelected;
   virtual Int_t IsSelected(TObject* obj,Int_t selectionLevel);
   Int_t IsD0FromDStarSelected(Double_t pt, TObject* obj,Int_t selectionLevel) const; 
-  
+  virtual Bool_t IsInFiducialAcceptance(Double_t pt,Double_t y) const;
   Float_t GetMassCut(Int_t iPtBin=0) const { return (GetCuts() ? fCutsRD[GetGlobalIndex(9,iPtBin)] : 1.e6);} // for the Dstar
   Float_t GetDCACut(Int_t iPtBin=0) const { return (GetCuts() ? fCutsRD[GetGlobalIndex(1,iPtBin)] : 1.e6);} // for the D0