|
|
[Rivet-svn] r3754 - in branches/2011-07-aida2yoda: data/anainfo data/plotinfo data/refdata src/Analyses
blackhole at projects.hepforge.org
blackhole at projects.hepforge.org
Fri Jun 8 15:17:44 BST 2012
Author: hoeth
Date: Fri Jun 8 15:17:44 2012
New Revision: 3754
Log:
merge c3716, 3717, 3742, 3753 from trunk
Added:
branches/2011-07-aida2yoda/data/anainfo/ATLAS_2012_I1093738.info
- copied unchanged from r3716, trunk/data/anainfo/ATLAS_2012_I1093738.info
branches/2011-07-aida2yoda/data/plotinfo/ATLAS_2012_I1093738.plot
- copied unchanged from r3717, trunk/data/plotinfo/ATLAS_2012_I1093738.plot
branches/2011-07-aida2yoda/data/refdata/ATLAS_2012_I1093738.aida
- copied unchanged from r3716, trunk/data/refdata/ATLAS_2012_I1093738.aida
branches/2011-07-aida2yoda/src/Analyses/ATLAS_2012_I1093738.cc
- copied, changed from r3716, trunk/src/Analyses/ATLAS_2012_I1093738.cc
Modified:
branches/2011-07-aida2yoda/data/anainfo/Makefile.am
branches/2011-07-aida2yoda/data/plotinfo/Makefile.am
branches/2011-07-aida2yoda/data/refdata/Makefile.am
branches/2011-07-aida2yoda/src/Analyses/Makefile.am
Copied: branches/2011-07-aida2yoda/data/anainfo/ATLAS_2012_I1093738.info (from r3716, trunk/data/anainfo/ATLAS_2012_I1093738.info)
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ branches/2011-07-aida2yoda/data/anainfo/ATLAS_2012_I1093738.info Fri Jun 8 15:17:44 2012 (r3754, copy of r3716, trunk/data/anainfo/ATLAS_2012_I1093738.info)
@@ -0,0 +1,41 @@
+Name: ATLAS_2012_I1093738
+Year: 2012
+Summary: Isolated prompt photon + jet xsection
+Experiment: ATLAS
+Collider: LHC
+SpiresID: 1093738
+Status: VALIDATED
+Authors:
+ - Giovanni Marchiori <giovanni.marchiori at cern.ch>
+References:
+ - arXiv:1203.3161 [hep-ex]
+RunInfo:
+ Inclusive photon+jet+X events at $\sqrt{s} = 7$~TeV.
+NumEvents: 1000000
+Beams: [p+, p+]
+Energies: [7000]
+Description:
+ 'A measurement of the production cross section for isolated photons in association with jets
+ in $pp$ collisions at $sqrt{s} = 7$ TeV.
+ Photons with $|\eta|<1.37$ and $E_T>25$~GeV and jets with $|y|<4.4$ and $p_T>20$ GeV
+ are selected.
+ The differential cross section as a function of the photon transverse energy is measured, for
+ three leading jet rapidity configurations, separately for the cases where the photon and
+ jet rapidities have the same or the opposite sign.
+ The measurement uses 37 pb$^{-1}$ of integrated luminosity collected with the ATLAS detector.'
+BibKey: Aad:2012cv
+BibTeX:
+'@article{Aad:2012cv,
+ author = "Aad, Georges and others",
+ title = "{Measurement of the production cross section of an
+ isolated photon associated with jets in proton-proton
+ collisions at sqrt(s) = 7 TeV with the ATLAS detector}",
+ collaboration = "ATLAS",
+ year = "2012",
+ eprint = "1203.3161",
+ archivePrefix = "arXiv",
+ primaryClass = "hep-ex",
+ reportNumber = "CERN-PH-EP-2012-009",
+ SLACcitation = "%%CITATION = ARXIV:1203.3161;%%",
+}'
+
Modified: branches/2011-07-aida2yoda/data/anainfo/Makefile.am
==============================================================================
--- branches/2011-07-aida2yoda/data/anainfo/Makefile.am Fri Jun 8 15:15:31 2012 (r3753)
+++ branches/2011-07-aida2yoda/data/anainfo/Makefile.am Fri Jun 8 15:17:44 2012 (r3754)
@@ -44,6 +44,7 @@
ATLAS_2011_S9212353.info \
ATLAS_2012_I1083318.info \
ATLAS_2012_I1084540.info \
+ ATLAS_2012_I1093738.info \
ATLAS_2012_I1094568.info \
ATLAS_2012_I1095236.info \
ATLAS_2012_CONF_2012_001.info \
Copied: branches/2011-07-aida2yoda/data/plotinfo/ATLAS_2012_I1093738.plot (from r3717, trunk/data/plotinfo/ATLAS_2012_I1093738.plot)
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ branches/2011-07-aida2yoda/data/plotinfo/ATLAS_2012_I1093738.plot Fri Jun 8 15:17:44 2012 (r3754, copy of r3717, trunk/data/plotinfo/ATLAS_2012_I1093738.plot)
@@ -0,0 +1,54 @@
+# BEGIN PLOT /ATLAS_2012_I1093738/d01-x01-y01
+Title=Leading photon $E_\perp$ (central jets, same-sign rapidity)
+XLabel=$E_\perp^\gamma$ [GeV]
+YLabel=$\mathrm{d}\sigma/\mathrm{d}E_\perp^\gamma$ [pb/GeV]
+LogX=1
+XMin=25.0
+XMax=400.0
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2012_I1093738/d02-x01-y01
+Title=Leading photon $E_\perp$ (forward jets, same-sign rapidity)
+XLabel=$E_\perp^\gamma$ [GeV]
+YLabel=$\mathrm{d}\sigma/\mathrm{d}E_\perp^\gamma$ [pb/GeV]
+LogX=1
+XMin=25.0
+XMax=400.0
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2012_I1093738/d03-x01-y01
+Title=Leading photon $E_\perp$ (very forward jets, same-sign rapidity)
+XLabel=$E_\perp^\gamma$ [GeV]
+YLabel=$\mathrm{d}\sigma/\mathrm{d}E_\perp^\gamma$ [pb/GeV]
+LogX=1
+XMin=25.0
+XMax=200.0
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2012_I1093738/d04-x01-y01
+Title=Leading photon $E_\perp$ (central jets, opposite-sign rapidity)
+XLabel=$E_\perp^\gamma$ [GeV]
+YLabel=$\mathrm{d}\sigma/\mathrm{d}E_\perp^\gamma$ [pb/GeV]
+LogX=1
+XMin=25.0
+XMax=400.0
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2012_I1093738/d05-x01-y01
+Title=Leading photon $E_\perp$ (forward jets, opposite-sign rapidity)
+XLabel=$E_\perp^\gamma$ [GeV]
+YLabel=$\mathrm{d}\sigma/\mathrm{d}E_\perp^\gamma$ [pb/GeV]
+LogX=1
+XMin=25.0
+XMax=400.0
+# END PLOT
+
+# BEGIN PLOT /ATLAS_2012_I1093738/d06-x01-y01
+Title=Leading photon $E_\perp$ (very forward jets, opposite-sign rapidity)
+XLabel=$E_\perp^\gamma$ [GeV]
+YLabel=$\mathrm{d}\sigma/\mathrm{d}E_\perp^\gamma$ [pb/GeV]
+LogX=1
+XMin=25.0
+XMax=200.0
+# END PLOT
+
Modified: branches/2011-07-aida2yoda/data/plotinfo/Makefile.am
==============================================================================
--- branches/2011-07-aida2yoda/data/plotinfo/Makefile.am Fri Jun 8 15:15:31 2012 (r3753)
+++ branches/2011-07-aida2yoda/data/plotinfo/Makefile.am Fri Jun 8 15:17:44 2012 (r3754)
@@ -41,6 +41,7 @@
ATLAS_2011_S9212353.plot \
ATLAS_2012_I1083318.plot \
ATLAS_2012_I1084540.plot \
+ ATLAS_2012_I1093738.plot \
ATLAS_2012_I1094568.plot \
ATLAS_2012_I1095236.plot \
ATLAS_2012_CONF_2012_001.plot \
Copied: branches/2011-07-aida2yoda/data/refdata/ATLAS_2012_I1093738.aida (from r3716, trunk/data/refdata/ATLAS_2012_I1093738.aida)
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ branches/2011-07-aida2yoda/data/refdata/ATLAS_2012_I1093738.aida Fri Jun 8 15:17:44 2012 (r3754, copy of r3716, trunk/data/refdata/ATLAS_2012_I1093738.aida)
@@ -0,0 +1,315 @@
+<?xml version="1.0" ?>
+<!DOCTYPE aida SYSTEM "http://aida.freehep.org/schemas/3.3/aida.dtd">
+<aida version="3.3">
+ <implementation version="1.1" package="Rivet"/>
+ <dataPointSet name="d01-x01-y01" dimension="2"
+ path="/REF/ATLAS_2012_I1093738" title="Leading photon $p_\perp$ (central jets, same-sign rapidity)">
+ <dimension dim="0" title="$p_\perp^\gamma$ [GeV]"/>
+ <dimension dim="1" title="$\mathrm{d}\sigma/\mathrm{d}p_\perp^\gamma$ [pb/GeV]"/>
+ <annotation>
+ <item key="Title" value="Leading photon $p_\perp$ (central jets, same-sign rapidity)" sticky="true"/>
+ <item key="AidaPath" value="/REF/ATLAS_2012_I1093738/d01-x01-y01" sticky="true"/>
+ </annotation>
+ <dataPoint>
+ <measurement errorPlus="2.500000e+00" value="2.750000e+01" errorMinus="2.500000e+00"/>
+ <measurement errorPlus="3.100000e+01" value="3.940000e+02" errorMinus="7.400000e+01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="2.500000e+00" value="3.250000e+01" errorMinus="2.500000e+00"/>
+ <measurement errorPlus="2.300000e+01" value="2.580000e+02" errorMinus="5.000000e+01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="5.000000e+00" value="4.000000e+01" errorMinus="5.000000e+00"/>
+ <measurement errorPlus="1.300000e+01" value="1.370000e+02" errorMinus="2.700000e+01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="5.000000e+00" value="5.000000e+01" errorMinus="5.000000e+00"/>
+ <measurement errorPlus="5.200000e+00" value="6.090000e+01" errorMinus="7.100000e+00"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="7.500000e+00" value="6.250000e+01" errorMinus="7.500000e+00"/>
+ <measurement errorPlus="2.400000e+00" value="2.480000e+01" errorMinus="3.000000e+00"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="7.500000e+00" value="7.750000e+01" errorMinus="7.500000e+00"/>
+ <measurement errorPlus="1.000000e+00" value="9.510000e+00" errorMinus="1.240000e+00"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="7.500000e+00" value="9.250000e+01" errorMinus="7.500000e+00"/>
+ <measurement errorPlus="5.000000e-01" value="4.400000e+00" errorMinus="5.700000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="1.250000e+01" value="1.125000e+02" errorMinus="1.250000e+01"/>
+ <measurement errorPlus="2.200000e-01" value="1.770000e+00" errorMinus="2.400000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="1.250000e+01" value="1.375000e+02" errorMinus="1.250000e+01"/>
+ <measurement errorPlus="9.300000e-02" value="6.980000e-01" errorMinus="1.030000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="2.500000e+01" value="1.750000e+02" errorMinus="2.500000e+01"/>
+ <measurement errorPlus="3.400000e-02" value="2.260000e-01" errorMinus="3.600000e-02"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="1.000000e+02" value="3.000000e+02" errorMinus="1.000000e+02"/>
+ <measurement errorPlus="4.800000e-03" value="2.830000e-02" errorMinus="5.000000e-03"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet name="d02-x01-y01" dimension="2"
+ path="/REF/ATLAS_2012_I1093738" title="Leading photon $p_\perp$ (forward jets, same-sign rapidity)">
+ <dimension dim="0" title="$p_\perp^\gamma$ [GeV]"/>
+ <dimension dim="1" title="$\mathrm{d}\sigma/\mathrm{d}p_\perp^\gamma$ [pb/GeV]"/>
+ <annotation>
+ <item key="Title" value="Leading photon $p_\perp$ (forward jets, same-sign rapidity)" sticky="true"/>
+ <item key="AidaPath" value="/REF/ATLAS_2012_I1093738/d02-x01-y01" sticky="true"/>
+ </annotation>
+ <dataPoint>
+ <measurement errorPlus="2.500000e+00" value="2.750000e+01" errorMinus="2.500000e+00"/>
+ <measurement errorPlus="3.100000e+01" value="3.160000e+02" errorMinus="5.500000e+01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="2.500000e+00" value="3.250000e+01" errorMinus="2.500000e+00"/>
+ <measurement errorPlus="2.300000e+01" value="2.100000e+02" errorMinus="3.700000e+01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="5.000000e+00" value="4.000000e+01" errorMinus="5.000000e+00"/>
+ <measurement errorPlus="1.200000e+01" value="1.050000e+02" errorMinus="1.900000e+01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="5.000000e+00" value="5.000000e+01" errorMinus="5.000000e+00"/>
+ <measurement errorPlus="3.800000e+00" value="4.360000e+01" errorMinus="5.100000e+00"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="7.500000e+00" value="6.250000e+01" errorMinus="7.500000e+00"/>
+ <measurement errorPlus="1.700000e+00" value="1.750000e+01" errorMinus="2.200000e+00"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="7.500000e+00" value="7.750000e+01" errorMinus="7.500000e+00"/>
+ <measurement errorPlus="6.800000e-01" value="6.390000e+00" errorMinus="8.400000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="7.500000e+00" value="9.250000e+01" errorMinus="7.500000e+00"/>
+ <measurement errorPlus="3.100000e-01" value="2.710000e+00" errorMinus="3.600000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="1.250000e+01" value="1.125000e+02" errorMinus="1.250000e+01"/>
+ <measurement errorPlus="1.600000e-01" value="1.270000e+00" errorMinus="1.800000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="1.250000e+01" value="1.375000e+02" errorMinus="1.250000e+01"/>
+ <measurement errorPlus="6.100000e-02" value="4.410000e-01" errorMinus="6.800000e-02"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="2.500000e+01" value="1.750000e+02" errorMinus="2.500000e+01"/>
+ <measurement errorPlus="1.800000e-02" value="1.020000e-01" errorMinus="1.900000e-02"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="1.000000e+02" value="3.000000e+02" errorMinus="1.000000e+02"/>
+ <measurement errorPlus="2.100000e-03" value="9.000000e-03" errorMinus="2.200000e-03"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet name="d03-x01-y01" dimension="2"
+ path="/REF/ATLAS_2012_I1093738" title="Leading photon $p_\perp$ (very forward jets, same-sign rapidity)">
+ <dimension dim="0" title="$p_\perp^\gamma$ [GeV]"/>
+ <dimension dim="1" title="$\mathrm{d}\sigma/\mathrm{d}p_\perp^\gamma$ [pb/GeV]"/>
+ <annotation>
+ <item key="Title" value="Leading photon $p_\perp$ (very forward jets, same-sign rapidity)" sticky="true"/>
+ <item key="AidaPath" value="/REF/ATLAS_2012_I1093738/d03-x01-y01" sticky="true"/>
+ </annotation>
+ <dataPoint>
+ <measurement errorPlus="2.500000e+00" value="2.750000e+01" errorMinus="2.500000e+00"/>
+ <measurement errorPlus="1.000000e+01" value="6.600000e+01" errorMinus="1.900000e+01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="2.500000e+00" value="3.250000e+01" errorMinus="2.500000e+00"/>
+ <measurement errorPlus="8.000000e+00" value="4.600000e+01" errorMinus="1.300000e+01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="5.000000e+00" value="4.000000e+01" errorMinus="5.000000e+00"/>
+ <measurement errorPlus="3.000000e+00" value="2.000000e+01" errorMinus="6.000000e+00"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="5.000000e+00" value="5.000000e+01" errorMinus="5.000000e+00"/>
+ <measurement errorPlus="8.000000e-01" value="8.100000e+00" errorMinus="1.400000e+00"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="7.500000e+00" value="6.250000e+01" errorMinus="7.500000e+00"/>
+ <measurement errorPlus="3.000000e-01" value="2.400000e+00" errorMinus="4.000000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="7.500000e+00" value="7.750000e+01" errorMinus="7.500000e+00"/>
+ <measurement errorPlus="1.100000e-01" value="8.600000e-01" errorMinus="1.700000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="7.500000e+00" value="9.250000e+01" errorMinus="7.500000e+00"/>
+ <measurement errorPlus="4.000000e-02" value="2.400000e-01" errorMinus="4.000000e-02"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="1.250000e+01" value="1.125000e+02" errorMinus="1.250000e+01"/>
+ <measurement errorPlus="2.000000e-02" value="7.000000e-02" errorMinus="2.000000e-02"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="1.250000e+01" value="1.375000e+02" errorMinus="1.250000e+01"/>
+ <measurement errorPlus="7.000000e-03" value="1.400000e-02" errorMinus="7.000000e-03"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="2.500000e+01" value="1.750000e+02" errorMinus="2.500000e+01"/>
+ <measurement errorPlus="2.000000e-03" value="2.800000e-03" errorMinus="1.900000e-03"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet name="d04-x01-y01" dimension="2"
+ path="/REF/ATLAS_2012_I1093738" title="Leading photon $p_\perp$ (central jets, opposite-sign rapidity)">
+ <dimension dim="0" title="$p_\perp^\gamma$ [GeV]"/>
+ <dimension dim="1" title="$\mathrm{d}\sigma/\mathrm{d}p_\perp^\gamma$ [pb/GeV]"/>
+ <annotation>
+ <item key="Title" value="Leading photon $p_\perp$ (central jets, opposite-sign rapidity)" sticky="true"/>
+ <item key="AidaPath" value="/REF/ATLAS_2012_I1093738/d04-x01-y01" sticky="true"/>
+ </annotation>
+ <dataPoint>
+ <measurement errorPlus="2.500000e+00" value="2.750000e+01" errorMinus="2.500000e+00"/>
+ <measurement errorPlus="3.000000e+01" value="3.240000e+02" errorMinus="6.500000e+01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="2.500000e+00" value="3.250000e+01" errorMinus="2.500000e+00"/>
+ <measurement errorPlus="2.000000e+01" value="2.010000e+02" errorMinus="4.100000e+01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="5.000000e+00" value="4.000000e+01" errorMinus="5.000000e+00"/>
+ <measurement errorPlus="1.200000e+01" value="1.120000e+02" errorMinus="2.300000e+01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="5.000000e+00" value="5.000000e+01" errorMinus="5.000000e+00"/>
+ <measurement errorPlus="3.900000e+00" value="4.550000e+01" errorMinus="5.600000e+00"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="7.500000e+00" value="6.250000e+01" errorMinus="7.500000e+00"/>
+ <measurement errorPlus="1.800000e+00" value="1.830000e+01" errorMinus="2.400000e+00"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="7.500000e+00" value="7.750000e+01" errorMinus="7.500000e+00"/>
+ <measurement errorPlus="7.600000e-01" value="7.180000e+00" errorMinus="9.900000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="7.500000e+00" value="9.250000e+01" errorMinus="7.500000e+00"/>
+ <measurement errorPlus="3.800000e-01" value="3.260000e+00" errorMinus="4.000000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="1.250000e+01" value="1.125000e+02" errorMinus="1.250000e+01"/>
+ <measurement errorPlus="1.700000e-01" value="1.360000e+00" errorMinus="1.700000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="1.250000e+01" value="1.375000e+02" errorMinus="1.250000e+01"/>
+ <measurement errorPlus="7.200000e-02" value="5.030000e-01" errorMinus="7.500000e-02"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="2.500000e+01" value="1.750000e+02" errorMinus="2.500000e+01"/>
+ <measurement errorPlus="2.500000e-02" value="1.560000e-01" errorMinus="2.700000e-02"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="1.000000e+02" value="3.000000e+02" errorMinus="1.000000e+02"/>
+ <measurement errorPlus="3.300000e-03" value="1.820000e-02" errorMinus="3.500000e-03"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet name="d05-x01-y01" dimension="2"
+ path="/REF/ATLAS_2012_I1093738" title="Leading photon $p_\perp$ (forward jets, opposite-sign rapidity)">
+ <dimension dim="0" title="$p_\perp^\gamma$ [GeV]"/>
+ <dimension dim="1" title="$\mathrm{d}\sigma/\mathrm{d}p_\perp^\gamma$ [pb/GeV]"/>
+ <annotation>
+ <item key="Title" value="Leading photon $p_\perp$ (forward jets, opposite-sign rapidity)" sticky="true"/>
+ <item key="AidaPath" value="/REF/ATLAS_2012_I1093738/d05-x01-y01" sticky="true"/>
+ </annotation>
+ <dataPoint>
+ <measurement errorPlus="2.500000e+00" value="2.750000e+01" errorMinus="2.500000e+00"/>
+ <measurement errorPlus="2.700000e+01" value="1.880000e+02" errorMinus="3.600000e+01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="2.500000e+00" value="3.250000e+01" errorMinus="2.500000e+00"/>
+ <measurement errorPlus="1.800000e+01" value="1.150000e+02" errorMinus="2.300000e+01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="5.000000e+00" value="4.000000e+01" errorMinus="5.000000e+00"/>
+ <measurement errorPlus="9.000000e+00" value="5.800000e+01" errorMinus="1.200000e+01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="5.000000e+00" value="5.000000e+01" errorMinus="5.000000e+00"/>
+ <measurement errorPlus="2.200000e+00" value="2.410000e+01" errorMinus="3.100000e+00"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="7.500000e+00" value="6.250000e+01" errorMinus="7.500000e+00"/>
+ <measurement errorPlus="9.000000e-01" value="8.800000e+00" errorMinus="1.200000e+00"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="7.500000e+00" value="7.750000e+01" errorMinus="7.500000e+00"/>
+ <measurement errorPlus="3.700000e-01" value="3.320000e+00" errorMinus="4.800000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="7.500000e+00" value="9.250000e+01" errorMinus="7.500000e+00"/>
+ <measurement errorPlus="1.800000e-01" value="1.490000e+00" errorMinus="1.900000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="1.250000e+01" value="1.125000e+02" errorMinus="1.250000e+01"/>
+ <measurement errorPlus="7.000000e-02" value="5.400000e-01" errorMinus="7.000000e-02"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="1.250000e+01" value="1.375000e+02" errorMinus="1.250000e+01"/>
+ <measurement errorPlus="3.100000e-02" value="1.750000e-01" errorMinus="3.100000e-02"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="2.500000e+01" value="1.750000e+02" errorMinus="2.500000e+01"/>
+ <measurement errorPlus="1.100000e-02" value="5.500000e-02" errorMinus="1.100000e-02"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="1.000000e+02" value="3.000000e+02" errorMinus="1.000000e+02"/>
+ <measurement errorPlus="1.200000e-03" value="4.100000e-03" errorMinus="1.100000e-03"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet name="d06-x01-y01" dimension="2"
+ path="/REF/ATLAS_2012_I1093738" title="Leading photon $p_\perp$ (very forward jets, opposite-sign rapidity)">
+ <dimension dim="0" title="$p_\perp^\gamma$ [GeV]"/>
+ <dimension dim="1" title="$\mathrm{d}\sigma/\mathrm{d}p_\perp^\gamma$ [pb/GeV]"/>
+ <annotation>
+ <item key="Title" value="Leading photon $p_\perp$ (very forward jets, opposite-sign rapidity)" sticky="true"/>
+ <item key="AidaPath" value="/REF/ATLAS_2012_I1093738/d06-x01-y01" sticky="true"/>
+ </annotation>
+ <dataPoint>
+ <measurement errorPlus="2.500000e+00" value="2.750000e+01" errorMinus="2.500000e+00"/>
+ <measurement errorPlus="6.000000e+00" value="3.100000e+01" errorMinus="1.300000e+01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="2.500000e+00" value="3.250000e+01" errorMinus="2.500000e+00"/>
+ <measurement errorPlus="4.000000e+00" value="2.100000e+01" errorMinus="9.000000e+00"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="5.000000e+00" value="4.000000e+01" errorMinus="5.000000e+00"/>
+ <measurement errorPlus="2.000000e+00" value="1.200000e+01" errorMinus="5.000000e+00"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="5.000000e+00" value="5.000000e+01" errorMinus="5.000000e+00"/>
+ <measurement errorPlus="6.000000e-01" value="3.500000e+00" errorMinus="1.100000e+00"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="7.500000e+00" value="6.250000e+01" errorMinus="7.500000e+00"/>
+ <measurement errorPlus="3.000000e-01" value="1.500000e+00" errorMinus="5.000000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="7.500000e+00" value="7.750000e+01" errorMinus="7.500000e+00"/>
+ <measurement errorPlus="8.000000e-02" value="3.800000e-01" errorMinus="1.200000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="7.500000e+00" value="9.250000e+01" errorMinus="7.500000e+00"/>
+ <measurement errorPlus="3.000000e-02" value="1.200000e-01" errorMinus="3.000000e-02"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="1.250000e+01" value="1.125000e+02" errorMinus="1.250000e+01"/>
+ <measurement errorPlus="1.000000e-02" value="4.000000e-02" errorMinus="1.000000e-02"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="1.250000e+01" value="1.375000e+02" errorMinus="1.250000e+01"/>
+ <measurement errorPlus="7.000000e-03" value="1.500000e-02" errorMinus="7.000000e-03"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorPlus="2.500000e+01" value="1.750000e+02" errorMinus="2.500000e+01"/>
+ <measurement errorPlus="1.900000e-03" value="2.300000e-03" errorMinus="1.900000e-03"/>
+ </dataPoint>
+ </dataPointSet>
+</aida>
Modified: branches/2011-07-aida2yoda/data/refdata/Makefile.am
==============================================================================
--- branches/2011-07-aida2yoda/data/refdata/Makefile.am Fri Jun 8 15:15:31 2012 (r3753)
+++ branches/2011-07-aida2yoda/data/refdata/Makefile.am Fri Jun 8 15:17:44 2012 (r3754)
@@ -43,6 +43,7 @@
ATLAS_2011_S9212183.aida \
ATLAS_2012_I1083318.aida \
ATLAS_2012_I1084540.aida \
+ ATLAS_2012_I1093738.aida \
ATLAS_2012_I1094568.aida \
ATLAS_2012_CONF_2012_001.aida \
BABAR_2005_S6181155.aida \
Copied and modified: branches/2011-07-aida2yoda/src/Analyses/ATLAS_2012_I1093738.cc (from r3716, trunk/src/Analyses/ATLAS_2012_I1093738.cc)
==============================================================================
--- trunk/src/Analyses/ATLAS_2012_I1093738.cc Tue May 1 20:31:32 2012 (r3716, copy source)
+++ branches/2011-07-aida2yoda/src/Analyses/ATLAS_2012_I1093738.cc Fri Jun 8 15:17:44 2012 (r3754)
@@ -4,7 +4,7 @@
#include <string>
#include "Rivet/Analysis.hh"
-#include "Rivet/RivetAIDA.hh"
+#include "Rivet/RivetYODA.hh"
#include "Rivet/Tools/Logging.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/LeadingParticlesFinalState.hh"
@@ -38,7 +38,6 @@
ATLAS_2012_I1093738()
: Analysis("ATLAS_2012_I1093738")
{
- setBeams(PROTON, PROTON);
setNeedsCrossSection(true);
_eta_bins_ph.push_back(0.0);
@@ -91,13 +90,13 @@
MYDEBUG << "\t... Booking Histograms " << std::endl;
// use autobook option, histo d??-x??-y??
- _h_phbarrel_jetcentral_SS = bookHistogram1D(1, 1, 1);
- _h_phbarrel_jetmedium_SS = bookHistogram1D(2, 1, 1);
- _h_phbarrel_jetforward_SS = bookHistogram1D(3, 1, 1);
-
- _h_phbarrel_jetcentral_OS = bookHistogram1D(4, 1, 1);
- _h_phbarrel_jetmedium_OS = bookHistogram1D(5, 1, 1);
- _h_phbarrel_jetforward_OS = bookHistogram1D(6, 1, 1);
+ _h_phbarrel_jetcentral_SS = bookHisto1D(1, 1, 1);
+ _h_phbarrel_jetmedium_SS = bookHisto1D(2, 1, 1);
+ _h_phbarrel_jetforward_SS = bookHisto1D(3, 1, 1);
+
+ _h_phbarrel_jetcentral_OS = bookHisto1D(4, 1, 1);
+ _h_phbarrel_jetmedium_OS = bookHisto1D(5, 1, 1);
+ _h_phbarrel_jetforward_OS = bookHisto1D(6, 1, 1);
MYDEBUG << "Exiting init." << std::endl;
}
@@ -108,22 +107,22 @@
int v_iter=0;
if (what==0) {
- for(v_iter=0; v_iter < (int)_eta_bins_ph.size()-1; v_iter++){
- if(eta >= _eta_bins_ph.at(v_iter) && eta < _eta_bins_ph.at(v_iter+1))
- break;
- }
+ for(v_iter=0; v_iter < (int)_eta_bins_ph.size()-1; v_iter++){
+ if(eta >= _eta_bins_ph.at(v_iter) && eta < _eta_bins_ph.at(v_iter+1))
+ break;
+ }
}
else if (what==1) {
- for (v_iter=0; v_iter < (int)_eta_bins_jet.size()-1; v_iter++){
- if(eta >= _eta_bins_jet.at(v_iter) && eta < _eta_bins_jet.at(v_iter+1))
- break;
- }
+ for (v_iter=0; v_iter < (int)_eta_bins_jet.size()-1; v_iter++){
+ if(eta >= _eta_bins_jet.at(v_iter) && eta < _eta_bins_jet.at(v_iter+1))
+ break;
+ }
}
else {
- for(v_iter=0; v_iter < (int)_eta_bins_areaoffset.size()-1; v_iter++){
- if(eta >= _eta_bins_areaoffset.at(v_iter) && eta < _eta_bins_areaoffset.at(v_iter+1))
- break;
- }
+ for(v_iter=0; v_iter < (int)_eta_bins_areaoffset.size()-1; v_iter++){
+ if(eta >= _eta_bins_areaoffset.at(v_iter) && eta < _eta_bins_areaoffset.at(v_iter+1))
+ break;
+ }
}
return v_iter;
@@ -143,7 +142,7 @@
// Get the photon
const FinalState& photonfs = applyProjection<FinalState>(event, "LeadingPhoton");
if (photonfs.particles().size() < 1) {
- MYDEBUG << "...Going to veto event(1)." << std::endl;
+ MYDEBUG << "...Going to veto event(1)." << std::endl;
vetoEvent;
}
MYDEBUG << "...Didn't veto event(1)." << std::endl;
@@ -155,7 +154,7 @@
// Get the jet
Jets jets = applyProjection<FastJets>(event, "Jets").jetsByPt(20.0*GeV);
if (jets.size()==0) {
- MYDEBUG << "...Going to veto event(3)." << std::endl;
+ MYDEBUG << "...Going to veto event(3)." << std::endl;
vetoEvent;
}
MYDEBUG << "...Didn't veto event(3)." << std::endl;
@@ -163,7 +162,7 @@
// Require jet separated from photon
if (deltaR(eta_P, phi_P, leadingJet.eta(), leadingJet.phi())<1.0) {
- MYDEBUG << "...Going to veto event(4)." << std::endl;
+ MYDEBUG << "...Going to veto event(4)." << std::endl;
vetoEvent;
}
MYDEBUG << "...Didn't veto event(4)." << std::endl;
@@ -171,57 +170,57 @@
// Veto if leading jet is outside plotted rapidity regions
const double abs_y1 = fabs(leadingJet.rapidity());
if (abs_y1 > 4.4) {
- MYDEBUG << "...Going to veto event(5)." << std::endl;
+ MYDEBUG << "...Going to veto event(5)." << std::endl;
vetoEvent;
}
MYDEBUG << "...Didn't veto event(5)." << std::endl;
// compute the median event energy density
- const unsigned int skipnhardjets = 0;
+ const unsigned int skipnhardjets = 0;
_ptDensity.clear();
_sigma.clear();
_Njets.clear();
std::vector< std::vector<double> > ptDensities;
std::vector<double> emptyVec;
ptDensities.assign(_eta_bins_areaoffset.size()-1,emptyVec);
-
+
const fastjet::ClusterSequenceArea* clust_seq_area = applyProjection<FastJets>(event, "KtJetsD05").clusterSeqArea();
foreach (const fastjet::PseudoJet& jet, applyProjection<FastJets>(event, "KtJetsD05").pseudoJets(0.0*GeV)) {
- double eta = fabs(jet.eta());
- double pt = fabs(jet.perp());
-
- /// get the cluster sequence
- double area = clust_seq_area->area(jet);
-
- if(area > 10e-4 && fabs(eta)<_eta_bins_areaoffset[_eta_bins_areaoffset.size()-1]){
- ptDensities.at(getEtaBin(fabs(eta),2)).push_back(pt/area);
- }
+ double eta = fabs(jet.eta());
+ double pt = fabs(jet.perp());
+
+ /// get the cluster sequence
+ double area = clust_seq_area->area(jet);
+
+ if(area > 10e-4 && fabs(eta)<_eta_bins_areaoffset[_eta_bins_areaoffset.size()-1]){
+ ptDensities.at(getEtaBin(fabs(eta),2)).push_back(pt/area);
+ }
}
-
+
for(int b=0; b<(int)_eta_bins_areaoffset.size()-1; b++){
- double median = 0.0;
- double sigma = 0.0;
- int Njets = 0;
- if(ptDensities[b].size() > skipnhardjets)
- {
- std::sort(ptDensities[b].begin(), ptDensities[b].end());
- int nDens = ptDensities[b].size() - skipnhardjets;
- if( nDens%2 == 0 )
- median = (ptDensities[b][nDens/2]+ptDensities[b][(nDens-2)/2])/2;
- else
- median = ptDensities[b][(nDens-1)/2];
- sigma = ptDensities[b][(int)(.15865*nDens)];
- Njets = nDens;
- }
- _ptDensity.push_back(median);
- _sigma.push_back(sigma);
- _Njets.push_back(Njets);
+ double median = 0.0;
+ double sigma = 0.0;
+ int Njets = 0;
+ if(ptDensities[b].size() > skipnhardjets)
+ {
+ std::sort(ptDensities[b].begin(), ptDensities[b].end());
+ int nDens = ptDensities[b].size() - skipnhardjets;
+ if( nDens%2 == 0 )
+ median = (ptDensities[b][nDens/2]+ptDensities[b][(nDens-2)/2])/2;
+ else
+ median = ptDensities[b][(nDens-1)/2];
+ sigma = ptDensities[b][(int)(.15865*nDens)];
+ Njets = nDens;
+ }
+ _ptDensity.push_back(median);
+ _sigma.push_back(sigma);
+ _Njets.push_back(Njets);
}
// compute photon isolation
-
+
// std EtCone
ParticleVector fs = applyProjection<FinalState>(event, "JetFS").particles();
FourMomentum mom_in_EtCone;
@@ -229,31 +228,30 @@
float cluster_eta_width = 0.25*7.0;
float cluster_phi_width = (PI/128.)*5.0;
foreach (const Particle& p, fs) {
- // check if it's in the cone of .4
- if (deltaR(eta_P, phi_P, p.momentum().eta(), p.momentum().phi()) >= iso_dR) continue;
-
- // check if it's in the 5x7 central core
- if (fabs(eta_P-p.momentum().eta()) < cluster_eta_width*0.5 &&
- fabs(phi_P-p.momentum().phi()) < cluster_phi_width*0.5) continue;
-
- mom_in_EtCone += p.momentum();
+ // check if it's in the cone of .4
+ if (deltaR(eta_P, phi_P, p.momentum().eta(), p.momentum().phi()) >= iso_dR) continue;
+
+ // check if it's in the 5x7 central core
+ if (fabs(eta_P-p.momentum().eta()) < cluster_eta_width*0.5 &&
+ fabs(phi_P-p.momentum().phi()) < cluster_phi_width*0.5) continue;
+
+ mom_in_EtCone += p.momentum();
}
MYDEBUG << "...Done with initial EtCone." << std::endl;
-
+
// now figure out the correction (area*density)
float EtCone_area = PI*iso_dR*iso_dR - cluster_eta_width*cluster_phi_width;
float correction = _ptDensity[getEtaBin(eta_P,2)]*EtCone_area;
MYDEBUG << "...Done with jet-area correction." << std::endl;
// require photon to be isolated
- if(mom_in_EtCone.Et()-correction > 4.0*GeV){
- MYDEBUG << "...Going to veto event(6)." << std::endl;
+ if(mom_in_EtCone.Et()-correction > 4.0*GeV){
+ MYDEBUG << "...Going to veto event(6)." << std::endl;
vetoEvent;
}
MYDEBUG << "...Didn't veto event(6)." << std::endl;
- //
int photon_jet_sign = sign( leadingJet.rapidity() * photon.rapidity() );
// Fill histos
@@ -263,30 +261,30 @@
if (abs_photon_eta<1.37) {
- if (abs_jet_rapidity < 1.2) {
+ if (abs_jet_rapidity < 1.2) {
- if (photon_jet_sign >= 1) {
- _h_phbarrel_jetcentral_SS->fill(photon_pt, weight);
- } else {
- _h_phbarrel_jetcentral_OS->fill(photon_pt, weight);
- }
-
- } else if (abs_jet_rapidity < 2.8) {
-
- if (photon_jet_sign >= 1) {
- _h_phbarrel_jetmedium_SS->fill(photon_pt, weight);
- } else {
- _h_phbarrel_jetmedium_OS->fill(photon_pt, weight);
- }
-
- } else if (abs_jet_rapidity < 4.4) {
-
- if (photon_jet_sign >= 1) {
- _h_phbarrel_jetforward_SS->fill(photon_pt, weight);
- } else {
- _h_phbarrel_jetforward_OS->fill(photon_pt, weight);
- }
- }
+ if (photon_jet_sign >= 1) {
+ _h_phbarrel_jetcentral_SS->fill(photon_pt, weight);
+ } else {
+ _h_phbarrel_jetcentral_OS->fill(photon_pt, weight);
+ }
+
+ } else if (abs_jet_rapidity < 2.8) {
+
+ if (photon_jet_sign >= 1) {
+ _h_phbarrel_jetmedium_SS->fill(photon_pt, weight);
+ } else {
+ _h_phbarrel_jetmedium_OS->fill(photon_pt, weight);
+ }
+
+ } else if (abs_jet_rapidity < 4.4) {
+
+ if (photon_jet_sign >= 1) {
+ _h_phbarrel_jetforward_SS->fill(photon_pt, weight);
+ } else {
+ _h_phbarrel_jetforward_OS->fill(photon_pt, weight);
+ }
+ }
}
@@ -311,13 +309,13 @@
private:
- AIDA::IHistogram1D* _h_phbarrel_jetcentral_SS;
- AIDA::IHistogram1D* _h_phbarrel_jetmedium_SS;
- AIDA::IHistogram1D* _h_phbarrel_jetforward_SS;
-
- AIDA::IHistogram1D* _h_phbarrel_jetcentral_OS;
- AIDA::IHistogram1D* _h_phbarrel_jetmedium_OS;
- AIDA::IHistogram1D* _h_phbarrel_jetforward_OS;
+ Histo1DPtr _h_phbarrel_jetcentral_SS;
+ Histo1DPtr _h_phbarrel_jetmedium_SS;
+ Histo1DPtr _h_phbarrel_jetforward_SS;
+
+ Histo1DPtr _h_phbarrel_jetcentral_OS;
+ Histo1DPtr _h_phbarrel_jetmedium_OS;
+ Histo1DPtr _h_phbarrel_jetforward_OS;
fastjet::AreaDefinition* _area_def;
@@ -330,8 +328,7 @@
std::vector<float> _Njets;
};
+ // The hook for the plugin system
+ DECLARE_RIVET_PLUGIN(ATLAS_2012_I1093738);
-
- // This global object acts as a hook for the plugin system
- AnalysisBuilder<ATLAS_2012_I1093738> plugin_ATLAS_2012_I1093738;
}
Modified: branches/2011-07-aida2yoda/src/Analyses/Makefile.am
==============================================================================
--- branches/2011-07-aida2yoda/src/Analyses/Makefile.am Fri Jun 8 15:15:31 2012 (r3753)
+++ branches/2011-07-aida2yoda/src/Analyses/Makefile.am Fri Jun 8 15:17:44 2012 (r3754)
@@ -67,6 +67,7 @@
ATLAS_2011_I926145.cc \
ATLAS_2011_I944826.cc \
ATLAS_2012_I1084540.cc \
+ ATLAS_2012_I1093738.cc \
ATLAS_2012_I1094568.cc
if ENABLE_PRELIMINARY
RivetATLASAnalyses_la_SOURCES += \
More information about the Rivet-svn
mailing list
|