Rivet analyses referenceJADE_1983_I190818Hadronic charged multiplicity measurement between 12 and 35 GeVExperiment: JADE (PETRA) Inspire ID: 190818 Status: VALIDATED Authors:
Beam energies: (6.0, 6.0); (15.0, 15.0); (17.5, 17.5) GeV Run details:
The charged particle multiplicity distribution of hadronic $e^+e^-$ events as measured between 12 and 35 GeV using the JADE detector at PETRA. Beam energy must be specified as analysis option "ENERGY" when rivet-merging samples. Source code: JADE_1983_I190818.cc 1// -*- C++ -*-
2#include "Rivet/Analysis.hh"
3#include "Rivet/Projections/ChargedFinalState.hh"
4
5namespace Rivet {
6
7
8 /// @brief Average multiplcity at a range of energies
9 class JADE_1983_I190818 : public Analysis {
10 public:
11
12 /// Constructor
13 RIVET_DEFAULT_ANALYSIS_CTOR(JADE_1983_I190818);
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 if( !(isCompatibleWithSqrtS(12.0*GeV) ||
24 isCompatibleWithSqrtS(30.0*GeV) ||
25 isCompatibleWithSqrtS(35.0*GeV) )) {
26 MSG_WARNING("CoM energy of events sqrt(s) = " << sqrtS()/GeV
27 << " doesn't match any available analysis energy .");
28 }
29 book(_mult, 1, 1, 1);
30 }
31
32
33 /// Perform the per-event analysis
34 void analyze(const Event& event) {
35 const FinalState& cfs = apply<FinalState>(event, "CFS");
36 MSG_DEBUG("Total charged multiplicity = " << cfs.size());
37 _mult->fill(int(sqrtS()),cfs.size());
38 }
39
40
41 /// Normalise histograms etc., after the run
42 void finalize() {
43 }
44 /// @}
45
46 private:
47
48 // Histogram
49 BinnedProfilePtr<int> _mult;
50
51 };
52
53 RIVET_DECLARE_PLUGIN(JADE_1983_I190818);
54
55}
|
