[Rivet-svn] r1936 - in trunk: . src/Analyses

[blackhole at projects.hepforge.org]

Author: hoeth
Date: Mon Oct 19 15:46:30 2009
New Revision: 1936

Log:
Fix most open issues in the STAR_2009_UE_HELEN analysis.

Modified:
   trunk/ChangeLog
   trunk/src/Analyses/STAR_2009_UE_HELEN.cc

Modified: trunk/ChangeLog
==============================================================================
--- trunk/ChangeLog	Mon Oct 19 15:15:31 2009	(r1935)
+++ trunk/ChangeLog	Mon Oct 19 15:46:30 2009	(r1936)
@@ -16,6 +16,9 @@
 	by looking at the genParticle(), so users need to take care of
 	any manually added particles themselves.
 
+	* Fixing most open issues with the STAR_2009_UE_HELEN analysis.
+	There is only one question left, regarding the away region.
+
 2009-10-17  Andy Buckley  <andy at insectnation.org>
 
 	* Adding parsing of units in cross-sections passed to the "-x"

Modified: trunk/src/Analyses/STAR_2009_UE_HELEN.cc
==============================================================================
--- trunk/src/Analyses/STAR_2009_UE_HELEN.cc	Mon Oct 19 15:15:31 2009	(r1935)
+++ trunk/src/Analyses/STAR_2009_UE_HELEN.cc	Mon Oct 19 15:46:30 2009	(r1936)
@@ -2,8 +2,10 @@
 #include "Rivet/Analysis.hh"
 #include "Rivet/RivetAIDA.hh"
 #include "Rivet/Tools/Logging.hh"
-#include "Rivet/Projections/FinalState.hh"
 #include "Rivet/Projections/ChargedFinalState.hh"
+#include "Rivet/Projections/NeutralFinalState.hh"
+#include "Rivet/Projections/MergedFinalState.hh"
+#include "Rivet/Projections/VetoedFinalState.hh"
 #include "Rivet/Projections/FastJets.hh"
 #include "fastjet/SISConePlugin.hh"
 
@@ -28,17 +30,30 @@
     //@{
 
     void init() {
-      // Final state for the jet finding
-      const FinalState fsj(-1, 1, 0.2*GeV);
-      addProjection(fsj, "FSJ");
-      // Split-merge is 0.75, so we need to initialize the plugin manually:
-      // R = 0.7, overlap_threshold = 0.75
-      addProjection(FastJets(fsj, fastjet::SISConePlugin(0.7, 0.75)), "AllJets");
-
-      // Charged final state for the distributions
+      // Charged final state, |eta|<1, pT>0.2GeV
       const ChargedFinalState cfs(-1.0, 1.0, 0.2*GeV);
       addProjection(cfs, "CFS");
 
+      // Neutral final state, |eta|<1, ET>0.2GeV (needed for the jets)
+      const NeutralFinalState nfs(-1.0, 1.0, 0.2*GeV);
+      addProjection(nfs, "NFS");
+
+      // STAR can't see neutrons and K^0_L
+      VetoedFinalState vfs(nfs);
+      vfs.vetoNeutrinos();
+      vfs.addVetoPairId(K0L);
+      vfs.addVetoPairId(NEUTRON);
+      addProjection(vfs, "VFS");
+
+      // Jets are reconstructed from charged and neutral particles,
+      // and the cuts are different (pT vs. ET), so we need to merge them.
+      const MergedFinalState jfs(cfs, vfs);
+      addProjection(jfs, "JFS");
+
+      // Split-merge is 0.75, so we need to initialize the plugin manually:
+      // R = 0.7, overlap_threshold = 0.75
+      addProjection(FastJets(jfs, fastjet::SISConePlugin(0.7, 0.75)), "AllJets");
+
       // Book histograms
       _hist_pmaxnchg   = bookProfile1D( 1, 1, 1);
       _hist_pminnchg   = bookProfile1D( 2, 1, 1);
@@ -48,8 +63,8 @@
 
     // Do the analysis
     void analyze(const Event& e) {
-      const FinalState& fsj = applyProjection<FinalState>(e, "FSJ");
-      if (fsj.particles().size() < 1) {
+      const FinalState& cfs = applyProjection<ChargedFinalState>(e, "CFS");
+      if (cfs.particles().size() < 1) {
         getLog() << Log::DEBUG << "Failed multiplicity cut" << endl;
         vetoEvent;
       }
@@ -75,9 +90,6 @@
       // Get the event weight
       const double weight = e.weight();
 
-      // Get the final states to work with for filling the distributions
-      const FinalState& cfs = applyProjection<ChargedFinalState>(e, "CFS");
-
       size_t numOverall(0),     numToward(0),     numTrans1(0),     numTrans2(0),     numAway(0)  ;
       double ptSumOverall(0.0), ptSumToward(0.0), ptSumTrans1(0.0), ptSumTrans2(0.0), ptSumAway(0.0);
       //double EtSumOverall(0.0), EtSumToward(0.0), EtSumTrans1(0.0), EtSumTrans2(0.0), EtSumAway(0.0);