Rivet analyses referenceATLAS_2013_I1234228High-mass Drell-Yan at 7 TeVExperiment: ATLAS (LHC) Inspire ID: 1234228 Status: VALIDATED Authors:
Beam energies: (3500.0, 3500.0) GeV Run details:
This Letter reports a measurement of the high-mass Drell-Yan differential cross-section in proton-proton collisions at a centre-of-mass energy of 7 TeV at the LHC. Based on an integrated luminosity of 4.9 fb$^{-1}$, the differential cross-section in the $Z/\gamma^\ast \rightarrow e^+e^-$ channel is measured with the ATLAS detector as a function of the invariant mass, $m_{ee}$, in the range $116< m_{ee} <1500$ GeV, for a fiducial region in which both the electron and the positron have transverse momentum $p_\text{T}>25$ GeV and pseudorapidity $|\eta|<2.5$. A comparison is made to various event generators and to the predictions of perturbative QCD calculations at next-to-next-to-leading order. Source code: ATLAS_2013_I1234228.cc 1// -*- C++ -*-
2#include "Rivet/Analysis.hh"
3#include "Rivet/Projections/FinalState.hh"
4#include "Rivet/Projections/DileptonFinder.hh"
5
6namespace Rivet {
7
8
9 /// @brief Add a short analysis description here
10 class ATLAS_2013_I1234228 : public Analysis {
11 public:
12
13
14 /// Constructor
15 RIVET_DEFAULT_ANALYSIS_CTOR(ATLAS_2013_I1234228);
16
17 /// @name Analysis methods
18 /// @{
19
20 /// Book histograms and initialise projections before the run
21 void init() {
22
23 const FinalState fs;
24 Cut cuts = Cuts::abseta < 2.5 && Cuts::pT > 25*GeV;
25 DileptonFinder zfinder(91.2*GeV, 0.1, cuts && Cuts::abspid == PID::ELECTRON, Cuts::massIn(116*GeV, 1500*GeV));
26 declare(zfinder, "DileptonFinder");
27
28 book(_hist_mll, 1, 1, 2);
29 }
30
31
32 /// Perform the per-event analysis
33 void analyze(const Event& event) {
34 const DileptonFinder& zfinder = apply<DileptonFinder>(event, "DileptonFinder");
35
36 if (zfinder.bosons().size() != 1) vetoEvent;
37
38 double mass = zfinder.bosons()[0].mass();
39 _hist_mll->fill(mass);
40 }
41
42
43 /// Normalise histograms etc., after the run
44 void finalize() {
45 const double sf = crossSection()/picobarn/sumOfWeights();
46 scale(_hist_mll, sf);
47 }
48
49 /// @}
50
51 private:
52
53
54 /// @name Histograms
55 /// @{
56 Histo1DPtr _hist_mll;
57 /// @}
58
59 };
60
61 RIVET_DECLARE_PLUGIN(ATLAS_2013_I1234228);
62
63}
|