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

[blackhole at projects.hepforge.org]

Author: hoeth
Date: Wed Jul 13 14:09:36 2011
New Revision: 3184

Log:
remove dangling comment

Modified:
   trunk/src/Analyses/DELPHI_2002_069_CONF_603.cc

Modified: trunk/src/Analyses/DELPHI_2002_069_CONF_603.cc
==============================================================================
--- trunk/src/Analyses/DELPHI_2002_069_CONF_603.cc	Wed Jul 13 14:01:39 2011	(r3183)
+++ trunk/src/Analyses/DELPHI_2002_069_CONF_603.cc	Wed Jul 13 14:09:36 2011	(r3184)
@@ -5,6 +5,7 @@
 #include "Rivet/Projections/Beam.hh"
 #include "Rivet/Projections/FinalState.hh"
 #include "Rivet/Projections/ChargedFinalState.hh"
+#include "Rivet/Projections/PartonFinalState.hh"
 
 
 /// @todo Use inline PID functions instead
@@ -34,11 +35,14 @@
     void init() {
       addProjection(Beam(), "Beams");
       addProjection(ChargedFinalState(), "FS");
+      addProjection(PartonFinalState(), "PFS");
 
       _histXbprim     = bookHistogram1D(1, 1, 1);
       _histXbweak     = bookHistogram1D(2, 1, 1);
       _histMeanXbprim = bookProfile1D(4, 1, 1);
       _histMeanXbweak = bookProfile1D(5, 1, 1);
+      _histXbpart     = bookHistogram1D("Xbpart", binEdges(1, 1, 1));
+      _histMeanXbpart = bookProfile1D("MeanXbpart", binEdges(5, 1, 1));
     }
 
 
@@ -62,7 +66,15 @@
                                    beams.second.momentum().vector3().mod() ) / 2.0;
       getLog() << Log::DEBUG << "Avg beam momentum = " << meanBeamMom << endl;
    
-   
+      const ParticleVector partons = applyProjection<PartonFinalState>(e, "PFS").partons();
+      foreach (Particle p, partons) {
+        if (abs(p.pdg_id()) == 5) {
+          const double xp = p.momentum().E()/meanBeamMom;
+          _histXbpart->fill(xp, weight);
+          _histMeanXbpart->fill(_histMeanXbpart->binMean(0), xp, weight);
+        }
+      }
+
       foreach (const GenParticle* p, particles(e.genEvent())) {
         const GenVertex* pv = p->production_vertex();
         const GenVertex* dv = p->end_vertex();
@@ -105,6 +117,7 @@
     void finalize() {
       normalize(_histXbprim);
       normalize(_histXbweak);
+      normalize(_histXbpart);
     }
 
 
@@ -116,9 +129,11 @@
 
     AIDA::IHistogram1D *_histXbprim;
     AIDA::IHistogram1D *_histXbweak;
+    AIDA::IHistogram1D *_histXbpart;
 
     AIDA::IProfile1D *_histMeanXbprim;
     AIDA::IProfile1D *_histMeanXbweak;
+    AIDA::IProfile1D *_histMeanXbpart;
 
     //@}