Rivet analyses referenceBABAR_2008_I767506$J/\psi\omega$ mass in $B\to J/\psi\omega K$Experiment: BABAR (PEP-II) Inspire ID: 767506 Status: VALIDATED NOHEPDATA Authors:
Beam energies: ANY Run details:
Measurement of the $J/\psi\omega$ mass distribution in $B\to J/\psi\omega K$. The background subtracted corrected data were read from figure 3 in the paper. Source code: BABAR_2008_I767506.cc 1// -*- C++ -*-
2#include "Rivet/Analysis.hh"
3#include "Rivet/Projections/UnstableParticles.hh"
4#include "Rivet/Projections/DecayedParticles.hh"
5
6namespace Rivet {
7
8
9 /// @brief B -> J/psi omega K
10 class BABAR_2008_I767506 : public Analysis {
11 public:
12
13 /// Constructor
14 RIVET_DEFAULT_ANALYSIS_CTOR(BABAR_2008_I767506);
15
16
17 /// @name Analysis methods
18 /// @{
19
20 /// Book histograms and initialise projections before the run
21 void init() {
22 // Initialise and register projections
23 UnstableParticles ufs = UnstableParticles(Cuts::abspid==511||
24 Cuts::abspid==521);
25 declare(ufs, "UFS");
26 DecayedParticles BB(ufs);
27 BB.addStable(310);
28 BB.addStable(223);
29 BB.addStable(443);
30 declare(BB, "BB");
31 // histos
32 for(unsigned int ix=0;ix<2;++ix)
33 book(_h[ix],1,1,1+ix);
34 }
35
36
37 /// Perform the per-event analysis
38 void analyze(const Event& event) {
39 static const map<PdgId,unsigned int> & mode1 = { { 443,1},{ 223,1}, { 321,1}};
40 static const map<PdgId,unsigned int> & mode1CC = { { 443,1},{ 223,1}, {-321,1}};
41 static const map<PdgId,unsigned int> & mode2 = { { 443,1},{ 223,1}, { 130,1}};
42 static const map<PdgId,unsigned int> & mode2CC = { { 443,1},{ 223,1}, { 310,1}};
43 DecayedParticles BB = apply<DecayedParticles>(event, "BB");
44 // loop over particles
45 for(unsigned int ix=0;ix<BB.decaying().size();++ix) {
46 unsigned int imode=0;
47 if (BB.modeMatches(ix,3,mode1) ||
48 BB.modeMatches(ix,3,mode1CC))
49 imode=0;
50 else if(BB.modeMatches(ix,3,mode2) ||
51 BB.modeMatches(ix,3,mode2CC))
52 imode=1;
53 else
54 continue;
55 const Particle & omega = BB.decayProducts()[ix].at(223)[0];
56 const Particle & psi = BB.decayProducts()[ix].at(443)[0];
57 _h[imode]->fill((omega.momentum()+psi.momentum()).mass());
58 }
59 }
60
61
62 /// Normalise histograms etc., after the run
63 void finalize() {
64 for(unsigned int ix=0;ix<2;++ix)
65 normalize(_h[ix],1.,false);
66 }
67
68 /// @}
69
70
71 /// @name Histograms
72 /// @{
73 Histo1DPtr _h[2];
74 /// @}
75
76
77 };
78
79
80 RIVET_DECLARE_PLUGIN(BABAR_2008_I767506);
81
82}
|