rivet is hosted by Hepforge, IPPP Durham

Rivet analyses reference

MC_HKTSPLITTINGS

Monte Carlo validation observables for $h[\tau^+ \, \tau^-]$ + jets production
Experiment: ()
Status: VALIDATED
Authors:
  • Frank Siegert
No references listed
Beams: * *
Beam energies: ANY
Run details:
  • $h [\to \tau^+ \tau^-]$ + jets.

Monte Carlo validation observables for $h[\tau^+ \, \tau^-]$ + jets production

Source code: MC_HKTSPLITTINGS.cc
 1// -*- C++ -*-
 2#include "Rivet/Analyses/MC_JetSplittings.hh"
 3#include "Rivet/Projections/ZFinder.hh"
 4#include "Rivet/Projections/FastJets.hh"
 5
 6namespace Rivet {
 7
 8
 9
10
11  /// @brief MC validation analysis for higgs [-> tau tau] + jets events
12  class MC_HKTSPLITTINGS : public MC_JetSplittings {
13  public:
14
15    /// Default constructor
16    MC_HKTSPLITTINGS()
17      : MC_JetSplittings("MC_HKTSPLITTINGS", 4, "Jets")
18    {    }
19
20
21    /// @name Analysis methods
22    //@{
23
24    /// Book histograms
25    void init() {
26      Cut cut = Cuts::abseta < 3.5 && Cuts::pT > 25*GeV;
27      /// @todo Urk, abuse! Need explicit HiggsFinder and TauFinder
28      ZFinder hfinder(FinalState(), cut, PID::TAU, 115*GeV, 135*GeV, 0.0, ZFinder::ClusterPhotons::NONE, ZFinder::AddPhotons::NO, 125*GeV);
29      declare(hfinder, "Hfinder");
30      FastJets jetpro(hfinder.remainingFinalState(), FastJets::KT, 0.6);
31      declare(jetpro, "Jets");
32
33      MC_JetSplittings::init();
34    }
35
36
37
38    /// Do the analysis
39    void analyze(const Event & e) {
40      const ZFinder& hfinder = apply<ZFinder>(e, "Hfinder");
41      if (hfinder.bosons().size() != 1) vetoEvent;
42      MC_JetSplittings::analyze(e);
43    }
44
45
46    /// Finalize
47    void finalize() {
48      MC_JetSplittings::finalize();
49    }
50
51    //@}
52
53  };
54
55
56
57  // The hook for the plugin system
58  RIVET_DECLARE_PLUGIN(MC_HKTSPLITTINGS);
59
60}