Rivet analyses

title: L3_1990_I298078

Jets Rates using the Jade Algorithm at 91 GeV

Experiment: L3 (LEP)

Inspire ID: 298078

Status: VALIDATED

Authors: - Peter Richardson

References: - Phys.Lett.B 248 (1990) 464-472, 1990.

Beams: e+ e-

Beam energies: (45.6, 45.6)GeV

Run details: - e+e- to hadrons

Measurement of the integrated jet rates using the Jade algorithm by the L3 experiment at LEP1.

Source code:L3_1990_I298078.cc

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

include "Rivet/Analysis.hh"

include "Rivet/Projections/FastJets.hh"

include "Rivet/Projections/FinalState.hh"

namespace Rivet {

/// @brief Add a short analysis description here class L3_1990_I298078 : public Analysis { public:

/// Constructor
RIVET_DEFAULT_ANALYSIS_CTOR(L3_1990_I298078);


/// @name Analysis methods
/// @{

/// Book histograms and initialise projections before the run
void init() {

  // Projections to use
  const FinalState FS;
  declare(FS, "FS");
  FastJets jadeJets = FastJets(FS, JetAlg::JADE, -1, JetMuons::ALL, JetInvisibles::DECAY);
  declare(jadeJets, "JadeJets");
  // book histos
  book(_h_y_2_JADE, 1, 1, 1);
  book(_h_y_3_JADE, 1, 1, 2);
  book(_h_y_4_JADE, 1, 1, 3);
}


/// Perform the per-event analysis
void analyze(const Event& event) {
  const FastJets& jadejet = apply<FastJets>(event, "JadeJets");
  if (jadejet.clusterSeq()) {
    const double y_23 = jadejet.clusterSeq()->exclusive_ymerge_max(2);
    const double y_34 = jadejet.clusterSeq()->exclusive_ymerge_max(3);
    const double y_45 = jadejet.clusterSeq()->exclusive_ymerge_max(4);
    for (const auto& b : _h_y_2_JADE->bins()) {
      const double ycut = std::stod(b.xEdge());
      if (y_23 < ycut) _h_y_2_JADE->fill(b.xEdge());
    }
    for (const auto& b : _h_y_3_JADE->bins()) {
      const double ycut = std::stod(b.xEdge());
      if (y_34 < ycut && y_23 > ycut) {
        _h_y_3_JADE->fill(b.xEdge());
      }
    }
    for (const auto& b : _h_y_4_JADE->bins()) {
      const double ycut = std::stod(b.xEdge());
      if (y_45 < ycut && y_34 > ycut) {
        _h_y_4_JADE->fill(b.xEdge());
      }
    }
  }
}


/// Normalise histograms etc., after the run
void finalize() {
  scale(_h_y_2_JADE, 100. / sumOfWeights());
  scale(_h_y_3_JADE, 100. / sumOfWeights());
  scale(_h_y_4_JADE, 100. / sumOfWeights());
}

/// @}


/// @name Histograms
/// @{
BinnedHistoPtr<string> _h_y_2_JADE, _h_y_3_JADE, _h_y_4_JADE;
/// @}

};

RIVET_DECLARE_PLUGIN(L3_1990_I298078);

} ```