Rivet analyses

title: MC_JETS

Monte Carlo validation observables for jet production

Experiment: ()

Status: VALIDATED

Authors: - Frank Siegert

References: none listed

Beams: * *

Beam energies: ANY

Run details: - Pure QCD jet production events at an arbitrary collider.

Jets with $p_\perp>20$ GeV are constructed with an anti-$k_\perp$ jet finder with $R=0.4$ and projected onto many different observables.

Source code:MC_JETS.cc

```c++ // -- C++ --

include "Rivet/Analyses/MC_JETS_BASE.hh"

include "Rivet/Projections/FastJets.hh"

include "Rivet/Projections/FinalState.hh"

include "fastjet/contrib/SoftDrop.hh"

namespace Rivet {

/// @brief MC validation analysis for jet events class MC_JETS : public MC_JETS_BASE { public:

MC_JETS()
    : MC_JETS_BASE("MC_JETS", 4, "Jets") { }


void init() {
  // 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
  JetAlg clusterAlgo;
  const string algoopt = getOption("ALGO", "ANTIKT");
  if (algoopt == "KT") {
    clusterAlgo = JetAlg::KT;
  }
  else if (algoopt == "CA") {
    clusterAlgo = JetAlg::CA;
  }
  else if (algoopt == "ANTIKT") {
    clusterAlgo = JetAlg::ANTIKT;
  }
  else {
    MSG_WARNING("Unknown jet clustering algorithm option " + algoopt + ". Defaulting to anti-kT");
    clusterAlgo = JetAlg::ANTIKT;
  }

  FinalState fs;
  FastJets jetpro(fs, clusterAlgo, R);

  const string groomopt = getOption("GROOM", "");

  if (groomopt == "SD") {
    jetpro.addTrf(new fastjet::contrib::SoftDrop(0.0, 0.1));
  }
  else if (groomopt == "TRIM") {
    jetpro.addTrf(new fastjet::Filter(fastjet::JetDefinition(fastjet::kt_algorithm, 0.2),
                                      fastjet::SelectorPtFractionMin(0.05)));
  }
  else if (groomopt != "") {
    MSG_WARNING("Unknown GROOM=" + groomopt + " option. Not applying jet grooming");
  }

  declare(jetpro, "Jets");
  MC_JETS_BASE::init();
}


void analyze(const Event& event) {
  MC_JETS_BASE::analyze(event);
}


void finalize() {
  MC_JETS_BASE::finalize();
}

};

RIVET_DECLARE_PLUGIN(MC_JETS);

} ```