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

blackhole at projects.hepforge.org blackhole at projects.hepforge.org
Wed Sep 14 09:55:10 BST 2011


Author: buckley
Date: Wed Sep 14 09:55:10 2011
New Revision: 3351

Log:
Adding a VH->bb MC analysis

Added:
   trunk/data/anainfo/MC_VH2BB.info
   trunk/data/plotinfo/MC_VH2BB.plot
   trunk/src/Analyses/MC_VH2BB.cc
Modified:
   trunk/data/anainfo/Makefile.am
   trunk/data/plotinfo/Makefile.am
   trunk/src/Analyses/Makefile.am

Added: trunk/data/anainfo/MC_VH2BB.info
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/anainfo/MC_VH2BB.info	Wed Sep 14 09:55:10 2011	(r3351)
@@ -0,0 +1,18 @@
+Name: MC_VH2BB
+Year: 2011
+Summary: MC unboosted VH2bb validation plots
+Experiment: ATLAS
+Collider: LHC
+Status: UNVALIDATED
+Authors:
+ - Ben Smart <bsmart at cern.ch>
+ - Andy Buckley <andy.buckley at cern.ch>
+RunInfo:
+  $VH$ with $H \to b \bar{b}$ and the vector boson decaying to electron or
+  muon channels.
+NumEvents: 100000
+Beams: [p+, p+]
+Energies: [7000]
+PtCuts: [0]
+Description:
+  'Various plots for characterising the process $V H \to b\bar{b}$'

Modified: trunk/data/anainfo/Makefile.am
==============================================================================
--- trunk/data/anainfo/Makefile.am	Tue Sep 13 18:00:19 2011	(r3350)
+++ trunk/data/anainfo/Makefile.am	Wed Sep 14 09:55:10 2011	(r3351)
@@ -110,6 +110,7 @@
   MC_TTBAR.info \
   MC_JETS.info \
   MC_PHOTONJETS.info \
+  MC_VH2BB.info \
   MC_HJETS.info \
   MC_WJETS.info \
   MC_WPOL.info \

