[Rivet-svn] r2050 - in trunk: data/anainfo src/Analyses

[blackhole at projects.hepforge.org]

Author: hoeth
Date: Tue Nov 10 01:13:37 2009
New Revision: 2050

Log:
STAR_2006_S6500200: added non-single diffractive trigger

Modified:
   trunk/data/anainfo/STAR_2006_S6500200.info
   trunk/src/Analyses/STAR_2006_S6500200.cc

Modified: trunk/data/anainfo/STAR_2006_S6500200.info
==============================================================================
--- trunk/data/anainfo/STAR_2006_S6500200.info	Tue Nov 10 00:57:24 2009	(r2049)
+++ trunk/data/anainfo/STAR_2006_S6500200.info	Tue Nov 10 01:13:37 2009	(r2050)
@@ -4,8 +4,9 @@
 Experiment: STAR
 Collider: RHIC pp 200 GeV
 SpiresID: 6500200
-Status: UNDER DEVELOPMENT, UNVALIDATED, DO NOT USE.
+Status: UNVALIDATED
 Authors:
+ - Bedanga Mohanty <bedanga at rcf.bnl.gov>
  - Hendrik Hoeth <hendrik.hoeth at cern.ch>
 References:
  - Phys. Lett. B637, 161

Modified: trunk/src/Analyses/STAR_2006_S6500200.cc
==============================================================================
--- trunk/src/Analyses/STAR_2006_S6500200.cc	Tue Nov 10 00:57:24 2009	(r2049)
+++ trunk/src/Analyses/STAR_2006_S6500200.cc	Tue Nov 10 01:13:37 2009	(r2050)
@@ -1,6 +1,7 @@
 // -*- C++ -*-
 #include "Rivet/Analysis.hh"
 #include "Rivet/Tools/Logging.hh"
+#include "Rivet/Projections/ChargedFinalState.hh"
 #include "Rivet/Projections/IdentifiedFinalState.hh"
 #include "Rivet/RivetAIDA.hh"
 
@@ -19,6 +20,10 @@
 
     /// Book projections and histograms
     void init() {
+      ChargedFinalState bbc1(-5.0,-3.3, 0.0*GeV); // beam-beam-counter trigger
+      ChargedFinalState bbc2( 3.3, 5.0, 0.0*GeV); // beam-beam-counter trigger
+      addProjection(bbc1, "BBC1");
+      addProjection(bbc2, "BBC2");
       IdentifiedFinalState pionfs(-2.5, 2.5, 0.3*GeV);
       IdentifiedFinalState protonfs(-2.5, 2.5, 0.4*GeV);
       pionfs.acceptIdPair(PIPLUS);
@@ -37,6 +42,13 @@
     void analyze(const Event& event) {
       const double weight = event.weight();
 
+      const ChargedFinalState& bbc1 = applyProjection<ChargedFinalState>(event, "BBC1");
+      const ChargedFinalState& bbc2 = applyProjection<ChargedFinalState>(event, "BBC2");
+      if (bbc1.size()<1 || bbc2.size()<1) {
+        getLog() << Log::DEBUG << "Failed beam-beam-counter trigger" << std::endl;
+        vetoEvent;
+      }
+
       const IdentifiedFinalState& pionfs = applyProjection<IdentifiedFinalState>(event, "PIONFS");
       const IdentifiedFinalState& protonfs = applyProjection<IdentifiedFinalState>(event, "PROTONFS");
       foreach (const Particle& p, pionfs.particles()) {