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

[blackhole at projects.hepforge.org]

Author: hoeth
Date: Fri Nov  6 12:46:33 2009
New Revision: 2032

Log:
Unless we know if STAR used eta or rapidity, I'll use eta
in the STAR_2006_S6500200 analysis, since this makes more sense.

Modified:
   trunk/src/Analyses/STAR_2006_S6500200.cc

Modified: trunk/src/Analyses/STAR_2006_S6500200.cc
==============================================================================
--- trunk/src/Analyses/STAR_2006_S6500200.cc	Fri Nov  6 10:45:07 2009	(r2031)
+++ trunk/src/Analyses/STAR_2006_S6500200.cc	Fri Nov  6 12:46:33 2009	(r2032)
@@ -19,8 +19,8 @@
 
     /// Book projections and histograms
     void init() {
-      IdentifiedFinalState pionfs(-3.0, 3.0, 0.3*GeV);
-      IdentifiedFinalState protonfs(-3.0, 3.0, 0.4*GeV);
+      IdentifiedFinalState pionfs(-0.5, 0.5, 0.3*GeV);
+      IdentifiedFinalState protonfs(-0.5, 0.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().rapidity()) < 0.5) {
+        if (fabs(p.momentum().eta()) < 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().rapidity()) < 0.5) {
+        if (fabs(p.momentum().eta()) < 0.5) {
           const double pT = p.momentum().pT() / GeV;
           if (p.pdgId()>0) {
             _h_pT_proton->fill(pT, weight/pT);