Added: trunk/data/plotinfo/MC_VH2BB.plot
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/data/plotinfo/MC_VH2BB.plot	Wed Sep 14 09:55:10 2011	(r3351)
@@ -0,0 +1,1165 @@
+# BEGIN PLOT /MC_VH2BB/pure_bb_Delta_eta
+Title= $\Delta\eta$ between b-hadron pair (without cuts)
+XLabel= $\Delta\eta$
+YLabel= $\text{d}\sigma/\text{d}\Delta\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_bb_Delta_phi
+Title= $\Delta\phi$ between b-hadron pair (without cuts)
+XLabel= $\Delta\phi$
+YLabel= $\text{d}\sigma/\text{d}\Delta\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_bb_Delta_pT
+Title= $\Delta{p_T}$ between b-hadron pair (without cuts)
+XLabel= $\Delta{p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}\Delta{p_T}$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_bb_Delta_R
+Title= $\Delta{R}$ between b-hadron pair (without cuts)
+XLabel= $\Delta{R}$
+YLabel= $\text{d}\sigma/\text{d}\Delta{R}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_b_eta
+Title= $\eta$ of b-hadron (without cuts)
+XLabel= $\eta$
+YLabel= $\text{d}\sigma/\text{d}\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_b_multiplicity
+Title= b-hadron multiplicity (without cuts)
+XLabel= $N_{b-Hadron}$
+YLabel= $\text{d}\sigma/\text{d}N_{b-Hadron}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_b_phi
+Title= $\phi$ of b-hadron (without cuts)
+XLabel= $\phi$
+YLabel= $\text{d}\sigma/\text{d}\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_b_pT
+Title= ${p_T}$ of b-hadron (without cuts)
+XLabel= ${p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}{p_T}$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_H_eta_using_bb
+Title= $\eta$ of Higgs candidate from b-hadron pair (without cuts)
+XLabel= $\eta$
+YLabel= $\text{d}\sigma/\text{d}\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_H_mass_using_bb
+Title= mass of Higgs candidate from b-hadron pair (without cuts)
+XLabel= $m$
+YLabel= $\text{d}\sigma/\text{d}m$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_H_phi_using_bb
+Title= $\phi$ of Higgs candidate from b-hadron pair (without cuts)
+XLabel= $\phi$
+YLabel= $\text{d}\sigma/\text{d}\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_H_pT_using_bb
+Title= ${p_T}$ of Higgs candidate from b-hadron pair (without cuts)
+XLabel= ${p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}p_T$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_VBbb_Delta_eta
+Title= $\Delta\eta$ between Higgs candidate and W (without cuts)
+XLabel= $\Delta\eta$
+YLabel= $\text{d}\sigma/\text{d}\Delta\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_VBbb_Delta_phi
+Title= $\Delta\phi$ between Higgs candidate and W (without cuts)
+XLabel= $\Delta\phi$
+YLabel= $\text{d}\sigma/\text{d}\Delta\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_VBbb_Delta_pT
+Title= $\Delta{p_T}$ between Higgs candidate and W (without cuts)
+XLabel= $\Delta{p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}\Delta{p_T}$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_VBbb_Delta_R
+Title= $\Delta{R}$ between Higgs candidate and W (without cuts)
+XLabel= $\Delta{R}$
+YLabel= $\text{d}\sigma/\text{d}\Delta{R}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_cuts_bb_Delta_eta
+Title= $\Delta\eta$ between b-hadron pair (with cuts)
+XLabel= $\Delta\eta$
+YLabel= $\text{d}\sigma/\text{d}\Delta\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_cuts_bb_Delta_phi
+Title= $\Delta\phi$ between b-hadron pair (with cuts)
+XLabel= $\Delta\phi$
+YLabel= $\text{d}\sigma/\text{d}\Delta\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_cuts_bb_Delta_pT
+Title= $\Delta{p_T}$ between b-hadron pair (with cuts)
+XLabel= $\Delta{p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}\Delta{p_T}$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_cuts_bb_Delta_R
+Title= $\Delta{R}$ between b-hadron pair (with cuts)
+XLabel= $\Delta{R}$
+YLabel= $\text{d}\sigma/\text{d}\Delta{R}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_cuts_b_eta
+Title= $\eta$ of b-hadron (with cuts)
+XLabel= $\eta$
+YLabel= $\text{d}\sigma/\text{d}\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_cuts_b_phi
+Title= $\phi$ of b-hadron (with cuts)
+XLabel= $\phi$
+YLabel= $\text{d}\sigma/\text{d}\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_cuts_b_pT
+Title= ${p_T}$ of b-hadron (with cuts)
+XLabel= ${p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}{p_T}$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_cuts_H_eta_using_bb
+Title= $\eta$ of Higgs candidate from b-hadron pair (with cuts)
+XLabel= $\eta$
+YLabel= $\text{d}\sigma/\text{d}\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_cuts_H_mass_using_bb
+Title= mass of Higgs candidate from b-hadron pair (with cuts)
+XLabel= $m$
+YLabel= $\text{d}\sigma/\text{d}m$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_cuts_H_phi_using_bb
+Title= $\phi$ of Higgs candidate from b-hadron pair (with cuts)
+XLabel= $\phi$
+YLabel= $\text{d}\sigma/\text{d}\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_cuts_H_pT_using_bb
+Title= ${p_T}$ of Higgs candidate from b-hadron pair (with cuts)
+XLabel= ${p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}p_T$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_cuts_VBbb_Delta_eta
+Title= $\Delta\eta$ between Higgs candidate and W (with cuts)
+XLabel= $\Delta\eta$
+YLabel= $\text{d}\sigma/\text{d}\Delta\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_cuts_VBbb_Delta_phi
+Title= $\Delta\phi$ between Higgs candidate and W (with cuts)
+XLabel= $\Delta\phi$
+YLabel= $\text{d}\sigma/\text{d}\Delta\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_cuts_VBbb_Delta_pT
+Title= $\Delta{p_T}$ between Higgs candidate and W (with cuts)
+XLabel= $\Delta{p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}\Delta{p_T}$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_cuts_VBbb_Delta_R
+Title= $\Delta{R}$ between Higgs candidate and W (with cuts)
+XLabel= $\Delta{R}$
+YLabel= $\text{d}\sigma/\text{d}\Delta{R}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_cuts_b_multiplicity
+Title= b-hadron multiplicity (with cuts)
+XLabel= $N_{b-Hadron}$
+YLabel= $\text{d}\sigma/\text{d}N_{b-Hadron}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_bb_Delta_eta
+Title= $\Delta\eta$ between b-jet pair (anti-$k_t$ $R=0.4$ jets without cuts)
+XLabel= $\Delta\eta$
+YLabel= $\text{d}\sigma/\text{d}\Delta\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_bb_Delta_phi
+Title= $\Delta\phi$ between b-jet pair (anti-$k_t$ $R=0.4$ jets without cuts)
+XLabel= $\Delta\phi$
+YLabel= $\text{d}\sigma/\text{d}\Delta\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_bb_Delta_pT
+Title= $\Delta{p_T}$ between b-jet pair (anti-$k_t$ $R=0.4$ jets without cuts)
+XLabel= $\Delta{p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}\Delta{p_T}$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_bb_Delta_R
+Title= $\Delta{R}$ between b-jet pair (anti-$k_t$ $R=0.4$ jets without cuts)
+XLabel= $\Delta{R}$
+YLabel= $\text{d}\sigma/\text{d}\Delta{R}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_b_jet_eta
+Title= $\eta$ of b-jet (anti-$k_t$ $R=0.4$ jet without cuts)
+XLabel= $\eta$
+YLabel= $\text{d}\sigma/\text{d}\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_b_jet_multiplicity
+Title= b-jet multiplicity (anti-$k_t$ $R=0.4$ jet without cuts)
+XLabel= $N_{b-jet}$
+YLabel= $\text{d}\sigma/\text{d}N_{b-jet}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_b_jet_phi
+Title= $\phi$ of b-jet (anti-$k_t$ $R=0.4$ jet without cuts)
+XLabel= $\phi$
+YLabel= $\text{d}\sigma/\text{d}\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_b_jet_pT
+Title= ${p_T}$ of b-jet (anti-$k_t$ $R=0.4$ jet without cuts)
+XLabel= ${p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}p_T$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_H_eta_using_bb
+Title= $\eta$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.4$ jets without cuts)
+XLabel= $\eta$
+YLabel= $\text{d}\sigma/\text{d}\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_H_mass_using_bb
+Title= $m$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.4$ jets without cuts)
+XLabel= $m [GeV]$
+YLabel= $\text{d}\sigma/\text{d}m$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_H_phi_using_bb
+Title= $\phi$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.4$ jets without cuts)
+XLabel= $\phi$
+YLabel= $\text{d}\sigma/\text{d}\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_H_pT_using_bb
+Title= ${p_T}$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.4$ jets without cuts)
+XLabel= $p_T [GeV]$
+YLabel= $\text{d}\sigma/\text{d}p_T$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_jet_eta
+Title= $\eta$ of jet (anti-$k_t$ $R=0.4$ jet without cuts)
+XLabel= $\eta$
+YLabel= $\text{d}\sigma/\text{d}\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_jet_multiplicity
+Title= jet multiplicity (anti-$k_t$ $R=0.4$ jet without cuts)
+XLabel= $N_{jet}$
+YLabel= $\text{d}\sigma/\text{d}N_{jet}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_jet_phi
+Title= $\phi$ of jet (anti-$k_t$ $R=0.4$ jet without cuts)
+XLabel= $\phi$
+YLabel= $\text{d}\sigma/\text{d}\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_jet_pT
+Title= ${p_T}$ of jet (anti-$k_t$ $R=0.4$ jet without cuts)
+XLabel= ${p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}p_T$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_VBbb_Delta_eta
+Title= $\Delta\eta$ between Higgs candidate and W (anti-$k_t$ $R=0.4$ jets without cuts)
+XLabel= $\Delta\eta$
+YLabel= $\text{d}\sigma/\text{d}\Delta\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_VBbb_Delta_phi
+Title= $\Delta\phi$ between Higgs candidate and W (anti-$k_t$ $R=0.4$ jets without cuts)
+XLabel= $\Delta\phi$
+YLabel= $\text{d}\sigma/\text{d}\Delta\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_VBbb_Delta_pT
+Title= $\Delta{p_T}$ between Higgs candidate and W (anti-$k_t$ $R=0.4$ jets without cuts)
+XLabel= $\Delta{p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}\Delta{p_T}$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_VBbb_Delta_R
+Title= $\Delta{R}$ between Higgs candidate and W (anti-$k_t$ $R=0.4$ jets without cuts)
+XLabel= $\Delta{R}$
+YLabel= $\text{d}\sigma/\text{d}\Delta{R}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_cuts_bb_Delta_eta
+Title= $\Delta\eta$ between b-jet pair (anti-$k_t$ $R=0.4$ jets with cuts)
+XLabel= $\Delta\eta$
+YLabel= $\text{d}\sigma/\text{d}\Delta\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_cuts_bb_Delta_phi
+Title= $\Delta\phi$ between b-jet pair (anti-$k_t$ $R=0.4$ jets with cuts)
+XLabel= $\Delta\phi$
+YLabel= $\text{d}\sigma/\text{d}\Delta\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_cuts_bb_Delta_pT
+Title= $\Delta{p_T}$ between b-jet pair (anti-$k_t$ $R=0.4$ jets with cuts)
+XLabel= $\Delta{p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}\Delta{p_T}$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_cuts_bb_Delta_R
+Title= $\Delta{R}$ between b-jet pair (anti-$k_t$ $R=0.4$ jets with cuts)
+XLabel= $\Delta{R}$
+YLabel= $\text{d}\sigma/\text{d}\Delta{R}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_cuts_b_jet_eta
+Title= $\eta$ of b-jet (anti-$k_t$ $R=0.4$ jet with cuts)
+XLabel= $\eta$
+YLabel= $\text{d}\sigma/\text{d}\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_cuts_b_jet_phi
+Title= $\phi$ of b-jet (anti-$k_t$ $R=0.4$ jet with cuts)
+XLabel= $\phi$
+YLabel= $\text{d}\sigma/\text{d}\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_cuts_b_jet_pT
+Title= ${p_T}$ of b-jet (anti-$k_t$ $R=0.4$ jet with cuts)
+XLabel= ${p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}p_T$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_cuts_H_eta_using_bb
+Title= $\eta$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.4$ jets with cuts)
+XLabel= $\eta$
+YLabel= $\text{d}\sigma/\text{d}\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_cuts_H_mass_using_bb
+Title= $m$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.4$ jets)
+XLabel= $m [GeV]$
+YLabel= $\text{d}\sigma/\text{d}m$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_cuts_H_phi_using_bb
+Title= $\phi$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.4$ jets with cuts)
+XLabel= $\phi$
+YLabel= $\text{d}\sigma/\text{d}\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_cuts_H_pT_using_bb
+Title= ${p_T}$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.4$ jets with cuts)
+XLabel= $p_T [GeV]$
+YLabel= $\text{d}\sigma/\text{d}p_T$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_cuts_VBbb_Delta_eta
+Title= $\Delta\eta$ between Higgs candidate and W (anti-$k_t$ $R=0.4$ jets with cuts)
+XLabel= $\Delta\eta$
+YLabel= $\text{d}\sigma/\text{d}\Delta\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_cuts_VBbb_Delta_phi
+Title= $\Delta\phi$ between Higgs candidate and W (anti-$k_t$ $R=0.4$ jets with cuts)
+XLabel= $\Delta\phi$
+YLabel= $\text{d}\sigma/\text{d}\Delta\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_cuts_VBbb_Delta_pT
+Title= $\Delta{p_T}$ between Higgs candidate and W (anti-$k_t$ $R=0.4$ jets with cuts)
+XLabel= $\Delta{p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}\Delta{p_T}$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_cuts_VBbb_Delta_R
+Title= $\Delta{R}$ between Higgs candidate and W (anti-$k_t$ $R=0.4$ jets with cuts)
+XLabel= $\Delta{R}$
+YLabel= $\text{d}\sigma/\text{d}\Delta{R}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_cuts_b_jet_multiplicity
+Title= b-jet multiplicity (anti-$k_t$ $R=0.4$ jet with cuts)
+XLabel= $N_{b-jet}$
+YLabel= $\text{d}\sigma/\text{d}N_{b-jet}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_bb_Delta_eta
+Title= $\Delta\eta$ between b-jet pair (anti-$k_t$ $R=0.5$ jets without cuts)
+XLabel= $\Delta\eta$
+YLabel= $\text{d}\sigma/\text{d}\Delta\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_bb_Delta_phi
+Title= $\Delta\phi$ between b-jet pair (anti-$k_t$ $R=0.5$ jets without cuts)
+XLabel= $\Delta\phi$
+YLabel= $\text{d}\sigma/\text{d}\Delta\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_bb_Delta_pT
+Title= $\Delta{p_T}$ between b-jet pair (anti-$k_t$ $R=0.5$ jets without cuts)
+XLabel= $\Delta{p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}\Delta{p_T}$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_bb_Delta_R
+Title= $\Delta{R}$ between b-jet pair (anti-$k_t$ $R=0.5$ jets without cuts)
+XLabel= $\Delta{R}$
+YLabel= $\text{d}\sigma/\text{d}\Delta{R}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_b_jet_eta
+Title= $\eta$ of b-jet (anti-$k_t$ $R=0.5$ jet without cuts)
+XLabel= $\eta$
+YLabel= $\text{d}\sigma/\text{d}\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_b_jet_multiplicity
+Title= b-jet multiplicity (anti-$k_t$ $R=0.5$ jet without cuts)
+XLabel= $N_{b-jet}$
+YLabel= $\text{d}\sigma/\text{d}N_{b-jet}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_b_jet_phi
+Title= $\phi$ of b-jet (anti-$k_t$ $R=0.5$ jet without cuts)
+XLabel= $\phi$
+YLabel= $\text{d}\sigma/\text{d}\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_b_jet_pT
+Title= ${p_T}$ of b-jet (anti-$k_t$ $R=0.5$ jet without cuts)
+XLabel= ${p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}p_T$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_H_eta_using_bb
+Title= $\eta$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.5$ jets without cuts)
+XLabel= $\eta$
+YLabel= $\text{d}\sigma/\text{d}\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_H_mass_using_bb
+Title= $m$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.5$ jets without cuts)
+XLabel= $m [GeV]$
+YLabel= $\text{d}\sigma/\text{d}m$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_H_phi_using_bb
+Title= $\phi$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.5$ jets without cuts)
+XLabel= $\phi$
+YLabel= $\text{d}\sigma/\text{d}\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_H_pT_using_bb
+Title= ${p_T}$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.5$ jets without cuts)
+XLabel= $p_T [GeV]$
+YLabel= $\text{d}\sigma/\text{d}p_T$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_jet_eta
+Title= $\eta$ of jet (anti-$k_t$ $R=0.5$ jet without cuts)
+XLabel= $\eta$
+YLabel= $\text{d}\sigma/\text{d}\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_jet_multiplicity
+Title= jet multiplicity (anti-$k_t$ $R=0.5$ jet without cuts)
+XLabel= $N_{jet}$
+YLabel= $\text{d}\sigma/\text{d}N_{jet}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_jet_phi
+Title= $\phi$ of jet (anti-$k_t$ $R=0.5$ jet without cuts)
+XLabel= $\phi$
+YLabel= $\text{d}\sigma/\text{d}\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_jet_pT
+Title= ${p_T}$ of jet (anti-$k_t$ $R=0.5$ jet without cuts)
+XLabel= ${p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}p_T$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_VBbb_Delta_eta
+Title= $\Delta\eta$ between Higgs candidate and W (anti-$k_t$ $R=0.5$ jets without cuts)
+XLabel= $\Delta\eta$
+YLabel= $\text{d}\sigma/\text{d}\Delta\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_VBbb_Delta_phi
+Title= $\Delta\phi$ between Higgs candidate and W (anti-$k_t$ $R=0.5$ jets without cuts)
+XLabel= $\Delta\phi$
+YLabel= $\text{d}\sigma/\text{d\Delta\phi}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_VBbb_Delta_pT
+Title= $\Delta{p_T}$ between Higgs candidate and W (anti-$k_t$ $R=0.5$ jets without cuts)
+XLabel= $\Delta{p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}\Delta{p_T}$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_VBbb_Delta_R
+Title= $\Delta{R}$ between Higgs candidate and W (anti-$k_t$ $R=0.5$ jets without cuts)
+XLabel= $\Delta{R}$
+YLabel= $\text{d}\sigma/\text{d}\Delta{R}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_cuts_bb_Delta_eta
+Title= $\Delta\eta$ between b-jet pair (anti-$k_t$ $R=0.5$ jets with cuts)
+XLabel= $\Delta\eta$
+YLabel= $\text{d}\sigma/\text{d}\Delta\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_cuts_bb_Delta_phi
+Title= $\Delta\phi$ between b-jet pair (anti-$k_t$ $R=0.5$ jets with cuts)
+XLabel= $\Delta\phi$
+YLabel= $\text{d}\sigma/\text{d}\Delta\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_cuts_bb_Delta_pT
+Title= $\Delta{p_T}$ between b-jet pair (anti-$k_t$ $R=0.5$ jets with cuts)
+XLabel= $\Delta{p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}\Delta{p_T}$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_cuts_bb_Delta_R
+Title= $\Delta{R}$ between b-jet pair (anti-$k_t$ $R=0.5$ jets with cuts)
+XLabel= $\Delta{R}$
+YLabel= $\text{d}\sigma/\text{d}\Delta{R}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_cuts_b_jet_eta
+Title= $\eta$ of b-jet (anti-$k_t$ $R=0.5$ jet with cuts)
+XLabel= $\eta$
+YLabel= $\text{d}\sigma/\text{d}\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_cuts_b_jet_phi
+Title= $\phi$ of b-jet (anti-$k_t$ $R=0.5$ jet with cuts)
+XLabel= $\phi$
+YLabel= $\text{d}\sigma/\text{d}\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_cuts_b_jet_pT
+Title= ${p_T}$ of b-jet (anti-$k_t$ $R=0.5$ jet with cuts)
+XLabel= ${p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}p_T$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_cuts_H_eta_using_bb
+Title= $\eta$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.5$ jets with cuts)
+XLabel= $\eta$
+YLabel= $\text{d}\sigma/\text{d}\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_cuts_H_mass_using_bb
+Title= $m$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.5$ jets with cuts)
+XLabel= $m [GeV]$
+YLabel= $\text{d}\sigma/\text{d}m$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_cuts_H_phi_using_bb
+Title= $\phi$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.5$ jets with cuts)
+XLabel= $\phi$
+YLabel= $\text{d}\sigma/\text{d}\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_cuts_H_pT_using_bb
+Title= ${p_T}$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.5$ jets with cuts)
+XLabel= $p_T [GeV]$
+YLabel= $\text{d}\sigma/\text{d}p_T$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_cuts_VBbb_Delta_eta
+Title= $\Delta\eta$ between Higgs candidate and W (anti-$k_t$ $R=0.5$ jets with cuts)
+XLabel= $\Delta\eta$
+YLabel= $\text{d}\sigma/\text{d}\Delta\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_cuts_VBbb_Delta_phi
+Title= $\Delta\phi$ between Higgs candidate and W (anti-$k_t$ $R=0.5$ jets with cuts)
+XLabel= $\Delta\phi$
+YLabel= $\text{d}\sigma/\text{d}\Delta\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_cuts_VBbb_Delta_pT
+Title= $\Delta{p_T}$ between Higgs candidate and W (anti-$k_t$ $R=0.5$ jets with cuts)
+XLabel= $\Delta{p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}\Delta{p_T}$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_cuts_VBbb_Delta_R
+Title= $\Delta{R}$ between Higgs candidate and W (anti-$k_t$ $R=0.5$ jets with cuts)
+XLabel= $\Delta{R}$
+YLabel= $\text{d}\sigma/\text{d}\Delta{R}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_cuts_b_jet_multiplicity
+Title= b-jet multiplicity (anti-$k_t$ $R=0.5$ jet with cuts)
+XLabel= $N_{b-jet}$
+YLabel= $\text{d}\sigma/\text{d}N_{b-jet}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_bb_Delta_eta
+Title= $\Delta\eta$ between b-jet pair (anti-$k_t$ $R=0.6$ jets without cuts)
+XLabel= $\Delta\eta$
+YLabel= $\text{d}\sigma/\text{d}\Delta\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_bb_Delta_phi
+Title= $\Delta\phi$ between b-jet pair (anti-$k_t$ $R=0.6$ jets without cuts)
+XLabel= $\Delta\phi$
+YLabel= $\text{d}\sigma/\text{d}\Delta\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_bb_Delta_pT
+Title= $\Delta{p_T}$ between b-jet pair (anti-$k_t$ $R=0.6$ jets without cuts)
+XLabel= $\Delta{p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}\Delta{p_T}$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_bb_Delta_R
+Title= $\Delta{R}$ between b-jet pair (anti-$k_t$ $R=0.6$ jets without cuts)
+XLabel= $\Delta{R}$
+YLabel= $\text{d}\sigma/\text{d}\Delta{R}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_b_jet_eta
+Title= $\eta$ of b-jet (anti-$k_t$ $R=0.6$ jet without cuts)
+XLabel= $\eta$
+YLabel= $\text{d}\sigma/\text{d}\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_b_jet_multiplicity
+Title= b-jet multiplicity (anti-$k_t$ $R=0.6$ jet without cuts)
+XLabel= $N_{jet}$
+YLabel= $\text{d}\sigma/\text{d}N_{jet}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_b_jet_phi
+Title= $\phi$ of b-jet (anti-$k_t$ $R=0.6$ jet without cuts)
+XLabel= $\phi$
+YLabel= $\text{d}\sigma/\text{d}\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_b_jet_pT
+Title= ${p_T}$ of b-jet (anti-$k_t$ $R=0.6$ jet without cuts)
+XLabel= ${p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}p_T$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_H_eta_using_bb
+Title= $\eta$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.6$ jets without cuts)
+XLabel= $\eta$
+YLabel= $\text{d}\sigma/\text{d}\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_H_mass_using_bb
+Title= $m$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.6$ jets without cuts)
+XLabel= $m [GeV]$
+YLabel= $\text{d}\sigma/\text{d}m$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_H_phi_using_bb
+Title= $\phi$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.6$ jets without cuts)
+XLabel= $\phi$
+YLabel= $\text{d}\sigma/\text{d}\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_H_pT_using_bb
+Title= ${p_T}$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.6$ jets without cuts)
+XLabel= $p_T [GeV]$
+YLabel= $\text{d}\sigma/\text{d}p_T$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_jet_eta
+Title= $\eta$ of jet (anti-$k_t$ $R=0.6$ jet without cuts)
+XLabel= $\eta$
+YLabel= $\text{d}\sigma/\text{d}\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_jet_multiplicity
+Title= jet multiplicity (anti-$k_t$ $R=0.6$ jet without cuts)
+XLabel= $N_{jet}$
+YLabel= $\text{d}\sigma/\text{d}N_{jet}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_jet_phi
+Title= $\phi$ of jet (anti-$k_t$ $R=0.6$ jet without cuts)
+XLabel= $\phi$
+YLabel= $\text{d}\sigma/\text{d}\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_jet_pT
+Title= ${p_T}$ of jet (anti-$k_t$ $R=0.6$ jet without cuts)
+XLabel= ${p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}p_T$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_VBbb_Delta_eta
+Title= $\Delta\eta$ between Higgs candidate and W (anti-$k_t$ $R=0.6$ jets without cuts)
+XLabel= $\Delta\eta$
+YLabel= $\text{d}\sigma/\text{d}\Delta\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_VBbb_Delta_phi
+Title= $\Delta\phi$ between Higgs candidate and W (anti-$k_t$ $R=0.6$ jets without cuts)
+XLabel= $\Delta\phi$
+YLabel= $\text{d}\sigma/\text{d}\Delta\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_VBbb_Delta_pT
+Title= $\Delta{p_T}$ between Higgs candidate and W (anti-$k_t$ $R=0.6$ jets without cuts)
+XLabel= $\Delta{p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}\Delta{p_T}$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_VBbb_Delta_R
+Title= $\Delta{R}$ between Higgs candidate and W (anti-$k_t$ $R=0.6$ jets without cuts)
+XLabel= $\Delta{R}$
+YLabel= $\text{d}\sigma/\text{d}\Delta{R}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_cuts_bb_Delta_eta
+Title= $\Delta\eta$ between b-jet pair (anti-$k_t$ $R=0.6$ jets with cuts)
+XLabel= $\Delta\eta$
+YLabel= $\text{d}\sigma/\text{d}\Delta\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_cuts_bb_Delta_phi
+Title= $\Delta\phi$ between b-jet pair (anti-$k_t$ $R=0.6$ jets with cuts)
+XLabel= $\Delta\phi$
+YLabel= $\text{d}\sigma/\text{d}\Delta\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_cuts_bb_Delta_pT
+Title= $\Delta{p_T}$ between b-jet pair (anti-$k_t$ $R=0.6$ jets with cuts)
+XLabel= $\Delta{p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}\Delta{p_T}$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_cuts_bb_Delta_R
+Title= $\Delta{R}$ between b-jet pair (anti-$k_t$ $R=0.6$ jets with cuts)
+XLabel= $\Delta{R}$
+YLabel= $\text{d}\sigma/\text{d}\Delta{R}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_cuts_b_jet_eta
+Title= $\eta$ of b-jet (anti-$k_t$ $R=0.6$ jet with cuts
+XLabel= $\eta$
+YLabel= $\text{d}\sigma/\text{d}\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_cuts_b_jet_phi
+Title= $\phi$ of b-jet (anti-$k_t$ $R=0.6$ jet with cuts)
+XLabel= $\phi$
+YLabel= $\text{d}\sigma/\text{d}\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_cuts_b_jet_pT
+Title= ${p_T}$ of b-jet (anti-$k_t$ $R=0.6$ jet with cuts)
+XLabel= ${p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}p_T$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_cuts_H_eta_using_bb
+Title= $\eta$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.6$ jets with cuts)
+XLabel= $\eta$
+YLabel= $\text{d}\sigma/\text{d}\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_cuts_H_mass_using_bb
+Title= $m$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.6$ jets with cuts)
+XLabel= $m [GeV]$
+YLabel= $\text{d}\sigma/\text{d}m$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_cuts_H_phi_using_bb
+Title= $\phi$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.6$ jets with cuts)
+XLabel= $\phi$
+YLabel= $\text{d}\sigma/\text{d}\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_cuts_H_pT_using_bb
+Title= ${p_T}$ of Higgs candidate from b-jet pair (anti-$k_t$ $R=0.6$ jets with cuts)
+XLabel= $p_T [GeV]$
+YLabel= $\text{d}\sigma/\text{d}p_T$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_cuts_VBbb_Delta_eta
+Title= $\Delta\eta$ between Higgs candidate and W (anti-$k_t$ $R=0.6$ jets with cuts)
+XLabel= $\Delta\eta$
+YLabel= $\text{d}\sigma/\text{d}\Delta\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_cuts_VBbb_Delta_phi
+Title= $\Delta\phi$ between Higgs candidate and W (anti-$k_t$ $R=0.6$ jets with cuts)
+XLabel= $\Delta\phi$
+YLabel= $\text{d}\sigma/\text{d}\Delta\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_cuts_VBbb_Delta_pT
+Title= $\Delta{p_T}$ between Higgs candidate and W (anti-$k_t$ $R=0.6$ jets with cuts)
+XLabel= $\Delta{p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}\Delta{p_T}$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_cuts_VBbb_Delta_R
+Title= $\Delta{R}$ between Higgs candidate and W (anti-$k_t$ $R=0.6$ jets with cuts)
+XLabel= $\Delta{R}$
+YLabel= $\text{d}\sigma/\text{d}\Delta{R}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_cuts_b_jet_multiplicity
+Title= b-jet multiplicity (anti-$k_t$ $R=0.6$ jet with cuts)
+XLabel= $N_{b-jet}$
+YLabel= $\text{d}\sigma/\text{d}N_{b-jet}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/VB_eta
+Title= $\eta$ of vector boson
+XLabel= $\eta$
+YLabel= $\text{d}\sigma/\text{d}\eta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/VB_mass
+Title= $m$ of vector boson
+XLabel= $m$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}m$ [pb/GeV]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/Z_multiplicity
+Title= Z boson multiplicity
+XLabel= $N_{Z}$
+YLabel= $\text{d}\sigma/\text{d}N_{Z}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/W_multiplicity
+Title= W boson multiplicity
+XLabel= $N_{W}$
+YLabel= $\text{d}\sigma/\text{d}N_{W}$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/VB_phi
+Title= $\phi$ of W
+XLabel= $\phi$
+YLabel= $\text{d}\sigma/\text{d}\phi$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/VB_pT
+Title= ${p_T}$ of W
+XLabel= ${p_T}$ [GeV]
+YLabel= $\text{d}\sigma/\text{d}{p_T}$ [pb/GeV]
+LogY=0
+# END PLOT
+
+
+# BEGIN PLOT /MC_VH2BB/pure_bVB_angle_Hframe
+Title= $\theta_{MIN}$ between b-hadron and W, in Higgs rest frame (without cuts)
+XLabel= $\theta$
+YLabel= $\text{d}\sigma/\text{d}\theta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_cuts_bVB_angle_Hframe
+Title= $\theta_{MIN}$ between b-hadron and W, in Higgs rest frame (with cuts)
+XLabel= $\theta$
+YLabel= $\text{d}\sigma/\text{d}\theta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_bVB_angle_Hframe
+Title= $\theta_{MIN}$ between b-jet and W, in Higgs rest frame (anti-$k_t$ $R=0.4$ jets without cuts)
+XLabel= $\theta$
+YLabel= $\text{d}\sigma/\text{d}\theta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_cuts_bVB_angle_Hframe
+Title= $\theta_{MIN}$ between b-jet and W, in Higgs rest frame (anti-$k_t$ $R=0.4$ jets with cuts)
+XLabel= $\theta$
+YLabel= $\text{d}\sigma/\text{d}\theta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_bVB_angle_Hframe
+Title= $\theta_{MIN}$ between b-jet and W, in Higgs rest frame (anti-$k_t$ $R=0.5$ jets without cuts)
+XLabel= $\theta$
+YLabel= $\text{d}\sigma/\text{d}\theta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_cuts_bVB_angle_Hframe
+Title= $\theta_{MIN}$ between b-jet and W, in Higgs rest frame (anti-$k_t$ $R=0.5$ jets with cuts)
+XLabel= $\theta$
+YLabel= $\text{d}\sigma/\text{d}\theta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_bVB_angle_Hframe
+Title= $\theta_{MIN}$ between b-jet and W, in Higgs rest frame (anti-$k_t$ $R=0.6$ jets without cuts)
+XLabel= $\theta$
+YLabel= $\text{d}\sigma/\text{d}\theta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_cuts_bVB_angle_Hframe
+Title= $\theta_{MIN}$ between b-jet and W, in Higgs rest frame (anti-$k_t$ $R=0.6$ jets with cuts)
+XLabel= $\theta$
+YLabel= $\text{d}\sigma/\text{d}\theta$ [pb]
+LogY=0
+# END PLOT
+
+
+# BEGIN PLOT /MC_VH2BB/pure_bb_angle_Hframe
+Title= $\theta$ between b-hadrons, in Higgs rest frame (without cuts)
+XLabel= $\theta$
+YLabel= $\text{d}\sigma/\text{d}\theta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_cuts_bb_angle_Hframe
+Title= $\theta$ between b-hadrons, in Higgs rest frame (with cuts)
+XLabel= $\theta$
+YLabel= $\text{d}\sigma/\text{d}\theta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_bb_angle_Hframe
+Title= $\theta$ between b-jets, in Higgs rest frame (anti-$k_t$ $R=0.4$ jets without cuts)
+XLabel= $\theta$
+YLabel= $\text{d}\sigma/\text{d}\theta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_cuts_bb_angle_Hframe
+Title= $\theta$ between b-jets, in Higgs rest frame (anti-$k_t$ $R=0.4$ jets with cuts)
+XLabel= $\theta$
+YLabel= $\text{d}\sigma/\text{d}\theta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_bb_angle_Hframe
+Title= $\theta$ between b-jets, in Higgs rest frame (anti-$k_t$ $R=0.5$ jets without cuts)
+XLabel= $\theta$
+YLabel= $\text{d}\sigma/\text{d}\theta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_cuts_bb_angle_Hframe
+Title= $\theta$ between b-jets, in Higgs rest frame (anti-$k_t$ $R=0.5$ jets with cuts)
+XLabel= $\theta$
+YLabel= $\text{d}\sigma/\text{d}\theta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_bb_angle_Hframe
+Title= $\theta$ between b-jets, in Higgs rest frame (anti-$k_t$ $R=0.6$ jets without cuts)
+XLabel= $\theta$
+YLabel= $\text{d}\sigma/\text{d}\theta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_cuts_bb_angle_Hframe
+Title= $\theta$ between b-jets, in Higgs rest frame (anti-$k_t$ $R=0.6$ jets with cuts)
+XLabel= $\theta$
+YLabel= $\text{d}\sigma/\text{d}\theta$ [pb]
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R04_cuts_bb_deltaR_v_HpT
+Title= $\Delta{R}$ between b-jets against $p_T$ of Higgs candidate (anti-$k_t$ $R=0.4$ jets with cuts)
+XLabel= $p_T$ [GeV]
+YLabel= $\Delta{R}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R05_cuts_bb_deltaR_v_HpT
+Title= $\Delta{R}$ between b-jets against $p_T$ of Higgs candidate (anti-$k_t$ $R=0.5$ jets with cuts)
+XLabel= $p_T$ [GeV]
+YLabel= $\Delta{R}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/R06_cuts_bb_deltaR_v_HpT
+Title= $\Delta{R}$ between b-jets against $p_T$ of Higgs candidate (anti-$k_t$ $R=0.6$ jets with cuts)
+XLabel= $p_T$ [GeV]
+YLabel= $\Delta{R}$
+LogY=0
+# END PLOT
+
+# BEGIN PLOT /MC_VH2BB/pure_cuts_bb_deltaR_v_HpT
+Title= $\Delta{R}$ between b-hadrons against $p_T$ of Higgs candidate (with cuts)
+XLabel= $p_T$ [GeV]
+YLabel= $\Delta{R}$
+LogY=0
+# END PLOT
+
+

