Rivet analyses

title: KLOE_2002_I585183

Mass spectrum of $\pi^0\pi^0$ in $\phi\to\pi^0\pi^0\gamma$ decays

Experiment: KLOE (DAPHNE)

Inspire ID: 585183

Status: VALIDATED

Authors: - Peter Richardson

References: - Phys.Lett. B485 (2000) 349-356

Beams: * *

Beam energies: ANY

Run details: - Any process producing phi mesons

Mass spectra for $\pi^0\pi^0$ in $\phi$ decays to $\pi^0\pi^0\gamma$ measured by KLOE.

Source code:KLOE_2002_I585183.cc

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

include "Rivet/Analysis.hh"

include "Rivet/Projections/UnstableParticles.hh"

namespace Rivet {

/// @brief phi -> pi0 pi0 gamma class KLOE_2002_I585183 : public Analysis { public:

/// Constructor
RIVET_DEFAULT_ANALYSIS_CTOR(KLOE_2002_I585183);


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

/// Book histograms and initialise projections before the run
void init() {
  declare(UnstableParticles(), "UFS");
  book(_h_pipi, 1, 1, 1);
  book(_nPhi, "TMP/nPhi");
}

void findDecayProducts(const Particle& mother,
                       unsigned int& nstable,
                       unsigned int& npi,
                       unsigned int& ngamma,
                       FourMomentum& ptot) {
  for (const Particle& p : mother.children()) {
    int id = p.pid();
    if (id == PID::PI0) {
      ++npi;
      ++nstable;
      ptot += p.momentum();
    }
    else if (id == PID::GAMMA) {
      ++ngamma;
      ++nstable;
    }
    else if (id == PID::PIPLUS || id == PID::PIMINUS) {
      ++nstable;
    }
    else if (!p.children().empty()) {
      findDecayProducts(p, nstable, npi, ngamma, ptot);
    }
    else
      ++nstable;
  }
}

/// Perform the per-event analysis
void analyze(const Event& event) {

  // Loop over phis
  for (const Particle& phi : apply<UnstableParticles>(event, "UFS").particles(Cuts::abspid == PID::PHI)) {
    _nPhi->fill();
    unsigned int nstable(0), npi(0), ngamma(0);
    FourMomentum p_tot(0, 0, 0, 0);
    findDecayProducts(phi, nstable, npi, ngamma, p_tot);
    if (nstable != 3) continue;
    if (npi == 2 && ngamma == 1) {
      _h_pipi->fill(p_tot.mass() / MeV);
    }
  }
}


/// Normalise histograms etc., after the run
void finalize() {

  scale(_h_pipi, 1. / *_nPhi);
}

/// @}


/// @name Histograms
/// @{
Histo1DPtr _h_pipi;
CounterPtr _nPhi;
/// @}

};

RIVET_DECLARE_PLUGIN(KLOE_2002_I585183);

} ```