|
|
[Rivet-svn] r3686 - in trunk: data/anainfo data/plotinfo src/Analyses
blackhole at projects.hepforge.org
blackhole at projects.hepforge.org
Mon Apr 16 13:16:53 BST 2012
Author: richardn
Date: Mon Apr 16 13:16:53 2012
New Revision: 3686
Log:
added ATLAS_2012_CONF_2012_023, latest 3 lepton + MET SUSY search
Added:
trunk/data/anainfo/ATLAS_2012_CONF_2012_023.info
- copied, changed from r3651, trunk/data/anainfo/ATLAS_2011_S9225137.info
trunk/data/plotinfo/ATLAS_2012_CONF_2012_023.plot
- copied, changed from r3651, trunk/data/plotinfo/ATLAS_2011_S9225137.plot
trunk/src/Analyses/ATLAS_2012_CONF_2012_023.cc
- copied, changed from r3647, trunk/src/Analyses/ATLAS_2011_S9225137.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_023.info (from r3651, trunk/data/anainfo/ATLAS_2011_S9225137.info)
==============================================================================
--- trunk/data/anainfo/ATLAS_2011_S9225137.info Sun Mar 25 11:21:04 2012 (r3651, copy source)
+++ trunk/data/anainfo/ATLAS_2012_CONF_2012_023.info Mon Apr 16 13:16:53 2012 (r3686)
@@ -1,14 +1,13 @@
-Name: ATLAS_2011_S9225137
+Name: ATLAS_2012_CONF_2012_023
Year: 2011
-Summary: High jet multiplicity squark and gluino search
+Summary: 3 lepton plus missing transverse energy SUSY search
Experiment: ATLAS
Collider: LHC
-SpiresID: 9225137
-Status: UNVALIDATED
+Status: PRELIMINARY
Authors:
- Peter Richardson <peter.richardson at durham.ac.uk>
References:
- - arXiv:1110.2299
+ - ATLAS-CONF-2012-023
RunInfo:
BSM signal events at 7000 GeV.
NumEvents: 25000 for BSM signals
@@ -16,28 +15,11 @@
Energies: [7000]
NeedCrossSection: yes
Description:
- 'Search for SUSY using events with 6 or more jets in
- association with missing transverse momentum produced in proton-proton
+ 'Search for SUSY using events with 3 leptons in association with
+ missing transverse energy in proton-proton
collisions at a centre-of-mass energy of 7 TeV. The data sample has a total
- integrated luminosity of 1.34 fb$^{-1}$. Distributions
- in the W and top control regions are not produced, while in addition to
- the plots from the paper the count of events in the different signal regions
- is included.'
-BibKey: Aad:2011qa
-BibTeX: '@article{Aad:2011qa,
- author = "Aad, Georges and others",
- title = "{Search for new phenomena in final states with large jet
- multiplicities and missing transverse momentum using
- sqrt(s)=7 TeV pp collisions with the ATLAS detector.}",
- collaboration = "Atlas Collaboration",
- journal = "JHEP",
- volume = "1111",
- pages = "099",
- doi = "10.1007/JHEP11(2011)099",
- year = "2011",
- eprint = "1110.2299",
- archivePrefix = "arXiv",
- primaryClass = "hep-ex",
- reportNumber = "CERN-PH-EP-2011-155",
- SLACcitation = "%%CITATION = ARXIV:1110.2299;%%",
-}'
+ integrated luminosity of 2.06 fb$^{-1}$. There is no reference data and
+ in addition to the control plots from the paper the number of events in
+ the two signal regions,
+ correctly normalized to an integrated luminosity 2.06 fb$^{-1}$,
+ are calculated.'
Modified: trunk/data/anainfo/Makefile.am
==============================================================================
--- trunk/data/anainfo/Makefile.am Mon Apr 16 13:14:25 2012 (r3685)
+++ trunk/data/anainfo/Makefile.am Mon Apr 16 13:16:53 2012 (r3686)
@@ -46,6 +46,7 @@
ATLAS_2012_I1094568.info \
ATLAS_2012_I1095236.info \
ATLAS_2012_CONF_2012_001.info \
+ ATLAS_2012_CONF_2012_023.info \
ATLAS_2012_CONF_2012_033.info \
ATLAS_2012_CONF_2012_037.info \
BABAR_2005_S6181155.info \
Copied and modified: trunk/data/plotinfo/ATLAS_2012_CONF_2012_023.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_023.plot Mon Apr 16 13:16:53 2012 (r3686)
@@ -1,241 +1,82 @@
-# 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_023/count_SR1
+Title=Number of events in the SR1 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_023/count_SR2
+Title=Number of events in the SR2 signal region
XLabel=
YLabel=Number of Events
-YLog=0
+LogY=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
+
+# BEGIN PLOT /ATLAS_2012_CONF_2012_023/hist_lepton_pT_1_SR1
+Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the leading lepton in SR1
+XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV]
+YLabel=$\mathrm{Events}/20 \mathrm{GeV}$
+LogY=1
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_023/hist_lepton_pT_2_SR1
+Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the $2^{\rm nd}$ leading lepton in SR1
+XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV]
+YLabel=$\mathrm{Events}/20 \mathrm{GeV}$
+LogY=1
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_023/hist_lepton_pT_3_SR1
+Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the $3^{\rm rd}$ leading lepton in SR1
+XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV]
+YLabel=$\mathrm{Events}/20 \mathrm{GeV}$
+LogY=1
# 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
+
+# BEGIN PLOT /ATLAS_2012_CONF_2012_023/hist_lepton_pT_1_SR2
+Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the leading lepton in SR2
+XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV]
+YLabel=$\mathrm{Events}/20 \mathrm{GeV}$
+LogY=1
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_023/hist_lepton_pT_2_SR2
+Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the $2^{\rm nd}$ leading lepton in SR2
+XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV]
+YLabel=$\mathrm{Events}/20 \mathrm{GeV}$
+LogY=1
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_023/hist_lepton_pT_3_SR2
+Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the $3^{\rm rd}$ leading lepton in SR2
+XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV]
+YLabel=$\mathrm{Events}/20 \mathrm{GeV}$
+LogY=1
# 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
+
+# BEGIN PLOT /ATLAS_2012_CONF_2012_023/hist_etmiss_SR1_A
+Title=Missing Transverse Energy in SR1
+XLabel=$E^{\rm miss}_{\mathrm{T}}$ [GeV]
+YLabel=$\mathrm{Events}/20$\,GeV
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_023/hist_etmiss_SR1_B
+Title=Missing Transverse Energy in SR1
+XLabel=$E^{\rm miss}_{\mathrm{T}}$ [GeV]
+YLabel=$\mathrm{Events}/20$\,GeV
+LogY=1
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_023/hist_etmiss_SR2_A
+Title=Missing Transverse Energy in SR2
+XLabel=$E^{\rm miss}_{\mathrm{T}}$ [GeV]
+YLabel=$\mathrm{Events}/20$\,GeV
+LogY=0
+# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_023/hist_etmiss_SR2_B
+Title=Missing Transverse Energy in SR2
+XLabel=$E^{\rm miss}_{\mathrm{T}}$ [GeV]
+YLabel=$\mathrm{Events}/20$\,GeV
+LogY=1
# END PLOT
+# BEGIN PLOT /ATLAS_2012_CONF_2012_023/hist_mSFOF
+Title=Mass of SFOF lepton pair closest to the $Z^0$ mass in SR1
+XLabel=$M_{\mathrm{SFOS}}$ [GeV]
+YLabel=$\mathrm{Events}/20$\,GeV
+LogY=1
+# END PLOT
Modified: trunk/data/plotinfo/Makefile.am
==============================================================================
--- trunk/data/plotinfo/Makefile.am Mon Apr 16 13:14:25 2012 (r3685)
+++ trunk/data/plotinfo/Makefile.am Mon Apr 16 13:16:53 2012 (r3686)
@@ -43,6 +43,7 @@
ATLAS_2012_I1094568.plot \
ATLAS_2012_I1095236.plot \
ATLAS_2012_CONF_2012_001.plot \
+ ATLAS_2012_CONF_2012_023.plot \
ATLAS_2012_CONF_2012_033.plot \
ATLAS_2012_CONF_2012_037.plot \
BABAR_2005_S6181155.plot \
Copied and modified: trunk/src/Analyses/ATLAS_2012_CONF_2012_023.cc (from r3647, trunk/src/Analyses/ATLAS_2011_S9225137.cc)
==============================================================================
--- trunk/src/Analyses/ATLAS_2011_S9225137.cc Thu Mar 22 20:35:24 2012 (r3647, copy source)
+++ trunk/src/Analyses/ATLAS_2012_CONF_2012_023.cc Mon Apr 16 13:16:53 2012 (r3686)
@@ -13,16 +13,16 @@
namespace Rivet {
-
- class ATLAS_2011_S9225137 : public Analysis {
+ /// @author Peter Richardson
+ class ATLAS_2012_CONF_2012_023 : public Analysis {
public:
/// @name Constructors etc.
//@{
/// Constructor
- ATLAS_2011_S9225137()
- : Analysis("ATLAS_2011_S9225137")
+ ATLAS_2012_CONF_2012_023()
+ : Analysis("ATLAS_2012_CONF_2012_023")
{ }
//@}
@@ -36,18 +36,10 @@
/// Book histograms and initialise projections before the run
void init() {
- // veto region electrons
- std::vector<std::pair<double, double> > eta_v_e;
- eta_v_e.push_back(make_pair(-1.52,-1.37));
- eta_v_e.push_back(make_pair( 1.37, 1.52));
- IdentifiedFinalState veto_elecs(eta_v_e, 10.0*GeV);
- veto_elecs.acceptIdPair(ELECTRON);
- addProjection(veto_elecs, "veto_elecs");
-
// projection to find the electrons
std::vector<std::pair<double, double> > eta_e;
eta_e.push_back(make_pair(-2.47,2.47));
- IdentifiedFinalState elecs(eta_e, 20.0*GeV);
+ IdentifiedFinalState elecs(eta_e, 10.0*GeV);
elecs.acceptIdPair(ELECTRON);
addProjection(elecs, "elecs");
@@ -71,69 +63,32 @@
// all tracks (to do deltaR with leptons)
addProjection(ChargedFinalState(-3.0,3.0),"cfs");
- /// Book histograms
- _etmisspT_55_NJ_6_obs = bookHistogram1D( 1,1,1);
- _etmisspT_55_NJ_6_bac = bookHistogram1D( 1,1,2);
- _etmisspT_55_NJ_6_sig = bookHistogram1D( 1,1,3);
- _etmisspT_55_NJ_7_obs = bookHistogram1D(13,1,1);
- _etmisspT_55_NJ_7_bac = bookHistogram1D(13,1,2);
- _etmisspT_55_NJ_7_sig = bookHistogram1D(13,1,3);
- _etmisspT_55_NJ_8_obs = bookHistogram1D(15,1,1);
- _etmisspT_55_NJ_8_bac = bookHistogram1D(15,1,2);
- _etmisspT_55_NJ_8_sig = bookHistogram1D(15,1,3);
- _etmisspT_80_NJ_5_obs = bookHistogram1D( 2,1,1);
- _etmisspT_80_NJ_5_bac = bookHistogram1D( 2,1,2);
- _etmisspT_80_NJ_5_sig = bookHistogram1D( 2,1,3);
- _etmisspT_80_NJ_6_obs = bookHistogram1D(14,1,1);
- _etmisspT_80_NJ_6_bac = bookHistogram1D(14,1,2);
- _etmisspT_80_NJ_6_sig = bookHistogram1D(14,1,3);
- _etmisspT_80_NJ_7_obs = bookHistogram1D(16,1,1);
- _etmisspT_80_NJ_7_bac = bookHistogram1D(16,1,2);
- _etmisspT_80_NJ_7_sig = bookHistogram1D(16,1,3);
-
- _njet55A_obs = bookHistogram1D( 3,1,1);
- _njet55A_bac = bookHistogram1D( 3,1,2);
- _njet55A_sig = bookHistogram1D( 3,1,3);
- _njet55B_obs = bookHistogram1D( 4,1,1);
- _njet55B_bac = bookHistogram1D( 4,1,2);
- _njet55B_sig = bookHistogram1D( 4,1,3);
- _njet55C_obs = bookHistogram1D(17,1,1);
- _njet55C_bac = bookHistogram1D(17,1,2);
- _njet55C_sig = bookHistogram1D(17,1,3);
- _njet80A_obs = bookHistogram1D( 5,1,1);
- _njet80A_bac = bookHistogram1D( 5,1,2);
- _njet80A_sig = bookHistogram1D( 5,1,3);
- _njet80B_obs = bookHistogram1D( 6,1,1);
- _njet80B_bac = bookHistogram1D( 6,1,2);
- _njet80B_sig = bookHistogram1D( 6,1,3);
- _njet80C_obs = bookHistogram1D(18,1,1);
- _njet80C_bac = bookHistogram1D(18,1,2);
- _njet80C_sig = bookHistogram1D(18,1,3);
-
- _count_7j55 = bookHistogram1D("count_7j55", 1, 0., 1.);
- _count_8j55 = bookHistogram1D("count_8j55", 1, 0., 1.);
- _count_6j80 = bookHistogram1D("count_6j80", 1, 0., 1.);
- _count_7j80 = bookHistogram1D("count_7j80", 1, 0., 1.);
+ // Book histograms
+ _hist_leptonpT_SR1.push_back(bookHistogram1D("hist_lepton_pT_1_SR1",11,0.,220.));
+ _hist_leptonpT_SR1.push_back(bookHistogram1D("hist_lepton_pT_2_SR1", 7,0.,140.));
+ _hist_leptonpT_SR1.push_back(bookHistogram1D("hist_lepton_pT_3_SR1", 8,0.,160.));
+ _hist_leptonpT_SR2.push_back(bookHistogram1D("hist_lepton_pT_1_SR2",11,0.,220.));
+ _hist_leptonpT_SR2.push_back(bookHistogram1D("hist_lepton_pT_2_SR2", 7,0.,140.));
+ _hist_leptonpT_SR2.push_back(bookHistogram1D("hist_lepton_pT_3_SR2", 8,0.,160.));
+ _hist_etmiss_SR1_A = bookHistogram1D("hist_etmiss_SR1_A",15,10.,310.);
+ _hist_etmiss_SR1_B = bookHistogram1D("hist_etmiss_SR1_B", 9,10.,190.);
+ _hist_etmiss_SR2_A = bookHistogram1D("hist_etmiss_SR2_A",15,10.,310.);
+ _hist_etmiss_SR2_B = bookHistogram1D("hist_etmiss_SR2_B", 9,10.,190.);
+ _hist_mSFOS= bookHistogram1D("hist_mSFOF",9,0.,180.);
+ _count_SR1 = bookHistogram1D("count_SR1", 1, 0., 1.);
+ _count_SR2 = bookHistogram1D("count_SR2", 1, 0., 1.);
}
/// Perform the per-event analysis
void analyze(const Event& event) {
const double weight = event.weight();
- // apply electron veto region
- ParticleVector veto_e
- = applyProjection<IdentifiedFinalState>(event, "veto_elecs").particles();
- if ( ! veto_e.empty() ) {
- MSG_DEBUG("electrons in veto region");
- vetoEvent;
- }
-
// get the jet candidates
Jets cand_jets;
foreach (const Jet& jet,
applyProjection<FastJets>(event, "AntiKtJets04").jetsByPt(20.0*GeV) ) {
- if ( fabs( jet.momentum().eta() ) < 4.9 ) {
+ if ( fabs( jet.momentum().eta() ) < 2.8 ) {
cand_jets.push_back(jet);
}
}
@@ -154,54 +109,82 @@
}
// candidate electrons
-
- ParticleVector cand_e =
- applyProjection<IdentifiedFinalState>(event, "elecs").particlesByPt();
+ ParticleVector cand_e;
+ foreach ( const Particle & e,
+ applyProjection<IdentifiedFinalState>(event, "elecs").particlesByPt() ) {
+ double eta = e.momentum().eta();
+ // remove electrons with pT<15 in old veto region
+ // (NOT EXPLICIT IN THIS PAPER BUT IN SIMILAR 4 LEPTON PAPER and THIS DESCRPITION
+ // IS MUCH WORSE SO ASSUME THIS IS DONE)
+ if( fabs(eta)>1.37 && fabs(eta) < 1.52 && e.momentum().perp()< 15.*GeV)
+ continue;
+ double pTinCone = -e.momentum().perp();
+ foreach ( const Particle & track, chg_tracks ) {
+ if ( deltaR(e.momentum(),track.momentum()) <= 0.2 )
+ pTinCone += track.momentum().pT();
+ }
+ if (pTinCone/e.momentum().perp()<0.1) {
+ cand_e.push_back(e);
+ }
+ }
// resolve jet/lepton ambiguity
- Jets cand_jets_2;
+ // (NOT EXPLICIT IN THIS PAPER BUT IN SIMILAR 4 LEPTON PAPER and THIS DESCRPITION
+ // IS MUCH WORSE SO ASSUME THIS IS DONE)
+ Jets recon_jets;
foreach ( const Jet& jet, cand_jets ) {
- // candidates above eta=2.8 are jets
- if ( fabs( jet.momentum().eta() ) >= 2.8 )
- cand_jets_2.push_back( jet );
- // otherwise more the R=0.2 from an electrons
- 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 );
- }
+ 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 );
}
// only keep electrons more than R=0.4 from jets
ParticleVector recon_e;
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;
}
}
+ // and 0.1 from any muons
+ if ( ! away ) {
+ foreach ( const Particle & mu, cand_e ) {
+ if ( deltaR(mu.momentum(),e.momentum()) < 0.1 ) {
+ away = false;
+ break;
+ }
+ }
+ }
if ( away )
recon_e.push_back( e );
}
-
// only keep muons more than R=0.4 from jets
ParticleVector recon_mu;
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;
}
}
+ // and 0.1 from any electrona
+ if ( ! away ) {
+ foreach ( const Particle & e, cand_e ) {
+ if ( deltaR(mu.momentum(),e.momentum()) < 0.1 ) {
+ away = false;
+ break;
+ }
+ }
+ }
if ( away )
recon_mu.push_back( mu );
}
@@ -215,245 +198,150 @@
}
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 );
- }
-
// now only use recon_jets, recon_mu, recon_e
- // reject events with electrons and muons
- if ( ! ( recon_mu.empty() && recon_e.empty() ) ) {
- MSG_DEBUG("Charged leptons left after selection");
+ // reject events with wrong number of leptons
+ if ( recon_mu.size() + recon_e.size() != 3 ) {
+ MSG_DEBUG("To few charged leptons left after selection");
vetoEvent;
}
- // calculate H_T
- double HT=0;
- foreach ( const Jet& jet, recon_jets ) {
- if ( jet.momentum().pT() > 40 * GeV )
- HT += jet.momentum().pT() ;
- }
-
- // number of jets and deltaR
- bool pass55DeltaR=true;
- unsigned int njet55=0;
- bool pass80DeltaR=true;
- unsigned int njet80=0;
- for (unsigned int ix=0;ix<recon_jets.size();++ix) {
- if(recon_jets[ix].momentum().pT()>80.*GeV) ++njet80;
- if(recon_jets[ix].momentum().pT()>55.*GeV) ++njet55;
-
- for (unsigned int iy=ix+1;iy<recon_jets.size();++iy) {
- if(recon_jets[ix].momentum().pT()>55.*GeV &&
- recon_jets[iy].momentum().pT()>55.*GeV &&
- deltaR(recon_jets[ix],recon_jets[iy]) <0.6 )
- pass55DeltaR = false;
- if(recon_jets[ix].momentum().pT()>80.*GeV &&
- recon_jets[iy].momentum().pT()>80.*GeV &&
- deltaR(recon_jets[ix],recon_jets[iy]) <0.6 )
- pass80DeltaR = false;
- }
- }
-
- // require at least four jets with et > 55
- if(njet55<=3) vetoEvent;
-
- // plots of etmiss/ht
- double etht = eTmiss/sqrt(HT);
- if(njet55==6) {
- _etmisspT_55_NJ_6_obs->fill(etht,weight);
- _etmisspT_55_NJ_6_bac->fill(etht,weight);
- _etmisspT_55_NJ_6_sig->fill(etht,weight);
- }
- else if(njet55==7) {
- _etmisspT_55_NJ_7_obs->fill(etht,weight);
- _etmisspT_55_NJ_7_bac->fill(etht,weight);
- _etmisspT_55_NJ_7_sig->fill(etht,weight);
- }
- else if(njet55==8) {
- _etmisspT_55_NJ_8_obs->fill(etht,weight);
- _etmisspT_55_NJ_8_bac->fill(etht,weight);
- _etmisspT_55_NJ_8_sig->fill(etht,weight);
- }
- if(njet80==5) {
- _etmisspT_80_NJ_5_obs->fill(etht,weight);
- _etmisspT_80_NJ_5_bac->fill(etht,weight);
- _etmisspT_80_NJ_5_sig->fill(etht,weight);
- }
- else if(njet80==6) {
- _etmisspT_80_NJ_6_obs->fill(etht,weight);
- _etmisspT_80_NJ_6_bac->fill(etht,weight);
- _etmisspT_80_NJ_6_sig->fill(etht,weight);
- }
- else if(njet80==7) {
- _etmisspT_80_NJ_7_obs->fill(etht,weight);
- _etmisspT_80_NJ_7_bac->fill(etht,weight);
- _etmisspT_80_NJ_7_sig->fill(etht,weight);
- }
-
- if(etht>1.5&ðt<2. ) {
- if(njet55>3) {
- _njet55A_obs->fill(njet55,weight);
- _njet55A_bac->fill(njet55,weight);
- _njet55A_sig->fill(njet55,weight);
+ // ATLAS calo problem
+ if(rand()/static_cast<double>(RAND_MAX)<=0.42) {
+ foreach ( const Particle & e, recon_e ) {
+ double eta = e.momentum().eta();
+ double phi = e.momentum().azimuthalAngle(MINUSPI_PLUSPI);
+ if(eta>-0.1&&eta<1.5&&phi>-0.9&&phi<-0.5)
+ vetoEvent;
}
- if(njet80>3) {
- _njet80A_obs->fill(njet80,weight);
- _njet80A_bac->fill(njet80,weight);
- _njet80A_sig->fill(njet80,weight);
+ foreach ( const Jet & jet, recon_jets ) {
+ double eta = jet.momentum().rapidity();
+ double phi = jet.momentum().azimuthalAngle(MINUSPI_PLUSPI);
+ if(jet.momentum().perp()>40 && eta>-0.1&&eta<1.5&&phi>-0.9&&phi<-0.5)
+ vetoEvent;
}
}
- else if(etht>2. &ðt<3. ) {
- if(njet55>3) {
- _njet55B_obs->fill(njet55,weight);
- _njet55B_bac->fill(njet55,weight);
- _njet55B_sig->fill(njet55,weight);
- }
- if(njet80>3) {
- _njet80B_obs->fill(njet80,weight);
- _njet80B_bac->fill(njet80,weight);
- _njet80B_sig->fill(njet80,weight);
+
+ // check at least one e/mu passing trigger
+ if( !( !recon_e .empty() && recon_e[0] .momentum().perp()>25.) &&
+ !( !recon_mu.empty() && recon_mu[0].momentum().perp()>20.) ) {
+ MSG_DEBUG("Hardest lepton fails trigger");
+ vetoEvent;
+ }
+
+ // eTmiss cut
+ if(eTmiss<50.) vetoEvent;
+
+ // check at least 1 SFOS pair
+ double mSFOS=1e30, mdiff=1e30;
+ unsigned int nSFOS=0;
+ for(unsigned int ix=0;ix<recon_e.size();++ix) {
+ for(unsigned int iy=ix+1;iy<recon_e.size();++iy) {
+ if(recon_e[ix].pdgId()*recon_e[iy].pdgId()>0) continue;
+ ++nSFOS;
+ double mtest = (recon_e[ix].momentum()+recon_e[iy].momentum()).mass();
+ // veto is mass<20
+ if(mtest<20.) vetoEvent;
+ if(fabs(mtest-90.)<mdiff) {
+ mSFOS = mtest;
+ mdiff = fabs(mtest-90.);
+ }
+ }
+ }
+ for(unsigned int ix=0;ix<recon_mu.size();++ix) {
+ for(unsigned int iy=ix+1;iy<recon_mu.size();++iy) {
+ if(recon_mu[ix].pdgId()*recon_mu[iy].pdgId()>0) continue;
+ ++nSFOS;
+ double mtest = (recon_mu[ix].momentum()+recon_mu[iy].momentum()).mass();
+ // veto is mass<20
+ if(mtest<20.) vetoEvent;
+ if(fabs(mtest-90.)<mdiff) {
+ mSFOS = mtest;
+ mdiff = fabs(mtest-90.);
+ }
+ }
+ }
+ // require at least 1 SFOS pair
+ if(nSFOS==0) vetoEvent;
+ // b-jet veto in SR!
+ if(mdiff>10.) {
+ foreach (const Jet & jet, recon_jets ) {
+ if(jet.containsBottom() &&
+ rand()/static_cast<double>(RAND_MAX)<=0.60)
+ vetoEvent;
}
}
+ // region SR1, Z depleted
+ if(mdiff>10.) {
+ _count_SR1->fill(0.5,weight);
+ _hist_etmiss_SR1_A->fill(eTmiss,weight);
+ _hist_etmiss_SR1_B->fill(eTmiss,weight);
+ _hist_mSFOS->fill(mSFOS,weight);
+ }
+ // region SR2, Z enriched
else {
- if(njet55>3) {
- _njet55C_obs->fill(njet55,weight);
- _njet55C_bac->fill(njet55,weight);
- _njet55C_sig->fill(njet55,weight);
+ _count_SR2->fill(0.5,weight);
+ _hist_etmiss_SR2_A->fill(eTmiss,weight);
+ _hist_etmiss_SR2_B->fill(eTmiss,weight);
+ }
+ // make the control plots
+ // lepton pT
+ unsigned int ie=0,imu=0;
+ for(unsigned int ix=0;ix<3;++ix) {
+ AIDA::IHistogram1D* hist = mdiff>10. ?
+ _hist_leptonpT_SR1[ix] : _hist_leptonpT_SR2[ix];
+ double pTe = ie <recon_e .size() ?
+ recon_e [ie ].momentum().perp() : -1*GeV;
+ double pTmu = imu<recon_mu.size() ?
+ recon_mu[imu].momentum().perp() : -1*GeV;
+ if(pTe>pTmu) {
+ hist->fill(pTe ,weight);
+ ++ie;
}
- if(njet80>3) {
- _njet80C_obs->fill(njet80,weight);
- _njet80C_bac->fill(njet80,weight);
- _njet80C_sig->fill(njet80,weight);
+ else {
+ hist->fill(pTmu,weight);
+ ++imu;
}
}
-
- // apply E_T/sqrt(H_T) cut
- if(etht<=3.5*GeV) {
- MSG_DEBUG("Fails ET/sqrt(HT) cut ");
- vetoEvent;
- }
-
- // check passes at least one delta5/ njet number cut
- if(!(pass55DeltaR && njet55 >= 7) &&
- !(pass80DeltaR && njet80 >= 6) ) {
- MSG_DEBUG("Fails DeltaR cut or jet number cuts");
- vetoEvent;
- }
-
- // 7j55
- if(njet55>=7&&pass55DeltaR)
- _count_7j55->fill( 0.5, weight) ;
- // 8j55
- if(njet55>=8&&pass55DeltaR)
- _count_8j55->fill( 0.5, weight) ;
- // 6j80
- if(njet80>=6&&pass80DeltaR)
- _count_6j80->fill( 0.5, weight) ;
- // 7j80
- if(njet80>=7&&pass80DeltaR)
- _count_7j80->fill( 0.5, weight) ;
-
}
//@}
void finalize() {
- double norm = crossSection()/femtobarn*1.34/sumOfWeights();
- scale(_etmisspT_55_NJ_6_obs,norm);
- scale(_etmisspT_55_NJ_6_bac,norm);
- scale(_etmisspT_55_NJ_6_sig,norm);
- scale(_etmisspT_55_NJ_7_obs,norm);
- scale(_etmisspT_55_NJ_7_bac,norm);
- scale(_etmisspT_55_NJ_7_sig,norm);
- scale(_etmisspT_55_NJ_8_obs,norm);
- scale(_etmisspT_55_NJ_8_bac,norm);
- scale(_etmisspT_55_NJ_8_sig,norm);
- scale(_etmisspT_80_NJ_5_obs,norm);
- scale(_etmisspT_80_NJ_5_bac,norm);
- scale(_etmisspT_80_NJ_5_sig,norm);
- scale(_etmisspT_80_NJ_6_obs,norm);
- scale(_etmisspT_80_NJ_6_bac,norm);
- scale(_etmisspT_80_NJ_6_sig,norm);
- scale(_etmisspT_80_NJ_7_obs,norm);
- scale(_etmisspT_80_NJ_7_bac,norm);
- scale(_etmisspT_80_NJ_7_sig,norm);
- scale(_njet55A_obs,norm);
- scale(_njet55A_bac,norm);
- scale(_njet55A_sig,norm);
- scale(_njet55B_obs,norm);
- scale(_njet55B_bac,norm);
- scale(_njet55B_sig,norm);
- scale(_njet55C_obs,norm);
- scale(_njet55C_bac,norm);
- scale(_njet55C_sig,norm);
- scale(_njet80A_obs,norm);
- scale(_njet80A_bac,norm);
- scale(_njet80A_sig,norm);
- scale(_njet80B_obs,norm);
- scale(_njet80B_bac,norm);
- scale(_njet80B_sig,norm);
- scale(_njet80C_obs,norm);
- scale(_njet80C_bac,norm);
- scale(_njet80C_sig,norm);
- scale(_count_7j55,norm);
- scale(_count_8j55,norm);
- scale(_count_6j80,norm);
- scale(_count_7j80,norm);
+ double norm = crossSection()/femtobarn*2.06/sumOfWeights();
+ // these are number of events at 2.06fb^-1 per 20 GeV
+ for(unsigned int ix=0;ix<3;++ix) {
+ scale(_hist_leptonpT_SR1[ix],norm*20.);
+ scale(_hist_leptonpT_SR2[ix],norm*20.);
+ }
+ scale(_hist_etmiss_SR1_A,norm*20.);
+ scale(_hist_etmiss_SR1_B,norm*20.);
+ scale(_hist_etmiss_SR2_A,norm*20.);
+ scale(_hist_etmiss_SR2_B,norm*20.);
+ scale(_hist_mSFOS ,norm*20.);
+ // these are number of events at 2.06fb^-1
+ scale(_count_SR1,norm);
+ scale(_count_SR2,norm);
}
private:
- /// @name Histograms
- //@{
- AIDA::IHistogram1D* _etmisspT_55_NJ_6_obs;
- AIDA::IHistogram1D* _etmisspT_55_NJ_6_bac;
- AIDA::IHistogram1D* _etmisspT_55_NJ_6_sig;
- AIDA::IHistogram1D* _etmisspT_55_NJ_7_obs;
- AIDA::IHistogram1D* _etmisspT_55_NJ_7_bac;
- AIDA::IHistogram1D* _etmisspT_55_NJ_7_sig;
- AIDA::IHistogram1D* _etmisspT_55_NJ_8_obs;
- AIDA::IHistogram1D* _etmisspT_55_NJ_8_bac;
- AIDA::IHistogram1D* _etmisspT_55_NJ_8_sig;
- AIDA::IHistogram1D* _etmisspT_80_NJ_5_obs;
- AIDA::IHistogram1D* _etmisspT_80_NJ_5_bac;
- AIDA::IHistogram1D* _etmisspT_80_NJ_5_sig;
- AIDA::IHistogram1D* _etmisspT_80_NJ_6_obs;
- AIDA::IHistogram1D* _etmisspT_80_NJ_6_bac;
- AIDA::IHistogram1D* _etmisspT_80_NJ_6_sig;
- AIDA::IHistogram1D* _etmisspT_80_NJ_7_obs;
- AIDA::IHistogram1D* _etmisspT_80_NJ_7_bac;
- AIDA::IHistogram1D* _etmisspT_80_NJ_7_sig;
- AIDA::IHistogram1D* _njet55A_obs;
- AIDA::IHistogram1D* _njet55A_bac;
- AIDA::IHistogram1D* _njet55A_sig;
- AIDA::IHistogram1D* _njet55B_obs;
- AIDA::IHistogram1D* _njet55B_bac;
- AIDA::IHistogram1D* _njet55B_sig;
- AIDA::IHistogram1D* _njet55C_obs;
- AIDA::IHistogram1D* _njet55C_bac;
- AIDA::IHistogram1D* _njet55C_sig;
- AIDA::IHistogram1D* _njet80A_obs;
- AIDA::IHistogram1D* _njet80A_bac;
- AIDA::IHistogram1D* _njet80A_sig;
- AIDA::IHistogram1D* _njet80B_obs;
- AIDA::IHistogram1D* _njet80B_bac;
- AIDA::IHistogram1D* _njet80B_sig;
- AIDA::IHistogram1D* _njet80C_obs;
- AIDA::IHistogram1D* _njet80C_bac;
- AIDA::IHistogram1D* _njet80C_sig;
- AIDA::IHistogram1D* _count_7j55;
- AIDA::IHistogram1D* _count_8j55;
- AIDA::IHistogram1D* _count_6j80;
- AIDA::IHistogram1D* _count_7j80;
- //@}
+ /// @name Histograms
+ //@{
+ vector<AIDA::IHistogram1D*> _hist_leptonpT_SR1;
+ vector<AIDA::IHistogram1D*> _hist_leptonpT_SR2;
+ AIDA::IHistogram1D* _hist_etmiss_SR1_A;
+ AIDA::IHistogram1D* _hist_etmiss_SR1_B;
+ AIDA::IHistogram1D* _hist_etmiss_SR2_A;
+ AIDA::IHistogram1D* _hist_etmiss_SR2_B;
+ AIDA::IHistogram1D* _hist_mSFOS;
+ AIDA::IHistogram1D* _count_SR1;
+ AIDA::IHistogram1D* _count_SR2;
+ //@}
};
// The hook for the plugin system
- DECLARE_RIVET_PLUGIN(ATLAS_2011_S9225137);
+ DECLARE_RIVET_PLUGIN(ATLAS_2012_CONF_2012_023);
}
Modified: trunk/src/Analyses/Makefile.am
==============================================================================
--- trunk/src/Analyses/Makefile.am Mon Apr 16 13:14:25 2012 (r3685)
+++ trunk/src/Analyses/Makefile.am Mon Apr 16 13:16:53 2012 (r3686)
@@ -79,6 +79,7 @@
ATLAS_2012_I1083318.cc \
ATLAS_2012_I1095236.cc \
ATLAS_2012_CONF_2012_001.cc \
+ ATLAS_2012_CONF_2012_023.cc \
ATLAS_2012_CONF_2012_033.cc \
ATLAS_2012_CONF_2012_037.cc
endif
More information about the Rivet-svn
mailing list
|