Added possibility to select on angle between D0 daughters and pi_s (Alessandro)

diff --git a/PWG3/vertexingHF/AliAODRecoCascadeHF.cxx b/PWG3/vertexingHF/AliAODRecoCascadeHF.cxx
index fc38929e083fd9fd1d00fbe2032f533eca4982aa..0608ae4a5dbdcdb817185f57fe48bf569fb4df08 100644 (file)
--- a/PWG3/vertexingHF/AliAODRecoCascadeHF.cxx
+++ b/PWG3/vertexingHF/AliAODRecoCascadeHF.cxx
@@ -275,3 +275,25 @@ Double_t AliAODRecoCascadeHF::AngleD0dkpPisoft() const {
 
   return theta;
 }
+//-----------------------------------------------------------------------------
+Bool_t AliAODRecoCascadeHF::TrigonometricalCut() const {
+  //  
+  // Trigonometrical constraint
+  //
+  TVector3 p3Trk0(Get2Prong()->PxProng(0),Get2Prong()->PyProng(0),Get2Prong()->PzProng(0)); // from D0
+  TVector3 p3Trk1(Get2Prong()->PxProng(1),Get2Prong()->PyProng(1),Get2Prong()->PzProng(1)); // from D0
+  TVector3 p3Trk2(PxProng(0),PyProng(0),PzProng(0)); // pi_s
+
+  Double_t alpha = p3Trk0.Angle(p3Trk2);
+  Double_t beta = p3Trk1.Angle(p3Trk2);
+
+  Double_t cosphi01 = TMath::Cos(alpha) / TMath::Cos(AngleD0dkpPisoft());
+  Double_t cosphi02 = TMath::Cos(beta) / TMath::Cos(AngleD0dkpPisoft());
+
+  Double_t phi01 = TMath::ACos(cosphi01);
+  Double_t phi02 = TMath::ACos(cosphi02);
+  Double_t phi00 = p3Trk0.Angle(p3Trk1);
+
+  if((phi01>phi00) || (phi02>phi00)) return kFALSE;
+  return kTRUE;
+}

diff --git a/PWG3/vertexingHF/AliAODRecoCascadeHF.h b/PWG3/vertexingHF/AliAODRecoCascadeHF.h
index 8f2ae31bb0204e23c2fc2cfca8e7e4ea8c39277d..c67a42940c6074db9c661ec8d8bc942dbf3eeb43 100644 (file)
--- a/PWG3/vertexingHF/AliAODRecoCascadeHF.h
+++ b/PWG3/vertexingHF/AliAODRecoCascadeHF.h
@@ -69,6 +69,7 @@ class AliAODRecoCascadeHF : public AliAODRecoDecayHF2Prong {
   Double_t InvMassDstarKpipi() const;
   Double_t DeltaInvMass() const {return (InvMassDstarKpipi()-InvMassD0());}
   Double_t AngleD0dkpPisoft() const;
+  Bool_t   TrigonometricalCut() const;
 
   // Lc invariant mass
   Double_t InvMassLctoK0sP() const {