|
[Rivet-svn] r2751 - in trunk: . data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgWed Nov 10 16:03:27 GMT 2010
Author: hoeth Date: Wed Nov 10 16:03:27 2010 New Revision: 2751 Log: Adding ATLAS-CONF-2010-083 conference note (azimuthal dijet decorrelation). Data is read from plots. When I run Pythia 6 the bins close to pi/2 are higher than in the note, so I call this "unvalidated". But then ... the note doesn't specify a tune or even just a version for the generators in the plots. Not even if they used Pythia 6 or Pythia 8. Probably 6, since they mention AGILe and cite the Py6 manual. Added: trunk/data/anainfo/ATLAS_2010_CONF_2010_083.info trunk/data/plotinfo/ATLAS_2010_CONF_2010_083.plot trunk/data/refdata/ATLAS_2010_CONF_2010_083.aida trunk/src/Analyses/ATLAS_2010_CONF_2010_083.cc Modified: trunk/ChangeLog trunk/data/anainfo/Makefile.am trunk/data/plotinfo/Makefile.am trunk/data/refdata/Makefile.am trunk/src/Analyses/Makefile.am Modified: trunk/ChangeLog ============================================================================== --- trunk/ChangeLog Wed Nov 10 15:59:14 2010 (r2750) +++ trunk/ChangeLog Wed Nov 10 16:03:27 2010 (r2751) @@ -1,3 +1,11 @@ +2010-11-10 Hendrik Hoeth <hendrik.hoeth at cern.ch> + + * Adding ATLAS-CONF-2010-083 conference note. Data is read from plots. + When I run Pythia 6 the bins close to pi/2 are higher than in the note, + so I call this "unvalidated". But then ... the note doesn't specify + a tune or even just a version for the generators in the plots. Not even + if they used Pythia 6 or Pythia 8. Probably 6, since they mention AGILe. + 2010-11-10 Andy Buckley <andy at insectnation.org> * Adding a JetAlg::useInvisibles(bool) mechanism to allow ATLAS Added: trunk/data/anainfo/ATLAS_2010_CONF_2010_083.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/anainfo/ATLAS_2010_CONF_2010_083.info Wed Nov 10 16:03:27 2010 (r2751) @@ -0,0 +1,21 @@ +Name: ATLAS_2010_CONF_2010_083 +Year: 2010 +Summary: Dijet azimuthal decorrelation at 7000 GeV in ATLAS +Experiment: ATLAS +Collider: LHC 7000GeV +SpiresID: None +Status: UNVALIDATED -- but PRELIMINARY DATA! +Authors: + - Hendrik Hoeth <hendrik.hoeth at cern.ch> +References: + - ATLAS-CONF-2010-083 +RunInfo: + pp QCD interactions at 7000 GeV including diffractive events. +NumEvents: 1000000 +Beams: [p+, p+] +Energies: [7000] +PtCuts: [0] +Description: + Dijet azimuthal decorrelation measured by ATLAS at 7 TeV. Jets are anti-kt with + $R=0.6$, $p_\perp>100$ GeV, $|\eta|<0.8$. The analysis is binned in four leading + jet $p_\perp$ bins. All data is fully corrected. ATTENTION - Data read from plots! Modified: trunk/data/anainfo/Makefile.am ============================================================================== --- trunk/data/anainfo/Makefile.am Wed Nov 10 15:59:14 2010 (r2750) +++ trunk/data/anainfo/Makefile.am Wed Nov 10 16:03:27 2010 (r2751) @@ -8,6 +8,7 @@ ATLAS_2010_CONF_2010_031.info \ ATLAS_2010_CONF_2010_049.info \ ATLAS_2010_CONF_2010_081.info \ + ATLAS_2010_CONF_2010_083.info \ BELLE_2006_S6265367.info \ CDF_1988_S1865951.info \ CDF_1990_S2089246.info \ Added: trunk/data/plotinfo/ATLAS_2010_CONF_2010_083.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/plotinfo/ATLAS_2010_CONF_2010_083.plot Wed Nov 10 16:03:27 2010 (r2751) @@ -0,0 +1,32 @@ +# BEGIN PLOT /ATLAS_2010_CONF_2010_083/d01-x01-y01 +Title=Jet--jet azimuthal angle, $p_\perp^\text{max} \in [110,160]\text{ GeV}$ +XLabel=$\Delta \phi_\text{dijet}$ / rad +YLabel=$1/\sigma_\text{dijet} \, \text{d}\sigma_\text{dijet}/\text{d}\Delta\phi_\text{dijet}$ +LegendXPos=0.50 +CustomLegend={\large ATLAS preliminary} +# END PLOT + +# BEGIN PLOT /ATLAS_2010_CONF_2010_083/d02-x01-y01 +Title=Jet--jet azimuthal angle, $p_\perp^\text{max} \in [160,210]\text{ GeV}$ +XLabel=$\Delta \phi_\text{dijet}$ / rad +YLabel=$1/\sigma_\text{dijet} \, \text{d}\sigma_\text{dijet}/\text{d}\Delta\phi_\text{dijet}$ +LegendXPos=0.50 +CustomLegend={\large ATLAS preliminary} +# END PLOT + +# BEGIN PLOT /ATLAS_2010_CONF_2010_083/d03-x01-y01 +Title=Jet--jet azimuthal angle, $p_\perp^\text{max} \in [210,310]\text{ GeV}$ +XLabel=$\Delta \phi_\text{dijet}$ / rad +YLabel=$1/\sigma_\text{dijet} \, \text{d}\sigma_\text{dijet}/\text{d}\Delta\phi_\text{dijet}$ +LegendXPos=0.50 +CustomLegend={\large ATLAS preliminary} +# END PLOT + +# BEGIN PLOT /ATLAS_2010_CONF_2010_083/d04-x01-y01 +Title=Jet--jet azimuthal angle, $p_\perp^\text{max} > 310\text{ GeV}$ +XLabel=$\Delta \phi_\text{dijet}$ / rad +YLabel=$1/\sigma_\text{dijet} \, \text{d}\sigma_\text{dijet}/\text{d}\Delta\phi_\text{dijet}$ +LegendXPos=0.50 +CustomLegend={\large ATLAS preliminary} +# END PLOT + Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Wed Nov 10 15:59:14 2010 (r2750) +++ trunk/data/plotinfo/Makefile.am Wed Nov 10 16:03:27 2010 (r2751) @@ -8,6 +8,7 @@ ATLAS_2010_CONF_2010_031.plot \ ATLAS_2010_CONF_2010_049.plot \ ATLAS_2010_CONF_2010_081.plot \ + ATLAS_2010_CONF_2010_083.plot \ CDF_1988_S1865951.plot \ CDF_1990_S2089246.plot \ CDF_1991_S2313472.plot \ Added: trunk/data/refdata/ATLAS_2010_CONF_2010_083.aida ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/refdata/ATLAS_2010_CONF_2010_083.aida Wed Nov 10 16:03:27 2010 (r2751) @@ -0,0 +1,137 @@ +<?xml version="1.0" encoding="ISO-8859-1" ?> +<!DOCTYPE aida SYSTEM "http://aida.freehep.org/schemas/3.3/aida.dtd"> +<aida version="3.3"> + <implementation version="1.1" package="FreeHEP"/> + <dataPointSet name="d01-x01-y01" dimension="2" + path="/REF/ATLAS_2010_CONF_2010_083" title=""> + <annotation> + <item key="Title" value="" sticky="true"/> + <item key="AidaPath" value="/REF/ATLAS_2010_CONF_2010_083/d01-x01-y01" sticky="true"/> + <item key="FullPath" value="/ATLAS_2010_CONF_2010_083.aida/REF/ATLAS_2010_ATLAS_CONF_2010_083/d01-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="1.884955e-01" value="1.759292e+00" errorMinus="1.884955e-01"/> + <measurement errorPlus="1.740000e-03" value="4.146000e-03" errorMinus="1.740000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.570797e-01" value="2.104867e+00" errorMinus="1.570797e-01"/> + <measurement errorPlus="6.520000e-03" value="2.801000e-02" errorMinus="6.520000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.413716e-01" value="2.403318e+00" errorMinus="1.413716e-01"/> + <measurement errorPlus="1.190000e-02" value="1.134000e-01" errorMinus="1.190000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.099557e-01" value="2.654646e+00" errorMinus="1.099557e-01"/> + <measurement errorPlus="2.640000e-02" value="3.532000e-01" errorMinus="2.640000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="7.853980e-02" value="2.843141e+00" errorMinus="7.853980e-02"/> + <measurement errorPlus="6.270000e-02" value="9.249000e-01" errorMinus="6.270000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="6.283185e-02" value="2.984513e+00" errorMinus="6.283185e-02"/> + <measurement errorPlus="7.870000e-02" value="2.315000e+00" errorMinus="7.870000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="4.712390e-02" value="3.094469e+00" errorMinus="4.712390e-02"/> + <measurement errorPlus="1.940000e-01" value="4.695000e+00" errorMinus="1.940000e-01"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d02-x01-y01" dimension="2" + path="/REF/ATLAS_2010_CONF_2010_083" title=""> + <annotation> + <item key="Title" value="" sticky="true"/> + <item key="AidaPath" value="/REF/ATLAS_2010_CONF_2010_083/d02-x01-y01" sticky="true"/> + <item key="FullPath" value="/ATLAS_2010_CONF_2010_083.aida/REF/ATLAS_2010_ATLAS_CONF_2010_083/d02-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="1.884955e-01" value="1.759292e+00" errorMinus="1.884955e-01"/> + <measurement errorPlus="3.060000e-03" value="4.923000e-03" errorMinus="3.060000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.570797e-01" value="2.104867e+00" errorMinus="1.570797e-01"/> + <measurement errorPlus="9.750000e-03" value="3.427000e-02" errorMinus="9.750000e-03"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.413716e-01" value="2.403318e+00" errorMinus="1.413716e-01"/> + <measurement errorPlus="2.700000e-02" value="1.779000e-01" errorMinus="2.700000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.099557e-01" value="2.654646e+00" errorMinus="1.099557e-01"/> + <measurement errorPlus="4.850000e-02" value="4.487000e-01" errorMinus="4.850000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="7.853980e-02" value="2.843141e+00" errorMinus="7.853980e-02"/> + <measurement errorPlus="7.260000e-02" value="1.026000e+00" errorMinus="7.260000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="6.283185e-02" value="2.984513e+00" errorMinus="6.283185e-02"/> + <measurement errorPlus="1.210000e-01" value="2.312000e+00" errorMinus="1.210000e-01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="4.712390e-02" value="3.094469e+00" errorMinus="4.712390e-02"/> + <measurement errorPlus="2.160000e-01" value="4.095000e+00" errorMinus="2.160000e-01"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d03-x01-y01" dimension="2" + path="/REF/ATLAS_2010_CONF_2010_083" title=""> + <annotation> + <item key="Title" value="" sticky="true"/> + <item key="AidaPath" value="/REF/ATLAS_2010_CONF_2010_083/d03-x01-y01" sticky="true"/> + <item key="FullPath" value="/ATLAS_2010_CONF_2010_083.aida/REF/ATLAS_2010_ATLAS_CONF_2010_083/d03-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="1.570797e-01" value="2.104867e+00" errorMinus="1.570797e-01"/> + <measurement errorPlus="1.220000e-02" value="2.923000e-02" errorMinus="1.220000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.413716e-01" value="2.403318e+00" errorMinus="1.413716e-01"/> + <measurement errorPlus="2.920000e-02" value="1.193000e-01" errorMinus="2.920000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.099557e-01" value="2.654646e+00" errorMinus="1.099557e-01"/> + <measurement errorPlus="6.320000e-02" value="4.724000e-01" errorMinus="6.320000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="7.853980e-02" value="2.843141e+00" errorMinus="7.853980e-02"/> + <measurement errorPlus="1.000000e-01" value="8.953000e-01" errorMinus="1.000000e-01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="6.283185e-02" value="2.984513e+00" errorMinus="6.283185e-02"/> + <measurement errorPlus="1.670000e-01" value="2.142000e+00" errorMinus="1.670000e-01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="4.712390e-02" value="3.094469e+00" errorMinus="4.712390e-02"/> + <measurement errorPlus="3.170000e-01" value="4.684000e+00" errorMinus="3.170000e-01"/> + </dataPoint> + </dataPointSet> + <dataPointSet name="d04-x01-y01" dimension="2" + path="/REF/ATLAS_2010_CONF_2010_083" title=""> + <annotation> + <item key="Title" value="" sticky="true"/> + <item key="AidaPath" value="/REF/ATLAS_2010_CONF_2010_083/d04-x01-y01" sticky="true"/> + <item key="FullPath" value="/ATLAS_2010_CONF_2010_083.aida/REF/ATLAS_2010_ATLAS_CONF_2010_083/d04-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="1.413716e-01" value="2.403318e+00" errorMinus="1.413716e-01"/> + <measurement errorPlus="6.230000e-02" value="1.219000e-01" errorMinus="6.230000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.099557e-01" value="2.654646e+00" errorMinus="1.099557e-01"/> + <measurement errorPlus="1.100000e-01" value="2.851000e-01" errorMinus="1.100000e-01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="7.853980e-02" value="2.843141e+00" errorMinus="7.853980e-02"/> + <measurement errorPlus="2.810000e-01" value="1.264000e+00" errorMinus="2.810000e-01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="6.283185e-02" value="2.984513e+00" errorMinus="6.283185e-02"/> + <measurement errorPlus="3.820000e-01" value="1.955000e+00" errorMinus="3.820000e-01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="4.712390e-02" value="3.094469e+00" errorMinus="4.712390e-02"/> + <measurement errorPlus="6.750000e-01" value="4.857000e+00" errorMinus="6.750000e-01"/> + </dataPoint> + </dataPointSet> +</aida> Modified: trunk/data/refdata/Makefile.am ============================================================================== --- trunk/data/refdata/Makefile.am Wed Nov 10 15:59:14 2010 (r2750) +++ trunk/data/refdata/Makefile.am Wed Nov 10 16:03:27 2010 (r2751) @@ -9,6 +9,7 @@ ATLAS_2010_CONF_2010_031.aida \ ATLAS_2010_CONF_2010_049.aida \ ATLAS_2010_CONF_2010_081.aida \ + ATLAS_2010_CONF_2010_083.aida \ DELPHI_1994_S3021912.aida \ DELPHI_1995_S3137023.aida \ DELPHI_1996_S3430090.aida \ Added: trunk/src/Analyses/ATLAS_2010_CONF_2010_083.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/src/Analyses/ATLAS_2010_CONF_2010_083.cc Wed Nov 10 16:03:27 2010 (r2751) @@ -0,0 +1,80 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/RivetAIDA.hh" +#include "Rivet/Projections/FinalState.hh" +#include "Rivet/Projections/FastJets.hh" +#include "Rivet/Tools/Logging.hh" +#include "LWH/Histogram1D.h" + +namespace Rivet { + + class ATLAS_2010_CONF_2010_083 : public Analysis { + public: + + ATLAS_2010_CONF_2010_083() : Analysis("ATLAS_2010_CONF_2010_083") { + setNeedsCrossSection(true); + } + + + void init() { + FinalState fs(-1.5, 1.5, 0.0*GeV); + FastJets jetsproj(fs, FastJets::ANTIKT, 0.6); + addProjection(jetsproj, "Jets"); + + for (size_t i=0 ; i<4 ; i++) { + _h_dphi[i] = bookHistogram1D(i+1, 1, 1); + } + } + + + void analyze(const Event& event) { + const double weight = event.weight(); + + const FastJets & jetsproj = applyProjection<FastJets>(event, "Jets"); + Jets alljets = jetsproj.jetsByPt(100.0*GeV); + Jets jets; + foreach (const Jet jet, alljets) { + if (fabs(jet.momentum().eta())<0.8) { + jets.push_back(jet); + } + if (jets.size()==2) break; + } + + if (jets.size() < 2 || jets[0].momentum().pT() < 110.0*GeV) { + vetoEvent; + } + + const double dphi = deltaPhi(jets[0].momentum().phi(), jets[1].momentum().phi()); + const double ptlead = jets[0].momentum().pT(); + + if (ptlead > 310.0*GeV) { + _h_dphi[3]->fill(dphi, weight); + } + else if (ptlead > 210.0*GeV) { + _h_dphi[2]->fill(dphi, weight); + } + else if (ptlead > 160.0*GeV) { + _h_dphi[1]->fill(dphi, weight); + } + else { + _h_dphi[0]->fill(dphi, weight); + } + } + + void finalize() { + for (size_t i=0 ; i<4 ; i++) { + normalize(_h_dphi[i], 1.0); + } + } + + + private: + + AIDA::IHistogram1D* _h_dphi[4]; + }; + + // This global object acts as a hook for the plugin system + AnalysisBuilder<ATLAS_2010_CONF_2010_083> plugin_ATLAS_2010_CONF_2010_083; + +} + Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Wed Nov 10 15:59:14 2010 (r2750) +++ trunk/src/Analyses/Makefile.am Wed Nov 10 16:03:27 2010 (r2751) @@ -50,7 +50,8 @@ ATLAS_2010_S8817804.cc \ ATLAS_2010_CONF_2010_031.cc \ ATLAS_2010_CONF_2010_049.cc \ - ATLAS_2010_CONF_2010_081.cc + ATLAS_2010_CONF_2010_081.cc \ + ATLAS_2010_CONF_2010_083.cc lib_LTLIBRARIES += RivetCDFAnalyses.la
More information about the Rivet-svn mailing list |