Rivet analyses

title: MC_REENTRANT

Experiment: ()

Status: UNVALIDATED

Authors: none listed

References: none listed

Beams: * *

Beam energies: ANY

Run details: none listed

Source code:MC_REENTRANT.cc

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

include "Rivet/Analysis.hh"

include "Rivet/Projections/ChargedFinalState.hh"

include "Rivet/Projections/FinalState.hh"

namespace Rivet {

/// Generic analysis looking att pp pseudorepidity distributions at /// two collision energies. It usess the possibility to read in a /// pre-exixsting yoda file, and if histograms for both energies are /// filled when finalize() is called, a ratio plot is produced. class MC_REENTRANT : public Analysis { public:

/// Constructor
RIVET_DEFAULT_ANALYSIS_CTOR(MC_REENTRANT);


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

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

  // Projections
  const FinalState fs(Cuts::abseta < 5 && Cuts::pT > 500 * MeV);
  declare(fs, "FS");
  declare(ChargedFinalState(fs), "CFS");

  // Histograms. Booked for both 900 GeV and 7 TeV and their ratio.
  book(_histEta70, "Eta70", 50, -5, 5);
  book(_histEta09, "Eta09", 50, -5, 5);
  book(_histEtaR, "EtaR", 50, -5, 5);

  if (isCompatibleWithSqrtS(900.))
    fill09 = true;
  else if (isCompatibleWithSqrtS(7000.))
    fill70 = true;
}


/// Perform the per-event analysis
void analyze(const Event& event) {
  const FinalState& cfs = apply<FinalState>(event, "CFS");
  for (const Particle& p : cfs.particles()) {
    if (fill09)
      _histEta09->fill(p.eta());
    else if (fill70)
      _histEta70->fill(p.eta());
  }
}


/// Finalize
void finalize() {
  if (fill70) scale(_histEta70, 1.0 / sumOfWeights());
  if (fill09) scale(_histEta09, 1.0 / sumOfWeights());
  if (_histEta70->numEntries() > 0 && _histEta09->numEntries() > 0) {
    divide(_histEta70, _histEta09, _histEtaR);
  }
}

/// @}

private:

/// @name Histograms
/// @{
Histo1DPtr _histEta09, _histEta70;
Estimate1DPtr _histEtaR;
/// @}

bool fill09 = false, fill70 = false;

};

RIVET_DECLARE_PLUGIN(MC_REENTRANT);

} ```