rivet is hosted by Hepforge, IPPP Durham

Rivet analyses reference

MC_HHJETS

Monte Carlo validation observables for $HH$ (stable) + jets production
Experiment: ()
Status: VALIDATED
Authors:
  • Andreas Papaefstathiou
No references listed
Beams: * *
Beam energies: ANY
Run details:
  • $HH$ production with stable Higgses

The available observables are the Higgs boson pair invariant mass, the separation between the two Higgs bosons, pT of the di-Higgs (any, hardest, second hardest), pT and pseudorapidity of Higgs bosons (any, hardest, second hardest), pT of jets 1--4, jet multiplicity, $\Delta\eta(h, \text{jet1})$, $\Delta R(\text{jet2}, \text{jet3})$, differential jet rates 0->1, 1->2, 2->3, 3->4, and integrated 0--4 jet rates.

Source code: MC_HHJETS.cc
  1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
// -*- C++ -*-
#include "Rivet/Analyses/MC_JetAnalysis.hh"
#include "Rivet/Projections/ZFinder.hh"
#include "Rivet/Projections/VetoedFinalState.hh"
#include "Rivet/Projections/IdentifiedFinalState.hh"
#include "Rivet/Projections/FastJets.hh"

namespace Rivet {

  /// @brief MC validation analysis for higgs pairs events (stable Higgses)
  class MC_HHJETS : public MC_JetAnalysis {
  public:

    /// Default constructor
    MC_HHJETS()
      : MC_JetAnalysis("MC_HHJETS", 4, "Jets")
    {    }


    /// @name Analysis methods
    //@{

    /// Book histograms
    void init() {
      IdentifiedFinalState ifs(Cuts::abseta < 10.0 && Cuts::pT > 0*GeV);
      ifs.acceptId(25);
      declare(ifs,"IFS");

      VetoedFinalState vfs;
      vfs.addVetoPairId(25);
      declare(FastJets(vfs, FastJets::ANTIKT, 0.4), "Jets");

      _h_HH_mass = bookHisto1D("HH_mass", 250, 240, 4000.0);
      _h_HH_dR = bookHisto1D("HH_dR", 25, 0.5, 10.0);
      _h_HH_dPhi = bookHisto1D("HH_dPhi", 64, 0, 3.2);
      _h_HH_deta= bookHisto1D("HH_deta", 50, -5, 5);
      _h_H_pT = bookHisto1D("H_pT", 50, 0, 2000.0);
      _h_HH_pT = bookHisto1D("HH_pT", 200, 0, 2000.0);
      _h_H_pT1 = bookHisto1D("H_pT1", 200, 0, 2000.0);
      _h_H_pT2 = bookHisto1D("H_pT2", 200, 0, 2000.0);
      _h_H_eta = bookHisto1D("H_eta", 50, -5.0, 5.0);
      _h_H_eta1 = bookHisto1D("H_eta1", 50, -5.0, 5.0);
      _h_H_eta2 = bookHisto1D("H_eta2", 50, -5.0, 5.0);
      _h_H_phi = bookHisto1D("H_phi", 25, 0.0, TWOPI);
      _h_H_jet1_deta = bookHisto1D("H_jet1_deta", 50, -5.0, 5.0);
      _h_H_jet1_dR = bookHisto1D("H_jet1_dR", 25, 0.5, 7.0);

      MC_JetAnalysis::init();
    }



    /// Do the analysis
    void analyze(const Event & e) {

      const IdentifiedFinalState& ifs = apply<IdentifiedFinalState>(e, "IFS");
      Particles allp = ifs.particlesByPt();
      if (allp.empty()) vetoEvent;

      const double weight = e.weight();

      FourMomentum hmom = allp[0].momentum();
      if (allp.size() > 1) {
        FourMomentum hmom2(allp[1].momentum());
        _h_HH_dR->fill(deltaR(hmom, hmom2), weight);
        _h_HH_dPhi->fill(deltaPhi(hmom, hmom2), weight);
        _h_HH_deta->fill(hmom.eta()-hmom2.eta(), weight);
        _h_HH_pT->fill((hmom+hmom2).pT(), weight);
        _h_HH_mass->fill((hmom+hmom2).mass(), weight);

        if (hmom.pT() > hmom2.pT()) {
          _h_H_pT1->fill(hmom.pT(), weight);
          _h_H_eta1->fill(hmom.eta(), weight);
          _h_H_pT2->fill(hmom2.pT(), weight);
          _h_H_eta2->fill(hmom2.eta(), weight);
        } else {
          _h_H_pT1->fill(hmom2.pT(), weight);
          _h_H_eta1->fill(hmom2.eta(), weight);
          _h_H_pT2->fill(hmom.pT(), weight);
          _h_H_eta2->fill(hmom.eta(), weight);
        }
      }
      _h_H_pT->fill(hmom.pT(), weight);
      _h_H_eta->fill(hmom.eta(), weight);
      _h_H_phi->fill(hmom.azimuthalAngle(), weight);


      // Get the jet candidates
      Jets jets = apply<FastJets>(e, "Jets").jetsByPt(20.0*GeV);
      if (!jets.empty()) {
        _h_H_jet1_deta->fill(deltaEta(hmom, jets[0]), weight);
        _h_H_jet1_dR->fill(deltaR(hmom, jets[0]), weight);
      }

      MC_JetAnalysis::analyze(e);
    }


    /// Finalize
    void finalize() {
      scale(_h_HH_mass, crossSection()/sumOfWeights());
      scale(_h_HH_dR, crossSection()/sumOfWeights());
      scale(_h_HH_deta, crossSection()/sumOfWeights());
      scale(_h_HH_dPhi, crossSection()/sumOfWeights());
      scale(_h_H_pT, crossSection()/sumOfWeights());
      scale(_h_H_pT1, crossSection()/sumOfWeights());
      scale(_h_H_pT2, crossSection()/sumOfWeights());
      scale(_h_HH_pT, crossSection()/sumOfWeights());
      scale(_h_H_eta, crossSection()/sumOfWeights());
      scale(_h_H_eta1, crossSection()/sumOfWeights());
      scale(_h_H_eta2, crossSection()/sumOfWeights());
      scale(_h_H_phi, crossSection()/sumOfWeights());
      scale(_h_H_jet1_deta, crossSection()/sumOfWeights());
      scale(_h_H_jet1_dR, crossSection()/sumOfWeights());

      MC_JetAnalysis::finalize();
    }

    //@}


  private:

    /// @name Histograms
    //@{
    Histo1DPtr _h_HH_mass, _h_HH_pT, _h_HH_dR, _h_HH_deta, _h_HH_dPhi;
    Histo1DPtr _h_H_pT, _h_H_pT1, _h_H_pT2, _h_H_eta, _h_H_eta1, _h_H_eta2, _h_H_phi;
    Histo1DPtr _h_H_jet1_deta, _h_H_jet1_dR;
    //@}

  };


  // The hook for the plugin system
  DECLARE_RIVET_PLUGIN(MC_HHJETS);

}