Rivet analyses

title: D0_2008_I779574

Measurement of D0 Run II differential jet cross sections

Experiment: D0 (Tevatron Run 2)

Inspire ID: 779574

Status: VALIDATED

Authors: - Andy Buckley - Gavin Hesketh

References: - Phys.Rev.Lett.101:062001,2008 - DOI: 10.1103/PhysRevLett.101.062001 - arXiv: 0802.2400

Beams: p- p+

Beam energies: (980.0, 980.0)GeV

Run details: - QCD events at $\sqrt{s} = 1960~\GeV$. A $\pT^\text{min} cut is probably necessary since the lowest jet $\pT$ bin is at 50~GeV

Measurement of the inclusive jet cross section in $p \bar{p}$ collisions at center-of-mass energy $\sqrt{s} = 1.96~\TeV$. The data cover jet transverse momenta from 50--600~GeV and jet rapidities in the range -2.4 to 2.4.

Source code:D0_2008_I779574.cc

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

include "Rivet/Analysis.hh"

include "Rivet/Projections/FastJets.hh"

include "Rivet/Projections/FinalState.hh"

include "Rivet/Projections/LeadingParticlesFinalState.hh"

include "Rivet/Projections/VetoedFinalState.hh"

namespace Rivet {

/// @brief D0 differential jet cross sections /// /// @author Andy Buckley /// @author Gavin Hesketh class D0_2008_I779574 : public Analysis { public:

RIVET_DEFAULT_ANALYSIS_CTOR(D0_2008_I779574);


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

void init() {
  // Full final state
  FinalState fs;
  declare(fs, "FS");

  // Jets
  FastJets jetpro(fs, JetAlg::D0ILCONE, 0.7);
  declare(jetpro, "Jets");

  // Book histograms
  book(_h_dsigdptdy_y00_04, 1, 1, 1);
  book(_h_dsigdptdy_y04_08, 2, 1, 1);
  book(_h_dsigdptdy_y08_12, 3, 1, 1);
  book(_h_dsigdptdy_y12_16, 4, 1, 1);
  book(_h_dsigdptdy_y16_20, 5, 1, 1);
  book(_h_dsigdptdy_y20_24, 6, 1, 1);
}


/// Do the analysis
void analyze(const Event& event) {
  // Skip if the event is empty
  const FinalState& fs = apply<FinalState>(event, "FS");
  if (fs.empty()) {
    MSG_DEBUG("Empty event!");
    vetoEvent;
  }

  // Find the jets
  const JetFinder& jetpro = apply<JetFinder>(event, "Jets");
  // Fill histo for each jet
  for (const Jet& j : jetpro.jets(Cuts::pT > 50 * GeV)) {
    const double pt = j.pT();
    const double y = j.absrap();
    MSG_TRACE("Filling histos: pT = " << pt / GeV << ", |y| = " << y);
    if (y < 0.4) {
      _h_dsigdptdy_y00_04->fill(pt / GeV);
    }
    else if (y < 0.8) {
      _h_dsigdptdy_y04_08->fill(pt / GeV);
    }
    else if (y < 1.2) {
      _h_dsigdptdy_y08_12->fill(pt / GeV);
    }
    else if (y < 1.6) {
      _h_dsigdptdy_y12_16->fill(pt / GeV);
    }
    else if (y < 2.0) {
      _h_dsigdptdy_y16_20->fill(pt / GeV);
    }
    else if (y < 2.4) {
      _h_dsigdptdy_y20_24->fill(pt / GeV);
    }
  }
}


/// Finalize
void finalize() {
  /// Scale by cross-section
  const double scalefactor = crossSection() / picobarn / sumOfWeights();
  scale(_h_dsigdptdy_y00_04, scalefactor);
  scale(_h_dsigdptdy_y04_08, scalefactor);
  scale(_h_dsigdptdy_y08_12, scalefactor);
  scale(_h_dsigdptdy_y12_16, scalefactor);
  scale(_h_dsigdptdy_y16_20, scalefactor);
  scale(_h_dsigdptdy_y20_24, scalefactor);
}

/// @}

private:

/// @name Histograms
/// @{
Histo1DPtr _h_dsigdptdy_y00_04;
Histo1DPtr _h_dsigdptdy_y04_08;
Histo1DPtr _h_dsigdptdy_y08_12;
Histo1DPtr _h_dsigdptdy_y12_16;
Histo1DPtr _h_dsigdptdy_y16_20;
Histo1DPtr _h_dsigdptdy_y20_24;
/// @}

};

RIVET_DECLARE_ALIASED_PLUGIN(D0_2008_I779574, D0_2008_S7662670);

} ```

Aliases: - D0_2008_S7662670