Rivet analyses

title: BABAR_2006_I724557

Measurement of $\Omega_c^{*0}$ production

Experiment: BABAR (PEP-II)

Inspire ID: 724557

Status: VALIDATED

Authors: - Peter Richardson

References: - Phys.Rev.Lett. 97 (2006) 232001

Beams: e+ e-

Beam energies: (5.3, 5.3)GeV

Run details: - e+e- to hadrons

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

```c++ // -- C++ --

include "Rivet/Analysis.hh"

include "Rivet/Projections/UnstableParticles.hh"

namespace Rivet {

/// @brief Omega_c*0 production class BABAR_2006_I724557 : public Analysis { public:

/// Constructor
RIVET_DEFAULT_ANALYSIS_CTOR(BABAR_2006_I724557);


/// @name Analysis methods
/// @{

/// Book histograms and initialise projections before the run
void init() {
  // projections
  declare(UnstableParticles(), "UFS");
  // histos
  book(_c_Omega, "TMP/Omega", refData<YODA::BinnedEstimate<string>>(1, 1, 1));
  book(_c_Omega_Star, "TMP/OmegaStar", refData<YODA::BinnedEstimate<string>>(1, 1, 1));
}


/// Perform the per-event analysis
void analyze(const Event& event) {
  for (const Particle& omega :
       apply<UnstableParticles>(event, "UFS").particles(Cuts::abspid == 4332 or Cuts::abspid == 4334)) {
    double xp = omega.momentum().p3().mod() / sqrt(0.25 * sqr(sqrtS()) - sqr(omega.mass()));
    if (xp < 0.5) continue;
    if (omega.abspid() == 4332)
      _c_Omega->fill("10.58"s);
    else
      _c_Omega_Star->fill("10.58"s);
  }
}


/// Normalise histograms etc., after the run
void finalize() {
  BinnedEstimatePtr<string> ratio;
  book(ratio, 1, 1, 1);
  divide(_c_Omega_Star, _c_Omega, ratio);
}

/// @}


/// @name Histograms
/// @{
BinnedHistoPtr<string> _c_Omega, _c_Omega_Star;
/// @}

};

RIVET_DECLARE_PLUGIN(BABAR_2006_I724557);

} ```