rivet is hosted by Hepforge, IPPP Durham

Rivet analyses reference

BESIII_2024_I2691955

$\psi(2S)\to K^-\Lambda^0\bar\Xi^++\mathrm{c.c.}$
Experiment: BESIII (BEPC)
Inspire ID: 2691955
Status: VALIDATED NOHEPDATA
Authors:
  • Peter Richardson
References:
  • Phys.Rev.D 109 (2024) 7, 072008
  • arXiv: 2308.15206
Beams: * *
Beam energies: ANY
Run details:
  • Any process producing psi(2S) originally e+e-

Measurement of the mass distributions in the decay $\psi(2S)\to K^-\Lambda^0\bar\Xi^++\mathrm{c.c.}$. The data were read from the plots in the paper. It is also not clear that any resolution effects have been unfolded.

Source code: BESIII_2024_I2691955.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 psi(2S) -> K- Lambda9 Xibar+
10  class BESIII_2024_I2691955 : public Analysis {
11  public:
12
13    /// Constructor
14    RIVET_DEFAULT_ANALYSIS_CTOR(BESIII_2024_I2691955);
15
16
17    /// @name Analysis methods
18    /// @{
19
20    /// Book histograms and initialise projections before the run
21    void init() {
22      UnstableParticles ufs = UnstableParticles(Cuts::pid==100443);
23      declare(ufs, "UFS");
24      DecayedParticles psi(ufs);
25      psi.addStable(PID::PI0);
26      psi.addStable( 3122);
27      psi.addStable(-3122);
28      psi.addStable( 3312);
29      psi.addStable(-3312);
30      declare(psi, "psi");
31      for (unsigned int iy=0; iy<3; ++iy)
32        book(_h[iy], 1, 1, 1+iy);
33    }
34
35
36    /// Perform the per-event analysis
37    void analyze(const Event& event) {
38      DecayedParticles psi = apply<DecayedParticles>(event, "psi");
39      // loop over particles
40      for (unsigned int ix=0;ix<psi.decaying().size();++ix) {
41        int isgn=1;
42       	if     (psi.modeMatches(ix,3,mode1)) isgn= 1;
43        else if(psi.modeMatches(ix,3,mode2)) isgn=-1;
44        else continue;
45       	const Particle& Km  = psi.decayProducts()[ix].at(-isgn* 321)[0];
46       	const Particle& lam = psi.decayProducts()[ix].at( isgn*3122)[0];
47       	const Particle& xi  = psi.decayProducts()[ix].at(-isgn*3312)[0];
48        _h[0]->fill((Km .mom()+xi .mom()).mass());
49        _h[1]->fill((lam.mom()+xi .mom()).mass());
50        _h[2]->fill((Km .mom()+lam.mom()).mass());
51      }
52    }
53
54
55    /// Normalise histograms etc., after the run
56    void finalize() {
57      normalize(_h);
58    }
59
60    /// @}
61
62
63    /// @name Histograms
64    /// @{
65    Histo1DPtr _h[3];
66    const map<PdgId,unsigned int> mode1 = { {-321,1}, { 3122,1}, {-3312,1} };
67    const map<PdgId,unsigned int> mode2 = { { 321,1}, {-3122,1}, { 3312,1} };
68    /// @}
69
70
71  };
72
73
74  RIVET_DECLARE_PLUGIN(BESIII_2024_I2691955);
75
76}