diff --git a/analyses/pluginNovosibirsk/SND_2001_I579319.cc b/analyses/pluginNovosibirsk/SND_2001_I579319.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2001_I579319.cc @@ -0,0 +1,97 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/Projections/FinalState.hh" + +namespace Rivet { + + + /// @brief Add a short analysis description here + class SND_2001_I579319 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(SND_2001_I579319); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + + // Initialise and register projections + declare(FinalState(), "FS"); + + // Book histograms + _npion1 = bookCounter("TMP/pion1"); + _npion2 = bookCounter("TMP/pion2"); + } + + + /// Perform the per-event analysis + void analyze(const Event& event) { + const FinalState& fs = apply(event, "FS"); + + map nCount; + int ntotal(0); + foreach (const Particle& p, fs.particles()) { + nCount[p.pdgId()] += 1; + ++ntotal; + } + if(ntotal!=4) vetoEvent; + if(nCount[-211]==2&&nCount[211]==2) + _npion1->fill(event.weight()); + else if(nCount[-211]==1&&nCount[211]==1&&nCount[111]==2) + _npion2->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 = _npion1->val(); + error = _npion1->err(); + } + else { + sigma = _npion2->val(); + error = _npion2->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 ex = temphisto.point(b).xErrs(); + pair 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 _npion1,_npion2; + //@} + + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(SND_2001_I579319); + + +} diff --git a/analyses/pluginNovosibirsk/SND_2001_I579319.info b/analyses/pluginNovosibirsk/SND_2001_I579319.info new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2001_I579319.info @@ -0,0 +1,27 @@ +Name: SND_2001_I579319 +Year: 2001 +Summary: Cross section for $e^+e^-2\pi^+2\pi^-$ and $\pi^+\pi^-2\pi^0$ between 1 and 1.4 GeV +Experiment: SND +Collider: VEPP-2M +InspireID: 579319 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - INP-2001-34 (2001) +RunInfo: e+e- to hadrons +NeedCrossSection: yes +Beams: [e+,e-] +Luminosity_fb: +Description: + 'Measurement of the cross section for $e^+e^-\to\pi^+2\pi^-$ and $\pi^+\pi^-2\pi^0$ between 1 and 1.4 GeV' +BibKey: Achasov:2001gp +BibTeX: '@article{Achasov:2001gp, + author = "Achasov, M. N. and others", + title = "{e+ e- --> 4pi processes investigation in the energy + range 0.98-GeV to 1.38-GeV with SND detector}", + collaboration = "SND", + year = "2001", + reportNumber = "BUDKER-INP-2001-34", + SLACcitation = "%%CITATION = BUDKER-INP-2001-34;%%" +}' diff --git a/analyses/pluginNovosibirsk/SND_2001_I579319.plot b/analyses/pluginNovosibirsk/SND_2001_I579319.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2001_I579319.plot @@ -0,0 +1,16 @@ +BEGIN PLOT /SND_2001_I579319/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 +BEGIN PLOT /SND_2001_I579319/d02-x01-y01 +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 +LegendXPos=0.8 +ConnectGaps=1 +END PLOT diff --git a/analyses/pluginNovosibirsk/SND_2001_I579319.yoda b/analyses/pluginNovosibirsk/SND_2001_I579319.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2001_I579319.yoda @@ -0,0 +1,99 @@ +BEGIN YODA_SCATTER2D_V2 /REF/SND_2001_I579319/d01-x01-y01 +IsRef: 1 +Path: /REF/SND_2001_I579319/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +0.98000 0.00000 0.00000 1.07000e+00 2.10000e-01 2.10000e-01 +0.98410 0.00000 0.00000 1.21000e+00 9.00000e-02 9.00000e-02 +1.00380 0.00000 0.00000 1.26000e+00 1.20000e-01 1.20000e-01 +1.00968 0.00000 0.00000 1.43000e+00 1.50000e-01 1.50000e-01 +1.01566 0.00000 0.00000 1.68000e+00 1.60000e-01 1.60000e-01 +1.01668 0.00000 0.00000 1.60000e+00 1.20000e-01 1.20000e-01 +1.01766 0.00000 0.00000 1.78000e+00 1.20000e-01 1.20000e-01 +1.01864 0.00000 0.00000 1.72000e+00 1.10000e-01 1.10000e-01 +1.01962 0.00000 0.00000 1.96000e+00 1.00000e-01 1.00000e-01 +1.02058 0.00000 0.00000 1.76000e+00 1.20000e-01 1.20000e-01 +1.02164 0.00000 0.00000 1.92000e+00 1.70000e-01 1.70000e-01 +1.02276 0.00000 0.00000 1.72000e+00 1.50000e-01 1.50000e-01 +1.02774 0.00000 0.00000 1.91000e+00 1.50000e-01 1.50000e-01 +1.03372 0.00000 0.00000 1.88000e+00 1.70000e-01 1.70000e-01 +1.04000 0.00000 0.00000 2.09000e+00 1.50000e-01 1.50000e-01 +1.05000 0.00000 0.00000 2.37000e+00 1.50000e-01 1.50000e-01 +1.06000 0.00000 0.00000 2.52000e+00 1.30000e-01 1.30000e-01 +1.07000 0.00000 0.00000 3.11000e+00 3.50000e-01 3.50000e-01 +1.08000 0.00000 0.00000 4.67000e+00 1.80000e-01 1.80000e-01 +1.09000 0.00000 0.00000 3.52000e+00 3.80000e-01 3.80000e-01 +1.10000 0.00000 0.00000 3.92000e+00 1.90000e-01 1.90000e-01 +1.11000 0.00000 0.00000 4.48000e+00 4.30000e-01 4.30000e-01 +1.12000 0.00000 0.00000 4.93000e+00 2.50000e-01 2.50000e-01 +1.13000 0.00000 0.00000 5.45000e+00 4.30000e-01 4.30000e-01 +1.14000 0.00000 0.00000 6.00000e+00 2.80000e-01 2.80000e-01 +1.15000 0.00000 0.00000 7.32000e+00 6.20000e-01 6.20000e-01 +1.16000 0.00000 0.00000 7.33000e+00 3.00000e-01 3.00000e-01 +1.18000 0.00000 0.00000 8.60000e+00 2.80000e-01 2.80000e-01 +1.19000 0.00000 0.00000 8.68000e+00 4.20000e-01 4.20000e-01 +1.20000 0.00000 0.00000 1.15700e+01 3.20000e-01 3.20000e-01 +1.21000 0.00000 0.00000 1.05800e+01 5.10000e-01 5.10000e-01 +1.22000 0.00000 0.00000 1.16600e+01 3.60000e-01 3.60000e-01 +1.23000 0.00000 0.00000 1.08400e+01 5.20000e-01 5.20000e-01 +1.24000 0.00000 0.00000 1.33600e+01 3.60000e-01 3.60000e-01 +1.25000 0.00000 0.00000 1.50300e+01 5.00000e-01 5.00000e-01 +1.26000 0.00000 0.00000 1.36900e+01 5.40000e-01 5.40000e-01 +1.27000 0.00000 0.00000 1.67400e+01 4.90000e-01 4.90000e-01 +1.28000 0.00000 0.00000 1.66600e+01 5.00000e-01 5.00000e-01 +1.29000 0.00000 0.00000 1.62100e+01 4.50000e-01 4.50000e-01 +1.30000 0.00000 0.00000 1.83200e+01 5.00000e-01 5.00000e-01 +1.31000 0.00000 0.00000 1.73800e+01 5.30000e-01 5.30000e-01 +1.32000 0.00000 0.00000 1.72800e+01 4.90000e-01 4.90000e-01 +1.33000 0.00000 0.00000 1.85500e+01 4.50000e-01 4.50000e-01 +1.34000 0.00000 0.00000 2.21900e+01 4.60000e-01 4.60000e-01 +1.35000 0.00000 0.00000 2.16600e+01 5.30000e-01 5.30000e-01 +1.36000 0.00000 0.00000 2.39900e+01 3.60000e-01 3.60000e-01 +1.37000 0.00000 0.00000 2.38600e+01 5.50000e-01 5.50000e-01 +1.38000 0.00000 0.00000 2.40900e+01 4.00000e-01 4.00000e-01 +END YODA_SCATTER2D_V2 +BEGIN YODA_SCATTER2D_V2 /REF/SND_2001_I579319/d02-x01-y01 +IsRef: 1 +Path: /REF/SND_2001_I579319/d02-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +0.98000 0.00000 0.00000 5.56000e+00 1.58000e+00 1.58000e+00 +1.01900 0.00000 0.00000 8.07000e+00 1.65000e+00 1.65000e+00 +1.04000 0.00000 0.00000 9.17000e+00 1.47000e+00 1.47000e+00 +1.05000 0.00000 0.00000 1.06600e+01 1.35000e+00 1.35000e+00 +1.06000 0.00000 0.00000 1.14500e+01 7.50000e-01 7.50000e-01 +1.07000 0.00000 0.00000 1.02600e+01 1.25000e+00 1.25000e+00 +1.08000 0.00000 0.00000 1.36700e+01 6.30000e-01 6.30000e-01 +1.09000 0.00000 0.00000 1.30100e+01 1.37000e+00 1.37000e+00 +1.10000 0.00000 0.00000 1.44100e+01 4.20000e-01 4.20000e-01 +1.11000 0.00000 0.00000 1.20200e+01 8.60000e-01 8.60000e-01 +1.12000 0.00000 0.00000 1.56700e+01 5.40000e-01 5.40000e-01 +1.13000 0.00000 0.00000 1.54100e+01 8.80000e-01 8.80000e-01 +1.14000 0.00000 0.00000 1.65000e+01 5.80000e-01 5.80000e-01 +1.15000 0.00000 0.00000 1.57800e+01 1.16000e+00 1.16000e+00 +1.16000 0.00000 0.00000 1.73400e+01 5.80000e-01 5.80000e-01 +1.18000 0.00000 0.00000 1.91400e+01 5.40000e-01 5.40000e-01 +1.20000 0.00000 0.00000 1.95400e+01 5.50000e-01 5.50000e-01 +1.21000 0.00000 0.00000 1.93600e+01 9.20000e-01 9.20000e-01 +1.22000 0.00000 0.00000 2.01000e+01 6.40000e-01 6.40000e-01 +1.23000 0.00000 0.00000 2.03400e+01 1.00000e+00 1.00000e+00 +1.24000 0.00000 0.00000 2.20300e+01 6.30000e-01 6.30000e-01 +1.25000 0.00000 0.00000 2.23500e+01 8.50000e-01 8.50000e-01 +1.26000 0.00000 0.00000 2.23300e+01 9.40000e-01 9.40000e-01 +1.27000 0.00000 0.00000 2.46900e+01 8.20000e-01 8.20000e-01 +1.28000 0.00000 0.00000 2.40300e+01 8.30000e-01 8.30000e-01 +1.29000 0.00000 0.00000 2.44600e+01 7.60000e-01 7.60000e-01 +1.30000 0.00000 0.00000 2.65700e+01 8.30000e-01 8.30000e-01 +1.31000 0.00000 0.00000 2.42900e+01 8.70000e-01 8.70000e-01 +1.32000 0.00000 0.00000 2.59200e+01 8.40000e-01 8.40000e-01 +1.33000 0.00000 0.00000 2.63100e+01 7.50000e-01 7.50000e-01 +1.34000 0.00000 0.00000 2.71900e+01 6.30000e-01 6.30000e-01 +1.35000 0.00000 0.00000 2.88100e+01 8.30000e-01 8.30000e-01 +1.36000 0.00000 0.00000 2.83500e+01 5.40000e-01 5.40000e-01 +1.37000 0.00000 0.00000 2.89900e+01 8.10000e-01 8.10000e-01 +1.38000 0.00000 0.00000 2.88000e+01 5.80000e-01 5.80000e-01 +END YODA_SCATTER2D_V2 diff --git a/analyses/pluginNovosibirsk/SND_2015_I1389908.cc b/analyses/pluginNovosibirsk/SND_2015_I1389908.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2015_I1389908.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 SND_2015_I1389908 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(SND_2015_I1389908); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + // Initialise and register projections + declare(FinalState(), "FS"); + _num3pi = bookCounter("TMP/num3"); + } + + + /// Perform the per-event analysis + void analyze(const Event& event) { + const FinalState& fs = apply(event, "FS"); + + map nCount; + int ntotal(0); + foreach (const Particle& p, fs.particles()) { + nCount[p.pdgId()] += 1; + ++ntotal; + } + if(ntotal!=3) vetoEvent; + if(nCount[-211]==1&&nCount[211]==1&&nCount[111]==1) + _num3pi->fill(event.weight()); + } + + + /// Normalise histograms etc., after the run + void finalize() { + double sigma = _num3pi->val(); + double error = _num3pi->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 ex = temphisto.point(b).xErrs(); + pair 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 _num3pi; + //@} + + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(SND_2015_I1389908); + + +} diff --git a/analyses/pluginNovosibirsk/SND_2015_I1389908.info b/analyses/pluginNovosibirsk/SND_2015_I1389908.info new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2015_I1389908.info @@ -0,0 +1,37 @@ +Name: SND_2015_I1389908 +Year: 2015 +Summary: Cross section for $e^+e^-\to \pi^+\pi^-\pi^0$ between 1.05 and 2 GeV. +Experiment: SND +Collider: VEPP-2M +InspireID: 1389908 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - J.Exp.Theor.Phys. 121 (2015) no.1, 27-34 +RunInfo: e+ e- to hadrons +NeedCrossSection: yes +Beams: [e+, e-] +Description: + 'Measurement of the cross section for $e^+e^-\to \pi^+\pi^-\pi^0$ by SND for energies between 1.05 and 2 GeV.' +BibKey: Aulchenko:2015mwt +BibTeX: '%%% contains utf-8, see: http://inspirehep.net/info/faq/general#utf8 +%%% add \usepackage[utf8]{inputenc} to your latex preamble + +@article{Aulchenko:2015mwt, + author = "Aulchenko, V. M. and others", + title = "{Study of the e$^{+}$ e$^{–}$ → + $\pi^+\pi^-\pi^0$ process in the energy range + 1.05–2.00 GeV}", + journal = "J. Exp. Theor. Phys.", + volume = "121", + year = "2015", + number = "1", + pages = "27-34", + doi = "10.1134/S1063776115060023", + note = "[Zh. Eksp. Teor. Fiz.148,no.1,34(2015)]", + SLACcitation = "%%CITATION = JTPHE,121,27;%%" +}' +ToDo: + - Implement the analysis, test it, remove this ToDo, and mark as VALIDATED :-) + diff --git a/analyses/pluginNovosibirsk/SND_2015_I1389908.plot b/analyses/pluginNovosibirsk/SND_2015_I1389908.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2015_I1389908.plot @@ -0,0 +1,7 @@ +BEGIN PLOT /SND_2015_I1389908/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 +ConnectGaps=1 +LegendXPos=0.1 +END PLOT diff --git a/analyses/pluginNovosibirsk/SND_2015_I1389908.yoda b/analyses/pluginNovosibirsk/SND_2015_I1389908.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2015_I1389908.yoda @@ -0,0 +1,48 @@ +BEGIN YODA_SCATTER2D_V2 /REF/SND_2015_I1389908/d01-x01-y01 +IsRef: 1 +Path: /REF/SND_2015_I1389908/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.05000 0.00000 0.00000 9.80000e-01 5.88218e-01 5.88218e-01 +1.07500 0.00000 0.00000 3.07000e+00 7.22357e-01 7.22357e-01 +1.10000 0.00000 0.00000 3.96000e+00 6.57647e-01 6.57647e-01 +1.12500 0.00000 0.00000 4.52000e+00 3.75366e-01 3.75366e-01 +1.15000 0.00000 0.00000 5.13000e+00 3.82099e-01 3.82099e-01 +1.17500 0.00000 0.00000 5.29000e+00 3.54683e-01 3.54683e-01 +1.20000 0.00000 0.00000 5.07000e+00 3.48281e-01 3.48281e-01 +1.22500 0.00000 0.00000 5.72000e+00 3.67967e-01 3.67967e-01 +1.25000 0.00000 0.00000 5.93000e+00 3.67967e-01 3.67967e-01 +1.27500 0.00000 0.00000 5.50000e+00 3.76431e-01 3.76431e-01 +1.30000 0.00000 0.00000 4.87000e+00 3.34215e-01 3.34215e-01 +1.32500 0.00000 0.00000 4.86000e+00 3.04138e-01 3.04138e-01 +1.35000 0.00000 0.00000 4.97000e+00 3.18904e-01 3.18904e-01 +1.37500 0.00000 0.00000 4.75000e+00 2.97321e-01 2.97321e-01 +1.40000 0.00000 0.00000 4.12000e+00 3.00000e-01 3.00000e-01 +1.42500 0.00000 0.00000 3.99000e+00 2.86007e-01 2.86007e-01 +1.45000 0.00000 0.00000 4.04000e+00 2.94109e-01 2.94109e-01 +1.47500 0.00000 0.00000 4.25000e+00 2.76586e-01 2.76586e-01 +1.50000 0.00000 0.00000 4.39000e+00 2.68701e-01 2.68701e-01 +1.52500 0.00000 0.00000 4.47000e+00 3.06105e-01 3.06105e-01 +1.55000 0.00000 0.00000 4.58000e+00 3.12410e-01 3.12410e-01 +1.57500 0.00000 0.00000 4.66000e+00 3.12410e-01 3.12410e-01 +1.60000 0.00000 0.00000 5.73000e+00 3.67967e-01 3.67967e-01 +1.62500 0.00000 0.00000 5.03000e+00 3.56090e-01 3.56090e-01 +1.65000 0.00000 0.00000 4.64000e+00 3.20156e-01 3.20156e-01 +1.67500 0.00000 0.00000 3.45000e+00 2.66271e-01 2.66271e-01 +1.70000 0.00000 0.00000 2.65000e+00 2.54951e-01 2.54951e-01 +1.72500 0.00000 0.00000 2.18000e+00 2.10238e-01 2.10238e-01 +1.75000 0.00000 0.00000 1.85000e+00 1.96977e-01 1.96977e-01 +1.77500 0.00000 0.00000 1.66000e+00 1.74642e-01 1.74642e-01 +1.80000 0.00000 0.00000 1.08000e+00 1.86815e-01 1.86815e-01 +1.82500 0.00000 0.00000 1.31000e+00 1.52315e-01 1.52315e-01 +1.85000 0.00000 0.00000 1.30000e+00 1.80278e-01 1.80278e-01 +1.87000 0.00000 0.00000 9.30000e-01 1.36015e-01 1.36015e-01 +1.89000 0.00000 0.00000 6.80000e-01 1.23693e-01 1.23693e-01 +1.90000 0.00000 0.00000 1.03000e+00 1.45602e-01 1.45602e-01 +1.92500 0.00000 0.00000 6.50000e-01 1.04403e-01 1.04403e-01 +1.95000 0.00000 0.00000 4.90000e-01 1.21655e-01 1.21655e-01 +1.97500 0.00000 0.00000 6.60000e-01 1.33417e-01 1.33417e-01 +2.00000 0.00000 0.00000 8.10000e-01 1.52971e-01 1.52971e-01 +END YODA_SCATTER2D_V2 diff --git a/analyses/pluginNovosibirsk/SND_2016_I1473343.cc b/analyses/pluginNovosibirsk/SND_2016_I1473343.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2016_I1473343.cc @@ -0,0 +1,115 @@ +// -*- 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 SND_2016_I1473343 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(SND_2016_I1473343); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + // Initialise and register projections + declare(FinalState(), "FS"); + declare(UnstableParticles(), "UFS"); + _nOmegaEta= bookCounter("/TMP/nOmegaEta"); + } + + void findChildren(const Particle & p,map & nRes, int &ncount) { + foreach(const Particle &child, p.children()) { + if(child.children().empty()) { + nRes[child.pdgId()]-=1; + --ncount; + } + else + findChildren(child,nRes,ncount); + } + } + + /// Perform the per-event analysis + void analyze(const Event& event) { + const FinalState& fs = apply(event, "FS"); + map nCount; + int ntotal(0); + foreach (const Particle& p, fs.particles()) { + nCount[p.pdgId()] += 1; + ++ntotal; + } + const FinalState& ufs = apply(event, "UFS"); + foreach (const Particle& p, ufs.particles()) { + if(p.children().empty()) continue; + if(p.pdgId()!=223) continue; + map nRes = nCount; + int ncount = ntotal; + findChildren(p,nRes,ncount); + foreach (const Particle& p2, ufs.particles()) { + if(p2.pdgId()!=221) continue; + map nResB = nRes; + int ncountB = ncount; + findChildren(p2,nResB,ncountB); + if(ncountB!=0) continue; + bool matched2 = true; + for(auto const & val : nResB) { + if(val.second!=0) { + matched2 = false; + break; + } + } + if(matched2) { + _nOmegaEta->fill(event.weight()); + break; + } + } + } + } + + + /// Normalise histograms etc., after the run + void finalize() { + double sigma = _nOmegaEta->val(); + double error = _nOmegaEta->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 ex = temphisto.point(b).xErrs(); + pair 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 _nOmegaEta; + //@} + + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(SND_2016_I1473343); + + +} diff --git a/analyses/pluginNovosibirsk/SND_2016_I1473343.info b/analyses/pluginNovosibirsk/SND_2016_I1473343.info new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2016_I1473343.info @@ -0,0 +1,33 @@ +Name: SND_2016_I1473343 +Year: 2016 +Summary: Cross section for $e^+e^-\to\eta\omega$ for energies below 2 GeV +Experiment: SND +Collider: VEPP-2M +InspireID: 1473343 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - Phys.Rev. D94 (2016) no.9, 092002 +RunInfo: e+ e- to hadrons +NeedCrossSection: yes +Beams: [e+, e-] +Description: + 'Cross section for $e^+e^-\to\eta\omega$ for energies below 2 GeV' +Keywords: [] +BibKey: Achasov:2016qvd +BibTeX: '@article{Achasov:2016qvd, + author = "Achasov, M. N. and others", + title = "{Measurement of the $e^+e^- \to \omega\eta$ cross section + below $\sqrt{s}=2$ GeV}", + journal = "Phys. Rev.", + volume = "D94", + year = "2016", + number = "9", + pages = "092002", + doi = "10.1103/PhysRevD.94.092002", + eprint = "1607.00371", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + SLACcitation = "%%CITATION = ARXIV:1607.00371;%%" +}' diff --git a/analyses/pluginNovosibirsk/SND_2016_I1473343.plot b/analyses/pluginNovosibirsk/SND_2016_I1473343.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2016_I1473343.plot @@ -0,0 +1,7 @@ +BEGIN PLOT /SND_2016_I1473343/d01-x01-y01 +Title=$\sigma(e^+e^-\to \eta\omega)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to \eta\omega)$/nb +ConnectGaps=1 +LogY=0 +END PLOT diff --git a/analyses/pluginNovosibirsk/SND_2016_I1473343.yoda b/analyses/pluginNovosibirsk/SND_2016_I1473343.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2016_I1473343.yoda @@ -0,0 +1,21 @@ +BEGIN YODA_SCATTER2D_V2 /REF/SND_2016_I1473343/d01-x01-y01 +IsRef: 1 +Path: /REF/SND_2016_I1473343/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.36700 0.02700 0.02700 -7.00000e-02 5.00000e-02 5.00000e-02 +1.41900 0.02500 0.02500 1.90000e-01 6.00000e-02 6.00000e-02 +1.46900 0.02500 0.02500 3.50000e-01 1.00000e-01 1.00000e-01 +1.51900 0.02500 0.02500 1.50000e-01 2.14709e-01 1.00000e-01 +1.56900 0.02500 0.02500 8.70000e-01 1.10454e-01 1.10454e-01 +1.61900 0.02500 0.02500 1.26000e+00 2.15870e-01 2.10238e-01 +1.66900 0.02500 0.02500 2.01000e+00 2.23607e-01 2.02237e-01 +1.71900 0.02500 0.02500 8.70000e-01 2.34307e-01 1.81108e-01 +1.76900 0.02500 0.02500 4.40000e-01 1.74642e-01 1.72627e-01 +1.81900 0.02500 0.02500 1.70000e-01 1.43178e-01 1.41421e-01 +1.86900 0.02500 0.02500 -3.00000e-02 1.51327e-01 1.50000e-01 +1.91900 0.02500 0.02500 -6.00000e-02 1.43178e-01 1.40357e-01 +1.97200 0.02800 0.02800 -6.00000e-02 1.72627e-01 1.70294e-01 +END YODA_SCATTER2D_V2 diff --git a/analyses/pluginNovosibirsk/SND_2016_I1484677.cc b/analyses/pluginNovosibirsk/SND_2016_I1484677.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2016_I1484677.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 SND_2016_I1484677 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(SND_2016_I1484677); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + + // Initialise and register projections + declare(FinalState(), "FS"); + + // Book histograms + _nkaon = bookCounter("TMP/kaon"); + } + + + /// Perform the per-event analysis + void analyze(const Event& event) { + const FinalState& fs = apply(event, "FS"); + if(fs.particles().size()!=2) vetoEvent; + foreach (const Particle& p, fs.particles()) { + if(abs(p.pdgId())!=PID::KPLUS) vetoEvent; + } + _nkaon->fill(event.weight()); + } + + + /// Normalise histograms etc., after the run + void finalize() { + double sigma = _nkaon->val(); + double error = _nkaon->err(); + sigma *= crossSection()/ sumOfWeights() /nanobarn; + error *= crossSection()/ sumOfWeights() /nanobarn; + for(unsigned int ix=1;ix<3;++ix) { + 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 ex = temphisto.point(b).xErrs(); + pair 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 _nkaon; + //@} + + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(SND_2016_I1484677); + + +} diff --git a/analyses/pluginNovosibirsk/SND_2016_I1484677.info b/analyses/pluginNovosibirsk/SND_2016_I1484677.info new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2016_I1484677.info @@ -0,0 +1,34 @@ +Name: SND_2016_I1484677 +Year: 2016 +Summary: Cross section for $e^+e^-\to K^+K^-$ between 1.05 and 2 GeV +Experiment: SND +Collider: VEPP-2M +InspireID: 1484677 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - Phys.Rev. D94 (2016) no.11, 112006 +RunInfo: e+ e- to hadrons +NeedCrossSection: yes +Beams: [e+, e-] +Description: + 'Measurement of the cross section for $e^+e^-\to K^+K^-$ at energies between 1.05 and 2 GeV' +Keywords: +BibKey: Achasov:2016lbc +BibTeX: '@article{Achasov:2016lbc, + author = "Achasov, M. N. and others", + title = "{Measurement of the $\mathbf{e^+e^-\to K^+K^-}$ cross + section in the energy range $\mathbf{\sqrt{s}=1.05-2.0}$ + GeV}", + journal = "Phys. Rev.", + volume = "D94", + year = "2016", + number = "11", + pages = "112006", + doi = "10.1103/PhysRevD.94.112006", + eprint = "1608.08757", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + SLACcitation = "%%CITATION = ARXIV:1608.08757;%%" +}' diff --git a/analyses/pluginNovosibirsk/SND_2016_I1484677.plot b/analyses/pluginNovosibirsk/SND_2016_I1484677.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2016_I1484677.plot @@ -0,0 +1,12 @@ +BEGIN PLOT /SND_2016_I1484677/d01-x01-y01 +Title=$\sigma(e^+e^-\to K^+K^-)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to K^+K^-)$/nb +ConnectGaps=1 +END PLOT +BEGIN PLOT /SND_2016_I1484677/d02-x01-y01 +Title=$\sigma(e^+e^-\to K^+K^-)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to K^+K^-)$/nb +ConnectGaps=1 +END PLOT diff --git a/analyses/pluginNovosibirsk/SND_2016_I1484677.yoda b/analyses/pluginNovosibirsk/SND_2016_I1484677.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2016_I1484677.yoda @@ -0,0 +1,72 @@ +BEGIN YODA_SCATTER2D_V2 /REF/SND_2016_I1484677/d01-x01-y01 +IsRef: 1 +Path: /REF/SND_2016_I1484677/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.04700 0.00000 0.00000 3.62430e+01 7.56508e-01 7.56508e-01 +1.07500 0.00000 0.00000 1.84830e+01 3.88104e-01 3.88104e-01 +1.09700 0.00000 0.00000 1.41840e+01 3.05601e-01 3.05601e-01 +1.12400 0.00000 0.00000 1.13070e+01 2.47407e-01 2.47407e-01 +1.15100 0.00000 0.00000 9.41500e+00 2.27482e-01 2.27482e-01 +1.17400 0.00000 0.00000 9.14800e+00 2.08331e-01 2.08331e-01 +1.19600 0.00000 0.00000 7.73700e+00 1.83330e-01 1.83330e-01 +1.22300 0.00000 0.00000 7.26900e+00 2.03571e-01 2.03571e-01 +1.24500 0.00000 0.00000 6.62400e+00 2.08406e-01 2.08406e-01 +1.27300 0.00000 0.00000 6.51800e+00 1.91888e-01 1.91888e-01 +1.29500 0.00000 0.00000 5.96400e+00 2.01040e-01 2.01040e-01 +1.32300 0.00000 0.00000 5.71800e+00 1.72177e-01 1.72177e-01 +1.34400 0.00000 0.00000 5.16000e+00 1.63294e-01 1.63294e-01 +1.37400 0.00000 0.00000 5.30800e+00 1.50360e-01 1.50360e-01 +1.39400 0.00000 0.00000 5.00100e+00 1.47139e-01 1.47139e-01 +1.42300 0.00000 0.00000 4.57900e+00 1.36418e-01 1.36418e-01 +1.44300 0.00000 0.00000 4.25400e+00 1.48385e-01 1.48385e-01 +1.47100 0.00000 0.00000 4.09300e+00 1.24326e-01 1.24326e-01 +1.49400 0.00000 0.00000 3.64600e+00 1.20420e-01 1.20420e-01 +1.52200 0.00000 0.00000 3.67900e+00 1.33701e-01 1.33701e-01 +1.54300 0.00000 0.00000 3.38200e+00 1.30035e-01 1.30035e-01 +1.57200 0.00000 0.00000 3.51800e+00 1.30173e-01 1.30173e-01 +1.59400 0.00000 0.00000 3.50700e+00 1.59101e-01 1.59101e-01 +1.62300 0.00000 0.00000 3.02200e+00 1.10693e-01 1.10693e-01 +1.64300 0.00000 0.00000 2.81500e+00 1.12947e-01 1.12947e-01 +1.66900 0.00000 0.00000 2.15500e+00 9.61353e-02 9.61353e-02 +1.69300 0.00000 0.00000 1.57000e+00 8.42021e-02 8.42021e-02 +1.72300 0.00000 0.00000 9.68000e-01 7.18401e-02 7.18401e-02 +1.74200 0.00000 0.00000 6.33000e-01 5.26308e-02 5.26308e-02 +1.77400 0.00000 0.00000 3.10000e-01 3.71080e-02 3.71080e-02 +1.79300 0.00000 0.00000 1.70000e-01 3.59026e-02 3.59026e-02 +1.82600 0.00000 0.00000 2.15000e-01 3.60555e-02 3.60555e-02 +1.84900 0.00000 0.00000 1.71000e-01 4.29418e-02 4.29418e-02 +1.87100 0.00000 0.00000 2.91000e-01 3.98121e-02 3.98121e-02 +1.89300 0.00000 0.00000 3.45000e-01 4.30813e-02 4.30813e-02 +1.90100 0.00000 0.00000 3.43000e-01 5.06952e-02 5.06952e-02 +1.92700 0.00000 0.00000 3.16000e-01 4.42719e-02 4.42719e-02 +1.95300 0.00000 0.00000 3.57000e-01 6.28172e-02 6.28172e-02 +1.97800 0.00000 0.00000 3.32000e-01 5.75674e-02 5.75674e-02 +2.00500 0.00000 0.00000 3.67000e-01 5.05964e-02 5.05964e-02 +END YODA_SCATTER2D_V2 +BEGIN YODA_SCATTER2D_V2 /REF/SND_2016_I1484677/d02-x01-y01 +IsRef: 1 +Path: /REF/SND_2016_I1484677/d02-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.27700 0.00000 0.00000 6.35800e+00 1.76604e-01 1.76604e-01 +1.35700 0.00000 0.00000 5.33900e+00 1.51013e-01 1.51013e-01 +1.43500 0.00000 0.00000 4.11400e+00 1.20673e-01 1.20673e-01 +1.51500 0.00000 0.00000 3.57400e+00 1.11517e-01 1.11517e-01 +1.59500 0.00000 0.00000 3.30200e+00 1.05475e-01 1.05475e-01 +1.67400 0.00000 0.00000 1.77000e+00 6.93181e-02 6.93181e-02 +1.71600 0.00000 0.00000 1.11500e+00 7.90569e-02 7.90569e-02 +1.75800 0.00000 0.00000 5.03000e-01 4.04969e-02 4.04969e-02 +1.79800 0.00000 0.00000 1.87000e-01 2.96816e-02 2.96816e-02 +1.84000 0.00000 0.00000 2.93000e-01 4.33244e-02 4.33244e-02 +1.87400 0.00000 0.00000 2.76000e-01 3.66742e-02 3.66742e-02 +1.90300 0.00000 0.00000 2.87000e-01 3.71214e-02 3.71214e-02 +1.92600 0.00000 0.00000 3.67000e-01 4.68188e-02 4.68188e-02 +1.94500 0.00000 0.00000 2.63000e-01 3.71618e-02 3.71618e-02 +1.96700 0.00000 0.00000 3.31000e-01 4.36578e-02 4.36578e-02 +1.98900 0.00000 0.00000 3.46000e-01 4.83839e-02 4.83839e-02 +END YODA_SCATTER2D_V2 diff --git a/analyses/pluginNovosibirsk/SND_2016_I1489182.cc b/analyses/pluginNovosibirsk/SND_2016_I1489182.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2016_I1489182.cc @@ -0,0 +1,111 @@ +// -*- 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 SND_2016_I1489182 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(SND_2016_I1489182); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + + // Initialise and register projections + declare(FinalState(), "FS"); + declare(UnstableParticles(), "UFS"); + _numOmegaPi = bookCounter("TMP/OmegaPi"); + } + + + + void findChildren(const Particle & p,map & nRes, int &ncount) { + foreach(const Particle &child, p.children()) { + if(child.children().empty()) { + nRes[child.pdgId()]+=1; + ++ncount; + } + else + findChildren(child,nRes,ncount); + } + } + /// Perform the per-event analysis + void analyze(const Event& event) { + const FinalState& fs = apply(event, "FS"); + + map nCount; + int ntotal(0); + foreach (const Particle& p, fs.particles()) { + nCount[p.pdgId()] += 1; + ++ntotal; + } + // three particles (pi0 pi0 gamma) + if(ntotal!=3) vetoEvent; + const FinalState& ufs = apply(event, "UFS"); + foreach (const Particle& p, ufs.particles()) { + if(p.children().empty()) continue; + // find the omega + if(p.pdgId()==223) { + map nRes; + int ncount(0); + findChildren(p,nRes,ncount); + // only omega to pi0 gamma mode + if(ncount!=2) continue; + if(nRes[111]!=1 || nRes[22]!=1) continue; + // omega pi0 + if(nCount[111]-nRes[111]==1) + _numOmegaPi->fill(event.weight()); + } + } + } + + + /// Normalise histograms etc., after the run + void finalize() { + + double sigma = _numOmegaPi->val(); + double error = _numOmegaPi->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 ex = temphisto.point(b).xErrs(); + pair 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 _numOmegaPi; + //@} + + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(SND_2016_I1489182); + + +} diff --git a/analyses/pluginNovosibirsk/SND_2016_I1489182.info b/analyses/pluginNovosibirsk/SND_2016_I1489182.info new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2016_I1489182.info @@ -0,0 +1,35 @@ +Name: SND_2016_I1489182 +Year: 2016 +Summary: Cross section for $e^+e^-\to\omega\pi^0\to\pi^0\pi^0\gamma$ between 1.05 and 2 GeV +Experiment: SND +Collider: VEPP-2M +InspireID: 1489182 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - Phys.Rev. D94 (2016) no.11, 112001 +RunInfo: e+e- to hadrons +NeedCrossSection: yes +Beams: [e+,e-] +Luminosity_fb: +Description: +Description: + 'Measurement of the cross section for $e^+e^-\to\omega\pi^0\to\pi^0\pi^0\gamma$ at energies between 1.05 and 2 GeV.' +Keywords: +BibKey: Achasov:2016zvn +BibTeX: '@article{Achasov:2016zvn, + author = "Achasov, M. N. and others", + title = "{Updated measurement of the $e^+e^- \to \omega \pi^0 \to + \pi^0\pi^0\gamma$ cross section with the SND detector}", + journal = "Phys. Rev.", + volume = "D94", + year = "2016", + number = "11", + pages = "112001", + doi = "10.1103/PhysRevD.94.112001", + eprint = "1610.00235", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + SLACcitation = "%%CITATION = ARXIV:1610.00235;%%" +}' diff --git a/analyses/pluginNovosibirsk/SND_2016_I1489182.plot b/analyses/pluginNovosibirsk/SND_2016_I1489182.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2016_I1489182.plot @@ -0,0 +1,7 @@ +BEGIN PLOT /SND_2016_I1489182/d01-x01-y01 +Title=$\sigma(e^+e^-\to \omega\pi^0\to \pi^0\pi^0\gamma)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to \omega\pi^0\to \pi^0\pi^0\gamma)$/nb +LogY=0 +ConnectGaps=1 +END PLOT diff --git a/analyses/pluginNovosibirsk/SND_2016_I1489182.yoda b/analyses/pluginNovosibirsk/SND_2016_I1489182.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2016_I1489182.yoda @@ -0,0 +1,64 @@ +BEGIN YODA_SCATTER2D_V2 /REF/SND_2016_I1489182/d01-x01-y01 +IsRef: 1 +Path: /REF/SND_2016_I1489182/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1047.00000 6.00000 6.00000 0.95000 0.11000 0.11000 +1075.00000 6.00000 6.00000 1.04000 0.09000 0.09000 +1097.00000 6.00000 6.00000 1.08000 0.07000 0.07000 +1124.00000 6.00000 6.00000 1.26000 0.11000 0.11000 +1151.00000 6.00000 6.00000 1.30000 0.10000 0.10000 +1174.00000 6.00000 6.00000 1.29000 0.12000 0.12000 +1196.00000 6.00000 6.00000 1.39000 0.07000 0.07000 +1223.00000 6.00000 6.00000 1.33000 0.11000 0.11000 +1245.00000 6.00000 6.00000 1.25000 0.12000 0.12000 +1273.00000 6.00000 2.00000 1.53000 0.13000 0.13000 +1277.00000 2.00000 2.00000 1.38000 0.09000 0.09000 +1295.00000 6.00000 6.00000 1.56000 0.07000 0.07000 +1323.00000 6.00000 6.00000 1.54000 0.12000 0.12000 +1344.00000 6.00000 6.00000 1.50000 0.13000 0.13000 +1357.00000 2.00000 2.00000 1.60000 0.11000 0.11000 +1374.00000 6.00000 6.00000 1.55000 0.12000 0.12000 +1394.00000 6.00000 6.00000 1.61000 0.09000 0.09000 +1423.00000 6.00000 6.00000 1.78000 0.10000 0.10000 +1435.00000 2.00000 2.00000 1.67000 0.09000 0.09000 +1443.00000 6.00000 6.00000 1.47000 0.12000 0.12000 +1471.00000 6.00000 6.00000 1.37000 0.10000 0.10000 +1494.00000 6.00000 6.00000 1.42000 0.06000 0.06000 +1515.00000 2.00000 2.00000 1.21000 0.08000 0.08000 +1522.00000 5.00000 6.00000 1.48000 0.10000 0.10000 +1543.00000 6.00000 6.00000 1.13000 0.08000 0.08000 +1572.00000 6.00000 6.00000 0.93000 0.08000 0.08000 +1594.00000 6.00000 0.50000 0.80000 0.06000 0.06000 +1595.00000 0.50000 2.00000 0.77000 0.07000 0.07000 +1623.00000 6.00000 6.00000 0.70000 0.08000 0.08000 +1643.00000 6.00000 6.00000 0.54000 0.07000 0.07000 +1669.00000 6.00000 3.00000 0.45000 0.07000 0.07000 +1674.00000 2.00000 2.00000 0.41000 0.04000 0.04000 +1693.00000 6.00000 6.00000 0.35000 0.04000 0.04000 +1716.00000 2.00000 2.00000 0.30000 0.05000 0.05000 +1723.00000 5.00000 6.00000 0.17000 0.04000 0.04000 +1742.00000 6.00000 6.00000 0.17000 0.04000 0.04000 +1758.00000 2.00000 2.00000 0.17000 0.04000 0.04000 +1774.00000 6.00000 6.00000 0.20000 0.04000 0.04000 +1793.00000 6.00000 3.00000 0.19000 0.03000 0.03000 +1798.00000 2.00000 2.00000 0.14000 0.03000 0.03000 +1826.00000 6.00000 6.00000 0.09000 0.03000 0.03000 +1840.00000 2.00000 2.00000 0.12000 0.03000 0.03000 +1849.00000 6.00000 6.00000 0.07000 0.04000 0.04000 +1871.00000 6.00000 1.50000 0.10000 0.03000 0.03000 +1874.00000 1.50000 2.00000 0.12000 0.02000 0.02000 +1893.00000 6.00000 4.00000 0.13000 0.03000 0.03000 +1901.00000 4.00000 1.00000 0.06000 0.02000 0.02000 +1903.00000 1.00000 2.00000 0.08000 0.03000 0.03000 +1926.00000 2.00000 0.50000 0.08000 0.03000 0.03000 +1927.00000 0.50000 6.00000 0.10000 0.03000 0.03000 +1945.00000 2.00000 2.00000 0.06000 0.02000 0.02000 +1953.00000 6.00000 6.00000 0.03000 0.03000 0.03000 +1967.00000 2.00000 2.00000 0.06000 0.03000 0.03000 +1978.00000 6.00000 6.00000 0.05000 0.05000 0.02000 +1989.00000 2.00000 2.00000 0.06000 0.02000 0.02000 +2005.00000 6.00000 6.00000 0.07000 0.03000 0.03000 +END YODA_SCATTER2D_V2 diff --git a/analyses/pluginNovosibirsk/SND_2018_I1637194.cc b/analyses/pluginNovosibirsk/SND_2018_I1637194.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2018_I1637194.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 SND_2018_I1637194 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(SND_2018_I1637194); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + // Initialise and register projections + declare(FinalState(), "FS"); + _nKKpi= bookCounter( "/TMP/nKKpi" ); + } + + + /// Perform the per-event analysis + void analyze(const Event& event) { + const FinalState& fs = apply(event, "FS"); + + map nCount; + int ntotal(0); + foreach (const Particle& p, fs.particles()) { + nCount[p.pdgId()] += 1; + ++ntotal; + } + + if(ntotal==3 && nCount[130]==1 && nCount[310]==1 && nCount[111] ==1) + _nKKpi->fill(event.weight()); + } + + + /// Normalise histograms etc., after the run + void finalize() { + double sigma = _nKKpi->val(); + double error = _nKKpi->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 ex = temphisto.point(b).xErrs(); + pair 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 _nKKpi; + //@} + + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(SND_2018_I1637194); + + +} diff --git a/analyses/pluginNovosibirsk/SND_2018_I1637194.info b/analyses/pluginNovosibirsk/SND_2018_I1637194.info new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2018_I1637194.info @@ -0,0 +1,35 @@ +Name: SND_2018_I1637194 +Year: 2018 +Summary: Cross section for $e^+e^-\to$ $K^0_SK^0_L\pi^0$ for energies between 1.3 and 2.0 +Experiment: SND +Collider: VEPP-2M +InspireID: 1637194 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - Phys.Rev. D97 (2018) no.3, 032011 +RunInfo: e+e- to hadrons +NeedCrossSection: yes +Beams: [e+,e-] +Luminosity_fb: +Description: + 'Measurement of the cross section for $e^+e^-\to$ $K^0_SK^0_L\pi^0$ for energies from 1.3 to 2.0 GeV.' +Keywords: +BibKey: Achasov:2017vaq +BibTeX: '@article{Achasov:2017vaq, + author = "Achasov, M. N. and others", + title = "{Measurement of the $\mathbf{e^+e^- \rightarrow K_S K_L + \pi^0}$ cross section in the energy range + $\mathbf{\sqrt{s}=1.3-2.0}$ GeV}", + journal = "Phys. Rev.", + volume = "D97", + year = "2018", + number = "3", + pages = "032011", + doi = "10.1103/PhysRevD.97.032011", + eprint = "1711.07143", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + SLACcitation = "%%CITATION = ARXIV:1711.07143;%%" +}' diff --git a/analyses/pluginNovosibirsk/SND_2018_I1637194.plot b/analyses/pluginNovosibirsk/SND_2018_I1637194.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2018_I1637194.plot @@ -0,0 +1,6 @@ +BEGIN PLOT /SND_2018_I1637194/d01-x01-y01 +Title=$\sigma(e^+e^-\to K^0_SK^0_L \pi^0)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to K^0_SK^0_L \pi^0)$/nb +ConnectGaps=1 +END PLOT diff --git a/analyses/pluginNovosibirsk/SND_2018_I1637194.yoda b/analyses/pluginNovosibirsk/SND_2018_I1637194.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2018_I1637194.yoda @@ -0,0 +1,23 @@ +BEGIN YODA_SCATTER2D_V2 /REF/SND_2018_I1637194/d01-x01-y01 +IsRef: 1 +Path: /REF/SND_2018_I1637194/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.32500 0.02500 0.02500 1.40000e-01 5.00000e-02 5.00000e-02 +1.36750 0.00750 0.00750 1.40000e-01 6.00000e-02 6.00000e-02 +1.42000 0.02000 0.02000 5.20000e-01 8.00000e-02 8.00000e-02 +1.46250 0.01250 0.01250 9.00000e-01 1.60000e-01 1.60000e-01 +1.50000 0.00000 0.00000 1.35000e+00 1.30000e-01 1.30000e-01 +1.52250 0.00250 0.00250 1.49000e+00 1.80000e-01 1.80000e-01 +1.56250 0.01250 0.01250 2.48000e+00 2.40000e-01 2.40000e-01 +1.62500 0.02500 0.02500 2.80000e+00 1.90000e-01 1.90000e-01 +1.68750 0.01250 0.01250 2.85000e+00 2.00000e-01 2.00000e-01 +1.73500 0.01500 0.01500 2.30000e+00 2.10000e-01 2.10000e-01 +1.78000 0.02000 0.02000 1.81000e+00 1.50000e-01 1.50000e-01 +1.83750 0.01250 0.01250 1.00000e+00 2.30000e-01 2.30000e-01 +1.88500 0.01500 0.01500 5.50000e-01 9.00000e-02 9.00000e-02 +1.93500 0.01500 0.01500 3.90000e-01 1.20000e-01 1.20000e-01 +1.98000 0.02000 0.02000 4.10000e-01 1.50000e-01 1.50000e-01 +END YODA_SCATTER2D_V2 diff --git a/analyses/pluginNovosibirsk/SND_2018_I1638368.cc b/analyses/pluginNovosibirsk/SND_2018_I1638368.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2018_I1638368.cc @@ -0,0 +1,120 @@ +// -*- 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 SND_2018_I1638368 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(SND_2018_I1638368); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + + // Initialise and register projections + declare(FinalState(), "FS"); + declare(UnstableParticles(), "UFS"); + _numEtaPiPi = bookCounter("TMP/EtaPiPi"); + } + + + void findChildren(const Particle & p,map & nRes, int &ncount) { + foreach(const Particle &child, p.children()) { + if(child.children().empty()) { + --nRes[child.pdgId()]; + --ncount; + } + else + findChildren(child,nRes,ncount); + } + } + + /// Perform the per-event analysis + void analyze(const Event& event) { + const FinalState& fs = apply(event, "FS"); + + map nCount; + int ntotal(0); + foreach (const Particle& p, fs.particles()) { + nCount[p.pdgId()] += 1; + ++ntotal; + } + const FinalState& ufs = apply(event, "UFS"); + foreach (const Particle& p, ufs.particles()) { + if(p.children().empty()) continue; + // find the omega + if(p.pdgId()==221) { + map nRes = nCount; + int ncount = ntotal; + findChildren(p,nRes,ncount); + // eta pi+pi- + if(ncount!=2) continue; + bool matched = true; + for(auto const & val : nRes) { + if(abs(val.first)==211) { + if(val.second !=1) { + matched = false; + break; + } + } + else if(val.second!=0) { + matched = false; + break; + } + } + if(matched) + _numEtaPiPi->fill(event.weight()); + } + } + } + + + /// Normalise histograms etc., after the run + void finalize() { + double sigma = _numEtaPiPi->val(); + double error = _numEtaPiPi->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 ex = temphisto.point(b).xErrs(); + pair 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 _numEtaPiPi; + //@} + + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(SND_2018_I1638368); + + +} diff --git a/analyses/pluginNovosibirsk/SND_2018_I1638368.info b/analyses/pluginNovosibirsk/SND_2018_I1638368.info new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2018_I1638368.info @@ -0,0 +1,35 @@ +Name: SND_2018_I1638368 +Year: 2018 +Summary: Cross section for $e^+e^-\to\eta\pi^+\pi^-$ between 1.08 and 2 GeV +Experiment: SND +Collider: VEPP-2M +InspireID: 1638368 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - Phys.Rev. D97 (2018) no.1, 012008 +RunInfo: e+e- to hadrons +NeedCrossSection: yes +Beams: [e+,e-] +Luminosity_fb: +Description: +Description: + 'Measurement of the cross section for $e^+e^-\to\eta\pi^+\pi^-$ at energies between 1.08 and 2 GeV.' +Keywords: +BibKey: Achasov:2017kqm +BibTeX: '@article{Achasov:2017kqm, + author = "Achasov, M. N. and others", + title = "{Measurement of the $e^+e^- \to \eta\pi^+\pi^-$ cross + section with the SND detector at the VEPP-2000 collider}", + journal = "Phys. Rev.", + volume = "D97", + year = "2018", + number = "1", + pages = "012008", + doi = "10.1103/PhysRevD.97.012008", + eprint = "1711.08862", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + SLACcitation = "%%CITATION = ARXIV:1711.08862;%%" +}' diff --git a/analyses/pluginNovosibirsk/SND_2018_I1638368.plot b/analyses/pluginNovosibirsk/SND_2018_I1638368.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2018_I1638368.plot @@ -0,0 +1,8 @@ +BEGIN PLOT /SND_2018_I1638368/d01-x01-y01 +Title=$\sigma(e^+e^-\to \eta\pi^+\pi^-)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to \eta\pi^+\pi^-)$/nb +LogY=0 +LegendXPos=0.7 +ConnectGaps=1 +END PLOT diff --git a/analyses/pluginNovosibirsk/SND_2018_I1638368.yoda b/analyses/pluginNovosibirsk/SND_2018_I1638368.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2018_I1638368.yoda @@ -0,0 +1,47 @@ +BEGIN YODA_SCATTER2D_V2 /REF/SND_2018_I1638368/d01-x01-y01 +IsRef: 1 +Path: /REF/SND_2018_I1638368/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.07500 0.00000 0.00000 3.70000e-01 2.50000e-01 2.50000e-01 +1.09700 0.00000 0.00000 9.00000e-02 1.50000e-01 1.50000e-01 +1.12400 0.00000 0.00000 1.30000e-01 2.00000e-01 2.00000e-01 +1.15100 0.00000 0.00000 1.70000e-01 1.60000e-01 1.60000e-01 +1.17400 0.00000 0.00000 8.00000e-02 1.90000e-01 1.90000e-01 +1.19600 0.00000 0.00000 2.10000e-01 1.60000e-01 1.60000e-01 +1.22300 0.00000 0.00000 3.30000e-01 1.30000e-01 1.30000e-01 +1.24500 0.00000 0.00000 1.50000e-01 1.20000e-01 1.20000e-01 +1.27500 0.00000 0.00000 3.30000e-01 9.00000e-02 9.00000e-02 +1.29500 0.00000 0.00000 5.10000e-01 1.50000e-01 1.50000e-01 +1.32300 0.00000 0.00000 7.10000e-01 1.60000e-01 1.60000e-01 +1.35100 0.00000 0.00000 1.02000e+00 1.20000e-01 1.20000e-01 +1.37400 0.00000 0.00000 1.22000e+00 1.90000e-01 1.90000e-01 +1.39400 0.00000 0.00000 1.86000e+00 2.00000e-01 2.00000e-01 +1.42300 0.00000 0.00000 2.06000e+00 2.10000e-01 2.10000e-01 +1.43800 0.00000 0.00000 3.03000e+00 1.60000e-01 1.60000e-01 +1.47100 0.00000 0.00000 3.29000e+00 2.50000e-01 2.50000e-01 +1.49400 0.00000 0.00000 3.82000e+00 2.40000e-01 2.40000e-01 +1.51700 0.00000 0.00000 4.44000e+00 1.90000e-01 1.90000e-01 +1.54300 0.00000 0.00000 4.55000e+00 2.80000e-01 2.80000e-01 +1.57200 0.00000 0.00000 3.94000e+00 3.00000e-01 3.00000e-01 +1.59400 0.00000 0.00000 3.34000e+00 3.10000e-01 3.10000e-01 +1.62300 0.00000 0.00000 3.12000e+00 2.80000e-01 2.80000e-01 +1.64300 0.00000 0.00000 2.45000e+00 2.80000e-01 2.80000e-01 +1.67200 0.00000 0.00000 2.30000e+00 1.60000e-01 1.60000e-01 +1.69300 0.00000 0.00000 2.67000e+00 2.70000e-01 2.70000e-01 +1.72000 0.00000 0.00000 2.23000e+00 1.70000e-01 1.70000e-01 +1.75100 0.00000 0.00000 2.36000e+00 1.70000e-01 1.70000e-01 +1.77400 0.00000 0.00000 1.96000e+00 2.50000e-01 2.50000e-01 +1.79700 0.00000 0.00000 2.00000e+00 1.60000e-01 1.60000e-01 +1.82600 0.00000 0.00000 1.44000e+00 2.20000e-01 2.20000e-01 +1.84300 0.00000 0.00000 1.31000e+00 1.40000e-01 1.40000e-01 +1.87300 0.00000 0.00000 9.70000e-01 1.30000e-01 1.30000e-01 +1.90000 0.00000 0.00000 7.90000e-01 1.00000e-01 1.00000e-01 +1.92700 0.00000 0.00000 7.60000e-01 1.30000e-01 1.30000e-01 +1.94700 0.00000 0.00000 7.10000e-01 1.20000e-01 1.20000e-01 +1.96700 0.00000 0.00000 7.40000e-01 1.80000e-01 1.80000e-01 +1.98400 0.00000 0.00000 7.40000e-01 1.50000e-01 1.50000e-01 +2.00500 0.00000 0.00000 7.80000e-01 2.20000e-01 2.20000e-01 +END YODA_SCATTER2D_V2 diff --git a/analyses/pluginNovosibirsk/SND_2018_I1693737.cc b/analyses/pluginNovosibirsk/SND_2018_I1693737.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2018_I1693737.cc @@ -0,0 +1,115 @@ +// -*- 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 SND_2018_I1693737 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(SND_2018_I1693737); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + // Initialise and register projections + declare(FinalState(), "FS"); + declare(UnstableParticles(), "UFS"); + _nKKEta= bookCounter("/TMP/nKKEta"); + } + + void findChildren(const Particle & p,map & nRes, int &ncount) { + foreach(const Particle &child, p.children()) { + if(child.children().empty()) { + nRes[child.pdgId()]-=1; + --ncount; + } + else + findChildren(child,nRes,ncount); + } + } + + /// Perform the per-event analysis + void analyze(const Event& event) { + const FinalState& fs = apply(event, "FS"); + map nCount; + int ntotal(0); + foreach (const Particle& p, fs.particles()) { + nCount[p.pdgId()] += 1; + ++ntotal; + } + const FinalState& ufs = apply(event, "UFS"); + foreach (const Particle& p, ufs.particles()) { + if(p.children().empty()) continue; + if(p.pdgId()!=221) continue; + map nRes = nCount; + int ncount = ntotal; + findChildren(p,nRes,ncount); + if(ncount!=2) continue; + bool matched = true; + for(auto const & val : nRes) { + if(abs(val.first)==321) { + if(val.second!=1) { + matched = false; + break; + } + } + else if(val.second!=0) { + matched = false; + break; + } + } + if(matched) { + _nKKEta->fill(event.weight()); + break; + } + } + } + + /// Normalise histograms etc., after the run + void finalize() { + double sigma = _nKKEta->val(); + double error = _nKKEta->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 ex = temphisto.point(b).xErrs(); + pair 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 _nKKEta; + //@} + + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(SND_2018_I1693737); + + +} diff --git a/analyses/pluginNovosibirsk/SND_2018_I1693737.info b/analyses/pluginNovosibirsk/SND_2018_I1693737.info new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2018_I1693737.info @@ -0,0 +1,34 @@ +Name: SND_2018_I1693737 +Year: 2018 +Summary: Cross section for $e^+e^-\to\eta K^+K^-$ for energies below 2 GeV +Experiment: SND +Collider: VEPP-2M +InspireID: 1693737 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - Phys.Atom.Nucl. 81 (2018) no.2, 205-213 +RunInfo: e+ e- to hadrons +NeedCrossSection: yes +Beams: [e+, e-] +Description: + 'Cross section for $e^+e^-\to\eta K^+K^-$ for energies below 2 GeV' +Keywords: [] +BibKey: Achasov:2018ygm +BibTeX: '%%% contains utf-8, see: http://inspirehep.net/info/faq/general#utf8 +%%% add \usepackage[utf8]{inputenc} to your latex preamble + +@article{Achasov:2018ygm, + author = "Achasov, M. N. and others", + title = "{Measurement of the $e^+e^−\to\eta K^+K^−$ Cross + Section by Means of the SND Detector}", + journal = "Phys. Atom. Nucl.", + volume = "81", + year = "2018", + number = "2", + pages = "205-213", + doi = "10.1134/S1063778818020023", + note = "[Yad. Fiz.81,no.2,195(2018)]", + SLACcitation = "%%CITATION = PANUE,81,205;%%" +}' diff --git a/analyses/pluginNovosibirsk/SND_2018_I1693737.plot b/analyses/pluginNovosibirsk/SND_2018_I1693737.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2018_I1693737.plot @@ -0,0 +1,7 @@ +BEGIN PLOT /SND_2018_I1693737/d01-x01-y01 +Title=$\sigma(e^+e^-\to \eta K^+K^-)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to \eta K^+K^-)$/nb +ConnectGaps=1 +LogY=0 +END PLOT diff --git a/analyses/pluginNovosibirsk/SND_2018_I1693737.yoda b/analyses/pluginNovosibirsk/SND_2018_I1693737.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2018_I1693737.yoda @@ -0,0 +1,20 @@ +BEGIN YODA_SCATTER2D_V2 /REF/SND_2018_I1693737/d01-x01-y01 +IsRef: 1 +Path: /REF/SND_2018_I1693737/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.59000 0.03000 0.03000 0.00000e+00 7.07107e-02 9.05539e-02 +1.64000 0.02000 0.02000 4.10000e-01 1.31529e-01 1.31529e-01 +1.68000 0.02000 0.02000 1.09000e+00 1.70880e-01 1.70880e-01 +1.72000 0.02000 0.02000 9.30000e-01 1.58114e-01 1.58114e-01 +1.75500 0.01500 0.01500 9.40000e-01 2.06155e-01 2.06155e-01 +1.79500 0.02500 0.02500 5.80000e-01 1.20830e-01 1.20830e-01 +1.83500 0.01500 0.01500 4.50000e-01 1.26491e-01 1.26491e-01 +1.87000 0.02000 0.02000 3.10000e-01 9.48683e-02 9.48683e-02 +1.89500 0.00500 0.00500 3.30000e-01 1.33417e-01 1.33417e-01 +1.92500 0.02500 0.02500 3.70000e-01 1.04403e-01 1.04403e-01 +1.96000 0.01000 0.01000 7.00000e-02 7.07107e-02 5.09902e-02 +1.98500 0.01500 0.01500 2.30000e-01 1.01119e-01 1.01119e-01 +END YODA_SCATTER2D_V2 diff --git a/analyses/pluginNovosibirsk/SND_2018_I1694988.cc b/analyses/pluginNovosibirsk/SND_2018_I1694988.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2018_I1694988.cc @@ -0,0 +1,81 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/Projections/FinalState.hh" + +namespace Rivet { + + + /// @brief Add a short analysis description here + class SND_2018_I1694988 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(SND_2018_I1694988); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + declare(FinalState(), "FS"); + _numPi0Gamma = bookCounter("TMP/Pi0Gamma"); + } + + + /// Perform the per-event analysis + void analyze(const Event& event) { + + const FinalState& fs = apply(event, "FS"); + + map nCount; + int ntotal(0); + foreach (const Particle& p, fs.particles()) { + nCount[p.pdgId()] += 1; + ++ntotal; + } + if(ntotal==2 && nCount[22]==1 && nCount[111]==1) + _numPi0Gamma->fill(event.weight()); + } + + + /// Normalise histograms etc., after the run + void finalize() { + double sigma = _numPi0Gamma->val(); + double error = _numPi0Gamma->err(); + sigma *= crossSection()/ sumOfWeights() /picobarn; + error *= crossSection()/ sumOfWeights() /picobarn; + 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 ex = temphisto.point(b).xErrs(); + pair 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 _numPi0Gamma; + //@} + + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(SND_2018_I1694988); + + +} diff --git a/analyses/pluginNovosibirsk/SND_2018_I1694988.info b/analyses/pluginNovosibirsk/SND_2018_I1694988.info new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2018_I1694988.info @@ -0,0 +1,28 @@ +Name: SND_2018_I1694988 +Year: 2018 +Summary: Cross section for $e^+e^-\to\pi\gamma$ for energies between 1.075 GeV and 2 GeV +Experiment: SND +Collider: VEPP-2M +InspireID: 1694988 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - arXiv:1809.07631 +RunInfo: e+ e- to hadrons +NeedCrossSection: yes +Beams: [e+, e-] +Description: + 'Cross section for $e^+e^-\to\pi\gamma$ for energies between 1.075 GeV and 2 GeV ' +Keywords: [] +BibKey: Achasov:2018ujw +BibTeX: '@article{Achasov:2018ujw, + author = "Achasov, M. N. and others", + title = "{Measurement of the $e^+e^- \to \pi^0\gamma$ cross + section in the energy range 1.075--2 GeV at SND}", + year = "2018", + eprint = "1809.07631", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + SLACcitation = "%%CITATION = ARXIV:1809.07631;%%" +}' diff --git a/analyses/pluginNovosibirsk/SND_2018_I1694988.plot b/analyses/pluginNovosibirsk/SND_2018_I1694988.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2018_I1694988.plot @@ -0,0 +1,7 @@ +BEGIN PLOT /SND_2018_I1694988/d01-x01-y01 +Title=$\sigma(e^+e^-\to \pi^0\gamma)$ +XLabel=$\sqrt{s}$/MeV +YLabel=$\sigma(e^+e^-\to \pi^0\gamma)$/pb +LogY=1 +ConnectGaps=1 +END PLOT diff --git a/analyses/pluginNovosibirsk/SND_2018_I1694988.yoda b/analyses/pluginNovosibirsk/SND_2018_I1694988.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginNovosibirsk/SND_2018_I1694988.yoda @@ -0,0 +1,19 @@ +BEGIN YODA_SCATTER2D_V2 /REF/SND_2018_I1694988/d01-x01-y01 +IsRef: 1 +Path: /REF/SND_2018_I1694988/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1075.00000 0.00000 0.00000 9.30000e+01 4.51774e+01 4.51774e+01 +1119.00000 0.00000 0.00000 4.80000e+01 2.91548e+01 2.91548e+01 +1200.00000 0.00000 0.00000 6.00000e+01 2.60768e+01 2.60768e+01 +1284.00000 0.00000 0.00000 5.00000e+01 2.10950e+01 2.10950e+01 +1353.00000 0.00000 0.00000 2.30000e+01 1.90263e+01 1.90263e+01 +1425.00000 0.00000 0.00000 -1.30000e+01 1.50333e+01 1.50333e+01 +1509.00000 0.00000 0.00000 -1.10000e+01 1.30384e+01 1.30384e+01 +1607.00000 0.00000 0.00000 1.20000e+01 1.61245e+01 1.61245e+01 +1705.00000 0.00000 0.00000 -1.00000e+00 1.31529e+01 1.31529e+01 +1804.00000 0.00000 0.00000 4.00000e+00 1.50333e+01 1.50333e+01 +1926.00000 0.00000 0.00000 -8.00000e+00 1.30384e+01 1.30384e+01 +END YODA_SCATTER2D_V2