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