rivet is hosted by Hepforge, IPPP Durham

Rivet analyses reference

MC_ZINC

Monte Carlo validation observables for $Z[\ell^+ \, \ell^-]$ production
Experiment: ()
Status: VALIDATED
Authors:
  • Frank Siegert
  • Christian Gutschow
No references listed
Beams: * *
Beam energies: ANY
Run details:
  • $\ell^+ \ell^-$ analysis. Needs mass cut on lepton pair to avoid photon singularity, e.g. a min range of $66 < m_{ll} < 116$ GeV

Monte Carlo validation observables for $Z[\ell^+ \, \ell^-]$ production

Source code: MC_ZINC.cc
  1// -*- C++ -*-
  2#include "Rivet/Analysis.hh"
  3#include "Rivet/Projections/ZFinder.hh"
  4
  5namespace Rivet {
  6
  7
  8
  9  /// @brief MC validation analysis for Z events
 10  class MC_ZINC : public Analysis {
 11  public:
 12
 13    /// Default constructor
 14    RIVET_DEFAULT_ANALYSIS_CTOR(MC_ZINC);
 15
 16    /// @name Analysis methods
 17    //@{
 18
 19    /// Book histograms
 20    void init() {
 21		  _dR=0.2;
 22      if (getOption("SCHEME") == "BARE")  _dR = 0.0;
 23		  _lepton=PID::ELECTRON;
 24      if (getOption("LMODE") == "MU")  _lepton = PID::MUON;
 25
 26      FinalState fs;
 27      Cut cut = Cuts::abseta < 3.5 && Cuts::pT > 25*GeV;
 28      ZFinder zfinder(fs, cut, _lepton, 65.0*GeV, 115.0*GeV, _dR, ZFinder::ClusterPhotons::NODECAY, ZFinder::AddPhotons::YES);
 29      declare(zfinder, "ZFinder");
 30
 31      book(_h_Z_mass ,"Z_mass", 50, 66.0, 116.0);
 32      book(_h_Z_pT ,"Z_pT", logspace(100, 1.0, 0.5*(sqrtS()>0.?sqrtS():14000.)/GeV));
 33      book(_h_Z_pT_peak ,"Z_pT_peak", 25, 0.0, 25.0);
 34      book(_h_Z_y ,"Z_y", 40, -4.0, 4.0);
 35      book(_h_Z_phi ,"Z_phi", 25, 0.0, TWOPI);
 36      book(_h_lepton_pT ,"lepton_pT", logspace(100, 10.0, 0.25*(sqrtS()>0.?sqrtS():14000.)/GeV));
 37      book(_h_lepton_eta ,"lepton_eta", 40, -4.0, 4.0);
 38
 39    }
 40
 41
 42
 43    /// Do the analysis
 44    void analyze(const Event & e) {
 45      const ZFinder& zfinder = apply<ZFinder>(e, "ZFinder");
 46      if (zfinder.bosons().size() != 1) vetoEvent;
 47
 48      FourMomentum zmom(zfinder.bosons()[0].momentum());
 49      _h_Z_mass->fill(zmom.mass()/GeV);
 50      _h_Z_pT->fill(zmom.pT()/GeV);
 51      _h_Z_pT_peak->fill(zmom.pT()/GeV);
 52      _h_Z_y->fill(zmom.rapidity());
 53      _h_Z_phi->fill(zmom.phi());
 54      for (const Particle& l : zfinder.constituents()) {
 55        _h_lepton_pT->fill(l.pT()/GeV);
 56        _h_lepton_eta->fill(l.eta());
 57      }
 58    }
 59
 60
 61    /// Finalize
 62    void finalize() {
 63      const double s = crossSection()/picobarn/sumOfWeights();
 64      scale(_h_Z_mass, s);
 65      scale(_h_Z_pT, s);
 66      scale(_h_Z_pT_peak, s);
 67      scale(_h_Z_y, s);
 68      scale(_h_Z_phi, s);
 69      scale(_h_lepton_pT, s);
 70      scale(_h_lepton_eta, s);
 71    }
 72
 73    //@}
 74
 75
 76  protected:
 77
 78    /// @name Parameters for specialised e/mu and dressed/bare subclassing
 79    //@{
 80    double _dR;
 81    PdgId _lepton;
 82    //@}
 83
 84
 85  private:
 86
 87    /// @name Histograms
 88    //@{
 89    Histo1DPtr _h_Z_mass;
 90    Histo1DPtr _h_Z_pT;
 91    Histo1DPtr _h_Z_pT_peak;
 92    Histo1DPtr _h_Z_y;
 93    Histo1DPtr _h_Z_phi;
 94    Histo1DPtr _h_lepton_pT;
 95    Histo1DPtr _h_lepton_eta;
 96    //@}
 97
 98  };
 99
100  // The hooks for the plugin system
101  RIVET_DECLARE_PLUGIN(MC_ZINC);
102}