Commit 69983f04 authored by Christian Gutschow's avatar Christian Gutschow Committed by Andy Buckley

ATLAS_2016_I1487726: collinear Wj at 8 TeV

parent d39f59ca
// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/DressedLeptons.hh"
#include "Rivet/Projections/FastJets.hh"
namespace Rivet {
class ATLAS_2016_I1487726 : public Analysis {
public:
/// Constructor
/// @brief Collinear W emissions at 8 TeV
DEFAULT_RIVET_ANALYSIS_CTOR(ATLAS_2016_I1487726);
public:
/// @name Analysis methods
//@{
/// Book histograms and initialise projections before the run
void init() {
_mode = 0;
if ( getOption("LMODE") == "EL" ) _mode = 1;
// These really should include non-prompt leptons/photons
FinalState mufs(Cuts::abspid == PID::MUON);
FinalState elfs(Cuts::abspid == PID::ELECTRON);
FinalState phs(Cuts::abspid == PID::PHOTON);
Cut lep_fid = (Cuts::abseta < 2.4 && Cuts::pT >= 25*GeV);
DressedLeptons dlep(phs, _mode? elfs : mufs, 0.1, lep_fid, true);
declare(dlep, "DressedLeptons");
FastJets fj(FinalState(), FastJets::ANTIKT, 0.4, JetAlg::Muons::NONE, JetAlg::Invisibles::NONE);
declare(fj, "AntiKt4Jets");
book(h_mu_jet_dr, 2, 1, 1);
book(h_mu_jet_dr_pt500600, 4, 1, 1);
book(h_mu_jet_dr_pt650, 5, 1, 1);
}
/// Perform the per-event analysis
void analyze(const Event& event) {
const vector<DressedLepton> leptons = apply<DressedLeptons>(event, "DressedLeptons").dressedLeptons();
const Jets jets = apply<FastJets>(event, "AntiKt4Jets").jetsByPt(Cuts::pT >= 100*GeV && Cuts::abseta <= 2.1);
if (leptons.size() != 1) vetoEvent;
if (jets.size() < 1) vetoEvent;
if (jets[0].pt() < 500.0*GeV) vetoEvent;
// find closest jet to the lepton.
Jet jet;
double drmin = 999;
for (const Jet &j : jets) {
double dr = deltaR(leptons[0], j);
if (dr < drmin) {
drmin = dr;
jet = j;
}
}
h_mu_jet_dr->fill(drmin);
if (jets[0].pT() > 650*GeV) h_mu_jet_dr_pt650->fill(drmin);
else if (jets[0].pT() > 500*GeV && jets[0].pT() < 600*GeV) {
h_mu_jet_dr_pt500600->fill(drmin);
}
}
/// Normalise histograms etc., after the run
void finalize() {
const double sf = crossSection() / femtobarn / sumOfWeights();
scale(h_mu_jet_dr, sf);
scale(h_mu_jet_dr_pt500600, sf);
scale(h_mu_jet_dr_pt650, sf);
}
//@}
protected:
size_t _mode;
private:
/// @name Histograms
//@{
Histo1DPtr h_mu_jet_dr;
Histo1DPtr h_mu_jet_dr_pt500600;
Histo1DPtr h_mu_jet_dr_pt650;
//@}
};
// The hook for the plugin system
DECLARE_RIVET_PLUGIN(ATLAS_2016_I1487726);
}
Name: ATLAS_2016_I1487726
Year: 2016
Summary: Collinear W emissions at 8 TeV
Experiment: ATLAS
Collider: LHC
InspireID: 1487726
Status: VALIDATED
Authors:
- marek.schoenherr@cern.ch
- mileswu@uchicago.edu
- chris.g@cern.ch
References:
- DOI:10.1016/j.physletb.2016.12.005
- arXiv:1609.07045 [hep-ex]
- Phys.Lett. B765 (2017) 132-153
Keywords:
- WJETS
- COLLINEAR
- VJETS
RunInfo: 8 TeV $pp \to W+\text{jets} \to \mu\nu+\text{jets}$.
Beams: [p+, p+]
Energies: [8000]
Luminosity_fb: 20.3
Options:
- LMODE=EL
Description:
'The $W$ boson angular distribution in events with high transverse momentum jets is measured using data collected by the ATLAS
experiment from proton-proton collisions at a centre-of-mass energy $\sqrt{s}=8$ TeV at the Large Hadron Collider, corresponding to an
integrated luminosity of 20.3 fb$^{-1}$. The focus is on the contributions to $W+$jets processes from real $W$ emission,
which is achieved by studying events where a muon is observed close to a high transverse momentum jet. At small angular separations,
these contributions are expected to be large. Various theoretical models of this process are compared to the data in terms of the
absolute cross-section and the angular distributions of the muon from the leptonic $W$ decay.'
BibKey: Aaboud:2016ylh
BibTeX: '@article{Aaboud:2016ylh,
author = "Aaboud, Morad and others",
title = "{Measurement of $W$ boson angular distributions in events
with high transverse momentum jets at $\sqrt{s}=$ 8 TeV
using the ATLAS detector}",
collaboration = "ATLAS",
journal = "Phys. Lett.",
volume = "B765",
year = "2017",
pages = "132-153",
doi = "10.1016/j.physletb.2016.12.005",
eprint = "1609.07045",
archivePrefix = "arXiv",
primaryClass = "hep-ex",
reportNumber = "CERN-EP-2016-182",
SLACcitation = "%%CITATION = ARXIV:1609.07045;%%"
}'
ReleaseTests:
- $A LHC-8-Wmv-boosted
- $A-2 LHC-8-Wev-boosted :LMODE=EL
# BEGIN PLOT /ATLAS_2016_I1487726/d02-x01-y01
LogY=0
LegendXPos=0.05
Title=$\Delta R$ between lepton and closest jet, $p_\mathrm{T}^\text{lead}>500\,\text{GeV}$
XLabel=$\Delta R(\mu,j)$
YLabel=d$\sigma$/d$\Delta R$ [fb]
# END PLOT
# BEGIN PLOT /ATLAS_2016_I1487726/d04-x01-y01
LogY=0
LegendXPos=0.05
Title=$\Delta R$ between lepton and closest jet, $500\,\text{GeV}<p_\mathrm{T}^\text{lead}<600\,\text{GeV}$
XLabel=$\Delta R(\mu,j)$
YLabel=d$\sigma$/d$\Delta R$ [fb]
# END PLOT
# BEGIN PLOT /ATLAS_2016_I1487726/d05-x01-y01
LogY=0
Title=$\Delta R$ between lepton and closest jet, $p_\mathrm{T}^\text{lead}>650\,\text{GeV}$
XLabel=$\Delta R(\mu,j)$
YLabel=d$\sigma$/d$\Delta R$ [fb]
LegendXPos=0.1
LegendYPos=0.25
# END PLOT
This diff is collapsed.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment