Rivet analyses referenceBABAR_2006_I724557Measurement of $\Omega_c^{*0}$ productionExperiment: BABAR (PEP-II) Inspire ID: 724557 Status: VALIDATED Authors:
Beam energies: (5.3, 5.3) GeV Run details:
Measurement of the ratio \[R=\sigma(e^+e^-\to\Omega^{*0}_cX, x_p(\Omega_c^{*0}) > 0.5)/\sigma(e^+e^-\to\Omega^0_cX, x_p(\Omega_c^0) > 0.5),\] by BABAR. This is currently the only measurement of the rate of $\Omega_c^{*0}$ production. , Source code: BABAR_2006_I724557.cc 1// -*- C++ -*-
2#include "Rivet/Analysis.hh"
3#include "Rivet/Projections/UnstableParticles.hh"
4
5namespace Rivet {
6
7
8 /// @brief Omega_c*0 production
9 class BABAR_2006_I724557 : public Analysis {
10 public:
11
12 /// Constructor
13 RIVET_DEFAULT_ANALYSIS_CTOR(BABAR_2006_I724557);
14
15
16 /// @name Analysis methods
17 /// @{
18
19 /// Book histograms and initialise projections before the run
20 void init() {
21 // projections
22 declare(UnstableParticles(),"UFS");
23 // histos
24 book(_c_Omega ,"TMP/Omega" ,refData<YODA::BinnedEstimate<string>>(1,1,1));
25 book(_c_Omega_Star,"TMP/OmegaStar",refData<YODA::BinnedEstimate<string>>(1,1,1));
26 }
27
28
29 /// Perform the per-event analysis
30 void analyze(const Event& event) {
31 for(const Particle & omega : apply<UnstableParticles>(event,"UFS").particles(Cuts::abspid==4332 or
32 Cuts::abspid==4334)) {
33 double xp = omega.momentum().p3().mod()/sqrt(0.25*sqr(sqrtS())-sqr(omega.mass()));
34 if(xp<0.5) continue;
35 if(omega.abspid()==4332)
36 _c_Omega->fill("10.58"s);
37 else
38 _c_Omega_Star->fill("10.58"s);
39 }
40 }
41
42
43 /// Normalise histograms etc., after the run
44 void finalize() {
45 BinnedEstimatePtr<string> ratio;
46 book(ratio,1,1,1);
47 divide(_c_Omega_Star,_c_Omega,ratio);
48 }
49
50 /// @}
51
52
53 /// @name Histograms
54 /// @{
55 BinnedHistoPtr<string> _c_Omega,_c_Omega_Star;
56 /// @}
57
58
59 };
60
61
62 RIVET_DECLARE_PLUGIN(BABAR_2006_I724557);
63
64}
|