[Rivet-svn] r3104 - in trunk: data/plotinfo src/Analyses

[blackhole at projects.hepforge.org]

Author: buckley
Date: Mon May 23 17:17:29 2011
New Revision: 3104

Log:
Adding plotting of the PID occurrence for *all* particles in the GenEvent, unphysical or otherwise. This shows up that UnstableFinalState is not finding certain decayed hadrons as written in the HERWIG event record. We need to fix this if there is going to be any tuning of HERWIG hadronisation flavour details using Rivet.

Modified:
   trunk/data/plotinfo/MC_GENERIC.plot
   trunk/src/Analyses/MC_GENERIC.cc

Modified: trunk/data/plotinfo/MC_GENERIC.plot
==============================================================================
--- trunk/data/plotinfo/MC_GENERIC.plot	Mon May 23 17:10:06 2011	(r3103)
+++ trunk/data/plotinfo/MC_GENERIC.plot	Mon May 23 17:17:29 2011	(r3104)
@@ -10,28 +10,34 @@
 YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}E$ [GeV$^{-1}$]
 # END PLOT
 
-# BEGIN PLOT /MC_GENERIC/EtaCh
-Title=Pseudorapidity of charged particles
+# BEGIN PLOT /MC_GENERIC/Eta.*
 XLabel=$\eta$
 YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\eta$
 LogY=0
 LegendYPos=0.5
 # END PLOT
 
-# BEGIN PLOT /MC_GENERIC/EtaChPMRatio
-Title=Pseudorapidity $+/-$ ratio of charged particles
+# BEGIN PLOT /MC_GENERIC/Eta$
+Title=Pseudorapidity of all particles
+# END PLOT
+
+# BEGIN PLOT /MC_GENERIC/EtaCh
+Title=Pseudorapidity of charged particles
+# END PLOT
+
+# BEGIN PLOT /MC_GENERIC/Eta.*PMRatio
 XLabel=$|\eta|$
 YLabel=$\mathrm{d}\sigma/\mathrm{d}\eta_+ \; / \; \mathrm{d}\sigma/\mathrm{d}\eta_-$
 LogY=0
 LegendYPos=0.5
 # END PLOT
 
-# BEGIN PLOT /MC_GENERIC/Eta
-Title=Pseudorapidity of all particles
-XLabel=$\eta$
-YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\eta$
-LogY=0
-LegendYPos=0.5
+# BEGIN PLOT /MC_GENERIC/EtaPMRatio
+Title=Pseudorapidity $+/-$ ratio of all particles
+# END PLOT
+
+# BEGIN PLOT /MC_GENERIC/EtaChPMRatio
+Title=Pseudorapidity $+/-$ ratio of charged particles
 # END PLOT
 
 # BEGIN PLOT /MC_GENERIC/EtaK
@@ -69,14 +75,6 @@
 LegendYPos=0.5
 # END PLOT
 
-# BEGIN PLOT /MC_GENERIC/EtaPMRatio
-Title=Pseudorapidity $+/-$ ratio of all particles
-XLabel=$|\eta|$
-YLabel=$\mathrm{d}\sigma/\mathrm{d}\eta_+ \; / \; \mathrm{d}\sigma/\mathrm{d}\eta_-$
-LogY=0
-LegendYPos=0.5
-# END PLOT
-
 # BEGIN PLOT /MC_GENERIC/MultCh
 Title=Charged multiplicity of events
 XLabel=$N_text{ch}$
@@ -91,18 +89,22 @@
 LogY=0
 # END PLOT
 
-# BEGIN PLOT /MC_GENERIC/MultsDecayedPIDs
-Title=Average per-event multiplicity of decayed particle IDs
+# BEGIN PLOT /MC_GENERIC/Mults.*PIDs
 XLabel=$|\text{PID}|$
 YLabel=$\langle N_\text{PID} \rangle$
-LogY=0
+LogY=1
+# END PLOT
+
+# BEGIN PLOT /MC_GENERIC/MultsDecayedPIDs
+Title=Average per-event multiplicity of decayed particle IDs
 # END PLOT
 
 # BEGIN PLOT /MC_GENERIC/MultsStablePIDs
 Title=Average per-event multiplicity of stable particle IDs
-XLabel=$|\text{PID}|$
-YLabel=$\langle N_\text{PID} \rangle$
-LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_GENERIC/MultsAllPIDs
+Title=Average per-event multiplicity of \emph{all} (unphysical) particle IDs
 # END PLOT
 
 # BEGIN PLOT /MC_GENERIC/PhiCh

Modified: trunk/src/Analyses/MC_GENERIC.cc
==============================================================================
--- trunk/src/Analyses/MC_GENERIC.cc	Mon May 23 17:10:06 2011	(r3103)
+++ trunk/src/Analyses/MC_GENERIC.cc	Mon May 23 17:17:29 2011	(r3104)
@@ -43,6 +43,7 @@
 
       _histStablePIDs  = bookHistogram1D("MultsStablePIDs", 3335, -0.5, 3334.5);
       _histDecayedPIDs = bookHistogram1D("MultsDecayedPIDs", 3335, -0.5, 3334.5);
+      _histAllPIDs  = bookHistogram1D("MultsAllPIDs", 3335, -0.5, 3334.5);
 
       _histPt    = bookHistogram1D("Pt", 300, 0, 30);
       _histPtCh  = bookHistogram1D("PtCh", 300, 0, 30);
@@ -79,7 +80,12 @@
     void analyze(const Event& event) {
       const double weight = event.weight();
 
-      // Analyse and print some info
+      // Unphysical (debug) plotting of all PIDs in the event, physical or otherwise
+      foreach (const GenParticle* gp, particles(event.genEvent())) {
+        _histAllPIDs->fill(abs(gp->pdg_id()), weight);
+      }
+
+      // Charged + neutral final state
       const FinalState& cnfs = applyProjection<FinalState>(event, "FS");
       MSG_DEBUG("Total multiplicity = " << cnfs.size());
       _histMult->fill(cnfs.size(), weight);
@@ -153,6 +159,7 @@
 
       scale(_histStablePIDs, 1/sumOfWeights());
       scale(_histDecayedPIDs, 1/sumOfWeights());
+      scale(_histAllPIDs, 1/sumOfWeights());
 
       scale(_histEta, 1/sumOfWeights());
       scale(_histEtaCh, 1/sumOfWeights());
@@ -193,7 +200,7 @@
     //@{
     /// Histograms
     AIDA::IHistogram1D *_histMult, *_histMultCh;
-    AIDA::IHistogram1D *_histStablePIDs, *_histDecayedPIDs;
+    AIDA::IHistogram1D *_histStablePIDs, *_histDecayedPIDs, *_histAllPIDs;
     AIDA::IHistogram1D *_histEtaPi, *_histEtaK, *_histEtaLambda;
     AIDA::IProfile1D   *_histEtaSumEt;
     AIDA::IHistogram1D *_histEta, *_histEtaCh;