|
[Rivet-svn] r3716 - in trunk: data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgTue May 1 20:31:32 BST 2012
Author: holsch Date: Tue May 1 20:31:32 2012 New Revision: 3716 Log: Add ATLAS photon+jets differential cross-section analysis Added: trunk/data/anainfo/ATLAS_2012_I1093738.info trunk/data/refdata/ATLAS_2012_I1093738.aida trunk/src/Analyses/ATLAS_2012_I1093738.cc Modified: trunk/data/anainfo/Makefile.am trunk/data/plotinfo/Makefile.am trunk/data/refdata/Makefile.am trunk/src/Analyses/Makefile.am Added: trunk/data/anainfo/ATLAS_2012_I1093738.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/anainfo/ATLAS_2012_I1093738.info Tue May 1 20:31:32 2012 (r3716) @@ -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: trunk/data/anainfo/Makefile.am ============================================================================== --- trunk/data/anainfo/Makefile.am Tue May 1 20:11:23 2012 (r3715) +++ trunk/data/anainfo/Makefile.am Tue May 1 20:31:32 2012 (r3716) @@ -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 \ Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Tue May 1 20:11:23 2012 (r3715) +++ trunk/data/plotinfo/Makefile.am Tue May 1 20:31:32 2012 (r3716) @@ -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 \ Added: trunk/data/refdata/ATLAS_2012_I1093738.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/ATLAS_2012_I1093738.aida Tue May 1 20:31:32 2012 (r3716) @@ -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: trunk/data/refdata/Makefile.am ============================================================================== --- trunk/data/refdata/Makefile.am Tue May 1 20:11:23 2012 (r3715) +++ trunk/data/refdata/Makefile.am Tue May 1 20:31:32 2012 (r3716) @@ -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 \ Added: trunk/src/Analyses/ATLAS_2012_I1093738.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/src/Analyses/ATLAS_2012_I1093738.cc Tue May 1 20:31:32 2012 (r3716) @@ -0,0 +1,337 @@ +// -*- C++ -*- +#include <iostream> +#include <sstream> +#include <string> + +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Tools/Logging.hh" +#include "Rivet/Projections/FinalState.hh" +#include "Rivet/Projections/LeadingParticlesFinalState.hh" +#include "Rivet/Projections/VetoedFinalState.hh" +#include "Rivet/Jet.hh" +#include "Rivet/Projections/FastJets.hh" + +#include "fastjet/internal/base.hh" +#include "fastjet/JetDefinition.hh" +#include "fastjet/AreaDefinition.hh" +#include "fastjet/ClusterSequence.hh" +#include "fastjet/ClusterSequenceArea.hh" +#include "fastjet/PseudoJet.hh" + +#define MYDEBUG if(false) getLog() << __LINE__ + + +namespace Rivet { + + /// @brief Measurement of isolated gamma + jet + X differential cross-sections + /// + /// Inclusive isolated gamma + jet cross-sections, differential in pT(gamma), for + /// various photon and jet rapidity configurations. + /// + /// @author Giovanni Marchiori + + class ATLAS_2012_I1093738 : public Analysis { + public: + + // Constructor + ATLAS_2012_I1093738() + : Analysis("ATLAS_2012_I1093738") + { + setBeams(PROTON, PROTON); + setNeedsCrossSection(true); + + _eta_bins_ph.push_back(0.0); + _eta_bins_ph.push_back(1.37); + _eta_bins_ph.push_back(1.52); + _eta_bins_ph.push_back(2.37); + + _eta_bins_jet.push_back(0.0); + _eta_bins_jet.push_back(1.2); + _eta_bins_jet.push_back(2.8); + _eta_bins_jet.push_back(4.4); + + _eta_bins_areaoffset.push_back(0.0); + _eta_bins_areaoffset.push_back(1.5); + _eta_bins_areaoffset.push_back(3.0); + } + + public: + + // Book histograms and initialise projections before the run + void init() { + MYDEBUG << "Entering init." << std::endl; + + // Final state + FinalState fs; + addProjection(fs, "FS"); + + // Voronoi eta-phi tassellation with KT jets, for ambient energy density calculation + FastJets fj(fs, FastJets::KT, 0.5); + _area_def = new fastjet::AreaDefinition(fastjet::VoronoiAreaSpec()); + fj.useJetArea(_area_def); + addProjection(fj, "KtJetsD05"); + + // Leading photon + LeadingParticlesFinalState photonfs(FinalState(-1.37, 1.37, 25.0*GeV)); + photonfs.addParticleId(PHOTON); + addProjection(photonfs, "LeadingPhoton"); + + // FS excluding the leading photon + VetoedFinalState vfs(fs); + vfs.addVetoOnThisFinalState(photonfs); + addProjection(vfs, "JetFS"); + + // Jets + FastJets jetpro(vfs, FastJets::ANTIKT, 0.4); + //FastJets jetpro(vfs, FastJets::ANTIKT, 0.6); + jetpro.useInvisibles(); + addProjection(jetpro, "Jets"); + + + 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); + + MYDEBUG << "Exiting init." << std::endl; + } + + + int getEtaBin(double eta_w, int what) const { + double eta = fabs(eta_w); + + 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; + } + } + 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; + } + } + 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; + } + } + + return v_iter; + } + + + // Perform the per-event analysis + void analyze(const Event& event) { + + MYDEBUG << "Entering Analyze." << std::endl; + + const double weight = event.weight(); + + MYDEBUG << "...First projections." << std::endl; + + + // Get the photon + const FinalState& photonfs = applyProjection<FinalState>(event, "LeadingPhoton"); + if (photonfs.particles().size() < 1) { + MYDEBUG << "...Going to veto event(1)." << std::endl; + vetoEvent; + } + MYDEBUG << "...Didn't veto event(1)." << std::endl; + + const FourMomentum photon = photonfs.particles().front().momentum(); + double eta_P = photon.eta(); + double phi_P = photon.phi(); + + // 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; + vetoEvent; + } + MYDEBUG << "...Didn't veto event(3)." << std::endl; + FourMomentum leadingJet = jets[0].momentum(); + + // 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; + vetoEvent; + } + MYDEBUG << "...Didn't veto event(4)." << std::endl; + + // 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; + vetoEvent; + } + MYDEBUG << "...Didn't veto event(5)." << std::endl; + + + // compute the median event energy density + 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); + } + } + + 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); + } + + + // compute photon isolation + + // std EtCone + ParticleVector fs = applyProjection<FinalState>(event, "JetFS").particles(); + FourMomentum mom_in_EtCone; + float iso_dR = 0.4; + 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(); + } + 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; + vetoEvent; + } + MYDEBUG << "...Didn't veto event(6)." << std::endl; + + + // + int photon_jet_sign = sign( leadingJet.rapidity() * photon.rapidity() ); + + // Fill histos + float abs_jet_rapidity = fabs(leadingJet.rapidity()); + float photon_pt = photon.pT()/GeV; + float abs_photon_eta = fabs(photon.eta()); + + if (abs_photon_eta<1.37) { + + 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); + } + } + + } + + MYDEBUG << "...Done with analyze." << std::endl; + } + + + /// Normalise histograms etc., after the run + void finalize() { + + MYDEBUG << "In finalize." << std::endl; + + scale(_h_phbarrel_jetcentral_SS, crossSection()/sumOfWeights()); + scale(_h_phbarrel_jetcentral_OS, crossSection()/sumOfWeights()); + scale(_h_phbarrel_jetmedium_SS, crossSection()/sumOfWeights()); + scale(_h_phbarrel_jetmedium_OS, crossSection()/sumOfWeights()); + scale(_h_phbarrel_jetforward_SS, crossSection()/sumOfWeights()); + scale(_h_phbarrel_jetforward_OS, crossSection()/sumOfWeights()); + + MYDEBUG << "Done with finalize." << std::endl; + } + + 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; + + fastjet::AreaDefinition* _area_def; + + std::vector<float> _eta_bins_ph; + std::vector<float> _eta_bins_jet; + std::vector<float> _eta_bins_areaoffset; + + std::vector<float> _ptDensity; + std::vector<float> _sigma; + std::vector<float> _Njets; + }; + + + + // This global object acts as a hook for the plugin system + AnalysisBuilder<ATLAS_2012_I1093738> plugin_ATLAS_2012_I1093738; +} Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Tue May 1 20:11:23 2012 (r3715) +++ trunk/src/Analyses/Makefile.am Tue May 1 20:31:32 2012 (r3716) @@ -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 |