Modified: trunk/data/plotinfo/Makefile.am
==============================================================================
--- trunk/data/plotinfo/Makefile.am	Tue Sep 13 18:00:19 2011	(r3350)
+++ trunk/data/plotinfo/Makefile.am	Wed Sep 14 09:55:10 2011	(r3351)
@@ -103,6 +103,7 @@
   MC_TTBAR.plot \
   MC_JETS.plot \
   MC_PHOTONJETS.plot \
+  MC_VH2BB.plot \
   MC_HJETS.plot \
   MC_WJETS.plot \
   MC_WPOL.plot \

Added: trunk/src/Analyses/MC_VH2BB.cc
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ trunk/src/Analyses/MC_VH2BB.cc	Wed Sep 14 09:55:10 2011	(r3351)
@@ -0,0 +1,298 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/RivetAIDA.hh"
+#include "Rivet/Tools/Logging.hh"
+#include "Rivet/Projections/FinalState.hh"
+#include "Rivet/Tools/ParticleIdUtils.hh"
+#include "Rivet/Projections/ZFinder.hh"
+#include "Rivet/Projections/WFinder.hh"
+#include "Rivet/Projections/UnstableFinalState.hh"
+#include "Rivet/Projections/FastJets.hh"
+#include "Rivet/Math/LorentzTrans.hh"
+#include "Rivet/Math/Constants.hh"
+#include <cmath>
+#include <vector>
+
+namespace Rivet {
+
+
+  class MC_VH2BB : public Analysis {
+  public:
+
+    /// @name Constructors etc.
+    //@{
+
+    /// Constructor
+    MC_VH2BB()
+      : Analysis("MC_VH2BB")
+    {
+      setNeedsCrossSection(true);
+    }
+
+    //@}
+
+
+  public:
+
+    /// @name Analysis methods
+    //@{
+
+    vector<double> boostAngles(const FourMomentum& b1, const FourMomentum& b2, const FourMomentum& vb){
+
+      // This should take in the four-momenta of two b's (jets/hadrons) and a vector boson, for the process VB*->VBH with H->bb
+      // It should return the smallest angle between the virtual vector boson and one of the b's, in the rest frame of the Higgs boson.
+      // It should also return (as the second element of the vector) the angle between the b's, in the rest frame of the Higgs boson.
+
+      FourMomentum higgsMomentum = b1 + b2;
+      FourMomentum virtualVBMomentum = higgsMomentum + vb;
+
+      LorentzTransform lt( -higgsMomentum.boostVector() );
+
+      FourMomentum virtualVBMomentumBOOSTED = lt.transform(virtualVBMomentum);
+      FourMomentum b1BOOSTED = lt.transform(b1);
+      FourMomentum b2BOOSTED = lt.transform(b2);
+
+      double angle1 = b1BOOSTED.angle(virtualVBMomentumBOOSTED);
+      double angle2 = b2BOOSTED.angle(virtualVBMomentumBOOSTED);
+
+      double anglebb = b1BOOSTED.angle(b2BOOSTED);
+
+      vector<double> toReturn;
+      toReturn.push_back(angle1 < angle2 ? angle1 : angle2);
+      toReturn.push_back(anglebb);
+
+      return toReturn;
+    }
+
+
+    /// @todo Replace with angle range mapper functions
+    double correctPhi(double phi) {
+      bool flag = false;
+
+      while(flag == false){
+        if(phi > M_PI){
+          phi = phi - (2*M_PI);
+        } else {
+          flag = true;
+        }
+      }
+
+      flag = false;
+      while (flag == false) {
+        if (phi <= -M_PI) {
+          phi = phi + (2*M_PI);
+        } else {
+          flag = true;
+        }
+      }
+
+      return phi;
+    }
+
+
+    /// Book histograms and initialise projections before the run
+    void init() {
+
+      ZFinder zeefinder(-3.5, 3.5, 25.0*GeV, ELECTRON, 65.0*GeV, 115.0*GeV, 0.2, true, true);
+      addProjection(zeefinder, "ZeeFinder");
+      ZFinder zmmfinder(-3.5, 3.5, 25.0*GeV, MUON, 65.0*GeV, 115.0*GeV, 0.2, true, true);
+      addProjection(zmmfinder, "ZmmFinder");
+
+      WFinder wefinder(-3.5, 3.5, 25.0*GeV, ELECTRON, 60.0*GeV, 100.0*GeV, 25.0*GeV, 0.2);
+      addProjection(wefinder, "WeFinder");
+      WFinder wmfinder(-3.5, 3.5, 25.0*GeV, MUON, 60.0*GeV, 100.0*GeV, 25.0*GeV, 0.2);
+      addProjection(wmfinder, "WmFinder");
+
+      FinalState fs;
+      addProjection(fs, "FinalState");
+      addProjection(FastJets(fs, FastJets::ANTIKT, 0.4), "AntiKT04");
+      addProjection(FastJets(fs, FastJets::ANTIKT, 0.5), "AntiKT05");
+      addProjection(FastJets(fs, FastJets::ANTIKT, 0.6), "AntiKT06");
+
+      /// Book histograms
+      _h_jet_bb_Delta_eta = bookHistogram1D("jet_bb_Delta_eta", 50, 0, 4);
+      _h_jet_bb_Delta_phi = bookHistogram1D("jet_bb_Delta_phi", 50, 0, M_PI);
+      _h_jet_bb_Delta_pT = bookHistogram1D("jet_bb_Delta_pT", 50,0, 500);
+      _h_jet_bb_Delta_R = bookHistogram1D("jet_bb_Delta_R", 50, 0, 5);
+      _h_jet_b_jet_eta = bookHistogram1D("jet_b_jet_eta", 50, -4, 4);
+      _h_jet_b_jet_multiplicity = bookHistogram1D("jet_b_jet_multiplicity", 11, -0.5, 10.5);
+      _h_jet_b_jet_phi = bookHistogram1D("jet_b_jet_phi", 50, -M_PI, M_PI);
+      _h_jet_b_jet_pT = bookHistogram1D("jet_b_jet_pT", 50, 0, 500);
+      _h_jet_H_eta_using_bb = bookHistogram1D("jet_H_eta_using_bb", 50, -4, 4);
+      _h_jet_H_mass_using_bb = bookHistogram1D("jet_H_mass_using_bb", 50, 50, 200);
+      _h_jet_H_phi_using_bb = bookHistogram1D("jet_H_phi_using_bb", 50, -M_PI, M_PI);
+      _h_jet_H_pT_using_bb = bookHistogram1D("jet_H_pT_using_bb", 50, 0, 500);
+      _h_jet_eta = bookHistogram1D("jet_eta", 50, -4, 4);
+      _h_jet_multiplicity = bookHistogram1D("jet_multiplicity", 11, -0.5, 10.5);
+      _h_jet_phi = bookHistogram1D("jet_phi", 50, -M_PI, M_PI);
+      _h_jet_pT = bookHistogram1D("jet_pT", 50, 0, 500);
+      _h_jet_VBbb_Delta_eta = bookHistogram1D("jet_VBbb_Delta_eta", 50, 0, 4);
+      _h_jet_VBbb_Delta_phi = bookHistogram1D("jet_VBbb_Delta_phi", 50, 0, M_PI);
+      _h_jet_VBbb_Delta_pT = bookHistogram1D("jet_VBbb_Delta_pT", 50, 0, 500);
+      _h_jet_VBbb_Delta_R = bookHistogram1D("jet_VBbb_Delta_R", 50, 0, 8);
+
+      _h_VB_eta = bookHistogram1D("VB_eta", 50, -4, 4);
+      _h_VB_mass = bookHistogram1D("VB_mass", 50, 60, 110);
+      _h_Z_multiplicity = bookHistogram1D("Z_multiplicity", 11, -0.5, 10.5);
+      _h_W_multiplicity = bookHistogram1D("W_multiplicity", 11, -0.5, 10.5);
+      _h_VB_phi = bookHistogram1D("VB_phi", 50, -M_PI, M_PI);
+      _h_VB_pT = bookHistogram1D("VB_pT", 50, 0, 500);
+
+      _h_jet_bVB_angle_Hframe = bookHistogram1D("jet_bVB_angle_Hframe", 50, 0, M_PI);
+      _h_jet_bb_angle_Hframe = bookHistogram1D("jet_bb_angle_Hframe", 50, 0, M_PI);
+    }
+
+
+    /// Perform the per-event analysis
+    void analyze(const Event& event) {
+      const double weight = event.weight();
+
+      const double JETPTCUT = 30*GeV;
+
+      const ZFinder& zeefinder = applyProjection<ZFinder>(event, "ZeeFinder");
+      const ZFinder& zmmfinder = applyProjection<ZFinder>(event, "ZmmFinder");
+
+      const WFinder& wefinder = applyProjection<WFinder>(event, "WeFinder");
+      const WFinder& wmfinder = applyProjection<WFinder>(event, "WmFinder");
+
+      Jets jets = applyProjection<FastJets>(event, "AntiKTJets").jetsByPt(JETPTCUT);
+
+      ParticleVector vectorBosons = zeefinder.particles();
+      /// @todo Don't we have a neater vector concatenation?
+      vectorBosons.insert(vectorBosons.end(), zeefinder.particles().begin(), zeefinder.particles().end());
+      vectorBosons.insert(vectorBosons.end(), wefinder.particles().begin(), wefinder.particles().end());
+      vectorBosons.insert(vectorBosons.end(), wmfinder.particles().begin(), wmfinder.particles().end());
+
+      _h_Z_multiplicity->fill(zeefinder.particles().size() + zmmfinder.particles().size(), weight);
+      _h_W_multiplicity->fill(wefinder.particles().size() + wmfinder.particles().size(), weight);
+      _h_jet_multiplicity->fill(jets.size(), weight);
+
+      // Identify the b-jets
+      Jets bjets;
+      foreach (const Jet& jet, jets) {
+        const double jetEta = jet.momentum().eta();
+        const double jetPhi = correctPhi(jet.momentum().phi());
+        const double jetPt = jet.momentum().pT();
+        _h_jet_eta->fill(jetEta, weight);
+        _h_jet_phi->fill(jetPhi, weight);
+        _h_jet_pT->fill(jetPt/GeV, weight);
+
+        if (jet.containsBottom() && jet.momentum().pT() > JETPTCUT) {
+          bjets.push_back(jet);
+          _h_jet_b_jet_eta->fill( jetEta , weight );
+          _h_jet_b_jet_phi->fill( jetPhi , weight );
+          _h_jet_b_jet_pT->fill( jetPt , weight );
+        }
+      }
+      _h_jet_b_jet_multiplicity->fill(bjets.size(), weight);
+
+      // Plot vector boson properties
+      foreach (const Particle& v, vectorBosons) {
+        _h_VB_phi->fill(v.momentum().phi(), weight);
+        _h_VB_pT->fill(v.momentum().pT(), weight);
+        _h_VB_eta->fill(v.momentum().eta(), weight);
+        _h_VB_mass->fill(v.momentum().mass(), weight);
+      }
+
+      // Construct Higgs candidates from pairs of b-jets
+      for (size_t i = 0; i < bjets.size()-1; ++i) {
+        for (size_t j = i+1; j < bjets.size(); ++j) {
+          const Jet& jet1 = bjets[i];
+          const Jet& jet2 = bjets[j];
+
+          const double deltaEtaJJ = fabs(jet1.momentum().eta() - jet2.momentum().eta());
+          const double deltaPhiJJ = deltaPhi(jet1.momentum(), jet2.momentum());
+          const double deltaRJJ = deltaR(jet1.momentum(), jet2.momentum());
+          const double deltaPtJJ = fabs(jet1.momentum().pT() - jet2.momentum().pT());
+          _h_jet_bb_Delta_eta->fill(deltaEtaJJ, weight);
+          _h_jet_bb_Delta_phi->fill(deltaPhiJJ, weight);
+          _h_jet_bb_Delta_pT->fill(deltaPtJJ, weight);
+          _h_jet_bb_Delta_R->fill(deltaRJJ, weight);
+
+          const FourMomentum phiggs = jet1.momentum() + jet2.momentum();
+          _h_jet_H_eta_using_bb->fill(phiggs.eta(), weight);
+          _h_jet_H_mass_using_bb->fill(phiggs.mass(), weight);
+          _h_jet_H_phi_using_bb->fill(phiggs.phi(), weight);
+          _h_jet_H_pT_using_bb->fill(phiggs.pT(), weight);
+
+          foreach (const Particle& v, vectorBosons) {
+            const double deltaEtaVH = fabs(phiggs.eta() - v.momentum().eta());
+            const double deltaPhiVH = deltaPhi(phiggs, v.momentum());
+            const double deltaRVH = deltaR(phiggs, v.momentum());
+            const double deltaPtVH = fabs(phiggs.pT() - v.momentum().pT());
+            _h_jet_VBbb_Delta_eta->fill(deltaEtaVH, weight);
+            _h_jet_VBbb_Delta_phi->fill(deltaPhiVH, weight);
+            _h_jet_VBbb_Delta_pT->fill(deltaPtVH, weight);
+            _h_jet_VBbb_Delta_R->fill(deltaRVH, weight);
+
+            // Calculate boost angles
+            const vector<double> boostAngles = angularCookbook(jet1.momentum(), jet2.momentum(), v.momentum());
+            _h_jet_bVB_angle_Hframe->fill(boostAngles[0], weight);
+            _h_jet_bb_angle_Hframe->fill(boostAngles[1], weight);
+          }
+        }
+      }
+    }
+
+
+    /// Normalise histograms etc., after the run
+    void finalize() {
+      scale(_h_jet_bb_Delta_eta, crossSection()/sumOfWeights());
+      scale(_h_jet_bb_Delta_phi, crossSection()/sumOfWeights());
+      scale(_h_jet_bb_Delta_pT, crossSection()/sumOfWeights());
+      scale(_h_jet_bb_Delta_R, crossSection()/sumOfWeights());
+      scale(_h_jet_b_jet_eta, crossSection()/sumOfWeights());
+      scale(_h_jet_b_jet_multiplicity, crossSection()/sumOfWeights());
+      scale(_h_jet_b_jet_phi, crossSection()/sumOfWeights());
+      scale(_h_jet_b_jet_pT, crossSection()/sumOfWeights());
+      scale(_h_jet_H_eta_using_bb, crossSection()/sumOfWeights());
+      scale(_h_jet_H_mass_using_bb, crossSection()/sumOfWeights());
+      scale(_h_jet_H_phi_using_bb, crossSection()/sumOfWeights());
+      scale(_h_jet_H_pT_using_bb, crossSection()/sumOfWeights());
+      scale(_h_jet_eta, crossSection()/sumOfWeights());
+      scale(_h_jet_multiplicity, crossSection()/sumOfWeights());
+      scale(_h_jet_phi, crossSection()/sumOfWeights());
+      scale(_h_jet_pT, crossSection()/sumOfWeights());
+      scale(_h_jet_VBbb_Delta_eta, crossSection()/sumOfWeights());
+      scale(_h_jet_VBbb_Delta_phi, crossSection()/sumOfWeights());
+      scale(_h_jet_VBbb_Delta_pT, crossSection()/sumOfWeights());
+      scale(_h_jet_VBbb_Delta_R, crossSection()/sumOfWeights());
+
+      scale(_h_VB_eta, crossSection()/sumOfWeights());
+      scale(_h_VB_mass, crossSection()/sumOfWeights());
+      scale(_h_Z_multiplicity, crossSection()/sumOfWeights());
+      scale(_h_W_multiplicity, crossSection()/sumOfWeights());
+      scale(_h_VB_phi, crossSection()/sumOfWeights());
+      scale(_h_VB_pT, crossSection()/sumOfWeights());
+
+      scale(_h_jet_bVB_angle_Hframe, crossSection()/sumOfWeights());
+      scale(_h_jet_bb_angle_Hframe, crossSection()/sumOfWeights());
+    }
+
+    //@}
+
+
+  private:
+
+    /// @name Histograms
+    //@{
+
+    AIDA::IHistogram1D *_h_Z_multiplicity, *_h_W_multiplicity;
+    AIDA::IHistogram1D *_h_jet_bb_Delta_eta, *_h_jet_bb_Delta_phi, *_h_jet_bb_Delta_pT, *_h_jet_bb_Delta_R;
+    AIDA::IHistogram1D *_h_jet_b_jet_eta, *_h_jet_b_jet_multiplicity, *_h_jet_b_jet_phi, *_h_jet_b_jet_pT;
+    AIDA::IHistogram1D *_h_jet_H_eta_using_bb, *_h_jet_H_mass_using_bb, *_h_jet_H_phi_using_bb, *_h_jet_H_pT_using_bb;
+    AIDA::IHistogram1D *_h_jet_eta, *_h_jet_multiplicity, *_h_jet_phi, *_h_jet_pT;
+    AIDA::IHistogram1D *_h_jet_VBbb_Delta_eta, *_h_jet_VBbb_Delta_phi, *_h_jet_VBbb_Delta_pT, *_h_jet_VBbb_Delta_R;
+    AIDA::IHistogram1D *_h_VB_eta, *_h_VB_mass, *_h_VB_phi, *_h_VB_pT;
+    AIDA::IHistogram1D *_h_jet_bVB_angle_Hframe, *_h_jet_bb_angle_Hframe;
+    //AIDA::IProfile1D *_h_jet_cuts_bb_deltaR_v_HpT;
+
+    //@}
+
+  };
+
+
+  // This global object acts as a hook for the plugin system
+  DECLARE_RIVET_PLUGIN(MC_VH2BB);
+
+}

Modified: trunk/src/Analyses/Makefile.am
==============================================================================
--- trunk/src/Analyses/Makefile.am	Tue Sep 13 18:00:19 2011	(r3350)
+++ trunk/src/Analyses/Makefile.am	Wed Sep 14 09:55:10 2011	(r3351)
@@ -243,6 +243,7 @@
     MC_GENERIC.cc \
     MC_JETS.cc \
     MC_PHOTONJETS.cc \
+    MC_VH2BB.cc \
     MC_HJETS.cc \
     MC_WJETS.cc \
     MC_WPOL.cc \


More information about the Rivet-svn mailing list