Rivet analyses

High-mass Drell-Yan at 8 TeV

Experiment: ATLAS (LHC)

Inspire ID: 1467454

Status: VALIDATED

Authors: - Markus Zinser

References: - Expt page: ATLAS-STDM-2014-06 - JHEP 1608 (2016) 009 - DOI: 10.1007/JHEP08(2016)009 - arXiv: 1606.01736

Beams: p+ p+

Beam energies: (4000.0, 4000.0)GeV

Run details: - Drell-Yan production in pp collisions at 8 TeV, leptonic Z decays

This paper presents a measurement of the double-differential cross section for the Drell-Yan Z/γ^* → ℓ⁺ℓ⁻ and photon-induced γγ → ℓ⁺ℓ⁻ processes where ℓ is an electron or muon. The measurement is performed for invariant masses of the lepton pairs, m_ℓℓ, between 116 GeV and 1500 GeV using a sample of 20.3 fb⁻¹ of pp collisions data at centre-of-mass energy of $\sqrt{s} = 8$ TeV collected by the ATLAS detector at the LHC in 2012. The data are presented double differentially in invariant mass and absolute dilepton rapidity as well as in invariant mass and absolute pseudorapidity separation of the lepton pair. The single-differential cross section as a function of m_ℓℓ is also reported. The electron and muon channel measurements are combined and a total experimental precision of better than 1% is achieved at low m_ℓℓ. A comparison to next-to-next-to-leading order perturbative QCD predictions using several recent parton distribution functions and including next-to-leading order electroweak effects indicates the potential of the data to constrain parton distribution functions. In particular, a large impact of the data on the photon PDF is demonstrated. Decay channel (e or mu) is selected by LMODE=EL,MU

Source code:`ATLAS_2016_I1467454.cc`

// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/DileptonFinder.hh"

namespace Rivet {


  /// @brief High-mass Drell-Yan at 8 TeV
  class ATLAS_2016_I1467454 : public Analysis {
  public:

    /// Constructor
    RIVET_DEFAULT_ANALYSIS_CTOR(ATLAS_2016_I1467454);
    /// @}


    /// @name Analysis methods
    /// @{

    /// Book histograms and initialise projections before the run
    void init() {

      // Get options from the new option system
      _mode = 0;
      if ( getOption("LMODE") == "EL" ) _mode = 0;
      if ( getOption("LMODE") == "MU" ) _mode = 1;

      Cut cuts = Cuts::abseta < 2.5 && Cuts::pT > 30*GeV;
      DileptonFinder zfinder(91.2*GeV, 0.1, cuts && Cuts::abspid == (_mode ? PID::MUON : PID::ELECTRON), Cuts::massIn(116*GeV, 1500*GeV));
      declare(zfinder, "DileptonFinder");

      size_t ch = _mode? 11 : 0; // offset
      book(_hist_mll, 18 + ch, 1, 1);

      const vector<double> mll_bins{ 116., 150., 200., 300., 500., 1500. };
      book(_hist_rap, mll_bins);
      book(_hist_deta, mll_bins);
      for (size_t i=0; i < _hist_rap->numBins(); ++i) {
        book(_hist_rap->bin(i+1),  19 + ch + i, 1, 1);
        book(_hist_deta->bin(i+1), 24 + ch + i, 1, 1);
      }

    }


    /// Perform the per-event analysis
    void analyze(const Event& event) {

      const DileptonFinder& zfinder = apply<DileptonFinder>(event, "DileptonFinder");
      if (zfinder.bosons().size() != 1)  vetoEvent;

      const Particle z0  = zfinder.bosons()[0];
      /// @todo Could use z0.constituents()
      const Particle el1 = zfinder.leptons()[0];
      const Particle el2 = zfinder.leptons()[1];

      if (el1.pT() > 40*GeV || el2.pT() > 40*GeV) {
        const double mass = z0.mass();
        _hist_mll->fill(mass/GeV);
        _hist_rap->fill(mass/GeV, z0.absrap());
        _hist_deta->fill(mass/GeV, deltaEta(el1,el2));
      }
    }


    /// Normalise histograms etc., after the run
    void finalize() {

      const double sf = crossSection()/picobarn/sumOfWeights();
      scale(_hist_mll, sf);
      scale(_hist_rap, sf*0.5);
      divByGroupWidth(_hist_rap);
      scale(_hist_deta, sf*0.5);
      divByGroupWidth(_hist_deta);

    }

    /// @}


    /// Choose to work in electron or muon mode
    size_t _mode;


    /// @name Histograms
    /// @{
    Histo1DPtr _hist_mll;
    Histo1DGroupPtr _hist_rap, _hist_deta;
    /// @}

  };

  RIVET_DECLARE_PLUGIN(ATLAS_2016_I1467454);

}

Rivet analyses

High-mass Drell-Yan at 8 TeV

Source code:ATLAS_2016_I1467454.cc

Source code:`ATLAS_2016_I1467454.cc`