|
[Rivet-svn] r3685 - in trunk: data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgMon Apr 16 13:14:25 BST 2012
Author: richardn Date: Mon Apr 16 13:14:25 2012 New Revision: 3685 Log: added ATLAS_2012_CONF_2012_001, latest 4 lepton + MET SUSY search Added: trunk/data/anainfo/ATLAS_2012_CONF_2012_001.info - copied, changed from r3651, trunk/data/anainfo/ATLAS_2011_S9225137.info trunk/data/plotinfo/ATLAS_2012_CONF_2012_001.plot - copied, changed from r3651, trunk/data/plotinfo/ATLAS_2011_S9225137.plot trunk/data/refdata/ATLAS_2012_CONF_2012_001.aida trunk/src/Analyses/ATLAS_2012_CONF_2012_001.cc - copied, changed from r3647, trunk/src/Analyses/ATLAS_2011_S9225137.cc Modified: trunk/data/anainfo/Makefile.am trunk/data/plotinfo/Makefile.am trunk/data/refdata/Makefile.am trunk/src/Analyses/Makefile.am Copied and modified: trunk/data/anainfo/ATLAS_2012_CONF_2012_001.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_001.info Mon Apr 16 13:14:25 2012 (r3685) @@ -1,14 +1,13 @@ -Name: ATLAS_2011_S9225137 +Name: ATLAS_2012_CONF_2012_001 Year: 2011 -Summary: High jet multiplicity squark and gluino search +Summary: 4 or more 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-001 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 4 or more 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:10:22 2012 (r3684) +++ trunk/data/anainfo/Makefile.am Mon Apr 16 13:14:25 2012 (r3685) @@ -45,6 +45,7 @@ ATLAS_2012_I1084540.info \ ATLAS_2012_I1094568.info \ ATLAS_2012_I1095236.info \ + ATLAS_2012_CONF_2012_001.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_001.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_001.plot Mon Apr 16 13:14:25 2012 (r3685) @@ -1,241 +1,62 @@ -# 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_001/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_001/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 -# 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 + +# BEGIN PLOT /ATLAS_2012_CONF_2012_001/d01-x01-y01 +Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the leading lepton +XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV] +YLabel=$\mathrm{Events}/10 \mathrm{GeV}$ +LogY=1 +# END PLOT +# BEGIN PLOT /ATLAS_2012_CONF_2012_001/d02-x01-y01 +Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the $2^{\rm nd}$ leading lepton +XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV] +YLabel=$\mathrm{Events}/10 \mathrm{GeV}$ +LogY=1 +# END PLOT +# BEGIN PLOT /ATLAS_2012_CONF_2012_001/d03-x01-y01 +Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the $3^{\rm rd}$ leading lepton +XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV] +YLabel=$\mathrm{Events}/5 \mathrm{GeV}$ +LogY=1 +# END PLOT +# BEGIN PLOT /ATLAS_2012_CONF_2012_001/d04-x01-y01 +Title=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ for the $4^{\rm th}$ leading lepton +XLabel=$E^e_{\mathrm{T}}(p^\mu_{\mathrm{T}})$ [GeV] +YLabel=$\mathrm{Events}/5 \mathrm{GeV}$ +LogY=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}$ + +# BEGIN PLOT /ATLAS_2012_CONF_2012_001/d05-x01-y01 +Title=Number of Jets XLabel=Number of Jets -YLabel=Number of events -YLog=1 +YLabel=Events +LogY=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_001/d06-x01-y01 +Title=Missing Transverse Energy +XLabel=$E^{\rm miss}_{\mathrm{T}}$ [GeV] +YLabel=$\mathrm{Events}/20$\,GeV +LogY=1 +# END PLOT +# BEGIN PLOT /ATLAS_2012_CONF_2012_001/d07-x01-y01 +Title=Mass of SFOF lepton pair closest to the $Z^0$ mass +XLabel=$M_{\mathrm{SFOS}}$ [GeV] +YLabel=$\mathrm{Events}/20$\,GeV +LogY=1 +# END PLOT +# BEGIN PLOT /ATLAS_2012_CONF_2012_001/d08-x01-y01 +Title=Effective Mass +XLabel=$M_{\mathrm{eff}}$ [GeV] +YLabel=$\mathrm{Events}/50$\,GeV +LogY=1 # END PLOT - Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Mon Apr 16 13:10:22 2012 (r3684) +++ trunk/data/plotinfo/Makefile.am Mon Apr 16 13:14:25 2012 (r3685) @@ -42,6 +42,7 @@ ATLAS_2012_I1084540.plot \ ATLAS_2012_I1094568.plot \ ATLAS_2012_I1095236.plot \ + ATLAS_2012_CONF_2012_001.plot \ ATLAS_2012_CONF_2012_033.plot \ ATLAS_2012_CONF_2012_037.plot \ BABAR_2005_S6181155.plot \ Added: trunk/data/refdata/ATLAS_2012_CONF_2012_001.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/ATLAS_2012_CONF_2012_001.aida Mon Apr 16 13:14:25 2012 (r3685) @@ -0,0 +1,709 @@ +<?xml version="1.0" encoding="utf-8" ?> +<!DOCTYPE aida SYSTEM "http://aida.freehep.org/schemas/3.3/aida.dtd"> +<aida version="3.3"> + <implementation version="1.0" package="HepData"/> + <dataPointSet name="d01-x01-y01" dimension="2" path="/REF/ATLAS_2012_CONF_2012_001" title="EVENTS" > + <dataPoint> + <measurement value="35.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="3.0" errorPlus="2.9" errorMinus="1.63"/> + </dataPoint> + <dataPoint> + <measurement value="45.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="5.0" errorPlus="3.37" errorMinus="2.15"/> + </dataPoint> + <dataPoint> + <measurement value="55.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="2.0" errorPlus="2.62" errorMinus="1.29"/> + </dataPoint> + <dataPoint> + <measurement value="65.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="1.0" errorPlus="2.29" errorMinus="0.83"/> + </dataPoint> + <dataPoint> + <measurement value="75.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="5.0" errorPlus="3.37" errorMinus="2.15"/> + </dataPoint> + <dataPoint> + <measurement value="85.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="3.0" errorPlus="2.9" errorMinus="1.63"/> + </dataPoint> + <dataPoint> + <measurement value="105.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="3.0" errorPlus="2.9" errorMinus="1.63"/> + </dataPoint> + <dataPoint> + <measurement value="115.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="1.0" errorPlus="2.29" errorMinus="0.83"/> + </dataPoint> + <dataPoint> + <measurement value="135.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="1.0" errorPlus="2.29" errorMinus="0.83"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d01-x01-y02" dimension="2" path="/REF/ATLAS_2012_CONF_2012_001" title="EVENTS" > + <dataPoint> + <measurement value="25.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.583" errorPlus="0.557" errorMinus="0.557"/> + </dataPoint> + <dataPoint> + <measurement value="35.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.678" errorPlus="0.24030189345904038" errorMinus="0.24030189345904038"/> + </dataPoint> + <dataPoint> + <measurement value="45.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="3.82" errorPlus="1.19" errorMinus="1.19"/> + </dataPoint> + <dataPoint> + <measurement value="55.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="4.66" errorPlus="1.1898319209031165" errorMinus="1.1898319209031165"/> + </dataPoint> + <dataPoint> + <measurement value="65.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="3.9" errorPlus="1.164817582284883" errorMinus="1.164817582284883"/> + </dataPoint> + <dataPoint> + <measurement value="75.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="2.96" errorPlus="0.724016574395918" errorMinus="0.724016574395918"/> + </dataPoint> + <dataPoint> + <measurement value="85.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="1.98" errorPlus="0.41182520563948" errorMinus="0.41182520563948"/> + </dataPoint> + <dataPoint> + <measurement value="95.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="1.85" errorPlus="0.604400529450463" errorMinus="0.604400529450463"/> + </dataPoint> + <dataPoint> + <measurement value="105.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="1.1" errorPlus="0.45880278987817846" errorMinus="0.45880278987817846"/> + </dataPoint> + <dataPoint> + <measurement value="115.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.495" errorPlus="0.13577186748365805" errorMinus="0.13577186748365805"/> + </dataPoint> + <dataPoint> + <measurement value="125.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.502" errorPlus="0.16324827717314508" errorMinus="0.16324827717314508"/> + </dataPoint> + <dataPoint> + <measurement value="135.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.497" errorPlus="0.1704230031421815" errorMinus="0.1704230031421815"/> + </dataPoint> + <dataPoint> + <measurement value="145.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.4" errorPlus="0.1492816130673835" errorMinus="0.1492816130673835"/> + </dataPoint> + <dataPoint> + <measurement value="155.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.206" errorPlus="0.09479451460923254" errorMinus="0.09479451460923254"/> + </dataPoint> + <dataPoint> + <measurement value="165.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.19" errorPlus="0.09990495483208028" errorMinus="0.09990495483208028"/> + </dataPoint> + <dataPoint> + <measurement value="175.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.361" errorPlus="0.5844287809476875" errorMinus="0.5844287809476875"/> + </dataPoint> + <dataPoint> + <measurement value="185.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.0474" errorPlus="0.4851555729866452" errorMinus="0.4851555729866452"/> + </dataPoint> + <dataPoint> + <measurement value="195.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.0191" errorPlus="0.01356244815658294" errorMinus="0.01356244815658294"/> + </dataPoint> + <dataPoint> + <measurement value="205.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.0646" errorPlus="0.03462080299473136" errorMinus="0.03462080299473136"/> + </dataPoint> + <dataPoint> + <measurement value="215.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.0392" errorPlus="0.030026821343592136" errorMinus="0.030026821343592136"/> + </dataPoint> + <dataPoint> + <measurement value="225.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.0127" errorPlus="0.032106852851065924" errorMinus="0.032106852851065924"/> + </dataPoint> + <dataPoint> + <measurement value="235.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.0298" errorPlus="0.029841246622753545" errorMinus="0.029841246622753545"/> + </dataPoint> + <dataPoint> + <measurement value="245.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.00572" errorPlus="0.006391189247706565" errorMinus="0.006391189247706565"/> + </dataPoint> + <dataPoint> + <measurement value="255.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.00215" errorPlus="0.002354591259645716" errorMinus="0.002354591259645716"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y01" dimension="2" path="/REF/ATLAS_2012_CONF_2012_001" title="EVENTS" > + <dataPoint> + <measurement value="15.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="3.0" errorPlus="2.9" errorMinus="1.63"/> + </dataPoint> + <dataPoint> + <measurement value="25.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="1.0" errorPlus="2.29" errorMinus="0.83"/> + </dataPoint> + <dataPoint> + <measurement value="35.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="1.0" errorPlus="2.29" errorMinus="0.83"/> + </dataPoint> + <dataPoint> + <measurement value="45.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="6.0" errorPlus="3.57" errorMinus="2.37"/> + </dataPoint> + <dataPoint> + <measurement value="55.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="4.0" errorPlus="3.15" errorMinus="1.91"/> + </dataPoint> + <dataPoint> + <measurement value="65.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="4.0" errorPlus="3.15" errorMinus="1.91"/> + </dataPoint> + <dataPoint> + <measurement value="75.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="4.0" errorPlus="3.15" errorMinus="1.91"/> + </dataPoint> + <dataPoint> + <measurement value="95.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="1.0" errorPlus="2.29" errorMinus="0.83"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y02" dimension="2" path="/REF/ATLAS_2012_CONF_2012_001" title="EVENTS" > + <dataPoint> + <measurement value="15.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.611" errorPlus="0.562594880886771" errorMinus="0.562594880886771"/> + </dataPoint> + <dataPoint> + <measurement value="25.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.995" errorPlus="0.3347312354710866" errorMinus="0.3347312354710866"/> + </dataPoint> + <dataPoint> + <measurement value="35.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="3.34" errorPlus="0.898276126811795" errorMinus="0.898276126811795"/> + </dataPoint> + <dataPoint> + <measurement value="45.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="7.27" errorPlus="1.61245154965971" errorMinus="1.61245154965971"/> + </dataPoint> + <dataPoint> + <measurement value="55.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="5.05" errorPlus="0.9538343671728337" errorMinus="0.9538343671728337"/> + </dataPoint> + <dataPoint> + <measurement value="65.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="2.56" errorPlus="0.5232590180780452" errorMinus="0.5232590180780452"/> + </dataPoint> + <dataPoint> + <measurement value="75.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="1.5" errorPlus="0.37589892258425006" errorMinus="0.37589892258425006"/> + </dataPoint> + <dataPoint> + <measurement value="85.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="1.18" errorPlus="0.4167733196834941" errorMinus="0.4167733196834941"/> + </dataPoint> + <dataPoint> + <measurement value="95.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.566" errorPlus="0.16606625183943907" errorMinus="0.16606625183943907"/> + </dataPoint> + <dataPoint> + <measurement value="105.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.432" errorPlus="0.1487346630748865" errorMinus="0.1487346630748865"/> + </dataPoint> + <dataPoint> + <measurement value="115.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.413" errorPlus="0.14354441821262157" errorMinus="0.14354441821262157"/> + </dataPoint> + <dataPoint> + <measurement value="125.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.173" errorPlus="0.07615773105863909" errorMinus="0.07615773105863909"/> + </dataPoint> + <dataPoint> + <measurement value="135.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.08" errorPlus="0.05817636977330228" errorMinus="0.05817636977330228"/> + </dataPoint> + <dataPoint> + <measurement value="145.0" errorPlus="5.0" errorMinus="5.0"/> + <measurement value="0.162" errorPlus="0.08563293758829017" errorMinus="0.08563293758829017"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d03-x01-y01" dimension="2" path="/REF/ATLAS_2012_CONF_2012_001" title="EVENTS" > + <dataPoint> + <measurement value="12.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="2.0" errorPlus="2.62" errorMinus="1.29"/> + </dataPoint> + <dataPoint> + <measurement value="17.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="1.0" errorPlus="2.29" errorMinus="0.83"/> + </dataPoint> + <dataPoint> + <measurement value="22.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="2.0" errorPlus="2.62" errorMinus="1.29"/> + </dataPoint> + <dataPoint> + <measurement value="27.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="3.0" errorPlus="2.9" errorMinus="1.63"/> + </dataPoint> + <dataPoint> + <measurement value="32.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="4.0" errorPlus="3.15" errorMinus="1.91"/> + </dataPoint> + <dataPoint> + <measurement value="37.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="4.0" errorPlus="3.15" errorMinus="1.91"/> + </dataPoint> + <dataPoint> + <measurement value="42.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="5.0" errorPlus="3.37" errorMinus="2.15"/> + </dataPoint> + <dataPoint> + <measurement value="47.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="2.0" errorPlus="2.62" errorMinus="1.29"/> + </dataPoint> + <dataPoint> + <measurement value="62.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="1.0" errorPlus="2.29" errorMinus="0.83"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d03-x01-y02" dimension="2" path="/REF/ATLAS_2012_CONF_2012_001" title="EVENTS" > + <dataPoint> + <measurement value="12.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.992" errorPlus="0.6235230549065528" errorMinus="0.6235230549065528"/> + </dataPoint> + <dataPoint> + <measurement value="17.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="1.54" errorPlus="0.7409453421137082" errorMinus="0.7409453421137082"/> + </dataPoint> + <dataPoint> + <measurement value="22.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="2.21" errorPlus="1.0183319694480772" errorMinus="1.0183319694480772"/> + </dataPoint> + <dataPoint> + <measurement value="27.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="3.72" errorPlus="1.3853880322855399" errorMinus="1.3853880322855399"/> + </dataPoint> + <dataPoint> + <measurement value="32.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="4.41" errorPlus="0.9060353193998565" errorMinus="0.9060353193998565"/> + </dataPoint> + <dataPoint> + <measurement value="37.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="3.22" errorPlus="0.6447480127925947" errorMinus="0.6447480127925947"/> + </dataPoint> + <dataPoint> + <measurement value="42.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="2.88" errorPlus="0.5769748694700663" errorMinus="0.5769748694700663"/> + </dataPoint> + <dataPoint> + <measurement value="47.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="1.79" errorPlus="0.37735924528226417" errorMinus="0.37735924528226417"/> + </dataPoint> + <dataPoint> + <measurement value="52.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="1.27" errorPlus="0.29410882339705485" errorMinus="0.29410882339705485"/> + </dataPoint> + <dataPoint> + <measurement value="57.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.714" errorPlus="0.18693849255838135" errorMinus="0.18693849255838135"/> + </dataPoint> + <dataPoint> + <measurement value="62.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.443" errorPlus="0.1492816130673835" errorMinus="0.1492816130673835"/> + </dataPoint> + <dataPoint> + <measurement value="67.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.256" errorPlus="0.09762171889492624" errorMinus="0.09762171889492624"/> + </dataPoint> + <dataPoint> + <measurement value="72.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.24" errorPlus="0.11715374513859982" errorMinus="0.11715374513859982"/> + </dataPoint> + <dataPoint> + <measurement value="77.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.117" errorPlus="0.4941831644238804" errorMinus="0.4941831644238804"/> + </dataPoint> + <dataPoint> + <measurement value="82.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.398" errorPlus="0.5850957186649036" errorMinus="0.5850957186649036"/> + </dataPoint> + <dataPoint> + <measurement value="87.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.00626" errorPlus="0.03811776488725434" errorMinus="0.03811776488725434"/> + </dataPoint> + <dataPoint> + <measurement value="92.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.0812" errorPlus="0.06215142798037709" errorMinus="0.06215142798037709"/> + </dataPoint> + <dataPoint> + <measurement value="97.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.00957" errorPlus="0.0099080018167136" errorMinus="0.0099080018167136"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d04-x01-y01" dimension="2" path="/REF/ATLAS_2012_CONF_2012_001" title="EVENTS" > + <dataPoint> + <measurement value="12.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="6.0" errorPlus="3.57" errorMinus="2.37"/> + </dataPoint> + <dataPoint> + <measurement value="17.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="6.0" errorPlus="3.57" errorMinus="2.37"/> + </dataPoint> + <dataPoint> + <measurement value="22.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="3.0" errorPlus="2.9" errorMinus="1.63"/> + </dataPoint> + <dataPoint> + <measurement value="27.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="3.0" errorPlus="2.9" errorMinus="1.63"/> + </dataPoint> + <dataPoint> + <measurement value="32.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="1.0" errorPlus="2.29" errorMinus="0.83"/> + </dataPoint> + <dataPoint> + <measurement value="37.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="5.0" errorPlus="3.37" errorMinus="2.15"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d04-x01-y02" dimension="2" path="/REF/ATLAS_2012_CONF_2012_001" title="EVENTS" > + <dataPoint> + <measurement value="12.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="5.64" errorPlus="1.7896647730790256" errorMinus="1.7896647730790256"/> + </dataPoint> + <dataPoint> + <measurement value="17.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="5.13" errorPlus="1.3675159962501353" errorMinus="1.3675159962501353"/> + </dataPoint> + <dataPoint> + <measurement value="22.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="4.05" errorPlus="0.8393449827097318" errorMinus="0.8393449827097318"/> + </dataPoint> + <dataPoint> + <measurement value="27.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="3.94" errorPlus="0.7849203781276163" errorMinus="0.7849203781276163"/> + </dataPoint> + <dataPoint> + <measurement value="32.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="2.48" errorPlus="0.4701063709417263" errorMinus="0.4701063709417263"/> + </dataPoint> + <dataPoint> + <measurement value="37.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="1.84" errorPlus="0.41182520563948" errorMinus="0.41182520563948"/> + </dataPoint> + <dataPoint> + <measurement value="42.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.521" errorPlus="0.14798986451781082" errorMinus="0.14798986451781082"/> + </dataPoint> + <dataPoint> + <measurement value="47.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.446" errorPlus="0.1380326048439281" errorMinus="0.1380326048439281"/> + </dataPoint> + <dataPoint> + <measurement value="52.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.164" errorPlus="0.08704596486914255" errorMinus="0.08704596486914255"/> + </dataPoint> + <dataPoint> + <measurement value="57.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.0134" errorPlus="0.04228853272460514" errorMinus="0.04228853272460514"/> + </dataPoint> + <dataPoint> + <measurement value="62.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.15" errorPlus="0.07951100552753687" errorMinus="0.07951100552753687"/> + </dataPoint> + <dataPoint> + <measurement value="67.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.00796" errorPlus="0.007420087600561061" errorMinus="0.007420087600561061"/> + </dataPoint> + <dataPoint> + <measurement value="72.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.0506" errorPlus="0.06631628759211419" errorMinus="0.06631628759211419"/> + </dataPoint> + <dataPoint> + <measurement value="77.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.00981" errorPlus="0.008057108662541421" errorMinus="0.008057108662541421"/> + </dataPoint> + <dataPoint> + <measurement value="87.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.00824" errorPlus="0.026284961479903293" errorMinus="0.026284961479903293"/> + </dataPoint> + <dataPoint> + <measurement value="97.5" errorPlus="2.5" errorMinus="2.5"/> + <measurement value="0.038" errorPlus="0.054166041022027815" errorMinus="0.054166041022027815"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d05-x01-y01" dimension="2" path="/REF/ATLAS_2012_CONF_2012_001" title="EVENTS" > + <dataPoint> + <measurement value="0.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="14.0" errorPlus="4.8" errorMinus="3.7"/> + </dataPoint> + <dataPoint> + <measurement value="1.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="8.0" errorPlus="3.93" errorMinus="2.76"/> + </dataPoint> + <dataPoint> + <measurement value="2.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="2.0" errorPlus="2.62" errorMinus="1.29"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d05-x01-y02" dimension="2" path="/REF/ATLAS_2012_CONF_2012_001" title="EVENTS" > + <dataPoint> + <measurement value="0.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="12.9" errorPlus="2.823118842698621" errorMinus="2.823118842698621"/> + </dataPoint> + <dataPoint> + <measurement value="1.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="6.84" errorPlus="1.6368567438844486" errorMinus="1.6368567438844486"/> + </dataPoint> + <dataPoint> + <measurement value="2.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="2.45" errorPlus="0.6711929677819933" errorMinus="0.6711929677819933"/> + </dataPoint> + <dataPoint> + <measurement value="3.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="0.956" errorPlus="0.6988569238406385" errorMinus="0.6988569238406385"/> + </dataPoint> + <dataPoint> + <measurement value="4.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="1.15" errorPlus="0.9475230867899738" errorMinus="0.9475230867899738"/> + </dataPoint> + <dataPoint> + <measurement value="5.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="0.197" errorPlus="0.1451206394693739" errorMinus="0.1451206394693739"/> + </dataPoint> + <dataPoint> + <measurement value="6.0" errorPlus="0.5" errorMinus="0.5"/> + <measurement value="0.0352" errorPlus="0.023754157530840782" errorMinus="0.023754157530840782"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d06-x01-y01" dimension="2" path="/REF/ATLAS_2012_CONF_2012_001" title="EVENTS" > + <dataPoint> + <measurement value="10.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="10.0" errorPlus="4.2" errorMinus="3.1"/> + </dataPoint> + <dataPoint> + <measurement value="30.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="9.0" errorPlus="4.09" errorMinus="2.93"/> + </dataPoint> + <dataPoint> + <measurement value="50.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="2.0" errorPlus="2.62" errorMinus="1.29"/> + </dataPoint> + <dataPoint> + <measurement value="70.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="1.0" errorPlus="2.29" errorMinus="0.83"/> + </dataPoint> + <dataPoint> + <measurement value="110.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="2.0" errorPlus="2.62" errorMinus="1.29"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d06-x01-y02" dimension="2" path="/REF/ATLAS_2012_CONF_2012_001" title="EVENTS" > + <dataPoint> + <measurement value="10.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="15.0" errorPlus="3.2893768406797057" errorMinus="3.2893768406797057"/> + </dataPoint> + <dataPoint> + <measurement value="30.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="7.39" errorPlus="2.3352944139872385" errorMinus="2.3352944139872385"/> + </dataPoint> + <dataPoint> + <measurement value="50.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="1.13" errorPlus="0.4326661530556787" errorMinus="0.4326661530556787"/> + </dataPoint> + <dataPoint> + <measurement value="70.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="0.301" errorPlus="0.15929846201391903" errorMinus="0.15929846201391903"/> + </dataPoint> + <dataPoint> + <measurement value="90.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="0.319" errorPlus="0.1345288073239334" errorMinus="0.1345288073239334"/> + </dataPoint> + <dataPoint> + <measurement value="110.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="0.138" errorPlus="0.06226556030423239" errorMinus="0.06226556030423239"/> + </dataPoint> + <dataPoint> + <measurement value="130.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="0.0517" errorPlus="0.05756735185849702" errorMinus="0.05756735185849702"/> + </dataPoint> + <dataPoint> + <measurement value="150.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="0.0689" errorPlus="0.06552312568856891" errorMinus="0.06552312568856891"/> + </dataPoint> + <dataPoint> + <measurement value="170.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="0.0391" errorPlus="0.023447174669882935" errorMinus="0.023447174669882935"/> + </dataPoint> + <dataPoint> + <measurement value="190.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="0.00278" errorPlus="0.007254157428674953" errorMinus="0.007254157428674953"/> + </dataPoint> + <dataPoint> + <measurement value="210.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="0.0117" errorPlus="0.010770329614269008" errorMinus="0.010770329614269008"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d07-x01-y01" dimension="2" path="/REF/ATLAS_2012_CONF_2012_001" title="EVENTS" > + <dataPoint> + <measurement value="50.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="2.0" errorPlus="2.62" errorMinus="1.29"/> + </dataPoint> + <dataPoint> + <measurement value="70.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="2.0" errorPlus="2.62" errorMinus="1.29"/> + </dataPoint> + <dataPoint> + <measurement value="90.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="20.0" errorPlus="5.5" errorMinus="4.4"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d07-x01-y02" dimension="2" path="/REF/ATLAS_2012_CONF_2012_001" title="EVENTS" > + <dataPoint> + <measurement value="30.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="0.633" errorPlus="0.5619368291899011" errorMinus="0.5619368291899011"/> + </dataPoint> + <dataPoint> + <measurement value="50.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="0.917" errorPlus="0.6867059050277637" errorMinus="0.6867059050277637"/> + </dataPoint> + <dataPoint> + <measurement value="70.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="0.919" errorPlus="0.4574330552113609" errorMinus="0.4574330552113609"/> + </dataPoint> + <dataPoint> + <measurement value="90.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="21.6" errorPlus="3.8275318418009276" errorMinus="3.8275318418009276"/> + </dataPoint> + <dataPoint> + <measurement value="110.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="0.308" errorPlus="0.12165936051122413" errorMinus="0.12165936051122413"/> + </dataPoint> + <dataPoint> + <measurement value="130.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="0.0382" errorPlus="0.032252286740632824" errorMinus="0.032252286740632824"/> + </dataPoint> + <dataPoint> + <measurement value="150.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="2.52E-4" errorPlus="1.7273389939441536E-4" errorMinus="1.7273389939441536E-4"/> + </dataPoint> + <dataPoint> + <measurement value="170.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="0.0159" errorPlus="0.01596558799418299" errorMinus="0.01596558799418299"/> + </dataPoint> + <dataPoint> + <measurement value="190.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="0.0023" errorPlus="0.0024781444671366517" errorMinus="0.0024781444671366517"/> + </dataPoint> + <dataPoint> + <measurement value="230.0" errorPlus="10.0" errorMinus="10.0"/> + <measurement value="0.00558" errorPlus="0.006024616170346456" errorMinus="0.006024616170346456"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d08-x01-y01" dimension="2" path="/REF/ATLAS_2012_CONF_2012_001" title="EVENTS" > + <dataPoint> + <measurement value="125.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="4.0" errorPlus="3.15" errorMinus="1.91"/> + </dataPoint> + <dataPoint> + <measurement value="175.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="8.0" errorPlus="3.93" errorMinus="2.76"/> + </dataPoint> + <dataPoint> + <measurement value="225.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="6.0" errorPlus="3.57" errorMinus="2.37"/> + </dataPoint> + <dataPoint> + <measurement value="275.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="1.0" errorPlus="2.29" errorMinus="0.83"/> + </dataPoint> + <dataPoint> + <measurement value="325.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="2.0" errorPlus="2.62" errorMinus="1.29"/> + </dataPoint> + <dataPoint> + <measurement value="375.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="2.0" errorPlus="2.62" errorMinus="1.29"/> + </dataPoint> + <dataPoint> + <measurement value="525.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="1.0" errorPlus="2.29" errorMinus="0.83"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d08-x01-y02" dimension="2" path="/REF/ATLAS_2012_CONF_2012_001" title="EVENTS" > + <dataPoint> + <measurement value="75.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="0.617" errorPlus="0.5625735507469224" errorMinus="0.5625735507469224"/> + </dataPoint> + <dataPoint> + <measurement value="125.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="3.27" errorPlus="1.1718361660232202" errorMinus="1.1718361660232202"/> + </dataPoint> + <dataPoint> + <measurement value="175.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="7.46" errorPlus="1.4560219778561039" errorMinus="1.4560219778561039"/> + </dataPoint> + <dataPoint> + <measurement value="225.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="5.69" errorPlus="1.0911003620199198" errorMinus="1.0911003620199198"/> + </dataPoint> + <dataPoint> + <measurement value="275.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="2.63" errorPlus="0.5546169849544819" errorMinus="0.5546169849544819"/> + </dataPoint> + <dataPoint> + <measurement value="325.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="2.05" errorPlus="0.7433034373659253" errorMinus="0.7433034373659253"/> + </dataPoint> + <dataPoint> + <measurement value="375.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="1.17" errorPlus="0.41182520563948" errorMinus="0.41182520563948"/> + </dataPoint> + <dataPoint> + <measurement value="425.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="0.506" errorPlus="0.1708010538609174" errorMinus="0.1708010538609174"/> + </dataPoint> + <dataPoint> + <measurement value="475.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="0.155" errorPlus="0.08905054744357274" errorMinus="0.08905054744357274"/> + </dataPoint> + <dataPoint> + <measurement value="525.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="0.187" errorPlus="0.08508818954473059" errorMinus="0.08508818954473059"/> + </dataPoint> + <dataPoint> + <measurement value="575.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="0.42" errorPlus="0.34988569562072697" errorMinus="0.34988569562072697"/> + </dataPoint> + <dataPoint> + <measurement value="625.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="0.12" errorPlus="0.0604400529450463" errorMinus="0.0604400529450463"/> + </dataPoint> + <dataPoint> + <measurement value="675.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="0.0901" errorPlus="0.05932663819904175" errorMinus="0.05932663819904175"/> + </dataPoint> + <dataPoint> + <measurement value="725.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="0.0391" errorPlus="0.036725059564281175" errorMinus="0.036725059564281175"/> + </dataPoint> + <dataPoint> + <measurement value="775.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="0.0172" errorPlus="0.015778466338652816" errorMinus="0.015778466338652816"/> + </dataPoint> + <dataPoint> + <measurement value="825.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="0.0217" errorPlus="0.014204224723651763" errorMinus="0.014204224723651763"/> + </dataPoint> + <dataPoint> + <measurement value="875.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="0.0231" errorPlus="0.016559287424282483" errorMinus="0.016559287424282483"/> + </dataPoint> + <dataPoint> + <measurement value="925.0" errorPlus="25.0" errorMinus="25.0"/> + <measurement value="0.00607" errorPlus="0.007376211764855995" errorMinus="0.007376211764855995"/> + </dataPoint> + </dataPointSet> +</aida> Modified: trunk/data/refdata/Makefile.am ============================================================================== --- trunk/data/refdata/Makefile.am Mon Apr 16 13:10:22 2012 (r3684) +++ trunk/data/refdata/Makefile.am Mon Apr 16 13:14:25 2012 (r3685) @@ -44,6 +44,7 @@ ATLAS_2012_I1083318.aida \ ATLAS_2012_I1084540.aida \ ATLAS_2012_I1094568.aida \ + ATLAS_2012_CONF_2012_001.aida \ BABAR_2005_S6181155.aida \ BABAR_2006_S6511112.aida \ BABAR_2007_S6895344.aida \ Copied and modified: trunk/src/Analyses/ATLAS_2012_CONF_2012_001.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_001.cc Mon Apr 16 13:14:25 2012 (r3685) @@ -13,16 +13,16 @@ namespace Rivet { - - class ATLAS_2011_S9225137 : public Analysis { + /// @author Peter Richardson + class ATLAS_2012_CONF_2012_001 : public Analysis { public: /// @name Constructors etc. //@{ /// Constructor - ATLAS_2011_S9225137() - : Analysis("ATLAS_2011_S9225137") + ATLAS_2012_CONF_2012_001() + : Analysis("ATLAS_2012_CONF_2012_001") { } //@} @@ -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,28 @@ // 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.push_back(bookHistogram1D(1,1,1)); + _hist_leptonpT.push_back(bookHistogram1D(2,1,1)); + _hist_leptonpT.push_back(bookHistogram1D(3,1,1)); + _hist_leptonpT.push_back(bookHistogram1D(4,1,1)); + _hist_njet = bookHistogram1D(5,1,1); + _hist_etmiss = bookHistogram1D(6,1,1); + _hist_mSFOS = bookHistogram1D(7,1,1); + _hist_meff = bookHistogram1D(8,1,1); + _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,40 +105,56 @@ } // 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 + 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; + 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 ); } @@ -196,12 +163,21 @@ 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 +191,139 @@ } 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 less than 4 electrons and muons + if ( recon_mu.size() + recon_e.size() < 4 ) { + 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); - } - 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); + // 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) { - _njet80B_obs->fill(njet80,weight); - _njet80B_bac->fill(njet80,weight); - _njet80B_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(njet55>3) { - _njet55C_obs->fill(njet55,weight); - _njet55C_bac->fill(njet55,weight); - _njet55C_sig->fill(njet55,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; + } + + // calculate meff + double meff = eTmiss; + foreach ( const Particle & e , recon_e ) + meff += e.momentum().perp(); + foreach ( const Particle & mu, recon_mu ) + meff += mu.momentum().perp(); + foreach ( const Jet & jet, recon_jets ) { + double pT = jet.momentum().perp(); + if(pT>40.) meff += pT; + } + + double mSFOS=1e30, mdiff=1e30; + // mass of SFOS pairs closest to the Z mass + 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; + double mtest = (recon_e[ix].momentum()+recon_e[iy].momentum()).mass(); + 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; + double mtest = (recon_mu[ix].momentum()+recon_mu[iy].momentum()).mass(); + if(fabs(mtest-90.)<mdiff) { + mSFOS = mtest; + mdiff = fabs(mtest-90.); + } + } + } + + // make the control plots + // lepton pT + unsigned int ie=0,imu=0; + for(unsigned int ix=0;ix<4;++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_leptonpT[ix]->fill(pTe ,weight); + ++ie; } - if(njet80>3) { - _njet80C_obs->fill(njet80,weight); - _njet80C_bac->fill(njet80,weight); - _njet80C_sig->fill(njet80,weight); + else { + _hist_leptonpT[ix]->fill(pTmu,weight); + ++imu; } } - - // apply E_T/sqrt(H_T) cut - if(etht<=3.5*GeV) { - MSG_DEBUG("Fails ET/sqrt(HT) cut "); - vetoEvent; + // njet + _hist_njet->fill(recon_jets.size(),weight); + // etmiss + _hist_etmiss->fill(eTmiss,weight); + if(mSFOS<1e30) _hist_mSFOS->fill(mSFOS,weight); + _hist_meff->fill(meff,weight); + + // finally the counts + if(eTmiss>50.) { + _count_SR1->fill(0.5,weight); + if(mdiff<10.) _count_SR2->fill(0.5,weight); } - - // 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 10 GeV + scale(_hist_leptonpT[0],norm*10.); + scale(_hist_leptonpT[1],norm*10.); + // these are number of events at 2.06fb^-1 per 5 GeV + scale(_hist_leptonpT[2],norm*5.); + scale(_hist_leptonpT[3],norm*5.); + // these are number of events at 2.06fb^-1 per 20 GeV + scale(_hist_etmiss ,norm*20.); + scale(_hist_mSFOS ,norm*20.); + // these are number of events at 2.06fb^-1 per 50 GeV + scale(_hist_meff ,norm*50.); + // these are number of events at 2.06fb^-1 + scale(_hist_njet ,norm); + 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; + vector<AIDA::IHistogram1D*> _hist_leptonpT; + AIDA::IHistogram1D* _hist_njet; + AIDA::IHistogram1D* _hist_etmiss; + AIDA::IHistogram1D* _hist_mSFOS; + AIDA::IHistogram1D* _hist_meff; + 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_001); } Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Mon Apr 16 13:10:22 2012 (r3684) +++ trunk/src/Analyses/Makefile.am Mon Apr 16 13:14:25 2012 (r3685) @@ -78,6 +78,7 @@ ATLAS_2011_S9225137.cc \ ATLAS_2012_I1083318.cc \ ATLAS_2012_I1095236.cc \ + ATLAS_2012_CONF_2012_001.cc \ ATLAS_2012_CONF_2012_033.cc \ ATLAS_2012_CONF_2012_037.cc endif
More information about the Rivet-svn mailing list |