1// -*- C++ -*-
 2#include "Rivet/Analysis.hh"
 3#include "Rivet/Projections/ChargedFinalState.hh"
 4
 5namespace Rivet {
 6
 7
 8  /// @brief Charged hadron multiplicity at 29 GeV from HRS experiment
 9  class HRS_1986_I18502 : public Analysis {
10  public:
11
12    /// Constructor
13    RIVET_DEFAULT_ANALYSIS_CTOR(HRS_1986_I18502);
14
15
16    /// @name Analysis methods
17    /// @{
18
19    /// Book histograms and initialise projections before the run
20    void init() {
21      const ChargedFinalState cfs;
22      declare(cfs, "CFS");
23
24      book(_histChTot, 1, 1, 1);
25      book(_histAver,  3, 1, 1);
26    }
27
28
29    /// Perform the per-event analysis
30    void analyze(const Event& event) {
31      const FinalState& cfs = apply<FinalState>(event, "CFS");
32      MSG_DEBUG("Total charged multiplicity = " << cfs.size());
33      _histChTot->fill(cfs.size());
34      _histAver->fill(29, cfs.size());
35    }
36
37
38    /// Normalise histograms etc., after the run
39    void finalize() {
40      scale(_histChTot, 100.0/sumOfWeights()); // and %age (100)
41    }
42
43    /// @}
44
45
46  private:
47
48    /// @name Histograms
49    /// @{
50    BinnedHistoPtr<int> _histChTot;
51    BinnedProfilePtr<int> _histAver;
52    /// @}
53  };
54
55
56  RIVET_DECLARE_PLUGIN(HRS_1986_I18502);
57
58
59}