Rivet analyses
Measurement of Rb 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 Rb 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
// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/UnstableParticles.hh"
namespace Rivet {
/// @brief measurement of R_b between 10.541 and 11.206 GeV
class BABAR_2009_I797507 : public Analysis {
public:
/// Constructor
RIVET_DEFAULT_ANALYSIS_CTOR(BABAR_2009_I797507);
/// @name Analysis methods
/// @{
/// Book histograms and initialise projections before the run
void init() {
// Initialise and register projections
declare(FinalState(), "FS");
declare(UnstableParticles(), "UFS");
// Book histograms
book(_c_hadrons, "/TMP/sigma_hadrons", refData<YODA::BinnedEstimate<string>>(1,1,1));
book(_c_muons, "/TMP/sigma_muons" , refData<YODA::BinnedEstimate<string>>(1,1,1));
for (const string& en : _c_hadrons.binning().edges<0>()) {
double eval = std::stod(en)*GeV;
if (isCompatibleWithSqrtS(eval)) {
_sqs = en;
break;
}
}
raiseBeamErrorIf(_sqs.empty());
}
/// Perform the per-event analysis
void analyze(const Event& event) {
const FinalState& fs = apply<FinalState>(event, "FS");
map<long,int> nCount;
int ntotal(0);
for (const Particle& p : fs.particles()) {
nCount[p.pid()] += 1;
++ntotal;
}
bool isBottom(false);
const UnstableParticles& ufs = apply<UnstableParticles>(event, "UFS");
for (const Particle& p : ufs.particles()) {
if (PID::isBottomHadron(p.pid())) {
isBottom = true;
break;
}
}
// mu+mu- + photons
if (nCount[-13]==1 and
nCount[13]==1 and
ntotal==2+nCount[22]) {
_c_muons->fill(_sqs);
}
else if (isBottom) { // everything else
_c_hadrons->fill(_sqs);
}
}
/// Normalise histograms etc., after the run
void finalize() {
BinnedEstimatePtr<string> mult;
book(mult, 1, 1, 1);
divide(_c_hadrons, _c_muons,mult);
}
/// @}
/// @name Histograms
/// @{
BinnedHistoPtr<string> _c_hadrons, _c_muons;
string _sqs = "";
/// @}
};
RIVET_DECLARE_PLUGIN(BABAR_2009_I797507);
}