|
|
[Rivet-svn] r3682 - in trunk: data/anainfo data/plotinfo src/Analyses
blackhole at projects.hepforge.org
blackhole at projects.hepforge.org
Mon Apr 16 12:57:09 BST 2012
Author: richardn
Date: Mon Apr 16 12:57:08 2012
New Revision: 3682
Log:
added ATLAS_2012_CONF_2012_033, latest 0 lepton + met SUSY search
Added:
trunk/data/anainfo/ATLAS_2012_CONF_2012_033.info
- copied, changed from r3651, trunk/data/anainfo/ATLAS_2011_S9212183.info
trunk/data/plotinfo/ATLAS_2012_CONF_2012_033.plot
- copied, changed from r3651, trunk/data/plotinfo/ATLAS_2011_S9225137.plot
trunk/src/Analyses/ATLAS_2012_CONF_2012_033.cc
- copied, changed from r3651, trunk/src/Analyses/ATLAS_2011_S9212183.cc
Modified:
trunk/data/anainfo/Makefile.am
trunk/data/plotinfo/Makefile.am
trunk/src/Analyses/Makefile.am
Copied and modified: trunk/data/anainfo/ATLAS_2012_CONF_2012_033.info (from r3651, trunk/data/anainfo/ATLAS_2011_S9212183.info)
==============================================================================
--- trunk/data/anainfo/ATLAS_2011_S9212183.info Sun Mar 25 11:21:04 2012 (r3651, copy source)
+++ trunk/data/anainfo/ATLAS_2012_CONF_2012_033.info Mon Apr 16 12:57:08 2012 (r3682)
@@ -1,15 +1,16 @@
-Name: ATLAS_2011_S9212183
-Year: 2011
+Name: ATLAS_2012_CONF_2012_033
+Year: 2012
Summary: 0-lepton squark and gluino search
Experiment: ATLAS
Collider: LHC
-SpiresID: 9212183
+SpiresID:
Status: VALIDATED
Authors:
+ - Peter Richardson <Peter.Richardson at durham.ac.uk>
- Chris Wymant <c.m.wymant at durham.ac.uk>
- David Grellscheid <david.grellscheid at durham.ac.uk>
References:
- - arXiv:1102.5290
+ - ATLAS-CONF-2012-033
RunInfo:
BSM signal events at 7000 GeV.
NumEvents: 10000 for BSM signals
@@ -19,17 +20,3 @@
'0-lepton search for squarks and gluinos by ATLAS at 7 TeV.
Event counts in five signal regions are implemented as one-bin
histograms.'
-BibKey: Aad:2011ib
-BibTeX: '@article{Aad:2011ib,
- author = "Aad, Georges and others",
- title = "{Search for squarks and gluinos using final states with
- jets and missing transverse momentum with the ATLAS
- detector in sqrt(s) = 7 TeV proton-proton collisions}",
- collaboration = "ATLAS Collaboration",
- year = "2011",
- note = "* Temporary entry *",
- eprint = "1109.6572",
- archivePrefix = "arXiv",
- primaryClass = "hep-ex",
-}'
-
Modified: trunk/data/anainfo/Makefile.am
==============================================================================
--- trunk/data/anainfo/Makefile.am Sun Apr 15 11:38:39 2012 (r3681)
+++ trunk/data/anainfo/Makefile.am Mon Apr 16 12:57:08 2012 (r3682)
@@ -44,6 +44,7 @@
ATLAS_2012_I1083318.info \
ATLAS_2012_I1084540.info \
ATLAS_2012_I1094568.info \
+ ATLAS_2012_CONF_2012_033.info \
BABAR_2005_S6181155.info \
BABAR_2007_S6895344.info \
BABAR_2007_S7266081.info \
Copied and modified: trunk/data/plotinfo/ATLAS_2012_CONF_2012_033.plot (from r3651, trunk/data/plotinfo/ATLAS_2011_S9225137.plot)
==============================================================================
--- trunk/data/plotinfo/ATLAS_2011_S9225137.plot Sun Mar 25 11:21:04 2012 (r3651, copy source)
+++ trunk/data/plotinfo/ATLAS_2012_CONF_2012_033.plot Mon Apr 16 12:57:08 2012 (r3682)
@@ -1,241 +1,109 @@
-# BEGIN PLOT /ATLAS_2011_S9225137/count_7j55
-Title=Number of events in 7 jet, $p_\perp>55$\,GeV signal region
+# BEGIN PLOT /ATLAS_2012_CONF_2012_033/count_A_tight
+Title=Number of events in the tight A signal region
XLabel=
YLabel=Number of Events
-YLog=0
+LogY=0
# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/count_8j55
-Title=Number of events in 8 jet, $p_\perp>55$\,GeV signal region
+# BEGIN PLOT /ATLAS_2012_CONF_2012_033/count_A_medium
+Title=Number of events in the medium A signal region
XLabel=
YLabel=Number of Events
-YLog=0
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/count_6j80
-Title=Number of events in 6 jet, $p_\perp>80$\,GeV signal region
-XLabel=
-YLabel=Number of Events
-YLog=0
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/count_7j80
-Title=Number of events in 7 jet, $p_\perp>80$\,GeV signal region
-XLabel=
-YLabel=Number of Events
-YLog=0
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d01-x01-y01
-Title=Observed $E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ for 6 jets with $p_\perp>55$\,GeV
-XLabel=$E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ [$\mathrm{GeV}^{\frac12}$]
-YLabel=$\mathrm{Events}/0.25$\,GeV
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d01-x01-y02
-Title=Background $E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ for 6 jets with $p_\perp>55$\,GeV
-XLabel=$E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ [$\mathrm{GeV}^{\frac12}$]
-YLabel=$\mathrm{Events}/0.25$\,GeV
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d01-x01-y03
-Title=Signal $E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ for 6 jets with $p_\perp>55$\,GeV
-XLabel=$E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ [$\mathrm{GeV}^{\frac12}$]
-YLabel=$\mathrm{Events}/0.25$\,GeV
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d13-x01-y01
-Title=Observed $E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ for 7 jets with $p_\perp>55$\,GeV
-XLabel=$E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ [$\mathrm{GeV}^{\frac12}$]
-YLabel=$\mathrm{Events}/0.25$\,GeV
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d13-x01-y02
-Title=Background $E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ for 7 jets with $p_\perp>55$\,GeV
-XLabel=$E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ [$\mathrm{GeV}^{\frac12}$]
-YLabel=$\mathrm{Events}/0.25$\,GeV
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d13-x01-y03
-Title=Signal $E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ for 7 jets with $p_\perp>55$\,GeV
-XLabel=$E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ [$\mathrm{GeV}^{\frac12}$]
-YLabel=$\mathrm{Events}/0.25$\,GeV
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d15-x01-y01
-Title=Observed $E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ for 8 jets with $p_\perp>55$\,GeV
-XLabel=$E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ [$\mathrm{GeV}^{\frac12}$]
-YLabel=$\mathrm{Events}/0.25$\,GeV
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d15-x01-y02
-Title=Background $E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ for 8 jets with $p_\perp>55$\,GeV
-XLabel=$E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ [$\mathrm{GeV}^{\frac12}$]
-YLabel=$\mathrm{Events}/0.25$\,GeV
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d15-x01-y03
-Title=Signal $E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ for 8 jets with $p_\perp>55$\,GeV
-XLabel=$E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ [$\mathrm{GeV}^{\frac12}$]
-YLabel=$\mathrm{Events}/0.25$\,GeV
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d02-x01-y01
-Title=Observed $E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ for 5 jets with $p_\perp>80$\,GeV
-XLabel=$E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ [$\mathrm{GeV}^{\frac12}$]
-YLabel=$\mathrm{Events}/0.25$\,GeV
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d02-x01-y02
-Title=Background $E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ for 5 jets with $p_\perp>80$\,GeV
-XLabel=$E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ [$\mathrm{GeV}^{\frac12}$]
-YLabel=$\mathrm{Events}/0.25$\,GeV
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d02-x01-y03
-Title=Signal $E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ for 5 jets with $p_\perp>80$\,GeV
-XLabel=$E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ [$\mathrm{GeV}^{\frac12}$]
-YLabel=$\mathrm{Events}/0.25$\,GeV
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d14-x01-y01
-Title=Observed $E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ for 6 jets with $p_\perp>80$\,GeV
-XLabel=$E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ [$\mathrm{GeV}^{\frac12}$]
-YLabel=$\mathrm{Events}/0.25$\,GeV
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d14-x01-y02
-Title=Background $E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ for 6 jets with $p_\perp>80$\,GeV
-XLabel=$E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ [$\mathrm{GeV}^{\frac12}$]
-YLabel=$\mathrm{Events}/0.25$\,GeV
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d14-x01-y03
-Title=Signal $E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ for 6 jets with $p_\perp>80$\,GeV
-XLabel=$E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ [$\mathrm{GeV}^{\frac12}$]
-YLabel=$\mathrm{Events}/0.25$\,GeV
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d16-x01-y01
-Title=Observed $E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ for 7 jets with $p_\perp>80$\,GeV
-XLabel=$E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ [$\mathrm{GeV}^{\frac12}$]
-YLabel=$\mathrm{Events}/0.25$\,GeV
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d16-x01-y02
-Title=Background $E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ for 7 jets with $p_\perp>80$\,GeV
-XLabel=$E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ [$\mathrm{GeV}^{\frac12}$]
-YLabel=$\mathrm{Events}/0.25$\,GeV
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d16-x01-y03
-Title=Signal $E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ for 7 jets with $p_\perp>80$\,GeV
-XLabel=$E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ [$\mathrm{GeV}^{\frac12}$]
-YLabel=$\mathrm{Events}/0.25$\,GeV
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d03-x01-y01
-Title=Observed number of jets with $p_\perp>55$\,GeV for $1.5<E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}<2\,\mathrm{GeV}^{\frac12}$
-XLabel=Number of Jets
-YLabel=Number of events
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d03-x01-y02
-Title=Background number of jets with $p_\perp>55$\,GeV for $1.5<E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}<2\,\mathrm{GeV}^{\frac12}$
-XLabel=Number of Jets
-YLabel=Number of events
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d03-x01-y03
-Title=Signal number of jets with $p_\perp>55$\,GeV for $1.5<E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}<2\,\mathrm{GeV}^{\frac12}$
-XLabel=Number of Jets
-YLabel=Number of events
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d04-x01-y01
-Title=Observed number of jets with $p_\perp>55$\,GeV for $2<E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}<3\,\mathrm{GeV}^{\frac12}$
-XLabel=Number of Jets
-YLabel=Number of events
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d04-x01-y02
-Title=Background number of jets with $p_\perp>55$\,GeV for $2<E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}<3\,\mathrm{GeV}^{\frac12}$
-XLabel=Number of Jets
-YLabel=Number of events
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d04-x01-y03
-Title=Signal number of jets with $p_\perp>55$\,GeV for $2<E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}<3\,\mathrm{GeV}^{\frac12}$
-XLabel=Number of Jets
-YLabel=Number of events
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d17-x01-y01
-Title=Observed number of jets with $p_\perp>55$\,GeV for $3<E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}\,\mathrm{GeV}^{\frac12}$
-XLabel=Number of Jets
-YLabel=Number of events
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d17-x01-y02
-Title=Background number of jets with $p_\perp>55$\,GeV for $3<E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}\,\mathrm{GeV}^{\frac12}$
-XLabel=Number of Jets
-YLabel=Number of events
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d17-x01-y03
-Title=Signal number of jets with $p_\perp>55$\,GeV for $3<E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}\,\mathrm{GeV}^{\frac12}$
-XLabel=Number of Jets
-YLabel=Number of events
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d05-x01-y01
-Title=Observed number of jets with $p_\perp>80$\,GeV for $1.5<E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}<2\,\mathrm{GeV}^{\frac12}$
-XLabel=Number of Jets
-YLabel=Number of events
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d05-x01-y02
-Title=Background number of jets with $p_\perp>80$\,GeV for $1.5<E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}<2\,\mathrm{GeV}^{\frac12}$
-XLabel=Number of Jets
-YLabel=Number of events
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d05-x01-y03
-Title=Signal number of jets with $p_\perp>80$\,GeV for $1.5<E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}<2\,\mathrm{GeV}^{\frac12}$
-XLabel=Number of Jets
-YLabel=Number of events
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d06-x01-y01
-Title=Observed number of jets with $p_\perp>80$\,GeV for $2<E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}<3\,\mathrm{GeV}^{\frac12}$
-XLabel=Number of Jets
-YLabel=Number of events
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d06-x01-y02
-Title=Background number of jets with $p_\perp>80$\,GeV for $2<E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}<3\,\mathrm{GeV}^{\frac12}$
-XLabel=Number of Jets
-YLabel=Number of events
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d06-x01-y03
-Title=Signal number of jets with $p_\perp>80$\,GeV for $2<E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}<3\,\mathrm{GeV}^{\frac12}$
-XLabel=Number of Jets
-YLabel=Number of events
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d18-x01-y01
-Title=Observed number of jets with $p_\perp>80$\,GeV for $3<E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}\,\mathrm{GeV}^{\frac12}$
-XLabel=Number of Jets
-YLabel=Number of events
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d18-x01-y02
-Title=Background number of jets with $p_\perp>80$\,GeV for $3<E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}\,\mathrm{GeV}^{\frac12}$
-XLabel=Number of Jets
-YLabel=Number of events
-YLog=1
-# END PLOT
-# BEGIN PLOT /ATLAS_2011_S9225137/d18-x01-y03
-Title=Signal number of jets with $p_\perp>80$\,GeV for $3<E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}\,\mathrm{GeV}^{\frac12}$
-XLabel=Number of Jets
-YLabel=Number of events
-YLog=1
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_033/count_A_loose
+Title=Number of events in the loose A signal region
+XLabel=
+YLabel=Number of Events
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_033/count_Ap_medium
+Title=Number of events in the medium ${\rm A}'$ signal region
+XLabel=
+YLabel=Number of Events
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_033/count_B_tight
+Title=Number of events in the tight B signal region
+XLabel=
+YLabel=Number of Events
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_033/count_C_tight
+Title=Number of events in the tight C signal region
+XLabel=
+YLabel=Number of Events
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_033/count_C_medium
+Title=Number of events in the medium C signal region
+XLabel=
+YLabel=Number of Events
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_033/count_C_loose
+Title=Number of events in the loose C signal region
+XLabel=
+YLabel=Number of Events
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_033/count_D_tight
+Title=Number of events in the tight D signal region
+XLabel=
+YLabel=Number of Events
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_033/count_E_tight
+Title=Number of events in the tight E signal region
+XLabel=
+YLabel=Number of Events
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_033/count_E_medium
+Title=Number of events in the medium E signal region
+XLabel=
+YLabel=Number of Events
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_033/count_E_loose
+Title=Number of events in the loose E signal region
+XLabel=
+YLabel=Number of Events
+LogY=0
# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_033/hist_m_eff_A
+Title=Effective Mass in the A Signal Region
+XLabel=$M_{\mathrm{eff}}$ [GeV]
+YLabel=$\mathrm{Events}/100 \mathrm{GeV}$
+LogY=1
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_033/hist_m_eff_Ap
+Title=Effective Mass in the ${\rm A}'$ Signal Region
+XLabel=$M_{\mathrm{eff}}$ [GeV]
+YLabel=$\mathrm{Events}/100 \mathrm{GeV}$
+LogY=1
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_033/hist_m_eff_B
+Title=Effective Mass in the B Signal Region
+XLabel=$M_{\mathrm{eff}}$ [GeV]
+YLabel=$\mathrm{Events}/100 \mathrm{GeV}$
+LogY=1
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_033/hist_m_eff_C
+Title=Effective Mass in the C Signal Region
+XLabel=$M_{\mathrm{eff}}$ [GeV]
+YLabel=$\mathrm{Events}/100 \mathrm{GeV}$
+LogY=1
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_033/hist_m_eff_D
+Title=Effective Mass in the D Signal Region
+XLabel=$M_{\mathrm{eff}}$ [GeV]
+YLabel=$\mathrm{Events}/100 \mathrm{GeV}$
+LogY=1
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_033/hist_m_eff_E
+Title=Effective Mass in the E Signal Region
+XLabel=$M_{\mathrm{eff}}$ [GeV]
+YLabel=$\mathrm{Events}/100 \mathrm{GeV}$
+LogY=1
+# END PLOT
Modified: trunk/data/plotinfo/Makefile.am
==============================================================================
--- trunk/data/plotinfo/Makefile.am Sun Apr 15 11:38:39 2012 (r3681)
+++ trunk/data/plotinfo/Makefile.am Mon Apr 16 12:57:08 2012 (r3682)
@@ -41,6 +41,7 @@
ATLAS_2012_I1083318.plot \
ATLAS_2012_I1084540.plot \
ATLAS_2012_I1094568.plot \
+ ATLAS_2012_CONF_2012_033.plot \
BABAR_2005_S6181155.plot \
BABAR_2007_S6895344.plot \
BABAR_2007_S7266081.plot \
Copied and modified: trunk/src/Analyses/ATLAS_2012_CONF_2012_033.cc (from r3651, trunk/src/Analyses/ATLAS_2011_S9212183.cc)
==============================================================================
--- trunk/src/Analyses/ATLAS_2011_S9212183.cc Sun Mar 25 11:21:04 2012 (r3651, copy source)
+++ trunk/src/Analyses/ATLAS_2012_CONF_2012_033.cc Mon Apr 16 12:57:08 2012 (r3682)
@@ -8,20 +8,20 @@
#include "Rivet/Projections/VisibleFinalState.hh"
#include "Rivet/Projections/IdentifiedFinalState.hh"
#include "Rivet/Projections/FastJets.hh"
+#include "Rivet/Projections/VetoedFinalState.hh"
namespace Rivet {
-
- /// @author Chris Wymant
- class ATLAS_2011_S9212183 : public Analysis {
+ /// @author Peter Richardson
+ class ATLAS_2012_CONF_2012_033 : public Analysis {
public:
/// @name Constructors etc.
//@{
/// Constructor
- ATLAS_2011_S9212183()
- : Analysis("ATLAS_2011_S9212183")
+ ATLAS_2012_CONF_2012_033()
+ : Analysis("ATLAS_2012_CONF_2012_033")
{ }
//@}
@@ -50,7 +50,9 @@
addProjection(muons, "muons");
// Jet finder
- addProjection(FastJets(FinalState(), FastJets::ANTIKT, 0.4), "AntiKtJets04");
+ VetoedFinalState vfs;
+ vfs.addVetoPairId(MUON);
+ addProjection(FastJets(vfs, FastJets::ANTIKT, 0.4), "AntiKtJets04");
// All tracks (to do deltaR with leptons)
addProjection(ChargedFinalState(-3.0,3.0),"cfs");
@@ -58,20 +60,26 @@
// Used for pTmiss (N.B. the real 'vfs' extends beyond 4.5 to |eta| = 4.9)
addProjection(VisibleFinalState(-4.5,4.5),"vfs");
-
// Book histograms
- _count_2j = bookHistogram1D("count_2j", 1, 0., 1.);
- _count_3j = bookHistogram1D("count_3j", 1, 0., 1.);
- _count_4j5 = bookHistogram1D("count_4j5", 1, 0., 1.);
- _count_4j10 = bookHistogram1D("count_4j10", 1, 0., 1.);
- _count_HM = bookHistogram1D("count_HM", 1, 0., 1.);
-
- _hist_meff_2j = bookHistogram1D("m_eff_2j", 30, 0., 3000.);
- _hist_meff_3j = bookHistogram1D("m_eff_3j", 30, 0., 3000.);
- _hist_meff_4j = bookHistogram1D("m_eff_4j", 30, 0., 3000.);
- _hist_meff_HM = bookHistogram1D("m_eff_HM", 20, 0., 3000.);
+ _count_A_tight = bookHistogram1D("count_A_tight" , 1, 0., 1.);
+ _count_A_medium = bookHistogram1D("count_A_medium" , 1, 0., 1.);
+ _count_Ap_medium = bookHistogram1D("count_Ap_medium" , 1, 0., 1.);
+ _count_B_tight = bookHistogram1D("count_B_tight" , 1, 0., 1.);
+ _count_C_tight = bookHistogram1D("count_C_tight" , 1, 0., 1.);
+ _count_C_medium = bookHistogram1D("count_C_medium" , 1, 0., 1.);
+ _count_C_loose = bookHistogram1D("count_C_loose" , 1, 0., 1.);
+ _count_D_tight = bookHistogram1D("count_D_tight" , 1, 0., 1.);
+ _count_E_tight = bookHistogram1D("count_E_tight" , 1, 0., 1.);
+ _count_E_medium = bookHistogram1D("count_E_medium" , 1, 0., 1.);
+ _count_E_loose = bookHistogram1D("count_E_loose" , 1, 0., 1.);
+
+ _hist_meff_A = bookHistogram1D("hist_m_eff_A" , 30, 0., 3000.);
+ _hist_meff_Ap = bookHistogram1D("hist_m_eff_Ap", 30, 0., 3000.);
+ _hist_meff_B = bookHistogram1D("hist_m_eff_B" , 30, 0., 3000.);
+ _hist_meff_C = bookHistogram1D("hist_m_eff_C" , 30, 0., 3000.);
+ _hist_meff_D = bookHistogram1D("hist_m_eff_D" , 30, 0., 3000.);
+ _hist_meff_E = bookHistogram1D("hist_m_eff_E" , 30, 0., 3000.);
- _hist_eTmiss = bookHistogram1D("Et_miss", 20, 0., 1000.);
}
@@ -82,15 +90,14 @@
Jets cand_jets;
const Jets jets = applyProjection<FastJets>(event, "AntiKtJets04").jetsByPt(20.0*GeV);
foreach (const Jet& jet, jets) {
- if ( fabs( jet.momentum().eta() ) < 4.9 ) {
- cand_jets.push_back(jet);
- }
+ if ( fabs( jet.momentum().eta() ) < 4.9 ) {
+ cand_jets.push_back(jet);
+ }
}
const ParticleVector cand_e = applyProjection<IdentifiedFinalState>(event, "elecs").particlesByPt();
- // Muon isolation not mentioned in hep-exp 1109.6572, unlike in 1102.5290,
- // but assumed to still be applicable
+ // Muon isolation not mentioned in hep-exp 1109.6572 but assumed to still be applicable
ParticleVector cand_mu;
const ParticleVector chg_tracks = applyProjection<ChargedFinalState>(event, "cfs").particles();
const ParticleVector muons = applyProjection<IdentifiedFinalState>(event, "muons").particlesByPt();
@@ -105,28 +112,24 @@
}
// Resolve jet-lepton overlap for jets with |eta| < 2.8
- Jets cand_jets_2;
+ Jets recon_jets;
foreach ( const Jet& jet, cand_jets ) {
- if ( fabs( jet.momentum().eta() ) >= 2.8 ) {
- cand_jets_2.push_back( jet );
- } else {
- bool away_from_e = true;
- foreach ( const Particle & e, cand_e ) {
- if ( deltaR(e.momentum(),jet.momentum()) <= 0.2 ) {
- away_from_e = false;
- break;
- }
- }
- if ( away_from_e ) cand_jets_2.push_back( jet );
- }
+ if ( fabs( jet.momentum().eta() ) >= 2.8 ) continue;
+ bool away_from_e = true;
+ foreach ( const Particle & e, cand_e ) {
+ if ( deltaR(e.momentum(),jet.momentum()) <= 0.2 ) {
+ away_from_e = false;
+ break;
+ }
+ }
+ if ( away_from_e ) recon_jets.push_back( jet );
}
-
ParticleVector recon_e, recon_mu;
foreach ( const Particle & e, cand_e ) {
bool away = true;
- foreach ( const Jet& jet, cand_jets_2 ) {
+ foreach ( const Jet& jet, recon_jets ) {
if ( deltaR(e.momentum(),jet.momentum()) < 0.4 ) {
away = false;
break;
@@ -137,7 +140,7 @@
foreach ( const Particle & mu, cand_mu ) {
bool away = true;
- foreach ( const Jet& jet, cand_jets_2 ) {
+ foreach ( const Jet& jet, recon_jets ) {
if ( deltaR(mu.momentum(),jet.momentum()) < 0.4 ) {
away = false;
break;
@@ -146,7 +149,6 @@
if ( away ) recon_mu.push_back( mu );
}
-
// pTmiss
// Based on all candidate electrons, muons and jets, plus everything else with |eta| < 4.5
// i.e. everything in our projection "vfs" plus the jets with |eta| > 4.5
@@ -155,38 +157,25 @@
foreach ( const Particle & p, vfs_particles ) {
pTmiss -= p.momentum();
}
- foreach ( const Jet& jet, cand_jets_2 ) {
+ foreach ( const Jet& jet, cand_jets ) {
if ( fabs( jet.momentum().eta() ) > 4.5 ) pTmiss -= jet.momentum();
}
double eTmiss = pTmiss.pT();
-
- // Final jet filter
- Jets recon_jets;
- foreach ( const Jet& jet, cand_jets_2 ) {
- if ( fabs( jet.momentum().eta() ) <= 2.8 ) recon_jets.push_back( jet );
- }
- // NB. It seems that jets with |eta| > 2.8 could have been thrown away at
- // the start; we don't do so, in order to follow both the structure of
- // the paper and the similar Rivet analysis ATLAS_2011_S8983313
-
- // 'candidate' muons needed only 10 GeV, to cause a veto they need 20 GeV
- ParticleVector veto_mu;
- foreach ( const Particle & mu, cand_mu ) {
- if ( mu.momentum().pT() >= 20.0*GeV ) veto_mu.push_back(mu);
- }
-
- if ( ! ( veto_mu.empty() && recon_e.empty() ) ) {
+ // no electron pT> 20 or muons pT>10
+ if ( !recon_mu.empty() || !recon_e.empty() ) {
MSG_DEBUG("Charged leptons left after selection");
vetoEvent;
}
-
- if ( eTmiss <= 130 * GeV ) {
+
+ if ( eTmiss <= 160 * GeV ) {
MSG_DEBUG("Not enough eTmiss: " << eTmiss << " < 130");
vetoEvent;
}
- if ( recon_jets.empty() || recon_jets[0].momentum().pT() <= 130.0 * GeV ) {
+ if ( recon_jets.size()<2 ||
+ recon_jets[0].momentum().pT() <= 130.0 * GeV ||
+ recon_jets[0].momentum().pT() <= 60.0 * GeV ) {
MSG_DEBUG("No hard leading jet in " << recon_jets.size() << " jets");
vetoEvent;
}
@@ -194,120 +183,115 @@
// ==================== observables ====================
int Njets = 0;
- double min_dPhi = 999.999;
+ double min_dPhi_All = 999.999;
+ double min_dPhi_2 = 999.999;
+ double min_dPhi_3 = 999.999;
double pTmiss_phi = pTmiss.phi();
foreach ( const Jet& jet, recon_jets ) {
- if ( jet.momentum().pT() > 40 * GeV ) {
- if ( Njets < 3 ) {
- min_dPhi = min( min_dPhi, deltaPhi( pTmiss_phi, jet.momentum().phi() ) );
- }
- ++Njets;
- }
+ if ( jet.momentum().pT() < 40 * GeV ) continue;
+ if ( Njets < 2 ) {
+ min_dPhi_2 = min( min_dPhi_2, deltaPhi( pTmiss_phi, jet.momentum().phi() ) );
+ }
+ if( Njets < 3) {
+ min_dPhi_3 = min( min_dPhi_3, deltaPhi( pTmiss_phi, jet.momentum().phi() ) );
+ }
+ min_dPhi_All = min( min_dPhi_All, deltaPhi( pTmiss_phi, jet.momentum().phi() ) );
+ ++Njets;
}
- int NjetsHighMass = 0;
+ // inclusive meff
+ double m_eff_inc = eTmiss;
foreach ( const Jet& jet, recon_jets ) {
- if ( jet.momentum().pT() > 80.0 * GeV ) {
- ++NjetsHighMass;
- }
+ double perp = jet.momentum().pT();
+ if(perp>40.) m_eff_inc += perp;
}
- if ( Njets < 2 ) {
- MSG_DEBUG("Only " << Njets << " >40 GeV jets left");
- vetoEvent;
- }
-
- if ( min_dPhi <= 0.4 ) {
- MSG_DEBUG("dPhi too small");
- vetoEvent;
- }
-
- // m_eff
- double m_eff_2j = eTmiss + recon_jets[0].momentum().pT() + recon_jets[1].momentum().pT();
- double m_eff_3j = recon_jets.size() < 3 ? -999.0 : m_eff_2j + recon_jets[2].momentum().pT();
- double m_eff_4j = recon_jets.size() < 4 ? -999.0 : m_eff_3j + recon_jets[3].momentum().pT();
- double m_eff_HM = eTmiss;
- foreach ( const Jet& jet, recon_jets ) {
- if ( jet.momentum().pT() > 40.0 * GeV ) m_eff_HM += jet.momentum().pT();
+ // region A
+ double m_eff_Nj = eTmiss + recon_jets[0].momentum().pT() + recon_jets[1].momentum().pT();
+ if( min_dPhi_2 > 0.4 && eTmiss/m_eff_Nj > 0.3 ) {
+ _hist_meff_A->fill(m_eff_inc,weight);
+ if(m_eff_inc>1900.) _count_A_tight ->fill(0.5,weight);
+ if(m_eff_inc>1400.) _count_A_medium->fill(0.5,weight);
+ }
+
+ // region A'
+ if( min_dPhi_2 > 0.4 && eTmiss/m_eff_Nj > 0.4 ) {
+ _hist_meff_Ap->fill(m_eff_inc,weight);
+ if(m_eff_inc>1200.) _count_Ap_medium->fill(0.5,weight);
+ }
+
+ // for rest of regions 3 jets pT> 60 needed
+ if(recon_jets.size()<3 || recon_jets[2].momentum().perp()<60.)
+ vetoEvent;
+
+ // region B
+ m_eff_Nj += recon_jets[2].momentum().perp();
+ if( min_dPhi_3 > 0.4 && eTmiss/m_eff_Nj > 0.25 ) {
+ _hist_meff_B->fill(m_eff_inc,weight);
+ if(m_eff_inc>1900.) _count_B_tight ->fill(0.5,weight);
+ }
+
+ // for rest of regions 4 jets pT> 60 needed
+ if(recon_jets.size()<4 || recon_jets[3].momentum().perp()<60.)
+ vetoEvent;
+
+ // region C
+ m_eff_Nj += recon_jets[3].momentum().perp();
+ if( min_dPhi_3 > 0.4 && min_dPhi_All > 0.2 && eTmiss/m_eff_Nj > 0.25 ) {
+ _hist_meff_C->fill(m_eff_inc,weight);
+ if(m_eff_inc>1500.) _count_C_tight ->fill(0.5,weight);
+ if(m_eff_inc>1200.) _count_C_medium->fill(0.5,weight);
+ if(m_eff_inc> 900.) _count_C_loose ->fill(0.5,weight);
+ }
+
+ // for rest of regions 5 jets pT> 40 needed
+ if(recon_jets.size()<5 || recon_jets[4].momentum().perp()<40.)
+ vetoEvent;
+
+ // region D
+ m_eff_Nj += recon_jets[4].momentum().perp();
+ if( min_dPhi_3 > 0.4 && min_dPhi_All > 0.2 && eTmiss/m_eff_Nj > 0.2 ) {
+ _hist_meff_D->fill(m_eff_inc,weight);
+ if(m_eff_inc>1500.) _count_D_tight ->fill(0.5,weight);
+ }
+
+ // for rest of regions 6 jets pT> 40 needed
+ if(recon_jets.size()<6 || recon_jets[5].momentum().perp()<40.)
+ vetoEvent;
+
+ // region E
+ m_eff_Nj += recon_jets[5].momentum().perp();
+ if( min_dPhi_3 > 0.4 && min_dPhi_All > 0.2 && eTmiss/m_eff_Nj > 0.15 ) {
+ _hist_meff_E->fill(m_eff_inc,weight);
+ if(m_eff_inc>1400.) _count_E_tight ->fill(0.5,weight);
+ if(m_eff_inc>1200.) _count_E_medium->fill(0.5,weight);
+ if(m_eff_inc> 900.) _count_E_loose ->fill(0.5,weight);
}
-
- double et_meff_2j = eTmiss / m_eff_2j;
- double et_meff_3j = eTmiss / m_eff_3j;
- double et_meff_4j = eTmiss / m_eff_4j;
- double et_meff_HM = eTmiss / m_eff_HM;
-
-
- // ==================== FILL ====================
-
- MSG_DEBUG( "Trying to fill "
- << Njets << ' '
- << m_eff_2j << ' '
- << et_meff_2j << ' '
- << m_eff_3j << ' '
- << et_meff_3j << ' '
- << m_eff_4j << ' '
- << et_meff_4j << ' '
- << m_eff_HM << ' '
- << et_meff_HM );
-
-
- _hist_eTmiss->fill(eTmiss, weight);
-
-
- // 2j region
- if ( et_meff_2j > 0.3 ) {
- _hist_meff_2j->fill(m_eff_2j, weight);
- if ( m_eff_2j > 1000 * GeV ) {
- MSG_DEBUG("Hits 2j");
- _count_2j->fill(0.5, weight);
- }
- }
-
-
- // 3j region
- if ( Njets >= 3 && et_meff_3j > 0.25 ) {
- _hist_meff_3j->fill(m_eff_3j, weight);
- if ( m_eff_3j > 1000 * GeV ) {
- MSG_DEBUG("Hits 3j");
- _count_3j->fill(0.5, weight);
- }
- }
-
-
- // 4j5 & 4j10 regions
- if ( Njets >= 4 && et_meff_4j > 0.25 ) {
- _hist_meff_4j->fill(m_eff_4j, weight);
- if ( m_eff_4j > 500 * GeV ) {
- MSG_DEBUG("Hits 4j5");
- _count_4j5->fill(0.5, weight);
- }
- if ( m_eff_4j > 1000 * GeV ) {
- MSG_DEBUG("Hits 4j10");
- _count_4j10->fill(0.5, weight);
- }
- }
-
-
- // High mass region
- if ( NjetsHighMass >= 4 && et_meff_HM > 0.2 ) {
- _hist_meff_HM->fill(m_eff_HM, weight);
- if ( m_eff_HM > 1100 * GeV ) {
- MSG_DEBUG("Hits HM");
- _count_HM->fill(0.5, weight);
- }
- }
-
}
void finalize() {
- // Two, three and four jet channels have bin width = 100 (GeV)
- // High mass channel has bin width = 150 (GeV)
- // Integrated luminosity = 1040 (pb)
- scale( _hist_meff_2j, 100. * 1040 * crossSection()/sumOfWeights() );
- scale( _hist_meff_3j, 100. * 1040 * crossSection()/sumOfWeights() );
- scale( _hist_meff_4j, 100. * 1040 * crossSection()/sumOfWeights() );
- scale( _hist_meff_HM, 150. * 1040 * crossSection()/sumOfWeights() );
+
+ double norm = crossSection()/femtobarn*4.7/sumOfWeights();
+ // these are number of events at 4.7fb^-1 per 100 GeV
+ scale( _hist_meff_A , 100. * norm );
+ scale( _hist_meff_Ap, 100. * norm );
+ scale( _hist_meff_B , 100. * norm );
+ scale( _hist_meff_C , 100. * norm );
+ scale( _hist_meff_D , 100. * norm );
+ scale( _hist_meff_E , 100. * norm );
+ // these are number of events at 4.7fb^-1
+ scale(_count_A_tight ,norm);
+ scale(_count_A_medium ,norm);
+ scale(_count_Ap_medium,norm);
+ scale(_count_B_tight ,norm);
+ scale(_count_C_tight ,norm);
+ scale(_count_C_medium ,norm);
+ scale(_count_C_loose ,norm);
+ scale(_count_D_tight ,norm);
+ scale(_count_E_tight ,norm);
+ scale(_count_E_medium ,norm);
+ scale(_count_E_loose ,norm);
}
//@}
@@ -315,21 +299,29 @@
private:
- AIDA::IHistogram1D* _count_2j;
- AIDA::IHistogram1D* _count_3j;
- AIDA::IHistogram1D* _count_4j5;
- AIDA::IHistogram1D* _count_4j10;
- AIDA::IHistogram1D* _count_HM;
- AIDA::IHistogram1D* _hist_meff_2j;
- AIDA::IHistogram1D* _hist_meff_3j;
- AIDA::IHistogram1D* _hist_meff_4j;
- AIDA::IHistogram1D* _hist_meff_HM;
- AIDA::IHistogram1D* _hist_eTmiss;
-
+ AIDA::IHistogram1D* _count_A_tight;
+ AIDA::IHistogram1D* _count_A_medium;
+ AIDA::IHistogram1D* _count_Ap_medium;
+ AIDA::IHistogram1D* _count_B_tight;
+ AIDA::IHistogram1D* _count_C_tight;
+ AIDA::IHistogram1D* _count_C_medium;
+ AIDA::IHistogram1D* _count_C_loose;
+ AIDA::IHistogram1D* _count_D_tight;
+ AIDA::IHistogram1D* _count_E_tight;
+ AIDA::IHistogram1D* _count_E_medium;
+ AIDA::IHistogram1D* _count_E_loose;
+
+ AIDA::IHistogram1D* _hist_meff_A ;
+ AIDA::IHistogram1D* _hist_meff_Ap;
+ AIDA::IHistogram1D* _hist_meff_B ;
+ AIDA::IHistogram1D* _hist_meff_C ;
+ AIDA::IHistogram1D* _hist_meff_D ;
+ AIDA::IHistogram1D* _hist_meff_E ;
+
};
// This global object acts as a hook for the plugin system
- DECLARE_RIVET_PLUGIN(ATLAS_2011_S9212183);
+ DECLARE_RIVET_PLUGIN(ATLAS_2012_CONF_2012_033);
}
Modified: trunk/src/Analyses/Makefile.am
==============================================================================
--- trunk/src/Analyses/Makefile.am Sun Apr 15 11:38:39 2012 (r3681)
+++ trunk/src/Analyses/Makefile.am Mon Apr 16 12:57:08 2012 (r3682)
@@ -76,7 +76,8 @@
ATLAS_2011_CONF_2011_098.cc \
ATLAS_2011_S9108483.cc \
ATLAS_2011_S9225137.cc \
- ATLAS_2012_I1083318.cc
+ ATLAS_2012_I1083318.cc \
+ ATLAS_2012_CONF_2012_033.cc
endif
More information about the Rivet-svn
mailing list
|