rivet is hosted by Hepforge, IPPP Durham

Rivet analyses reference

CDF_1990_I283352

CDF pseudorapidity distributions at 630 and 1800 GeV
Experiment: CDF (Tevatron Run 0)
Inspire ID: 283352
Status: VALIDATED
Authors:
  • Andy Buckley
References: Beams: p- p+
Beam energies: (315.0, 315.0); (900.0, 900.0) GeV
Run details:
  • QCD min bias events at $\sqrt{s} = 630$ and 1800 GeV. Particles with $c \tau > 10$ mm should be set stable. Beam energy must be specified as analysis option "ENERGY" when rivet-merge'ing samples.

Pseudorapidity distributions based on the CDF 630 and 1800 GeV runs from 1987. All data is detector corrected. The data confirms the UA5 measurement of a $\mathrm{d}{N}/\mathrm{d}{\eta}$ rise with energy faster than $\ln{\sqrt{s}}$, and as such this analysis is important for constraining the energy evolution of minimum bias and underlying event characteristics in MC simulations. Beam energy must be specified (in GeV) as analysis option "ENERGY" when rivet-merging samples.

Source code: CDF_1990_I283352.cc
 1// -*- C++ -*-
 2#include "Rivet/Analysis.hh"
 3#include "Rivet/Projections/ChargedFinalState.hh"
 4#include "Rivet/Projections/TriggerCDFRun0Run1.hh"
 5
 6namespace Rivet {
 7
 8
 9  /// @brief CDF pseudorapidity analysis at 630 and 1800 GeV
10  ///
11  /// @author Andy Buckley
12  class CDF_1990_I283352 : public Analysis {
13  public:
14
15    RIVET_DEFAULT_ANALYSIS_CTOR(CDF_1990_I283352);
16
17
18    /// @name Analysis methods
19    /// @{
20
21    void init() {
22      // Setup projections
23      declare(TriggerCDFRun0Run1(), "Trigger");
24      declare(ChargedFinalState((Cuts::etaIn(-3.5, 3.5))), "CFS");
25
26      // Book histo
27      if (isCompatibleWithSqrtS(1800*GeV)) {
28        book(_hist_eta ,3, 1, 1);
29      } else if (isCompatibleWithSqrtS(630*GeV)) {
30        book(_hist_eta ,4, 1, 1);
31      }
32      book(_sumWTrig, "sumWTrig");
33    }
34
35
36    /// Do the analysis
37    void analyze(const Event& event) {
38      // Trigger
39      const bool trigger = apply<TriggerCDFRun0Run1>(event, "Trigger").minBiasDecision();
40      if (!trigger) vetoEvent;
41      _sumWTrig->fill();
42
43      // Loop over final state charged particles to fill eta histos
44      const FinalState& fs = apply<FinalState>(event, "CFS");
45      for (const Particle& p : fs.particles()) {
46        const double eta = p.eta();
47        _hist_eta->fill(fabs(eta));
48      }
49    }
50
51
52    /// Finalize
53    void finalize() {
54      // Divide through by num events to get d<N>/d(eta) in bins
55      // Factor of 1/2 for |eta| -> eta
56      scale(_hist_eta, 0.5/ *_sumWTrig);
57    }
58
59    /// @}
60
61
62  private:
63
64    /// Counter
65    CounterPtr _sumWTrig;
66
67    /// Histogram
68    Histo1DPtr _hist_eta;
69
70  };
71
72
73
74  RIVET_DECLARE_ALIASED_PLUGIN(CDF_1990_I283352, CDF_1990_S2089246);
75
76}