ATLAS_2014_I1268975.aida 0000644 0001750 0001750 00000103340 12357717236 013314 0 ustar lysak lysak
ATLAS_2014_I1268975.cc 0000644 0001750 0001750 00000006214 12357720137 012777 0 ustar lysak lysak // -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/RivetAIDA.hh"
#include "Rivet/Tools/Logging.hh"
#include "Rivet/Projections/FastJets.hh"
#include "Rivet/Tools/BinnedHistogram.hh"
#include "Rivet/Projections/FinalState.hh"
namespace Rivet {
class ATLAS_2014_I1268975 : public Analysis {
public:
/// @name Constructors etc.
//@{
/// Constructor
ATLAS_2014_I1268975()
: Analysis("ATLAS_2014_I1268975")
{
}
//@}
public:
/// @name Analysis methods
//@{
/// Book histograms and initialise projections before the run
void init() {
const FinalState fs;
addProjection(fs,"FinalState");
FastJets fj04(fs, FastJets::ANTIKT, 0.4);
fj04.useInvisibles();
addProjection(fj04, "AntiKT04");
FastJets fj06(fs, FastJets::ANTIKT, 0.6);
fj06.useInvisibles();
addProjection(fj06, "AntiKT06");
double ystarbins[] = { 0.0, 0.5, 1.0, 1.5, 2.0, 2.5, 3.0};
size_t massDsOffset(0);
for (size_t alg = 0; alg < 2; ++alg) {
for (size_t i = 0; i < 6; ++i) {
_mass[alg].addHistogram(ystarbins[i], ystarbins[i+1], bookHistogram1D(1 + massDsOffset, 1, i+1));
}
massDsOffset += 1;
}
}
/// Perform the per-event analysis
void analyze(const Event& event) {
const double weight = event.weight();
Jets jetAr[2];
jetAr[AKT4] = applyProjection(event, "AntiKT04").jetsByPt(50*GeV);
jetAr[AKT6] = applyProjection(event, "AntiKT06").jetsByPt(50*GeV);
// Loop over jet "radii" used in analysis
for (size_t alg = 0; alg < 2; ++alg) {
// Identify dijets
vector leadjets;
foreach (const Jet& jet, jetAr[alg]) {
const double absy = fabs(jet.momentum().rapidity());
const double pT = jet.momentum().pT();
if (absy < 3.0 && leadjets.size() < 2) {
if ( leadjets.empty() && pT < 100*GeV) continue;
leadjets.push_back(jet.momentum());
}
}
// Make sure we have a leading jet with pT > 100 GeV and a second to leading jet with pT > 50 GeV
if (leadjets.size() < 2) {
MSG_DEBUG("Could not find two suitable leading jets");
continue;
}
const double y1 = leadjets[0].rapidity();
const double y2 = leadjets[1].rapidity();
const double ystar = fabs( y1-y2 ) / 2.;
const double m = (leadjets[0] + leadjets[1]).mass();
// Fill mass histogram
_mass[alg].fill(ystar, m/TeV, weight);
}
}
/// Normalise histograms etc., after the run
void finalize() {
for (size_t alg = 0; alg < 2; ++alg) {
_mass[alg].scale(crossSectionPerEvent()/picobarn, this);
}
}
//@}
private:
// Data members like post-cuts event weight counters go here
enum Alg { AKT4=0, AKT6=1 };
private:
/// The di-jet mass spectrum binned in rapidity for akt6 and akt4 jets (array index is jet type from enum above)
BinnedHistogram _mass[2];
};
// The hook for the plugin system
DECLARE_RIVET_PLUGIN(ATLAS_2014_I1268975);
}
ATLAS_2014_I1268975.info 0000644 0001750 0001750 00000005127 12357724153 013351 0 ustar lysak lysak Name: ATLAS_2014_I1268975
Year: 2014
Summary: high mass dijet cross section
Experiment: ATLAS
Collider: LHC
SpiresID: 1268975
Status: VALIDATED
Authors:
- Christopher Meyer
References:
- arXiv:1312.3524 [hep-ex]
- JHEP 1405 (2014) 059
RunInfo:
QCD jet production with a minimum leading jet pT of 100 GeV and minimum second jet pT of 50 GeV at 7 TeV.
NeedCrossSection: True
NumEvents: 10000000
Beams: [p+, p+]
Energies: [7000]
PtCuts: [100,50]
Description:
'Double-differential dijet cross sections measured in pp collisions
at the LHC with a 7 TeV centre-of-mass energy are presented as
functions of dijet mass and rapidity separation of the two highest-pT
jets. These measurements are obtained using data corresponding to an
integrated luminosity of 4.5/fb, recorded by the ATLAS detector in
2011. The data are corrected for detector effects so that cross
sections are presented at the particle level. Cross sections are
measured up to 5 TeV dijet mass using jets reconstructed with the
anti-kt algorithm for values of the jet radius parameter of 0.4 and
0.6. The cross sections are compared with next-to-leading-order
perturbative QCD calculations by NLOJET++ corrected to account for
non-perturbative effects. Comparisons with POWHEG predictions, using a
next-to-leading-order matrix element calculation interfaced to a
parton-shower Monte Carlo simulation, are also shown. Electroweak
effects are accounted for in both cases. The quantitative comparison
of data and theoretical predictions obtained using various
parameterizations of the parton distribution functions is performed
using a frequentist method. An example setting a lower limit on the
compositeness scale for a model of contact interactions is presented,
showing that the unfolded results can be used to constrain
contributions to dijet production beyond that predicted by the
Standard Model.'
BibKey: Aad:2013tea
BibTeX: '@article{Aad:2013tea,
author = "Aad, Georges and others",
title = "{Measurement of dijet cross sections in $pp$ collisions
at 7 TeV centre-of-mass energy using the ATLAS detector}",
collaboration = "ATLAS Collaboration",
journal = "JHEP",
volume = "1405",
pages = "059",
doi = "10.1007/JHEP05(2014)059",
year = "2014",
eprint = "1312.3524",
archivePrefix = "arXiv",
primaryClass = "hep-ex",
reportNumber = "CERN-PH-EP-2013-192",
SLACcitation = "%%CITATION = ARXIV:1312.3524;%%",
}'
ATLAS_2014_I1268975.plot 0000644 0001750 0001750 00000003277 12357714344 013401 0 ustar lysak lysak # BEGIN PLOT /ATLAS_2014_I1268975/d.*-x01-y.*
XLabel=$m_{12} [TeV]$
YLabel=$d^2\sigma/dm_{12}dy^{*}$ [pb/TeV]
LogY=1
# END PLOT
# BEGIN PLOT /ATLAS_2014_I1268975/d01-x01-y01
Title=Dijet double-differential cross sections ($y^{*} < 0.5$), R=0.4
# END PLOT
# BEGIN PLOT /ATLAS_2014_I1268975/d01-x01-y02
Title=Dijet double-differential cross sections ($0.5 < y^{*} < 1.0$), R=0.4
# END PLOT
# BEGIN PLOT /ATLAS_2014_I1268975/d01-x01-y03
Title=Dijet double-differential cross sections ($1.0 < y^{*} < 1.5$), R=0.4
# END PLOT
# BEGIN PLOT /ATLAS_2014_I1268975/d01-x01-y04
Title=Dijet double-differential cross sections ($1.5 < y^{*} < 2.0$), R=0.4
# END PLOT
# BEGIN PLOT /ATLAS_2014_I1268975/d01-x01-y05
Title=Dijet double-differential cross sections ($2.0 < y^{*} < 2.5$), R=0.4
# END PLOT
# BEGIN PLOT /ATLAS_2014_I1268975/d01-x01-y06
Title=Dijet double-differential cross sections ($2.5 < y^{*} < 3.0$), R=0.4
# END PLOT
# BEGIN PLOT /ATLAS_2014_I1268975/d02-x01-y01
Title=Dijet double-differential cross sections ($y^{*} < 0.5$), R=0.6
# END PLOT
# BEGIN PLOT /ATLAS_2014_I1268975/d02-x01-y02
Title=Dijet double-differential cross sections ($0.5 < y^{*} < 1.0$), R=0.6
# END PLOT
# BEGIN PLOT /ATLAS_2014_I1268975/d02-x01-y03
Title=Dijet double-differential cross sections ($1.0 < y^{*} < 1.5$), R=0.6
# END PLOT
# BEGIN PLOT /ATLAS_2014_I1268975/d02-x01-y04
Title=Dijet double-differential cross sections ($1.5 < y^{*} < 2.0$), R=0.6
# END PLOT
# BEGIN PLOT /ATLAS_2014_I1268975/d02-x01-y05
Title=Dijet double-differential cross sections ($2.0 < y^{*} < 2.5$), R=0.6
# END PLOT
# BEGIN PLOT /ATLAS_2014_I1268975/d02-x01-y06
Title=Dijet double-differential cross sections ($2.5 < y^{*} < 3.0$), R=0.6
# END PLOT