1// -*- C++ -*-
 2#include "Rivet/Analysis.hh"
 3#include "Rivet/Projections/DileptonFinder.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      // set FS cuts from input options
24      const double etacut = getOption<double>("ABSETATAUMAX", 3.5);
25      const double ptcut = getOption<double>("PTTAUMIN", 25.);
26
27      Cut cut = Cuts::abseta < etacut && Cuts::pT > ptcut*GeV;
28      /// @todo Er, FS taus?
29      DileptonFinder hfinder(125*GeV, 0.0, cut && Cuts::abspid == PID::TAU, Cuts::massIn(115*GeV, 135*GeV));
30
31      declare(hfinder, "Hfinder");
32      book(_h_H_mass ,"H_mass", 50, 119.7, 120.3);
33      book(_h_H_pT ,"H_pT", logspace(100, 1.0, 0.5*(sqrtS()>0.?sqrtS():14000.)/GeV));
34      book(_h_H_pT_peak ,"H_pT_peak", 25, 0.0, 25.0);
35      book(_h_H_y ,"H_y", 40, -4, 4);
36      book(_h_H_phi ,"H_phi", 25, 0.0, TWOPI);
37      book(_h_lepton_pT ,"lepton_pT", logspace(100, 10.0, 0.25*(sqrtS()>0.?sqrtS():14000.)/GeV));
38      book(_h_lepton_eta ,"lepton_eta", 40, -4, 4);
39    }
40
41
42    /// Do the analysis
43    void analyze(const Event & e) {
44      const DileptonFinder& hfinder = apply<DileptonFinder>(e, "Hfinder");
45      if (hfinder.bosons().size() != 1) vetoEvent;
46
47      FourMomentum hmom(hfinder.bosons()[0].momentum());
48      _h_H_mass->fill(hmom.mass()/GeV);
49      _h_H_pT->fill(hmom.pT()/GeV);
50      _h_H_pT_peak->fill(hmom.pT()/GeV);
51      _h_H_y->fill(hmom.rapidity());
52      _h_H_phi->fill(hmom.phi());
53      for (const Particle& l : hfinder.constituents()) {
54        _h_lepton_pT->fill(l.pT()/GeV);
55        _h_lepton_eta->fill(l.eta());
56      }
57    }
58
59
60    /// Finalize
61    void finalize() {
62      const double xsec = crossSection()/picobarn;
63      normalize(_h_H_mass, xsec);
64      normalize(_h_H_pT, xsec);
65      normalize(_h_H_pT_peak, xsec);
66      normalize(_h_H_y, xsec);
67      normalize(_h_H_phi, xsec);
68      normalize(_h_lepton_pT, xsec);
69      normalize(_h_lepton_eta, xsec);
70    }
71
72    /// @}
73
74
75  private:
76
77    /// @name Histograms
78    /// @{
79    Histo1DPtr _h_H_mass;
80    Histo1DPtr _h_H_pT;
81    Histo1DPtr _h_H_pT_peak;
82    Histo1DPtr _h_H_y;
83    Histo1DPtr _h_H_phi;
84    Histo1DPtr _h_lepton_pT;
85    Histo1DPtr _h_lepton_eta;
86    /// @}
87
88  };
89
90
91
92  RIVET_DECLARE_PLUGIN(MC_HINC);
93
94}