Rivet analyses reference

BABAR_2009_I797507

Measurement of $R_b$ for energies between 10.541 and 11.206 GeV
Experiment: BABAR (PEP-II)
Inspire ID: 797507
Status: VALIDATED
Authors:

Peter Richardson

References:

PRL 102,012001

Beams: e- e+
Beam energies: (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.3, 5.3); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.4, 5.4); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.5, 5.5); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6); (5.6, 5.6) GeV
Run details:

e+ e- to hadrons and e+ e- to mu+ mu- (for normalization)

Measurement of $R_b$ in $e^+e^-$ collisions by BaBar for energies between 10.541 and 11.206 GeV. The individual bottom 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: BABAR_2009_I797507.cc

 1// -*- C++ -*-
 2#include "Rivet/Analysis.hh"
 3#include "Rivet/Projections/FinalState.hh"
 4#include "Rivet/Projections/UnstableParticles.hh"
 5
 6namespace Rivet {
 7
 8
 9  /// @brief measurement of R_b between 10.541 and 11.206 GeV
10  class BABAR_2009_I797507 : public Analysis {
11  public:
12
13    /// Constructor
14    RIVET_DEFAULT_ANALYSIS_CTOR(BABAR_2009_I797507);
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      declare(UnstableParticles(), "UFS");
26
27      // Book histograms
28      book(_c_hadrons, "/TMP/sigma_hadrons", refData<YODA::BinnedEstimate<string>>(1,1,1));
29      book(_c_muons,   "/TMP/sigma_muons"  , refData<YODA::BinnedEstimate<string>>(1,1,1));
30      for (const string& en : _c_hadrons.binning().edges<0>()) {
31        double end = std::stod(en)*GeV;
32        if(isCompatibleWithSqrtS(end)) {
33          _ecms = en;
34          break;
35        }
36      }
37      if(_ecms.empty()) MSG_ERROR("Beam energy incompatible with analysis.");
38    }
39
40
41    /// Perform the per-event analysis
42    void analyze(const Event& event) {
43      const FinalState& fs = apply<FinalState>(event, "FS");
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      bool isBottom(false);
51      const UnstableParticles& ufs = apply<UnstableParticles>(event, "UFS");
52      for (const Particle& p : ufs.particles()) {
53        if (PID::isBottomHadron(p.pid())) {
54	  isBottom = true;
55	  break;
56	}
57      }
58
59      // mu+mu- + photons
60      if(nCount[-13]==1 and nCount[13]==1 &&
61	 ntotal==2+nCount[22])
62	_c_muons->fill(_ecms);
63      // everything else
64      else if(isBottom) {
65	_c_hadrons->fill(_ecms);
66      }
67    }
68
69
70    /// Normalise histograms etc., after the run
71    void finalize() {
72      BinnedEstimatePtr<string> mult;
73      book(mult, 1, 1, 1);
74      divide(_c_hadrons,_c_muons,mult);
75    }
76
77    /// @}
78
79
80    /// @name Histograms
81    /// @{
82    BinnedHistoPtr<string> _c_hadrons, _c_muons;
83    string _ecms;
84    /// @}
85
86
87  };
88
89
90  RIVET_DECLARE_PLUGIN(BABAR_2009_I797507);
91
92
93}