Page Menu
Home
HEPForge
Search
Configure Global Search
Log In
Files
F7879358
No One
Temporary
Actions
View File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Flag For Later
Size
142 KB
Subscribers
None
View Options
diff --git a/analyses/Makefile.am b/analyses/Makefile.am
--- a/analyses/Makefile.am
+++ b/analyses/Makefile.am
@@ -1,60 +1,61 @@
AUTOMAKE_OPTIONS = -Wno-portability
SHELL := bash -O extglob
if ENABLE_ANALYSES
PLUGIN_DIRS = $(shell ls -d $(srcdir)/plugin*)
PLUGIN_NAMES = $(notdir $(subst plugin,,$(PLUGIN_DIRS)))
PLUGIN_LIBS = $(patsubst %,Rivet%Analyses.so,$(PLUGIN_NAMES))
PLUGIN_DATAFILES = $(shell ls $(abs_srcdir)/plugin*/*.{info,plot,yoda})
CLEANFILES = $(PLUGIN_LIBS)
EXTRA_DIST = $(PLUGIN_DIRS)
%.so:
@+echo && RIVET_BUILDPLUGIN_BEFORE_INSTALL=1 bash $(top_builddir)/bin/rivet-buildplugin -j2 $@ $^ -I$(top_builddir)/include
RivetALICEAnalyses.so: $(filter-out $(srcdir)/pluginALICE/tmp*.cc, $(wildcard $(srcdir)/pluginALICE/*.cc))
RivetATLASAnalyses.so: $(filter-out $(srcdir)/pluginATLAS/tmp*.cc, $(wildcard $(srcdir)/pluginATLAS/*.cc))
RivetBABARAnalyses.so: $(filter-out $(srcdir)/pluginBABAR/tmp*.cc, $(wildcard $(srcdir)/pluginBABAR/*.cc))
RivetBELLEAnalyses.so: $(filter-out $(srcdir)/pluginBELLE/tmp*.cc, $(wildcard $(srcdir)/pluginBELLE/*.cc))
RivetBESAnalyses.so: $(filter-out $(srcdir)/pluginBES/tmp*.cc, $(wildcard $(srcdir)/pluginBES/*.cc))
RivetCDFAnalyses.so: $(filter-out $(srcdir)/pluginCDF/tmp*.cc, $(wildcard $(srcdir)/pluginCDF/*.cc))
RivetCMSAnalyses.so: $(filter-out $(srcdir)/pluginCMS/tmp*.cc, $(wildcard $(srcdir)/pluginCMS/*.cc))
RivetD0Analyses.so: $(filter-out $(srcdir)/pluginD0/tmp*.cc, $(wildcard $(srcdir)/pluginD0/*.cc))
+RivetFrascatiAnalyses.so: $(filter-out $(srcdir)/pluginFrascati/tmp*.cc, $(wildcard $(srcdir)/pluginFrascati/*.cc))
RivetHERAAnalyses.so: $(filter-out $(srcdir)/pluginHERA/tmp*.cc, $(wildcard $(srcdir)/pluginHERA/*.cc))
RivetLEPAnalyses.so: $(filter-out $(srcdir)/pluginLEP/tmp*.cc, $(wildcard $(srcdir)/pluginLEP/*.cc))
RivetLHCbAnalyses.so: $(filter-out $(srcdir)/pluginLHCb/tmp*.cc, $(wildcard $(srcdir)/pluginLHCb/*.cc))
RivetLHCfAnalyses.so: $(filter-out $(srcdir)/pluginLHCf/tmp*.cc, $(wildcard $(srcdir)/pluginLHCf/*.cc))
RivetMCAnalyses.so: $(filter-out $(srcdir)/pluginMC/tmp*.cc, $(wildcard $(srcdir)/pluginMC/*.cc))
RivetNovosibirskAnalyses.so: $(filter-out $(srcdir)/pluginNovosibirsk/tmp*.cc, $(wildcard $(srcdir)/pluginNovosibirsk/*.cc))
RivetOrsayAnalyses.so: $(filter-out $(srcdir)/pluginOrsay/tmp*.cc, $(wildcard $(srcdir)/pluginOrsay/*.cc))
RivetMiscAnalyses.so: $(filter-out $(srcdir)/pluginMisc/tmp*.cc, $(wildcard $(srcdir)/pluginMisc/*.cc))
RivetPetraAnalyses.so: $(filter-out $(srcdir)/pluginPetra/tmp*.cc, $(wildcard $(srcdir)/pluginPetra/*.cc))
RivetRHICAnalyses.so: $(filter-out $(srcdir)/pluginRHIC/tmp*.cc, $(wildcard $(srcdir)/pluginRHIC/*.cc))
RivetSPSAnalyses.so: $(filter-out $(srcdir)/pluginSPS/tmp*.cc, $(wildcard $(srcdir)/pluginSPS/*.cc))
RivetTOTEMAnalyses.so: $(filter-out $(srcdir)/pluginTOTEM/tmp*.cc, $(wildcard $(srcdir)/pluginTOTEM/*.cc))
RivetTristanAnalyses.so: $(filter-out $(srcdir)/pluginTristan/tmp*.cc, $(wildcard $(srcdir)/pluginTristan/*.cc))
all-local: $(PLUGIN_LIBS) $(PLUGIN_DATAFILES)
mkdir -p $(builddir)/data
$(LN_S) -f $(abs_srcdir)/plugin*/*.{info,plot,yoda} $(builddir)/data
clean-local:
rm -rf data
install-exec-local: $(PLUGIN_LIBS)
$(INSTALL) $^ $(DESTDIR)$(libdir)
install-data-local: $(PLUGIN_DATAFILES)
@echo "Installing analysis data files..."
$(mkdir_p) $(DESTDIR)$(pkgdatadir)
cp $(abs_srcdir)/plugin*/*.{info,plot,yoda} $(DESTDIR)$(pkgdatadir)/
#@$(install_sh_DATA) $^ $(DESTDIR)$(pkgdatadir)
uninstall-local:
cd $(DESTDIR)$(libdir) && rm -f $(PLUGIN_LIBS)
@echo "Uninstalling analysis data files..."
rm -f $(DESTDIR)$(pkgdatadir)/*.{info,plot,yoda}
endif
diff --git a/analyses/pluginFrascati/BBAR_1980_I152630.cc b/analyses/pluginFrascati/BBAR_1980_I152630.cc
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/BBAR_1980_I152630.cc
@@ -0,0 +1,99 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/Projections/FinalState.hh"
+
+namespace Rivet {
+
+
+ /// @brief Add a short analysis description here
+ class BBAR_1980_I152630 : public Analysis {
+ public:
+
+ /// Constructor
+ DEFAULT_RIVET_ANALYSIS_CTOR(BBAR_1980_I152630);
+
+
+ /// @name Analysis methods
+ //@{
+
+ /// Book histograms and initialise projections before the run
+ void init() {
+ // Initialise and register projections
+ declare(FinalState(), "FS");
+
+ // Book histograms
+ _c_hadrons = bookCounter("/TMP/sigma_hadrons");
+ _c_muons = bookCounter("/TMP/sigma_muons");
+ }
+
+
+ /// Perform the per-event analysis
+ void analyze(const Event& event) {
+ const FinalState& fs = apply<FinalState>(event, "FS");
+
+ map<long,int> nCount;
+ int ntotal(0);
+ foreach (const Particle& p, fs.particles()) {
+ nCount[p.pdgId()] += 1;
+ ++ntotal;
+ }
+ // mu+mu- + photons
+ if(nCount[-13]==1 and nCount[13]==1 &&
+ ntotal==2+nCount[22])
+ _c_muons->fill(event.weight());
+ // everything else
+ else
+ _c_hadrons->fill(event.weight());
+ }
+
+
+ /// Normalise histograms etc., after the run
+ void finalize() {
+ Scatter1D R = *_c_hadrons/ *_c_muons;
+ double rval = R.point(0).x();
+ pair<double,double> rerr = R.point(0).xErrs();
+ double fact = crossSection()/ sumOfWeights() /picobarn;
+ double sig_h = _c_hadrons->val()*fact;
+ double err_h = _c_hadrons->err()*fact;
+ double sig_m = _c_muons ->val()*fact;
+ double err_m = _c_muons ->err()*fact;
+ Scatter2D temphisto(refData(1, 1, 1));
+ Scatter2DPtr hadrons = bookScatter2D("sigma_hadrons");
+ Scatter2DPtr muons = bookScatter2D("sigma_muons" );
+ Scatter2DPtr mult = bookScatter2D(1, 1, 1);
+ for (size_t b = 0; b < temphisto.numPoints(); b++) {
+ const double x = temphisto.point(b).x();
+ pair<double,double> ex = temphisto.point(b).xErrs();
+ pair<double,double> ex2 = ex;
+ if(ex2.first ==0.) ex2. first=0.0001;
+ if(ex2.second==0.) ex2.second=0.0001;
+ if (inRange(sqrtS()/GeV, x-ex2.first, x+ex2.second)) {
+ mult ->addPoint(x, rval, ex, rerr);
+ hadrons->addPoint(x, sig_h, ex, make_pair(err_h,err_h));
+ muons ->addPoint(x, sig_m, ex, make_pair(err_m,err_m));
+ }
+ else {
+ mult ->addPoint(x, 0., ex, make_pair(0.,.0));
+ hadrons->addPoint(x, 0., ex, make_pair(0.,.0));
+ muons ->addPoint(x, 0., ex, make_pair(0.,.0));
+ }
+ }
+ }
+
+ //@}
+
+
+ /// @name Histograms
+ //@{
+ CounterPtr _c_hadrons, _c_muons;
+ //@}
+
+
+ };
+
+
+ // The hook for the plugin system
+ DECLARE_RIVET_PLUGIN(BBAR_1980_I152630);
+
+
+}
diff --git a/analyses/pluginFrascati/BBAR_1980_I152630.info b/analyses/pluginFrascati/BBAR_1980_I152630.info
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/BBAR_1980_I152630.info
@@ -0,0 +1,31 @@
+Name: BBAR_1980_I152630
+Year: 1980
+Summary: Measurement of $R$ between 1.427 and 1.946
+Experiment: BBAR
+Collider: ADONE
+InspireID: 152630
+Status: VALIDATED
+Authors:
+ - Peter Richardson <peter.richardson@durham.ac.uk>
+References:
+ - Phys.Lett. B91 (1980) 155-160, 1980
+RunInfo: e+ e- to hadrons and e+ e- to mu+ mu- (for normalization)
+NeedCrossSection: no
+Beams: [e-, e+]
+Description:
+ 'Measurement of $R$ in $e^+e^-$ collisions by the baryon-antibaryon experiment between 1.427 and 1.946 GeV.
+ The individual hadronic and muonic cross sections are also outputted to the yoda file
+ so that ratio $R$ can be recalcuated if runs are combined.'
+BibKey: Ambrosio:1980mf
+BibTeX: '@article{Ambrosio:1980mf,
+ author = "Ambrosio, M. and others",
+ title = "{Total Cross-section for Hadron Production by Electron -
+ Positron Annihilation at Adone}",
+ journal = "Phys. Lett.",
+ volume = "91B",
+ year = "1980",
+ pages = "155-160",
+ doi = "10.1016/0370-2693(80)90683-8",
+ reportNumber = "LNF-80-4-P",
+ SLACcitation = "%%CITATION = PHLTA,91B,155;%%"
+}'
diff --git a/analyses/pluginFrascati/BBAR_1980_I152630.plot b/analyses/pluginFrascati/BBAR_1980_I152630.plot
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/BBAR_1980_I152630.plot
@@ -0,0 +1,21 @@
+BEGIN PLOT /BBAR_1980_I152630/d01-x01-y01
+Title=$R=\sigma(e^+e^-\to \text{hadrons})/\sigma(e^+e^-\to \mu^+\mu^-)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$R$
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /BBAR_1980_I152630/sigma_hadrons
+Title=$\sigma(e^+e^-\to \text{hadrons})$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to \text{hadrons})/pb$
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /BBAR_1980_I152630/sigma_muons
+Title=$\sigma(e^+e^-\to \mu^+\mu^-)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to \mu^+\mu^-)/pb$
+LogY=0
+ConnectGaps=1
+END PLOT
diff --git a/analyses/pluginFrascati/BBAR_1980_I152630.yoda b/analyses/pluginFrascati/BBAR_1980_I152630.yoda
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/BBAR_1980_I152630.yoda
@@ -0,0 +1,27 @@
+BEGIN YODA_SCATTER2D_V2 /REF/BBAR_1980_I152630/d01-x01-y01
+IsRef: 1
+Path: /REF/BBAR_1980_I152630/d01-x01-y01
+Title: doi:10.17182/hepdata.27239.v1/t1
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.427000e+00 1.500000e-02 1.000000e-02 1.840000e+00 3.190235e-01 3.190235e-01
+1.452000e+00 1.500000e-02 1.000000e-02 1.800000e+00 2.996665e-01 2.996665e-01
+1.472000e+00 1.000000e-02 1.500000e-02 2.410000e+00 3.953255e-01 3.953255e-01
+1.501000e+00 1.400000e-02 1.100000e-02 2.390000e+00 3.696245e-01 3.696245e-01
+1.523000e+00 1.100000e-02 1.400000e-02 2.660000e+00 4.508891e-01 4.508891e-01
+1.550000e+00 1.300000e-02 1.200000e-02 1.990000e+00 3.768319e-01 3.768319e-01
+1.574000e+00 1.200000e-02 1.300000e-02 1.730000e+00 3.276282e-01 3.276282e-01
+1.599000e+00 1.200000e-02 1.300000e-02 1.970000e+00 4.070875e-01 4.070875e-01
+1.624000e+00 1.200000e-02 1.300000e-02 1.940000e+00 3.902320e-01 3.902320e-01
+1.650000e+00 1.300000e-02 1.200000e-02 2.470000e+00 4.469567e-01 4.469567e-01
+1.673000e+00 1.100000e-02 1.400000e-02 2.100000e+00 3.900320e-01 3.900320e-01
+1.700000e+00 1.300000e-02 1.200000e-02 2.280000e+00 4.178086e-01 4.178086e-01
+1.722000e+00 1.000000e-02 1.500000e-02 2.310000e+00 4.518432e-01 4.518432e-01
+1.746000e+00 9.000000e-03 1.600000e-02 2.370000e+00 5.059449e-01 5.059449e-01
+1.799000e+00 1.200000e-02 1.300000e-02 2.440000e+00 4.322684e-01 4.322684e-01
+1.826000e+00 1.400000e-02 1.100000e-02 2.710000e+00 4.720617e-01 4.720617e-01
+1.845000e+00 8.000000e-03 1.700000e-02 2.820000e+00 6.473245e-01 6.473245e-01
+1.874000e+00 1.200000e-02 1.300000e-02 1.460000e+00 4.213799e-01 4.213799e-01
+1.946000e+00 9.000000e-03 1.600000e-02 1.410000e+00 3.509021e-01 3.509021e-01
+END YODA_SCATTER2D_V2
diff --git a/analyses/pluginFrascati/FENICE_1994_I377833.cc b/analyses/pluginFrascati/FENICE_1994_I377833.cc
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/FENICE_1994_I377833.cc
@@ -0,0 +1,82 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/Projections/FinalState.hh"
+#include "Rivet/Projections/UnstableParticles.hh"
+
+namespace Rivet {
+
+
+ /// @brief Add a short analysis description here
+ class FENICE_1994_I377833 : public Analysis {
+ public:
+
+ /// Constructor
+ DEFAULT_RIVET_ANALYSIS_CTOR(FENICE_1994_I377833);
+
+
+ /// @name Analysis methods
+ //@{
+
+ /// Book histograms and initialise projections before the run
+ void init() {
+ // Initialise and register projections
+ declare(FinalState(), "FS");
+ declare(UnstableParticles(), "UFS");
+ _nProton= bookCounter( "/TMP/nProton" );
+ }
+
+
+ /// Perform the per-event analysis
+ void analyze(const Event& event) {
+ const FinalState& fs = apply<FinalState>(event, "FS");
+ // total hadronic and muonic cross sections
+ map<long,int> nCount;
+ int ntotal(0);
+ foreach (const Particle& p, fs.particles()) {
+ nCount[p.pdgId()] += 1;
+ ++ntotal;
+ }
+ if(ntotal==2 && nCount[2212]==1 && nCount[-2212]==1)
+ _nProton->fill(event.weight());
+ }
+
+
+ /// Normalise histograms etc., after the run
+ void finalize() {
+ double fact = crossSection()/ sumOfWeights() /nanobarn;
+ double sigma = _nProton->val()*fact;
+ double error = _nProton->err()*fact;
+ Scatter2D temphisto(refData(1, 1, 1));
+ Scatter2DPtr mult = bookScatter2D(1, 1, 1);
+ for (size_t b = 0; b < temphisto.numPoints(); b++) {
+ const double x = temphisto.point(b).x();
+ pair<double,double> ex = temphisto.point(b).xErrs();
+ pair<double,double> ex2 = ex;
+ if(ex2.first ==0.) ex2. first=0.0001;
+ if(ex2.second==0.) ex2.second=0.0001;
+ if (inRange(sqr(sqrtS()/GeV), x-ex2.first, x+ex2.second)) {
+ mult->addPoint(x, sigma, ex, make_pair(error,error));
+ }
+ else {
+ mult->addPoint(x, 0., ex, make_pair(0.,.0));
+ }
+ }
+ }
+
+ //@}
+
+
+ /// @name Histograms
+ //@{
+ CounterPtr _nProton;
+ //@}
+
+
+ };
+
+
+ // The hook for the plugin system
+ DECLARE_RIVET_PLUGIN(FENICE_1994_I377833);
+
+
+}
diff --git a/analyses/pluginFrascati/FENICE_1994_I377833.info b/analyses/pluginFrascati/FENICE_1994_I377833.info
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/FENICE_1994_I377833.info
@@ -0,0 +1,30 @@
+Name: FENICE_1994_I377833
+Year: 1994
+Summary: Cross section for $e^+e^-\to p\bar{p}$ between 1.9 and 2.44 GeV
+Experiment: FENICE
+Collider: ADONE
+InspireID: 377833
+Status: VALIDATED
+Authors:
+ - Peter Richardson <peter.richardson@durham.ac.uk>
+References:
+ - Phys.Lett. B334 (1994) 431-434, 1994
+RunInfo: e+e- to hadrons
+NeedCrossSection: yes
+Beams: [e+,e-]
+Description:
+ 'Measurement of the cross section for $e^+e^-\to p\bar{p}$ between 1.9 and 2.44 GeV.'
+Keywords:
+BibKey: Antonelli:1994kq
+BibTeX: '@article{Antonelli:1994kq,
+ author = "Antonelli, A. and others",
+ title = "{Measurement of the electromagnetic form-factor of the
+ proton in the timelike region}",
+ journal = "Phys. Lett.",
+ volume = "B334",
+ year = "1994",
+ pages = "431-434",
+ doi = "10.1016/0370-2693(94)90710-2",
+ reportNumber = "LNF-94-032-P",
+ SLACcitation = "%%CITATION = PHLTA,B334,431;%%"
+}'
diff --git a/analyses/pluginFrascati/FENICE_1994_I377833.plot b/analyses/pluginFrascati/FENICE_1994_I377833.plot
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/FENICE_1994_I377833.plot
@@ -0,0 +1,8 @@
+BEGIN PLOT /FENICE_1994_I377833/d01-x01-y01
+Title=$\sigma(e^+e^-\to p\bar{p})$
+XLabel=$s/\text{GeV}^2$
+YLabel=$\sigma(e^+e^-\to p\bar{p})$/nb
+ConnectGaps=1
+XMin=3.5
+XMax=6.0
+END PLOT
diff --git a/analyses/pluginFrascati/FENICE_1994_I377833.yoda b/analyses/pluginFrascati/FENICE_1994_I377833.yoda
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/FENICE_1994_I377833.yoda
@@ -0,0 +1,24 @@
+BEGIN YODA_SCATTER2D_V2 /REF/FENICE_1994_I377833/d01-x01-y01
+IsRef: 1
+Path: /REF/FENICE_1994_I377833/d01-x01-y01
+Title: doi:10.17182/hepdata.28572.v1/t1
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+3.690000e+00 0.000000e+00 0.000000e+00 9.600000e-01 2.731300e-01 2.731300e-01
+4.000000e+00 0.000000e+00 0.000000e+00 6.200000e-01 1.708801e-01 1.708801e-01
+4.410000e+00 0.000000e+00 0.000000e+00 6.200000e-01 1.208305e-01 1.208305e-01
+5.950000e+00 0.000000e+00 0.000000e+00 2.800000e-01 1.208305e-01 1.208305e-01
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/FENICE_1994_I377833/d02-x01-y01
+IsRef: 1
+Path: /REF/FENICE_1994_I377833/d02-x01-y01
+Title: doi:10.17182/hepdata.28572.v1/t2
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+3.690000e+00 0.000000e+00 0.000000e+00 3.600000e-01 5.385165e-02 5.385165e-02
+4.000000e+00 0.000000e+00 0.000000e+00 2.400000e-01 3.162278e-02 3.162278e-02
+4.410000e+00 0.000000e+00 0.000000e+00 2.200000e-01 2.236068e-02 2.236068e-02
+5.950000e+00 0.000000e+00 0.000000e+00 1.500000e-01 3.162278e-02 3.162278e-02
+END YODA_SCATTER2D_V2
diff --git a/analyses/pluginFrascati/FENICE_1998_I471263.cc b/analyses/pluginFrascati/FENICE_1998_I471263.cc
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/FENICE_1998_I471263.cc
@@ -0,0 +1,82 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/Projections/FinalState.hh"
+#include "Rivet/Projections/UnstableParticles.hh"
+
+namespace Rivet {
+
+
+ /// @brief Add a short analysis description here
+ class FENICE_1998_I471263 : public Analysis {
+ public:
+
+ /// Constructor
+ DEFAULT_RIVET_ANALYSIS_CTOR(FENICE_1998_I471263);
+
+
+ /// @name Analysis methods
+ //@{
+
+ /// Book histograms and initialise projections before the run
+ void init() {
+ // Initialise and register projections
+ declare(FinalState(), "FS");
+ declare(UnstableParticles(), "UFS");
+ _nNeutron= bookCounter( "/TMP/nNeutron" );
+ }
+
+
+ /// Perform the per-event analysis
+ void analyze(const Event& event) {
+ const FinalState& fs = apply<FinalState>(event, "FS");
+ // total hadronic and muonic cross sections
+ map<long,int> nCount;
+ int ntotal(0);
+ foreach (const Particle& p, fs.particles()) {
+ nCount[p.pdgId()] += 1;
+ ++ntotal;
+ }
+ if(ntotal==2 && nCount[2112]==1 && nCount[-2112]==1)
+ _nNeutron->fill(event.weight());
+ }
+
+
+ /// Normalise histograms etc., after the run
+ void finalize() {
+ double fact = crossSection()/ sumOfWeights() /nanobarn;
+ double sigma = _nNeutron->val()*fact;
+ double error = _nNeutron->err()*fact;
+ Scatter2D temphisto(refData(1, 1, 1));
+ Scatter2DPtr mult = bookScatter2D(1, 1, 1);
+ for (size_t b = 0; b < temphisto.numPoints(); b++) {
+ const double x = temphisto.point(b).x();
+ pair<double,double> ex = temphisto.point(b).xErrs();
+ pair<double,double> ex2 = ex;
+ if(ex2.first ==0.) ex2. first=0.0001;
+ if(ex2.second==0.) ex2.second=0.0001;
+ if (inRange(sqrtS()/GeV, x-ex2.first, x+ex2.second)) {
+ mult->addPoint(x, sigma, ex, make_pair(error,error));
+ }
+ else {
+ mult->addPoint(x, 0., ex, make_pair(0.,.0));
+ }
+ }
+ }
+
+ //@}
+
+
+ /// @name Histograms
+ //@{
+ CounterPtr _nNeutron;
+ //@}
+
+
+ };
+
+
+ // The hook for the plugin system
+ DECLARE_RIVET_PLUGIN(FENICE_1998_I471263);
+
+
+}
diff --git a/analyses/pluginFrascati/FENICE_1998_I471263.info b/analyses/pluginFrascati/FENICE_1998_I471263.info
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/FENICE_1998_I471263.info
@@ -0,0 +1,29 @@
+Name: FENICE_1998_I471263
+Year: 1998
+Summary: Cross section for $e^+e^-\to n\bar{n}$ between 1.9 and 2.44 GeV
+Experiment: FENICE
+Collider: ADONE
+InspireID: 471263
+Status: VALIDATED
+Authors:
+ - Peter Richardson <peter.richardson@durham.ac.uk>
+References:
+ - Nucl.Phys. B517 (1998) 3-35, 1998
+RunInfo: e+e- to hadrons
+NeedCrossSection: yes
+Beams: [e+,e-]
+Description:
+ 'Measurement of the cross section for $e^+e^-\to n\bar{n}$ between 1.9 and 2.44 GeV.'
+Keywords:
+BibKey: Antonelli:1998fv
+BibTeX: '@article{Antonelli:1998fv,
+ author = "Antonelli, A. and others",
+ title = "{The first measurement of the neutron electromagnetic
+ form-factors in the timelike region}",
+ journal = "Nucl. Phys.",
+ volume = "B517",
+ year = "1998",
+ pages = "3-35",
+ doi = "10.1016/S0550-3213(98)00083-2",
+ SLACcitation = "%%CITATION = NUPHA,B517,3;%%"
+}'
diff --git a/analyses/pluginFrascati/FENICE_1998_I471263.plot b/analyses/pluginFrascati/FENICE_1998_I471263.plot
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/FENICE_1998_I471263.plot
@@ -0,0 +1,8 @@
+BEGIN PLOT /FENICE_1998_I471263/d01-x01-y01
+Title=$\sigma(e^+e^-\to n\bar{n})$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to n\bar{n})$/nb
+ConnectGaps=1
+XMin=1.8
+XMax=2.5
+END PLOT
diff --git a/analyses/pluginFrascati/FENICE_1998_I471263.yoda b/analyses/pluginFrascati/FENICE_1998_I471263.yoda
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/FENICE_1998_I471263.yoda
@@ -0,0 +1,65 @@
+BEGIN YODA_SCATTER2D_V2 /REF/FENICE_1998_I471263/d01-x01-y01
+IsRef: 1
+Path: /REF/FENICE_1998_I471263/d01-x01-y01
+Title: doi:10.17182/hepdata.32681.v1/t1
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.900000e+00 0.000000e+00 0.000000e+00 1.100000e+00 8.000000e-01 8.000000e-01
+1.920000e+00 0.000000e+00 0.000000e+00 1.060000e+00 4.400000e-01 4.400000e-01
+2.000000e+00 0.000000e+00 0.000000e+00 8.800000e-01 3.400000e-01 3.400000e-01
+2.100000e+00 0.000000e+00 0.000000e+00 1.000000e+00 2.800000e-01 2.800000e-01
+2.440000e+00 0.000000e+00 0.000000e+00 6.500000e-01 2.900000e-01 2.900000e-01
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/FENICE_1998_I471263/d02-x01-y01
+IsRef: 1
+Path: /REF/FENICE_1998_I471263/d02-x01-y01
+Title: doi:10.17182/hepdata.32681.v1/t2
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+3.640000e+00 0.000000e+00 0.000000e+00 4.500000e-01 1.600000e-01 1.600000e-01
+3.690000e+00 0.000000e+00 0.000000e+00 3.900000e-01 8.000000e-02 8.000000e-02
+4.000000e+00 0.000000e+00 0.000000e+00 2.900000e-01 5.000000e-02 5.000000e-02
+4.410000e+00 0.000000e+00 0.000000e+00 2.800000e-01 4.000000e-02 4.000000e-02
+5.950000e+00 0.000000e+00 0.000000e+00 2.300000e-01 5.000000e-02 5.000000e-02
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/FENICE_1998_I471263/d02-x01-y02
+IsRef: 1
+Path: /REF/FENICE_1998_I471263/d02-x01-y02
+Title: doi:10.17182/hepdata.32681.v1/t2
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+3.640000e+00 0.000000e+00 0.000000e+00 5.500000e-01 2.000000e-01 2.000000e-01
+3.690000e+00 0.000000e+00 0.000000e+00 4.700000e-01 1.000000e-01 1.000000e-01
+4.000000e+00 0.000000e+00 0.000000e+00 3.400000e-01 6.000000e-02 6.000000e-02
+4.410000e+00 0.000000e+00 0.000000e+00 3.400000e-01 5.000000e-02 5.000000e-02
+5.950000e+00 0.000000e+00 0.000000e+00 2.600000e-01 6.000000e-02 6.000000e-02
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/FENICE_1998_I471263/d03-x01-y01
+IsRef: 1
+Path: /REF/FENICE_1998_I471263/d03-x01-y01
+Title: doi:10.17182/hepdata.32681.v1/t3
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+3.610000e+00 0.000000e+00 0.000000e+00 9.700000e-01 4.400000e-01 5.300000e-01
+3.690000e+00 0.000000e+00 0.000000e+00 9.600000e-01 2.700000e-01 2.700000e-01
+4.000000e+00 0.000000e+00 0.000000e+00 6.200000e-01 1.700000e-01 1.700000e-01
+4.420000e+00 0.000000e+00 0.000000e+00 6.200000e-01 1.200000e-01 1.200000e-01
+5.950000e+00 0.000000e+00 0.000000e+00 2.800000e-01 1.200000e-01 1.200000e-01
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/FENICE_1998_I471263/d03-x01-y02
+IsRef: 1
+Path: /REF/FENICE_1998_I471263/d03-x01-y02
+Title: doi:10.17182/hepdata.32681.v1/t3
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+3.610000e+00 0.000000e+00 0.000000e+00 4.100000e-01 9.000000e-02 1.100000e-01
+3.690000e+00 0.000000e+00 0.000000e+00 3.600000e-01 5.000000e-02 5.000000e-02
+4.000000e+00 0.000000e+00 0.000000e+00 2.400000e-01 3.000000e-02 3.000000e-02
+4.420000e+00 0.000000e+00 0.000000e+00 2.200000e-01 2.000000e-02 2.000000e-02
+5.950000e+00 0.000000e+00 0.000000e+00 1.500000e-01 3.000000e-02 3.000000e-02
+END YODA_SCATTER2D_V2
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1973_I84794.cc b/analyses/pluginFrascati/GAMMAGAMMA_1973_I84794.cc
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1973_I84794.cc
@@ -0,0 +1,99 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/Projections/FinalState.hh"
+
+namespace Rivet {
+
+
+ /// @brief Add a short analysis description here
+ class GAMMAGAMMA_1973_I84794 : public Analysis {
+ public:
+
+ /// Constructor
+ DEFAULT_RIVET_ANALYSIS_CTOR(GAMMAGAMMA_1973_I84794);
+
+
+ /// @name Analysis methods
+ //@{
+
+ /// Book histograms and initialise projections before the run
+ void init() {
+ // Initialise and register projections
+ declare(FinalState(), "FS");
+
+ // Book histograms
+ _c_hadrons = bookCounter("/TMP/sigma_hadrons");
+ _c_muons = bookCounter("/TMP/sigma_muons");
+ }
+
+
+ /// Perform the per-event analysis
+ void analyze(const Event& event) {
+ const FinalState& fs = apply<FinalState>(event, "FS");
+
+ map<long,int> nCount;
+ int ntotal(0);
+ foreach (const Particle& p, fs.particles()) {
+ nCount[p.pdgId()] += 1;
+ ++ntotal;
+ }
+ // mu+mu- + photons
+ if(nCount[-13]==1 and nCount[13]==1 &&
+ ntotal==2+nCount[22])
+ _c_muons->fill(event.weight());
+ // everything else
+ else
+ _c_hadrons->fill(event.weight());
+ }
+
+
+ /// Normalise histograms etc., after the run
+ void finalize() {
+ Scatter1D R = *_c_hadrons/ *_c_muons;
+ double rval = R.point(0).x();
+ pair<double,double> rerr = R.point(0).xErrs();
+ double fact = crossSection()/ sumOfWeights() /nanobarn;
+ double sig_h = _c_hadrons->val()*fact;
+ double err_h = _c_hadrons->err()*fact;
+ double sig_m = _c_muons ->val()*fact;
+ double err_m = _c_muons ->err()*fact;
+ Scatter2D temphisto(refData(3, 1, 1));
+ Scatter2DPtr hadrons = bookScatter2D(3,1,1);
+ Scatter2DPtr muons = bookScatter2D("sigma_muons" );
+ Scatter2DPtr mult = bookScatter2D(4,1,1);
+ for (size_t b = 0; b < temphisto.numPoints(); b++) {
+ const double x = temphisto.point(b).x();
+ pair<double,double> ex = temphisto.point(b).xErrs();
+ pair<double,double> ex2 = ex;
+ if(ex2.first ==0.) ex2. first=0.0001;
+ if(ex2.second==0.) ex2.second=0.0001;
+ if (inRange(sqrtS()/GeV, x-ex2.first, x+ex2.second)) {
+ mult ->addPoint(x, rval, ex, rerr);
+ hadrons->addPoint(x, sig_h, ex, make_pair(err_h,err_h));
+ muons ->addPoint(x, sig_m, ex, make_pair(err_m,err_m));
+ }
+ else {
+ mult ->addPoint(x, 0., ex, make_pair(0.,.0));
+ hadrons->addPoint(x, 0., ex, make_pair(0.,.0));
+ muons ->addPoint(x, 0., ex, make_pair(0.,.0));
+ }
+ }
+ }
+
+ //@}
+
+
+ /// @name Histograms
+ //@{
+ CounterPtr _c_hadrons, _c_muons;
+ //@}
+
+
+ };
+
+
+ // The hook for the plugin system
+ DECLARE_RIVET_PLUGIN(GAMMAGAMMA_1973_I84794);
+
+
+}
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1973_I84794.info b/analyses/pluginFrascati/GAMMAGAMMA_1973_I84794.info
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1973_I84794.info
@@ -0,0 +1,32 @@
+Name: GAMMAGAMMA_1973_I84794
+Year: 1973
+Summary: Measurement of $R$ and the hadronic cross section for energies between 1.35 and 3 GeV
+Experiment: GAMMAGAMMA
+Collider: ADONE
+InspireID: 84794
+Status: VALIDATED
+Authors:
+ - Peter Richardson <peter.richardson@durham.ac.uk>
+References:
+ - Phys.Lett. B44 (1973) 533-536, 1973
+RunInfo: e+ e- to hadrons and e+ e- to mu+ mu- (for normalization)
+NeedCrossSection: yes
+Beams: [e-, e+]
+Description:
+ 'Measurement of $R$ and the hadronic cross section in $e^+e^-$ collisions by the GAMMA-GAMMA group for energies between 1.35 and 3 GeV.
+ The muonic cross section is also outputted to the yoda file
+ so that ratio $R$ can be recalcuated if runs are combined.'
+BibKey: Bacci:1973fb
+BibTeX: '@article{Bacci:1973fb,
+ author = "Bacci, C. and Penson, G. and Salvini, G. and Stella, B.
+ and Baldini-Celio, R. and Capon, G. and Mencuccini, C. and
+ Murtas, G. P. and Spinetti, M. and Zallo, A.",
+ title = "{Multihadronic cross-sections from e+ e- annihilation up
+ to 3 gev center-of-mass energy}",
+ journal = "Phys. Lett.",
+ volume = "44B",
+ year = "1973",
+ pages = "533-536",
+ doi = "10.1016/0370-2693(73)90016-6",
+ SLACcitation = "%%CITATION = PHLTA,44B,533;%%"
+}'
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1973_I84794.plot b/analyses/pluginFrascati/GAMMAGAMMA_1973_I84794.plot
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1973_I84794.plot
@@ -0,0 +1,21 @@
+BEGIN PLOT /GAMMAGAMMA_1973_I84794/d04-x01-y01
+Title=$R=\sigma(e^+e^-\to \text{hadrons})/\sigma(e^+e^-\to \mu^+\mu^-)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$R$
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /GAMMAGAMMA_1973_I84794/d03-x01-y01
+Title=$\sigma(e^+e^-\to \text{hadrons})$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to \text{hadrons})/nb$
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /GAMMAGAMMA_1973_I84794/sigma_muons
+Title=$\sigma(e^+e^-\to \mu^+\mu^-)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to \mu^+\mu^-)/nb$
+LogY=0
+ConnectGaps=1
+END PLOT
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1973_I84794.yoda b/analyses/pluginFrascati/GAMMAGAMMA_1973_I84794.yoda
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1973_I84794.yoda
@@ -0,0 +1,60 @@
+BEGIN YODA_SCATTER2D_V2 /REF/GAMMAGAMMA_1973_I84794/d01-x01-y01
+IsRef: 1
+Path: /REF/GAMMAGAMMA_1973_I84794/d01-x01-y01
+Title: doi:10.17182/hepdata.6496.v1/t1
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.350000e+00 0.000000e+00 0.000000e+00 4.500000e+01 1.800000e+01 1.800000e+01
+1.450000e+00 0.000000e+00 0.000000e+00 5.500000e+01 1.500000e+01 1.500000e+01
+1.650000e+00 0.000000e+00 0.000000e+00 3.600000e+01 7.000000e+00 7.000000e+00
+1.900000e+00 0.000000e+00 0.000000e+00 1.550000e+01 4.000000e+00 4.000000e+00
+1.980000e+00 0.000000e+00 0.000000e+00 3.000000e+01 1.000000e+01 1.000000e+01
+2.100000e+00 0.000000e+00 0.000000e+00 1.700000e+01 4.000000e+00 4.000000e+00
+2.400000e+00 0.000000e+00 0.000000e+00 3.400000e+01 1.200000e+01 1.200000e+01
+2.800000e+00 0.000000e+00 0.000000e+00 1.500000e+01 3.500000e+00 3.500000e+00
+3.000000e+00 0.000000e+00 0.000000e+00 2.800000e+01 4.500000e+00 4.500000e+00
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/GAMMAGAMMA_1973_I84794/d02-x01-y01
+IsRef: 1
+Path: /REF/GAMMAGAMMA_1973_I84794/d02-x01-y01
+Title: doi:10.17182/hepdata.6496.v1/t2
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+3.000000e+00 0.000000e+00 0.000000e+00 5.000000e+00 3.000000e+00 3.000000e+00
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/GAMMAGAMMA_1973_I84794/d03-x01-y01
+IsRef: 1
+Path: /REF/GAMMAGAMMA_1973_I84794/d03-x01-y01
+Title: doi:10.17182/hepdata.6496.v1/t3
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.350000e+00 0.000000e+00 0.000000e+00 4.500000e+01 2.163331e+01 2.163331e+01
+1.450000e+00 0.000000e+00 0.000000e+00 5.500000e+01 2.051828e+01 2.051828e+01
+1.650000e+00 0.000000e+00 0.000000e+00 3.600000e+01 1.140175e+01 1.140175e+01
+1.900000e+00 0.000000e+00 0.000000e+00 1.550000e+01 5.656854e+00 5.656854e+00
+1.980000e+00 0.000000e+00 0.000000e+00 3.000000e+01 1.280625e+01 1.280625e+01
+2.100000e+00 0.000000e+00 0.000000e+00 1.700000e+01 5.656854e+00 5.656854e+00
+2.400000e+00 0.000000e+00 0.000000e+00 3.400000e+01 1.500000e+01 1.500000e+01
+2.800000e+00 0.000000e+00 0.000000e+00 1.500000e+01 4.609772e+00 4.609772e+00
+3.000000e+00 0.000000e+00 0.000000e+00 2.800000e+01 8.746428e+00 8.746428e+00
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/GAMMAGAMMA_1973_I84794/d04-x01-y01
+IsRef: 1
+Path: /REF/GAMMAGAMMA_1973_I84794/d04-x01-y01
+Title: doi:10.17182/hepdata.6496.v1/t4
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.350000e+00 0.000000e+00 0.000000e+00 9.400000e-01 4.548626e-01 4.548626e-01
+1.450000e+00 0.000000e+00 0.000000e+00 1.330000e+00 4.951767e-01 4.951767e-01
+1.650000e+00 0.000000e+00 0.000000e+00 1.130000e+00 3.560899e-01 3.560899e-01
+1.900000e+00 0.000000e+00 0.000000e+00 6.400000e-01 2.404163e-01 2.404163e-01
+1.980000e+00 0.000000e+00 0.000000e+00 1.350000e+00 5.762812e-01 5.762812e-01
+2.100000e+00 0.000000e+00 0.000000e+00 8.600000e-01 2.828427e-01 2.828427e-01
+2.400000e+00 0.000000e+00 0.000000e+00 2.250000e+00 1.000000e+00 1.000000e+00
+2.800000e+00 0.000000e+00 0.000000e+00 1.350000e+00 4.186884e-01 4.186884e-01
+3.000000e+00 0.000000e+00 0.000000e+00 2.900000e+00 9.106591e-01 9.106591e-01
+END YODA_SCATTER2D_V2
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1975_I100016.cc b/analyses/pluginFrascati/GAMMAGAMMA_1975_I100016.cc
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1975_I100016.cc
@@ -0,0 +1,99 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/Projections/FinalState.hh"
+
+namespace Rivet {
+
+
+ /// @brief Add a short analysis description here
+ class GAMMAGAMMA_1975_I100016 : public Analysis {
+ public:
+
+ /// Constructor
+ DEFAULT_RIVET_ANALYSIS_CTOR(GAMMAGAMMA_1975_I100016);
+
+
+ /// @name Analysis methods
+ //@{
+
+ /// Book histograms and initialise projections before the run
+ void init() {
+ // Initialise and register projections
+ declare(FinalState(), "FS");
+
+ // Book histograms
+ _c_hadrons = bookCounter("/TMP/sigma_hadrons");
+ _c_muons = bookCounter("/TMP/sigma_muons");
+ }
+
+
+ /// Perform the per-event analysis
+ void analyze(const Event& event) {
+ const FinalState& fs = apply<FinalState>(event, "FS");
+
+ map<long,int> nCount;
+ int ntotal(0);
+ foreach (const Particle& p, fs.particles()) {
+ nCount[p.pdgId()] += 1;
+ ++ntotal;
+ }
+ // mu+mu- + photons
+ if(nCount[-13]==1 and nCount[13]==1 &&
+ ntotal==2+nCount[22])
+ _c_muons->fill(event.weight());
+ // everything else
+ else
+ _c_hadrons->fill(event.weight());
+ }
+
+
+ /// Normalise histograms etc., after the run
+ void finalize() {
+ Scatter1D R = *_c_hadrons/ *_c_muons;
+ double rval = R.point(0).x();
+ pair<double,double> rerr = R.point(0).xErrs();
+ double fact = crossSection()/ sumOfWeights() /nanobarn;
+ double sig_h = _c_hadrons->val()*fact;
+ double err_h = _c_hadrons->err()*fact;
+ double sig_m = _c_muons ->val()*fact;
+ double err_m = _c_muons ->err()*fact;
+ Scatter2D temphisto(refData(1, 1, 1));
+ Scatter2DPtr hadrons = bookScatter2D(1,1,1);
+ Scatter2DPtr muons = bookScatter2D("sigma_muons" );
+ Scatter2DPtr mult = bookScatter2D(1, 1, 2);
+ for (size_t b = 0; b < temphisto.numPoints(); b++) {
+ const double x = temphisto.point(b).x();
+ pair<double,double> ex = temphisto.point(b).xErrs();
+ pair<double,double> ex2 = ex;
+ if(ex2.first ==0.) ex2. first=0.0001;
+ if(ex2.second==0.) ex2.second=0.0001;
+ if (inRange(sqrtS()/GeV, x-ex2.first, x+ex2.second)) {
+ mult ->addPoint(x, rval, ex, rerr);
+ hadrons->addPoint(x, sig_h, ex, make_pair(err_h,err_h));
+ muons ->addPoint(x, sig_m, ex, make_pair(err_m,err_m));
+ }
+ else {
+ mult ->addPoint(x, 0., ex, make_pair(0.,.0));
+ hadrons->addPoint(x, 0., ex, make_pair(0.,.0));
+ muons ->addPoint(x, 0., ex, make_pair(0.,.0));
+ }
+ }
+ }
+
+ //@}
+
+
+ /// @name Histograms
+ //@{
+ CounterPtr _c_hadrons, _c_muons;
+ //@}
+
+
+ };
+
+
+ // The hook for the plugin system
+ DECLARE_RIVET_PLUGIN(GAMMAGAMMA_1975_I100016);
+
+
+}
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1975_I100016.info b/analyses/pluginFrascati/GAMMAGAMMA_1975_I100016.info
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1975_I100016.info
@@ -0,0 +1,32 @@
+Name: GAMMAGAMMA_1975_I100016
+Year: 1975
+Summary: Measurement of $R$ and the hadronic cross section for energies between 1.91 and 2.97 GeV
+Experiment: GAMMAGAMMA
+Collider: ADONE
+InspireID: 100016
+Status: VALIDATED
+Authors:
+ - Peter Richardson <peter.richardson@durham.ac.uk>
+References:
+ - Phys.Lett. B58 (1975) 481-483, 1975
+RunInfo: e+ e- to hadrons and e+ e- to mu+ mu- (for normalization)
+NeedCrossSection: yes
+Beams: [e-, e+]
+Description:
+ 'Measurement of $R$ and the hadronic cross section in $e^+e^-$ collisions by the GAMMA-GAMMA group for energies between 1.91 and 2.97 GeV.
+ The muonic cross section is also outputted to the yoda file
+ so that ratio $R$ can be recalcuated if runs are combined.'
+BibKey: Bacci:1975eb
+BibTeX: '@article{Bacci:1975eb,
+ author = "Bacci, C. and others",
+ title = "{Search for Narrow Resonances in e+ e- Annihilation Into
+ Hadrons in the Mass Regions 1910-MeV-2545-MeV and
+ 2970-MeV-3090-MeV}",
+ journal = "Phys. Lett.",
+ volume = "58B",
+ year = "1975",
+ pages = "481-483",
+ doi = "10.1016/0370-2693(75)90707-8",
+ reportNumber = "LNF-75/38-P",
+ SLACcitation = "%%CITATION = PHLTA,58B,481;%%"
+}'
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1975_I100016.plot b/analyses/pluginFrascati/GAMMAGAMMA_1975_I100016.plot
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1975_I100016.plot
@@ -0,0 +1,21 @@
+BEGIN PLOT /GAMMAGAMMA_1975_I100016/d01-x01-y02
+Title=$R=\sigma(e^+e^-\to \text{hadrons})/\sigma(e^+e^-\to \mu^+\mu^-)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$R$
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /GAMMAGAMMA_1975_I100016/d01-x01-y01
+Title=$\sigma(e^+e^-\to \text{hadrons})$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to \text{hadrons})/nb$
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /GAMMAGAMMA_1975_I100016/sigma_muons
+Title=$\sigma(e^+e^-\to \mu^+\mu^-)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to \mu^+\mu^-)/nb$
+LogY=0
+ConnectGaps=1
+END PLOT
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1975_I100016.yoda b/analyses/pluginFrascati/GAMMAGAMMA_1975_I100016.yoda
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1975_I100016.yoda
@@ -0,0 +1,22 @@
+BEGIN YODA_SCATTER2D_V2 /REF/GAMMAGAMMA_1975_I100016/d01-x01-y01
+IsRef: 1
+Path: /REF/GAMMAGAMMA_1975_I100016/d01-x01-y01
+Title: doi:10.17182/hepdata.27776.v1/t1
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+2.055000e+00 1.450000e-01 1.450000e-01 2.900000e+01 1.040829e+01 1.040829e+01
+2.372500e+00 1.725000e-01 1.725000e-01 3.000000e+01 1.050000e+01 1.050000e+01
+3.030000e+00 6.000000e-02 6.000000e-02 2.500000e+01 9.182184e+00 9.182184e+00
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/GAMMAGAMMA_1975_I100016/d01-x01-y02
+IsRef: 1
+Path: /REF/GAMMAGAMMA_1975_I100016/d01-x01-y02
+Title: doi:10.17182/hepdata.27776.v1/t1
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+2.055000e+00 1.450000e-01 1.450000e-01 1.400000e+00 4.729693e-01 4.729693e-01
+2.372500e+00 1.725000e-01 1.725000e-01 1.900000e+00 6.908871e-01 6.908871e-01
+3.030000e+00 6.000000e-02 6.000000e-02 2.600000e+00 9.272001e-01 9.272001e-01
+END YODA_SCATTER2D_V2
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1979_I133588.cc b/analyses/pluginFrascati/GAMMAGAMMA_1979_I133588.cc
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1979_I133588.cc
@@ -0,0 +1,99 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/Projections/FinalState.hh"
+
+namespace Rivet {
+
+
+ /// @brief Add a short analysis description here
+ class GAMMAGAMMA_1979_I133588 : public Analysis {
+ public:
+
+ /// Constructor
+ DEFAULT_RIVET_ANALYSIS_CTOR(GAMMAGAMMA_1979_I133588);
+
+
+ /// @name Analysis methods
+ //@{
+
+ /// Book histograms and initialise projections before the run
+ void init() {
+ // Initialise and register projections
+ declare(FinalState(), "FS");
+
+ // Book histograms
+ _c_hadrons = bookCounter("/TMP/sigma_hadrons");
+ _c_muons = bookCounter("/TMP/sigma_muons");
+ }
+
+
+ /// Perform the per-event analysis
+ void analyze(const Event& event) {
+ const FinalState& fs = apply<FinalState>(event, "FS");
+
+ map<long,int> nCount;
+ int ntotal(0);
+ foreach (const Particle& p, fs.particles()) {
+ nCount[p.pdgId()] += 1;
+ ++ntotal;
+ }
+ // mu+mu- + photons
+ if(nCount[-13]==1 and nCount[13]==1 &&
+ ntotal==2+nCount[22])
+ _c_muons->fill(event.weight());
+ // everything else
+ else
+ _c_hadrons->fill(event.weight());
+ }
+
+
+ /// Normalise histograms etc., after the run
+ void finalize() {
+ Scatter1D R = *_c_hadrons/ *_c_muons;
+ double rval = R.point(0).x();
+ pair<double,double> rerr = R.point(0).xErrs();
+ double fact = crossSection()/ sumOfWeights() /nanobarn;
+ double sig_h = _c_hadrons->val()*fact;
+ double err_h = _c_hadrons->err()*fact;
+ double sig_m = _c_muons ->val()*fact;
+ double err_m = _c_muons ->err()*fact;
+ Scatter2D temphisto(refData(1, 1, 1));
+ Scatter2DPtr hadrons = bookScatter2D(1,1,2);
+ Scatter2DPtr muons = bookScatter2D("sigma_muons" );
+ Scatter2DPtr mult = bookScatter2D(1, 1, 1);
+ for (size_t b = 0; b < temphisto.numPoints(); b++) {
+ const double x = temphisto.point(b).x();
+ pair<double,double> ex = temphisto.point(b).xErrs();
+ pair<double,double> ex2 = ex;
+ if(ex2.first ==0.) ex2. first=0.0001;
+ if(ex2.second==0.) ex2.second=0.0001;
+ if (inRange(sqrtS()/GeV, x-ex2.first, x+ex2.second)) {
+ mult ->addPoint(x, rval, ex, rerr);
+ hadrons->addPoint(x, sig_h, ex, make_pair(err_h,err_h));
+ muons ->addPoint(x, sig_m, ex, make_pair(err_m,err_m));
+ }
+ else {
+ mult ->addPoint(x, 0., ex, make_pair(0.,.0));
+ hadrons->addPoint(x, 0., ex, make_pair(0.,.0));
+ muons ->addPoint(x, 0., ex, make_pair(0.,.0));
+ }
+ }
+ }
+
+ //@}
+
+
+ /// @name Histograms
+ //@{
+ CounterPtr _c_hadrons, _c_muons;
+ //@}
+
+
+ };
+
+
+ // The hook for the plugin system
+ DECLARE_RIVET_PLUGIN(GAMMAGAMMA_1979_I133588);
+
+
+}
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1979_I133588.info b/analyses/pluginFrascati/GAMMAGAMMA_1979_I133588.info
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1979_I133588.info
@@ -0,0 +1,32 @@
+Name: GAMMAGAMMA_1979_I133588
+Year: 1979
+Summary: Measurement of $R$ and the hadronic cross section for energies between 1.92 and 3.09 GeV
+Experiment: GAMMAGAMMA
+Collider: ADONE
+InspireID: 133588
+Status: VALIDATED
+Authors:
+ - Peter Richardson <peter.richardson@durham.ac.uk>
+References:
+ - Lett.Nuovo Cim. 24 (1979) 324, 1979
+RunInfo: e+ e- to hadrons and e+ e- to mu+ mu- (for normalization)
+NeedCrossSection: yes
+Beams: [e-, e+]
+Description:
+ 'Measurement of $R$ and the hadronic cross section in $e^+e^-$ collisions by the GAMMA-GAMMA group for energies between 1.92 and 3.09 GeV.
+ The muonic cross section is also outputted to the yoda file
+ so that ratio $R$ can be recalcuated if runs are combined.'
+BibKey: Bacci:1978mq
+BibTeX: '@article{Bacci:1978mq,
+ author = "Bacci, C. and others",
+ title = "{Total Hadronic Cross-section From $e^+ e^-$ Annihilation
+ in the Total Center-of-mass Energy Range 1920-{MeV} -
+ 3090-{MeV}}",
+ journal = "Lett. Nuovo Cim.",
+ volume = "24",
+ year = "1979",
+ pages = "324",
+ doi = "10.1007/BF02724851",
+ reportNumber = "LNF-78/56-P",
+ SLACcitation = "%%CITATION = NCLTA,24,324;%%"
+}'
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1979_I133588.plot b/analyses/pluginFrascati/GAMMAGAMMA_1979_I133588.plot
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1979_I133588.plot
@@ -0,0 +1,21 @@
+BEGIN PLOT /GAMMAGAMMA_1979_I133588/d01-x01-y01
+Title=$R=\sigma(e^+e^-\to \text{hadrons})/\sigma(e^+e^-\to \mu^+\mu^-)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$R$
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /GAMMAGAMMA_1979_I133588/d01-x01-y02
+Title=$\sigma(e^+e^-\to \text{hadrons})$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to \text{hadrons})/nb$
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /GAMMAGAMMA_1979_I133588/sigma_muons
+Title=$\sigma(e^+e^-\to \mu^+\mu^-)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to \mu^+\mu^-)/nb$
+LogY=0
+ConnectGaps=1
+END PLOT
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1979_I133588.yoda b/analyses/pluginFrascati/GAMMAGAMMA_1979_I133588.yoda
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1979_I133588.yoda
@@ -0,0 +1,32 @@
+BEGIN YODA_SCATTER2D_V2 /REF/GAMMAGAMMA_1979_I133588/d01-x01-y01
+IsRef: 1
+Path: /REF/GAMMAGAMMA_1979_I133588/d01-x01-y01
+Title: doi:10.17182/hepdata.37911.v1/t1
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.960000e+00 4.000000e-02 4.000000e-02 1.740000e+00 1.700000e-01 1.700000e-01
+2.050000e+00 5.000000e-02 5.000000e-02 1.860000e+00 1.900000e-01 1.900000e-01
+2.150000e+00 5.000000e-02 5.000000e-02 2.090000e+00 1.300000e-01 1.300000e-01
+2.370000e+00 1.700000e-01 1.700000e-01 2.690000e+00 2.700000e-01 2.700000e-01
+2.590000e+00 5.000000e-02 5.000000e-02 2.890000e+00 2.900000e-01 2.900000e-01
+2.700000e+00 6.000000e-02 6.000000e-02 2.990000e+00 3.000000e-01 3.000000e-01
+2.870000e+00 1.100000e-01 1.100000e-01 2.380000e+00 1.900000e-01 1.900000e-01
+3.035000e+00 5.500000e-02 5.500000e-02 2.460000e+00 2.700000e-01 2.700000e-01
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/GAMMAGAMMA_1979_I133588/d01-x01-y02
+IsRef: 1
+Path: /REF/GAMMAGAMMA_1979_I133588/d01-x01-y02
+Title: doi:10.17182/hepdata.37911.v1/t1
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.960000e+00 4.000000e-02 4.000000e-02 3.930000e+01 3.800000e+00 3.800000e+00
+2.050000e+00 5.000000e-02 5.000000e-02 3.840000e+01 3.900000e+00 3.900000e+00
+2.150000e+00 5.000000e-02 5.000000e-02 3.930000e+01 2.400000e+00 2.400000e+00
+2.370000e+00 1.700000e-01 1.700000e-01 4.160000e+01 4.200000e+00 4.200000e+00
+2.590000e+00 5.000000e-02 5.000000e-02 3.740000e+01 3.800000e+00 3.800000e+00
+2.700000e+00 6.000000e-02 6.000000e-02 3.560000e+01 3.600000e+00 3.600000e+00
+2.870000e+00 1.100000e-01 1.100000e-01 2.510000e+01 2.000000e+00 2.000000e+00
+3.035000e+00 5.500000e-02 5.500000e-02 2.320000e+01 2.500000e+00 2.500000e+00
+END YODA_SCATTER2D_V2
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1979_I141722.cc b/analyses/pluginFrascati/GAMMAGAMMA_1979_I141722.cc
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1979_I141722.cc
@@ -0,0 +1,140 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/Projections/FinalState.hh"
+
+namespace Rivet {
+
+
+ /// @brief Add a short analysis description here
+ class GAMMAGAMMA_1979_I141722 : public Analysis {
+ public:
+
+ /// Constructor
+ DEFAULT_RIVET_ANALYSIS_CTOR(GAMMAGAMMA_1979_I141722);
+
+
+ /// @name Analysis methods
+ //@{
+
+ /// Book histograms and initialise projections before the run
+ void init() {
+ // Initialise and register projections
+ declare(FinalState(), "FS");
+
+ // Book histograms
+ _c_hadrons = bookCounter("/TMP/sigma_hadrons");
+ _c_muons = bookCounter("/TMP/sigma_muons");
+ _c_charged = bookCounter("/TMP/Ncharged");
+ _c_neutral = bookCounter("/TMP/Nneutral");
+ _nHadrons = 0.;
+ }
+
+
+ /// Perform the per-event analysis
+ void analyze(const Event& event) {
+ const FinalState& fs = apply<FinalState>(event, "FS");
+
+ map<long,int> nCount;
+ int ntotal(0),ncharged(0),nneutral(0);
+ foreach (const Particle& p, fs.particles()) {
+ nCount[p.pdgId()] += 1;
+ ++ntotal;
+ if(PID::isCharged(p.pdgId()))
+ ncharged += 1;
+ else
+ nneutral += 1;
+ }
+ // mu+mu- + photons
+ if(nCount[-13]==1 and nCount[13]==1 &&
+ ntotal==2+nCount[22])
+ _c_muons->fill(event.weight());
+ // everything else
+ else {
+ if(ntotal==2) vetoEvent;
+ _c_hadrons->fill(event.weight());
+ _c_charged->fill(ncharged*event.weight());
+ _c_neutral->fill(nneutral*event.weight());
+ _nHadrons += event.weight();
+ }
+ }
+
+
+ /// Normalise histograms etc., after the run
+ void finalize() {
+ Scatter1D R = *_c_hadrons/ *_c_muons;
+ double rval = R.point(0).x();
+ pair<double,double> rerr = R.point(0).xErrs();
+ double fact = crossSection()/ sumOfWeights() /picobarn;
+ double sig_h = _c_hadrons->val()*fact;
+ double err_h = _c_hadrons->err()*fact;
+ double sig_m = _c_muons ->val()*fact;
+ double err_m = _c_muons ->err()*fact;
+ Scatter2D temphisto(refData(1, 1, 1));
+ Scatter2DPtr hadrons = bookScatter2D("sigma_hadrons");
+ Scatter2DPtr muons = bookScatter2D("sigma_muons" );
+ Scatter2DPtr mult = bookScatter2D(1, 1, 1);
+ for (size_t b = 0; b < temphisto.numPoints(); b++) {
+ const double x = temphisto.point(b).x();
+ pair<double,double> ex = temphisto.point(b).xErrs();
+ pair<double,double> ex2 = ex;
+ if(ex2.first ==0.) ex2. first=0.0001;
+ if(ex2.second==0.) ex2.second=0.0001;
+ if (inRange(sqrtS()/GeV, x-ex2.first, x+ex2.second)) {
+ mult ->addPoint(x, rval, ex, rerr);
+ hadrons->addPoint(x, sig_h, ex, make_pair(err_h,err_h));
+ muons ->addPoint(x, sig_m, ex, make_pair(err_m,err_m));
+ }
+ else {
+ mult ->addPoint(x, 0., ex, make_pair(0.,.0));
+ hadrons->addPoint(x, 0., ex, make_pair(0.,.0));
+ muons ->addPoint(x, 0., ex, make_pair(0.,.0));
+ }
+ }
+ scale(_c_charged, 1./_nHadrons);
+ scale(_c_neutral, 1./_nHadrons);
+ for(unsigned int iy=1; iy<3;++iy) {
+ double aver(0.),error(0.);
+ if(iy==1) {
+ aver = _c_charged->val();
+ error = _c_charged->err();
+ }
+ else {
+ aver = _c_neutral->val();
+ error = _c_neutral->err();
+ }
+ Scatter2D temphisto(refData(2, 1, iy));
+ Scatter2DPtr mult = bookScatter2D(2, 1, iy);
+ for (size_t b = 0; b < temphisto.numPoints(); b++) {
+ const double x = temphisto.point(b).x();
+ pair<double,double> ex = temphisto.point(b).xErrs();
+ pair<double,double> ex2 = ex;
+ if(ex2.first ==0.) ex2. first=0.0001;
+ if(ex2.second==0.) ex2.second=0.0001;
+ if (inRange(sqrtS()/GeV, x-ex2.first, x+ex2.second)) {
+ mult ->addPoint(x, aver, ex, make_pair(error,error));
+ }
+ else {
+ mult ->addPoint(x, 0., ex, make_pair(0.,.0));
+ }
+ }
+ }
+ }
+
+ //@}
+
+
+ /// @name Histograms
+ //@{
+ CounterPtr _c_hadrons, _c_muons,_c_neutral,_c_charged;
+ double _nHadrons;
+ //@}
+
+
+ };
+
+
+ // The hook for the plugin system
+ DECLARE_RIVET_PLUGIN(GAMMAGAMMA_1979_I141722);
+
+
+}
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1979_I141722.info b/analyses/pluginFrascati/GAMMAGAMMA_1979_I141722.info
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1979_I141722.info
@@ -0,0 +1,33 @@
+Name: GAMMAGAMMA_1979_I141722
+Year: 1979
+Summary: Measurement of $R$ and the hadron multiplicity between 1.42 and 3.09 GeV
+Experiment: GAMMAGAMMA
+Collider: ADONE
+InspireID: 141722
+Status: VALIDATED
+Authors:
+ - Peter Richardson <peter.richardson@durham.ac.uk>
+References:
+ - Phys.Lett. B86 (1979) 234-238, 1979
+RunInfo: e+ e- to hadrons and e+ e- to mu+ mu- (for normalization)
+NeedCrossSection: no
+Beams: [e-, e+]
+Description:
+ 'Measurement of $R$ in $e^+e^-$ collisions by Gamma-Gamma-2 for energies between 1.42 and 3.09 GeV.
+ The average charged and neutral particle multiplicity is also measured.
+ The individual hadronic and muonic cross sections are also outputted to the yoda file
+ so that ratio $R$ can be recalcuated if runs are combined.'
+BibKey: Bacci:1979ab
+BibTeX: '@article{Bacci:1979ab,
+ author = "Bacci, C. and others",
+ title = "{Total Cross-Section for Hadronic Production by e+ e-
+ Annihilation in the Total Center-Of-Mass Energy Range
+ 1.42-GeV - 3.09-GeV}",
+ journal = "Phys. Lett.",
+ volume = "86B",
+ year = "1979",
+ pages = "234-238",
+ doi = "10.1016/0370-2693(79)90828-1",
+ reportNumber = "LNF-79/34-P",
+ SLACcitation = "%%CITATION = PHLTA,86B,234;%%"
+}'
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1979_I141722.plot b/analyses/pluginFrascati/GAMMAGAMMA_1979_I141722.plot
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1979_I141722.plot
@@ -0,0 +1,35 @@
+BEGIN PLOT /GAMMAGAMMA_1979_I141722/d01-x01-y01
+Title=$R=\sigma(e^+e^-\to \text{hadrons})/\sigma(e^+e^-\to \mu^+\mu^-)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$R$
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /GAMMAGAMMA_1979_I141722/sigma_hadrons
+Title=$\sigma(e^+e^-\to \text{hadrons})$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to \text{hadrons})/pb$
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /GAMMAGAMMA_1979_I141722/sigma_muons
+Title=$\sigma(e^+e^-\to \mu^+\mu^-)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to \mu^+\mu^-)/pb$
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /GAMMAGAMMA_1979_I141722/d02-x01-y01
+Title=Average Charged particle Multiplicity
+XLabel=$\sqrt{s}$/GeV
+YLabel=$N_{\text{charged}}$
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /GAMMAGAMMA_1979_I141722/d02-x01-y02
+Title=Average Neutral particle Multiplicity
+XLabel=$\sqrt{s}$/GeV
+YLabel=$N_{\text{neutral}}$
+LogY=0
+ConnectGaps=1
+END PLOT
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1979_I141722.yoda b/analyses/pluginFrascati/GAMMAGAMMA_1979_I141722.yoda
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1979_I141722.yoda
@@ -0,0 +1,91 @@
+BEGIN YODA_SCATTER2D_V2 /REF/GAMMAGAMMA_1979_I141722/d01-x01-y01
+IsRef: 1
+Path: /REF/GAMMAGAMMA_1979_I141722/d01-x01-y01
+Title: doi:10.17182/hepdata.27311.v1/t1
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.435000e+00 1.600000e-02 1.500000e-02 1.460000e+00 3.456776e-01 3.456776e-01
+1.464000e+00 1.400000e-02 1.100000e-02 1.580000e+00 3.721519e-01 3.721519e-01
+1.491000e+00 1.600000e-02 9.000000e-03 1.930000e+00 4.102539e-01 4.102539e-01
+1.506000e+00 6.000000e-03 1.800000e-02 2.040000e+00 4.322823e-01 4.322823e-01
+1.525000e+00 0.000000e+00 0.000000e+00 2.380000e+00 5.750974e-01 5.750974e-01
+1.539000e+00 1.300000e-02 1.000000e-02 2.280000e+00 5.576128e-01 5.576128e-01
+1.550000e+00 0.000000e+00 0.000000e+00 2.200000e+00 5.550676e-01 5.550676e-01
+1.560000e+00 9.000000e-03 1.400000e-02 2.200000e+00 6.181424e-01 6.181424e-01
+1.575000e+00 0.000000e+00 0.000000e+00 2.090000e+00 5.138524e-01 5.138524e-01
+1.587000e+00 1.100000e-02 1.300000e-02 1.860000e+00 4.994327e-01 4.994327e-01
+1.610000e+00 1.000000e-02 1.500000e-02 2.070000e+00 4.945789e-01 4.945789e-01
+1.638000e+00 1.300000e-02 1.200000e-02 1.780000e+00 4.330785e-01 4.330785e-01
+1.663000e+00 1.300000e-02 1.200000e-02 2.350000e+00 5.698300e-01 5.698300e-01
+1.687000e+00 1.200000e-02 1.300000e-02 2.100000e+00 5.267115e-01 5.267115e-01
+1.713000e+00 1.300000e-02 1.200000e-02 2.100000e+00 5.210806e-01 5.210806e-01
+1.738000e+00 1.300000e-02 1.200000e-02 1.940000e+00 5.060168e-01 5.060168e-01
+1.764000e+00 1.400000e-02 1.100000e-02 1.680000e+00 4.271440e-01 4.271440e-01
+1.790000e+00 1.500000e-02 1.000000e-02 1.820000e+00 4.803925e-01 4.803925e-01
+1.812000e+00 1.200000e-02 1.300000e-02 1.700000e+00 4.035158e-01 4.035158e-01
+1.834000e+00 9.000000e-03 1.600000e-02 1.910000e+00 4.611337e-01 4.611337e-01
+1.859000e+00 9.000000e-03 1.600000e-02 1.550000e+00 4.382993e-01 4.382993e-01
+1.886000e+00 1.100000e-02 1.400000e-02 1.740000e+00 4.874146e-01 4.874146e-01
+1.930000e+00 3.000000e-02 4.200000e-02 1.580000e+00 3.987443e-01 3.987443e-01
+1.960000e+00 4.000000e-02 4.000000e-02 1.610000e+00 3.684620e-01 3.684620e-01
+2.050000e+00 5.000000e-02 5.000000e-02 1.760000e+00 4.050284e-01 4.050284e-01
+2.150000e+00 5.000000e-02 5.000000e-02 2.000000e+00 4.294182e-01 4.294182e-01
+2.370000e+00 1.700000e-01 1.700000e-01 2.610000e+00 5.975904e-01 5.975904e-01
+2.590000e+00 5.000000e-02 5.000000e-02 2.770000e+00 6.360018e-01 6.360018e-01
+2.700000e+00 6.000000e-02 6.000000e-02 2.850000e+00 6.552147e-01 6.552147e-01
+2.870000e+00 1.100000e-01 1.100000e-01 2.260000e+00 4.994727e-01 4.994727e-01
+3.035000e+00 5.500000e-02 5.500000e-02 2.320000e+00 5.396777e-01 5.396777e-01
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/GAMMAGAMMA_1979_I141722/d02-x01-y01
+IsRef: 1
+Path: /REF/GAMMAGAMMA_1979_I141722/d02-x01-y01
+Title: doi:10.17182/hepdata.27311.v1/t2
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.435000e+00 1.600000e-02 1.500000e-02 3.050000e+00 1.300000e-01 1.300000e-01
+1.485000e+00 3.500000e-02 1.500000e-02 3.040000e+00 5.000000e-02 5.000000e-02
+1.511000e+00 1.100000e-02 3.900000e-02 3.150000e+00 5.000000e-02 5.000000e-02
+1.570000e+00 2.000000e-02 3.000000e-02 3.120000e+00 9.000000e-02 9.000000e-02
+1.624000e+00 2.400000e-02 2.600000e-02 3.250000e+00 1.000000e-01 1.000000e-01
+1.675000e+00 2.500000e-02 2.500000e-02 3.070000e+00 1.100000e-01 1.100000e-01
+1.725000e+00 2.500000e-02 2.500000e-02 3.140000e+00 1.200000e-01 1.200000e-01
+1.776000e+00 2.600000e-02 2.400000e-02 3.000000e+00 1.300000e-01 1.300000e-01
+1.822000e+00 2.200000e-02 2.800000e-02 3.240000e+00 1.000000e-01 1.000000e-01
+1.872000e+00 2.200000e-02 2.800000e-02 3.180000e+00 1.500000e-01 1.500000e-01
+1.930000e+00 3.000000e-02 4.200000e-02 3.090000e+00 1.700000e-01 1.700000e-01
+1.960000e+00 4.000000e-02 4.000000e-02 3.200000e+00 1.100000e-01 1.100000e-01
+2.050000e+00 5.000000e-02 5.000000e-02 3.350000e+00 1.100000e-01 1.100000e-01
+2.150000e+00 5.000000e-02 5.000000e-02 3.530000e+00 7.000000e-02 7.000000e-02
+2.370000e+00 1.700000e-01 1.700000e-01 3.380000e+00 1.100000e-01 1.100000e-01
+2.590000e+00 5.000000e-02 5.000000e-02 3.400000e+00 1.200000e-01 1.200000e-01
+2.700000e+00 6.000000e-02 5.000000e-02 3.630000e+00 1.200000e-01 1.200000e-01
+2.870000e+00 1.100000e-01 1.100000e-01 3.600000e+00 1.100000e-01 1.100000e-01
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/GAMMAGAMMA_1979_I141722/d02-x01-y02
+IsRef: 1
+Path: /REF/GAMMAGAMMA_1979_I141722/d02-x01-y02
+Title: doi:10.17182/hepdata.27311.v1/t2
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.435000e+00 1.600000e-02 1.500000e-02 1.340000e+00 1.300000e-01 1.300000e-01
+1.485000e+00 3.500000e-02 1.500000e-02 1.170000e+00 5.000000e-02 5.000000e-02
+1.511000e+00 1.100000e-02 3.900000e-02 1.310000e+00 5.000000e-02 5.000000e-02
+1.570000e+00 2.000000e-02 3.000000e-02 1.340000e+00 9.000000e-02 9.000000e-02
+1.624000e+00 2.400000e-02 2.600000e-02 1.280000e+00 1.000000e-01 1.000000e-01
+1.675000e+00 2.500000e-02 2.500000e-02 1.340000e+00 1.100000e-01 1.100000e-01
+1.725000e+00 2.500000e-02 2.500000e-02 1.620000e+00 1.200000e-01 1.200000e-01
+1.776000e+00 2.600000e-02 2.400000e-02 1.440000e+00 1.300000e-01 1.300000e-01
+1.822000e+00 2.200000e-02 2.800000e-02 1.610000e+00 1.000000e-01 1.000000e-01
+1.872000e+00 2.200000e-02 2.800000e-02 1.420000e+00 1.500000e-01 1.500000e-01
+1.930000e+00 3.000000e-02 4.200000e-02 1.730000e+00 1.700000e-01 1.700000e-01
+1.960000e+00 4.000000e-02 4.000000e-02 1.750000e+00 1.500000e-01 1.500000e-01
+2.050000e+00 5.000000e-02 5.000000e-02 1.850000e+00 1.500000e-01 1.500000e-01
+2.150000e+00 5.000000e-02 5.000000e-02 1.900000e+00 1.000000e-01 1.000000e-01
+2.370000e+00 1.700000e-01 1.700000e-01 1.900000e+00 1.500000e-01 1.500000e-01
+2.590000e+00 5.000000e-02 5.000000e-02 2.350000e+00 1.600000e-01 1.600000e-01
+2.700000e+00 6.000000e-02 5.000000e-02 2.400000e+00 1.600000e-01 1.600000e-01
+2.870000e+00 1.100000e-01 1.100000e-01 2.250000e+00 1.500000e-01 1.500000e-01
+END YODA_SCATTER2D_V2
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1980_I153382.cc b/analyses/pluginFrascati/GAMMAGAMMA_1980_I153382.cc
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1980_I153382.cc
@@ -0,0 +1,78 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/Projections/FinalState.hh"
+
+namespace Rivet {
+
+
+ /// @brief Add a short analysis description here
+ class GAMMAGAMMA_1980_I153382 : public Analysis {
+ public:
+
+ /// Constructor
+ DEFAULT_RIVET_ANALYSIS_CTOR(GAMMAGAMMA_1980_I153382);
+
+
+ /// @name Analysis methods
+ //@{
+
+ /// Book histograms and initialise projections before the run
+ void init() {
+ // Initialise and register projections
+ declare(FinalState(), "FS");
+ // Book histograms
+ _npion = bookCounter("TMP/pion");
+ }
+
+
+ /// Perform the per-event analysis
+ void analyze(const Event& event) {
+ const FinalState& fs = apply<FinalState>(event, "FS");
+ if(fs.particles().size()!=4) vetoEvent;
+ foreach (const Particle& p, fs.particles()) {
+ if(abs(p.pdgId())!=PID::PIPLUS) vetoEvent;
+ }
+ _npion->fill(event.weight());
+ }
+
+
+ /// Normalise histograms etc., after the run
+ void finalize() {
+ double sigma = _npion->val();
+ double error = _npion->err();
+ sigma *= crossSection()/ sumOfWeights() /nanobarn;
+ error *= crossSection()/ sumOfWeights() /nanobarn;
+ Scatter2D temphisto(refData(1, 1, 1));
+ Scatter2DPtr mult = bookScatter2D(1, 1, 1);
+ for (size_t b = 0; b < temphisto.numPoints(); b++) {
+ const double x = temphisto.point(b).x();
+ pair<double,double> ex = temphisto.point(b).xErrs();
+ pair<double,double> ex2 = ex;
+ if(ex2.first ==0.) ex2. first=0.0001;
+ if(ex2.second==0.) ex2.second=0.0001;
+ if (inRange(sqrtS()/GeV, x-ex2.first, x+ex2.second)) {
+ mult->addPoint(x, sigma, ex, make_pair(error,error));
+ }
+ else {
+ mult->addPoint(x, 0., ex, make_pair(0.,.0));
+ }
+ }
+ }
+
+ //@}
+
+
+ /// @name Histograms
+ //@{
+ CounterPtr _npion;
+ //@}
+
+
+ };
+
+
+ // The hook for the plugin system
+ DECLARE_RIVET_PLUGIN(GAMMAGAMMA_1980_I153382);
+
+
+}
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1980_I153382.info b/analyses/pluginFrascati/GAMMAGAMMA_1980_I153382.info
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1980_I153382.info
@@ -0,0 +1,31 @@
+Name: GAMMAGAMMA_1980_I153382
+Year: 1980
+Summary: Cross section for $e^+e^-\to\pi^+\pi^+\pi^-\pi^+$ between 1.456 and 2.15 GeV
+Experiment: GAMMAGAMMA
+Collider: ADONE
+InspireID: 153382
+Status: VALIDATED
+Authors:
+ - Peter Richardson <peter.richardson@durham.ac.uk>
+References:
+ - Phys.Lett. B95 (1980) 139-142, 1980
+RunInfo: e+e- to hadrons
+NeedCrossSection: yes
+Beams: [e+,e-]
+Luminosity_fb:
+Description:
+ 'Measurement of the cross section for $e^+e^-\to\pi^+\pi^+\pi^-\pi^+$ for energies between 1.456 and 2.15 GeV'
+Keywords:
+BibKey: Bacci:1980ru
+BibTeX: '@article{Bacci:1980ru,
+ author = "Bacci, C. and others",
+ title = "{Measurement of the $e^+ e^- \to \pi^+ \pi^- \pi^+ \pi^-$
+ Cross-section in the $\rho^\prime$ (1600) Energy Region}",
+ journal = "Phys. Lett.",
+ volume = "95B",
+ year = "1980",
+ pages = "139-142",
+ doi = "10.1016/0370-2693(80)90418-9",
+ reportNumber = "LNF-80/25-P",
+ SLACcitation = "%%CITATION = PHLTA,95B,139;%%"
+}'
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1980_I153382.plot b/analyses/pluginFrascati/GAMMAGAMMA_1980_I153382.plot
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1980_I153382.plot
@@ -0,0 +1,8 @@
+BEGIN PLOT /GAMMAGAMMA_1980_I153382/d01-x01-y01
+Title=$\sigma(e^+e^-\to 2\pi^+2\pi^-)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to 2\pi^+2\pi^-)$/nb
+LogY=0
+LegendXPos=0.8
+ConnectGaps=1
+END PLOT
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1980_I153382.yoda b/analyses/pluginFrascati/GAMMAGAMMA_1980_I153382.yoda
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1980_I153382.yoda
@@ -0,0 +1,23 @@
+BEGIN YODA_SCATTER2D_V2 /REF/GAMMAGAMMA_1980_I153382/d01-x01-y01
+IsRef: 1
+Path: /REF/GAMMAGAMMA_1980_I153382/d01-x01-y01
+Title: doi:10.17182/hepdata.27154.v1/t1
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.456000e+00 3.700000e-02 1.900000e-02 2.260000e+01 5.319972e+00 5.319972e+00
+1.491000e+00 1.600000e-02 9.000000e-03 3.060000e+01 5.134014e+00 5.134014e+00
+1.507000e+00 7.000000e-03 1.800000e-02 2.340000e+01 4.490000e+00 4.490000e+00
+1.546000e+00 2.100000e-02 2.900000e-02 2.630000e+01 5.273806e+00 5.273806e+00
+1.600000e+00 2.500000e-02 1.500000e-02 2.440000e+01 5.064148e+00 5.064148e+00
+1.633000e+00 1.800000e-02 1.700000e-02 2.240000e+01 4.438423e+00 4.438423e+00
+1.660000e+00 1.000000e-02 1.000000e-02 2.930000e+01 6.010493e+00 6.010493e+00
+1.701000e+00 3.100000e-02 2.400000e-02 1.730000e+01 3.402062e+00 3.402062e+00
+1.752000e+00 2.700000e-02 2.300000e-02 1.680000e+01 3.480000e+00 3.480000e+00
+1.794000e+00 1.900000e-02 1.100000e-02 1.320000e+01 3.189106e+00 3.189106e+00
+1.821000e+00 1.600000e-02 1.400000e-02 1.010000e+01 2.589445e+00 2.589445e+00
+1.851000e+00 1.600000e-02 2.400000e-02 1.190000e+01 2.911396e+00 2.911396e+00
+1.935000e+00 6.000000e-02 6.500000e-02 8.200000e+00 1.863572e+00 1.863572e+00
+2.050000e+00 5.000000e-02 5.000000e-02 8.300000e+00 2.441316e+00 2.441316e+00
+2.150000e+00 5.000000e-02 5.000000e-02 5.300000e+00 1.523819e+00 1.523819e+00
+END YODA_SCATTER2D_V2
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1981_I158474.cc b/analyses/pluginFrascati/GAMMAGAMMA_1981_I158474.cc
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1981_I158474.cc
@@ -0,0 +1,177 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/Projections/FinalState.hh"
+
+namespace Rivet {
+
+
+ /// @brief Add a short analysis description here
+ class GAMMAGAMMA_1981_I158474 : public Analysis {
+ public:
+
+ /// Constructor
+ DEFAULT_RIVET_ANALYSIS_CTOR(GAMMAGAMMA_1981_I158474);
+
+
+ /// @name Analysis methods
+ //@{
+
+ /// Book histograms and initialise projections before the run
+ void init() {
+ // Initialise and register projections
+ declare(FinalState(), "FS");
+ _n3pi = bookCounter("TMP/n3pi");
+ _n4pi = bookCounter("TMP/n4pi");
+ _n5pi = bookCounter("TMP/n5pi");
+ _n6pi = bookCounter("TMP/n6pi");
+ _n35pi = bookCounter("TMP/n35pi");
+ _n46pi = bookCounter("TMP/n46pi");
+ _nC2 = bookCounter("TMP/nC2");
+ _nC4 = bookCounter("TMP/nC4");
+ _nmu = bookCounter("TMP/nmu");
+ }
+
+
+ /// Perform the per-event analysis
+ void analyze(const Event& event) {
+ const FinalState& fs = apply<FinalState>(event, "FS");
+
+ map<long,int> nCount;
+ int ntotal(0);
+ foreach (const Particle& p, fs.particles()) {
+ nCount[p.pdgId()] += 1;
+ ++ntotal;
+ }
+ // mu+mu- + photons
+ if(nCount[-13]==1 and nCount[13]==1 &&
+ ntotal==2+nCount[22])
+ _nmu->fill(event.weight());
+ else {
+ if(ntotal==3 && nCount[211] == 1 && nCount[-211]==1 && nCount[111]==1 ) {
+ _n3pi->fill(event.weight());
+ }
+ if(ntotal==4 && nCount[211] == 1 && nCount[-211]==1 && nCount[111]==2 ) {
+ _n4pi->fill(event.weight());
+ }
+ if(ntotal==5 && nCount[211] == 2 && nCount[-211]==2 && nCount[111]==1 ) {
+ _n5pi->fill(event.weight());
+ }
+ if(ntotal==6 && nCount[211] == 2 && nCount[-211]==2 && nCount[111]==2 ) {
+ _n6pi->fill(event.weight());
+ }
+ if(nCount[211] == 1 && nCount[-211]==1 && ntotal == 2+nCount[111]) {
+ _nC2->fill(event.weight());
+ }
+ if(nCount[211] == 2 && nCount[-211]==2 && ntotal == 4+nCount[111]) {
+ _nC4->fill(event.weight());
+ }
+ if((nCount[211]+nCount[-211]+nCount[111])==ntotal ) {
+ if(ntotal==3 || ntotal ==5)
+ _n35pi->fill(event.weight());
+ else if(ntotal==4 || ntotal==6)
+ _n46pi ->fill(event.weight());
+ }
+ }
+ }
+
+
+ /// Normalise histograms etc., after the run
+ void finalize() {
+ double fact = crossSection()/ sumOfWeights() /nanobarn;
+ for(unsigned int ix=1;ix<7;++ix) {
+ double sigma,error;
+ if(ix==1) {
+ sigma = _n3pi->val()*fact;
+ error = _n3pi->err()*fact;
+ }
+ else if(ix==2) {
+ sigma = _n4pi->val()*fact;
+ error = _n4pi->err()*fact;
+ }
+ else if(ix==3) {
+ sigma = _n5pi->val()*fact;
+ error = _n5pi->err()*fact;
+ }
+ else if(ix==4) {
+ sigma = _n6pi->val()*fact;
+ error = _n6pi->err()*fact;
+ }
+ else if(ix==5) {
+ sigma = _n35pi->val()*fact;
+ error = _n35pi->err()*fact;
+ }
+ else if(ix==6) {
+ sigma = _n46pi->val()*fact;
+ error = _n46pi->err()*fact;
+ }
+ Scatter2D temphisto(refData(1, 1, ix));
+ Scatter2DPtr mult = bookScatter2D(1, 1, ix);
+ for (size_t b = 0; b < temphisto.numPoints(); b++) {
+ const double x = temphisto.point(b).x();
+ pair<double,double> ex = temphisto.point(b).xErrs();
+ pair<double,double> ex2 = ex;
+ if(ex2.first ==0.) ex2. first=0.0001;
+ if(ex2.second==0.) ex2.second=0.0001;
+ if (inRange(sqrtS()/GeV, x-ex2.first, x+ex2.second)) {
+ mult->addPoint(x, sigma, ex, make_pair(error,error));
+ }
+ else {
+ mult->addPoint(x, 0., ex, make_pair(0.,.0));
+ }
+ }
+ }
+ for(unsigned int ix=1;ix<3;++ix) {
+ Scatter1D R = (ix==1? *_nC2 : *_nC4)/ *_nmu;
+ double rval = R.point(0).x();
+ pair<double,double> rerr = R.point(0).xErrs();
+ double sig_h = (ix ==1 ? _nC2 : _nC4)->val()*fact;
+ double err_h = (ix ==1 ? _nC2 : _nC4)->err()*fact;
+ double sig_m = _nmu->val()*fact;
+ double err_m = _nmu->err()*fact;
+ Scatter2D temphisto(refData(2, 1, ix));
+ ostringstream title;
+ if(ix==1)
+ title << "sigma_2pi";
+ else
+ title << "sigma_4pi";
+ Scatter2DPtr hadrons = bookScatter2D(title.str());
+ Scatter2DPtr muons;
+ if(ix==1) muons = bookScatter2D("sigma_muons");
+ Scatter2DPtr mult = bookScatter2D(2,1,ix);
+ for (size_t b = 0; b < temphisto.numPoints(); b++) {
+ const double x = temphisto.point(b).x();
+ pair<double,double> ex = temphisto.point(b).xErrs();
+ pair<double,double> ex2 = ex;
+ if(ex2.first ==0.) ex2. first=0.0001;
+ if(ex2.second==0.) ex2.second=0.0001;
+ if (inRange(sqrtS()/GeV, x-ex2.first, x+ex2.second)) {
+ mult ->addPoint(x, rval, ex, rerr);
+ hadrons->addPoint(x, sig_h, ex, make_pair(err_h,err_h));
+ if(ix==1) muons ->addPoint(x, sig_m, ex, make_pair(err_m,err_m));
+ }
+ else {
+ mult ->addPoint(x, 0., ex, make_pair(0.,.0));
+ hadrons->addPoint(x, 0., ex, make_pair(0.,.0));
+ if(ix==1) muons ->addPoint(x, 0., ex, make_pair(0.,.0));
+ }
+ }
+ }
+ }
+
+ //@}
+
+
+ /// @name Histograms
+ //@{
+ CounterPtr _n3pi,_n4pi,_n5pi,_n6pi,_n35pi,_n46pi,_nC2,_nC4,_nmu;
+ //@}
+
+
+ };
+
+
+ // The hook for the plugin system
+ DECLARE_RIVET_PLUGIN(GAMMAGAMMA_1981_I158474);
+
+
+}
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1981_I158474.info b/analyses/pluginFrascati/GAMMAGAMMA_1981_I158474.info
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1981_I158474.info
@@ -0,0 +1,31 @@
+Name: GAMMAGAMMA_1981_I158474
+Year: 1981
+Summary: Multi pion cross sections for energies between 1.44 and 2.15 GeV
+Experiment: GAMMAGAMMA
+Collider: ADONE
+InspireID: 158474
+Status: VALIDATED
+Authors:
+ - Peter Richardson <peter.richardson@durham.ac.uk>
+References:
+ - Nucl.Phys. B184 (1981) 31-39, 1981
+RunInfo: e+e- to hadrons
+NeedCrossSection: yes
+Beams: [e+,e-]
+Luminosity_fb:
+Description:
+ 'Measurement of the cross section for $e^+e^-\to$pions for energies between 1.44 and 2.15 GeV'
+Keywords:
+BibKey: Bacci:1980zs
+BibTeX: '@article{Bacci:1980zs,
+ author = "Bacci, C. and others",
+ title = "{Measurement of Hadronic Exclusive Cross-sections in $e^+
+ e^-$ Annihilation From 1.42-{GeV} to 2.20-{GeV}}",
+ journal = "Nucl. Phys.",
+ volume = "B184",
+ year = "1981",
+ pages = "31-39",
+ doi = "10.1016/0550-3213(81)90208-X",
+ reportNumber = "LNF-80/72-P",
+ SLACcitation = "%%CITATION = NUPHA,B184,31;%%"
+}'
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1981_I158474.plot b/analyses/pluginFrascati/GAMMAGAMMA_1981_I158474.plot
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1981_I158474.plot
@@ -0,0 +1,77 @@
+BEGIN PLOT /GAMMAGAMMA_1981_I158474/d01-x01-y01
+Title=$\sigma(e^+e^-\to \pi^+\pi^-\pi^0)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to \pi^+\pi^-\pi^0)$/nb
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /GAMMAGAMMA_1981_I158474/d01-x01-y02
+Title=$\sigma(e^+e^-\to \pi^+\pi^-2\pi^0)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to \pi^+\pi^-2\pi^0)$/nb
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /GAMMAGAMMA_1981_I158474/d01-x01-y03
+Title=$\sigma(e^+e^-\to 2\pi^+2\pi^-\pi^0)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to 2\pi^+2\pi^-\pi^0)$/nb
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /GAMMAGAMMA_1981_I158474/d01-x01-y04
+Title=$\sigma(e^+e^-\to 2\pi^+2\pi^-2\pi^0)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to 2\pi^+2\pi^-2\pi^0)$/nb
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /GAMMAGAMMA_1981_I158474/d01-x01-y05
+Title=$\sigma(e^+e^-\to 3\pi+5\pi)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to 3\pi+5\pi)$/nb
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /GAMMAGAMMA_1981_I158474/d01-x01-y06
+Title=$\sigma(e^+e^-\to 4\pi+6\pi)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to 4\pi^+6\pi^-)$/nb
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /GAMMAGAMMA_1981_I158474/d02-x01-y01
+Title==$R=\sigma(e^+e^-\to \pi^+\pi^-+n\pi^0)/\sigma(e^+e^-\to \mu^+\mu^-)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$R$
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /GAMMAGAMMA_1981_I158474/d02-x01-y02
+Title==$R=\sigma(e^+e^-\to 2\pi^+2\pi^-+n\pi^0)/\sigma(e^+e^-\to \mu^+\mu^-)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$R$
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /GAMMAGAMMA_1981_I158474/sigma_2pi
+Title=$\sigma(e^+e^-\to \pi^+\pi^-+n\pi^0)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to \pi^+\pi^-+n\pi^0)/nb$
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /GAMMAGAMMA_1981_I158474/sigma_4pi
+Title=$\sigma(e^+e^-\to 2\pi^+2\pi^-+n\pi^0)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to 2\pi^+2\pi^-+n\pi^0)/nb$
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /GAMMAGAMMA_1981_I158474/sigma_muons
+Title=$\sigma(e^+e^-\to \mu^+\mu^-)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to \mu^+\mu^-)/nb$
+LogY=0
+ConnectGaps=1
+END PLOT
\ No newline at end of file
diff --git a/analyses/pluginFrascati/GAMMAGAMMA_1981_I158474.yoda b/analyses/pluginFrascati/GAMMAGAMMA_1981_I158474.yoda
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/GAMMAGAMMA_1981_I158474.yoda
@@ -0,0 +1,192 @@
+BEGIN YODA_SCATTER2D_V2 /REF/GAMMAGAMMA_1981_I158474/d01-x01-y01
+IsRef: 1
+Path: /REF/GAMMAGAMMA_1981_I158474/d01-x01-y01
+Title: doi:10.17182/hepdata.34278.v1/t1
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.437000e+00 1.700000e-02 3.800000e-02 5.200000e+00 3.400000e+00 3.400000e+00
+1.491000e+00 1.600000e-02 9.000000e-03 0.000000e+00 1.200000e+00 1.200000e+00
+1.505000e+00 5.000000e-03 2.000000e-02 4.000000e-01 1.100000e+00 1.100000e+00
+1.549000e+00 2.400000e-02 2.600000e-02 1.000000e-01 1.300000e+00 1.300000e+00
+1.600000e+00 2.500000e-02 1.500000e-02 0.000000e+00 2.600000e+00 2.600000e+00
+1.633000e+00 1.800000e-02 1.700000e-02 0.000000e+00 3.600000e+00 3.600000e+00
+1.660000e+00 1.000000e-02 1.000000e-02 1.220000e+01 8.000000e+00 8.000000e+00
+1.701000e+00 3.100000e-02 2.400000e-02 0.000000e+00 1.500000e+00 1.500000e+00
+1.752000e+00 2.700000e-02 2.300000e-02 0.000000e+00 2.400000e+00 2.400000e+00
+1.794000e+00 1.900000e-02 1.100000e-02 0.000000e+00 3.200000e+00 3.200000e+00
+1.821000e+00 1.600000e-02 1.400000e-02 6.000000e-01 4.800000e+00 4.800000e+00
+1.851000e+00 1.600000e-02 2.400000e-02 4.900000e+00 5.000000e+00 5.000000e+00
+1.935000e+00 6.000000e-02 6.500000e-02 0.000000e+00 1.700000e+00 1.700000e+00
+2.050000e+00 5.000000e-02 5.000000e-02 0.000000e+00 3.000000e+00 3.000000e+00
+2.150000e+00 5.000000e-02 5.000000e-02 0.000000e+00 2.000000e+00 2.000000e+00
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/GAMMAGAMMA_1981_I158474/d01-x01-y02
+IsRef: 1
+Path: /REF/GAMMAGAMMA_1981_I158474/d01-x01-y02
+Title: doi:10.17182/hepdata.34278.v1/t1
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.437000e+00 1.700000e-02 3.800000e-02 1.670000e+01 4.700000e+00 4.700000e+00
+1.491000e+00 1.600000e-02 9.000000e-03 3.160000e+01 2.700000e+00 2.700000e+00
+1.505000e+00 5.000000e-03 2.000000e-02 2.800000e+01 3.100000e+00 3.100000e+00
+1.549000e+00 2.400000e-02 2.600000e-02 3.120000e+01 3.400000e+00 3.400000e+00
+1.600000e+00 2.500000e-02 1.500000e-02 1.690000e+01 5.800000e+00 5.800000e+00
+1.633000e+00 1.800000e-02 1.700000e-02 2.000000e+01 4.900000e+00 4.900000e+00
+1.660000e+00 1.000000e-02 1.000000e-02 2.010000e+01 8.100000e+00 8.100000e+00
+1.701000e+00 3.100000e-02 2.400000e-02 2.730000e+01 4.400000e+00 4.400000e+00
+1.752000e+00 2.700000e-02 2.300000e-02 2.210000e+01 4.500000e+00 4.500000e+00
+1.794000e+00 1.900000e-02 1.100000e-02 2.710000e+01 5.400000e+00 5.400000e+00
+1.821000e+00 1.600000e-02 1.400000e-02 1.630000e+01 4.900000e+00 4.900000e+00
+1.851000e+00 1.600000e-02 2.400000e-02 1.040000e+01 5.500000e+00 5.500000e+00
+1.935000e+00 6.000000e-02 6.500000e-02 1.690000e+01 3.200000e+00 3.200000e+00
+2.050000e+00 5.000000e-02 5.000000e-02 1.660000e+01 4.000000e+00 4.000000e+00
+2.150000e+00 5.000000e-02 5.000000e-02 1.500000e+01 2.700000e+00 2.700000e+00
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/GAMMAGAMMA_1981_I158474/d01-x01-y03
+IsRef: 1
+Path: /REF/GAMMAGAMMA_1981_I158474/d01-x01-y03
+Title: doi:10.17182/hepdata.34278.v1/t1
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.437000e+00 1.700000e-02 3.800000e-02 4.400000e+00 2.400000e+00 2.400000e+00
+1.491000e+00 1.600000e-02 9.000000e-03 6.200000e+00 1.500000e+00 1.500000e+00
+1.505000e+00 5.000000e-03 2.000000e-02 6.100000e+00 1.900000e+00 1.900000e+00
+1.549000e+00 2.400000e-02 2.600000e-02 7.100000e+00 2.400000e+00 2.400000e+00
+1.600000e+00 2.500000e-02 1.500000e-02 6.700000e+00 4.900000e+00 4.900000e+00
+1.633000e+00 1.800000e-02 1.700000e-02 3.900000e+00 3.900000e+00 3.900000e+00
+1.660000e+00 1.000000e-02 1.000000e-02 0.000000e+00 5.000000e+00 5.000000e+00
+1.701000e+00 3.100000e-02 2.400000e-02 0.000000e+00 4.000000e+00 4.000000e+00
+1.752000e+00 2.700000e-02 2.300000e-02 3.400000e+00 3.400000e+00 3.400000e+00
+1.794000e+00 1.900000e-02 1.100000e-02 4.100000e+00 2.800000e+00 2.800000e+00
+1.821000e+00 1.600000e-02 1.400000e-02 6.700000e+00 3.400000e+00 3.400000e+00
+1.851000e+00 1.600000e-02 2.400000e-02 1.000000e-01 3.200000e+00 3.200000e+00
+1.935000e+00 6.000000e-02 6.500000e-02 2.500000e+00 2.300000e+00 2.300000e+00
+2.050000e+00 5.000000e-02 5.000000e-02 2.400000e+00 3.200000e+00 3.200000e+00
+2.150000e+00 5.000000e-02 5.000000e-02 2.700000e+00 1.600000e+00 1.600000e+00
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/GAMMAGAMMA_1981_I158474/d01-x01-y04
+IsRef: 1
+Path: /REF/GAMMAGAMMA_1981_I158474/d01-x01-y04
+Title: doi:10.17182/hepdata.34278.v1/t1
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.437000e+00 1.700000e-02 3.800000e-02 6.600000e+00 2.900000e+00 2.900000e+00
+1.491000e+00 1.600000e-02 9.000000e-03 2.900000e+00 1.800000e+00 1.800000e+00
+1.505000e+00 5.000000e-03 2.000000e-02 8.700000e+00 2.100000e+00 2.100000e+00
+1.549000e+00 2.400000e-02 2.600000e-02 1.150000e+01 3.000000e+00 3.000000e+00
+1.600000e+00 2.500000e-02 1.500000e-02 1.210000e+01 6.300000e+00 6.300000e+00
+1.633000e+00 1.800000e-02 1.700000e-02 8.700000e+00 4.600000e+00 4.600000e+00
+1.660000e+00 1.000000e-02 1.000000e-02 1.350000e+01 4.600000e+00 4.600000e+00
+1.701000e+00 3.100000e-02 2.400000e-02 1.380000e+01 2.600000e+00 2.600000e+00
+1.752000e+00 2.700000e-02 2.300000e-02 6.500000e+00 3.200000e+00 3.200000e+00
+1.794000e+00 1.900000e-02 1.100000e-02 2.400000e+00 2.500000e+00 2.500000e+00
+1.821000e+00 1.600000e-02 1.400000e-02 1.260000e+01 3.200000e+00 3.200000e+00
+1.851000e+00 1.600000e-02 2.400000e-02 1.020000e+01 2.900000e+00 2.900000e+00
+1.935000e+00 6.000000e-02 6.500000e-02 1.020000e+01 2.300000e+00 2.300000e+00
+2.050000e+00 5.000000e-02 5.000000e-02 9.700000e+00 2.800000e+00 2.800000e+00
+2.150000e+00 5.000000e-02 5.000000e-02 1.340000e+01 1.700000e+00 1.700000e+00
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/GAMMAGAMMA_1981_I158474/d01-x01-y05
+IsRef: 1
+Path: /REF/GAMMAGAMMA_1981_I158474/d01-x01-y05
+Title: doi:10.17182/hepdata.34278.v1/t1
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.437000e+00 1.700000e-02 3.800000e-02 1.170000e+01 4.900000e+00 4.900000e+00
+1.491000e+00 1.600000e-02 9.000000e-03 9.300000e+00 2.600000e+00 2.600000e+00
+1.505000e+00 5.000000e-03 2.000000e-02 9.300000e+00 3.100000e+00 3.100000e+00
+1.549000e+00 2.400000e-02 2.600000e-02 1.030000e+01 3.800000e+00 3.800000e+00
+1.600000e+00 2.500000e-02 1.500000e-02 1.000000e+01 7.800000e+00 7.800000e+00
+1.633000e+00 1.800000e-02 1.700000e-02 5.900000e+00 6.900000e+00 6.900000e+00
+1.660000e+00 1.000000e-02 1.000000e-02 8.900000e+00 1.100000e+01 1.100000e+01
+1.701000e+00 3.100000e-02 2.400000e-02 0.000000e+00 6.200000e+00 6.200000e+00
+1.752000e+00 2.700000e-02 2.300000e-02 5.000000e+00 5.600000e+00 5.600000e+00
+1.794000e+00 1.900000e-02 1.100000e-02 6.100000e+00 5.300000e+00 5.300000e+00
+1.821000e+00 1.600000e-02 1.400000e-02 1.060000e+01 7.000000e+00 7.000000e+00
+1.851000e+00 1.600000e-02 2.400000e-02 5.100000e+00 6.900000e+00 6.900000e+00
+1.935000e+00 6.000000e-02 6.500000e-02 3.600000e+00 4.000000e+00 4.000000e+00
+2.050000e+00 5.000000e-02 5.000000e-02 3.600000e+00 5.700000e+00 5.700000e+00
+2.150000e+00 5.000000e-02 5.000000e-02 4.100000e+00 3.100000e+00 3.100000e+00
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/GAMMAGAMMA_1981_I158474/d01-x01-y06
+IsRef: 1
+Path: /REF/GAMMAGAMMA_1981_I158474/d01-x01-y06
+Title: doi:10.17182/hepdata.34278.v1/t1
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.437000e+00 1.700000e-02 3.800000e-02 4.590000e+01 5.600000e+00 5.600000e+00
+1.491000e+00 1.600000e-02 9.000000e-03 6.510000e+01 3.300000e+00 3.300000e+00
+1.505000e+00 5.000000e-03 2.000000e-02 6.070000e+01 3.800000e+00 3.800000e+00
+1.549000e+00 2.400000e-02 2.600000e-02 7.010000e+01 4.700000e+00 4.700000e+00
+1.600000e+00 2.500000e-02 1.500000e-02 5.340000e+01 8.900000e+00 8.900000e+00
+1.633000e+00 1.800000e-02 1.700000e-02 5.110000e+01 7.000000e+00 7.000000e+00
+1.660000e+00 1.000000e-02 1.000000e-02 6.290000e+01 9.500000e+00 9.500000e+00
+1.701000e+00 3.100000e-02 2.400000e-02 5.840000e+01 5.200000e+00 5.200000e+00
+1.752000e+00 2.700000e-02 2.300000e-02 4.540000e+01 5.700000e+00 5.700000e+00
+1.794000e+00 1.900000e-02 1.100000e-02 4.270000e+01 6.000000e+00 6.000000e+00
+1.821000e+00 1.600000e-02 1.400000e-02 3.990000e+01 5.900000e+00 5.900000e+00
+1.851000e+00 1.600000e-02 2.400000e-02 3.250000e+01 6.200000e+00 6.200000e+00
+1.935000e+00 6.000000e-02 6.500000e-02 3.560000e+01 4.000000e+00 4.000000e+00
+2.050000e+00 5.000000e-02 5.000000e-02 3.460000e+01 5.200000e+00 5.200000e+00
+2.150000e+00 5.000000e-02 5.000000e-02 3.510000e+01 3.300000e+00 3.300000e+00
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/GAMMAGAMMA_1981_I158474/d02-x01-y01
+IsRef: 1
+Path: /REF/GAMMAGAMMA_1981_I158474/d02-x01-y01
+Title: doi:10.17182/hepdata.34278.v1/t2
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.437000e+00 1.700000e-02 3.800000e-02 5.800000e-01 1.300000e-01 1.300000e-01
+1.491000e+00 1.600000e-02 9.000000e-03 9.000000e-01 7.000000e-02 7.000000e-02
+1.505000e+00 5.000000e-03 2.000000e-02 8.300000e-01 7.000000e-02 7.000000e-02
+1.549000e+00 2.400000e-02 2.600000e-02 9.700000e-01 1.000000e-01 1.000000e-01
+1.600000e+00 2.500000e-02 1.500000e-02 6.000000e-01 1.700000e-01 1.700000e-01
+1.633000e+00 1.800000e-02 1.700000e-02 6.800000e-01 1.900000e-01 1.900000e-01
+1.660000e+00 1.000000e-02 1.000000e-02 1.030000e+00 2.200000e-01 2.200000e-01
+1.701000e+00 3.100000e-02 2.400000e-02 9.200000e-01 1.400000e-01 1.400000e-01
+1.752000e+00 2.700000e-02 2.300000e-02 8.500000e-01 1.700000e-01 1.700000e-01
+1.794000e+00 1.900000e-02 1.100000e-02 1.100000e+00 2.200000e-01 2.200000e-01
+1.821000e+00 1.600000e-02 1.400000e-02 7.800000e-01 2.400000e-01 2.400000e-01
+1.851000e+00 1.600000e-02 2.400000e-02 6.200000e-01 2.900000e-01 2.900000e-01
+1.935000e+00 6.000000e-02 6.500000e-02 7.900000e-01 1.500000e-01 1.500000e-01
+2.050000e+00 5.000000e-02 5.000000e-02 8.800000e-01 2.200000e-01 2.200000e-01
+2.150000e+00 5.000000e-02 5.000000e-02 8.800000e-01 1.600000e-01 1.600000e-01
+2.370000e+00 1.700000e-01 1.700000e-01 1.230000e+00 3.600000e-01 3.600000e-01
+2.590000e+00 5.000000e-02 5.000000e-02 1.300000e+00 3.800000e-01 3.800000e-01
+2.700000e+00 6.000000e-02 6.000000e-02 9.600000e-01 2.000000e-01 2.000000e-01
+2.870000e+00 1.100000e-01 1.100000e-01 8.300000e-01 1.900000e-01 1.900000e-01
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/GAMMAGAMMA_1981_I158474/d02-x01-y02
+IsRef: 1
+Path: /REF/GAMMAGAMMA_1981_I158474/d02-x01-y02
+Title: doi:10.17182/hepdata.34278.v1/t2
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.437000e+00 1.700000e-02 3.800000e-02 7.400000e-01 6.000000e-02 6.000000e-02
+1.491000e+00 1.600000e-02 9.000000e-03 9.100000e-01 4.000000e-02 4.000000e-02
+1.505000e+00 5.000000e-03 2.000000e-02 9.900000e-01 5.000000e-02 5.000000e-02
+1.549000e+00 2.400000e-02 2.600000e-02 1.130000e+00 7.000000e-02 7.000000e-02
+1.600000e+00 2.500000e-02 1.500000e-02 1.280000e+00 1.400000e-01 1.400000e-01
+1.633000e+00 1.800000e-02 1.700000e-02 1.080000e+00 1.100000e-01 1.100000e-01
+1.660000e+00 1.000000e-02 1.000000e-02 1.370000e+00 1.300000e-01 1.300000e-01
+1.701000e+00 3.100000e-02 2.400000e-02 1.050000e+00 9.000000e-02 9.000000e-02
+1.752000e+00 2.700000e-02 2.300000e-02 9.500000e-01 9.000000e-02 9.000000e-02
+1.794000e+00 1.900000e-02 1.100000e-02 7.400000e-01 9.000000e-02 9.000000e-02
+1.821000e+00 1.600000e-02 1.400000e-02 1.130000e+00 9.000000e-02 9.000000e-02
+1.851000e+00 1.600000e-02 2.400000e-02 8.900000e-01 1.000000e-01 1.000000e-01
+1.935000e+00 6.000000e-02 6.500000e-02 9.500000e-01 7.000000e-02 7.000000e-02
+2.050000e+00 5.000000e-02 5.000000e-02 9.900000e-01 1.100000e-01 1.100000e-01
+2.150000e+00 5.000000e-02 5.000000e-02 1.140000e+00 8.000000e-02 8.000000e-02
+2.370000e+00 1.700000e-01 1.700000e-01 1.160000e+00 1.500000e-01 1.500000e-01
+2.590000e+00 5.000000e-02 5.000000e-02 1.070000e+00 1.400000e-01 1.400000e-01
+2.700000e+00 6.000000e-02 6.000000e-02 1.630000e+00 2.500000e-01 2.500000e-01
+2.870000e+00 1.100000e-01 1.100000e-01 1.090000e+00 1.000000e-01 1.000000e-01
+END YODA_SCATTER2D_V2
diff --git a/analyses/pluginFrascati/KLOE_2005_I655225.cc b/analyses/pluginFrascati/KLOE_2005_I655225.cc
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/KLOE_2005_I655225.cc
@@ -0,0 +1,79 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/Projections/FinalState.hh"
+
+namespace Rivet {
+
+
+ /// @brief Add a short analysis description here
+ class KLOE_2005_I655225 : public Analysis {
+ public:
+
+ /// Constructor
+ DEFAULT_RIVET_ANALYSIS_CTOR(KLOE_2005_I655225);
+
+
+ /// @name Analysis methods
+ //@{
+
+ /// Book histograms and initialise projections before the run
+ void init() {
+ // Initialise and register projections
+ declare(FinalState(), "FS");
+
+ // Book histograms
+ _npion = bookCounter("TMP/pion");
+ }
+
+
+ /// Perform the per-event analysis
+ void analyze(const Event& event) {
+ const FinalState& fs = apply<FinalState>(event, "FS");
+ if(fs.particles().size()!=2) vetoEvent;
+ foreach (const Particle& p, fs.particles()) {
+ if(abs(p.pdgId())!=PID::PIPLUS) vetoEvent;
+ }
+ _npion->fill(event.weight());
+ }
+
+
+ /// Normalise histograms etc., after the run
+ void finalize() {
+ double sigma = _npion->val();
+ double error = _npion->err();
+ sigma *= crossSection()/ sumOfWeights() /nanobarn;
+ error *= crossSection()/ sumOfWeights() /nanobarn;
+ Scatter2D temphisto(refData(2, 1, 1));
+ Scatter2DPtr mult = bookScatter2D(2, 1, 1);
+ for (size_t b = 0; b < temphisto.numPoints(); b++) {
+ const double x = temphisto.point(b).x();
+ pair<double,double> ex = temphisto.point(b).xErrs();
+ pair<double,double> ex2 = ex;
+ if(ex2.first ==0.) ex2. first=0.0001;
+ if(ex2.second==0.) ex2.second=0.0001;
+ if (inRange(sqr(sqrtS()/GeV), x-ex2.first, x+ex2.second)) {
+ mult->addPoint(x, sigma, ex, make_pair(error,error));
+ }
+ else {
+ mult->addPoint(x, 0., ex, make_pair(0.,.0));
+ }
+ }
+ }
+
+ //@}
+
+
+ /// @name Histograms
+ //@{
+ CounterPtr _npion;
+ //@}
+
+
+ };
+
+
+ // The hook for the plugin system
+ DECLARE_RIVET_PLUGIN(KLOE_2005_I655225);
+
+
+}
diff --git a/analyses/pluginFrascati/KLOE_2005_I655225.info b/analyses/pluginFrascati/KLOE_2005_I655225.info
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/KLOE_2005_I655225.info
@@ -0,0 +1,34 @@
+Name: KLOE_2005_I655225
+Year: 2005
+Summary: Cross section for $e^+e^-\to\pi^+\pi^-$ below 1 GeV
+Experiment: KLOE
+Collider: DAPHNE
+InspireID: 655225
+Status: VALIDATED
+Authors:
+ - Peter Richardson <peter.richardson@durham.ac.uk>
+References:
+ - Phys.Lett. B606 (2005) 12-24, 2005
+RunInfo: e+e- to hadrons
+NeedCrossSection: yes
+Beams: [e+,e-]
+Luminosity_fb:
+Description:
+ 'Measurement of the cross section for $e^+e^-\to\pi^+\pi^-$ below 1 GeV'
+BibKey: Aloisio:2004bu
+BibTeX: '@article{Aloisio:2004bu,
+ author = "Aloisio, A. and others",
+ title = "{Measurement of $\sigma(e^+e^- \to \pi^+ \pi^- \gamma$)
+ and extraction of $\sigma(e^+e^- \to \pi^+ \pi^-$) below
+ 1-GeV with the KLOE detector}",
+ collaboration = "KLOE",
+ journal = "Phys. Lett.",
+ volume = "B606",
+ year = "2005",
+ pages = "12-24",
+ doi = "10.1016/j.physletb.2004.11.068",
+ eprint = "hep-ex/0407048",
+ archivePrefix = "arXiv",
+ primaryClass = "hep-ex",
+ SLACcitation = "%%CITATION = HEP-EX/0407048;%%"
+}'
diff --git a/analyses/pluginFrascati/KLOE_2005_I655225.plot b/analyses/pluginFrascati/KLOE_2005_I655225.plot
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/KLOE_2005_I655225.plot
@@ -0,0 +1,6 @@
+BEGIN PLOT /KLOE_2005_I655225/d02-x01-y01
+Title=$\sigma(e^+e^-\to \pi^+\pi^-)$
+XLabel=$s/\text{GeV}^2$
+YLabel=$\sigma(e^+e^-\to \pi^+\pi^-)$/nb
+ConnectGaps=1
+END PLOT
diff --git a/analyses/pluginFrascati/KLOE_2005_I655225.yoda b/analyses/pluginFrascati/KLOE_2005_I655225.yoda
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/KLOE_2005_I655225.yoda
@@ -0,0 +1,204 @@
+BEGIN YODA_SCATTER2D_V2 /REF/KLOE_2005_I655225/d01-x01-y01
+IsRef: 1
+Path: /REF/KLOE_2005_I655225/d01-x01-y01
+Title: doi:10.17182/hepdata.41901.v1/t1
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+3.500000e-01 0.000000e+00 0.000000e+00 1.340000e+01 2.944295e-01 2.944295e-01
+3.600000e-01 0.000000e+00 0.000000e+00 1.459000e+01 3.034545e-01 3.034545e-01
+3.700000e-01 0.000000e+00 0.000000e+00 1.578000e+01 3.129526e-01 3.129526e-01
+3.800000e-01 0.000000e+00 0.000000e+00 1.704000e+01 3.234789e-01 3.234789e-01
+3.900000e-01 0.000000e+00 0.000000e+00 1.863000e+01 3.303429e-01 3.303429e-01
+4.000000e-01 0.000000e+00 0.000000e+00 2.034000e+01 3.740614e-01 3.740614e-01
+4.100000e-01 0.000000e+00 0.000000e+00 2.264000e+01 3.750150e-01 3.750150e-01
+4.200000e-01 0.000000e+00 0.000000e+00 2.456000e+01 4.130585e-01 4.130585e-01
+4.300000e-01 0.000000e+00 0.000000e+00 2.707000e+01 4.439720e-01 4.439720e-01
+4.400000e-01 0.000000e+00 0.000000e+00 2.999000e+01 4.675498e-01 4.675498e-01
+4.500000e-01 0.000000e+00 0.000000e+00 3.265000e+01 5.010944e-01 5.010944e-01
+4.600000e-01 0.000000e+00 0.000000e+00 3.624000e+01 5.344723e-01 5.344723e-01
+4.700000e-01 0.000000e+00 0.000000e+00 4.010000e+01 5.870244e-01 5.870244e-01
+4.800000e-01 0.000000e+00 0.000000e+00 4.434000e+01 6.438927e-01 6.438927e-01
+4.900000e-01 0.000000e+00 0.000000e+00 4.894000e+01 6.829422e-01 6.829422e-01
+5.000000e-01 0.000000e+00 0.000000e+00 5.410000e+01 7.963311e-01 7.963311e-01
+5.100000e-01 0.000000e+00 0.000000e+00 5.977000e+01 8.253105e-01 8.253105e-01
+5.200000e-01 0.000000e+00 0.000000e+00 6.493000e+01 8.862148e-01 8.862148e-01
+5.300000e-01 0.000000e+00 0.000000e+00 7.024000e+01 9.600794e-01 9.600794e-01
+5.400000e-01 0.000000e+00 0.000000e+00 7.560000e+01 1.042060e+00 1.042060e+00
+5.500000e-01 0.000000e+00 0.000000e+00 8.020000e+01 1.096353e+00 1.096353e+00
+5.600000e-01 0.000000e+00 0.000000e+00 8.347000e+01 1.118567e+00 1.118567e+00
+5.700000e-01 0.000000e+00 0.000000e+00 8.606000e+01 1.146920e+00 1.146920e+00
+5.800000e-01 0.000000e+00 0.000000e+00 8.785000e+01 1.168698e+00 1.168698e+00
+5.900000e-01 0.000000e+00 0.000000e+00 8.950000e+01 1.207336e+00 1.207336e+00
+6.000000e-01 0.000000e+00 0.000000e+00 9.031000e+01 1.201564e+00 1.201564e+00
+6.100000e-01 0.000000e+00 0.000000e+00 7.420000e+01 1.007181e+00 1.007181e+00
+6.200000e-01 0.000000e+00 0.000000e+00 6.549000e+01 8.793226e-01 8.793226e-01
+6.300000e-01 0.000000e+00 0.000000e+00 6.414000e+01 8.630517e-01 8.630517e-01
+6.400000e-01 0.000000e+00 0.000000e+00 6.209000e+01 8.319479e-01 8.319479e-01
+6.500000e-01 0.000000e+00 0.000000e+00 5.940000e+01 8.062223e-01 8.062223e-01
+6.600000e-01 0.000000e+00 0.000000e+00 5.638000e+01 7.566704e-01 7.566704e-01
+6.700000e-01 0.000000e+00 0.000000e+00 5.304000e+01 7.131936e-01 7.131936e-01
+6.800000e-01 0.000000e+00 0.000000e+00 4.987000e+01 6.859276e-01 6.859276e-01
+6.900000e-01 0.000000e+00 0.000000e+00 4.698000e+01 6.371448e-01 6.371448e-01
+7.000000e-01 0.000000e+00 0.000000e+00 4.416000e+01 6.000143e-01 6.000143e-01
+7.100000e-01 0.000000e+00 0.000000e+00 4.154000e+01 5.618208e-01 5.618208e-01
+7.200000e-01 0.000000e+00 0.000000e+00 3.905000e+01 5.395685e-01 5.395685e-01
+7.300000e-01 0.000000e+00 0.000000e+00 3.687000e+01 4.991215e-01 4.991215e-01
+7.400000e-01 0.000000e+00 0.000000e+00 3.520000e+01 4.828297e-01 4.828297e-01
+7.500000e-01 0.000000e+00 0.000000e+00 3.322000e+01 4.520819e-01 4.520819e-01
+7.600000e-01 0.000000e+00 0.000000e+00 3.199000e+01 4.374750e-01 4.374750e-01
+7.700000e-01 0.000000e+00 0.000000e+00 3.051000e+01 4.239096e-01 4.239096e-01
+7.800000e-01 0.000000e+00 0.000000e+00 2.960000e+01 4.093140e-01 4.093140e-01
+7.900000e-01 0.000000e+00 0.000000e+00 2.852000e+01 3.855765e-01 3.855765e-01
+8.000000e-01 0.000000e+00 0.000000e+00 2.753000e+01 3.773327e-01 3.773327e-01
+8.100000e-01 0.000000e+00 0.000000e+00 2.700000e+01 3.710768e-01 3.710768e-01
+8.200000e-01 0.000000e+00 0.000000e+00 2.648000e+01 3.612379e-01 3.612379e-01
+8.300000e-01 0.000000e+00 0.000000e+00 2.584000e+01 3.614807e-01 3.614807e-01
+8.400000e-01 0.000000e+00 0.000000e+00 2.545000e+01 3.490384e-01 3.490384e-01
+8.500000e-01 0.000000e+00 0.000000e+00 2.516000e+01 3.456156e-01 3.456156e-01
+8.600000e-01 0.000000e+00 0.000000e+00 2.496000e+01 3.395972e-01 3.395972e-01
+8.700000e-01 0.000000e+00 0.000000e+00 2.481000e+01 3.495953e-01 3.495953e-01
+8.800000e-01 0.000000e+00 0.000000e+00 2.509000e+01 3.486837e-01 3.486837e-01
+8.900000e-01 0.000000e+00 0.000000e+00 2.517000e+01 3.420989e-01 3.420989e-01
+9.000000e-01 0.000000e+00 0.000000e+00 2.537000e+01 3.480936e-01 3.480936e-01
+9.100000e-01 0.000000e+00 0.000000e+00 2.586000e+01 3.503367e-01 3.503367e-01
+9.200000e-01 0.000000e+00 0.000000e+00 2.687000e+01 3.695450e-01 3.695450e-01
+9.300000e-01 0.000000e+00 0.000000e+00 2.794000e+01 3.821835e-01 3.821835e-01
+9.400000e-01 0.000000e+00 0.000000e+00 2.949000e+01 4.080257e-01 4.080257e-01
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/KLOE_2005_I655225/d02-x01-y01
+IsRef: 1
+Path: /REF/KLOE_2005_I655225/d02-x01-y01
+Title: doi:10.17182/hepdata.41901.v1/t2
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+3.500000e-01 0.000000e+00 0.000000e+00 3.300000e+02 8.163443e+00 8.163443e+00
+3.600000e-01 0.000000e+00 0.000000e+00 3.490000e+02 8.290462e+00 8.290462e+00
+3.700000e-01 0.000000e+00 0.000000e+00 3.700000e+02 8.436694e+00 8.436694e+00
+3.800000e-01 0.000000e+00 0.000000e+00 3.920000e+02 7.803433e+00 7.803433e+00
+3.900000e-01 0.000000e+00 0.000000e+00 4.160000e+02 8.002192e+00 8.002192e+00
+4.000000e-01 0.000000e+00 0.000000e+00 4.500000e+02 9.044612e+00 9.044612e+00
+4.100000e-01 0.000000e+00 0.000000e+00 4.890000e+02 8.645091e+00 8.645091e+00
+4.200000e-01 0.000000e+00 0.000000e+00 5.210000e+02 9.642274e+00 9.642274e+00
+4.300000e-01 0.000000e+00 0.000000e+00 5.640000e+02 1.002654e+01 1.002654e+01
+4.400000e-01 0.000000e+00 0.000000e+00 6.080000e+02 1.043482e+01 1.043482e+01
+4.500000e-01 0.000000e+00 0.000000e+00 6.490000e+02 1.082749e+01 1.082749e+01
+4.600000e-01 0.000000e+00 0.000000e+00 7.100000e+02 1.143084e+01 1.143084e+01
+4.700000e-01 0.000000e+00 0.000000e+00 7.690000e+02 1.203331e+01 1.203331e+01
+4.800000e-01 0.000000e+00 0.000000e+00 8.300000e+02 1.267288e+01 1.267288e+01
+4.900000e-01 0.000000e+00 0.000000e+00 8.950000e+02 1.337034e+01 1.337034e+01
+5.000000e-01 0.000000e+00 0.000000e+00 9.670000e+02 1.467939e+01 1.467939e+01
+5.100000e-01 0.000000e+00 0.000000e+00 1.041000e+03 1.498520e+01 1.498520e+01
+5.200000e-01 0.000000e+00 0.000000e+00 1.102000e+03 1.567589e+01 1.567589e+01
+5.300000e-01 0.000000e+00 0.000000e+00 1.171000e+03 1.691570e+01 1.691570e+01
+5.400000e-01 0.000000e+00 0.000000e+00 1.226000e+03 1.753563e+01 1.753563e+01
+5.500000e-01 0.000000e+00 0.000000e+00 1.279000e+03 1.813853e+01 1.813853e+01
+5.600000e-01 0.000000e+00 0.000000e+00 1.288000e+03 1.782551e+01 1.782551e+01
+5.700000e-01 0.000000e+00 0.000000e+00 1.302000e+03 1.798953e+01 1.798953e+01
+5.800000e-01 0.000000e+00 0.000000e+00 1.297000e+03 1.793092e+01 1.793092e+01
+5.900000e-01 0.000000e+00 0.000000e+00 1.282000e+03 1.775531e+01 1.775531e+01
+6.000000e-01 0.000000e+00 0.000000e+00 1.266000e+03 1.756834e+01 1.756834e+01
+6.100000e-01 0.000000e+00 0.000000e+00 1.006000e+03 1.414036e+01 1.414036e+01
+6.200000e-01 0.000000e+00 0.000000e+00 8.570000e+02 1.199920e+01 1.199920e+01
+6.300000e-01 0.000000e+00 0.000000e+00 8.170000e+02 1.153834e+01 1.153834e+01
+6.400000e-01 0.000000e+00 0.000000e+00 7.720000e+02 1.060893e+01 1.060893e+01
+6.500000e-01 0.000000e+00 0.000000e+00 7.140000e+02 9.929096e+00 9.929096e+00
+6.600000e-01 0.000000e+00 0.000000e+00 6.570000e+02 9.269689e+00 9.269689e+00
+6.700000e-01 0.000000e+00 0.000000e+00 5.950000e+02 8.564581e+00 8.564581e+00
+6.800000e-01 0.000000e+00 0.000000e+00 5.430000e+02 7.985333e+00 7.985333e+00
+6.900000e-01 0.000000e+00 0.000000e+00 4.932000e+02 7.001135e+00 7.001135e+00
+7.000000e-01 0.000000e+00 0.000000e+00 4.470000e+02 6.385848e+00 6.385848e+00
+7.100000e-01 0.000000e+00 0.000000e+00 4.050000e+02 5.773392e+00 5.773392e+00
+7.200000e-01 0.000000e+00 0.000000e+00 3.671000e+02 5.347103e+00 5.347103e+00
+7.300000e-01 0.000000e+00 0.000000e+00 3.333000e+02 4.778745e+00 4.778745e+00
+7.400000e-01 0.000000e+00 0.000000e+00 3.046000e+02 4.362402e+00 4.362402e+00
+7.500000e-01 0.000000e+00 0.000000e+00 2.778000e+02 3.967619e+00 3.967619e+00
+7.600000e-01 0.000000e+00 0.000000e+00 2.572000e+02 3.688712e+00 3.688712e+00
+7.700000e-01 0.000000e+00 0.000000e+00 2.338000e+02 3.427144e+00 3.427144e+00
+7.800000e-01 0.000000e+00 0.000000e+00 2.177000e+02 3.199643e+00 3.199643e+00
+7.900000e-01 0.000000e+00 0.000000e+00 2.003000e+02 2.861722e+00 2.861722e+00
+8.000000e-01 0.000000e+00 0.000000e+00 1.845000e+02 2.684124e+00 2.684124e+00
+8.100000e-01 0.000000e+00 0.000000e+00 1.721000e+02 2.497635e+00 2.497635e+00
+8.200000e-01 0.000000e+00 0.000000e+00 1.600000e+02 2.314563e+00 2.314563e+00
+8.300000e-01 0.000000e+00 0.000000e+00 1.484000e+02 2.185785e+00 2.185785e+00
+8.400000e-01 0.000000e+00 0.000000e+00 1.385000e+02 2.026703e+00 2.026703e+00
+8.500000e-01 0.000000e+00 0.000000e+00 1.292000e+02 1.874622e+00 1.874622e+00
+8.600000e-01 0.000000e+00 0.000000e+00 1.203000e+02 1.727564e+00 1.727564e+00
+8.700000e-01 0.000000e+00 0.000000e+00 1.118000e+02 1.683709e+00 1.683709e+00
+8.800000e-01 0.000000e+00 0.000000e+00 1.063000e+02 1.571798e+00 1.571798e+00
+8.900000e-01 0.000000e+00 0.000000e+00 9.950000e+01 1.447011e+00 1.447011e+00
+9.000000e-01 0.000000e+00 0.000000e+00 9.310000e+01 1.328214e+00 1.328214e+00
+9.100000e-01 0.000000e+00 0.000000e+00 8.760000e+01 1.266155e+00 1.266155e+00
+9.200000e-01 0.000000e+00 0.000000e+00 8.300000e+01 1.214915e+00 1.214915e+00
+9.300000e-01 0.000000e+00 0.000000e+00 7.910000e+01 1.124101e+00 1.124101e+00
+9.400000e-01 0.000000e+00 0.000000e+00 7.530000e+01 1.080997e+00 1.080997e+00
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/KLOE_2005_I655225/d03-x01-y01
+IsRef: 1
+Path: /REF/KLOE_2005_I655225/d03-x01-y01
+Title: doi:10.17182/hepdata.41901.v1/t3
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+3.500000e-01 0.000000e+00 0.000000e+00 7.680000e+00 1.600000e-01 1.600000e-01
+3.600000e-01 0.000000e+00 0.000000e+00 8.260000e+00 1.600000e-01 1.600000e-01
+3.700000e-01 0.000000e+00 0.000000e+00 8.920000e+00 1.600000e-01 1.600000e-01
+3.800000e-01 0.000000e+00 0.000000e+00 9.600000e+00 1.600000e-01 1.600000e-01
+3.900000e-01 0.000000e+00 0.000000e+00 1.035000e+01 1.500000e-01 1.500000e-01
+4.000000e-01 0.000000e+00 0.000000e+00 1.140000e+01 1.700000e-01 1.700000e-01
+4.100000e-01 0.000000e+00 0.000000e+00 1.259000e+01 1.600000e-01 1.600000e-01
+4.200000e-01 0.000000e+00 0.000000e+00 1.363000e+01 1.800000e-01 1.800000e-01
+4.300000e-01 0.000000e+00 0.000000e+00 1.501000e+01 1.800000e-01 1.800000e-01
+4.400000e-01 0.000000e+00 0.000000e+00 1.643000e+01 1.800000e-01 1.800000e-01
+4.500000e-01 0.000000e+00 0.000000e+00 1.782000e+01 1.900000e-01 1.900000e-01
+4.600000e-01 0.000000e+00 0.000000e+00 1.979000e+01 1.800000e-01 1.800000e-01
+4.700000e-01 0.000000e+00 0.000000e+00 2.178000e+01 2.000000e-01 2.000000e-01
+4.800000e-01 0.000000e+00 0.000000e+00 2.386000e+01 2.000000e-01 2.000000e-01
+4.900000e-01 0.000000e+00 0.000000e+00 2.611000e+01 2.000000e-01 2.000000e-01
+5.000000e-01 0.000000e+00 0.000000e+00 2.860000e+01 2.300000e-01 2.300000e-01
+5.100000e-01 0.000000e+00 0.000000e+00 3.123000e+01 2.200000e-01 2.200000e-01
+5.200000e-01 0.000000e+00 0.000000e+00 3.350000e+01 2.200000e-01 2.200000e-01
+5.300000e-01 0.000000e+00 0.000000e+00 3.605000e+01 2.300000e-01 2.300000e-01
+5.400000e-01 0.000000e+00 0.000000e+00 3.820000e+01 2.600000e-01 2.600000e-01
+5.500000e-01 0.000000e+00 0.000000e+00 4.032000e+01 2.400000e-01 2.400000e-01
+5.600000e-01 0.000000e+00 0.000000e+00 4.107000e+01 2.400000e-01 2.400000e-01
+5.700000e-01 0.000000e+00 0.000000e+00 4.198000e+01 2.300000e-01 2.300000e-01
+5.800000e-01 0.000000e+00 0.000000e+00 4.236000e+01 2.300000e-01 2.300000e-01
+5.900000e-01 0.000000e+00 0.000000e+00 4.246000e+01 2.400000e-01 2.400000e-01
+6.000000e-01 0.000000e+00 0.000000e+00 4.258000e+01 2.300000e-01 2.300000e-01
+6.100000e-01 0.000000e+00 0.000000e+00 3.243000e+01 2.000000e-01 2.000000e-01
+6.200000e-01 0.000000e+00 0.000000e+00 2.799000e+01 1.600000e-01 1.600000e-01
+6.300000e-01 0.000000e+00 0.000000e+00 2.732000e+01 1.600000e-01 1.600000e-01
+6.400000e-01 0.000000e+00 0.000000e+00 2.627000e+01 1.500000e-01 1.500000e-01
+6.500000e-01 0.000000e+00 0.000000e+00 2.469000e+01 1.500000e-01 1.500000e-01
+6.600000e-01 0.000000e+00 0.000000e+00 2.305000e+01 1.400000e-01 1.400000e-01
+6.700000e-01 0.000000e+00 0.000000e+00 2.118000e+01 1.300000e-01 1.300000e-01
+6.800000e-01 0.000000e+00 0.000000e+00 1.957000e+01 1.300000e-01 1.300000e-01
+6.900000e-01 0.000000e+00 0.000000e+00 1.802000e+01 1.100000e-01 1.100000e-01
+7.000000e-01 0.000000e+00 0.000000e+00 1.654000e+01 1.100000e-01 1.100000e-01
+7.100000e-01 0.000000e+00 0.000000e+00 1.517000e+01 1.000000e-01 1.000000e-01
+7.200000e-01 0.000000e+00 0.000000e+00 1.392000e+01 1.000000e-01 1.000000e-01
+7.300000e-01 0.000000e+00 0.000000e+00 1.278000e+01 8.000000e-02 8.000000e-02
+7.400000e-01 0.000000e+00 0.000000e+00 1.181000e+01 8.000000e-02 8.000000e-02
+7.500000e-01 0.000000e+00 0.000000e+00 1.089000e+01 7.000000e-02 7.000000e-02
+7.600000e-01 0.000000e+00 0.000000e+00 1.019000e+01 7.000000e-02 7.000000e-02
+7.700000e-01 0.000000e+00 0.000000e+00 9.370000e+00 7.000000e-02 7.000000e-02
+7.800000e-01 0.000000e+00 0.000000e+00 8.820000e+00 6.000000e-02 6.000000e-02
+7.900000e-01 0.000000e+00 0.000000e+00 8.200000e+00 5.000000e-02 5.000000e-02
+8.000000e-01 0.000000e+00 0.000000e+00 7.630000e+00 5.000000e-02 5.000000e-02
+8.100000e-01 0.000000e+00 0.000000e+00 7.200000e+00 5.000000e-02 5.000000e-02
+8.200000e-01 0.000000e+00 0.000000e+00 6.760000e+00 5.000000e-02 5.000000e-02
+8.300000e-01 0.000000e+00 0.000000e+00 6.330000e+00 5.000000e-02 5.000000e-02
+8.400000e-01 0.000000e+00 0.000000e+00 5.970000e+00 4.000000e-02 4.000000e-02
+8.500000e-01 0.000000e+00 0.000000e+00 5.630000e+00 4.000000e-02 4.000000e-02
+8.600000e-01 0.000000e+00 0.000000e+00 5.290000e+00 4.000000e-02 4.000000e-02
+8.700000e-01 0.000000e+00 0.000000e+00 4.970000e+00 4.000000e-02 4.000000e-02
+8.800000e-01 0.000000e+00 0.000000e+00 4.774000e+00 3.500000e-02 3.500000e-02
+8.900000e-01 0.000000e+00 0.000000e+00 4.516000e+00 3.000000e-02 3.000000e-02
+9.000000e-01 0.000000e+00 0.000000e+00 4.269000e+00 3.000000e-02 3.000000e-02
+9.100000e-01 0.000000e+00 0.000000e+00 4.059000e+00 2.700000e-02 2.700000e-02
+9.200000e-01 0.000000e+00 0.000000e+00 3.886000e+00 2.600000e-02 2.600000e-02
+9.300000e-01 0.000000e+00 0.000000e+00 3.741000e+00 2.500000e-02 2.500000e-02
+9.400000e-01 0.000000e+00 0.000000e+00 3.599000e+00 2.500000e-02 2.500000e-02
+END YODA_SCATTER2D_V2
diff --git a/analyses/pluginFrascati/KLOE_2008_I791841.cc b/analyses/pluginFrascati/KLOE_2008_I791841.cc
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/KLOE_2008_I791841.cc
@@ -0,0 +1,95 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/Projections/FinalState.hh"
+
+namespace Rivet {
+
+
+ /// @brief Add a short analysis description here
+ class KLOE_2008_I791841 : public Analysis {
+ public:
+
+ /// Constructor
+ DEFAULT_RIVET_ANALYSIS_CTOR(KLOE_2008_I791841);
+
+
+ /// @name Analysis methods
+ //@{
+
+ /// Book histograms and initialise projections before the run
+ void init() {
+ // Initialise and register projections
+ declare(FinalState(), "FS");
+ _n4pi = bookCounter("TMP/4pi");
+ _n2pigamma = bookCounter("TMP/2pigamma");
+ }
+
+
+ /// Perform the per-event analysis
+ void analyze(const Event& event) {
+ const FinalState& fs = apply<FinalState>(event, "FS");
+
+ map<long,int> nCount;
+ int ntotal(0);
+ foreach (const Particle& p, fs.particles()) {
+ nCount[p.pdgId()] += 1;
+ ++ntotal;
+ }
+ if(nCount[111]==2) {
+ if( nCount[211] == 1 && nCount[-211] == 1 )
+ _n4pi->fill(event.weight());
+ else if( nCount[22] == 1)
+ _n2pigamma->fill(event.weight());
+ }
+ }
+
+
+ /// Normalise histograms etc., after the run
+ void finalize() {
+ for(unsigned int ix=1;ix<3;++ix) {
+ double sigma,error;
+ if(ix==1) {
+ sigma = _n4pi->val();
+ error = _n4pi->err();
+ }
+ else if(ix==2) {
+ sigma = _n2pigamma->val();
+ error = _n2pigamma->err();
+ }
+ sigma *= crossSection()/ sumOfWeights() /nanobarn;
+ error *= crossSection()/ sumOfWeights() /nanobarn;
+ Scatter2D temphisto(refData(ix, 1, 1));
+ Scatter2DPtr mult = bookScatter2D(ix, 1, 1);
+ for (size_t b = 0; b < temphisto.numPoints(); b++) {
+ const double x = temphisto.point(b).x();
+ pair<double,double> ex = temphisto.point(b).xErrs();
+ pair<double,double> ex2 = ex;
+ if(ex2.first ==0.) ex2. first=0.0001;
+ if(ex2.second==0.) ex2.second=0.0001;
+ if (inRange(sqrtS()/MeV, x-ex2.first, x+ex2.second)) {
+ mult->addPoint(x, sigma, ex, make_pair(error,error));
+ }
+ else {
+ mult->addPoint(x, 0., ex, make_pair(0.,.0));
+ }
+ }
+ }
+ }
+
+ //@}
+
+
+ /// @name Histograms
+ //@{
+ CounterPtr _n4pi,_n2pigamma;
+ //@}
+
+
+ };
+
+
+ // The hook for the plugin system
+ DECLARE_RIVET_PLUGIN(KLOE_2008_I791841);
+
+
+}
diff --git a/analyses/pluginFrascati/KLOE_2008_I791841.info b/analyses/pluginFrascati/KLOE_2008_I791841.info
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/KLOE_2008_I791841.info
@@ -0,0 +1,34 @@
+Name: KLOE_2008_I791841
+Year: 2008
+Summary: Cross section for $\pi^+\pi^-2\pi^0$ and $2\pi^0\gamma$ near the $\phi$ mass
+Experiment: KLOE
+Collider: DAPHNE
+InspireID: 791841
+Status: VALIDATED
+Authors:
+ - Peter Richardson <peter.richardson@durham.ac.uk>
+References:
+ - Phys.Lett. B669 (2008) 223-228, 2008
+RunInfo: e+e- to hadrons
+NeedCrossSection: yes
+Beams: [e+,e-]
+Luminosity_fb:
+Description:
+ 'Measurement of the cross section for $\pi^+\pi^-2\pi^0$ and $2\pi^0\gamma$ near the $\phi$ mass'
+Keywords: []
+BibKey: Ambrosino:2008gb
+BibTeX: '@article{Ambrosino:2008gb,
+ author = "Ambrosino, F. and others",
+ title = "{Study of the process $e^+ e^- \to \omega \pi^0$ in the
+ $\phi$-meson mass region with the KLOE detector}",
+ collaboration = "KLOE",
+ journal = "Phys. Lett.",
+ volume = "B669",
+ year = "2008",
+ pages = "223-228",
+ doi = "10.1016/j.physletb.2008.09.056",
+ eprint = "0807.4909",
+ archivePrefix = "arXiv",
+ primaryClass = "hep-ex",
+ SLACcitation = "%%CITATION = ARXIV:0807.4909;%%"
+}'
diff --git a/analyses/pluginFrascati/KLOE_2008_I791841.plot b/analyses/pluginFrascati/KLOE_2008_I791841.plot
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/KLOE_2008_I791841.plot
@@ -0,0 +1,16 @@
+BEGIN PLOT /KLOE_2008_I791841/d01-x01-y01
+Title=$\sigma(e^+e^-\to \pi^+\pi^-2\pi^0)$
+XLabel=$\sqrt{s}$/MeV
+YLabel=$\sigma(e^+e^-\to \pi^+\pi^-2\pi^0)$/nb
+LogY=0
+LegendXPos=0.8
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /KLOE_2008_I791841/d02-x01-y01
+Title=$\sigma(e^+e^-\to 2\pi^0\gamma)$
+XLabel=$\sqrt{s}$/MeV
+YLabel=$\sigma(e^+e^-\to 2\pi^0\gamma)$/nb
+LogY=0
+LegendXPos=0.8
+ConnectGaps=1
+END PLOT
diff --git a/analyses/pluginFrascati/KLOE_2008_I791841.yoda b/analyses/pluginFrascati/KLOE_2008_I791841.yoda
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/KLOE_2008_I791841.yoda
@@ -0,0 +1,50 @@
+BEGIN YODA_SCATTER2D_V2 /REF/KLOE_2008_I791841/d01-x01-y01
+IsRef: 1
+Path: /REF/KLOE_2008_I791841/d01-x01-y01
+Title: doi:10.17182/hepdata.50498.v1/t1
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.000100e+03 0.000000e+00 0.000000e+00 5.720000e+00 5.000000e-02 5.000000e-02
+1.009900e+03 0.000000e+00 0.000000e+00 6.200000e+00 6.000000e-02 6.000000e-02
+1.017200e+03 0.000000e+00 0.000000e+00 5.710000e+00 8.000000e-02 8.000000e-02
+1.018150e+03 0.000000e+00 0.000000e+00 5.600000e+00 6.000000e-02 6.000000e-02
+1.019300e+03 0.000000e+00 0.000000e+00 5.880000e+00 1.200000e-01 1.200000e-01
+1.019450e+03 0.000000e+00 0.000000e+00 5.890000e+00 6.000000e-02 6.000000e-02
+1.019550e+03 0.000000e+00 0.000000e+00 5.930000e+00 5.000000e-02 5.000000e-02
+1.019650e+03 0.000000e+00 0.000000e+00 5.980000e+00 5.000000e-02 5.000000e-02
+1.019750e+03 0.000000e+00 0.000000e+00 6.040000e+00 5.000000e-02 5.000000e-02
+1.019850e+03 0.000000e+00 0.000000e+00 6.080000e+00 5.000000e-02 5.000000e-02
+1.019950e+03 0.000000e+00 0.000000e+00 6.200000e+00 7.000000e-02 7.000000e-02
+1.020050e+03 0.000000e+00 0.000000e+00 6.210000e+00 8.000000e-02 8.000000e-02
+1.020150e+03 0.000000e+00 0.000000e+00 6.230000e+00 1.100000e-01 1.100000e-01
+1.020450e+03 0.000000e+00 0.000000e+00 6.410000e+00 9.000000e-02 9.000000e-02
+1.022300e+03 0.000000e+00 0.000000e+00 7.240000e+00 8.000000e-02 8.000000e-02
+1.023000e+03 0.000000e+00 0.000000e+00 7.410000e+00 7.000000e-02 7.000000e-02
+1.029950e+03 0.000000e+00 0.000000e+00 7.840000e+00 7.000000e-02 7.000000e-02
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/KLOE_2008_I791841/d02-x01-y01
+IsRef: 1
+Path: /REF/KLOE_2008_I791841/d02-x01-y01
+Title: doi:10.17182/hepdata.50498.v1/t2
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.000100e+03 0.000000e+00 0.000000e+00 5.500000e-01 5.000000e-03 5.000000e-03
+1.009900e+03 0.000000e+00 0.000000e+00 5.810000e-01 1.200000e-02 1.200000e-02
+1.017200e+03 0.000000e+00 0.000000e+00 5.640000e-01 1.800000e-02 1.800000e-02
+1.018150e+03 0.000000e+00 0.000000e+00 5.410000e-01 1.400000e-02 1.400000e-02
+1.019300e+03 0.000000e+00 0.000000e+00 4.800000e-01 3.400000e-02 3.400000e-02
+1.019450e+03 0.000000e+00 0.000000e+00 4.970000e-01 9.000000e-03 9.000000e-03
+1.019550e+03 0.000000e+00 0.000000e+00 4.870000e-01 8.000000e-03 8.000000e-03
+1.019650e+03 0.000000e+00 0.000000e+00 5.090000e-01 7.000000e-03 7.000000e-03
+1.019750e+03 0.000000e+00 0.000000e+00 5.050000e-01 6.000000e-03 6.000000e-03
+1.019850e+03 0.000000e+00 0.000000e+00 5.080000e-01 6.000000e-03 6.000000e-03
+1.019950e+03 0.000000e+00 0.000000e+00 4.910000e-01 1.100000e-02 1.100000e-02
+1.020050e+03 0.000000e+00 0.000000e+00 5.160000e-01 1.600000e-02 1.600000e-02
+1.020150e+03 0.000000e+00 0.000000e+00 5.010000e-01 2.400000e-02 2.400000e-02
+1.020450e+03 0.000000e+00 0.000000e+00 4.880000e-01 2.400000e-02 2.400000e-02
+1.022300e+03 0.000000e+00 0.000000e+00 6.120000e-01 1.800000e-02 1.800000e-02
+1.023000e+03 0.000000e+00 0.000000e+00 6.190000e-01 1.300000e-02 1.300000e-02
+1.029950e+03 0.000000e+00 0.000000e+00 6.890000e-01 1.300000e-02 1.300000e-02
+END YODA_SCATTER2D_V2
diff --git a/analyses/pluginFrascati/KLOE_2009_I797438.cc b/analyses/pluginFrascati/KLOE_2009_I797438.cc
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/KLOE_2009_I797438.cc
@@ -0,0 +1,82 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/Projections/FinalState.hh"
+
+namespace Rivet {
+
+
+ /// @brief Add a short analysis description here
+ class KLOE_2009_I797438 : public Analysis {
+ public:
+
+ /// Constructor
+ DEFAULT_RIVET_ANALYSIS_CTOR(KLOE_2009_I797438);
+
+
+ /// @name Analysis methods
+ //@{
+
+ /// Book histograms and initialise projections before the run
+ void init() {
+
+ // Initialise and register projections
+ declare(FinalState(), "FS");
+
+ // Book histograms
+ _npion = bookCounter("TMP/pion");
+
+ }
+
+
+ /// Perform the per-event analysis
+ void analyze(const Event& event) {
+ const FinalState& fs = apply<FinalState>(event, "FS");
+ if(fs.particles().size()!=2) vetoEvent;
+ foreach (const Particle& p, fs.particles()) {
+ if(abs(p.pdgId())!=PID::PIPLUS) vetoEvent;
+ }
+ _npion->fill(event.weight());
+ }
+
+
+ /// Normalise histograms etc., after the run
+ void finalize() {
+ double sigma = _npion->val();
+ double error = _npion->err();
+ sigma *= crossSection()/ sumOfWeights() /nanobarn;
+ error *= crossSection()/ sumOfWeights() /nanobarn;
+ Scatter2D temphisto(refData(2, 1, 1));
+ Scatter2DPtr mult = bookScatter2D(2, 1, 1);
+ for (size_t b = 0; b < temphisto.numPoints(); b++) {
+ const double x = temphisto.point(b).x();
+ pair<double,double> ex = temphisto.point(b).xErrs();
+ pair<double,double> ex2 = ex;
+ if(ex2.first ==0.) ex2. first=0.0001;
+ if(ex2.second==0.) ex2.second=0.0001;
+ if (inRange(sqr(sqrtS()/GeV), x-ex2.first, x+ex2.second)) {
+ mult->addPoint(x, sigma, ex, make_pair(error,error));
+ }
+ else {
+ mult->addPoint(x, 0., ex, make_pair(0.,.0));
+ }
+ }
+
+ }
+
+ //@}
+
+
+ /// @name Histograms
+ //@{
+ CounterPtr _npion;
+ //@}
+
+
+ };
+
+
+ // The hook for the plugin system
+ DECLARE_RIVET_PLUGIN(KLOE_2009_I797438);
+
+
+}
diff --git a/analyses/pluginFrascati/KLOE_2009_I797438.info b/analyses/pluginFrascati/KLOE_2009_I797438.info
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/KLOE_2009_I797438.info
@@ -0,0 +1,35 @@
+Name: KLOE_2009_I797438
+Year: 2009
+Summary: Cross section for $e^+e^-\to\pi^+\pi^-$ below 0.93 GeV
+Experiment: KLOE
+Collider: DAPHNE
+InspireID: 797438
+Status: VALIDATED
+Authors:
+ - Peter Richardson <peter.richardson@durham.ac.uk>
+References:
+ - Phys.Lett. B670 (2009) 285-291, 2009
+RunInfo: e+e- to hadrons
+NeedCrossSection: yes
+Beams: [e+,e-]
+Luminosity_fb:
+Description:
+ 'Measurement of the cross section for $e^+e^-\to\pi^+\pi^-$ below 0.93 GeV'
+Keywords: []
+BibKey: Ambrosino:2008aa
+BibTeX: '@article{Ambrosino:2008aa,
+ author = "Ambrosino, F. and others",
+ title = "{Measurement of $\sigma(e^+ e^- \to \pi^+ \pi^-
+ \gamma(\gamma)$ and the dipion contribution to the muon
+ anomaly with the KLOE detector}",
+ collaboration = "KLOE",
+ journal = "Phys. Lett.",
+ volume = "B670",
+ year = "2009",
+ pages = "285-291",
+ doi = "10.1016/j.physletb.2008.10.060",
+ eprint = "0809.3950",
+ archivePrefix = "arXiv",
+ primaryClass = "hep-ex",
+ SLACcitation = "%%CITATION = ARXIV:0809.3950;%%"
+}'
diff --git a/analyses/pluginFrascati/KLOE_2009_I797438.plot b/analyses/pluginFrascati/KLOE_2009_I797438.plot
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/KLOE_2009_I797438.plot
@@ -0,0 +1,7 @@
+BEGIN PLOT /KLOE_2009_I797438/d02-x01-y01
+Title=$\sigma(e^+e^-\to \pi^+\pi^-)$
+XLabel=$s/\text{GeV}^2$
+YLabel=$\sigma(e^+e^-\to \pi^+\pi^-)$/nb
+ConnectGaps=1
+LogY=0
+END PLOT
diff --git a/analyses/pluginFrascati/KLOE_2009_I797438.yoda b/analyses/pluginFrascati/KLOE_2009_I797438.yoda
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/KLOE_2009_I797438.yoda
@@ -0,0 +1,204 @@
+BEGIN YODA_SCATTER2D_V2 /REF/KLOE_2009_I797438/d01-x01-y01
+IsRef: 1
+Path: /REF/KLOE_2009_I797438/d01-x01-y01
+Title: doi:10.17182/hepdata.57088.v1/t1
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+3.550000e-01 5.000000e-03 5.000000e-03 1.307000e+01 1.937982e-01 1.937982e-01
+3.650000e-01 5.000000e-03 5.000000e-03 1.421000e+01 1.983200e-01 1.983200e-01
+3.750000e-01 5.000000e-03 5.000000e-03 1.520000e+01 1.963116e-01 1.963116e-01
+3.850000e-01 5.000000e-03 5.000000e-03 1.660000e+01 2.005110e-01 2.005110e-01
+3.950000e-01 5.000000e-03 5.000000e-03 1.823000e+01 2.094158e-01 2.094158e-01
+4.050000e-01 5.000000e-03 5.000000e-03 1.997000e+01 2.057903e-01 2.057903e-01
+4.150000e-01 5.000000e-03 5.000000e-03 2.200000e+01 2.251222e-01 2.251222e-01
+4.250000e-01 5.000000e-03 5.000000e-03 2.409000e+01 2.308111e-01 2.308111e-01
+4.350000e-01 5.000000e-03 5.000000e-03 2.657000e+01 2.222003e-01 2.222003e-01
+4.450000e-01 5.000000e-03 5.000000e-03 2.926000e+01 2.337944e-01 2.337944e-01
+4.550000e-01 5.000000e-03 5.000000e-03 3.256000e+01 2.437691e-01 2.437691e-01
+4.650000e-01 5.000000e-03 5.000000e-03 3.560000e+01 2.529465e-01 2.529465e-01
+4.750000e-01 5.000000e-03 5.000000e-03 3.918000e+01 2.784725e-01 2.784725e-01
+4.850000e-01 5.000000e-03 5.000000e-03 4.428000e+01 2.779454e-01 2.779454e-01
+4.950000e-01 5.000000e-03 5.000000e-03 4.973000e+01 2.935000e-01 2.935000e-01
+5.050000e-01 5.000000e-03 5.000000e-03 5.417000e+01 3.360901e-01 3.360901e-01
+5.150000e-01 5.000000e-03 5.000000e-03 5.920000e+01 3.390990e-01 3.390990e-01
+5.250000e-01 5.000000e-03 5.000000e-03 6.390000e+01 3.127544e-01 3.127544e-01
+5.350000e-01 5.000000e-03 5.000000e-03 6.982000e+01 3.185490e-01 3.185490e-01
+5.450000e-01 5.000000e-03 5.000000e-03 7.468000e+01 4.044316e-01 4.044316e-01
+5.550000e-01 5.000000e-03 5.000000e-03 7.920000e+01 3.813358e-01 3.813358e-01
+5.650000e-01 5.000000e-03 5.000000e-03 8.379000e+01 4.009872e-01 4.009872e-01
+5.750000e-01 5.000000e-03 5.000000e-03 8.579000e+01 3.588026e-01 3.588026e-01
+5.850000e-01 5.000000e-03 5.000000e-03 8.866000e+01 5.089351e-01 5.089351e-01
+5.950000e-01 5.000000e-03 5.000000e-03 9.024000e+01 4.997949e-01 4.997949e-01
+6.050000e-01 5.000000e-03 5.000000e-03 9.138000e+01 1.948098e+00 1.948098e+00
+6.150000e-01 5.000000e-03 5.000000e-03 7.010000e+01 2.824060e+00 2.824060e+00
+6.250000e-01 5.000000e-03 5.000000e-03 6.502000e+01 4.084458e-01 4.084458e-01
+6.350000e-01 5.000000e-03 5.000000e-03 6.492000e+01 2.791621e-01 2.791621e-01
+6.450000e-01 5.000000e-03 5.000000e-03 6.240000e+01 2.667397e-01 2.667397e-01
+6.550000e-01 5.000000e-03 5.000000e-03 5.962000e+01 2.868956e-01 2.868956e-01
+6.650000e-01 5.000000e-03 5.000000e-03 5.628000e+01 2.712503e-01 2.712503e-01
+6.750000e-01 5.000000e-03 5.000000e-03 5.343000e+01 2.636512e-01 2.636512e-01
+6.850000e-01 5.000000e-03 5.000000e-03 4.984000e+01 2.208443e-01 2.208443e-01
+6.950000e-01 5.000000e-03 5.000000e-03 4.722000e+01 2.336375e-01 2.336375e-01
+7.050000e-01 5.000000e-03 5.000000e-03 4.465000e+01 2.269738e-01 2.269738e-01
+7.150000e-01 5.000000e-03 5.000000e-03 4.140000e+01 1.990970e-01 1.990970e-01
+7.250000e-01 5.000000e-03 5.000000e-03 3.940000e+01 2.144531e-01 2.144531e-01
+7.350000e-01 5.000000e-03 5.000000e-03 3.780000e+01 1.840880e-01 1.840880e-01
+7.450000e-01 5.000000e-03 5.000000e-03 3.605000e+01 1.805437e-01 1.805437e-01
+7.550000e-01 5.000000e-03 5.000000e-03 3.413000e+01 1.689632e-01 1.689632e-01
+7.650000e-01 5.000000e-03 5.000000e-03 3.250000e+01 1.866983e-01 1.866983e-01
+7.750000e-01 5.000000e-03 5.000000e-03 3.114000e+01 1.552321e-01 1.552321e-01
+7.850000e-01 5.000000e-03 5.000000e-03 3.001000e+01 1.529902e-01 1.529902e-01
+7.950000e-01 5.000000e-03 5.000000e-03 2.923000e+01 1.436800e-01 1.436800e-01
+8.050000e-01 5.000000e-03 5.000000e-03 2.846000e+01 1.608400e-01 1.608400e-01
+8.150000e-01 5.000000e-03 5.000000e-03 2.779000e+01 1.588359e-01 1.588359e-01
+8.250000e-01 5.000000e-03 5.000000e-03 2.706000e+01 1.393644e-01 1.393644e-01
+8.350000e-01 5.000000e-03 5.000000e-03 2.643000e+01 1.678159e-01 1.678159e-01
+8.450000e-01 5.000000e-03 5.000000e-03 2.602000e+01 1.466618e-01 1.466618e-01
+8.550000e-01 5.000000e-03 5.000000e-03 2.563000e+01 1.645580e-01 1.645580e-01
+8.650000e-01 5.000000e-03 5.000000e-03 2.563000e+01 1.454897e-01 1.454897e-01
+8.750000e-01 5.000000e-03 5.000000e-03 2.563000e+01 1.454897e-01 1.454897e-01
+8.850000e-01 5.000000e-03 5.000000e-03 2.563000e+01 1.645580e-01 1.645580e-01
+8.950000e-01 5.000000e-03 5.000000e-03 2.563000e+01 1.645580e-01 1.645580e-01
+9.050000e-01 5.000000e-03 5.000000e-03 2.563000e+01 1.645580e-01 1.645580e-01
+9.150000e-01 5.000000e-03 5.000000e-03 2.563000e+01 1.852166e-01 1.852166e-01
+9.250000e-01 5.000000e-03 5.000000e-03 2.563000e+01 1.852166e-01 1.852166e-01
+9.350000e-01 5.000000e-03 5.000000e-03 2.857000e+01 2.005020e-01 2.005020e-01
+9.450000e-01 5.000000e-03 5.000000e-03 2.986000e+01 2.336257e-01 2.336257e-01
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/KLOE_2009_I797438/d02-x01-y01
+IsRef: 1
+Path: /REF/KLOE_2009_I797438/d02-x01-y01
+Title: doi:10.17182/hepdata.57088.v1/t2
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+3.550000e-01 5.000000e-03 5.000000e-03 3.090000e+02 4.762738e+00 4.762738e+00
+3.650000e-01 5.000000e-03 5.000000e-03 3.350000e+02 4.861204e+00 4.861204e+00
+3.750000e-01 5.000000e-03 5.000000e-03 3.540000e+02 4.797676e+00 4.797676e+00
+3.850000e-01 5.000000e-03 5.000000e-03 3.800000e+02 4.863456e+00 4.863456e+00
+3.950000e-01 5.000000e-03 5.000000e-03 4.050000e+02 4.835403e+00 4.835403e+00
+4.050000e-01 5.000000e-03 5.000000e-03 4.390000e+02 4.908593e+00 4.908593e+00
+4.150000e-01 5.000000e-03 5.000000e-03 4.720000e+02 5.101498e+00 5.101498e+00
+4.250000e-01 5.000000e-03 5.000000e-03 5.110000e+02 5.192984e+00 5.192984e+00
+4.350000e-01 5.000000e-03 5.000000e-03 5.480000e+02 4.970796e+00 4.970796e+00
+4.450000e-01 5.000000e-03 5.000000e-03 5.920000e+02 5.011194e+00 5.011194e+00
+4.550000e-01 5.000000e-03 5.000000e-03 6.480000e+02 5.023732e+00 5.023732e+00
+4.650000e-01 5.000000e-03 5.000000e-03 6.950000e+02 5.160092e+00 5.160092e+00
+4.750000e-01 5.000000e-03 5.000000e-03 7.490000e+02 5.580952e+00 5.580952e+00
+4.850000e-01 5.000000e-03 5.000000e-03 8.260000e+02 6.161432e+00 6.161432e+00
+4.950000e-01 5.000000e-03 5.000000e-03 9.080000e+02 6.246270e+00 6.246270e+00
+5.050000e-01 5.000000e-03 5.000000e-03 9.630000e+02 6.738109e+00 6.738109e+00
+5.150000e-01 5.000000e-03 5.000000e-03 1.035000e+03 6.734484e+00 6.734484e+00
+5.250000e-01 5.000000e-03 5.000000e-03 1.085000e+03 6.160315e+00 6.160315e+00
+5.350000e-01 5.000000e-03 5.000000e-03 1.158000e+03 6.088405e+00 6.088405e+00
+5.450000e-01 5.000000e-03 5.000000e-03 1.209000e+03 7.264430e+00 7.264430e+00
+5.550000e-01 5.000000e-03 5.000000e-03 1.242000e+03 6.826119e+00 6.826119e+00
+5.650000e-01 5.000000e-03 5.000000e-03 1.289000e+03 6.947035e+00 6.947035e+00
+5.750000e-01 5.000000e-03 5.000000e-03 1.276000e+03 6.297109e+00 6.297109e+00
+5.850000e-01 5.000000e-03 5.000000e-03 1.285000e+03 8.141291e+00 8.141291e+00
+5.950000e-01 5.000000e-03 5.000000e-03 1.282000e+03 7.924711e+00 7.924711e+00
+6.050000e-01 5.000000e-03 5.000000e-03 1.262000e+03 2.714612e+01 2.714612e+01
+6.150000e-01 5.000000e-03 5.000000e-03 8.981000e+02 3.625019e+01 3.625019e+01
+6.250000e-01 5.000000e-03 5.000000e-03 8.017000e+02 5.433386e+00 5.433386e+00
+6.350000e-01 5.000000e-03 5.000000e-03 7.857000e+02 3.894345e+00 3.894345e+00
+6.450000e-01 5.000000e-03 5.000000e-03 7.342000e+02 3.566847e+00 3.566847e+00
+6.550000e-01 5.000000e-03 5.000000e-03 6.838000e+02 3.656306e+00 3.656306e+00
+6.650000e-01 5.000000e-03 5.000000e-03 6.269000e+02 3.370319e+00 3.370319e+00
+6.750000e-01 5.000000e-03 5.000000e-03 5.735000e+02 3.167922e+00 3.167922e+00
+6.850000e-01 5.000000e-03 5.000000e-03 5.208000e+02 2.647614e+00 2.647614e+00
+6.950000e-01 5.000000e-03 5.000000e-03 4.760000e+02 2.635429e+00 2.635429e+00
+7.050000e-01 5.000000e-03 5.000000e-03 4.358000e+02 2.465559e+00 2.465559e+00
+7.150000e-01 5.000000e-03 5.000000e-03 3.895000e+02 2.099310e+00 2.099310e+00
+7.250000e-01 5.000000e-03 5.000000e-03 3.607000e+02 2.184440e+00 2.184440e+00
+7.350000e-01 5.000000e-03 5.000000e-03 3.311000e+02 1.829282e+00 1.829282e+00
+7.450000e-01 5.000000e-03 5.000000e-03 3.026000e+02 1.695779e+00 1.695779e+00
+7.550000e-01 5.000000e-03 5.000000e-03 2.760000e+02 1.565810e+00 1.565810e+00
+7.650000e-01 5.000000e-03 5.000000e-03 2.514000e+02 1.585696e+00 1.585696e+00
+7.750000e-01 5.000000e-03 5.000000e-03 2.302000e+02 1.319060e+00 1.319060e+00
+7.850000e-01 5.000000e-03 5.000000e-03 2.123000e+02 1.204455e+00 1.204455e+00
+7.950000e-01 5.000000e-03 5.000000e-03 1.974000e+02 1.095293e+00 1.095293e+00
+8.050000e-01 5.000000e-03 5.000000e-03 1.837000e+02 1.176298e+00 1.176298e+00
+8.150000e-01 5.000000e-03 5.000000e-03 1.713000e+02 1.067166e+00 1.067166e+00
+8.250000e-01 5.000000e-03 5.000000e-03 1.584000e+02 9.438780e-01 9.438780e-01
+8.350000e-01 5.000000e-03 5.000000e-03 1.470000e+02 1.025590e+00 1.025590e+00
+8.450000e-01 5.000000e-03 5.000000e-03 1.375000e+02 8.254733e-01 8.254733e-01
+8.550000e-01 5.000000e-03 5.000000e-03 1.274000e+02 8.843075e-01 8.843075e-01
+8.650000e-01 5.000000e-03 5.000000e-03 1.192000e+02 7.755945e-01 7.755945e-01
+8.750000e-01 5.000000e-03 5.000000e-03 1.115000e+02 6.792262e-01 6.792262e-01
+8.850000e-01 5.000000e-03 5.000000e-03 1.049000e+02 7.321914e-01 7.321914e-01
+8.950000e-01 5.000000e-03 5.000000e-03 9.870000e+01 6.428716e-01 6.428716e-01
+9.050000e-01 5.000000e-03 5.000000e-03 9.310000e+01 6.207720e-01 6.207720e-01
+9.150000e-01 5.000000e-03 5.000000e-03 8.760000e+01 6.662801e-01 6.662801e-01
+9.250000e-01 5.000000e-03 5.000000e-03 8.280000e+01 6.431688e-01 6.431688e-01
+9.350000e-01 5.000000e-03 5.000000e-03 7.874000e+01 5.816352e-01 5.816352e-01
+9.450000e-01 5.000000e-03 5.000000e-03 7.474000e+01 6.127017e-01 6.127017e-01
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/KLOE_2009_I797438/d03-x01-y01
+IsRef: 1
+Path: /REF/KLOE_2009_I797438/d03-x01-y01
+Title: doi:10.17182/hepdata.57088.v1/t3
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+3.550000e-01 5.000000e-03 5.000000e-03 7.350000e+00 1.260221e-01 1.260221e-01
+3.650000e-01 5.000000e-03 5.000000e-03 8.090000e+00 1.286486e-01 1.286486e-01
+3.750000e-01 5.000000e-03 5.000000e-03 8.680000e+00 1.277465e-01 1.277465e-01
+3.850000e-01 5.000000e-03 5.000000e-03 9.450000e+00 1.297422e-01 1.297422e-01
+3.950000e-01 5.000000e-03 5.000000e-03 1.023000e+01 1.296510e-01 1.296510e-01
+4.050000e-01 5.000000e-03 5.000000e-03 1.128000e+01 1.320758e-01 1.320758e-01
+4.150000e-01 5.000000e-03 5.000000e-03 1.230000e+01 1.375065e-01 1.375065e-01
+4.250000e-01 5.000000e-03 5.000000e-03 1.351000e+01 1.405911e-01 1.405911e-01
+4.350000e-01 5.000000e-03 5.000000e-03 1.470000e+01 1.355235e-01 1.355235e-01
+4.450000e-01 5.000000e-03 5.000000e-03 1.613000e+01 1.454806e-01 1.454806e-01
+4.550000e-01 5.000000e-03 5.000000e-03 1.791000e+01 1.464817e-01 1.464817e-01
+4.650000e-01 5.000000e-03 5.000000e-03 1.949000e+01 1.508540e-01 1.508540e-01
+4.750000e-01 5.000000e-03 5.000000e-03 2.131000e+01 1.707663e-01 1.707663e-01
+4.850000e-01 5.000000e-03 5.000000e-03 2.385000e+01 1.664561e-01 1.664561e-01
+4.950000e-01 5.000000e-03 5.000000e-03 2.661000e+01 1.778695e-01 1.778695e-01
+5.050000e-01 5.000000e-03 5.000000e-03 2.865000e+01 1.940569e-01 1.940569e-01
+5.150000e-01 5.000000e-03 5.000000e-03 3.125000e+01 2.026196e-01 2.026196e-01
+5.250000e-01 5.000000e-03 5.000000e-03 3.325000e+01 1.861752e-01 1.861752e-01
+5.350000e-01 5.000000e-03 5.000000e-03 3.605000e+01 1.931228e-01 1.931228e-01
+5.450000e-01 5.000000e-03 5.000000e-03 3.822000e+01 2.309861e-01 2.309861e-01
+5.550000e-01 5.000000e-03 5.000000e-03 3.988000e+01 2.187825e-01 2.187825e-01
+5.650000e-01 5.000000e-03 5.000000e-03 4.206000e+01 2.244251e-01 2.244251e-01
+5.750000e-01 5.000000e-03 5.000000e-03 4.227000e+01 2.041587e-01 2.041587e-01
+5.850000e-01 5.000000e-03 5.000000e-03 4.318000e+01 2.687244e-01 2.687244e-01
+5.950000e-01 5.000000e-03 5.000000e-03 4.361000e+01 2.633290e-01 2.633290e-01
+6.050000e-01 5.000000e-03 5.000000e-03 4.337000e+01 9.307995e-01 9.307995e-01
+6.150000e-01 5.000000e-03 5.000000e-03 3.303000e+01 1.333320e+00 1.333320e+00
+6.250000e-01 5.000000e-03 5.000000e-03 2.984000e+01 2.027628e-01 2.027628e-01
+6.350000e-01 5.000000e-03 5.000000e-03 2.931000e+01 1.487672e-01 1.487672e-01
+6.450000e-01 5.000000e-03 5.000000e-03 2.757000e+01 1.348363e-01 1.348363e-01
+6.550000e-01 5.000000e-03 5.000000e-03 2.590000e+01 1.368230e-01 1.368230e-01
+6.650000e-01 5.000000e-03 5.000000e-03 2.398000e+01 1.321950e-01 1.321950e-01
+6.750000e-01 5.000000e-03 5.000000e-03 2.216000e+01 1.203489e-01 1.203489e-01
+6.850000e-01 5.000000e-03 5.000000e-03 2.033000e+01 1.068011e-01 1.068011e-01
+6.950000e-01 5.000000e-03 5.000000e-03 1.878000e+01 1.048091e-01 1.048091e-01
+7.050000e-01 5.000000e-03 5.000000e-03 1.738000e+01 1.016210e-01 1.016210e-01
+7.150000e-01 5.000000e-03 5.000000e-03 1.570000e+01 8.581900e-02 8.581900e-02
+7.250000e-01 5.000000e-03 5.000000e-03 1.469000e+01 9.256637e-02 9.256637e-02
+7.350000e-01 5.000000e-03 5.000000e-03 1.363000e+01 7.387671e-02 7.387671e-02
+7.450000e-01 5.000000e-03 5.000000e-03 1.260000e+01 7.202500e-02 7.202500e-02
+7.550000e-01 5.000000e-03 5.000000e-03 1.163000e+01 6.206907e-02 6.206907e-02
+7.650000e-01 5.000000e-03 5.000000e-03 1.070000e+01 6.668081e-02 6.668081e-02
+7.750000e-01 5.000000e-03 5.000000e-03 9.910000e+00 5.900916e-02 5.900916e-02
+7.850000e-01 5.000000e-03 5.000000e-03 9.240000e+00 4.953560e-02 4.953560e-02
+7.950000e-01 5.000000e-03 5.000000e-03 8.680000e+00 4.851210e-02 4.851210e-02
+8.050000e-01 5.000000e-03 5.000000e-03 8.160000e+00 5.226811e-02 5.226811e-02
+8.150000e-01 5.000000e-03 5.000000e-03 7.690000e+00 5.104227e-02 5.104227e-02
+8.250000e-01 5.000000e-03 5.000000e-03 7.180000e+00 4.171959e-02 4.171959e-02
+8.350000e-01 5.000000e-03 5.000000e-03 6.732000e+00 4.692883e-02 4.692883e-02
+8.450000e-01 5.000000e-03 5.000000e-03 6.358000e+00 3.983981e-02 3.983981e-02
+8.550000e-01 5.000000e-03 5.000000e-03 5.948000e+00 4.127767e-02 4.127767e-02
+8.650000e-01 5.000000e-03 5.000000e-03 5.621000e+00 3.482996e-02 3.482996e-02
+8.750000e-01 5.000000e-03 5.000000e-03 5.304000e+00 3.321522e-02 3.321522e-02
+8.850000e-01 5.000000e-03 5.000000e-03 5.038000e+00 3.448068e-02 3.448068e-02
+8.950000e-01 5.000000e-03 5.000000e-03 4.784000e+00 3.284894e-02 3.284894e-02
+9.050000e-01 5.000000e-03 5.000000e-03 4.550000e+00 3.063111e-02 3.063111e-02
+9.150000e-01 5.000000e-03 5.000000e-03 4.322000e+00 3.243684e-02 3.243684e-02
+9.250000e-01 5.000000e-03 5.000000e-03 4.117000e+00 3.084053e-02 3.084053e-02
+9.350000e-01 5.000000e-03 5.000000e-03 3.950000e+00 2.943285e-02 2.943285e-02
+9.450000e-01 5.000000e-03 5.000000e-03 3.780000e+00 3.115156e-02 3.115156e-02
+END YODA_SCATTER2D_V2
diff --git a/analyses/pluginFrascati/MUPI_1972_I84978.cc b/analyses/pluginFrascati/MUPI_1972_I84978.cc
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/MUPI_1972_I84978.cc
@@ -0,0 +1,99 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/Projections/FinalState.hh"
+
+namespace Rivet {
+
+
+ /// @brief Add a short analysis description here
+ class MUPI_1972_I84978 : public Analysis {
+ public:
+
+ /// Constructor
+ DEFAULT_RIVET_ANALYSIS_CTOR(MUPI_1972_I84978);
+
+
+ /// @name Analysis methods
+ //@{
+
+ /// Book histograms and initialise projections before the run
+ void init() {
+ // Initialise and register projections
+ declare(FinalState(), "FS");
+
+ // Book histograms
+ _c_hadrons = bookCounter("/TMP/sigma_hadrons");
+ _c_muons = bookCounter("/TMP/sigma_muons");
+ }
+
+
+ /// Perform the per-event analysis
+ void analyze(const Event& event) {
+ const FinalState& fs = apply<FinalState>(event, "FS");
+
+ map<long,int> nCount;
+ int ntotal(0);
+ foreach (const Particle& p, fs.particles()) {
+ nCount[p.pdgId()] += 1;
+ ++ntotal;
+ }
+ // mu+mu- + photons
+ if(nCount[-13]==1 and nCount[13]==1 &&
+ ntotal==2+nCount[22])
+ _c_muons->fill(event.weight());
+ // everything else
+ else
+ _c_hadrons->fill(event.weight());
+ }
+
+
+ /// Normalise histograms etc., after the run
+ void finalize() {
+ Scatter1D R = *_c_hadrons/ *_c_muons;
+ double rval = R.point(0).x();
+ pair<double,double> rerr = R.point(0).xErrs();
+ double fact = crossSection()/ sumOfWeights() /nanobarn;
+ double sig_h = _c_hadrons->val()*fact;
+ double err_h = _c_hadrons->err()*fact;
+ double sig_m = _c_muons ->val()*fact;
+ double err_m = _c_muons ->err()*fact;
+ Scatter2D temphisto(refData(1, 1, 1));
+ Scatter2DPtr hadrons = bookScatter2D(1, 1, 1);
+ Scatter2DPtr muons = bookScatter2D("sigma_muons" );
+ Scatter2DPtr mult = bookScatter2D(2, 1, 1);
+ for (size_t b = 0; b < temphisto.numPoints(); b++) {
+ const double x = temphisto.point(b).x();
+ pair<double,double> ex = temphisto.point(b).xErrs();
+ pair<double,double> ex2 = ex;
+ if(ex2.first ==0.) ex2. first=0.0001;
+ if(ex2.second==0.) ex2.second=0.0001;
+ if (inRange(sqrtS()/GeV, x-ex2.first, x+ex2.second)) {
+ mult ->addPoint(x, rval, ex, rerr);
+ hadrons->addPoint(x, sig_h, ex, make_pair(err_h,err_h));
+ muons ->addPoint(x, sig_m, ex, make_pair(err_m,err_m));
+ }
+ else {
+ mult ->addPoint(x, 0., ex, make_pair(0.,.0));
+ hadrons->addPoint(x, 0., ex, make_pair(0.,.0));
+ muons ->addPoint(x, 0., ex, make_pair(0.,.0));
+ }
+ }
+ }
+
+ //@}
+
+
+ /// @name Histograms
+ //@{
+ CounterPtr _c_hadrons, _c_muons;
+ //@}
+
+
+ };
+
+
+ // The hook for the plugin system
+ DECLARE_RIVET_PLUGIN(MUPI_1972_I84978);
+
+
+}
diff --git a/analyses/pluginFrascati/MUPI_1972_I84978.info b/analyses/pluginFrascati/MUPI_1972_I84978.info
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/MUPI_1972_I84978.info
@@ -0,0 +1,36 @@
+Name: MUPI_1972_I84978
+Year: 1972
+Summary: Measurement of $R$ and the hadronic cross section for energies between 1.2 and 2.1 GeV
+Experiment: MUPI
+Collider: ADONE
+InspireID: 84978
+Status: VALIDATED
+Authors:
+ - Peter Richardson <peter.richardson@durham.ac.uk>
+References:
+ - Phys.Lett. B218 (1989) 499-507, 1989
+RunInfo: e+ e- to hadrons and e+ e- to mu+ mu- (for normalization)
+NeedCrossSection: yes
+Beams: [e-, e+]
+Description:
+ 'Measurement of $R$ and the hadronic cross seection in $e^+e^-$ collisions by
+ the Mu-Pi group for energies between 1.2 and 2.1 GeV.
+ The muonic cross section is also outputted to the yoda file
+ so that ratio $R$ can be recalcuated if runs are combined.'
+Keywords: []
+InspireID: 84978
+BibKey: Ceradini:1973pr
+BibTeX: '@article{Ceradini:1973pr,
+ author = "Ceradini, F. and Santonico, R. and Conversi, M. and
+ Paoluzi, L. and Grilli, M. and Iarocci, E. and
+ Spillantini, P. and Valente, V. and Visentin, R. and
+ Nigro, M.",
+ title = "{Multiplicity in hadron production by e+ e- colliding
+ beams}",
+ journal = "Phys. Lett.",
+ volume = "42B",
+ year = "1972",
+ pages = "501-503",
+ doi = "10.1016/0370-2693(72)90116-5",
+ SLACcitation = "%%CITATION = PHLTA,42B,501;%%"
+}'
diff --git a/analyses/pluginFrascati/MUPI_1972_I84978.plot b/analyses/pluginFrascati/MUPI_1972_I84978.plot
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/MUPI_1972_I84978.plot
@@ -0,0 +1,21 @@
+BEGIN PLOT /MUPI_1972_I84978/d02-x01-y01
+Title=$R=\sigma(e^+e^-\to \text{hadrons})/\sigma(e^+e^-\to \mu^+\mu^-)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$R$
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /MUPI_1972_I84978/d01-x01-y01
+Title=$\sigma(e^+e^-\to \text{hadrons})$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to \text{hadrons})/pb$
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /MUPI_1972_I84978/sigma_muons
+Title=$\sigma(e^+e^-\to \mu^+\mu^-)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to \mu^+\mu^-)/pb$
+LogY=0
+ConnectGaps=1
+END PLOT
diff --git a/analyses/pluginFrascati/MUPI_1972_I84978.yoda b/analyses/pluginFrascati/MUPI_1972_I84978.yoda
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/MUPI_1972_I84978.yoda
@@ -0,0 +1,24 @@
+BEGIN YODA_SCATTER2D_V2 /REF/MUPI_1972_I84978/d01-x01-y01
+IsRef: 1
+Path: /REF/MUPI_1972_I84978/d01-x01-y01
+Title: doi:10.17182/hepdata.28188.v1/t1
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.200000e+00 0.000000e+00 0.000000e+00 3.400000e+01 1.612452e+01 1.612452e+01
+1.500000e+00 0.000000e+00 0.000000e+00 6.600000e+01 9.000000e+00 9.000000e+00
+1.900000e+00 0.000000e+00 0.000000e+00 3.600000e+01 5.000000e+00 5.000000e+00
+2.100000e+00 0.000000e+00 0.000000e+00 3.600000e+01 4.000000e+00 4.000000e+00
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/MUPI_1972_I84978/d02-x01-y01
+IsRef: 1
+Path: /REF/MUPI_1972_I84978/d02-x01-y01
+Title: doi:10.17182/hepdata.28188.v1/t2
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+1.200000e+00 0.000000e+00 0.000000e+00 5.600000e-01 2.716616e-01 2.716616e-01
+1.500000e+00 0.000000e+00 0.000000e+00 1.710000e+00 2.300000e-01 2.300000e-01
+1.900000e+00 0.000000e+00 0.000000e+00 1.500000e+00 2.100000e-01 2.100000e-01
+2.100000e+00 0.000000e+00 0.000000e+00 1.830000e+00 2.000000e-01 2.000000e-01
+END YODA_SCATTER2D_V2
diff --git a/analyses/pluginFrascati/MUPI_1973_I95215.cc b/analyses/pluginFrascati/MUPI_1973_I95215.cc
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/MUPI_1973_I95215.cc
@@ -0,0 +1,99 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/Projections/FinalState.hh"
+
+namespace Rivet {
+
+
+ /// @brief Add a short analysis description here
+ class MUPI_1973_I95215 : public Analysis {
+ public:
+
+ /// Constructor
+ DEFAULT_RIVET_ANALYSIS_CTOR(MUPI_1973_I95215);
+
+
+ /// @name Analysis methods
+ //@{
+
+ /// Book histograms and initialise projections before the run
+ void init() {
+ // Initialise and register projections
+ declare(FinalState(), "FS");
+
+ // Book histograms
+ _c_hadrons = bookCounter("/TMP/sigma_hadrons");
+ _c_muons = bookCounter("/TMP/sigma_muons");
+ }
+
+
+ /// Perform the per-event analysis
+ void analyze(const Event& event) {
+ const FinalState& fs = apply<FinalState>(event, "FS");
+
+ map<long,int> nCount;
+ int ntotal(0);
+ foreach (const Particle& p, fs.particles()) {
+ nCount[p.pdgId()] += 1;
+ ++ntotal;
+ }
+ // mu+mu- + photons
+ if(nCount[-13]==1 and nCount[13]==1 &&
+ ntotal==2+nCount[22])
+ _c_muons->fill(event.weight());
+ // everything else
+ else
+ _c_hadrons->fill(event.weight());
+ }
+
+
+ /// Normalise histograms etc., after the run
+ void finalize() {
+ Scatter1D R = *_c_hadrons/ *_c_muons;
+ double rval = R.point(0).x();
+ pair<double,double> rerr = R.point(0).xErrs();
+ double fact = crossSection()/ sumOfWeights() /nanobarn;
+ double sig_h = _c_hadrons->val()*fact;
+ double err_h = _c_hadrons->err()*fact;
+ double sig_m = _c_muons ->val()*fact;
+ double err_m = _c_muons ->err()*fact;
+ Scatter2D temphisto(refData(1, 1, 1));
+ Scatter2DPtr hadrons = bookScatter2D(1, 1, 1);
+ Scatter2DPtr muons = bookScatter2D("sigma_muons" );
+ Scatter2DPtr mult = bookScatter2D(1, 1, 2);
+ for (size_t b = 0; b < temphisto.numPoints(); b++) {
+ const double x = temphisto.point(b).x();
+ pair<double,double> ex = temphisto.point(b).xErrs();
+ pair<double,double> ex2 = ex;
+ if(ex2.first ==0.) ex2. first=0.0001;
+ if(ex2.second==0.) ex2.second=0.0001;
+ if (inRange(sqrtS()/GeV, x-ex2.first, x+ex2.second)) {
+ mult ->addPoint(x, rval, ex, rerr);
+ hadrons->addPoint(x, sig_h, ex, make_pair(err_h,err_h));
+ muons ->addPoint(x, sig_m, ex, make_pair(err_m,err_m));
+ }
+ else {
+ mult ->addPoint(x, 0., ex, make_pair(0.,.0));
+ hadrons->addPoint(x, 0., ex, make_pair(0.,.0));
+ muons ->addPoint(x, 0., ex, make_pair(0.,.0));
+ }
+ }
+ }
+
+ //@}
+
+
+ /// @name Histograms
+ //@{
+ CounterPtr _c_hadrons, _c_muons;
+ //@}
+
+
+ };
+
+
+ // The hook for the plugin system
+ DECLARE_RIVET_PLUGIN(MUPI_1973_I95215);
+
+
+}
diff --git a/analyses/pluginFrascati/MUPI_1973_I95215.info b/analyses/pluginFrascati/MUPI_1973_I95215.info
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/MUPI_1973_I95215.info
@@ -0,0 +1,33 @@
+Name: MUPI_1973_I95215
+Year: 1973
+Summary: Measurement of $R$ and the hadronic cross section for energies betwee 2.6 and 3 GeV
+Experiment: MUPI
+Collider: ADONE
+InspireID: 95215
+Status: VALIDATED
+Authors:
+ - Peter Richardson <peter.richardson@durham.ac.uk>
+References:
+ - Phys.Lett. B47 (1973) 80-84, 1973
+RunInfo: e+ e- to hadrons and e+ e- to mu+ mu- (for normalization)
+NeedCrossSection: yes
+Beams: [e-, e+]
+Description:
+ 'Measurement of $R$ and the hadronic cross seection in $e^+e^-$ collisions by
+ the Mu-Pi group for energies between 2.6 and 3 GeV.
+ The muonic cross section is also outputted to the yoda file
+ so that ratio $R$ can be recalcuated if runs are combined.'
+Keywords: []
+BibKey: Ceradini:1974nt
+BibTeX: '@article{Ceradini:1974nt,
+ author = "Ceradini, F. and Conversi, M. and D Angelo, S. and
+ Paoluzi, L. and Santonico, R. and Visentin, R.",
+ title = "{Multihadron production in e+ e- collisions up to 3 gev
+ total center-of-mass energy}",
+ journal = "Phys. Lett.",
+ volume = "47B",
+ year = "1973",
+ pages = "80-84",
+ doi = "10.1016/0370-2693(73)90574-1",
+ SLACcitation = "%%CITATION = PHLTA,47B,80;%%"
+}'
diff --git a/analyses/pluginFrascati/MUPI_1973_I95215.plot b/analyses/pluginFrascati/MUPI_1973_I95215.plot
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/MUPI_1973_I95215.plot
@@ -0,0 +1,21 @@
+BEGIN PLOT /MUPI_1973_I95215/d01-x01-y02
+Title=$R=\sigma(e^+e^-\to \text{hadrons})/\sigma(e^+e^-\to \mu^+\mu^-)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$R$
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /MUPI_1973_I95215/d01-x01-y01
+Title=$\sigma(e^+e^-\to \text{hadrons})$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to \text{hadrons})/$nb
+LogY=0
+ConnectGaps=1
+END PLOT
+BEGIN PLOT /MUPI_1973_I95215/sigma_muons
+Title=$\sigma(e^+e^-\to \mu^+\mu^-)$
+XLabel=$\sqrt{s}$/GeV
+YLabel=$\sigma(e^+e^-\to \mu^+\mu^-)/$nb
+LogY=0
+ConnectGaps=1
+END PLOT
diff --git a/analyses/pluginFrascati/MUPI_1973_I95215.yoda b/analyses/pluginFrascati/MUPI_1973_I95215.yoda
new file mode 100644
--- /dev/null
+++ b/analyses/pluginFrascati/MUPI_1973_I95215.yoda
@@ -0,0 +1,22 @@
+BEGIN YODA_SCATTER2D_V2 /REF/MUPI_1973_I95215/d01-x01-y01
+IsRef: 1
+Path: /REF/MUPI_1973_I95215/d01-x01-y01
+Title: doi:10.17182/hepdata.28037.v1/t1
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+2.600000e+00 0.000000e+00 0.000000e+00 1.830000e+01 4.700000e+00 4.700000e+00
+2.800000e+00 0.000000e+00 0.000000e+00 1.700000e+01 5.000000e+00 5.000000e+00
+3.000000e+00 0.000000e+00 0.000000e+00 1.390000e+01 5.000000e+00 5.000000e+00
+END YODA_SCATTER2D_V2
+BEGIN YODA_SCATTER2D_V2 /REF/MUPI_1973_I95215/d01-x01-y02
+IsRef: 1
+Path: /REF/MUPI_1973_I95215/d01-x01-y02
+Title: doi:10.17182/hepdata.28037.v1/t1
+Type: Scatter2D
+---
+# xval xerr- xerr+ yval yerr- yerr+
+2.600000e+00 0.000000e+00 0.000000e+00 1.470000e+00 3.800000e-01 3.800000e-01
+2.800000e+00 0.000000e+00 0.000000e+00 1.580000e+00 4.700000e-01 4.700000e-01
+3.000000e+00 0.000000e+00 0.000000e+00 1.490000e+00 5.400000e-01 5.400000e-01
+END YODA_SCATTER2D_V2
File Metadata
Details
Attached
Mime Type
text/x-diff
Expires
Tue, Nov 19, 8:04 PM (1 d, 7 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
3805966
Default Alt Text
(142 KB)
Attached To
rRIVETHG rivethg
Event Timeline
Log In to Comment