[Rivet-svn] r3362 - trunk/src/Analyses

[blackhole at projects.hepforge.org]

Author: hoeth
Date: Sat Sep 17 11:45:10 2011
New Revision: 3362

Log:
Fix for decay chain downfeed in ALICE_2011_S8945144 (patch provided by ALICE)

Modified:
   trunk/src/Analyses/ALICE_2011_S8945144.cc

Modified: trunk/src/Analyses/ALICE_2011_S8945144.cc
==============================================================================
--- trunk/src/Analyses/ALICE_2011_S8945144.cc	Fri Sep 16 15:41:43 2011	(r3361)
+++ trunk/src/Analyses/ALICE_2011_S8945144.cc	Sat Sep 17 11:45:10 2011	(r3362)
@@ -45,12 +45,22 @@
               _histAveragePt->fill(p.mass()/GeV, p.momentum().pT()/GeV, weight);
               break;
             case 2212:
-              _histPtProtons->fill(p.momentum().pT()/GeV, weight);
-              _histAveragePt->fill(p.mass()/GeV, p.momentum().pT()/GeV, weight);
+              if ( !(p.hasAncestor(3322) ||                             // Xi0
+                     p.hasAncestor(3122) || p.hasAncestor(-3122) ||     // Lambda
+                     p.hasAncestor(3222) || p.hasAncestor(-3222) ||     // Sigma+/-
+                     p.hasAncestor(3312) || p.hasAncestor(-3312) ) ) {  // Xi-/+
+                _histPtProtons->fill(p.momentum().pT()/GeV, weight);
+                _histAveragePt->fill(p.mass()/GeV, p.momentum().pT()/GeV, weight);
+              }
               break;
             case -2212:
-              _histPtAntiProtons->fill(p.momentum().pT()/GeV, weight);
-              _histAveragePt->fill(p.mass()/GeV, p.momentum().pT()/GeV, weight);
+              if ( !(p.hasAncestor(3322) ||                             // Xi0
+                     p.hasAncestor(3122) || p.hasAncestor(-3122) ||     // Lambda
+                     p.hasAncestor(3222) || p.hasAncestor(-3222) ||     // Sigma+/-
+                     p.hasAncestor(3312) || p.hasAncestor(-3312) ) ) {  // Xi-/+
+                _histPtAntiProtons->fill(p.momentum().pT()/GeV, weight);
+                _histAveragePt->fill(p.mass()/GeV, p.momentum().pT()/GeV, weight);
+              }
               break;
             case 321:
               _histPtKaons->fill(p.momentum().pT()/GeV, weight);