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

[blackhole at projects.hepforge.org]

Author: hoeth
Date: Tue Nov 10 00:57:24 2009
New Revision: 2049

Log:
STAR_2006_S6500200: They really cut on rapidity.
Those heavy-ion people are always good for a surprise.

Modified:
   trunk/src/Analyses/STAR_2006_S6500200.cc

Modified: trunk/src/Analyses/STAR_2006_S6500200.cc
==============================================================================
--- trunk/src/Analyses/STAR_2006_S6500200.cc	Tue Nov 10 00:54:56 2009	(r2048)
+++ trunk/src/Analyses/STAR_2006_S6500200.cc	Tue Nov 10 00:57:24 2009	(r2049)
@@ -19,8 +19,8 @@
 
     /// Book projections and histograms
     void init() {
-      IdentifiedFinalState pionfs(-0.5, 0.5, 0.3*GeV);
-      IdentifiedFinalState protonfs(-0.5, 0.5, 0.4*GeV);
+      IdentifiedFinalState pionfs(-2.5, 2.5, 0.3*GeV);
+      IdentifiedFinalState protonfs(-2.5, 2.5, 0.4*GeV);
       pionfs.acceptIdPair(PIPLUS);
       protonfs.acceptIdPair(PROTON);
       addProjection(pionfs, "PIONFS");
@@ -40,7 +40,7 @@
       const IdentifiedFinalState& pionfs = applyProjection<IdentifiedFinalState>(event, "PIONFS");
       const IdentifiedFinalState& protonfs = applyProjection<IdentifiedFinalState>(event, "PROTONFS");
       foreach (const Particle& p, pionfs.particles()) {
-        if (fabs(p.momentum().eta()) < 0.5) {
+        if (fabs(p.momentum().rapidity()) < 0.5) {
           const double pT = p.momentum().pT() / GeV;
           if (p.pdgId()>0) {
             _h_pT_piplus->fill(pT, weight/pT);
@@ -51,7 +51,7 @@
         }
       }
       foreach (const Particle& p, protonfs.particles()) {
-        if (fabs(p.momentum().eta()) < 0.5) {
+        if (fabs(p.momentum().rapidity()) < 0.5) {
           const double pT = p.momentum().pT() / GeV;
           if (p.pdgId()>0) {
             _h_pT_proton->fill(pT, weight/pT);