Rivet is hosted by Hepforge, IPPP Durham

Rivet analyses reference

CDF_2009_NOTE_9936

CDF Run 2 min bias charged multiplicity analysis
Experiment: CDF (Tevatron Run 2)
Spires ID: None
Status: OBSOLETE
Authors:
  • Holger Schulz
References:
  • CDF public note 9936
  • http://www-cdf.fnal.gov/physics/new/qcd/minbias_mult09/multpage.html
Beams: p- p+
Beam energies: (980.0, 980.0) GeV
Run details:
  • $p\bar{p}$ QCD interactions at 1960 GeV. Particles with $c \tau > 10$ mm should be set stable.

Niccolo Moggi's min bias analysis. Minimum bias events are used to measure the charged multiplicity distribution. The multiplicity distribution was not published in S8233977 but the numbers and a public note are available from the CDF website given above. Note: the systematic and statistical errors in Rivet were added in quadrature.

Source code: CDF_2009_NOTE_9936.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
// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/ChargedFinalState.hh"
#include "Rivet/Projections/TriggerCDFRun2.hh"

namespace Rivet {


  class CDF_2009_NOTE_9936 : public Analysis {
  public:

    /// @name Constructors etc.
    //@{

    /// Constructor
    CDF_2009_NOTE_9936()
      : Analysis("CDF_2009_NOTE_9936")
    {    }

    //@}


  public:

    /// @name Analysis methods
    //@{

    /// Book histograms and initialise projections before the run
    void init() {

      declare(TriggerCDFRun2(), "Trigger");

      declare(ChargedFinalState(-1.0, 1.0, 0.4*GeV), "CFS");

      _hist_nch = bookHisto1D(1, 1, 1);
    }


    /// Perform the per-event analysis
    void analyze(const Event& event) {
      // MinBias Trigger
      const bool trigger = apply<TriggerCDFRun2>(event, "Trigger").minBiasDecision();
      if (!trigger) vetoEvent;
      //_sumWPassed += event.weight();
      const double weight = event.weight();

      // Get events charged multiplicity and fill histogram
      const ChargedFinalState& cfs = apply<ChargedFinalState>(event, "CFS");
      _hist_nch->fill(cfs.size(), weight);

    }


    /// Normalise histograms etc., after the run
    void finalize() {
      normalize(_hist_nch);
    }

    //@}

  private:

    Histo1DPtr _hist_nch;

  };



  // The hook for the plugin system
  DECLARE_RIVET_PLUGIN(CDF_2009_NOTE_9936);

}