Rivet analyses referenceBABAR_2004_I656680Cross section for $e^+e^-\to$ $\pi^+\pi^-\pi^0$ between 1.05 and 2.325 GeVExperiment: BABAR (PEP-II) Inspire ID: 656680 Status: VALIDATED Authors:
Beam energies: (0.5, 0.5); (0.5, 0.5); (0.6, 0.6); (0.6, 0.6); (0.6, 0.6); (0.6, 0.6); (0.6, 0.6); (0.6, 0.6); (0.6, 0.6); (0.6, 0.6); (0.7, 0.7); (0.7, 0.7); (0.7, 0.7); (0.7, 0.7); (0.7, 0.7); (0.7, 0.7); (0.7, 0.7); (0.7, 0.7); (0.8, 0.8); (0.8, 0.8); (0.8, 0.8); (0.8, 0.8); (0.8, 0.8); (0.8, 0.8); (0.8, 0.8); (0.8, 0.8); (0.9, 0.9); (0.9, 0.9); (0.9, 0.9); (0.9, 0.9); (0.9, 0.9); (0.9, 0.9); (0.9, 0.9); (0.9, 0.9); (1.0, 1.0); (1.0, 1.0); (1.0, 1.0); (1.0, 1.0); (1.0, 1.0); (1.0, 1.0); (1.0, 1.0); (1.0, 1.0); (1.1, 1.1); (1.1, 1.1); (1.1, 1.1); (1.1, 1.1); (1.1, 1.1); (1.1, 1.1); (1.1, 1.1); (1.1, 1.1); (1.2, 1.2); (1.2, 1.2); (1.2, 1.2); (1.2, 1.2); (1.2, 1.2); (1.2, 1.2); (1.2, 1.2); (1.2, 1.2); (1.3, 1.3); (1.3, 1.3); (1.3, 1.3); (1.3, 1.3); (1.3, 1.3); (1.3, 1.3); (1.3, 1.3); (1.3, 1.3); (1.4, 1.4); (1.4, 1.4); (1.4, 1.4); (1.4, 1.4); (1.4, 1.4); (1.4, 1.4); (1.4, 1.4); (1.4, 1.4); (1.5, 1.5); (1.5, 1.5); (1.5, 1.5); (1.5, 1.5) GeV Run details:
Measurement of the cross section for $e^+e^-\to \pi^+\pi^-\pi^0$ by BaBar using radiative return for energies between 1.05 and 2.325 GeV. Beam energy must be specified as analysis option "ENERGY" when rivet-merging samples. Source code: BABAR_2004_I656680.cc 1// -*- C++ -*-
2#include "Rivet/Analysis.hh"
3#include "Rivet/Projections/FinalState.hh"
4
5
6namespace Rivet {
7
8
9 /// @brief e+e- -> pi+pi-pi0
10 class BABAR_2004_I656680 : public Analysis {
11 public:
12
13 /// Constructor
14 RIVET_DEFAULT_ANALYSIS_CTOR(BABAR_2004_I656680);
15
16
17 /// @name Analysis methods
18 /// @{
19
20 /// Book histograms and initialise projections before the run
21 void init() {
22
23 // Initialise and register projections
24 declare(FinalState(), "FS");
25
26 book(_num3pi, 1, 1, 1);
27 for (const string& en : _num3pi.binning().edges<0>()) {
28 const double end = std::stod(en)*GeV;
29 if (isCompatibleWithSqrtS(end)) {
30 _ecms = en;
31 break;
32 }
33 }
34 if (_ecms.empty())
35 MSG_ERROR("Beam energy incompatible with analysis.");
36
37 }
38
39
40 /// Perform the per-event analysis
41 void analyze(const Event& event) {
42 const FinalState& fs = apply<FinalState>(event, "FS");
43
44 map<long,int> nCount;
45 int ntotal(0);
46 for (const Particle& p : fs.particles()) {
47 nCount[p.pid()] += 1;
48 ++ntotal;
49 }
50 if(ntotal!=3) vetoEvent;
51 if(nCount[-211]==1&&nCount[211]==1&&nCount[111]==1)
52 _num3pi->fill(_ecms);
53
54 }
55
56
57 /// Normalise histograms etc., after the run
58 void finalize() {
59 scale(_num3pi, crossSection()/ sumOfWeights() /nanobarn);
60 }
61
62 /// @}
63
64
65 /// @name Histograms
66 /// @{
67 BinnedHistoPtr<string> _num3pi;
68 string _ecms;
69 /// @}
70
71
72 };
73
74
75 RIVET_DECLARE_PLUGIN(BABAR_2004_I656680);
76
77
78}
|
