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
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
// -*- 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);

      // set ptcut from input option
      const double jetptcut = getOption<double>("PTJMIN", 20.0);
      _jetptcut = jetptcut * GeV;

      // set clustering radius from input option
      const double R = getOption<double>("R", 0.4);

      // set clustering algorithm from input option
      FastJets::Algo clusterAlgo;
      const string algoopt = getOption("ALGO", "ANTIKT");

      if ( algoopt == "KT" ) {
	clusterAlgo = FastJets::KT;
      } else if ( algoopt == "CA" ) {
	clusterAlgo = FastJets::CA;
      } else if ( algoopt == "ANTIKT" ) {
	clusterAlgo = FastJets::ANTIKT;
      } else {
	MSG_WARNING("Unknown jet clustering algorithm option " + algoopt + ". "
		    "Defaulting to anti-kT");
	clusterAlgo = FastJets::ANTIKT;
      }

      FastJets jetpro(vfs, clusterAlgo, R);
      declare(jetpro, "Jets");

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

      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
  RIVET_DECLARE_PLUGIN(MC_HHJETS);

}