|
[Rivet-svn] r3683 - in trunk: data/anainfo data/plotinfo src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgMon Apr 16 13:00:38 BST 2012
Author: richardn Date: Mon Apr 16 13:00:38 2012 New Revision: 3683 Log: added ATLAS_2012_CONF_2012_037, latest 0 lepton + met + high jet multiplicity SUSY search Added: trunk/data/anainfo/ATLAS_2012_CONF_2012_037.info - copied, changed from r3651, trunk/data/anainfo/ATLAS_2011_S9225137.info trunk/data/plotinfo/ATLAS_2012_CONF_2012_037.plot - copied, changed from r3651, trunk/data/plotinfo/ATLAS_2011_S9225137.plot trunk/src/Analyses/ATLAS_2012_CONF_2012_037.cc - copied, changed from r3651, 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_037.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_037.info Mon Apr 16 13:00:38 2012 (r3683) @@ -1,14 +1,13 @@ -Name: ATLAS_2011_S9225137 -Year: 2011 +Name: ATLAS_2012_CONF_2012_037 +Year: 2012 Summary: High jet multiplicity squark and gluino search Experiment: ATLAS Collider: LHC -SpiresID: 9225137 Status: UNVALIDATED Authors: - Peter Richardson <peter.richardson at durham.ac.uk> References: - - arXiv:1110.2299 + - ATLAS-CONF-2012-037 RunInfo: BSM signal events at 7000 GeV. NumEvents: 25000 for BSM signals @@ -19,25 +18,7 @@ 'Search for SUSY using events with 6 or more jets in association with missing transverse momentum produced 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 + integrated luminosity of 4.7 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;%%", -}' Modified: trunk/data/anainfo/Makefile.am ============================================================================== --- trunk/data/anainfo/Makefile.am Mon Apr 16 12:57:08 2012 (r3682) +++ trunk/data/anainfo/Makefile.am Mon Apr 16 13:00:38 2012 (r3683) @@ -45,6 +45,7 @@ ATLAS_2012_I1084540.info \ ATLAS_2012_I1094568.info \ ATLAS_2012_CONF_2012_033.info \ + ATLAS_2012_CONF_2012_037.info \ BABAR_2005_S6181155.info \ BABAR_2007_S6895344.info \ BABAR_2007_S7266081.info \ Copied and modified: trunk/data/plotinfo/ATLAS_2012_CONF_2012_037.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_037.plot Mon Apr 16 13:00:38 2012 (r3683) @@ -1,241 +1,87 @@ -# BEGIN PLOT /ATLAS_2011_S9225137/count_7j55 +# BEGIN PLOT /ATLAS_2012_CONF_2012_037/count_7j55 Title=Number of events in 7 jet, $p_\perp>55$\,GeV signal region XLabel= YLabel=Number of Events -YLog=0 +LogY=0 # END PLOT -# BEGIN PLOT /ATLAS_2011_S9225137/count_8j55 +# BEGIN PLOT /ATLAS_2012_CONF_2012_037/count_8j55 Title=Number of events in 8 jet, $p_\perp>55$\,GeV signal region XLabel= YLabel=Number of Events -YLog=0 +LogY=0 # END PLOT -# BEGIN PLOT /ATLAS_2011_S9225137/count_6j80 +# BEGIN PLOT /ATLAS_2012_CONF_2012_037/count_9j55 +Title=Number of events in 9 jet, $p_\perp>55$\,GeV signal region +XLabel= +YLabel=Number of Events +LogY=0 +# END PLOT +# BEGIN PLOT /ATLAS_2012_CONF_2012_037/count_6j80 Title=Number of events in 6 jet, $p_\perp>80$\,GeV signal region XLabel= YLabel=Number of Events -YLog=0 +LogY=0 # END PLOT -# BEGIN PLOT /ATLAS_2011_S9225137/count_7j80 +# BEGIN PLOT /ATLAS_2012_CONF_2012_037/count_7j80 Title=Number of events in 7 jet, $p_\perp>80$\,GeV signal region XLabel= YLabel=Number of Events -YLog=0 +LogY=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 +# BEGIN PLOT /ATLAS_2012_CONF_2012_038/count_8j80 +Title=Number of events in 8 jet, $p_\perp>80$\,GeV signal region +XLabel= +YLabel=Number of Events +LogY=0 # 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 + +# BEGIN PLOT /ATLAS_2012_CONF_2012_037/hist_njet55 +Title=Observed number of jets with $p_\perp>55$\,GeV for $E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}>4\,\mathrm{GeV}^{\frac12}$ +XLabel=Number of Jets +YLabel=Number of events +LogY=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 +# BEGIN PLOT /ATLAS_2012_CONF_2012_037/hist_njet80 +Title=Observed number of jets with $p_\perp>80$\,GeV for $E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}>4\,\mathrm{GeV}^{\frac12}$ +XLabel=Number of Jets +YLabel=Number of events +LogY=1 # END PLOT -# BEGIN PLOT /ATLAS_2011_S9225137/d13-x01-y01 + +# BEGIN PLOT /ATLAS_2012_CONF_2012_037/etmiss_HT_7j55 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 +YLabel=$\mathrm{Events}/2$\,GeV +LogY=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 +# BEGIN PLOT /ATLAS_2012_CONF_2012_037/etmiss_HT_8j55 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 +YLabel=$\mathrm{Events}/2$\,GeV +LogY=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 +# BEGIN PLOT /ATLAS_2012_CONF_2012_037/etmiss_HT_9j55 +Title=Observed $E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ for 9 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 +YLabel=$\mathrm{Events}/2$\,GeV +LogY=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 + +# BEGIN PLOT /ATLAS_2012_CONF_2012_037/etmiss_HT_6j80 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 +YLabel=$\mathrm{Events}/2$\,GeV +LogY=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 +# BEGIN PLOT /ATLAS_2012_CONF_2012_037/etmiss_HT_7j80 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 +YLabel=$\mathrm{Events}/2$\,GeV +LogY=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 +# BEGIN PLOT /ATLAS_2012_CONF_2012_037/etmiss_HT_8j80 +Title=Observed $E_{\mathrm{T}}/\sqrt{H_{\mathrm{T}}}$ for 8 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 +YLabel=$\mathrm{Events}/2$\,GeV +LogY=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 -# END PLOT - Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Mon Apr 16 12:57:08 2012 (r3682) +++ trunk/data/plotinfo/Makefile.am Mon Apr 16 13:00:38 2012 (r3683) @@ -42,6 +42,7 @@ ATLAS_2012_I1084540.plot \ ATLAS_2012_I1094568.plot \ ATLAS_2012_CONF_2012_033.plot \ + ATLAS_2012_CONF_2012_037.plot \ BABAR_2005_S6181155.plot \ BABAR_2007_S6895344.plot \ BABAR_2007_S7266081.plot \ Copied and modified: trunk/src/Analyses/ATLAS_2012_CONF_2012_037.cc (from r3651, trunk/src/Analyses/ATLAS_2011_S9225137.cc) ============================================================================== --- trunk/src/Analyses/ATLAS_2011_S9225137.cc Sun Mar 25 11:21:04 2012 (r3651, copy source) +++ trunk/src/Analyses/ATLAS_2012_CONF_2012_037.cc Mon Apr 16 13:00:38 2012 (r3683) @@ -14,15 +14,15 @@ namespace Rivet { - class ATLAS_2011_S9225137 : public Analysis { + class ATLAS_2012_CONF_2012_037 : public Analysis { public: /// @name Constructors etc. //@{ /// Constructor - ATLAS_2011_S9225137() - : Analysis("ATLAS_2011_S9225137") + ATLAS_2012_CONF_2012_037() + : Analysis("ATLAS_2012_CONF_2012_037") { } //@} @@ -36,14 +36,6 @@ /// 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)); @@ -72,48 +64,22 @@ 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); + _etmiss_HT_7j55 = bookHistogram1D("etmiss_HT_7j55", 8, 0., 16.); + _etmiss_HT_8j55 = bookHistogram1D("etmiss_HT_8j55", 8, 0., 16.); + _etmiss_HT_9j55 = bookHistogram1D("etmiss_HT_9j55", 8, 0., 16.); + _etmiss_HT_6j80 = bookHistogram1D("etmiss_HT_6j80", 8, 0., 16.); + _etmiss_HT_7j80 = bookHistogram1D("etmiss_HT_7j80", 8, 0., 16.); + _etmiss_HT_8j80 = bookHistogram1D("etmiss_HT_8j80", 8, 0., 16.); + + _hist_njet55 = bookHistogram1D("hist_njet55", 11, 2.5, 13.5); + _hist_njet80 = bookHistogram1D("hist_njet80", 11, 2.5, 13.5); _count_7j55 = bookHistogram1D("count_7j55", 1, 0., 1.); _count_8j55 = bookHistogram1D("count_8j55", 1, 0., 1.); + _count_9j55 = bookHistogram1D("count_9j55", 1, 0., 1.); _count_6j80 = bookHistogram1D("count_6j80", 1, 0., 1.); _count_7j80 = bookHistogram1D("count_7j80", 1, 0., 1.); + _count_8j80 = bookHistogram1D("count_8j80", 1, 0., 1.); } @@ -121,19 +87,12 @@ /// 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,35 +113,29 @@ } // candidate electrons - ParticleVector cand_e = applyProjection<IdentifiedFinalState>(event, "elecs").particlesByPt(); // resolve jet/lepton ambiguity - Jets cand_jets_2; + 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 ); - } + // candidates after |eta| < 2.8 + 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 ); } // 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; @@ -196,7 +149,7 @@ 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; @@ -215,13 +168,6 @@ } 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 @@ -237,223 +183,106 @@ HT += jet.momentum().pT() ; } - // number of jets and deltaR - bool pass55DeltaR=true; - unsigned int njet55=0; - bool pass80DeltaR=true; - unsigned int njet80=0; + // number of jets + unsigned int njet55=0, 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; + double ratio = eTmiss/sqrt(HT); - // 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(ratio>4.) { + _hist_njet55->fill(njet55,weight); + _hist_njet80->fill(njet80,weight); + // 7j55 + if(njet55>=7) + _count_7j55->fill( 0.5, weight); + // 8j55 + if(njet55>=8) + _count_8j55->fill( 0.5, weight) ; + // 8j55 + if(njet55>=9) + _count_9j55->fill( 0.5, weight) ; + // 6j80 + if(njet80>=6) + _count_6j80->fill( 0.5, weight) ; + // 7j80 + if(njet80>=7) + _count_7j80->fill( 0.5, weight) ; + // 8j80 + if(njet80>=7) + _count_8j80->fill( 0.5, 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); - } - if(njet80>3) { - _njet80A_obs->fill(njet80,weight); - _njet80A_bac->fill(njet80,weight); - _njet80A_sig->fill(njet80,weight); - } - } - 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); - } - } - else { - if(njet55>3) { - _njet55C_obs->fill(njet55,weight); - _njet55C_bac->fill(njet55,weight); - _njet55C_sig->fill(njet55,weight); - } - if(njet80>3) { - _njet80C_obs->fill(njet80,weight); - _njet80C_bac->fill(njet80,weight); - _njet80C_sig->fill(njet80,weight); - } - } - - // 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) ; + if(njet55>=7) + _etmiss_HT_7j55->fill( ratio, weight); + // 8j55 + if(njet55>=8) + _etmiss_HT_8j55->fill( ratio, weight) ; // 8j55 - if(njet55>=8&&pass55DeltaR) - _count_8j55->fill( 0.5, weight) ; + if(njet55>=9) + _etmiss_HT_9j55->fill( ratio, weight) ; // 6j80 - if(njet80>=6&&pass80DeltaR) - _count_6j80->fill( 0.5, weight) ; + if(njet80>=6) + _etmiss_HT_6j80->fill( ratio, weight) ; // 7j80 - if(njet80>=7&&pass80DeltaR) - _count_7j80->fill( 0.5, weight) ; + if(njet80>=7) + _etmiss_HT_7j80->fill( ratio, weight) ; + // 8j80 + if(njet80>=7) + _etmiss_HT_8j80->fill( ratio, 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); + double norm = crossSection()/femtobarn*4.7/sumOfWeights(); + + scale(_etmiss_HT_7j55,2.*norm); + scale(_etmiss_HT_8j55,2.*norm); + scale(_etmiss_HT_9j55,2.*norm); + scale(_etmiss_HT_6j80,2.*norm); + scale(_etmiss_HT_7j80,2.*norm); + scale(_etmiss_HT_8j80,2.*norm); + + scale(_hist_njet55,norm); + scale(_hist_njet80,norm); + scale(_count_7j55,norm); scale(_count_8j55,norm); + scale(_count_9j55,norm); scale(_count_6j80,norm); scale(_count_7j80,norm); + scale(_count_8j80,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* _etmiss_HT_7j55; + AIDA::IHistogram1D* _etmiss_HT_8j55; + AIDA::IHistogram1D* _etmiss_HT_9j55; + AIDA::IHistogram1D* _etmiss_HT_6j80; + AIDA::IHistogram1D* _etmiss_HT_7j80; + AIDA::IHistogram1D* _etmiss_HT_8j80; + + AIDA::IHistogram1D* _hist_njet55; + AIDA::IHistogram1D* _hist_njet80; + AIDA::IHistogram1D* _count_7j55; AIDA::IHistogram1D* _count_8j55; + AIDA::IHistogram1D* _count_9j55; AIDA::IHistogram1D* _count_6j80; AIDA::IHistogram1D* _count_7j80; + AIDA::IHistogram1D* _count_8j80; //@} }; // The hook for the plugin system - DECLARE_RIVET_PLUGIN(ATLAS_2011_S9225137); + DECLARE_RIVET_PLUGIN(ATLAS_2012_CONF_2012_037); } Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Mon Apr 16 12:57:08 2012 (r3682) +++ trunk/src/Analyses/Makefile.am Mon Apr 16 13:00:38 2012 (r3683) @@ -77,7 +77,8 @@ ATLAS_2011_S9108483.cc \ ATLAS_2011_S9225137.cc \ ATLAS_2012_I1083318.cc \ - ATLAS_2012_CONF_2012_033.cc + ATLAS_2012_CONF_2012_033.cc \ + ATLAS_2012_CONF_2012_037.cc endif
More information about the Rivet-svn mailing list |