Rivet analyses referenceCDF_1990_I283352CDF pseudorapidity distributions at 630 and 1800 GeVExperiment: CDF (Tevatron Run 0) Inspire ID: 283352 Status: VALIDATED Authors:
Beam energies: (315.0, 315.0); (900.0, 900.0) GeV Run details:
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}
|
