Rivet analyses referenceJADE_1979_I142874Charged hadron multiplicityExperiment: JADE (PETRA) Inspire ID: 142874 Status: VALIDATED Authors:
Beam energies: (11.0, 11.0); (13.8, 13.8); (15.0, 15.0); (15.8, 15.8) GeV Run details:
Charged hadron multiplicity at a range of energies, use with care as its not clear if $K_S^0$ decays are included of not, within the systematic error of 1.5 Source code: JADE_1979_I142874.cc 1// -*- C++ -*-
2#include "Rivet/Analysis.hh"
3#include "Rivet/Projections/ChargedFinalState.hh"
4
5namespace Rivet {
6
7
8 /// @brief Add a short analysis description here
9 class JADE_1979_I142874 : public Analysis {
10 public:
11
12 /// Constructor
13 RIVET_DEFAULT_ANALYSIS_CTOR(JADE_1979_I142874);
14
15
16 /// @name Analysis methods
17 /// @{
18
19 /// Book histograms and initialise projections before the run
20 void init() {
21
22 // Initialise and register projections
23 declare(ChargedFinalState(), "FS");
24 if ( isCompatibleWithSqrtS(22.0*GeV) ) _sqs = "22.0";
25 else if( isCompatibleWithSqrtS(27.7*GeV) ) _sqs = "27.7";
26 else if( isCompatibleWithSqrtS(30.0*GeV) ) _sqs = "30.0";
27 else if( isCompatibleWithSqrtS(31.6*GeV) ) _sqs = "31.6";
28 else
29 MSG_WARNING("CoM energy of events sqrt(s) = " << sqrtS()/GeV
30 << " doesn't match any available analysis energy .");
31 // Book histograms
32 book(_mult, 2, 1, 1);
33
34 }
35
36
37 /// Perform the per-event analysis
38 void analyze(const Event& event) {
39 const ChargedFinalState& fs = apply<ChargedFinalState>(event, "FS");
40 _mult->fill(_sqs,fs.particles().size());
41 }
42
43
44 /// Normalise histograms etc., after the run
45 void finalize() {
46 }
47 /// @}
48
49
50 /// @name Histograms
51 /// @{
52 BinnedProfilePtr<string> _mult;
53 string _sqs;
54 /// @}
55
56
57 };
58
59
60 RIVET_DECLARE_PLUGIN(JADE_1979_I142874);
61
62
63}
|