Rivet analyses
Differential Branching ratio and asymmetries in B → K(*)ℓ+ℓ−
Experiment: BABAR (PEP-II)
Inspire ID: 1111233
Status: VALIDATED NOHEPDATA SINGLEWEIGHT
Authors: - Peter Richardson
References: - Phys.Rev.D 86 (2012) 032012
Beams: * *
Beam energies: ANY
Run details: - Any proces producing B0 and B+, originally Upsilon(4S) decays
Measurement of the differential Branching ratio and asymmetries in B → K(*)ℓ+ℓ−. As well as the differential branching ratio as a function of q2 the ratio of e+e− and μ+μ− and isospin and CP asymmetries are measured.
Source
code:BABAR_2012_I1111233.cc
// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/DecayedParticles.hh"
#include "Rivet/Projections/UnstableParticles.hh"
namespace Rivet {
/// @brief B -> K(*) l+ l-
class BABAR_2012_I1111233 : public Analysis {
public:
/// Constructor
RIVET_DEFAULT_ANALYSIS_CTOR(BABAR_2012_I1111233);
/// @name Analysis methods
/// @{
/// Book histograms and initialise projections before the run
void init() {
// Initialise and register projections
UnstableParticles ufs = UnstableParticles(Cuts::abspid == 511 or Cuts::abspid == 521);
declare(ufs, "UFS");
DecayedParticles BB(ufs);
BB.addStable(443);
BB.addStable(100443);
BB.addStable(310);
BB.addStable(313);
BB.addStable(323);
BB.addStable(-313);
BB.addStable(-323);
declare(BB, "BB");
// histograms
for (unsigned int ix = 0; ix < 2; ++ix) {
for (unsigned int iy = 0; iy < 2; ++iy) {
book(_h_br[ix][iy], 1, 1 + ix, 1 + iy);
book(_p_CP[ix][iy], 2, 1 + ix, 1 + iy);
for (unsigned int il = 0; il < 2; ++il) {
book(_h_br_l[il][ix][iy], "TMP/h_br_l_" + toString(il) + "_" + toString(ix) + "_" + toString(iy),
refData(3, 1 + ix, 1 + iy));
book(_h_br_I[il][ix][iy], "TMP/h_br_I_" + toString(il) + "_" + toString(ix) + "_" + toString(iy),
refData(4, 1 + ix, 1 + iy));
}
}
}
for (unsigned int ix = 0; ix < 3; ++ix) book(_c[ix], "TMP/c_" + toString(ix + 1));
}
/// Perform the per-event analysis
void analyze(const Event& event) {
// kaon modes
static const map<PdgId, unsigned int>& mode1 = {{321, 1}, {13, 1}, {-13, 1}};
static const map<PdgId, unsigned int>& mode1CC = {{-321, 1}, {13, 1}, {-13, 1}};
static const map<PdgId, unsigned int>& mode2 = {{310, 1}, {13, 1}, {-13, 1}};
static const map<PdgId, unsigned int>& mode3 = {{321, 1}, {11, 1}, {-11, 1}};
static const map<PdgId, unsigned int>& mode3CC = {{-321, 1}, {11, 1}, {-11, 1}};
static const map<PdgId, unsigned int>& mode4 = {{310, 1}, {11, 1}, {-11, 1}};
// K* modes
static const map<PdgId, unsigned int>& mode5 = {{323, 1}, {13, 1}, {-13, 1}};
static const map<PdgId, unsigned int>& mode5CC = {{-323, 1}, {13, 1}, {-13, 1}};
static const map<PdgId, unsigned int>& mode6 = {{313, 1}, {13, 1}, {-13, 1}};
static const map<PdgId, unsigned int>& mode6CC = {{-313, 1}, {13, 1}, {-13, 1}};
static const map<PdgId, unsigned int>& mode7 = {{323, 1}, {11, 1}, {-11, 1}};
static const map<PdgId, unsigned int>& mode7CC = {{-323, 1}, {11, 1}, {-11, 1}};
static const map<PdgId, unsigned int>& mode8 = {{313, 1}, {11, 1}, {-11, 1}};
static const map<PdgId, unsigned int>& mode8CC = {{-313, 1}, {11, 1}, {-11, 1}};
DecayedParticles BB = apply<DecayedParticles>(event, "BB");
// loop over particles
for (unsigned int ix = 0; ix < BB.decaying().size(); ++ix) {
_c[0]->fill();
if (BB.decaying()[ix].abspid() == 521)
_c[1]->fill();
else
_c[2]->fill();
int imode = 0;
if ((BB.decaying()[ix].pid() > 0 && BB.modeMatches(ix, 3, mode1))
|| (BB.decaying()[ix].pid() < 0 && BB.modeMatches(ix, 3, mode1CC)))
imode = 0;
else if (BB.modeMatches(ix, 3, mode2))
imode = 1;
else if ((BB.decaying()[ix].pid() > 0 && BB.modeMatches(ix, 3, mode3))
|| (BB.decaying()[ix].pid() < 0 && BB.modeMatches(ix, 3, mode3CC)))
imode = 2;
else if (BB.modeMatches(ix, 3, mode4))
imode = 3;
else if ((BB.decaying()[ix].pid() > 0 && BB.modeMatches(ix, 3, mode5))
|| (BB.decaying()[ix].pid() < 0 && BB.modeMatches(ix, 3, mode5CC)))
imode = 4;
else if ((BB.decaying()[ix].pid() > 0 && BB.modeMatches(ix, 3, mode6))
|| (BB.decaying()[ix].pid() < 0 && BB.modeMatches(ix, 3, mode6CC)))
imode = 5;
else if ((BB.decaying()[ix].pid() > 0 && BB.modeMatches(ix, 3, mode7))
|| (BB.decaying()[ix].pid() < 0 && BB.modeMatches(ix, 3, mode7CC)))
imode = 6;
else if ((BB.decaying()[ix].pid() > 0 && BB.modeMatches(ix, 3, mode8))
|| (BB.decaying()[ix].pid() < 0 && BB.modeMatches(ix, 3, mode8CC)))
imode = 7;
else
continue;
int il = imode < 2 || imode == 4 || imode == 5 ? 13 : 11;
const Particle& lp = BB.decayProducts()[ix].at(-il)[0];
const Particle& lm = BB.decayProducts()[ix].at(il)[0];
double qq = (lp.momentum() + lm.momentum()).mass2();
double ACP = BB.decaying()[ix].pid() > 0 ? -1. : 1.;
if (imode < 4) {
double wgt = (imode == 1 || imode == 3) ? 2 : 1;
for (unsigned int iy = 0; iy < 2; ++iy) {
_h_br[iy][0]->fill(qq, wgt);
_p_CP[iy][0]->fill(qq, ACP, wgt);
if (il == 13)
_h_br_l[0][iy][0]->fill(qq, wgt);
else
_h_br_l[1][iy][0]->fill(qq, wgt);
if (BB.decaying()[ix].abspid() == 521)
_h_br_I[0][iy][0]->fill(qq, wgt);
else
_h_br_I[1][iy][0]->fill(qq, wgt);
}
}
else {
for (unsigned int iy = 0; iy < 2; ++iy) {
_h_br[iy][1]->fill(qq);
_p_CP[iy][1]->fill(qq, ACP);
if (BB.decaying()[ix].abspid() == 521)
_h_br_I[0][iy][1]->fill(qq);
else
_h_br_I[1][iy][1]->fill(qq);
if (il == 13)
_h_br_l[0][iy][1]->fill(qq);
else
_h_br_l[1][iy][1]->fill(qq);
}
}
}
}
/// Normalise histograms etc., after the run
void finalize() {
// ratio of lifetimes
double rLife = 1. / 1.078;
for (unsigned int ix = 0; ix < 2; ++ix) {
for (unsigned int iy = 0; iy < 2; ++iy) {
scale(_h_br[ix][iy], 1e7 / *_c[0]);
for (unsigned int il = 0; il < 2; ++il) {
scale(_h_br_l[il][ix][iy], 1e7 / *_c[0]);
scale(_h_br_I[il][ix][iy], 1e7 / *_c[il + 1]);
if (il == 0) scale(_h_br_I[il][ix][iy], rLife);
}
// RK plots
Estimate1DPtr RK;
book(RK, 3, 1 + ix, 1 + iy);
divide(_h_br_l[0][ix][iy], _h_br_l[1][ix][iy], RK);
book(RK, 4, 1 + ix, 1 + iy);
asymm(_h_br_I[1][ix][iy], _h_br_I[0][ix][iy], RK);
}
}
}
/// @}
/// @name Histograms
/// @{
Histo1DPtr _h_br[2][2], _h_br_l[2][2][2], _h_br_I[2][2][2];
Profile1DPtr _p_CP[2][2];
CounterPtr _c[3];
/// @}
};
RIVET_DECLARE_PLUGIN(BABAR_2012_I1111233);
}