Rivet analyses referenceLENA_1982_I179431Measurement of $R$ between 7.4 and 9.4 GeVExperiment: LENA (DORIS) Inspire ID: 179431 Status: VALIDATED Authors:
Beam energies: ANY Run details:
Measurement of $R$ in $e^+e^-$ collisions by LENA for energies between .4 and 9.4 GeV. The individual hadronic and muonic cross sections are also outputted to the yoda file so that ratio $R$ can be recalculated if runs are combined. Source code: LENA_1982_I179431.cc 1// -*- C++ -*-
2#include "Rivet/Analysis.hh"
3#include "Rivet/Projections/FinalState.hh"
4
5namespace Rivet {
6
7
8 /// @brief Add a short analysis description here
9 class LENA_1982_I179431 : public Analysis {
10 public:
11
12 /// Constructor
13 RIVET_DEFAULT_ANALYSIS_CTOR(LENA_1982_I179431);
14
15
16 /// @name Analysis methods
17 /// @{
18
19 /// Book histograms and initialise projections before the run
20 void init() {
21 // Initialise and register projections
22 declare(FinalState(), "FS");
23
24 // Book histograms
25 for (size_t i=0; i<2; ++i) {
26 const auto& ref = refData(i+1, 1, 1);
27 book(_h_hadrons[i], "d0"+to_string(i+1)+"_sigma_hadrons", ref);
28 book(_h_muons[i], "d0"+to_string(i+1)+"_sigma_muons", ref);
29 book(_ratios[i], i+1, 1, 1);
30 }
31 book(_c_hadrons, "d03_sigma_hadrons");
32 book(_c_muons, "d03_sigma_muons");
33 }
34
35
36 /// Perform the per-event analysis
37 void analyze(const Event& event) {
38 const FinalState& fs = apply<FinalState>(event, "FS");
39
40 map<long,int> nCount;
41 int ntotal(0);
42 for (const Particle& p : fs.particles()) {
43 nCount[p.pid()] += 1;
44 ++ntotal;
45 }
46 // mu+mu- + photons
47 if (nCount[-13]==1 and nCount[13]==1 && ntotal==2+nCount[22]) {
48 _h_muons[0]->fill(sqrtS()/GeV);
49 _h_muons[1]->fill(sqrtS()/GeV);
50 _c_muons->fill();
51 }
52 // everything else
53 else {
54 _h_hadrons[0]->fill(sqrtS()/GeV);
55 _h_hadrons[1]->fill(sqrtS()/GeV);
56 _c_hadrons->fill();
57 }
58 }
59
60
61 /// Normalise histograms etc., after the run
62 void finalize() {
63 const double fact = crossSection() / sumOfWeights() / picobarn;
64 scale(_h_hadrons, fact);
65 scale(_h_muons, fact);
66 scale({_c_hadrons, _c_muons}, fact);
67 divide(_h_hadrons[0], _h_muons[0], _ratios[0]);
68 divide(_h_hadrons[1], _h_muons[1], _ratios[1]);
69 const Estimate0D R = *_c_hadrons / *_c_muons;
70 BinnedEstimatePtr<string> mult;
71 book(mult, 3, 1, 1);
72 for (auto& b : mult->bins()) {
73 const double Ecm = std::stod(b.xEdge());
74 if (isCompatibleWithSqrtS(Ecm/GeV)) {
75 b.set(R.val(), R.errPos());
76 }
77 }
78 }
79
80 /// @}
81
82
83 /// @name Histograms
84 /// @{
85 CounterPtr _c_hadrons, _c_muons;
86 Histo1DPtr _h_hadrons[2], _h_muons[2];
87 Estimate1DPtr _ratios[2];
88 /// @}
89
90
91 };
92
93
94 RIVET_DECLARE_PLUGIN(LENA_1982_I179431);
95
96
97}
|