Rivet analyses referenceMC_HINCMonte Carlo validation observables for $h[\tau^+ \, \tau^-]$ productionExperiment: () Status: VALIDATED Authors:
Beams: * * Beam energies: ANY Run details:
Monte Carlo validation observables for $h[\tau^+ \, \tau^-]$ production Source code: MC_HINC.cc 1// -*- C++ -*-
2#include "Rivet/Analysis.hh"
3#include "Rivet/Projections/ZFinder.hh"
4
5namespace Rivet {
6
7
8 /// @brief MC validation analysis for higgs [-> tau tau] events
9 class MC_HINC : public Analysis {
10 public:
11
12 /// Default constructor
13 MC_HINC()
14 : Analysis("MC_HINC")
15 { }
16
17
18 /// @name Analysis methods
19 //@{
20
21 /// Book histograms
22 void init() {
23 Cut cut = Cuts::abseta < 3.5 && Cuts::pT > 25*GeV;
24 /// @todo Urk, abuse! Need explicit HiggsFinder and TauFinder?
25 ZFinder hfinder(FinalState(), cut, PID::TAU, 115*GeV, 135*GeV, 0.0, ZFinder::ClusterPhotons::NONE, ZFinder::AddPhotons::NO, 125*GeV);
26 declare(hfinder, "Hfinder");
27 book(_h_H_mass ,"H_mass", 50, 119.7, 120.3);
28 book(_h_H_pT ,"H_pT", logspace(100, 1.0, 0.5*(sqrtS()>0.?sqrtS():14000.)/GeV));
29 book(_h_H_pT_peak ,"H_pT_peak", 25, 0.0, 25.0);
30 book(_h_H_y ,"H_y", 40, -4, 4);
31 book(_h_H_phi ,"H_phi", 25, 0.0, TWOPI);
32 book(_h_lepton_pT ,"lepton_pT", logspace(100, 10.0, 0.25*(sqrtS()>0.?sqrtS():14000.)/GeV));
33 book(_h_lepton_eta ,"lepton_eta", 40, -4, 4);
34 }
35
36
37 /// Do the analysis
38 void analyze(const Event & e) {
39 const ZFinder& hfinder = apply<ZFinder>(e, "Hfinder");
40 if (hfinder.bosons().size() != 1) vetoEvent;
41 const double weight = 1.0;
42
43 FourMomentum hmom(hfinder.bosons()[0].momentum());
44 _h_H_mass->fill(hmom.mass()/GeV, weight);
45 _h_H_pT->fill(hmom.pT()/GeV, weight);
46 _h_H_pT_peak->fill(hmom.pT()/GeV, weight);
47 _h_H_y->fill(hmom.rapidity(), weight);
48 _h_H_phi->fill(hmom.phi(), weight);
49 for (const Particle& l : hfinder.constituents()) {
50 _h_lepton_pT->fill(l.pT()/GeV, weight);
51 _h_lepton_eta->fill(l.eta(), weight);
52 }
53 }
54
55
56 /// Finalize
57 void finalize() {
58 const double xsec = crossSection()/picobarn;
59 normalize(_h_H_mass, xsec);
60 normalize(_h_H_pT, xsec);
61 normalize(_h_H_pT_peak, xsec);
62 normalize(_h_H_y, xsec);
63 normalize(_h_H_phi, xsec);
64 normalize(_h_lepton_pT, xsec);
65 normalize(_h_lepton_eta, xsec);
66 }
67
68 //@}
69
70
71 private:
72
73 /// @name Histograms
74 //@{
75 Histo1DPtr _h_H_mass;
76 Histo1DPtr _h_H_pT;
77 Histo1DPtr _h_H_pT_peak;
78 Histo1DPtr _h_H_y;
79 Histo1DPtr _h_H_phi;
80 Histo1DPtr _h_lepton_pT;
81 Histo1DPtr _h_lepton_eta;
82 //@}
83
84 };
85
86
87
88 // The hook for the plugin system
89 RIVET_DECLARE_PLUGIN(MC_HINC);
90
91}
|