Adding cuts on impact parameters

diff --git a/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoESDTrackCut.cxx b/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoESDTrackCut.cxx
index 3841c33406178f23b8e12e448a0a2909d3439e8f..c1298c133691dcbe1b98ef79ec1cc909860550ea 100644 (file)
--- a/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoESDTrackCut.cxx
+++ b/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoESDTrackCut.cxx
@@ -95,7 +95,9 @@ AliFemtoESDTrackCut::AliFemtoESDTrackCut() :
     fNTracksPassed(0),
     fNTracksFailed(0),
     fRemoveKinks(kFALSE),
-    fMostProbable(0)
+    fMostProbable(0), 
+    fMaxImpactXY(1000.0),
+    fMaxImpactZ(1000.0)
 {
   // Default constructor
   fNTracksPassed = fNTracksFailed = 0;
@@ -156,6 +158,12 @@ bool AliFemtoESDTrackCut::Pass(const AliFemtoTrack* track)
       return false;
     }
        
+  if (fMaxImpactXY < track->ImpactD())
+    return false;
+
+  if (fMaxImpactZ < track->ImpactZ())
+    return false;
+  
   if (fMaxSigmaToVertex < track->SigmaToVertex()) {
     return false;
   }
@@ -352,6 +360,10 @@ TList *AliFemtoESDTrackCut::ListSettings()
   tListSetttings->AddLast(new TObjString(buf));
   snprintf(buf, 200, "AliFemtoESDTrackCut.maxsigmatovertex=%f", fMaxSigmaToVertex);
   tListSetttings->AddLast(new TObjString(buf));
+  snprintf(buf, 200, "AliFemtoESDTrackCut.maximpactxy=%f", fMaxImpactXY);
+  tListSetttings->AddLast(new TObjString(buf));
+  snprintf(buf, 200, "AliFemtoESDTrackCut.maximpactz=%f", fMaxImpactZ);
+  tListSetttings->AddLast(new TObjString(buf));
   if (fMostProbable) {
     if (fMostProbable == 2)
       snprintf(buf, 200, "AliFemtoESDTrackCut.mostprobable=%s", "Pion");

diff --git a/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoESDTrackCut.h b/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoESDTrackCut.h
index 83398a3c5fe763944ce0c302a742dcbd82de2614..d1cd78966f536282d7c6f73fa4dc80dc2e847fb2 100644 (file)
--- a/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoESDTrackCut.h
+++ b/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoESDTrackCut.h
@@ -47,6 +47,8 @@ class AliFemtoESDTrackCut : public AliFemtoTrackCut
   void SetMaxITSChiNdof(const float& maxchi);
   void SetMaxTPCChiNdof(const float& maxchi);
   void SetMaxSigmaToVertex(const float& maxsig);
+  void SetMaxImpactXY(const float& maximpxy);
+  void SetMaxImpactZ(const float& maximpz);
   void SetMostProbablePion();
   void SetMostProbableKaon();
   void SetMostProbableProton();
@@ -76,6 +78,9 @@ class AliFemtoESDTrackCut : public AliFemtoTrackCut
   bool              fRemoveKinks;        // if true particles with any kink label will not pass
   int               fMostProbable;       // this particle type is required to be most probable
 
+  float             fMaxImpactXY;        // Max XY impact parameter
+  float             fMaxImpactZ;         // Max Z impact parameter
+
   float PidFractionElectron(float mom) const;
   float PidFractionPion(float mom) const;
   float PidFractionKaon(float mom) const;
@@ -107,6 +112,8 @@ inline void AliFemtoESDTrackCut::SetNoMostProbable() { fMostProbable = 0; }
 inline void AliFemtoESDTrackCut::SetMaxITSChiNdof(const float& maxchi) { fMaxITSchiNdof = maxchi; }
 inline void AliFemtoESDTrackCut::SetMaxTPCChiNdof(const float& maxchi) { fMaxTPCchiNdof = maxchi; }
 inline void AliFemtoESDTrackCut::SetMaxSigmaToVertex(const float& maxsig) { fMaxSigmaToVertex = maxsig; }
+inline void AliFemtoESDTrackCut::SetMaxImpactXY(const float& maximpxy) { fMaxImpactXY = maximpxy; }
+inline void AliFemtoESDTrackCut::SetMaxImpactZ(const float& maximpz) { fMaxImpactZ = maximpz; }
 
 #endif