[Rivet-svn] rivet: Add aggregate signal regions to CMS_2016_PAS_SUS_16_14.

[Rivet Mercurial]

details:   https://rivet.hepforge.org/hg/rivet/rev/3eba456fb5b8
branches:  release-2-5-x
changeset: 5697:3eba456fb5b8
user:      Andy Buckley <andy at insectnation.org>
date:      Wed Feb 22 22:24:54 2017 +0000
description:
Add aggregate signal regions to CMS_2016_PAS_SUS_16_14.

diffs (truncated from 84 to 50 lines):

--- a/ChangeLog	Mon Feb 20 23:53:59 2017 +0000
+++ b/ChangeLog	Wed Feb 22 22:24:54 2017 +0000
@@ -1,3 +1,7 @@
+2017-02-22  Andy Buckley  <andy.buckley at cern.ch>
+
+	* Add aggregate signal regions to CMS_2016_PAS_SUS_16_14.
+
 2017-02-18  Andy Buckley  <andy.buckley at cern.ch>
 
 	* Add getEnvParam function, for neater use of environment variable
--- a/src/Analyses/CMS_2016_PAS_SUS_16_14.cc	Mon Feb 20 23:53:59 2017 +0000
+++ b/src/Analyses/CMS_2016_PAS_SUS_16_14.cc	Wed Feb 22 22:24:54 2017 +0000
@@ -58,6 +58,10 @@
           }
         }
       }
+      _h_srcountsagg.resize(12);
+      for (size_t ia = 0; ia < 12; ++ia) {
+        _h_srcountsagg[ia] = bookCounter("agg-" + toString(ia));
+      }
 
     }
 
@@ -134,20 +138,35 @@
       if (deltaPhi(htmissvec, jets24[2]) < 0.3) vetoEvent;
       if (jets24.size() >= 4 && deltaPhi(htmissvec, jets24[3]) < 0.3) vetoEvent;
 
+      // Count jet and b-jet multiplicities
+      const size_t nj = jets24.size();
+      size_t nbj = 0;
+      for (const Jet& j : jets24)
+        if (j.bTagged()) nbj += 1;
+
 
       ////////
 
 
-      // Calculate a bin index for this event
-      // Nj bin
+      // Fill the aggregate signal regions first
+      if (nj >= 3 && nbj == 0 && ht >  500*GeV && htmiss > 500*GeV) _h_srcountsagg[ 0]->fill(event.weight());
+      if (nj >= 3 && nbj == 0 && ht > 1500*GeV && htmiss > 750*GeV) _h_srcountsagg[ 1]->fill(event.weight());
+      if (nj >= 5 && nbj == 0 && ht >  500*GeV && htmiss > 500*GeV) _h_srcountsagg[ 2]->fill(event.weight());
+      if (nj >= 5 && nbj == 0 && ht > 1500*GeV && htmiss > 750*GeV) _h_srcountsagg[ 3]->fill(event.weight());
+      if (nj >= 9 && nbj == 0 && ht > 1500*GeV && htmiss > 750*GeV) _h_srcountsagg[ 4]->fill(event.weight());
+      if (nj >= 3 && nbj >= 2 && ht >  500*GeV && htmiss > 500*GeV) _h_srcountsagg[ 5]->fill(event.weight());
+      if (nj >= 3 && nbj >= 1 && ht >  750*GeV && htmiss > 750*GeV) _h_srcountsagg[ 6]->fill(event.weight());
+      if (nj >= 5 && nbj >= 3 && ht >  500*GeV && htmiss > 500*GeV) _h_srcountsagg[ 7]->fill(event.weight());
+      if (nj >= 5 && nbj >= 2 && ht > 1500*GeV && htmiss > 750*GeV) _h_srcountsagg[ 8]->fill(event.weight());
+      if (nj >= 9 && nbj >= 3 && ht >  750*GeV && htmiss > 750*GeV) _h_srcountsagg[ 9]->fill(event.weight());