Rivet Analyses Reference

CDF_1990_S2089246

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_S2089246.cc
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76

// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/ChargedFinalState.hh"
#include "Rivet/Projections/TriggerCDFRun0Run1.hh"

namespace Rivet {


  /// @brief CDF pseudorapidity analysis at 630 and 1800 GeV
  ///
  /// @author Andy Buckley
  class CDF_1990_S2089246 : public Analysis {
  public:

    RIVET_DEFAULT_ANALYSIS_CTOR(CDF_1990_S2089246);


    /// @name Analysis methods
    //@{

    void init() {
      // Setup projections
      declare(TriggerCDFRun0Run1(), "Trigger");
      declare(ChargedFinalState((Cuts::etaIn(-3.5, 3.5))), "CFS");

      // Book histo
      if (isCompatibleWithSqrtS(1800)) {
        book(_hist_eta ,3, 1, 1);
      } else if (isCompatibleWithSqrtS(630)) {
        book(_hist_eta ,4, 1, 1);
      }
      book(_sumWTrig, "sumWTrig");
    }


    /// Do the analysis
    void analyze(const Event& event) {
      // Trigger
      const bool trigger = apply<TriggerCDFRun0Run1>(event, "Trigger").minBiasDecision();
      if (!trigger) vetoEvent;
      _sumWTrig->fill();

      // Loop over final state charged particles to fill eta histos
      const FinalState& fs = apply<FinalState>(event, "CFS");
      for (const Particle& p : fs.particles()) {
        const double eta = p.eta();
        _hist_eta->fill(fabs(eta));
      }
    }


    /// Finalize
    void finalize() {
      // Divide through by num events to get d<N>/d(eta) in bins
      // Factor of 1/2 for |eta| -> eta
      scale(_hist_eta, 0.5/ *_sumWTrig);
    }

    //@}


  private:

    /// Counter
    CounterPtr _sumWTrig;

    /// Histogram
    Histo1DPtr _hist_eta;

  };



  RIVET_DECLARE_ALIASED_PLUGIN(CDF_1990_S2089246, CDF_1990_I283352);

}

Edit this page on GitLab