Rivet analyses

title: D0_2009_I824127

Dijet angular distributions

Experiment: D0 (Tevatron Run 2)

Inspire ID: 824127

Status: VALIDATED

Authors: - Frank Siegert

References: - arXiv: 0906.4819

Beams: p- p+

Beam energies: (980.0, 980.0)GeV

Run details: - $p \bar{p} \to$ jets at 1960 GeV

Dijet angular distributions in different bins of dijet mass from 0.25 TeV to above 1.1 TeV in $p \bar{p}$ collisions at $\sqrt{s} = 1.96$ TeV, based on an integrated luminosity of 0.7 fb$^{-1}$.

Source code:D0_2009_I824127.cc

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

include "Rivet/Analysis.hh"

include "Rivet/Projections/ChargedFinalState.hh"

include "Rivet/Projections/FastJets.hh"

include "Rivet/Projections/FinalState.hh"

namespace Rivet {

/// D0 dijet angular distributions class D0_2009_I824127 : public Analysis { public:

RIVET_DEFAULT_ANALYSIS_CTOR(D0_2009_I824127);


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

// Book histograms
void init() {
  FinalState fs;
  FastJets conefinder(fs, JetAlg::D0ILCONE, 0.7);
  declare(conefinder, "ConeFinder");

  book(_h_chi_dijet, {250., 300., 400., 500., 600., 700., 800., 900., 1000., 1100., 1960.});
  for (auto& b : _h_chi_dijet->bins()) {
    book(b, b.index(), 1, 1);
  }
}


/// Do the analysis
void analyze(const Event& e) {

  const Jets& jets = apply<JetFinder>(e, "ConeFinder").jetsByPt();
  if (jets.size() < 2) vetoEvent;

  FourMomentum j0(jets[0].momentum());
  FourMomentum j1(jets[1].momentum());
  double y0 = j0.rapidity();
  double y1 = j1.rapidity();

  if (fabs(y0 + y1) > 2) vetoEvent;

  double mjj = FourMomentum(j0 + j1).mass();
  double chi = exp(fabs(y0 - y1));
  if (chi < 16.) _h_chi_dijet->fill(mjj, chi);
}


/// Finalize
void finalize() {
  normalize(_h_chi_dijet);
}

/// @}

private:

/// @name Histograms
/// @{
Histo1DGroupPtr _h_chi_dijet;
/// @}

};

RIVET_DECLARE_ALIASED_PLUGIN(D0_2009_I824127, D0_2009_S8320160);

} ```

Aliases: - D0_2009_S8320160