Added the possibility to cut in dip angle, to reject electrons close to the kinematic...

diff --git a/PWG2/RESONANCES/AliRsnCutStd.cxx b/PWG2/RESONANCES/AliRsnCutStd.cxx
index c90c62bce1dad3e21da33ce694e4fa7f6be43f79..6f26217199ff7153151a7ff27eea9dfc2e5c9ce4 100644 (file)
--- a/PWG2/RESONANCES/AliRsnCutStd.cxx
+++ b/PWG2/RESONANCES/AliRsnCutStd.cxx
@@ -68,6 +68,7 @@ AliRsnCutStd::AliRsnCutStd
     case kPtLeading:
     case kEta:
     case kY:
+    case kDipAngle:
     case kThetaDeg:
       if (fVarType != kDouble) 
       {
@@ -115,6 +116,7 @@ AliRsnCutStd::AliRsnCutStd
     case kPtLeading:
     case kEta:
     case kY:
+    case kDipAngle:
     case kThetaDeg:
       break;
     // other cuts are not based on a value, so no problem
@@ -142,6 +144,7 @@ AliRsnCut::EVarType AliRsnCutStd::CheckType()
     case kPtLeading:
     case kEta:
     case kY:
+    case kDipAngle:
     case kThetaDeg:
       return kDouble;
     // other cuts are not based on a value, so no problem
@@ -239,6 +242,12 @@ Bool_t AliRsnCutStd::IsMotherSelected(AliRsnMother * const mother)
     case kY:
       fCutValueD = ref.Rapidity();
       return OkRange();
+    case kDipAngle:
+      fCutValueD  = mother->GetDaughter(0)->P().Perp() * mother->GetDaughter(1)->P().Perp();
+      fCutValueD += mother->GetDaughter(0)->P().Pz() * mother->GetDaughter(1)->P().Pz();
+      fCutValueD += mother->GetDaughter(0)->P().Mag() * mother->GetDaughter(1)->P().Mag();
+      fCutValueD  = TMath::ACos(fCutValueD);
+      return OkRange();
     case kSameLabel:
       return mother->IsLabelEqual();
     default:

diff --git a/PWG2/RESONANCES/AliRsnCutStd.h b/PWG2/RESONANCES/AliRsnCutStd.h
index a96a067d548cf2744a3b3eb98e4c3c99797433bf..500722c25c8c74834c8826e20db04a3dfac78453 100644 (file)
--- a/PWG2/RESONANCES/AliRsnCutStd.h
+++ b/PWG2/RESONANCES/AliRsnCutStd.h
@@ -33,6 +33,7 @@ class AliRsnCutStd : public AliRsnCut
       kThetaDeg,
       kMult,
       kPtLeading,
+      kDipAngle,
       // cut without reference values
       kCharge,
       kSameLabel,