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}