rivet is hosted by Hepforge, IPPP Durham

## Rivet analyses reference

### BELLE_2004_I668024

Rate and spectra for the production of $\Sigma_c(2800)$
Experiment: BELLE (KEKB)
Inspire ID: 668024
Status: VALIDATED
Authors:
• Peter Richardson
References:
• Phys.Rev.Lett. 94 (2005) 122002
Beams: e+ e-
Beam energies: (5.3, 5.3) GeV
Run details:
Measurement of the rate and scaled momentum spectra for the production of $\Sigma_c(2800)^0$, $\Sigma_c(2800)^+$ and $\Sigma_c(2800)^{++}$ baryons by BELLE.
  1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 `// -*- C++ -*- #include "Rivet/Analysis.hh" #include "Rivet/Projections/UnstableParticles.hh" namespace Rivet { /// @brief Sigma_c(2800) rate and spectra class BELLE_2004_I668024 : public Analysis { public: /// Constructor DEFAULT_RIVET_ANALYSIS_CTOR(BELLE_2004_I668024); /// @name Analysis methods ///@{ /// Book histograms and initialise projections before the run void init() { // Initialise and register projections declare(UnstableParticles(), "UFS"); // Book histograms for(unsigned int ix=0;ix<3;++ix) { book(_r[ix], 1, 1, 1+ix); book(_h[ix], 2, 1, 1+ix); } } /// Perform the per-event analysis void analyze(const Event& event) { const UnstableParticles& ufs = apply(event, "UFS"); for (const Particle& p : ufs.particles()) { const Vector3 mom3 = p.p3(); double pp = mom3.mod(); double xp = pp/sqrt(0.25*sqr(sqrtS())-sqr(p.mass())); int id = abs(p.pid()); if(id==_sid+110 && isDecay(p,{4122,-211})) { _h[0]->fill(xp); _r[0]->fill(0.5); } else if(id==_sid+210&& isDecay(p,{4122,111})) { _h[1]->fill(xp); _r[1]->fill(0.5); } else if(id==_sid+220&& isDecay(p,{4122,211})) { _h[2]->fill(xp); _r[2]->fill(0.5); } } } // Check for explicit decay into pdgids bool isDecay(const Particle& mother, vector ids) { if(mother.pid()<0) { for(unsigned int ix=0;ix