rivet is hosted by Hepforge, IPPP Durham

Rivet analyses reference

MC_ZZKTSPLITTINGS

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

Monte Carlo validation observables for $Z[e^+ \, e^-]Z[\mu^+ \, \mu^-]$ + jets production

Source code: MC_ZZKTSPLITTINGS.cc
 1// -*- C++ -*-
 2#include "Rivet/Analyses/MC_JetSplittings.hh"
 3#include "Rivet/Projections/ZFinder.hh"
 4#include "Rivet/Projections/FastJets.hh"
 5#include "Rivet/Projections/VetoedFinalState.hh"
 6
 7namespace Rivet {
 8
 9
10  /// @brief MC validation analysis for Z[ee]Z[mumu] + jets events
11  class MC_ZZKTSPLITTINGS : public MC_JetSplittings {
12  public:
13
14    /// Default constructor
15    MC_ZZKTSPLITTINGS()
16      : MC_JetSplittings("MC_ZZKTSPLITTINGS", 4, "Jets")
17    {    }
18
19
20    /// @name Analysis methods
21    //@{
22
23    /// Book histograms
24    void init() {
25      Cut cut = Cuts::abseta < 3.5 && Cuts::pT > 25*GeV;
26      ZFinder zeefinder(FinalState(), cut, PID::ELECTRON, 65*GeV, 115*GeV,
27                        0.2, ZFinder::ClusterPhotons::NODECAY, ZFinder::AddPhotons::YES);
28      declare(zeefinder, "ZeeFinder");
29
30      VetoedFinalState zmminput;
31      zmminput.addVetoOnThisFinalState(zeefinder);
32      ZFinder zmmfinder(zmminput, cut, PID::MUON, 65*GeV, 115*GeV,
33                        0.2, ZFinder::ClusterPhotons::NODECAY, ZFinder::AddPhotons::YES);
34      declare(zmmfinder, "ZmmFinder");
35
36      VetoedFinalState jetinput;
37      jetinput
38          .addVetoOnThisFinalState(zeefinder)
39          .addVetoOnThisFinalState(zmmfinder);
40      FastJets jetpro(jetinput, FastJets::KT, 0.6);
41      declare(jetpro, "Jets");
42
43      MC_JetSplittings::init();
44    }
45
46
47    /// Do the analysis
48    void analyze(const Event & e) {
49      const ZFinder& zeefinder = apply<ZFinder>(e, "ZeeFinder");
50      if (zeefinder.bosons().size() != 1) vetoEvent;
51      const ZFinder& zmmfinder = apply<ZFinder>(e, "ZmmFinder");
52      if (zmmfinder.bosons().size() != 1) vetoEvent;
53      MC_JetSplittings::analyze(e);
54    }
55
56
57    /// Finalize
58    void finalize() {
59      MC_JetSplittings::finalize();
60    }
61
62    //@}
63
64  };
65
66
67
68  // The hook for the plugin system
69  RIVET_DECLARE_PLUGIN(MC_ZZKTSPLITTINGS);
70
71}