Rivet analyses

Differential meson cross sections in e+e− at $\sqrt{s}=10.58$GeV

Experiment: BELLE (KEKB)

Inspire ID: 2849895

Status: VALIDATED NOHEPDATA

Authors: - Peter Richardson

References: - Phys.Rev.D 111 (2025) 5, 052003 - arXiv: 2411.12216

Beams: e+ e-

Beam energies: (5.3, 5.3)GeV

Run details: - e+ e- to hadrons

Differential cross sections vs scaled momentum fraction for D+, D0, Ds+, D*+, D*0, Ds*+, ρ+, ρ0, ω, K*+, K*0, ϕ, η, KS0, f0 mesons at $\sqrt{s}=10.58$GeV by the BELLE collaboration.

Source code:BELLE_2025_I2849895.cc

// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/UnstableParticles.hh"

namespace Rivet {


  /// @brief meson spectra at 10.58 GeV
  class BELLE_2025_I2849895 : public Analysis {
  public:

    /// Constructor
    RIVET_DEFAULT_ANALYSIS_CTOR(BELLE_2025_I2849895);


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

    /// Book histograms and initialise projections before the run
    void init() {
      declare(UnstableParticles(), "UFS");
      for (unsigned int ix = 0; ix < 6; ++ix) book(_h[ix], 1, 1, 1 + ix);
      for (unsigned int ix = 0; ix < 9; ++ix) book(_h[6 + ix], 2, 1, 1 + ix);
    }


    /// Perform the per-event analysis
    void analyze(const Event& event) {
      // apply projection
      const UnstableParticles& ufs = apply<UnstableParticles>(event, "UFS");
      // fill single particle histos
      for (const Particle& p : ufs.particles()) {
        // remove particles from radiation/mixing
        bool self = false;
        for (const Particle& child : p.children()) self |= p.abspid() == child.pid();
        if (self) continue;
        const double xp = 2. * p.momentum().p3().mod() / sqrtS();
        switch (p.abspid()) {
          case 411: _h[0]->fill(xp); break;
          case 421: _h[1]->fill(xp); break;
          case 431: _h[2]->fill(xp); break;
          case 413: _h[3]->fill(xp); break;
          case 423: _h[4]->fill(xp); break;
          case 433: _h[5]->fill(xp); break;
          case 213:
            if (p.pid() > 0) _h[6]->fill(xp);
            break;
          case 113: _h[7]->fill(xp); break;
          case 223: _h[8]->fill(xp); break;
          case 323: _h[9]->fill(xp); break;
          case 313: _h[10]->fill(xp); break;
          case 333: _h[11]->fill(xp); break;
          case 221: _h[12]->fill(xp); break;
          case PID::K0S: _h[13]->fill(xp); break;
          case 9010221: _h[14]->fill(xp); break;
          default: break;
        }
      }
    }


    /// Normalise histograms etc., after the run
    void finalize() {
      scale(_h, crossSection() / picobarn / sumOfWeights());
    }

    /// @}


    /// @name Histograms
    /// @{
    Histo1DPtr _h[15];
    /// @}
  };


  RIVET_DECLARE_PLUGIN(BELLE_2025_I2849895);

}