diff --git a/analyses/pluginFrascati/KLOE_2009_I818106.cc b/analyses/pluginFrascati/KLOE_2009_I818106.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginFrascati/KLOE_2009_I818106.cc @@ -0,0 +1,98 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/Projections/UnstableParticles.hh" + +namespace Rivet { + + + /// @brief Add a short analysis description here + class KLOE_2009_I818106 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(KLOE_2009_I818106); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + declare(UnstableParticles(), "UFS"); + _h_etapi = bookHisto1D( 1, 1, 1); + _nPhi=0.; + } + + void findDecayProducts(const Particle & mother, unsigned int & nstable, + unsigned int & neta, unsigned int & npi, + unsigned int & ngamma, FourMomentum & ptot) { + for(const Particle & p : mother.children()) { + int id = p.pdgId(); + if ( id == PID::ETA ) { + ++neta; + ++nstable; + ptot += p.momentum(); + } + else 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, neta, npi, ngamma, ptot); + } + else + ++nstable; + } + } + + /// Perform the per-event analysis + void analyze(const Event& event) { + + // Loop over phis + for(const Particle& phi : apply(event, "UFS").particles(Cuts::abspid==PID::PHI)) { + _nPhi+=event.weight(); + unsigned int nstable(0),neta(0),npi(0),ngamma(0); + FourMomentum p_tot(0,0,0,0); + findDecayProducts(phi, nstable, neta, npi, ngamma, p_tot); + if(nstable!=3) continue; + if(neta==1 && npi==1 && ngamma==1 ) { + _h_etapi->fill(p_tot.mass()/MeV, event.weight()); + } + } + + } + + + /// Normalise histograms etc., after the run + void finalize() { + // normalise to total no of phi mesons + // and mult by 10^7 due normalisation in paper + scale( _h_etapi, 1./_nPhi*1e7); + } + + //@} + + + /// @name Histograms + //@{ + Histo1DPtr _h_etapi; + double _nPhi; + //@} + + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(KLOE_2009_I818106); + + +} diff --git a/analyses/pluginFrascati/KLOE_2009_I818106.info b/analyses/pluginFrascati/KLOE_2009_I818106.info new file mode 100644 --- /dev/null +++ b/analyses/pluginFrascati/KLOE_2009_I818106.info @@ -0,0 +1,32 @@ +Name: KLOE_2009_I818106 +Year: 2009 +Summary: Mass spectrum of $\eta\pi$ in $\phi\to\eta\pi^0\gamma$ decays +Experiment: KLOE +Collider: DAPHNE +InspireID: 818106 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - Phys.Lett. B681 (2009) 5-13 +RunInfo: Any process producing phi mesons +Description: + 'Mass spectra for $\eta\pi$ in $\phi$ decays to $\eta\pi^0\gamma$ + measured by KLOE. Useful for teting he treatment of the a_0(980) meson.' +Keywords: [] +BibKey: Ambrosino:2009py +BibTeX: '@article{Ambrosino:2009py, + author = "Ambrosino, F. and others", + title = "{Study of the a(0)(980) meson via the radiative decay + $\phi \to \eta \pi_0 \gamma$ with the KLOE detector}", + collaboration = "KLOE", + journal = "Phys. Lett.", + volume = "B681", + year = "2009", + pages = "5-13", + doi = "10.1016/j.physletb.2009.09.022", + eprint = "0904.2539", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + SLACcitation = "%%CITATION = ARXIV:0904.2539;%%" +}' diff --git a/analyses/pluginFrascati/KLOE_2009_I818106.plot b/analyses/pluginFrascati/KLOE_2009_I818106.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginFrascati/KLOE_2009_I818106.plot @@ -0,0 +1,6 @@ +BEGIN PLOT /KLOE_2009_I818106/d01-x01-y01 +Title=$\eta\pi^0$ mass in $\phi^0\to\eta\pi^0\gamma$ decays +XLabel=$m_{\eta\pi^0}$ [$\text{MeV}$] +YLabel=$1/\Gamma\text{d}\Gamma/\text{d}m_{\eta\pi^0}\times10^7$ [$\text{MeV}^{-1}$] +LogY=0 +END PLOT diff --git a/analyses/pluginFrascati/KLOE_2009_I818106.yoda b/analyses/pluginFrascati/KLOE_2009_I818106.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginFrascati/KLOE_2009_I818106.yoda @@ -0,0 +1,57 @@ +BEGIN YODA_SCATTER2D_V2 /REF/KLOE_2009_I818106/d01-x01-y01 +IsRef: 1 +Path: /REF/KLOE_2009_I818106/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +691.53000 3.17500 3.17500 6.00000e-02 7.00000e-02 7.00000e-02 +697.88000 3.17500 3.17500 1.80000e-01 1.00000e-01 1.00000e-01 +704.24000 3.17500 3.17500 1.80000e-01 1.20000e-01 1.20000e-01 +710.59000 3.17500 3.17500 3.10000e-01 1.30000e-01 1.30000e-01 +716.94000 3.17500 3.17500 3.00000e-01 8.00000e-02 8.00000e-02 +723.29000 3.17500 3.17500 3.80000e-01 1.10000e-01 1.10000e-01 +729.65000 3.17500 3.17500 5.30000e-01 1.70000e-01 1.70000e-01 +736.00000 3.17500 3.17500 5.10000e-01 1.30000e-01 1.30000e-01 +742.35000 3.17500 3.17500 5.30000e-01 5.00000e-02 5.00000e-02 +748.71000 3.17500 3.17500 6.70000e-01 7.00000e-02 7.00000e-02 +755.06000 3.17500 3.17500 8.10000e-01 7.00000e-02 7.00000e-02 +761.41000 3.17500 3.17500 9.40000e-01 1.00000e-01 1.00000e-01 +767.76000 3.17500 3.17500 9.90000e-01 1.10000e-01 1.10000e-01 +774.12000 3.17500 3.17500 9.90000e-01 8.00000e-02 8.00000e-02 +780.47000 3.17500 3.17500 1.08000e+00 9.00000e-02 9.00000e-02 +786.82000 3.17500 3.17500 1.30000e+00 1.00000e-01 1.00000e-01 +793.18000 3.17500 3.17500 1.27000e+00 1.30000e-01 1.30000e-01 +799.53000 3.17500 3.17500 1.42000e+00 2.80000e-01 2.80000e-01 +805.88000 3.17500 3.17500 1.63000e+00 2.80000e-01 2.80000e-01 +812.24000 3.17500 3.17500 1.71000e+00 1.40000e-01 1.40000e-01 +818.59000 3.17500 3.17500 1.79000e+00 1.60000e-01 1.60000e-01 +824.94000 3.17500 3.17500 1.66000e+00 1.80000e-01 1.80000e-01 +831.29000 3.17500 3.17500 1.82000e+00 1.50000e-01 1.50000e-01 +837.65000 3.17500 3.17500 1.96000e+00 1.20000e-01 1.20000e-01 +844.00000 3.17500 3.17500 2.13000e+00 1.30000e-01 1.30000e-01 +850.35000 3.17500 3.17500 2.25000e+00 1.30000e-01 1.30000e-01 +856.71000 3.17500 3.17500 2.35000e+00 1.40000e-01 1.40000e-01 +863.06000 3.17500 3.17500 2.27000e+00 1.30000e-01 1.30000e-01 +869.41000 3.17500 3.17500 2.35000e+00 1.30000e-01 1.30000e-01 +875.76000 3.17500 3.17500 2.42000e+00 1.60000e-01 1.60000e-01 +882.12000 3.17500 3.17500 2.59000e+00 1.60000e-01 1.60000e-01 +888.47000 3.17500 3.17500 2.80000e+00 1.40000e-01 1.40000e-01 +894.82000 3.17500 3.17500 2.92000e+00 1.90000e-01 1.90000e-01 +901.18000 3.17500 3.17500 3.18000e+00 2.00000e-01 2.00000e-01 +907.53000 3.17500 3.17500 3.37000e+00 1.70000e-01 1.70000e-01 +913.88000 3.17500 3.17500 3.48000e+00 1.70000e-01 1.70000e-01 +920.24000 3.17500 3.17500 3.67000e+00 1.70000e-01 1.70000e-01 +926.59000 3.17500 3.17500 3.94000e+00 1.70000e-01 1.70000e-01 +932.94000 3.17500 3.17500 4.29000e+00 2.50000e-01 2.50000e-01 +939.29000 3.17500 3.17500 4.63000e+00 2.50000e-01 2.50000e-01 +945.65000 3.17500 3.17500 4.89000e+00 2.10000e-01 2.10000e-01 +952.00000 3.17500 3.17500 5.20000e+00 2.20000e-01 2.20000e-01 +958.35000 3.17500 3.17500 5.40000e+00 2.80000e-01 2.80000e-01 +964.71000 3.17500 3.17500 5.44000e+00 3.30000e-01 3.30000e-01 +971.06000 3.17500 3.17500 5.35000e+00 2.20000e-01 2.20000e-01 +977.41000 3.17500 3.17500 4.94000e+00 2.10000e-01 2.10000e-01 +983.76000 3.17500 3.17500 4.02000e+00 1.90000e-01 1.90000e-01 +990.12000 3.17500 3.17500 2.80000e+00 2.70000e-01 2.70000e-01 +996.47000 3.17500 3.17500 1.51000e+00 3.20000e-01 3.20000e-01 +END YODA_SCATTER2D_V2