diff --git a/analyses/pluginBABAR/BABAR_2007_I760730.cc b/analyses/pluginBABAR/BABAR_2007_I760730.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2007_I760730.cc @@ -0,0 +1,159 @@ +// -*- 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 BABAR_2007_I760730 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(BABAR_2007_I760730); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + // Initialise and register projections + declare(FinalState(), "FS"); + declare(UnstableParticles(), "UFS"); + _nLL = bookCounter( "/TMP/nLL" ); + _nSS = bookCounter( "/TMP/nSS" ); + _nLS = bookCounter( "/TMP/nLS" ); + } + + 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"); + // total hadronic and muonic cross sections + map nCount; + int ntotal(0); + foreach (const Particle& p, fs.particles()) { + nCount[p.pdgId()] += 1; + ++ntotal; + } + + // find the Lambdas and Sigmas + const FinalState& ufs = apply(event, "UFS"); + for(unsigned int ix=0;ix nRes = nCount; + int ncount = ntotal; + findChildren(p1,nRes,ncount); + for(unsigned int iy=ix+1;iy nRes2 = nRes; + int ncount2 = ncount; + findChildren(p2,nRes2,ncount2); + if(ncount2!=0) continue; + matched=true; + for(auto const & val : nRes2) { + if(val.second!=0) { + matched = false; + break; + } + } + if(matched) { + if(abs(p1.pdgId())==3122 && abs(p2.pdgId())==3122) + _nLL->fill(event.weight()); + else if(abs(p1.pdgId())==3212 && abs(p2.pdgId())==3212) + _nSS->fill(event.weight()); + else + _nLS->fill(event.weight()); + break; + } + } + if(matched) break; + } + } + + /// Normalise histograms etc., after the run + void finalize() { + double fact = crossSection()/ sumOfWeights() /picobarn; + for(unsigned int ix=1;ix<4;++ix) { + double sigma,error; + if(ix==1) { + sigma = _nLL->val()*fact; + error = _nLL->err()*fact; + } + else if(ix==2) { + sigma = _nSS->val()*fact; + error = _nSS->err()*fact; + } + else { + sigma = _nLS->val()*fact; + error = _nLS->err()*fact; + } + 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 _nLL,_nSS,_nLS; + //@} + + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(BABAR_2007_I760730); + + +} diff --git a/analyses/pluginBABAR/BABAR_2007_I760730.info b/analyses/pluginBABAR/BABAR_2007_I760730.info new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2007_I760730.info @@ -0,0 +1,35 @@ +Name: BABAR_2007_I760730 +Year: 2007 +Summary: Cross section for $e^+e^-\to\Lambda\bar{\Lambda}$, $\Sigma^0\bar{\Sigma}^0$ and $\Lambda\bar{\Sigma}^0$ between threshold and 3 GeV +Experiment: BABAR +Collider: PEP-II +InspireID: 760730 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - Phys.Rev. D76 (2007) 092006 +RunInfo: e+e- to hadrons +NeedCrossSection: yes +Beams: [e+,e-] +Description: + 'Measurement of the cross section for $e^+e^-\to\Lambda\bar{\Lambda}$ between threshold and 3 GeV.' +Keywords: +BibKey: Aubert:2007uf +BibTeX: '@article{Aubert:2007uf, + author = "Aubert, Bernard and others", + title = "{Study of $e^{+} e^{-} \to \Lambda \bar{\Lambda}$, + $\Lambda \bar{\Sigma}^0$, $\Sigma^0 \bar{\Sigma}^0$ using + initial state radiation with BABAR}", + collaboration = "BaBar", + journal = "Phys. Rev.", + volume = "D76", + year = "2007", + pages = "092006", + doi = "10.1103/PhysRevD.76.092006", + eprint = "0709.1988", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + reportNumber = "SLAC-PUB-12806, BABAR-PUB-07-054", + SLACcitation = "%%CITATION = ARXIV:0709.1988;%%" +}' diff --git a/analyses/pluginBABAR/BABAR_2007_I760730.plot b/analyses/pluginBABAR/BABAR_2007_I760730.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2007_I760730.plot @@ -0,0 +1,18 @@ +BEGIN PLOT /BABAR_2007_I760730/d01-x01-y01 +Title=$\sigma(e^+e^-\to \Lambda\bar{\Lambda})$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to \Lambda\bar{\Lambda})$/pb +ConnectGaps=1 +END PLOT +BEGIN PLOT /BABAR_2007_I760730/d02-x01-y01 +Title=$\sigma(e^+e^-\to \Sigma\bar{\Sigma}^0)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to \Sigma\bar{\Sigma}^0)$/pb +ConnectGaps=1 +END PLOT +BEGIN PLOT /BABAR_2007_I760730/d03-x01-y01 +Title=$\sigma(e^+e^-\to \Lambda\bar{\Sigma}^0+\text{c.c})$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to \Lambda\bar{\Sigma}^0+\text{c.c.})$/pb +ConnectGaps=1 +END PLOT diff --git a/analyses/pluginBABAR/BABAR_2007_I760730.yoda b/analyses/pluginBABAR/BABAR_2007_I760730.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2007_I760730.yoda @@ -0,0 +1,41 @@ +BEGIN YODA_SCATTER2D_V2 /REF/BABAR_2007_I760730/d01-x01-y01 +IsRef: 1 +Path: /REF/BABAR_2007_I760730/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +2.25000 0.02000 0.02000 2.04000e+02 6.39062e+01 6.57875e+01 +2.28500 0.01500 0.01500 2.02000e+02 5.12640e+01 5.31413e+01 +2.32500 0.02500 0.02500 1.55000e+02 2.77308e+01 3.04631e+01 +2.37500 0.02500 0.02500 1.76000e+02 3.44384e+01 3.44384e+01 +2.42500 0.02500 0.02500 9.00000e+01 3.10483e+01 3.20156e+01 +2.47500 0.02500 0.02500 1.04000e+02 2.32594e+01 2.60000e+01 +2.55000 0.05000 0.05000 4.60000e+01 9.84886e+00 1.17047e+01 +2.65000 0.05000 0.05000 2.50000e+01 1.31529e+01 1.21655e+01 +2.75000 0.05000 0.05000 1.00000e+01 8.06226e+00 9.05539e+00 +2.90000 0.10000 0.10000 2.10000e+00 3.50571e+00 3.20624e+00 +END YODA_SCATTER2D_V2 +BEGIN YODA_SCATTER2D_V2 /REF/BABAR_2007_I760730/d02-x01-y01 +IsRef: 1 +Path: /REF/BABAR_2007_I760730/d02-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +2.49250 0.10750 0.10750 3.00000e+01 1.33417e+01 1.33417e+01 +2.70000 0.10000 0.10000 1.70000e+01 1.01980e+01 8.24621e+00 +2.90000 0.10000 0.10000 3.40000e+00 7.81025e+00 8.50941e+00 +END YODA_SCATTER2D_V2 +BEGIN YODA_SCATTER2D_V2 /REF/BABAR_2007_I760730/d03-x01-y01 +IsRef: 1 +Path: /REF/BABAR_2007_I760730/d03-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +2.35400 0.04600 0.04600 4.70000e+01 2.15870e+01 2.35372e+01 +2.45000 0.05000 0.05000 3.20000e+01 1.84391e+01 1.94165e+01 +2.60000 0.10000 0.10000 7.10000e+00 7.53260e+00 8.22982e+00 +2.80000 0.10000 0.10000 2.90000e+00 3.71214e+00 4.31045e+00 +END YODA_SCATTER2D_V2 diff --git a/analyses/pluginBABAR/BABAR_2009_I829441.cc b/analyses/pluginBABAR/BABAR_2009_I829441.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2009_I829441.cc @@ -0,0 +1,80 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/Projections/FinalState.hh" + +namespace Rivet { + + + /// @brief Add a short analysis description here + class BABAR_2009_I829441 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(BABAR_2009_I829441); + + + /// @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(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(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 _npion; + //@} + + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(BABAR_2009_I829441); + + +} diff --git a/analyses/pluginBABAR/BABAR_2009_I829441.info b/analyses/pluginBABAR/BABAR_2009_I829441.info new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2009_I829441.info @@ -0,0 +1,38 @@ +Name: BABAR_2009_I829441 +Year: 2009 +Summary: Cross section for $e^+e^-\to \pi^+\pi^-$ below 3 GeV +Experiment: BABAR +Collider: PEP-II +InspireID: 829441 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - Phys.Rev.Lett. 103 (2009) 231801 + - Phys.Rev. D86 (2012) 032013 +RunInfo: e+e- to hadrons +NeedCrossSection: yes +Beams: [e+,e-] +Luminosity_fb: +Description: +Description: + 'Measurement of the cross section for $e^+e^-\to \pi^+\pi^-$ at energies near the $\rho$ mass.' +Keywords: [] +BibKey: Aubert:2009ad +BibTeX: '@article{Aubert:2009ad, + author = "Aubert, Bernard and others", + title = "{Precise measurement of the e+ e- ---> pi+ pi- (gamma) + cross section with the Initial State Radiation method at + BABAR}", + collaboration = "BaBar", + journal = "Phys. Rev. Lett.", + volume = "103", + year = "2009", + pages = "231801", + doi = "10.1103/PhysRevLett.103.231801", + eprint = "0908.3589", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + reportNumber = "SLAC-PUB-13754, BABAR-PUB-09-027", + SLACcitation = "%%CITATION = ARXIV:0908.3589;%%" +}' diff --git a/analyses/pluginBABAR/BABAR_2009_I829441.plot b/analyses/pluginBABAR/BABAR_2009_I829441.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2009_I829441.plot @@ -0,0 +1,6 @@ +BEGIN PLOT /BABAR_2009_I829441/d01-x01-y01 +Title=$\sigma(e^+e^-\to \pi^+\pi^-)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to \pi^+\pi^-)$/nb +ConnectGaps=1 +END PLOT diff --git a/analyses/pluginBABAR/BABAR_2009_I829441.yoda b/analyses/pluginBABAR/BABAR_2009_I829441.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2009_I829441.yoda @@ -0,0 +1,345 @@ +BEGIN YODA_SCATTER2D_V2 /REF/BABAR_2009_I829441/d01-x01-y01 +IsRef: 1 +Path: /REF/BABAR_2009_I829441/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +0.30500 0.00500 0.00500 2.54904e+01 2.69943e+00 2.69943e+00 +0.31500 0.00500 0.00500 3.54801e+01 2.91464e+00 2.91464e+00 +0.32500 0.00500 0.00500 4.54858e+01 3.04669e+00 3.04669e+00 +0.33500 0.00500 0.00500 5.17825e+01 3.13355e+00 3.13355e+00 +0.34500 0.00500 0.00500 6.44156e+01 3.49953e+00 3.49953e+00 +0.35500 0.00500 0.00500 7.52105e+01 3.86247e+00 3.86247e+00 +0.36500 0.00500 0.00500 8.01998e+01 3.87376e+00 3.87376e+00 +0.37500 0.00500 0.00500 8.63805e+01 4.02805e+00 4.02805e+00 +0.38500 0.00500 0.00500 1.02438e+02 4.31976e+00 4.31976e+00 +0.39500 0.00500 0.00500 1.04710e+02 3.96803e+00 3.96803e+00 +0.40500 0.00500 0.00500 1.18252e+02 4.09898e+00 4.09898e+00 +0.41500 0.00500 0.00500 1.18000e+02 3.90209e+00 3.90209e+00 +0.42500 0.00500 0.00500 1.27675e+02 4.17054e+00 4.17054e+00 +0.43500 0.00500 0.00500 1.29278e+02 3.83703e+00 3.83703e+00 +0.44500 0.00500 0.00500 1.37056e+02 4.08884e+00 4.08884e+00 +0.45500 0.00500 0.00500 1.38780e+02 3.80708e+00 3.80708e+00 +0.46500 0.00500 0.00500 1.49933e+02 4.05805e+00 4.05805e+00 +0.47500 0.00500 0.00500 1.61274e+02 4.29024e+00 4.29024e+00 +0.48500 0.00500 0.00500 1.67659e+02 4.29975e+00 4.29975e+00 +0.49500 0.00500 0.00500 1.70777e+02 4.09567e+00 4.09567e+00 +0.50100 0.00100 0.00100 1.70490e+02 4.94587e+00 4.94587e+00 +0.50300 0.00100 0.00100 1.73419e+02 4.53398e+00 4.53398e+00 +0.50500 0.00100 0.00100 1.77172e+02 4.50925e+00 4.50925e+00 +0.50700 0.00100 0.00100 1.75923e+02 4.79673e+00 4.79673e+00 +0.50900 0.00100 0.00100 1.76444e+02 5.32704e+00 5.32704e+00 +0.51100 0.00100 0.00100 1.78940e+02 5.10733e+00 5.10733e+00 +0.51300 0.00100 0.00100 1.88560e+02 5.72383e+00 5.72383e+00 +0.51500 0.00100 0.00100 1.85026e+02 5.34796e+00 5.34796e+00 +0.51700 0.00100 0.00100 1.87884e+02 4.52075e+00 4.52075e+00 +0.51900 0.00100 0.00100 1.91173e+02 4.83832e+00 4.83832e+00 +0.52100 0.00100 0.00100 1.86696e+02 4.97632e+00 4.97632e+00 +0.52300 0.00100 0.00100 1.93177e+02 5.58947e+00 5.58947e+00 +0.52500 0.00100 0.00100 1.99509e+02 5.91304e+00 5.91304e+00 +0.52700 0.00100 0.00100 1.97814e+02 4.73803e+00 4.73803e+00 +0.52900 0.00100 0.00100 2.02529e+02 4.70092e+00 4.70092e+00 +0.53100 0.00100 0.00100 2.01567e+02 5.89086e+00 5.89086e+00 +0.53300 0.00100 0.00100 2.05391e+02 5.87140e+00 5.87140e+00 +0.53500 0.00100 0.00100 2.12441e+02 5.04889e+00 5.04889e+00 +0.53700 0.00100 0.00100 2.17834e+02 4.88096e+00 4.88096e+00 +0.53900 0.00100 0.00100 2.16395e+02 5.23306e+00 5.23306e+00 +0.54100 0.00100 0.00100 2.22313e+02 6.02197e+00 6.02197e+00 +0.54300 0.00100 0.00100 2.27740e+02 5.73655e+00 5.73655e+00 +0.54500 0.00100 0.00100 2.25003e+02 5.73778e+00 5.73778e+00 +0.54700 0.00100 0.00100 2.34542e+02 6.23250e+00 6.23250e+00 +0.54900 0.00100 0.00100 2.35788e+02 5.23052e+00 5.23052e+00 +0.55100 0.00100 0.00100 2.36245e+02 5.65545e+00 5.65545e+00 +0.55300 0.00100 0.00100 2.35052e+02 6.10188e+00 6.10188e+00 +0.55500 0.00100 0.00100 2.44258e+02 6.03505e+00 6.03505e+00 +0.55700 0.00100 0.00100 2.41296e+02 6.26508e+00 6.26508e+00 +0.55900 0.00100 0.00100 2.38944e+02 4.88960e+00 4.88960e+00 +0.56100 0.00100 0.00100 2.47561e+02 5.95113e+00 5.95113e+00 +0.56300 0.00100 0.00100 2.52729e+02 6.22876e+00 6.22876e+00 +0.56500 0.00100 0.00100 2.54064e+02 6.29356e+00 6.29356e+00 +0.56700 0.00100 0.00100 2.51711e+02 5.89993e+00 5.89993e+00 +0.56900 0.00100 0.00100 2.63511e+02 6.23737e+00 6.23737e+00 +0.57100 0.00100 0.00100 2.58408e+02 5.95967e+00 5.95967e+00 +0.57300 0.00100 0.00100 2.57723e+02 6.36086e+00 6.36086e+00 +0.57500 0.00100 0.00100 2.66283e+02 5.16458e+00 5.16458e+00 +0.57700 0.00100 0.00100 2.76459e+02 6.53556e+00 6.53556e+00 +0.57900 0.00100 0.00100 2.68516e+02 6.16405e+00 6.16405e+00 +0.58100 0.00100 0.00100 2.78692e+02 6.61913e+00 6.61913e+00 +0.58300 0.00100 0.00100 2.79647e+02 5.81224e+00 5.81224e+00 +0.58500 0.00100 0.00100 2.87129e+02 6.07523e+00 6.07523e+00 +0.58700 0.00100 0.00100 2.90007e+02 5.74166e+00 5.74166e+00 +0.58900 0.00100 0.00100 2.89457e+02 6.42293e+00 6.42293e+00 +0.59100 0.00100 0.00100 2.88048e+02 5.95134e+00 5.95134e+00 +0.59300 0.00100 0.00100 3.01113e+02 5.52231e+00 5.52231e+00 +0.59500 0.00100 0.00100 3.09435e+02 6.95643e+00 6.95643e+00 +0.59700 0.00100 0.00100 3.21707e+02 6.80183e+00 6.80183e+00 +0.59900 0.00100 0.00100 3.30973e+02 7.60966e+00 7.60966e+00 +0.60100 0.00100 0.00100 3.25413e+02 6.06285e+00 6.06285e+00 +0.60300 0.00100 0.00100 3.32822e+02 5.21301e+00 5.21301e+00 +0.60500 0.00100 0.00100 3.36632e+02 5.88944e+00 5.88944e+00 +0.60700 0.00100 0.00100 3.41428e+02 6.67362e+00 6.67362e+00 +0.60900 0.00100 0.00100 3.44518e+02 6.71287e+00 6.71287e+00 +0.61100 0.00100 0.00100 3.49659e+02 5.70493e+00 5.70493e+00 +0.61300 0.00100 0.00100 3.53204e+02 5.30912e+00 5.30912e+00 +0.61500 0.00100 0.00100 3.57844e+02 6.82366e+00 6.82366e+00 +0.61700 0.00100 0.00100 3.66571e+02 7.12967e+00 7.12967e+00 +0.61900 0.00100 0.00100 3.73152e+02 5.35323e+00 5.35323e+00 +0.62100 0.00100 0.00100 3.77674e+02 6.75055e+00 6.75055e+00 +0.62300 0.00100 0.00100 3.92495e+02 6.66409e+00 6.66409e+00 +0.62500 0.00100 0.00100 3.89935e+02 6.72768e+00 6.72768e+00 +0.62700 0.00100 0.00100 4.03266e+02 7.19893e+00 7.19893e+00 +0.62900 0.00100 0.00100 4.09291e+02 6.64353e+00 6.64353e+00 +0.63100 0.00100 0.00100 4.32872e+02 7.16007e+00 7.16007e+00 +0.63300 0.00100 0.00100 4.32317e+02 6.74057e+00 6.74057e+00 +0.63500 0.00100 0.00100 4.41693e+02 7.32202e+00 7.32202e+00 +0.63700 0.00100 0.00100 4.47307e+02 7.55775e+00 7.55775e+00 +0.63900 0.00100 0.00100 4.40543e+02 7.26787e+00 7.26787e+00 +0.64100 0.00100 0.00100 4.62622e+02 6.96407e+00 6.96407e+00 +0.64300 0.00100 0.00100 4.69842e+02 7.75651e+00 7.75651e+00 +0.64500 0.00100 0.00100 4.73894e+02 6.26555e+00 6.26555e+00 +0.64700 0.00100 0.00100 4.81965e+02 7.71686e+00 7.71686e+00 +0.64900 0.00100 0.00100 5.06772e+02 7.42557e+00 7.42557e+00 +0.65100 0.00100 0.00100 5.15120e+02 7.32130e+00 7.32130e+00 +0.65300 0.00100 0.00100 5.21883e+02 7.97169e+00 7.97169e+00 +0.65500 0.00100 0.00100 5.27421e+02 8.10339e+00 8.10339e+00 +0.65700 0.00100 0.00100 5.38889e+02 8.17683e+00 8.17683e+00 +0.65900 0.00100 0.00100 5.49541e+02 8.30160e+00 8.30160e+00 +0.66100 0.00100 0.00100 5.56211e+02 7.72668e+00 7.72668e+00 +0.66300 0.00100 0.00100 5.70704e+02 8.09063e+00 8.09063e+00 +0.66500 0.00100 0.00100 5.86913e+02 7.50127e+00 7.50127e+00 +0.66700 0.00100 0.00100 5.89595e+02 7.81293e+00 7.81293e+00 +0.66900 0.00100 0.00100 6.02804e+02 8.54057e+00 8.54057e+00 +0.67100 0.00100 0.00100 6.19429e+02 8.63669e+00 8.63669e+00 +0.67300 0.00100 0.00100 6.46475e+02 8.35179e+00 8.35179e+00 +0.67500 0.00100 0.00100 6.53007e+02 9.22731e+00 9.22731e+00 +0.67700 0.00100 0.00100 6.70383e+02 9.25269e+00 9.25269e+00 +0.67900 0.00100 0.00100 6.88002e+02 9.06100e+00 9.06100e+00 +0.68100 0.00100 0.00100 7.04235e+02 9.05242e+00 9.05242e+00 +0.68300 0.00100 0.00100 7.14122e+02 9.62278e+00 9.62278e+00 +0.68500 0.00100 0.00100 7.38793e+02 9.93198e+00 9.93198e+00 +0.68700 0.00100 0.00100 7.56886e+02 9.10639e+00 9.10639e+00 +0.68900 0.00100 0.00100 7.62360e+02 9.25970e+00 9.25970e+00 +0.69100 0.00100 0.00100 7.98399e+02 1.04788e+01 1.04788e+01 +0.69300 0.00100 0.00100 8.02306e+02 8.97791e+00 8.97791e+00 +0.69500 0.00100 0.00100 8.30193e+02 9.59223e+00 9.59223e+00 +0.69700 0.00100 0.00100 8.52331e+02 1.04668e+01 1.04668e+01 +0.69900 0.00100 0.00100 8.58301e+02 1.04834e+01 1.04834e+01 +0.70100 0.00100 0.00100 8.75293e+02 1.07817e+01 1.07817e+01 +0.70300 0.00100 0.00100 8.82796e+02 1.10894e+01 1.10894e+01 +0.70500 0.00100 0.00100 9.05022e+02 1.04656e+01 1.04656e+01 +0.70700 0.00100 0.00100 9.44689e+02 1.14607e+01 1.14607e+01 +0.70900 0.00100 0.00100 9.43541e+02 1.10681e+01 1.10681e+01 +0.71100 0.00100 0.00100 9.74289e+02 1.17619e+01 1.17619e+01 +0.71300 0.00100 0.00100 1.00055e+03 1.17577e+01 1.17577e+01 +0.71500 0.00100 0.00100 1.01121e+03 1.22018e+01 1.22018e+01 +0.71700 0.00100 0.00100 1.04300e+03 1.18594e+01 1.18594e+01 +0.71900 0.00100 0.00100 1.06149e+03 1.15454e+01 1.15454e+01 +0.72100 0.00100 0.00100 1.07526e+03 1.27249e+01 1.27249e+01 +0.72300 0.00100 0.00100 1.10393e+03 1.29544e+01 1.29544e+01 +0.72500 0.00100 0.00100 1.10533e+03 1.26579e+01 1.26579e+01 +0.72700 0.00100 0.00100 1.13328e+03 1.16202e+01 1.16202e+01 +0.72900 0.00100 0.00100 1.14954e+03 1.31084e+01 1.31084e+01 +0.73100 0.00100 0.00100 1.15838e+03 1.24490e+01 1.24490e+01 +0.73300 0.00100 0.00100 1.19523e+03 1.25723e+01 1.25723e+01 +0.73500 0.00100 0.00100 1.20812e+03 1.36549e+01 1.36549e+01 +0.73700 0.00100 0.00100 1.23571e+03 1.23235e+01 1.23235e+01 +0.73900 0.00100 0.00100 1.26629e+03 1.42119e+01 1.42119e+01 +0.74100 0.00100 0.00100 1.26664e+03 1.24516e+01 1.24516e+01 +0.74300 0.00100 0.00100 1.27859e+03 1.42056e+01 1.42056e+01 +0.74500 0.00100 0.00100 1.28014e+03 1.25097e+01 1.25097e+01 +0.74700 0.00100 0.00100 1.29324e+03 1.43521e+01 1.43521e+01 +0.74900 0.00100 0.00100 1.31638e+03 1.42986e+01 1.42986e+01 +0.75100 0.00100 0.00100 1.31756e+03 1.30731e+01 1.30731e+01 +0.75300 0.00100 0.00100 1.32269e+03 1.36937e+01 1.36937e+01 +0.75500 0.00100 0.00100 1.33774e+03 1.37788e+01 1.37788e+01 +0.75700 0.00100 0.00100 1.31038e+03 1.44103e+01 1.44103e+01 +0.75900 0.00100 0.00100 1.31385e+03 1.32756e+01 1.32756e+01 +0.76100 0.00100 0.00100 1.32091e+03 1.38624e+01 1.38624e+01 +0.76300 0.00100 0.00100 1.32731e+03 1.33224e+01 1.33224e+01 +0.76500 0.00100 0.00100 1.31828e+03 1.37332e+01 1.37332e+01 +0.76700 0.00100 0.00100 1.31243e+03 1.29789e+01 1.29789e+01 +0.76900 0.00100 0.00100 1.32067e+03 1.45973e+01 1.45973e+01 +0.77100 0.00100 0.00100 1.31279e+03 1.44817e+01 1.44817e+01 +0.77300 0.00100 0.00100 1.32546e+03 1.44727e+01 1.44727e+01 +0.77500 0.00100 0.00100 1.33186e+03 1.47100e+01 1.47100e+01 +0.77700 0.00100 0.00100 1.33377e+03 1.48083e+01 1.48083e+01 +0.77900 0.00100 0.00100 1.31142e+03 1.29713e+01 1.29713e+01 +0.78100 0.00100 0.00100 1.19470e+03 1.35010e+01 1.35010e+01 +0.78300 0.00100 0.00100 9.85965e+02 1.14479e+01 1.14479e+01 +0.78500 0.00100 0.00100 8.66014e+02 1.01592e+01 1.01592e+01 +0.78700 0.00100 0.00100 8.31708e+02 1.00297e+01 1.00297e+01 +0.78900 0.00100 0.00100 8.31511e+02 9.75744e+00 9.75744e+00 +0.79100 0.00100 0.00100 8.22419e+02 8.80198e+00 8.80198e+00 +0.79300 0.00100 0.00100 8.22363e+02 8.95940e+00 8.95940e+00 +0.79500 0.00100 0.00100 8.29933e+02 1.01655e+01 1.01655e+01 +0.79700 0.00100 0.00100 8.14409e+02 9.85160e+00 9.85160e+00 +0.79900 0.00100 0.00100 7.86562e+02 8.90542e+00 8.90542e+00 +0.80100 0.00100 0.00100 7.75412e+02 8.71936e+00 8.71936e+00 +0.80300 0.00100 0.00100 7.61926e+02 9.46256e+00 9.46256e+00 +0.80500 0.00100 0.00100 7.47172e+02 8.48462e+00 8.48462e+00 +0.80700 0.00100 0.00100 7.21769e+02 7.90131e+00 7.90131e+00 +0.80900 0.00100 0.00100 7.00670e+02 8.51080e+00 8.51080e+00 +0.81100 0.00100 0.00100 6.77105e+02 8.68613e+00 8.68613e+00 +0.81300 0.00100 0.00100 6.57361e+02 8.32101e+00 8.32101e+00 +0.81500 0.00100 0.00100 6.51412e+02 8.35477e+00 8.35477e+00 +0.81700 0.00100 0.00100 6.24082e+02 8.29506e+00 8.29506e+00 +0.81900 0.00100 0.00100 5.99225e+02 7.26060e+00 7.26060e+00 +0.82100 0.00100 0.00100 5.94326e+02 7.00189e+00 7.00189e+00 +0.82300 0.00100 0.00100 5.80226e+02 7.80643e+00 7.80643e+00 +0.82500 0.00100 0.00100 5.47810e+02 7.62455e+00 7.62455e+00 +0.82700 0.00100 0.00100 5.44839e+02 7.34869e+00 7.34869e+00 +0.82900 0.00100 0.00100 5.29460e+02 6.88818e+00 6.88818e+00 +0.83100 0.00100 0.00100 5.06543e+02 7.08034e+00 7.08034e+00 +0.83300 0.00100 0.00100 4.92976e+02 6.91965e+00 6.91965e+00 +0.83500 0.00100 0.00100 4.76862e+02 6.69314e+00 6.69314e+00 +0.83700 0.00100 0.00100 4.65248e+02 6.80354e+00 6.80354e+00 +0.83900 0.00100 0.00100 4.51144e+02 5.84514e+00 5.84514e+00 +0.84100 0.00100 0.00100 4.35738e+02 6.29826e+00 6.29826e+00 +0.84300 0.00100 0.00100 4.20266e+02 6.34874e+00 6.34874e+00 +0.84500 0.00100 0.00100 4.10063e+02 5.43764e+00 5.43764e+00 +0.84700 0.00100 0.00100 4.01241e+02 6.16992e+00 6.16992e+00 +0.84900 0.00100 0.00100 3.78963e+02 5.07024e+00 5.07024e+00 +0.85100 0.00100 0.00100 3.72188e+02 5.54278e+00 5.54278e+00 +0.85300 0.00100 0.00100 3.65819e+02 5.73991e+00 5.73991e+00 +0.85500 0.00100 0.00100 3.38705e+02 5.44493e+00 5.44493e+00 +0.85700 0.00100 0.00100 3.32300e+02 4.78480e+00 4.78480e+00 +0.85900 0.00100 0.00100 3.29412e+02 4.76065e+00 4.76065e+00 +0.86100 0.00100 0.00100 3.15291e+02 4.84528e+00 4.84528e+00 +0.86300 0.00100 0.00100 3.15040e+02 5.25050e+00 5.25050e+00 +0.86500 0.00100 0.00100 3.00625e+02 5.13391e+00 5.13391e+00 +0.86700 0.00100 0.00100 2.95889e+02 4.33473e+00 4.33473e+00 +0.86900 0.00100 0.00100 2.86377e+02 4.71602e+00 4.71602e+00 +0.87100 0.00100 0.00100 2.73246e+02 4.84049e+00 4.84049e+00 +0.87300 0.00100 0.00100 2.62975e+02 4.61824e+00 4.61824e+00 +0.87500 0.00100 0.00100 2.62874e+02 4.05910e+00 4.05910e+00 +0.87700 0.00100 0.00100 2.52677e+02 3.81778e+00 3.81778e+00 +0.87900 0.00100 0.00100 2.42882e+02 4.54679e+00 4.54679e+00 +0.88100 0.00100 0.00100 2.36673e+02 3.65542e+00 3.65542e+00 +0.88300 0.00100 0.00100 2.34518e+02 4.43607e+00 4.43607e+00 +0.88500 0.00100 0.00100 2.28634e+02 3.51008e+00 3.51008e+00 +0.88700 0.00100 0.00100 2.24883e+02 3.69221e+00 3.69221e+00 +0.88900 0.00100 0.00100 2.10195e+02 3.33349e+00 3.33349e+00 +0.89100 0.00100 0.00100 2.11491e+02 4.01305e+00 4.01305e+00 +0.89300 0.00100 0.00100 2.03657e+02 4.05319e+00 4.05319e+00 +0.89500 0.00100 0.00100 1.97638e+02 3.59671e+00 3.59671e+00 +0.89700 0.00100 0.00100 1.99534e+02 3.84660e+00 3.84660e+00 +0.89900 0.00100 0.00100 1.92686e+02 3.29938e+00 3.29938e+00 +0.90100 0.00100 0.00100 1.86710e+02 3.30402e+00 3.30402e+00 +0.90300 0.00100 0.00100 1.83346e+02 3.86473e+00 3.86473e+00 +0.90500 0.00100 0.00100 1.75526e+02 3.42085e+00 3.42085e+00 +0.90700 0.00100 0.00100 1.68741e+02 3.02536e+00 3.02536e+00 +0.90900 0.00100 0.00100 1.60017e+02 3.52612e+00 3.52612e+00 +0.91100 0.00100 0.00100 1.59971e+02 3.33427e+00 3.33427e+00 +0.91300 0.00100 0.00100 1.55649e+02 3.44126e+00 3.44126e+00 +0.91500 0.00100 0.00100 1.48082e+02 3.10591e+00 3.10591e+00 +0.91700 0.00100 0.00100 1.47357e+02 2.76511e+00 2.76511e+00 +0.91900 0.00100 0.00100 1.45097e+02 3.42667e+00 3.42667e+00 +0.92100 0.00100 0.00100 1.35646e+02 2.57948e+00 2.57948e+00 +0.92300 0.00100 0.00100 1.39373e+02 3.22257e+00 3.22257e+00 +0.92500 0.00100 0.00100 1.32803e+02 3.18767e+00 3.18767e+00 +0.92700 0.00100 0.00100 1.28585e+02 3.14486e+00 3.14486e+00 +0.92900 0.00100 0.00100 1.27816e+02 2.50977e+00 2.50977e+00 +0.93100 0.00100 0.00100 1.23242e+02 3.00352e+00 3.00352e+00 +0.93300 0.00100 0.00100 1.22269e+02 3.02470e+00 3.02470e+00 +0.93500 0.00100 0.00100 1.15368e+02 2.84473e+00 2.84473e+00 +0.93700 0.00100 0.00100 1.12606e+02 2.69351e+00 2.69351e+00 +0.93900 0.00100 0.00100 1.09421e+02 2.43689e+00 2.43689e+00 +0.94100 0.00100 0.00100 1.11148e+02 2.33278e+00 2.33278e+00 +0.94300 0.00100 0.00100 1.07676e+02 2.63038e+00 2.63038e+00 +0.94500 0.00100 0.00100 1.03334e+02 2.16834e+00 2.16834e+00 +0.94700 0.00100 0.00100 1.05928e+02 2.52356e+00 2.52356e+00 +0.94900 0.00100 0.00100 9.92949e+01 2.58985e+00 2.58985e+00 +0.95100 0.00100 0.00100 9.98899e+01 2.63741e+00 2.63741e+00 +0.95300 0.00100 0.00100 9.89337e+01 2.60154e+00 2.60154e+00 +0.95500 0.00100 0.00100 9.46227e+01 2.54042e+00 2.54042e+00 +0.95700 0.00100 0.00100 9.49455e+01 2.57963e+00 2.57963e+00 +0.95900 0.00100 0.00100 8.92567e+01 2.13248e+00 2.13248e+00 +0.96100 0.00100 0.00100 8.84068e+01 2.54816e+00 2.54816e+00 +0.96300 0.00100 0.00100 8.63153e+01 1.95674e+00 1.95674e+00 +0.96500 0.00100 0.00100 8.23814e+01 2.09390e+00 2.09390e+00 +0.96700 0.00100 0.00100 8.08113e+01 2.13539e+00 2.13539e+00 +0.96900 0.00100 0.00100 7.95135e+01 2.32598e+00 2.32598e+00 +0.97100 0.00100 0.00100 7.91463e+01 2.38683e+00 2.38683e+00 +0.97300 0.00100 0.00100 8.20929e+01 2.39002e+00 2.39002e+00 +0.97500 0.00100 0.00100 7.61146e+01 1.74979e+00 1.74979e+00 +0.97700 0.00100 0.00100 7.31396e+01 2.07533e+00 2.07533e+00 +0.97900 0.00100 0.00100 7.39941e+01 2.12724e+00 2.12724e+00 +0.98100 0.00100 0.00100 7.07019e+01 1.96335e+00 1.96335e+00 +0.98300 0.00100 0.00100 7.19296e+01 1.81812e+00 1.81812e+00 +0.98500 0.00100 0.00100 6.71707e+01 2.18207e+00 2.18207e+00 +0.98700 0.00100 0.00100 6.88309e+01 1.97739e+00 1.97739e+00 +0.98900 0.00100 0.00100 6.43775e+01 1.90782e+00 1.90782e+00 +0.99100 0.00100 0.00100 6.65032e+01 2.15347e+00 2.15347e+00 +0.99300 0.00100 0.00100 6.36086e+01 1.71685e+00 1.71685e+00 +0.99500 0.00100 0.00100 6.34289e+01 1.69409e+00 1.69409e+00 +0.99700 0.00100 0.00100 6.18220e+01 1.81624e+00 1.81624e+00 +0.99900 0.00100 0.00100 6.13907e+01 2.06360e+00 2.06360e+00 +1.00500 0.00500 0.00500 5.94632e+01 1.57942e+00 1.57942e+00 +1.01500 0.00500 0.00500 5.54884e+01 1.48186e+00 1.48186e+00 +1.02500 0.00500 0.00500 4.26740e+01 1.16288e+00 1.16288e+00 +1.03500 0.00500 0.00500 4.27609e+01 1.29947e+00 1.29947e+00 +1.04500 0.00500 0.00500 3.97599e+01 1.19043e+00 1.19043e+00 +1.05500 0.00500 0.00500 3.67208e+01 1.08106e+00 1.08106e+00 +1.06500 0.00500 0.00500 3.48831e+01 1.15828e+00 1.15828e+00 +1.07500 0.00500 0.00500 3.33657e+01 1.15134e+00 1.15134e+00 +1.08500 0.00500 0.00500 3.00331e+01 1.00221e+00 1.00221e+00 +1.09500 0.00500 0.00500 2.96114e+01 9.78244e-01 9.78244e-01 +1.10500 0.00500 0.00500 2.63388e+01 9.45589e-01 9.45589e-01 +1.11500 0.00500 0.00500 2.52242e+01 9.59383e-01 9.59383e-01 +1.12500 0.00500 0.00500 2.28780e+01 9.06155e-01 9.06155e-01 +1.13500 0.00500 0.00500 2.12351e+01 7.96671e-01 7.96671e-01 +1.14500 0.00500 0.00500 2.10231e+01 8.44067e-01 8.44067e-01 +1.15500 0.00500 0.00500 1.97813e+01 8.24146e-01 8.24146e-01 +1.16500 0.00500 0.00500 1.76864e+01 7.76961e-01 7.76961e-01 +1.17500 0.00500 0.00500 1.63228e+01 6.77793e-01 6.77793e-01 +1.18500 0.00500 0.00500 1.58713e+01 7.34172e-01 7.34172e-01 +1.19500 0.00500 0.00500 1.46102e+01 6.71846e-01 6.71846e-01 +1.20500 0.00500 0.00500 1.32862e+01 6.18357e-01 6.18357e-01 +1.21500 0.00500 0.00500 1.24045e+01 6.28415e-01 6.28415e-01 +1.22500 0.00500 0.00500 1.13420e+01 5.59182e-01 5.59182e-01 +1.23500 0.00500 0.00500 1.08011e+01 5.35139e-01 5.35139e-01 +1.24500 0.00500 0.00500 1.08246e+01 5.99093e-01 5.99093e-01 +1.25500 0.00500 0.00500 9.84104e+00 5.21959e-01 5.21959e-01 +1.26500 0.00500 0.00500 9.83737e+00 5.68676e-01 5.68676e-01 +1.27500 0.00500 0.00500 8.58472e+00 5.13665e-01 5.13665e-01 +1.28500 0.00500 0.00500 8.06022e+00 5.04002e-01 5.04002e-01 +1.29500 0.00500 0.00500 7.27264e+00 4.54801e-01 4.54801e-01 +1.30500 0.00500 0.00500 6.81295e+00 4.47716e-01 4.47716e-01 +1.31500 0.00500 0.00500 6.74482e+00 4.30431e-01 4.30431e-01 +1.32500 0.00500 0.00500 5.62010e+00 3.71815e-01 3.71815e-01 +1.33500 0.00500 0.00500 5.68195e+00 3.90390e-01 3.90390e-01 +1.34500 0.00500 0.00500 4.80148e+00 3.75290e-01 3.75290e-01 +1.35500 0.00500 0.00500 4.92940e+00 3.54625e-01 3.54625e-01 +1.36500 0.00500 0.00500 4.22120e+00 3.56537e-01 3.56537e-01 +1.37500 0.00500 0.00500 3.86059e+00 3.14442e-01 3.14442e-01 +1.38500 0.00500 0.00500 2.79557e+00 2.98523e-01 2.98523e-01 +1.39500 0.00500 0.00500 2.37278e+00 2.32561e-01 2.32561e-01 +1.42500 0.02500 0.02500 1.91735e+00 1.46111e-01 1.46111e-01 +1.47500 0.02500 0.02500 7.13561e-01 1.02489e-01 1.02489e-01 +1.52500 0.02500 0.02500 1.29645e-01 6.35916e-02 6.35916e-02 +1.57500 0.02500 0.02500 4.41370e-02 5.98525e-02 5.98525e-02 +1.62500 0.02500 0.02500 1.62515e-01 6.06476e-02 6.06476e-02 +1.67500 0.02500 0.02500 6.15551e-01 7.83663e-02 7.83663e-02 +1.72500 0.02500 0.02500 1.08982e+00 9.16829e-02 9.16829e-02 +1.77500 0.02500 0.02500 1.29044e+00 9.77287e-02 9.77287e-02 +1.82500 0.02500 0.02500 1.31079e+00 1.02564e-01 1.02564e-01 +1.87500 0.02500 0.02500 1.21806e+00 9.53119e-02 9.53119e-02 +1.92500 0.02500 0.02500 8.61740e-01 7.95639e-02 7.95639e-02 +1.97500 0.02500 0.02500 6.18126e-01 6.72227e-02 6.72227e-02 +2.02500 0.02500 0.02500 2.88671e-01 4.66405e-02 4.66405e-02 +2.07500 0.02500 0.02500 2.09002e-01 4.32806e-02 4.32806e-02 +2.12500 0.02500 0.02500 1.32147e-01 3.98234e-02 3.98234e-02 +2.17500 0.02500 0.02500 8.86300e-02 3.73953e-02 3.73953e-02 +2.22500 0.02500 0.02500 1.13896e-01 3.44593e-02 3.44593e-02 +2.27500 0.02500 0.02500 2.02709e-01 4.36286e-02 4.36286e-02 +2.32500 0.02500 0.02500 1.04124e-01 3.10663e-02 3.10663e-02 +2.37500 0.02500 0.02500 1.28319e-01 3.56242e-02 3.56242e-02 +2.42500 0.02500 0.02500 8.12150e-02 2.90157e-02 2.90157e-02 +2.47500 0.02500 0.02500 4.30790e-02 3.22583e-02 3.22583e-02 +2.55000 0.05000 0.05000 4.76500e-02 1.86345e-02 1.86345e-02 +2.65000 0.05000 0.05000 2.42110e-02 1.36672e-02 1.36672e-02 +2.75000 0.05000 0.05000 1.39450e-02 1.41180e-02 1.41180e-02 +2.85000 0.05000 0.05000 9.18100e-03 1.32598e-02 1.32598e-02 +2.95000 0.05000 0.05000 1.02280e-02 1.23726e-02 1.23726e-02 +END YODA_SCATTER2D_V2 diff --git a/analyses/pluginBABAR/BABAR_2013_I1238807.cc b/analyses/pluginBABAR/BABAR_2013_I1238807.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2013_I1238807.cc @@ -0,0 +1,84 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/Projections/FinalState.hh" + + +namespace Rivet { + + + /// @brief Add a short analysis description here + class BABAR_2013_I1238807 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(BABAR_2013_I1238807); + + + /// @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; + 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 _nkaon; + //@} + + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(BABAR_2013_I1238807); + + +} diff --git a/analyses/pluginBABAR/BABAR_2013_I1238807.info b/analyses/pluginBABAR/BABAR_2013_I1238807.info new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2013_I1238807.info @@ -0,0 +1,39 @@ +Name: BABAR_2013_I1238807 +Year: 2013 +Summary: Cross section for $e^+e^-\to K^+K^-$ between threshold and 5 GeV +Experiment: BABAR +Collider: PEP-II +InspireID: 1238807 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - Phys.Rev. D88 (2013) no.3, 032013 +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 threshold and 5 GeV using radiative return' +Keywords: +BibKey: Lees:2013gzt +BibTeX: '%%% contains utf-8, see: http://inspirehep.net/info/faq/general#utf8 +%%% add \usepackage[utf8]{inputenc} to your latex preamble + +@article{Lees:2013gzt, + author = "Lees, J. P. and others", + title = "{Precision measurement of the $e^+e^- → K^+K^-(\gamma)$ + cross section with the initial-state radiation method at + BABAR}", + collaboration = "BaBar", + journal = "Phys. Rev.", + volume = "D88", + year = "2013", + number = "3", + pages = "032013", + doi = "10.1103/PhysRevD.88.032013", + eprint = "1306.3600", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + reportNumber = "BABAR-PUB-13-006, SLAC-PUB-15487", + SLACcitation = "%%CITATION = ARXIV:1306.3600;%%" +}' diff --git a/analyses/pluginBABAR/BABAR_2013_I1238807.plot b/analyses/pluginBABAR/BABAR_2013_I1238807.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2013_I1238807.plot @@ -0,0 +1,6 @@ +BEGIN PLOT /BABAR_2013_I1238807/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 diff --git a/analyses/pluginBABAR/BABAR_2013_I1238807.yoda b/analyses/pluginBABAR/BABAR_2013_I1238807.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2013_I1238807.yoda @@ -0,0 +1,166 @@ +BEGIN YODA_SCATTER2D_V2 /REF/BABAR_2013_I1238807/d01-x01-y01IsRef: 1 +Path: /REF/BABAR_2013_I1238807/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +0.98750 0.00250 0.00250 0.03404 0.05198 0.05198 +0.99250 0.00250 0.00250 0.52023 0.20054 0.20054 +0.99750 0.00250 0.00250 3.19217 0.40650 0.40650 +1.00050 0.00050 0.00050 7.76484 1.51642 1.51642 +1.00150 0.00050 0.00050 5.94779 1.38601 1.38601 +1.00250 0.00050 0.00050 10.59020 1.80259 1.80259 +1.00350 0.00050 0.00050 9.63146 1.86215 1.86215 +1.00450 0.00050 0.00050 13.94398 2.23428 2.23428 +1.00550 0.00050 0.00050 20.33949 2.70745 2.70745 +1.00650 0.00050 0.00050 29.08580 3.12644 3.12644 +1.00750 0.00050 0.00050 30.13583 3.33408 3.33408 +1.00850 0.00050 0.00050 40.40224 3.86399 3.86399 +1.00950 0.00050 0.00050 52.32349 4.47783 4.47783 +1.01050 0.00050 0.00050 67.42239 4.78165 4.78165 +1.01150 0.00050 0.00050 104.99093 5.90879 5.90879 +1.01250 0.00050 0.00050 132.71958 7.17124 7.17124 +1.01350 0.00050 0.00050 182.32391 8.32592 8.32592 +1.01450 0.00050 0.00050 274.11728 10.53492 10.53492 +1.01550 0.00050 0.00050 444.48373 13.64209 13.64209 +1.01650 0.00050 0.00050 695.80222 17.75332 17.75332 +1.01750 0.00050 0.00050 1211.07223 21.79731 21.79731 +1.01850 0.00050 0.00050 1815.48296 28.84015 28.84015 +1.01950 0.00050 0.00050 2009.06013 29.20881 29.20881 +1.02050 0.00050 0.00050 1510.64308 24.19325 24.19325 +1.02150 0.00050 0.00050 949.83482 16.44565 16.44565 +1.02250 0.00050 0.00050 675.07793 13.31370 13.31370 +1.02350 0.00050 0.00050 455.07878 12.23909 12.23909 +1.02450 0.00050 0.00050 363.97604 10.76803 10.76803 +1.02550 0.00050 0.00050 278.69647 9.80299 9.80299 +1.02650 0.00050 0.00050 230.22939 8.37779 8.37779 +1.02750 0.00050 0.00050 171.82932 6.73530 6.73530 +1.02850 0.00050 0.00050 155.89575 7.28194 7.28194 +1.02950 0.00050 0.00050 133.72340 6.77856 6.77856 +1.03050 0.00050 0.00050 126.68812 5.96282 5.96282 +1.03150 0.00050 0.00050 101.37972 5.81777 5.81777 +1.03250 0.00050 0.00050 94.47374 5.60552 5.60552 +1.03350 0.00050 0.00050 78.97748 5.15112 5.15112 +1.03450 0.00050 0.00050 74.62534 5.07906 5.07906 +1.03550 0.00050 0.00050 65.20443 4.63823 4.63823 +1.03650 0.00050 0.00050 60.89312 4.45762 4.45762 +1.03750 0.00050 0.00050 64.81491 4.18438 4.18438 +1.03850 0.00050 0.00050 51.66602 4.08711 4.08711 +1.03950 0.00050 0.00050 45.79948 3.76907 3.76907 +1.04250 0.00250 0.00250 43.40466 1.57036 1.57036 +1.04750 0.00250 0.00250 34.99578 1.45252 1.45252 +1.05250 0.00250 0.00250 27.42988 1.27134 1.27134 +1.05750 0.00250 0.00250 22.91198 1.17089 1.17089 +1.06250 0.00250 0.00250 21.09122 1.09872 1.09872 +1.06750 0.00250 0.00250 19.21957 1.03872 1.03872 +1.07250 0.00250 0.00250 16.65472 0.95991 0.95991 +1.07750 0.00250 0.00250 16.67370 0.97234 0.97234 +1.08250 0.00250 0.00250 16.07066 0.93845 0.93845 +1.08750 0.00250 0.00250 13.80418 0.86376 0.86376 +1.09250 0.00250 0.00250 13.87001 0.88342 0.88342 +1.09750 0.00250 0.00250 12.74153 0.85121 0.85121 +1.11000 0.01000 0.01000 11.50355 0.36345 0.36345 +1.13000 0.01000 0.01000 9.76316 0.33553 0.33553 +1.15000 0.01000 0.01000 8.88819 0.32304 0.32304 +1.17000 0.01000 0.01000 7.91165 0.30661 0.30661 +1.19000 0.01000 0.01000 7.61149 0.30212 0.30212 +1.21000 0.01000 0.01000 7.27651 0.30303 0.30303 +1.23000 0.01000 0.01000 6.32698 0.28422 0.28422 +1.25000 0.01000 0.01000 7.25539 0.29315 0.29315 +1.27000 0.01000 0.01000 6.24005 0.27201 0.27201 +1.29000 0.01000 0.01000 6.02363 0.26125 0.26125 +1.31000 0.01000 0.01000 5.66601 0.25118 0.25118 +1.33000 0.01000 0.01000 5.50158 0.24554 0.24554 +1.35000 0.01000 0.01000 5.00275 0.22984 0.22984 +1.37000 0.01000 0.01000 4.75062 0.22128 0.22128 +1.39000 0.01000 0.01000 5.06118 0.22529 0.22529 +1.41000 0.01000 0.01000 4.33639 0.20842 0.20842 +1.43000 0.01000 0.01000 4.25246 0.20244 0.20244 +1.45000 0.01000 0.01000 4.10792 0.19887 0.19887 +1.47000 0.01000 0.01000 3.75452 0.18989 0.18989 +1.49000 0.01000 0.01000 3.46905 0.18099 0.18099 +1.51000 0.01000 0.01000 3.61010 0.18348 0.18348 +1.53000 0.01000 0.01000 3.22218 0.17507 0.17507 +1.55000 0.01000 0.01000 3.06929 0.17213 0.17213 +1.57000 0.01000 0.01000 3.30626 0.17399 0.17399 +1.59000 0.01000 0.01000 3.23367 0.17013 0.17013 +1.61000 0.01000 0.01000 2.91618 0.16424 0.16424 +1.63000 0.01000 0.01000 2.69413 0.15741 0.15741 +1.65000 0.01000 0.01000 2.52585 0.15017 0.15017 +1.67000 0.01000 0.01000 1.90007 0.13284 0.13284 +1.69000 0.01000 0.01000 1.18224 0.10676 0.10676 +1.71000 0.01000 0.01000 0.84654 0.08840 0.08840 +1.73000 0.01000 0.01000 0.86783 0.08836 0.08836 +1.75000 0.01000 0.01000 0.54285 0.07274 0.07274 +1.77000 0.01000 0.01000 0.33428 0.05854 0.05854 +1.79000 0.01000 0.01000 0.13643 0.04163 0.04163 +1.81000 0.01000 0.01000 0.21617 0.04525 0.04525 +1.83000 0.01000 0.01000 0.30144 0.05383 0.05383 +1.85000 0.01000 0.01000 0.20391 0.04581 0.04581 +1.87000 0.01000 0.01000 0.32506 0.05411 0.05411 +1.89000 0.01000 0.01000 0.23759 0.04843 0.04843 +1.91000 0.01000 0.01000 0.33583 0.05767 0.05767 +1.93000 0.01000 0.01000 0.33025 0.05690 0.05690 +1.95000 0.01000 0.01000 0.33111 0.05679 0.05679 +1.97000 0.01000 0.01000 0.34416 0.05801 0.05801 +1.99000 0.01000 0.01000 0.28873 0.05412 0.05412 +2.01000 0.01000 0.01000 0.37396 0.05773 0.05773 +2.03000 0.01000 0.01000 0.32173 0.05253 0.05253 +2.05000 0.01000 0.01000 0.19299 0.04382 0.04382 +2.07000 0.01000 0.01000 0.29156 0.04940 0.04940 +2.09000 0.01000 0.01000 0.15410 0.04044 0.04044 +2.11000 0.01000 0.01000 0.19067 0.04199 0.04199 +2.13000 0.01000 0.01000 0.22421 0.04437 0.04437 +2.15000 0.01000 0.01000 0.14127 0.03491 0.03491 +2.17000 0.01000 0.01000 0.09976 0.03196 0.03196 +2.19000 0.01000 0.01000 0.22013 0.04426 0.04426 +2.21000 0.01000 0.01000 0.23219 0.04380 0.04380 +2.23000 0.01000 0.01000 0.25413 0.04675 0.04675 +2.25000 0.01000 0.01000 0.31236 0.04700 0.04700 +2.27000 0.01000 0.01000 0.13028 0.03197 0.03197 +2.29000 0.01000 0.01000 0.26799 0.04314 0.04314 +2.31000 0.01000 0.01000 0.21411 0.04027 0.04027 +2.33000 0.01000 0.01000 0.15053 0.03497 0.03497 +2.35000 0.01000 0.01000 0.10188 0.02936 0.02936 +2.37000 0.01000 0.01000 0.15979 0.03388 0.03388 +2.39000 0.01000 0.01000 0.07784 0.02430 0.02430 +2.41000 0.01000 0.01000 0.11294 0.02826 0.02826 +2.43000 0.01000 0.01000 0.11989 0.03027 0.03027 +2.45000 0.01000 0.01000 0.05867 0.02199 0.02199 +2.47000 0.01000 0.01000 0.05324 0.02115 0.02115 +2.49000 0.01000 0.01000 0.10564 0.02690 0.02690 +2.51000 0.01000 0.01000 0.07906 0.02565 0.02565 +2.53000 0.01000 0.01000 0.09094 0.02659 0.02659 +2.55000 0.01000 0.01000 0.10364 0.02810 0.02810 +2.57000 0.01000 0.01000 0.05557 0.02203 0.02203 +2.59000 0.01000 0.01000 0.08155 0.02418 0.02418 +2.62000 0.02000 0.02000 0.06594 0.01464 0.01464 +2.66000 0.02000 0.02000 0.04467 0.01276 0.01276 +2.70000 0.02000 0.02000 0.04615 0.01199 0.01199 +2.74000 0.02000 0.02000 0.03700 0.01072 0.01072 +2.78000 0.02000 0.02000 0.03738 0.01095 0.01095 +2.82000 0.02000 0.02000 0.03908 0.01044 0.01044 +2.86000 0.02000 0.02000 0.02904 0.00925 0.00925 +2.90000 0.02000 0.02000 0.03853 0.01094 0.01094 +2.94000 0.02000 0.02000 0.04467 0.01061 0.01061 +2.98000 0.02000 0.02000 0.02110 0.00791 0.00791 +3.02000 0.02000 0.02000 0.02358 0.00794 0.00794 +3.06000 0.02000 0.02000 0.03378 0.00928 0.00928 +3.10000 0.02000 0.02000 0.02013 0.00692 0.00692 +3.14000 0.02000 0.02000 0.03313 0.00927 0.00927 +3.18000 0.02000 0.02000 0.01784 0.00722 0.00722 +3.22000 0.02000 0.02000 0.02287 0.00757 0.00757 +3.26000 0.02000 0.02000 0.02440 0.00793 0.00793 +3.30000 0.02000 0.02000 0.01005 0.00547 0.00547 +3.34000 0.02000 0.02000 0.03010 0.00775 0.00775 +3.38000 0.02000 0.02000 0.01215 0.00584 0.00584 +3.44000 0.04000 0.04000 0.01232 0.00357 0.00357 +3.52000 0.04000 0.04000 0.00768 0.00275 0.00275 +3.60000 0.04000 0.04000 0.00823 0.00293 0.00293 +3.68000 0.04000 0.04000 0.01103 0.00323 0.00323 +3.76000 0.04000 0.04000 0.00536 0.00246 0.00246 +3.95000 0.15000 0.15000 0.00250 0.00096 0.00096 +4.25000 0.15000 0.15000 0.00390 0.00088 0.00088 +4.55000 0.15000 0.15000 0.00158 0.00062 0.00062 +4.85000 0.15000 0.15000 0.00116 0.00053 0.00053 +END YODA_SCATTER2D_V2 diff --git a/analyses/pluginBABAR/BABAR_2017_I1511276.cc b/analyses/pluginBABAR/BABAR_2017_I1511276.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2017_I1511276.cc @@ -0,0 +1,164 @@ +// -*- 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 BABAR_2017_I1511276 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(BABAR_2017_I1511276); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + + // Initialise and register projections + declare(FinalState(), "FS"); + declare(UnstableParticles(), "UFS"); + + _nKKpi = bookCounter("TMP/KKpi"); + _nPhipi = bookCounter("TMP/Phipi"); + _nKKeta = bookCounter("TMP/KKeta"); + _nKKpipi = bookCounter("TMP/KKpipi");; + + } + + 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; + } + // stable histos + if( ntotal == 3 && nCount[130] == 1 && + nCount[310] == 1 && nCount[111] == 1) + _nKKpi->fill(event.weight()); + else if( ntotal == 4 && nCount[130] == 1 && + nCount[310] == 1 && nCount[111] == 2) + _nKKpipi->fill(event.weight()); + // unstable particles + const FinalState& ufs = apply(event, "UFS"); + foreach (const Particle& p, ufs.particles()) { + if(p.children().empty()) continue; + if(p.pdgId()!=333 && p.pdgId()!=221) continue; + map nRes=nCount; + int ncount = ntotal; + findChildren(p,nRes,ncount); + bool matched = true; + if(p.pdgId()==333 && ncount==1) { + for(auto const & val : nRes) { + if(val.first==111) { + if(val.second!=1) { + matched = false; + break; + } + } + else if(val.second!=0) { + matched = false; + break; + } + } + if(matched) + _nPhipi->fill(event.weight()); + } + else if(p.pdgId()==221 && ncount==2) { + for(auto const & val : nRes) { + if(val.first==130 || val.first==310) { + if(val.second!=1) { + matched = false; + break; + } + } + else if(val.second!=0) { + matched = false; + break; + } + } + if(matched) + _nKKeta->fill(event.weight()); + } + } + } + + + /// Normalise histograms etc., after the run + void finalize() { + for(unsigned int ix=1;ix<5;++ix) { + double sigma,error; + if(ix==1) { + sigma = _nKKpi->val(); + error = _nKKpi->err(); + } + else if (ix==2) { + sigma = _nPhipi->val(); + error = _nPhipi->err(); + } + else if (ix==3) { + sigma = _nKKeta->val(); + error = _nKKeta->err(); + } + else if (ix==4) { + sigma = _nKKpipi->val(); + error = _nKKpipi->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 _nKKpi,_nPhipi,_nKKeta,_nKKpipi; + //@} + + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(BABAR_2017_I1511276); + + +} diff --git a/analyses/pluginBABAR/BABAR_2017_I1511276.info b/analyses/pluginBABAR/BABAR_2017_I1511276.info new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2017_I1511276.info @@ -0,0 +1,37 @@ +Name: BABAR_2017_I1511276 +Year: 2017 +Summary: Cross section for $e^+e^-\to K_S^0K_L^0$ with $\pi^0$, $\eta$ and $\pi^0\pi^0$ between threshold and 4 GeV +Experiment: BABAR +Collider: PEP-II +InspireID: 1511276 +Status: VALIDATED +Authors: + - Peter Richardson +References: +- Phys.Rev. D95 (2017) no.5, 052001 +RunInfo: e+e- to hadrons +NeedCrossSection: yes +Beams: [e+,e-] +Description: + 'Cross section for $e^+e^-\to K_S^0K_L^0$ with $\pi^0$, $\eta$ and $\pi^0\pi^0$ between threshold and 4 GeV using radiative events.' +Keywords: +BibKey: TheBABAR:2017vgl +BibTeX: '@article{TheBABAR:2017vgl, + author = "Lees, J. P. and others", + title = "{Cross sections for the reactions $e^+ e^-\to K^0_S + K^0_L\pi^0$, $K^0_S K^0_L\eta$, and $K^0_S + K^0_L\pi^0\pi^0$ from events with initial-state + radiation}", + collaboration = "BaBar", + journal = "Phys. Rev.", + volume = "D95", + year = "2017", + number = "5", + pages = "052001", + doi = "10.1103/PhysRevD.95.052001", + eprint = "1701.08297", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + reportNumber = "BABAR-PUB-16-004, SLAC-PUB-16918", + SLACcitation = "%%CITATION = ARXIV:1701.08297;%%" +}' diff --git a/analyses/pluginBABAR/BABAR_2017_I1511276.plot b/analyses/pluginBABAR/BABAR_2017_I1511276.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2017_I1511276.plot @@ -0,0 +1,25 @@ +BEGIN PLOT /BABAR_2017_I1511276/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 +BEGIN PLOT /BABAR_2017_I1511276/d02-x01-y01 +Title=$\sigma(e^+e^-\to \phi\pi^0)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to \phi\pi^0)$/nb +ConnectGaps=1 +END PLOT +BEGIN PLOT /BABAR_2017_I1511276/d03-x01-y01 +Title=$\sigma(e^+e^-\to K^0_SK^0_L\eta)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to K^0_SK^0_L\eta)$/nb +ConnectGaps=1 +LogY=0 +END PLOT +BEGIN PLOT /BABAR_2017_I1511276/d04-x01-y01 +Title=$\sigma(e^+e^-\to K^0_SK^0_L\pi^0\pi^0)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to K^0_SK^0_L\pi^0\pi^0)$/nb +ConnectGaps=1 +END PLOT diff --git a/analyses/pluginBABAR/BABAR_2017_I1511276.yoda b/analyses/pluginBABAR/BABAR_2017_I1511276.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2017_I1511276.yoda @@ -0,0 +1,177 @@ +BEGIN YODA_SCATTER2D_V2 /REF/BABAR_2017_I1511276/d01-x01-y01 +IsRef: 1 +Path: /REF/BABAR_2017_I1511276/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.42500 0.00000 0.00000 2.80000e-01 7.00000e-02 7.00000e-02 +1.47500 0.00000 0.00000 9.40000e-01 1.50000e-01 1.50000e-01 +1.52500 0.00000 0.00000 1.68000e+00 2.20000e-01 2.20000e-01 +1.57500 0.00000 0.00000 2.60000e+00 2.80000e-01 2.80000e-01 +1.62500 0.00000 0.00000 2.89000e+00 2.60000e-01 2.60000e-01 +1.67500 0.00000 0.00000 3.15000e+00 3.00000e-01 3.00000e-01 +1.72500 0.00000 0.00000 2.79000e+00 2.90000e-01 2.90000e-01 +1.77500 0.00000 0.00000 1.96000e+00 2.30000e-01 2.30000e-01 +1.82500 0.00000 0.00000 1.30000e+00 1.80000e-01 1.80000e-01 +1.87500 0.00000 0.00000 1.12000e+00 1.80000e-01 1.80000e-01 +1.92500 0.00000 0.00000 7.90000e-01 1.20000e-01 1.20000e-01 +1.97500 0.00000 0.00000 5.50000e-01 1.00000e-01 1.00000e-01 +2.02500 0.00000 0.00000 5.80000e-01 1.10000e-01 1.10000e-01 +2.07500 0.00000 0.00000 5.50000e-01 1.10000e-01 1.10000e-01 +2.12500 0.00000 0.00000 4.40000e-01 9.00000e-02 9.00000e-02 +2.17500 0.00000 0.00000 4.20000e-01 9.00000e-02 9.00000e-02 +2.22500 0.00000 0.00000 3.60000e-01 8.00000e-02 8.00000e-02 +2.27500 0.00000 0.00000 4.80000e-01 9.00000e-02 9.00000e-02 +2.32500 0.00000 0.00000 3.10000e-01 7.00000e-02 7.00000e-02 +2.37500 0.00000 0.00000 1.80000e-01 7.00000e-02 7.00000e-02 +2.42500 0.00000 0.00000 2.60000e-01 7.00000e-02 7.00000e-02 +2.47500 0.00000 0.00000 2.50000e-01 8.00000e-02 8.00000e-02 +2.52500 0.00000 0.00000 1.90000e-01 6.00000e-02 6.00000e-02 +2.57500 0.00000 0.00000 9.00000e-02 4.00000e-02 4.00000e-02 +2.62500 0.00000 0.00000 1.40000e-01 5.00000e-02 5.00000e-02 +2.67500 0.00000 0.00000 7.00000e-02 3.00000e-02 3.00000e-02 +2.72500 0.00000 0.00000 1.20000e-01 5.00000e-02 5.00000e-02 +2.77500 0.00000 0.00000 1.20000e-01 8.00000e-02 8.00000e-02 +2.82500 0.00000 0.00000 1.00000e-01 6.00000e-02 6.00000e-02 +2.87500 0.00000 0.00000 1.50000e-01 1.00000e-01 1.00000e-01 +2.92500 0.00000 0.00000 4.00000e-02 3.00000e-02 3.00000e-02 +2.97500 0.00000 0.00000 1.00000e-01 1.00000e-01 1.00000e-01 +3.02500 0.00000 0.00000 6.00000e-02 3.00000e-02 3.00000e-02 +3.07500 0.00000 0.00000 3.10000e-01 9.00000e-02 9.00000e-02 +3.12500 0.00000 0.00000 2.40000e-01 1.20000e-01 1.20000e-01 +3.17500 0.00000 0.00000 5.00000e-02 4.00000e-02 4.00000e-02 +3.22500 0.00000 0.00000 2.00000e-02 3.00000e-02 3.00000e-02 +3.27500 0.00000 0.00000 4.00000e-02 3.00000e-02 3.00000e-02 +3.32500 0.00000 0.00000 5.00000e-02 4.00000e-02 4.00000e-02 +3.37500 0.00000 0.00000 3.00000e-02 4.00000e-02 4.00000e-02 +3.42500 0.00000 0.00000 2.00000e-02 1.00000e-02 1.00000e-02 +3.47500 0.00000 0.00000 5.00000e-02 3.00000e-02 3.00000e-02 +3.52500 0.00000 0.00000 1.00000e-02 3.00000e-02 3.00000e-02 +3.57500 0.00000 0.00000 4.00000e-02 3.00000e-02 3.00000e-02 +3.62500 0.00000 0.00000 3.00000e-02 3.00000e-02 3.00000e-02 +3.67500 0.00000 0.00000 0.00000e+00 1.00000e-02 1.00000e-02 +3.72500 0.00000 0.00000 4.00000e-02 4.00000e-02 4.00000e-02 +3.77500 0.00000 0.00000 1.00000e-02 1.00000e-02 1.00000e-02 +3.82500 0.00000 0.00000 3.00000e-02 2.00000e-02 2.00000e-02 +3.87500 0.00000 0.00000 0.00000e+00 1.00000e-02 1.00000e-02 +3.92500 0.00000 0.00000 2.00000e-02 2.00000e-02 2.00000e-02 +3.97500 0.00000 0.00000 1.00000e-02 1.00000e-02 1.00000e-02 +END YODA_SCATTER2D_V2 +BEGIN YODA_SCATTER2D_V2 /REF/BABAR_2017_I1511276/d02-x01-y01 +IsRef: 1 +Path: /REF/BABAR_2017_I1511276/d02-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.25000 0.00000 0.00000 0.00000e+00 2.00000e-02 2.00000e-02 +1.35000 0.00000 0.00000 0.00000e+00 2.00000e-02 2.00000e-02 +1.45000 0.00000 0.00000 6.00000e-02 6.00000e-02 6.00000e-02 +1.55000 0.00000 0.00000 1.70000e-01 6.00000e-02 6.00000e-02 +1.65000 0.00000 0.00000 4.00000e-02 3.00000e-02 3.00000e-02 +1.75000 0.00000 0.00000 5.00000e-02 3.00000e-02 3.00000e-02 +1.85000 0.00000 0.00000 2.00000e-02 2.00000e-02 2.00000e-02 +1.95000 0.00000 0.00000 0.00000e+00 1.00000e-02 1.00000e-02 +2.05000 0.00000 0.00000 3.00000e-02 2.00000e-02 2.00000e-02 +2.15000 0.00000 0.00000 0.00000e+00 1.00000e-02 1.00000e-02 +2.25000 0.00000 0.00000 1.00000e-02 1.00000e-02 1.00000e-02 +2.35000 0.00000 0.00000 2.00000e-02 2.00000e-02 2.00000e-02 +2.45000 0.00000 0.00000 2.00000e-02 1.00000e-02 1.00000e-02 +2.55000 0.00000 0.00000 0.00000e+00 1.00000e-02 1.00000e-02 +2.65000 0.00000 0.00000 1.00000e-02 1.00000e-02 1.00000e-02 +2.75000 0.00000 0.00000 0.00000e+00 1.00000e-02 1.00000e-02 +2.85000 0.00000 0.00000 0.00000e+00 1.00000e-02 1.00000e-02 +END YODA_SCATTER2D_V2 +BEGIN YODA_SCATTER2D_V2 /REF/BABAR_2017_I1511276/d03-x01-y01 +IsRef: 1 +Path: /REF/BABAR_2017_I1511276/d03-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.57500 0.00000 0.00000 1.20000e-01 1.40000e-01 1.40000e-01 +1.62500 0.00000 0.00000 5.20000e-01 1.40000e-01 1.40000e-01 +1.67500 0.00000 0.00000 7.80000e-01 1.80000e-01 1.80000e-01 +1.72500 0.00000 0.00000 1.01000e+00 2.30000e-01 2.30000e-01 +1.77500 0.00000 0.00000 5.00000e-01 1.40000e-01 1.40000e-01 +1.82500 0.00000 0.00000 5.30000e-01 1.40000e-01 1.40000e-01 +1.87500 0.00000 0.00000 4.10000e-01 1.20000e-01 1.20000e-01 +1.92500 0.00000 0.00000 3.60000e-01 1.20000e-01 1.20000e-01 +1.97500 0.00000 0.00000 2.80000e-01 9.00000e-02 9.00000e-02 +2.02500 0.00000 0.00000 4.60000e-01 1.10000e-01 1.10000e-01 +2.07500 0.00000 0.00000 4.10000e-01 1.20000e-01 1.20000e-01 +2.12500 0.00000 0.00000 4.50000e-01 1.20000e-01 1.20000e-01 +2.17500 0.00000 0.00000 4.40000e-01 1.10000e-01 1.10000e-01 +2.22500 0.00000 0.00000 4.90000e-01 1.20000e-01 1.20000e-01 +2.27500 0.00000 0.00000 3.00000e-01 1.10000e-01 1.10000e-01 +2.32500 0.00000 0.00000 3.10000e-01 9.00000e-02 9.00000e-02 +2.37500 0.00000 0.00000 3.00000e-01 1.00000e-01 1.00000e-01 +2.42500 0.00000 0.00000 3.20000e-01 1.00000e-01 1.00000e-01 +2.47500 0.00000 0.00000 1.90000e-01 7.00000e-02 7.00000e-02 +2.52500 0.00000 0.00000 1.20000e-01 6.00000e-02 6.00000e-02 +2.57500 0.00000 0.00000 2.90000e-01 9.00000e-02 9.00000e-02 +2.62500 0.00000 0.00000 1.10000e-01 6.00000e-02 6.00000e-02 +2.67500 0.00000 0.00000 1.90000e-01 8.00000e-02 8.00000e-02 +2.72500 0.00000 0.00000 1.50000e-01 6.00000e-02 6.00000e-02 +2.77500 0.00000 0.00000 1.00000e-01 5.00000e-02 5.00000e-02 +2.82500 0.00000 0.00000 8.00000e-02 5.00000e-02 5.00000e-02 +2.87500 0.00000 0.00000 1.40000e-01 6.00000e-02 6.00000e-02 +2.92500 0.00000 0.00000 5.00000e-02 4.00000e-02 4.00000e-02 +2.97500 0.00000 0.00000 5.00000e-02 6.00000e-02 6.00000e-02 +3.02500 0.00000 0.00000 1.40000e-01 6.00000e-02 6.00000e-02 +3.07500 0.00000 0.00000 2.20000e-01 6.00000e-02 6.00000e-02 +3.12500 0.00000 0.00000 9.00000e-02 5.00000e-02 5.00000e-02 +3.17500 0.00000 0.00000 3.00000e-02 4.00000e-02 4.00000e-02 +3.22500 0.00000 0.00000 -1.00000e-02 3.00000e-02 3.00000e-02 +3.27500 0.00000 0.00000 4.00000e-02 4.00000e-02 4.00000e-02 +3.32500 0.00000 0.00000 5.00000e-02 3.00000e-02 3.00000e-02 +3.37500 0.00000 0.00000 1.00000e-02 2.00000e-02 2.00000e-02 +3.42500 0.00000 0.00000 2.00000e-02 3.00000e-02 3.00000e-02 +3.47500 0.00000 0.00000 1.00000e-02 2.00000e-02 2.00000e-02 +3.52500 0.00000 0.00000 2.00000e-02 2.00000e-02 2.00000e-02 +3.57500 0.00000 0.00000 -1.00000e-02 2.00000e-02 2.00000e-02 +3.62500 0.00000 0.00000 -1.00000e-02 1.00000e-02 1.00000e-02 +3.67500 0.00000 0.00000 7.00000e-02 3.00000e-02 3.00000e-02 +3.72500 0.00000 0.00000 4.00000e-02 3.00000e-02 3.00000e-02 +3.77500 0.00000 0.00000 7.00000e-02 4.00000e-02 4.00000e-02 +3.82500 0.00000 0.00000 1.00000e-02 2.00000e-02 2.00000e-02 +3.87500 0.00000 0.00000 1.00000e-02 2.00000e-02 2.00000e-02 +3.92500 0.00000 0.00000 1.00000e-02 2.00000e-02 2.00000e-02 +3.97500 0.00000 0.00000 2.00000e-02 2.00000e-02 2.00000e-02 +END YODA_SCATTER2D_V2 +BEGIN YODA_SCATTER2D_V2 /REF/BABAR_2017_I1511276/d04-x01-y01 +IsRef: 1 +Path: /REF/BABAR_2017_I1511276/d04-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.35000 0.00000 0.00000 1.80000e-02 2.90000e-02 2.90000e-02 +1.45000 0.00000 0.00000 5.90000e-02 7.00000e-02 7.00000e-02 +1.55000 0.00000 0.00000 3.34000e-01 1.08000e-01 1.08000e-01 +1.65000 0.00000 0.00000 3.50000e-01 1.15000e-01 1.15000e-01 +1.75000 0.00000 0.00000 5.49000e-01 1.62000e-01 1.62000e-01 +1.85000 0.00000 0.00000 2.90000e-01 1.13000e-01 1.13000e-01 +1.95000 0.00000 0.00000 1.13000e-01 1.03000e-01 1.03000e-01 +2.05000 0.00000 0.00000 1.85000e-01 1.04000e-01 1.04000e-01 +2.15000 0.00000 0.00000 1.41000e-01 9.80000e-02 9.80000e-02 +2.25000 0.00000 0.00000 3.80000e-02 8.80000e-02 8.80000e-02 +2.35000 0.00000 0.00000 8.00000e-03 8.10000e-02 8.10000e-02 +2.45000 0.00000 0.00000 -2.20000e-02 7.30000e-02 7.30000e-02 +2.55000 0.00000 0.00000 8.00000e-02 7.40000e-02 7.40000e-02 +2.65000 0.00000 0.00000 -1.90000e-02 6.20000e-02 6.20000e-02 +2.75000 0.00000 0.00000 2.90000e-02 5.90000e-02 5.90000e-02 +2.85000 0.00000 0.00000 2.80000e-02 5.40000e-02 5.40000e-02 +2.95000 0.00000 0.00000 5.10000e-02 8.30000e-02 8.30000e-02 +3.05000 0.00000 0.00000 2.13000e-01 8.10000e-02 8.10000e-02 +3.15000 0.00000 0.00000 -2.00000e-03 7.10000e-02 7.10000e-02 +3.25000 0.00000 0.00000 2.00000e-03 3.30000e-02 3.30000e-02 +3.35000 0.00000 0.00000 -2.00000e-03 2.90000e-02 2.90000e-02 +3.45000 0.00000 0.00000 3.10000e-02 2.90000e-02 2.90000e-02 +3.55000 0.00000 0.00000 1.50000e-02 2.50000e-02 2.50000e-02 +3.65000 0.00000 0.00000 -1.50000e-02 2.00000e-02 2.00000e-02 +3.75000 0.00000 0.00000 5.20000e-02 2.40000e-02 2.40000e-02 +3.85000 0.00000 0.00000 2.00000e-03 1.80000e-02 1.80000e-02 +3.95000 0.00000 0.00000 -1.60000e-02 1.60000e-02 1.60000e-02 +END YODA_SCATTER2D_V2 diff --git a/analyses/pluginBABAR/BABAR_2017_I1591716.cc b/analyses/pluginBABAR/BABAR_2017_I1591716.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2017_I1591716.cc @@ -0,0 +1,143 @@ +// -*- 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 BABAR_2017_I1591716 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(BABAR_2017_I1591716); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + + // Initialise and register projections + declare(FinalState(), "FS"); + declare(UnstableParticles(), "UFS"); + + _nKKpipi = bookCounter("TMP/KKpipi"); + _nKKpieta = bookCounter("TMP/KKpieta"); + + } + + 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; + } + // stable histos + if( ntotal == 4 && nCount[310] == 1 && nCount[111] == 1 && + ( (nCount[ 321]==1 && nCount[-211]==1) || + (nCount[-321]==1 && nCount[ 211]==1))) + _nKKpipi->fill(event.weight()); + + // unstable particles + 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); + bool matched = true; + if(p.pdgId()==221 && ncount==3) { + for(auto const & val : nRes) { + if(val.first==310 || val.first==111) { + if(val.second!=1) { + matched = false; + break; + } + } + else if(abs(val.first)==321 || abs(val.first)==211) + continue; + else if(val.second!=0) { + matched = false; + break; + } + } + if(matched) { + if((nRes[321]==1 && nRes[-211]==1 && nRes[-321]==0 && nRes[211]==0) || + (nRes[321]==0 && nRes[-211]==0 && nRes[-321]==1 && nRes[211]==1)) + _nKKpieta->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 = _nKKpipi->val(); + error = _nKKpipi->err(); + } + else if (ix==2) { + sigma = _nKKpieta->val(); + error = _nKKpieta->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 _nKKpipi,_nKKpieta; + //@} + + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(BABAR_2017_I1591716); + + +} diff --git a/analyses/pluginBABAR/BABAR_2017_I1591716.info b/analyses/pluginBABAR/BABAR_2017_I1591716.info new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2017_I1591716.info @@ -0,0 +1,40 @@ +Name: BABAR_2017_I1591716 +Year: 2017 +Summary: Cross section for $e^+e^-\to K_S^0K^\pm\pi^\mp$ with $\pi^0$, $\eta$ between threshold and 4 GeV +Experiment: BABAR +Collider: PEP-II +InspireID: 1591716 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - Phys.Rev. D95 (2017) no.9, 092005 +RunInfo: e+e- to hadrons +NeedCrossSection: yes +Beams: [e+,e-] +Description: + 'Cross section for $e^+e^-\to K_S^0K^\pm\pi^\mp$ with $\pi^0$, $\eta$ between threshold and 4 GeV using radiative events.' +Keywords: +BibKey: TheBABAR:2017aph +BibTeX: '@article{TheBABAR:2017aph, + author = "Lees, J. P. and others", + title = "{Measurement of the $e^+e^-\to K^0_{\scriptscriptstyle + S}K^\pm\pi^{\mp}\pi^0$ and $K^0_{\scriptscriptstyle + S}K^\pm\pi^\mp\eta$ cross sections using initial-state + radiation}", + collaboration = "BaBar", + journal = "Phys. Rev.", + volume = "D95", + year = "2017", + number = "9", + pages = "092005", + doi = "10.1103/PhysRevD.95.092005", + eprint = "1704.05009", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + reportNumber = "BABAR-PUB-15-005, SLAC-PUB-16940", + SLACcitation = "%%CITATION = ARXIV:1704.05009;%%" +}' +ToDo: + - Implement the analysis, test it, remove this ToDo, and mark as VALIDATED :-) + diff --git a/analyses/pluginBABAR/BABAR_2017_I1591716.plot b/analyses/pluginBABAR/BABAR_2017_I1591716.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2017_I1591716.plot @@ -0,0 +1,12 @@ +BEGIN PLOT /BABAR_2017_I1591716/d01-x01-y01 +Title=$\sigma(e^+e^-\to K^0_SK^\pm\pi^\mp\pi^0)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to K^0_SK^\pm\pi^\mp\pi^0)$/nb +ConnectGaps=1 +END PLOT +BEGIN PLOT /BABAR_2017_I1591716/d02-x01-y01 +Title=$\sigma(e^+e^-\to K^0_SK^\pm\pi^\mp\eta)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to K^0_SK^\pm\pi^\mp\eta)$/nb +ConnectGaps=1 +END PLOT diff --git a/analyses/pluginBABAR/BABAR_2017_I1591716.yoda b/analyses/pluginBABAR/BABAR_2017_I1591716.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2017_I1591716.yoda @@ -0,0 +1,241 @@ +BEGIN YODA_SCATTER2D_V2 /REF/BABAR_2017_I1591716/d01-x01-y01 +IsRef: 1 +Path: /REF/BABAR_2017_I1591716/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.51000 0.00000 0.00000 5.00000e-02 3.00000e-02 3.00000e-02 +1.53000 0.00000 0.00000 5.00000e-02 3.00000e-02 3.00000e-02 +1.55000 0.00000 0.00000 2.00000e-02 2.00000e-02 2.00000e-02 +1.57000 0.00000 0.00000 6.00000e-02 4.00000e-02 4.00000e-02 +1.59000 0.00000 0.00000 1.90000e-01 6.00000e-02 6.00000e-02 +1.61000 0.00000 0.00000 1.60000e-01 6.00000e-02 6.00000e-02 +1.63000 0.00000 0.00000 3.60000e-01 9.00000e-02 9.00000e-02 +1.65000 0.00000 0.00000 5.30000e-01 1.00000e-01 1.00000e-01 +1.67000 0.00000 0.00000 5.20000e-01 1.00000e-01 1.00000e-01 +1.69000 0.00000 0.00000 7.20000e-01 1.20000e-01 1.20000e-01 +1.71000 0.00000 0.00000 7.00000e-01 1.20000e-01 1.20000e-01 +1.73000 0.00000 0.00000 1.09000e+00 1.40000e-01 1.40000e-01 +1.75000 0.00000 0.00000 9.10000e-01 1.30000e-01 1.30000e-01 +1.77000 0.00000 0.00000 1.11000e+00 1.40000e-01 1.40000e-01 +1.79000 0.00000 0.00000 1.48000e+00 1.60000e-01 1.60000e-01 +1.81000 0.00000 0.00000 1.35000e+00 1.50000e-01 1.50000e-01 +1.83000 0.00000 0.00000 1.67000e+00 1.70000e-01 1.70000e-01 +1.85000 0.00000 0.00000 1.73000e+00 1.70000e-01 1.70000e-01 +1.87000 0.00000 0.00000 1.98000e+00 1.80000e-01 1.80000e-01 +1.89000 0.00000 0.00000 2.12000e+00 1.90000e-01 1.90000e-01 +1.91000 0.00000 0.00000 1.99000e+00 1.80000e-01 1.80000e-01 +1.93000 0.00000 0.00000 2.31000e+00 1.90000e-01 1.90000e-01 +1.95000 0.00000 0.00000 2.05000e+00 1.80000e-01 1.80000e-01 +1.97000 0.00000 0.00000 2.32000e+00 1.90000e-01 1.90000e-01 +1.99000 0.00000 0.00000 2.00000e+00 1.80000e-01 1.80000e-01 +2.01000 0.00000 0.00000 1.65000e+00 1.60000e-01 1.60000e-01 +2.03000 0.00000 0.00000 1.67000e+00 1.60000e-01 1.60000e-01 +2.05000 0.00000 0.00000 1.62000e+00 1.60000e-01 1.60000e-01 +2.07000 0.00000 0.00000 1.91000e+00 1.70000e-01 1.70000e-01 +2.09000 0.00000 0.00000 1.44000e+00 1.50000e-01 1.50000e-01 +2.11000 0.00000 0.00000 1.90000e+00 1.70000e-01 1.70000e-01 +2.13000 0.00000 0.00000 1.78000e+00 1.60000e-01 1.60000e-01 +2.15000 0.00000 0.00000 1.73000e+00 1.60000e-01 1.60000e-01 +2.17000 0.00000 0.00000 1.36000e+00 1.40000e-01 1.40000e-01 +2.19000 0.00000 0.00000 1.49000e+00 1.40000e-01 1.40000e-01 +2.21000 0.00000 0.00000 1.42000e+00 1.40000e-01 1.40000e-01 +2.23000 0.00000 0.00000 1.36000e+00 1.40000e-01 1.40000e-01 +2.25000 0.00000 0.00000 1.36000e+00 1.40000e-01 1.40000e-01 +2.27000 0.00000 0.00000 1.15000e+00 1.20000e-01 1.20000e-01 +2.29000 0.00000 0.00000 9.90000e-01 1.20000e-01 1.20000e-01 +2.31000 0.00000 0.00000 9.50000e-01 1.10000e-01 1.10000e-01 +2.33000 0.00000 0.00000 1.25000e+00 1.30000e-01 1.30000e-01 +2.35000 0.00000 0.00000 9.80000e-01 1.10000e-01 1.10000e-01 +2.37000 0.00000 0.00000 9.80000e-01 1.10000e-01 1.10000e-01 +2.39000 0.00000 0.00000 6.10000e-01 9.00000e-02 9.00000e-02 +2.41000 0.00000 0.00000 1.08000e+00 1.20000e-01 1.20000e-01 +2.43000 0.00000 0.00000 8.40000e-01 1.00000e-01 1.00000e-01 +2.45000 0.00000 0.00000 1.03000e+00 1.10000e-01 1.10000e-01 +2.47000 0.00000 0.00000 9.30000e-01 1.10000e-01 1.10000e-01 +2.49000 0.00000 0.00000 7.70000e-01 1.00000e-01 1.00000e-01 +2.51000 0.00000 0.00000 6.50000e-01 9.00000e-02 9.00000e-02 +2.53000 0.00000 0.00000 7.70000e-01 1.00000e-01 1.00000e-01 +2.55000 0.00000 0.00000 8.30000e-01 1.00000e-01 1.00000e-01 +2.57000 0.00000 0.00000 7.10000e-01 9.00000e-02 9.00000e-02 +2.59000 0.00000 0.00000 8.50000e-01 1.00000e-01 1.00000e-01 +2.61000 0.00000 0.00000 5.60000e-01 8.00000e-02 8.00000e-02 +2.63000 0.00000 0.00000 4.30000e-01 7.00000e-02 7.00000e-02 +2.65000 0.00000 0.00000 5.60000e-01 8.00000e-02 8.00000e-02 +2.67000 0.00000 0.00000 6.40000e-01 9.00000e-02 9.00000e-02 +2.69000 0.00000 0.00000 4.60000e-01 7.00000e-02 7.00000e-02 +2.71000 0.00000 0.00000 6.30000e-01 8.00000e-02 8.00000e-02 +2.73000 0.00000 0.00000 4.90000e-01 7.00000e-02 7.00000e-02 +2.75000 0.00000 0.00000 5.90000e-01 8.00000e-02 8.00000e-02 +2.77000 0.00000 0.00000 3.70000e-01 6.00000e-02 6.00000e-02 +2.79000 0.00000 0.00000 5.10000e-01 7.00000e-02 7.00000e-02 +2.81000 0.00000 0.00000 3.50000e-01 6.00000e-02 6.00000e-02 +2.83000 0.00000 0.00000 3.00000e-01 6.00000e-02 6.00000e-02 +2.85000 0.00000 0.00000 3.60000e-01 6.00000e-02 6.00000e-02 +2.87000 0.00000 0.00000 4.20000e-01 7.00000e-02 7.00000e-02 +2.89000 0.00000 0.00000 2.80000e-01 5.00000e-02 5.00000e-02 +2.91000 0.00000 0.00000 4.40000e-01 7.00000e-02 7.00000e-02 +2.93000 0.00000 0.00000 3.70000e-01 6.00000e-02 6.00000e-02 +2.95000 0.00000 0.00000 2.30000e-01 5.00000e-02 5.00000e-02 +2.97000 0.00000 0.00000 2.90000e-01 6.00000e-02 6.00000e-02 +2.99000 0.00000 0.00000 4.20000e-01 7.00000e-02 7.00000e-02 +3.01000 0.00000 0.00000 4.70000e-01 7.00000e-02 7.00000e-02 +3.03000 0.00000 0.00000 2.60000e-01 5.00000e-02 5.00000e-02 +3.05000 0.00000 0.00000 3.30000e-01 6.00000e-02 6.00000e-02 +3.07000 0.00000 0.00000 3.90000e-01 6.00000e-02 6.00000e-02 +3.09000 0.00000 0.00000 2.69000e+00 1.60000e-01 1.60000e-01 +3.11000 0.00000 0.00000 1.61000e+00 1.30000e-01 1.30000e-01 +3.13000 0.00000 0.00000 3.80000e-01 6.00000e-02 6.00000e-02 +3.15000 0.00000 0.00000 3.00000e-01 5.00000e-02 5.00000e-02 +3.17000 0.00000 0.00000 2.50000e-01 5.00000e-02 5.00000e-02 +3.19000 0.00000 0.00000 1.60000e-01 4.00000e-02 4.00000e-02 +3.21000 0.00000 0.00000 2.10000e-01 4.00000e-02 4.00000e-02 +3.23000 0.00000 0.00000 1.80000e-01 4.00000e-02 4.00000e-02 +3.25000 0.00000 0.00000 1.90000e-01 4.00000e-02 4.00000e-02 +3.27000 0.00000 0.00000 2.30000e-01 5.00000e-02 5.00000e-02 +3.29000 0.00000 0.00000 1.60000e-01 4.00000e-02 4.00000e-02 +3.31000 0.00000 0.00000 1.90000e-01 4.00000e-02 4.00000e-02 +3.33000 0.00000 0.00000 7.00000e-02 3.00000e-02 3.00000e-02 +3.35000 0.00000 0.00000 1.50000e-01 4.00000e-02 4.00000e-02 +3.37000 0.00000 0.00000 1.30000e-01 3.00000e-02 3.00000e-02 +3.39000 0.00000 0.00000 1.20000e-01 3.00000e-02 3.00000e-02 +3.41000 0.00000 0.00000 1.40000e-01 3.00000e-02 3.00000e-02 +3.43000 0.00000 0.00000 1.50000e-01 4.00000e-02 4.00000e-02 +3.45000 0.00000 0.00000 1.80000e-01 4.00000e-02 4.00000e-02 +3.47000 0.00000 0.00000 9.00000e-02 3.00000e-02 3.00000e-02 +3.49000 0.00000 0.00000 1.40000e-01 4.00000e-02 4.00000e-02 +3.61000 0.00000 0.00000 1.40000e-01 3.00000e-02 3.00000e-02 +3.63000 0.00000 0.00000 7.00000e-02 2.00000e-02 2.00000e-02 +3.65000 0.00000 0.00000 1.50000e-01 4.00000e-02 4.00000e-02 +3.67000 0.00000 0.00000 1.10000e-01 3.00000e-02 3.00000e-02 +3.69000 0.00000 0.00000 1.70000e-01 4.00000e-02 4.00000e-02 +3.71000 0.00000 0.00000 1.60000e-01 4.00000e-02 4.00000e-02 +3.73000 0.00000 0.00000 7.00000e-02 2.00000e-02 2.00000e-02 +3.75000 0.00000 0.00000 8.00000e-02 2.00000e-02 2.00000e-02 +3.77000 0.00000 0.00000 8.00000e-02 3.00000e-02 3.00000e-02 +3.79000 0.00000 0.00000 5.00000e-02 2.00000e-02 2.00000e-02 +3.81000 0.00000 0.00000 9.00000e-02 3.00000e-02 3.00000e-02 +3.83000 0.00000 0.00000 7.00000e-02 2.00000e-02 2.00000e-02 +3.85000 0.00000 0.00000 4.00000e-02 2.00000e-02 2.00000e-02 +3.87000 0.00000 0.00000 4.00000e-02 2.00000e-02 2.00000e-02 +3.89000 0.00000 0.00000 1.10000e-01 3.00000e-02 3.00000e-02 +3.51000 0.00000 0.00000 5.00000e-02 2.00000e-02 2.00000e-02 +3.53000 0.00000 0.00000 1.70000e-01 4.00000e-02 4.00000e-02 +3.55000 0.00000 0.00000 9.00000e-02 3.00000e-02 3.00000e-02 +3.57000 0.00000 0.00000 8.00000e-02 3.00000e-02 3.00000e-02 +3.59000 0.00000 0.00000 1.30000e-01 3.00000e-02 3.00000e-02 +3.91000 0.00000 0.00000 8.00000e-02 2.00000e-02 2.00000e-02 +3.93000 0.00000 0.00000 8.00000e-02 3.00000e-02 3.00000e-02 +3.95000 0.00000 0.00000 5.00000e-02 2.00000e-02 2.00000e-02 +3.97000 0.00000 0.00000 1.00000e-01 3.00000e-02 3.00000e-02 +3.99000 0.00000 0.00000 8.00000e-02 2.00000e-02 2.00000e-02 +END YODA_SCATTER2D_V2 +BEGIN YODA_SCATTER2D_V2 /REF/BABAR_2017_I1591716/d02-x01-y01 +IsRef: 1 +Path: /REF/BABAR_2017_I1591716/d02-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +2.01000 0.00000 0.00000 1.00000e-02 3.00000e-02 3.00000e-02 +2.03000 0.00000 0.00000 4.00000e-02 4.00000e-02 4.00000e-02 +2.05000 0.00000 0.00000 8.00000e-02 4.00000e-02 4.00000e-02 +2.07000 0.00000 0.00000 5.00000e-02 4.00000e-02 4.00000e-02 +2.09000 0.00000 0.00000 -4.00000e-02 3.00000e-02 3.00000e-02 +2.11000 0.00000 0.00000 5.00000e-02 3.00000e-02 3.00000e-02 +2.13000 0.00000 0.00000 2.00000e-02 4.00000e-02 4.00000e-02 +2.15000 0.00000 0.00000 5.00000e-02 5.00000e-02 5.00000e-02 +2.17000 0.00000 0.00000 2.00000e-02 4.00000e-02 4.00000e-02 +2.19000 0.00000 0.00000 1.90000e-01 7.00000e-02 7.00000e-02 +2.21000 0.00000 0.00000 4.00000e-02 4.00000e-02 4.00000e-02 +2.23000 0.00000 0.00000 4.00000e-02 4.00000e-02 4.00000e-02 +2.25000 0.00000 0.00000 4.00000e-02 6.00000e-02 6.00000e-02 +2.27000 0.00000 0.00000 1.00000e-01 6.00000e-02 6.00000e-02 +2.29000 0.00000 0.00000 2.30000e-01 7.00000e-02 7.00000e-02 +2.31000 0.00000 0.00000 1.40000e-01 7.00000e-02 7.00000e-02 +2.33000 0.00000 0.00000 4.00000e-02 4.00000e-02 4.00000e-02 +2.35000 0.00000 0.00000 7.00000e-02 5.00000e-02 5.00000e-02 +2.37000 0.00000 0.00000 1.10000e-01 5.00000e-02 5.00000e-02 +2.39000 0.00000 0.00000 8.00000e-02 6.00000e-02 6.00000e-02 +2.41000 0.00000 0.00000 1.70000e-01 7.00000e-02 7.00000e-02 +2.43000 0.00000 0.00000 9.00000e-02 6.00000e-02 6.00000e-02 +2.45000 0.00000 0.00000 1.20000e-01 7.00000e-02 7.00000e-02 +2.47000 0.00000 0.00000 5.00000e-02 5.00000e-02 5.00000e-02 +2.49000 0.00000 0.00000 1.00000e-01 7.00000e-02 7.00000e-02 +2.51000 0.00000 0.00000 9.00000e-02 5.00000e-02 5.00000e-02 +2.53000 0.00000 0.00000 2.00000e-01 7.00000e-02 7.00000e-02 +2.55000 0.00000 0.00000 1.10000e-01 5.00000e-02 5.00000e-02 +2.57000 0.00000 0.00000 6.00000e-02 5.00000e-02 5.00000e-02 +2.59000 0.00000 0.00000 1.20000e-01 6.00000e-02 6.00000e-02 +2.61000 0.00000 0.00000 7.00000e-02 6.00000e-02 6.00000e-02 +2.63000 0.00000 0.00000 2.20000e-01 7.00000e-02 7.00000e-02 +2.65000 0.00000 0.00000 9.00000e-02 4.00000e-02 4.00000e-02 +2.67000 0.00000 0.00000 7.00000e-02 5.00000e-02 5.00000e-02 +2.69000 0.00000 0.00000 2.00000e-02 2.00000e-02 2.00000e-02 +2.71000 0.00000 0.00000 1.80000e-01 7.00000e-02 7.00000e-02 +2.73000 0.00000 0.00000 8.00000e-02 4.00000e-02 4.00000e-02 +2.75000 0.00000 0.00000 1.10000e-01 5.00000e-02 5.00000e-02 +2.77000 0.00000 0.00000 9.00000e-02 6.00000e-02 6.00000e-02 +2.79000 0.00000 0.00000 5.00000e-02 4.00000e-02 4.00000e-02 +2.81000 0.00000 0.00000 1.60000e-01 6.00000e-02 6.00000e-02 +2.83000 0.00000 0.00000 8.00000e-02 4.00000e-02 4.00000e-02 +2.85000 0.00000 0.00000 1.90000e-01 7.00000e-02 7.00000e-02 +2.87000 0.00000 0.00000 9.00000e-02 5.00000e-02 5.00000e-02 +2.89000 0.00000 0.00000 3.00000e-02 3.00000e-02 3.00000e-02 +2.91000 0.00000 0.00000 5.00000e-02 4.00000e-02 4.00000e-02 +2.93000 0.00000 0.00000 8.00000e-02 4.00000e-02 4.00000e-02 +2.95000 0.00000 0.00000 7.00000e-02 4.00000e-02 4.00000e-02 +2.97000 0.00000 0.00000 1.00000e-01 5.00000e-02 5.00000e-02 +2.99000 0.00000 0.00000 1.00000e-02 2.00000e-02 2.00000e-02 +3.01000 0.00000 0.00000 1.40000e-01 5.00000e-02 5.00000e-02 +3.03000 0.00000 0.00000 8.00000e-02 5.00000e-02 5.00000e-02 +3.05000 0.00000 0.00000 3.00000e-02 3.00000e-02 3.00000e-02 +3.07000 0.00000 0.00000 0.00000e+00 3.00000e-02 3.00000e-02 +3.09000 0.00000 0.00000 6.80000e-01 1.10000e-01 1.10000e-01 +3.11000 0.00000 0.00000 3.30000e-01 8.00000e-02 8.00000e-02 +3.13000 0.00000 0.00000 8.00000e-02 4.00000e-02 4.00000e-02 +3.15000 0.00000 0.00000 7.00000e-02 6.00000e-02 6.00000e-02 +3.17000 0.00000 0.00000 1.30000e-01 6.00000e-02 6.00000e-02 +3.19000 0.00000 0.00000 1.00000e-01 6.00000e-02 6.00000e-02 +3.21000 0.00000 0.00000 3.00000e-02 4.00000e-02 4.00000e-02 +3.23000 0.00000 0.00000 7.00000e-02 4.00000e-02 4.00000e-02 +3.25000 0.00000 0.00000 0.00000e+00 0.00000e+00 0.00000e+00 +3.27000 0.00000 0.00000 8.00000e-02 5.00000e-02 5.00000e-02 +3.29000 0.00000 0.00000 3.00000e-02 4.00000e-02 4.00000e-02 +3.31000 0.00000 0.00000 3.00000e-02 3.00000e-02 3.00000e-02 +3.33000 0.00000 0.00000 0.00000e+00 3.00000e-02 3.00000e-02 +3.35000 0.00000 0.00000 2.00000e-02 2.00000e-02 2.00000e-02 +3.37000 0.00000 0.00000 0.00000e+00 0.00000e+00 0.00000e+00 +3.39000 0.00000 0.00000 1.20000e-01 5.00000e-02 5.00000e-02 +3.41000 0.00000 0.00000 2.00000e-02 3.00000e-02 3.00000e-02 +3.43000 0.00000 0.00000 3.00000e-02 2.00000e-02 2.00000e-02 +3.45000 0.00000 0.00000 4.00000e-02 3.00000e-02 3.00000e-02 +3.47000 0.00000 0.00000 3.00000e-02 4.00000e-02 4.00000e-02 +3.49000 0.00000 0.00000 7.00000e-02 4.00000e-02 4.00000e-02 +3.51000 0.00000 0.00000 -2.00000e-02 2.00000e-02 2.00000e-02 +3.53000 0.00000 0.00000 7.00000e-02 4.00000e-02 4.00000e-02 +3.55000 0.00000 0.00000 7.00000e-02 4.00000e-02 4.00000e-02 +3.57000 0.00000 0.00000 2.00000e-02 2.00000e-02 2.00000e-02 +3.59000 0.00000 0.00000 -2.00000e-02 2.00000e-02 2.00000e-02 +3.61000 0.00000 0.00000 3.00000e-02 3.00000e-02 3.00000e-02 +3.63000 0.00000 0.00000 9.00000e-02 4.00000e-02 4.00000e-02 +3.65000 0.00000 0.00000 -1.00000e-02 3.00000e-02 3.00000e-02 +3.67000 0.00000 0.00000 2.00000e-02 2.00000e-02 2.00000e-02 +3.69000 0.00000 0.00000 3.00000e-02 3.00000e-02 3.00000e-02 +3.71000 0.00000 0.00000 9.00000e-02 4.00000e-02 4.00000e-02 +3.73000 0.00000 0.00000 5.00000e-02 4.00000e-02 4.00000e-02 +3.75000 0.00000 0.00000 3.00000e-02 2.00000e-02 2.00000e-02 +3.77000 0.00000 0.00000 0.00000e+00 1.00000e-02 1.00000e-02 +3.79000 0.00000 0.00000 3.00000e-02 2.00000e-02 2.00000e-02 +3.81000 0.00000 0.00000 4.00000e-02 2.00000e-02 2.00000e-02 +3.83000 0.00000 0.00000 0.00000e+00 1.00000e-02 1.00000e-02 +3.85000 0.00000 0.00000 1.00000e-02 2.00000e-02 2.00000e-02 +3.87000 0.00000 0.00000 3.00000e-02 2.00000e-02 2.00000e-02 +3.89000 0.00000 0.00000 8.00000e-02 4.00000e-02 4.00000e-02 +3.91000 0.00000 0.00000 0.00000e+00 1.00000e-02 1.00000e-02 +3.93000 0.00000 0.00000 3.00000e-02 2.00000e-02 2.00000e-02 +3.95000 0.00000 0.00000 1.00000e-02 1.00000e-02 1.00000e-02 +3.97000 0.00000 0.00000 5.00000e-02 3.00000e-02 3.00000e-02 +3.99000 0.00000 0.00000 0.00000e+00 0.00000e+00 0.00000e+00 +END YODA_SCATTER2D_V2 diff --git a/analyses/pluginBABAR/BABAR_2017_I1621593.cc b/analyses/pluginBABAR/BABAR_2017_I1621593.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2017_I1621593.cc @@ -0,0 +1,87 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/Projections/FinalState.hh" + + +namespace Rivet { + + + /// @brief Add a short analysis description here + class BABAR_2017_I1621593 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(BABAR_2017_I1621593); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + + // Initialise and register projections + declare(FinalState(), "FS"); + + _num4pi = bookCounter("TMP/num4"); + + } + + + /// 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]==1&&nCount[211]==1&&nCount[111]==2) + _num4pi->fill(event.weight()); + + } + + + /// Normalise histograms etc., after the run + void finalize() { + + double sigma = _num4pi->val(); + double error = _num4pi->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 _num4pi; + //@} + + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(BABAR_2017_I1621593); + + +} diff --git a/analyses/pluginBABAR/BABAR_2017_I1621593.info b/analyses/pluginBABAR/BABAR_2017_I1621593.info new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2017_I1621593.info @@ -0,0 +1,36 @@ +Name: BABAR_2017_I1621593 +Year: 2017 +Summary: $e^+e^-\to \pi^+\pi^-\pi^0\pi^0$ cross section from 0.85 to 4.5 GeV +Experiment: BABAR +Collider: PEP-II +InspireID: 1621593 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - Phys.Rev. D96 (2017) no.9, 092009 +RunInfo: e+ e- to hadrons +NeedCrossSection: yes +Beams: [e+, e-] +Description: + 'Measurement of the cross section for $e^+e^-\to \pi^+\pi^-\pi^0\pi^0$ between 0.85 to 4.5 GeV using radiative return.' +Keywords: [] +BibKey: TheBaBar:2017vzo +BibTeX: '@article{TheBaBar:2017vzo, + author = "Lees, J. P. and others", + title = "{Measurement of the + ${e}^{+}{e}^{{-}}{\rightarrow}{{\pi}}^{+}{{\pi}}^{{-}}{{\pi}}^{0}{{\pi}}^{0}$ + cross section using initial-state radiation at BABAR}", + collaboration = "BaBar", + journal = "Phys. Rev.", + volume = "D96", + year = "2017", + number = "9", + pages = "092009", + doi = "10.1103/PhysRevD.96.092009", + eprint = "1709.01171", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + reportNumber = "SLAC-PUB-17147, BABAR-PUB-17-002", + SLACcitation = "%%CITATION = ARXIV:1709.01171;%%" +}' diff --git a/analyses/pluginBABAR/BABAR_2017_I1621593.plot b/analyses/pluginBABAR/BABAR_2017_I1621593.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2017_I1621593.plot @@ -0,0 +1,7 @@ +BEGIN PLOT /BABAR_2017_I1621593/d01-x01-y01 +Title=$\sigma(e^+e^-\to \pi^+\pi^-\pi^0\pi^0)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to \pi^+\pi^-\pi^0\pi^0)$/nb +LogY=0 +ConnectGaps=1 +END PLOT diff --git a/analyses/pluginBABAR/BABAR_2017_I1621593.yoda b/analyses/pluginBABAR/BABAR_2017_I1621593.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2017_I1621593.yoda @@ -0,0 +1,382 @@ +BEGIN YODA_SCATTER2D_V2 /REF/BABAR_2017_I1621593/d01-x01-y01 +IsRef: 1 +Path: /REF/BABAR_2017_I1621593/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +0.85000 0.01000 0.01000 0.05000 0.11000 0.11000 +0.87000 0.01000 0.01000 0.23000 0.07000 0.07000 +0.89000 0.01000 0.01000 0.22000 0.12000 0.12000 +0.91000 0.01000 0.01000 0.30000 0.07000 0.07000 +0.93000 0.01000 0.01000 0.95000 0.16000 0.16000 +0.95000 0.01000 0.01000 2.38000 0.23000 0.23000 +0.97000 0.01000 0.01000 3.86000 0.30000 0.30000 +0.99000 0.01000 0.01000 4.75000 0.32000 0.32000 +1.01000 0.01000 0.01000 6.41000 0.39000 0.39000 +1.03000 0.01000 0.01000 7.43000 0.37000 0.37000 +1.05000 0.01000 0.01000 9.32000 0.40000 0.40000 +1.07000 0.01000 0.01000 9.59000 0.40000 0.40000 +1.09000 0.01000 0.01000 11.56000 0.43000 0.43000 +1.11000 0.01000 0.01000 12.10000 0.43000 0.43000 +1.13000 0.01000 0.01000 13.47000 0.45000 0.45000 +1.15000 0.01000 0.01000 14.67000 0.47000 0.47000 +1.17000 0.01000 0.01000 15.77000 0.47000 0.47000 +1.19000 0.01000 0.01000 16.69000 0.47000 0.47000 +1.21000 0.01000 0.01000 17.98000 0.51000 0.51000 +1.23000 0.01000 0.01000 19.89000 0.51000 0.51000 +1.25000 0.01000 0.01000 19.93000 0.50000 0.50000 +1.27000 0.01000 0.01000 21.00000 0.53000 0.53000 +1.29000 0.01000 0.01000 22.81000 0.52000 0.52000 +1.31000 0.01000 0.01000 23.68000 0.53000 0.53000 +1.33000 0.01000 0.01000 24.57000 0.53000 0.53000 +1.35000 0.01000 0.01000 25.25000 0.54000 0.54000 +1.37000 0.01000 0.01000 27.54000 0.56000 0.56000 +1.39000 0.01000 0.01000 27.55000 0.55000 0.55000 +1.41000 0.01000 0.01000 28.58000 0.55000 0.55000 +1.43000 0.01000 0.01000 30.41000 0.57000 0.57000 +1.45000 0.01000 0.01000 30.62000 0.57000 0.57000 +1.47000 0.01000 0.01000 30.75000 0.57000 0.57000 +1.49000 0.01000 0.01000 31.00000 0.56000 0.56000 +1.51000 0.01000 0.01000 30.57000 0.56000 0.56000 +1.53000 0.01000 0.01000 29.48000 0.54000 0.54000 +1.55000 0.01000 0.01000 28.21000 0.53000 0.53000 +1.57000 0.01000 0.01000 28.23000 0.53000 0.53000 +1.59000 0.01000 0.01000 26.05000 0.49000 0.49000 +1.61000 0.01000 0.01000 26.06000 0.49000 0.49000 +1.63000 0.01000 0.01000 25.26000 0.48000 0.48000 +1.65000 0.01000 0.01000 23.91000 0.46000 0.46000 +1.67000 0.01000 0.01000 23.71000 0.46000 0.46000 +1.69000 0.01000 0.01000 21.73000 0.44000 0.44000 +1.71000 0.01000 0.01000 21.07000 0.43000 0.43000 +1.73000 0.01000 0.01000 20.14000 0.43000 0.43000 +1.75000 0.01000 0.01000 18.47000 0.40000 0.40000 +1.77000 0.01000 0.01000 17.08000 0.39000 0.39000 +1.79000 0.01000 0.01000 15.33000 0.36000 0.36000 +1.81000 0.01000 0.01000 14.37000 0.35000 0.35000 +1.83000 0.01000 0.01000 12.53000 0.32000 0.32000 +1.85000 0.01000 0.01000 11.99000 0.32000 0.32000 +1.87000 0.01000 0.01000 10.53000 0.29000 0.29000 +1.89000 0.01000 0.01000 10.34000 0.29000 0.29000 +1.91000 0.01000 0.01000 9.02000 0.27000 0.27000 +1.93000 0.01000 0.01000 8.81000 0.26000 0.26000 +1.95000 0.01000 0.01000 8.68000 0.26000 0.26000 +1.97000 0.01000 0.01000 8.40000 0.26000 0.26000 +1.99000 0.01000 0.01000 8.35000 0.26000 0.26000 +2.01000 0.01000 0.01000 8.59000 0.26000 0.26000 +2.03000 0.01000 0.01000 8.84000 0.27000 0.27000 +2.05000 0.01000 0.01000 7.76000 0.24000 0.24000 +2.07000 0.01000 0.01000 7.84000 0.25000 0.25000 +2.09000 0.01000 0.01000 7.31000 0.23000 0.23000 +2.11000 0.01000 0.01000 7.00000 0.23000 0.23000 +2.13000 0.01000 0.01000 6.42000 0.22000 0.22000 +2.15000 0.01000 0.01000 6.92000 0.23000 0.23000 +2.17000 0.01000 0.01000 5.93000 0.21000 0.21000 +2.19000 0.01000 0.01000 6.26000 0.22000 0.22000 +2.21000 0.01000 0.01000 5.85000 0.21000 0.21000 +2.23000 0.01000 0.01000 5.08000 0.19000 0.19000 +2.25000 0.01000 0.01000 4.77000 0.18000 0.18000 +2.27000 0.01000 0.01000 4.72000 0.19000 0.19000 +2.29000 0.01000 0.01000 4.56000 0.18000 0.18000 +2.31000 0.01000 0.01000 4.03000 0.16000 0.16000 +2.33000 0.01000 0.01000 4.16000 0.17000 0.17000 +2.35000 0.01000 0.01000 3.82000 0.16000 0.16000 +2.37000 0.01000 0.01000 3.36000 0.15000 0.15000 +2.39000 0.01000 0.01000 3.24000 0.15000 0.15000 +2.41000 0.01000 0.01000 3.28000 0.15000 0.15000 +2.43000 0.01000 0.01000 3.45000 0.15000 0.15000 +2.45000 0.01000 0.01000 3.20000 0.15000 0.15000 +2.47000 0.01000 0.01000 2.97000 0.14000 0.14000 +2.49000 0.01000 0.01000 2.90000 0.14000 0.14000 +2.51000 0.01000 0.01000 2.87000 0.14000 0.14000 +2.53000 0.01000 0.01000 2.62000 0.13000 0.13000 +2.55000 0.01000 0.01000 2.51000 0.13000 0.13000 +2.57000 0.01000 0.01000 2.44000 0.13000 0.13000 +2.59000 0.01000 0.01000 2.27000 0.12000 0.12000 +2.61000 0.01000 0.01000 2.17000 0.12000 0.12000 +2.63000 0.01000 0.01000 1.94000 0.11000 0.11000 +2.65000 0.01000 0.01000 2.25000 0.12000 0.12000 +2.67000 0.01000 0.01000 2.13000 0.12000 0.12000 +2.69000 0.01000 0.01000 1.70000 0.10000 0.10000 +2.71000 0.01000 0.01000 1.61000 0.10000 0.10000 +2.73000 0.01000 0.01000 1.58000 0.10000 0.10000 +2.75000 0.01000 0.01000 1.92000 0.12000 0.12000 +2.77000 0.01000 0.01000 1.44000 0.10000 0.10000 +2.79000 0.01000 0.01000 1.34000 0.09000 0.09000 +2.81000 0.01000 0.01000 1.41000 0.09000 0.09000 +2.83000 0.01000 0.01000 1.44000 0.10000 0.10000 +2.85000 0.01000 0.01000 1.07000 0.08000 0.08000 +2.87000 0.01000 0.01000 1.22000 0.09000 0.09000 +2.89000 0.01000 0.01000 1.23000 0.09000 0.09000 +2.91000 0.01000 0.01000 1.11000 0.09000 0.09000 +2.93000 0.01000 0.01000 1.19000 0.09000 0.09000 +2.95000 0.01000 0.01000 1.04000 0.09000 0.09000 +2.97000 0.01000 0.01000 1.09000 0.09000 0.09000 +2.99000 0.01000 0.01000 0.97000 0.08000 0.08000 +3.01000 0.01000 0.01000 0.90000 0.08000 0.08000 +3.03000 0.01000 0.01000 0.90000 0.08000 0.08000 +3.05000 0.01000 0.01000 0.99000 0.08000 0.08000 +3.07000 0.01000 0.01000 1.21000 0.09000 0.09000 +3.09000 0.01000 0.01000 3.14000 0.19000 0.19000 +3.11000 0.01000 0.01000 1.83000 0.12000 0.12000 +3.13000 0.01000 0.01000 0.85000 0.07000 0.07000 +3.15000 0.01000 0.01000 0.80000 0.07000 0.07000 +3.17000 0.01000 0.01000 0.71000 0.06000 0.06000 +3.19000 0.01000 0.01000 0.64000 0.06000 0.06000 +3.21000 0.01000 0.01000 0.77000 0.07000 0.07000 +3.23000 0.01000 0.01000 0.61000 0.06000 0.06000 +3.25000 0.01000 0.01000 0.61000 0.06000 0.06000 +3.27000 0.01000 0.01000 0.55000 0.05000 0.05000 +3.29000 0.01000 0.01000 0.59000 0.06000 0.06000 +3.31000 0.01000 0.01000 0.51000 0.05000 0.05000 +3.33000 0.01000 0.01000 0.56000 0.06000 0.06000 +3.35000 0.01000 0.01000 0.45000 0.05000 0.05000 +3.37000 0.01000 0.01000 0.56000 0.06000 0.06000 +3.39000 0.01000 0.01000 0.51000 0.05000 0.05000 +3.41000 0.01000 0.01000 0.56000 0.06000 0.06000 +3.43000 0.01000 0.01000 0.45000 0.05000 0.05000 +3.45000 0.01000 0.01000 0.37000 0.04000 0.04000 +3.47000 0.01000 0.01000 0.39000 0.05000 0.05000 +3.49000 0.01000 0.01000 0.48000 0.05000 0.05000 +3.51000 0.01000 0.01000 0.53000 0.06000 0.06000 +3.53000 0.01000 0.01000 0.50000 0.06000 0.06000 +3.55000 0.01000 0.01000 0.44000 0.05000 0.05000 +3.57000 0.01000 0.01000 0.40000 0.05000 0.05000 +3.59000 0.01000 0.01000 0.35000 0.04000 0.04000 +3.61000 0.01000 0.01000 0.35000 0.05000 0.05000 +3.63000 0.01000 0.01000 0.37000 0.05000 0.05000 +3.65000 0.01000 0.01000 0.31000 0.04000 0.04000 +3.67000 0.01000 0.01000 0.35000 0.05000 0.05000 +3.69000 0.01000 0.01000 0.27000 0.03000 0.03000 +3.71000 0.01000 0.01000 0.28000 0.04000 0.04000 +3.73000 0.01000 0.01000 0.27000 0.04000 0.04000 +3.75000 0.01000 0.01000 0.30000 0.04000 0.04000 +3.77000 0.01000 0.01000 0.20000 0.03000 0.03000 +3.79000 0.01000 0.01000 0.26000 0.04000 0.04000 +3.81000 0.01000 0.01000 0.25000 0.04000 0.04000 +3.83000 0.01000 0.01000 0.17000 0.03000 0.03000 +3.85000 0.01000 0.01000 0.22000 0.03000 0.03000 +3.87000 0.01000 0.01000 0.24000 0.04000 0.04000 +3.89000 0.01000 0.01000 0.20000 0.03000 0.03000 +3.91000 0.01000 0.01000 0.20000 0.03000 0.03000 +3.93000 0.01000 0.01000 0.24000 0.04000 0.04000 +3.95000 0.01000 0.01000 0.14000 0.02000 0.02000 +3.97000 0.01000 0.01000 0.19000 0.03000 0.03000 +3.99000 0.01000 0.01000 0.13000 0.02000 0.02000 +4.01000 0.01000 0.01000 0.18000 0.03000 0.03000 +4.03000 0.01000 0.01000 0.17000 0.03000 0.03000 +4.05000 0.01000 0.01000 0.16000 0.03000 0.03000 +4.07000 0.01000 0.01000 0.17000 0.03000 0.03000 +4.09000 0.01000 0.01000 0.10000 0.02000 0.02000 +4.11000 0.01000 0.01000 0.14000 0.02000 0.02000 +4.13000 0.01000 0.01000 0.18000 0.03000 0.03000 +4.15000 0.01000 0.01000 0.15000 0.03000 0.03000 +4.17000 0.01000 0.01000 0.18000 0.03000 0.03000 +4.19000 0.01000 0.01000 0.15000 0.03000 0.03000 +4.21000 0.01000 0.01000 0.08000 0.02000 0.02000 +4.23000 0.01000 0.01000 0.10000 0.02000 0.02000 +4.25000 0.01000 0.01000 0.11000 0.02000 0.02000 +4.27000 0.01000 0.01000 0.15000 0.03000 0.03000 +4.29000 0.01000 0.01000 0.07000 0.02000 0.02000 +4.31000 0.01000 0.01000 0.12000 0.02000 0.02000 +4.33000 0.01000 0.01000 0.12000 0.02000 0.02000 +4.35000 0.01000 0.01000 0.13000 0.03000 0.03000 +4.37000 0.01000 0.01000 0.11000 0.02000 0.02000 +4.39000 0.01000 0.01000 0.08000 0.02000 0.02000 +4.41000 0.01000 0.01000 0.10000 0.02000 0.02000 +4.43000 0.01000 0.01000 0.11000 0.02000 0.02000 +4.45000 0.01000 0.01000 0.07000 0.02000 0.02000 +4.47000 0.01000 0.01000 0.07000 0.02000 0.02000 +4.49000 0.01000 0.01000 0.08000 0.02000 0.02000 +END YODA_SCATTER2D_V2 +BEGIN YODA_SCATTER2D_V2 /REF/BABAR_2017_I1621593/d01-x01-y02 +IsRef: 1 +Path: /REF/BABAR_2017_I1621593/d01-x01-y02 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +0.85000 0.01000 0.01000 0.05000 0.11000 0.11000 +0.87000 0.01000 0.01000 0.23000 0.07000 0.07000 +0.89000 0.01000 0.01000 0.22000 0.12000 0.12000 +0.91000 0.01000 0.01000 0.30000 0.07000 0.07000 +0.93000 0.01000 0.01000 0.95000 0.16000 0.16000 +0.95000 0.01000 0.01000 2.38000 0.23000 0.23000 +0.97000 0.01000 0.01000 3.86000 0.30000 0.30000 +0.99000 0.01000 0.01000 4.75000 0.32000 0.32000 +1.01000 0.01000 0.01000 6.41000 0.39000 0.39000 +1.03000 0.01000 0.01000 7.43000 0.37000 0.37000 +1.05000 0.01000 0.01000 9.32000 0.40000 0.40000 +1.07000 0.01000 0.01000 9.59000 0.40000 0.40000 +1.09000 0.01000 0.01000 11.56000 0.43000 0.43000 +1.11000 0.01000 0.01000 12.10000 0.43000 0.43000 +1.13000 0.01000 0.01000 13.47000 0.45000 0.45000 +1.15000 0.01000 0.01000 14.67000 0.47000 0.47000 +1.17000 0.01000 0.01000 15.77000 0.47000 0.47000 +1.19000 0.01000 0.01000 16.69000 0.47000 0.47000 +1.21000 0.01000 0.01000 17.98000 0.51000 0.51000 +1.23000 0.01000 0.01000 19.89000 0.51000 0.51000 +1.25000 0.01000 0.01000 19.93000 0.50000 0.50000 +1.27000 0.01000 0.01000 21.00000 0.53000 0.53000 +1.29000 0.01000 0.01000 22.81000 0.52000 0.52000 +1.31000 0.01000 0.01000 23.68000 0.53000 0.53000 +1.33000 0.01000 0.01000 24.57000 0.53000 0.53000 +1.35000 0.01000 0.01000 25.25000 0.54000 0.54000 +1.37000 0.01000 0.01000 27.54000 0.56000 0.56000 +1.39000 0.01000 0.01000 27.55000 0.55000 0.55000 +1.41000 0.01000 0.01000 28.58000 0.55000 0.55000 +1.43000 0.01000 0.01000 30.41000 0.57000 0.57000 +1.45000 0.01000 0.01000 30.62000 0.57000 0.57000 +1.47000 0.01000 0.01000 30.75000 0.57000 0.57000 +1.49000 0.01000 0.01000 31.00000 0.56000 0.56000 +1.51000 0.01000 0.01000 30.57000 0.56000 0.56000 +1.53000 0.01000 0.01000 29.48000 0.54000 0.54000 +1.55000 0.01000 0.01000 28.21000 0.53000 0.53000 +1.57000 0.01000 0.01000 28.23000 0.53000 0.53000 +1.59000 0.01000 0.01000 26.05000 0.49000 0.49000 +1.61000 0.01000 0.01000 26.06000 0.49000 0.49000 +1.63000 0.01000 0.01000 25.26000 0.48000 0.48000 +1.65000 0.01000 0.01000 23.91000 0.46000 0.46000 +1.67000 0.01000 0.01000 23.71000 0.46000 0.46000 +1.69000 0.01000 0.01000 21.73000 0.44000 0.44000 +1.71000 0.01000 0.01000 21.07000 0.43000 0.43000 +1.73000 0.01000 0.01000 20.14000 0.43000 0.43000 +1.75000 0.01000 0.01000 18.47000 0.40000 0.40000 +1.77000 0.01000 0.01000 17.08000 0.39000 0.39000 +1.79000 0.01000 0.01000 15.33000 0.36000 0.36000 +1.81000 0.01000 0.01000 14.37000 0.35000 0.35000 +1.83000 0.01000 0.01000 12.53000 0.32000 0.32000 +1.85000 0.01000 0.01000 11.99000 0.32000 0.32000 +1.87000 0.01000 0.01000 10.53000 0.29000 0.29000 +1.89000 0.01000 0.01000 10.34000 0.29000 0.29000 +1.91000 0.01000 0.01000 9.02000 0.27000 0.27000 +1.93000 0.01000 0.01000 8.81000 0.26000 0.26000 +1.95000 0.01000 0.01000 8.68000 0.26000 0.26000 +1.97000 0.01000 0.01000 8.40000 0.26000 0.26000 +1.99000 0.01000 0.01000 8.35000 0.26000 0.26000 +2.01000 0.01000 0.01000 8.59000 0.26000 0.26000 +2.03000 0.01000 0.01000 8.84000 0.27000 0.27000 +2.05000 0.01000 0.01000 7.76000 0.24000 0.24000 +2.07000 0.01000 0.01000 7.84000 0.25000 0.25000 +2.09000 0.01000 0.01000 7.31000 0.23000 0.23000 +2.11000 0.01000 0.01000 7.00000 0.23000 0.23000 +2.13000 0.01000 0.01000 6.42000 0.22000 0.22000 +2.15000 0.01000 0.01000 6.92000 0.23000 0.23000 +2.17000 0.01000 0.01000 5.93000 0.21000 0.21000 +2.19000 0.01000 0.01000 6.26000 0.22000 0.22000 +2.21000 0.01000 0.01000 5.85000 0.21000 0.21000 +2.23000 0.01000 0.01000 5.08000 0.19000 0.19000 +2.25000 0.01000 0.01000 4.77000 0.18000 0.18000 +2.27000 0.01000 0.01000 4.72000 0.19000 0.19000 +2.29000 0.01000 0.01000 4.56000 0.18000 0.18000 +2.31000 0.01000 0.01000 4.03000 0.16000 0.16000 +2.33000 0.01000 0.01000 4.16000 0.17000 0.17000 +2.35000 0.01000 0.01000 3.82000 0.16000 0.16000 +2.37000 0.01000 0.01000 3.36000 0.15000 0.15000 +2.39000 0.01000 0.01000 3.24000 0.15000 0.15000 +2.41000 0.01000 0.01000 3.28000 0.15000 0.15000 +2.43000 0.01000 0.01000 3.45000 0.15000 0.15000 +2.45000 0.01000 0.01000 3.20000 0.15000 0.15000 +2.47000 0.01000 0.01000 2.97000 0.14000 0.14000 +2.49000 0.01000 0.01000 2.90000 0.14000 0.14000 +2.51000 0.01000 0.01000 2.87000 0.14000 0.14000 +2.53000 0.01000 0.01000 2.62000 0.13000 0.13000 +2.55000 0.01000 0.01000 2.51000 0.13000 0.13000 +2.57000 0.01000 0.01000 2.44000 0.13000 0.13000 +2.59000 0.01000 0.01000 2.27000 0.12000 0.12000 +2.61000 0.01000 0.01000 2.17000 0.12000 0.12000 +2.63000 0.01000 0.01000 1.94000 0.11000 0.11000 +2.65000 0.01000 0.01000 2.25000 0.12000 0.12000 +2.67000 0.01000 0.01000 2.13000 0.12000 0.12000 +2.69000 0.01000 0.01000 1.70000 0.10000 0.10000 +2.71000 0.01000 0.01000 1.61000 0.10000 0.10000 +2.73000 0.01000 0.01000 1.58000 0.10000 0.10000 +2.75000 0.01000 0.01000 1.92000 0.12000 0.12000 +2.77000 0.01000 0.01000 1.44000 0.10000 0.10000 +2.79000 0.01000 0.01000 1.34000 0.09000 0.09000 +2.81000 0.01000 0.01000 1.41000 0.09000 0.09000 +2.83000 0.01000 0.01000 1.44000 0.10000 0.10000 +2.85000 0.01000 0.01000 1.07000 0.08000 0.08000 +2.87000 0.01000 0.01000 1.22000 0.09000 0.09000 +2.89000 0.01000 0.01000 1.23000 0.09000 0.09000 +2.91000 0.01000 0.01000 1.11000 0.09000 0.09000 +2.93000 0.01000 0.01000 1.19000 0.09000 0.09000 +2.95000 0.01000 0.01000 1.04000 0.09000 0.09000 +2.97000 0.01000 0.01000 1.09000 0.09000 0.09000 +2.99000 0.01000 0.01000 0.97000 0.08000 0.08000 +3.01000 0.01000 0.01000 0.90000 0.08000 0.08000 +3.03000 0.01000 0.01000 0.90000 0.08000 0.08000 +3.05000 0.01000 0.01000 0.99000 0.08000 0.08000 +3.07000 0.01000 0.01000 1.21000 0.09000 0.09000 +3.09000 0.01000 0.01000 3.14000 0.19000 0.19000 +3.11000 0.01000 0.01000 1.83000 0.12000 0.12000 +3.13000 0.01000 0.01000 0.85000 0.07000 0.07000 +3.15000 0.01000 0.01000 0.80000 0.07000 0.07000 +3.17000 0.01000 0.01000 0.71000 0.06000 0.06000 +3.19000 0.01000 0.01000 0.64000 0.06000 0.06000 +3.21000 0.01000 0.01000 0.77000 0.07000 0.07000 +3.23000 0.01000 0.01000 0.61000 0.06000 0.06000 +3.25000 0.01000 0.01000 0.61000 0.06000 0.06000 +3.27000 0.01000 0.01000 0.55000 0.05000 0.05000 +3.29000 0.01000 0.01000 0.59000 0.06000 0.06000 +3.31000 0.01000 0.01000 0.51000 0.05000 0.05000 +3.33000 0.01000 0.01000 0.56000 0.06000 0.06000 +3.35000 0.01000 0.01000 0.45000 0.05000 0.05000 +3.37000 0.01000 0.01000 0.56000 0.06000 0.06000 +3.39000 0.01000 0.01000 0.51000 0.05000 0.05000 +3.41000 0.01000 0.01000 0.56000 0.06000 0.06000 +3.43000 0.01000 0.01000 0.45000 0.05000 0.05000 +3.45000 0.01000 0.01000 0.37000 0.04000 0.04000 +3.47000 0.01000 0.01000 0.39000 0.05000 0.05000 +3.49000 0.01000 0.01000 0.48000 0.05000 0.05000 +3.51000 0.01000 0.01000 0.53000 0.06000 0.06000 +3.53000 0.01000 0.01000 0.50000 0.06000 0.06000 +3.55000 0.01000 0.01000 0.44000 0.05000 0.05000 +3.57000 0.01000 0.01000 0.40000 0.05000 0.05000 +3.59000 0.01000 0.01000 0.35000 0.04000 0.04000 +3.61000 0.01000 0.01000 0.35000 0.05000 0.05000 +3.63000 0.01000 0.01000 0.37000 0.05000 0.05000 +3.65000 0.01000 0.01000 0.31000 0.04000 0.04000 +3.67000 0.01000 0.01000 0.35000 0.05000 0.05000 +3.69000 0.01000 0.01000 0.27000 0.03000 0.03000 +3.71000 0.01000 0.01000 0.28000 0.04000 0.04000 +3.73000 0.01000 0.01000 0.27000 0.04000 0.04000 +3.75000 0.01000 0.01000 0.30000 0.04000 0.04000 +3.77000 0.01000 0.01000 0.20000 0.03000 0.03000 +3.79000 0.01000 0.01000 0.26000 0.04000 0.04000 +3.81000 0.01000 0.01000 0.25000 0.04000 0.04000 +3.83000 0.01000 0.01000 0.17000 0.03000 0.03000 +3.85000 0.01000 0.01000 0.22000 0.03000 0.03000 +3.87000 0.01000 0.01000 0.24000 0.04000 0.04000 +3.89000 0.01000 0.01000 0.20000 0.03000 0.03000 +3.91000 0.01000 0.01000 0.20000 0.03000 0.03000 +3.93000 0.01000 0.01000 0.24000 0.04000 0.04000 +3.95000 0.01000 0.01000 0.14000 0.02000 0.02000 +3.97000 0.01000 0.01000 0.19000 0.03000 0.03000 +3.99000 0.01000 0.01000 0.13000 0.02000 0.02000 +4.01000 0.01000 0.01000 0.18000 0.03000 0.03000 +4.03000 0.01000 0.01000 0.17000 0.03000 0.03000 +4.05000 0.01000 0.01000 0.16000 0.03000 0.03000 +4.07000 0.01000 0.01000 0.17000 0.03000 0.03000 +4.09000 0.01000 0.01000 0.10000 0.02000 0.02000 +4.11000 0.01000 0.01000 0.14000 0.02000 0.02000 +4.13000 0.01000 0.01000 0.18000 0.03000 0.03000 +4.15000 0.01000 0.01000 0.15000 0.03000 0.03000 +4.17000 0.01000 0.01000 0.18000 0.03000 0.03000 +4.19000 0.01000 0.01000 0.15000 0.03000 0.03000 +4.21000 0.01000 0.01000 0.08000 0.02000 0.02000 +4.23000 0.01000 0.01000 0.10000 0.02000 0.02000 +4.25000 0.01000 0.01000 0.11000 0.02000 0.02000 +4.27000 0.01000 0.01000 0.15000 0.03000 0.03000 +4.29000 0.01000 0.01000 0.07000 0.02000 0.02000 +4.31000 0.01000 0.01000 0.12000 0.02000 0.02000 +4.33000 0.01000 0.01000 0.12000 0.02000 0.02000 +4.35000 0.01000 0.01000 0.13000 0.03000 0.03000 +4.37000 0.01000 0.01000 0.11000 0.02000 0.02000 +4.39000 0.01000 0.01000 0.08000 0.02000 0.02000 +4.41000 0.01000 0.01000 0.10000 0.02000 0.02000 +4.43000 0.01000 0.01000 0.11000 0.02000 0.02000 +4.45000 0.01000 0.01000 0.07000 0.02000 0.02000 +4.47000 0.01000 0.01000 0.07000 0.02000 0.02000 +4.49000 0.01000 0.01000 0.08000 0.02000 0.02000 +END YODA_SCATTER2D_V2 diff --git a/analyses/pluginBABAR/BABAR_2018_I1647139.cc b/analyses/pluginBABAR/BABAR_2018_I1647139.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2018_I1647139.cc @@ -0,0 +1,122 @@ +// -*- 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 BABAR_2018_I1647139 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(BABAR_2018_I1647139); + + + /// @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(BABAR_2018_I1647139); + + +} diff --git a/analyses/pluginBABAR/BABAR_2018_I1647139.info b/analyses/pluginBABAR/BABAR_2018_I1647139.info new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2018_I1647139.info @@ -0,0 +1,34 @@ +Name: BABAR_2018_I1647139 +Year: 2018 +Summary: Cross section for $e^+e^-\to\eta\pi^+\pi^-$ between 1.2 and 3.5 GeV +Experiment: BABAR +Collider: PEP-II +InspireID: 1647139 +Status: VALIDATED +Authors: + - Peter Richardson +References: +- Phys.Rev. D97 (2018) 052007 +RunInfo: e+e- to hadrons +NeedCrossSection: yes +Beams: [e+,e-] +Description: + 'Measurement of the cross section for $e^+e^-\to\eta\pi^+\pi^-$ at energies between 1.2 and 3.5 GeV using radiative events.' +Keywords: +BibKey: TheBABAR:2018vvb +BibTeX: '@article{TheBABAR:2018vvb, + author = "Lees, J. P. and others", + title = "{Study of the process $e^+e^- \to \pi^+\pi^-\eta $ using + initial state radiation}", + collaboration = "BaBar", + journal = "Phys. Rev.", + volume = "D97", + year = "2018", + pages = "052007", + doi = "10.1103/PhysRevD.97.052007", + eprint = "1801.02960", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + reportNumber = "BABAR-PUB-17-003, SLAC-PUB-17204", + SLACcitation = "%%CITATION = ARXIV:1801.02960;%%" +}' diff --git a/analyses/pluginBABAR/BABAR_2018_I1647139.plot b/analyses/pluginBABAR/BABAR_2018_I1647139.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2018_I1647139.plot @@ -0,0 +1,8 @@ +BEGIN PLOT /BABAR_2018_I1647139/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 +ConnectGaps=1 +LegendXPos=0.7 +END PLOT diff --git a/analyses/pluginBABAR/BABAR_2018_I1647139.yoda b/analyses/pluginBABAR/BABAR_2018_I1647139.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2018_I1647139.yoda @@ -0,0 +1,65 @@ +BEGIN YODA_SCATTER2D_V2 /REF/BABAR_2018_I1647139/d01-x01-y01 +IsRef: 1 +Path: /REF/BABAR_2018_I1647139/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.21250 0.01250 0.01250 0.23100 0.08300 0.08300 +1.23750 0.01250 0.01250 0.05800 0.05200 0.05200 +1.26250 0.01250 0.01250 0.22800 0.08100 0.08100 +1.28750 0.01250 0.01250 0.58300 0.11200 0.11200 +1.31250 0.01250 0.01250 0.44400 0.10300 0.10300 +1.33750 0.01250 0.01250 0.91400 0.13400 0.13400 +1.36250 0.01250 0.01250 1.28000 0.15400 0.15400 +1.38750 0.01250 0.01250 1.62800 0.17000 0.17000 +1.41250 0.01250 0.01250 2.10300 0.18900 0.18900 +1.43750 0.01250 0.01250 2.92000 0.21600 0.21600 +1.46250 0.01250 0.01250 3.58200 0.23500 0.23500 +1.48750 0.01250 0.01250 3.69900 0.23700 0.23700 +1.51250 0.01250 0.01250 3.97000 0.24100 0.24100 +1.53750 0.01250 0.01250 4.03500 0.24000 0.24000 +1.56250 0.01250 0.01250 3.82600 0.23100 0.23100 +1.58750 0.01250 0.01250 3.47600 0.21800 0.21800 +1.61250 0.01250 0.01250 3.06500 0.20300 0.20300 +1.63750 0.01250 0.01250 2.73200 0.18900 0.18900 +1.66250 0.01250 0.01250 2.52800 0.17900 0.17900 +1.68750 0.01250 0.01250 2.16100 0.16600 0.16600 +1.71250 0.01250 0.01250 2.23300 0.16400 0.16400 +1.73750 0.01250 0.01250 2.13000 0.15900 0.15900 +1.76250 0.01250 0.01250 2.09100 0.15500 0.15500 +1.78750 0.01250 0.01250 1.96500 0.14900 0.14900 +1.81250 0.01250 0.01250 2.00500 0.14600 0.14600 +1.83750 0.01250 0.01250 1.26200 0.11800 0.11800 +1.86250 0.01250 0.01250 0.98700 0.10100 0.10100 +1.88750 0.01250 0.01250 0.57500 0.08100 0.08100 +1.91250 0.01250 0.01250 0.88800 0.09200 0.09200 +1.93750 0.01250 0.01250 0.72800 0.08600 0.08600 +1.96250 0.01250 0.01250 0.73600 0.08500 0.08500 +1.98750 0.01250 0.01250 0.64800 0.08100 0.08100 +2.02500 0.02500 0.02500 0.63200 0.05700 0.05700 +2.07500 0.02500 0.02500 0.54400 0.05000 0.05000 +2.12500 0.02500 0.02500 0.33700 0.03900 0.03900 +2.17500 0.02500 0.02500 0.31700 0.03800 0.03800 +2.22500 0.02500 0.02500 0.23600 0.03300 0.03300 +2.27500 0.02500 0.02500 0.20900 0.03100 0.03100 +2.32500 0.02500 0.02500 0.31700 0.03600 0.03600 +2.37500 0.02500 0.02500 0.28100 0.03300 0.03300 +2.42500 0.02500 0.02500 0.22300 0.03000 0.03000 +2.47500 0.02500 0.02500 0.28700 0.03200 0.03200 +2.52500 0.02500 0.02500 0.17300 0.02600 0.02600 +2.57500 0.02500 0.02500 0.09600 0.01900 0.01900 +2.62500 0.02500 0.02500 0.14700 0.02300 0.02300 +2.67500 0.02500 0.02500 0.09500 0.01800 0.01800 +2.72500 0.02500 0.02500 0.09700 0.01800 0.01800 +2.77500 0.02500 0.02500 0.08800 0.01700 0.01700 +2.82500 0.02500 0.02500 0.10100 0.01800 0.01800 +2.87500 0.02500 0.02500 0.07900 0.01500 0.01500 +2.92500 0.02500 0.02500 0.04400 0.01200 0.01200 +2.97500 0.02500 0.02500 0.06300 0.01400 0.01400 +3.02500 0.02500 0.02500 0.05800 0.01400 0.01400 +3.17500 0.02500 0.02500 0.03000 0.01000 0.01000 +3.25000 0.05000 0.05000 0.03300 0.00700 0.00700 +3.35000 0.05000 0.05000 0.01700 0.00500 0.00500 +3.45000 0.05000 0.05000 0.00800 0.00300 0.00300 +END YODA_SCATTER2D_V2 diff --git a/analyses/pluginBABAR/BABAR_2018_I1700745.cc b/analyses/pluginBABAR/BABAR_2018_I1700745.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2018_I1700745.cc @@ -0,0 +1,215 @@ +// -*- 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 BABAR_2018_I1700745 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(BABAR_2018_I1700745); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + + + // Initialise and register projections + declare(FinalState(), "FS"); + declare(UnstableParticles(), "UFS"); + _num5Pi = bookCounter("TMP/5Pi"); + _num2PiEta = bookCounter("TMP/2PiEta"); + _numOmegaPiPi = bookCounter("TMP/OmegaPiPi"); + _num4PiEta = bookCounter("TMP/4PiEta"); + _numOmegaPiEta = bookCounter("TMP/OmegaPiEta"); + } + + 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; + } + if(ntotal==5 && nCount[211]==1 && nCount[-211]==1 && nCount[111]==3) + _num5Pi->fill(event.weight()); + const FinalState& ufs = apply(event, "UFS"); + foreach (const Particle& p, ufs.particles()) { + if(p.children().empty()) continue; + // find eta/omegas + if(p.pdgId()==221 || p.pdgId()==223 ) { + map nRes = nCount; + int ncount = ntotal; + findChildren(p,nRes,ncount); + // eta + if(p.pdgId()==221) { + // 2 pi eta + if(ncount==2) { + 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) + _num2PiEta->fill(event.weight()); + } + // 4 pi eta + else if(ncount==4) { + bool matched = true; + for(auto const & val : nRes) { + if(abs(val.first)==211) { + if(val.second !=1) { + matched = false; + break; + } + } + else if(abs(val.first)==111) { + if(val.second !=2) { + matched = false; + break; + } + } + else if(val.second!=0) { + matched = false; + break; + } + } + if(matched) + _num4PiEta->fill(event.weight()); + } + // pi0 omega eta + foreach (const Particle& p2, ufs.particles()) { + if(p2.pdgId()!=223) continue; + map nResB = nRes; + int ncountB = ncount; + findChildren(p2,nResB,ncountB); + if(ncountB!=1) continue; + bool matched = true; + for(auto const & val : nResB) { + if(abs(val.first)==111) { + if(val.second !=1) { + matched = false; + break; + } + } + else if(val.second!=0) { + matched = false; + break; + } + } + if(matched) { + _numOmegaPiEta->fill(event.weight()); + break; + } + } + } + else { + if(ncount!=2) continue; + bool matched = true; + for(auto const & val : nRes) { + if(abs(val.first)==111) { + if(val.second !=2) { + matched = false; + break; + } + } + else if(val.second!=0) { + matched = false; + break; + } + } + if(matched) + _numOmegaPiPi->fill(event.weight()); + } + } + } + } + + /// Normalise histograms etc., after the run + void finalize() { + for(unsigned int ix=1;ix<6;++ix) { + double sigma,error; + if(ix==1) { + sigma = _num5Pi->val(); + error = _num5Pi->err(); + } + else if(ix==2) { + sigma = _num2PiEta->val(); + error = _num2PiEta->err(); + } + else if(ix==3) { + sigma = _numOmegaPiPi->val(); + error = _numOmegaPiPi->err(); + } + else if(ix==4) { + sigma = _num4PiEta->val(); + error = _num4PiEta->err(); + } + else if(ix==5) { + sigma = _numOmegaPiEta->val(); + error = _numOmegaPiEta->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 _num5Pi,_num2PiEta,_numOmegaPiPi, + _num4PiEta,_numOmegaPiEta; + //@} + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(BABAR_2018_I1700745); + + +} diff --git a/analyses/pluginBABAR/BABAR_2018_I1700745.info b/analyses/pluginBABAR/BABAR_2018_I1700745.info new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2018_I1700745.info @@ -0,0 +1,34 @@ +Name: BABAR_2018_I1700745 +Year: 2018 +Summary: Cross section for $e^+e^-\to\pi^+\pi^-3\pi^0$ and $\pi^+\pi^-2\pi^0\eta$ between threshold and 4.35 GeV +Experiment: BABAR +Collider: PEP-II +InspireID: 1700745 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - arXiv:1810.11962 +RunInfo: e+e- to hadrons +NeedCrossSection: yes +Beams: [e+,e-] +Description: + 'Measurement of the cross section for $e^+e^-\to\pi^+\pi^-3\pi^0$ and $\pi^+\pi^-2\pi^0\eta$ between threshold and 4.35 GeV. + The cross sections for the $\pi^+\pi^-\eta$, $\omega\pi^0\pi^0$ and $\omega\pi^0\eta$ resonant contributions are also measured' +Keywords: +BibKey: Lees:2018dnv +BibTeX: '@article{Lees:2018dnv, + author = "Lees, J. P.", + title = "{Study of the reactions + $e^+e^-\to\pi^+\pi^-\pi^0\pi^0\pi^0\gamma$ and + $\pi^+\pi^-\pi^0\pi^0\eta\gamma$ at center-of-mass + energies from threshold to 4.35 GeV using initial-state + radiation}", + collaboration = "BaBar", + year = "2018", + eprint = "1810.11962", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + reportNumber = "SLAC-PUB-17344", + SLACcitation = "%%CITATION = ARXIV:1810.11962;%%" +}' diff --git a/analyses/pluginBABAR/BABAR_2018_I1700745.plot b/analyses/pluginBABAR/BABAR_2018_I1700745.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2018_I1700745.plot @@ -0,0 +1,40 @@ +BEGIN PLOT /BABAR_2018_I1700745/d01-x01-y01 +Title=$\sigma(e^+e^-\to \pi^+\pi^-3\pi^0)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to \pi^+\pi^-3\pi^0)$/nb +LogY=0 +ConnectGaps=1 +LegendXPos=0.6 +END PLOT +BEGIN PLOT /BABAR_2018_I1700745/d02-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 +ConnectGaps=1 +LegendXPos=0.6 +END PLOT +BEGIN PLOT /BABAR_2018_I1700745/d03-x01-y01 +Title=$\sigma(e^+e^-\to \omega\pi^0\pi^0)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to \omega\pi^0\pi^0)$/nb +LogY=0 +ConnectGaps=1 +LegendXPos=0.6 +END PLOT +BEGIN PLOT /BABAR_2018_I1700745/d04-x01-y01 +Title=$\sigma(e^+e^-\to \eta\pi^+\pi^-2\pi^0)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to \eta\pi^+\pi^-2\pi^0)$/nb +LogY=0 +ConnectGaps=1 +LegendXPos=0.6 +END PLOT +BEGIN PLOT /BABAR_2018_I1700745/d05-x01-y01 +Title=$\sigma(e^+e^-\to \eta\omega\pi^0)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to \eta\omega\pi^0)$/nb +LogY=0 +ConnectGaps=1 +LegendXPos=0.6 +END PLOT diff --git a/analyses/pluginBABAR/BABAR_2018_I1700745.yoda b/analyses/pluginBABAR/BABAR_2018_I1700745.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginBABAR/BABAR_2018_I1700745.yoda @@ -0,0 +1,322 @@ +BEGIN YODA_SCATTER2D_V2 /REF/BABAR_2018_I1700745/d01-x01-y01 +IsRef: 1 +Path: /REF/BABAR_2018_I1700745/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.12500 0.00000 0.00000 0.00000e+00 2.00000e-02 2.00000e-02 +1.17500 0.00000 0.00000 0.00000e+00 3.00000e-02 3.00000e-02 +1.22500 0.00000 0.00000 -3.00000e-02 5.00000e-02 5.00000e-02 +1.27500 0.00000 0.00000 2.10000e-01 1.20000e-01 1.20000e-01 +1.32500 0.00000 0.00000 5.10000e-01 1.20000e-01 1.20000e-01 +1.37500 0.00000 0.00000 1.17000e+00 2.00000e-01 2.00000e-01 +1.42500 0.00000 0.00000 1.68000e+00 1.50000e-01 1.50000e-01 +1.47500 0.00000 0.00000 2.10000e+00 2.60000e-01 2.60000e-01 +1.52500 0.00000 0.00000 1.92000e+00 2.80000e-01 2.80000e-01 +1.57500 0.00000 0.00000 2.49000e+00 2.70000e-01 2.70000e-01 +1.62500 0.00000 0.00000 2.36000e+00 2.70000e-01 2.70000e-01 +1.67500 0.00000 0.00000 2.81000e+00 2.00000e-01 2.00000e-01 +1.72500 0.00000 0.00000 2.20000e+00 2.50000e-01 2.50000e-01 +1.77500 0.00000 0.00000 2.20000e+00 2.30000e-01 2.30000e-01 +1.82500 0.00000 0.00000 2.03000e+00 1.70000e-01 1.70000e-01 +1.87500 0.00000 0.00000 1.65000e+00 1.50000e-01 1.50000e-01 +1.92500 0.00000 0.00000 1.23000e+00 1.50000e-01 1.50000e-01 +1.97500 0.00000 0.00000 1.46000e+00 1.90000e-01 1.90000e-01 +2.02500 0.00000 0.00000 1.41000e+00 1.40000e-01 1.40000e-01 +2.07500 0.00000 0.00000 1.42000e+00 1.40000e-01 1.40000e-01 +2.12500 0.00000 0.00000 1.30000e+00 1.20000e-01 1.20000e-01 +2.17500 0.00000 0.00000 1.12000e+00 1.30000e-01 1.30000e-01 +2.22500 0.00000 0.00000 1.16000e+00 1.30000e-01 1.30000e-01 +2.27500 0.00000 0.00000 1.03000e+00 1.20000e-01 1.20000e-01 +2.32500 0.00000 0.00000 8.20000e-01 1.10000e-01 1.10000e-01 +2.37500 0.00000 0.00000 6.80000e-01 1.00000e-01 1.00000e-01 +2.42500 0.00000 0.00000 9.20000e-01 1.00000e-01 1.00000e-01 +2.47500 0.00000 0.00000 6.10000e-01 9.00000e-02 9.00000e-02 +2.52500 0.00000 0.00000 4.50000e-01 8.00000e-02 8.00000e-02 +2.57500 0.00000 0.00000 7.10000e-01 1.00000e-01 1.00000e-01 +2.62500 0.00000 0.00000 4.50000e-01 8.00000e-02 8.00000e-02 +2.67500 0.00000 0.00000 5.60000e-01 9.00000e-02 9.00000e-02 +2.72500 0.00000 0.00000 2.20000e-01 8.00000e-02 8.00000e-02 +2.77500 0.00000 0.00000 4.00000e-01 8.00000e-02 8.00000e-02 +2.82500 0.00000 0.00000 2.90000e-01 8.00000e-02 8.00000e-02 +2.87500 0.00000 0.00000 6.20000e-01 8.00000e-02 8.00000e-02 +2.92500 0.00000 0.00000 5.50000e-01 8.00000e-02 8.00000e-02 +2.97500 0.00000 0.00000 6.00000e-01 9.00000e-02 9.00000e-02 +3.02500 0.00000 0.00000 8.50000e-01 1.00000e-01 1.00000e-01 +3.07500 0.00000 0.00000 4.36000e+00 1.30000e-01 1.30000e-01 +3.12500 0.00000 0.00000 2.66000e+00 1.10000e-01 1.10000e-01 +3.17500 0.00000 0.00000 6.00000e-01 6.00000e-02 6.00000e-02 +3.22500 0.00000 0.00000 3.30000e-01 5.00000e-02 5.00000e-02 +3.27500 0.00000 0.00000 3.10000e-01 5.00000e-02 5.00000e-02 +3.32500 0.00000 0.00000 2.00000e-01 5.00000e-02 5.00000e-02 +3.37500 0.00000 0.00000 3.50000e-01 5.00000e-02 5.00000e-02 +3.42500 0.00000 0.00000 2.20000e-01 5.00000e-02 5.00000e-02 +3.47500 0.00000 0.00000 1.90000e-01 5.00000e-02 5.00000e-02 +3.52500 0.00000 0.00000 2.60000e-01 5.00000e-02 5.00000e-02 +3.57500 0.00000 0.00000 1.20000e-01 5.00000e-02 5.00000e-02 +3.62500 0.00000 0.00000 3.80000e-01 5.00000e-02 5.00000e-02 +3.67500 0.00000 0.00000 4.10000e-01 6.00000e-02 6.00000e-02 +3.72500 0.00000 0.00000 2.90000e-01 5.00000e-02 5.00000e-02 +3.77500 0.00000 0.00000 1.50000e-01 4.00000e-02 4.00000e-02 +3.82500 0.00000 0.00000 2.00000e-01 4.00000e-02 4.00000e-02 +3.87500 0.00000 0.00000 1.80000e-01 4.00000e-02 4.00000e-02 +3.92500 0.00000 0.00000 1.40000e-01 4.00000e-02 4.00000e-02 +3.97500 0.00000 0.00000 2.20000e-01 4.00000e-02 4.00000e-02 +4.02500 0.00000 0.00000 1.40000e-01 4.00000e-02 4.00000e-02 +4.07500 0.00000 0.00000 1.40000e-01 3.00000e-02 3.00000e-02 +4.12500 0.00000 0.00000 4.00000e-02 3.00000e-02 3.00000e-02 +4.17500 0.00000 0.00000 8.00000e-02 3.00000e-02 3.00000e-02 +4.22500 0.00000 0.00000 9.00000e-02 3.00000e-02 3.00000e-02 +4.27500 0.00000 0.00000 1.20000e-01 3.00000e-02 3.00000e-02 +4.32500 0.00000 0.00000 9.00000e-02 3.00000e-02 3.00000e-02 +END YODA_SCATTER2D_V2 +BEGIN YODA_SCATTER2D_V2 /REF/BABAR_2018_I1700745/d02-x01-y01 +IsRef: 1 +Path: /REF/BABAR_2018_I1700745/d02-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.07500 0.00000 0.00000 6.00000e-02 3.00000e-02 3.00000e-02 +1.12500 0.00000 0.00000 2.90000e-01 2.30000e-01 2.30000e-01 +1.17500 0.00000 0.00000 0.00000e+00 1.20000e-01 1.20000e-01 +1.22500 0.00000 0.00000 2.30000e-01 2.50000e-01 2.50000e-01 +1.27500 0.00000 0.00000 5.70000e-01 2.70000e-01 2.70000e-01 +1.32500 0.00000 0.00000 1.15000e+00 3.40000e-01 3.40000e-01 +1.37500 0.00000 0.00000 1.83000e+00 3.60000e-01 3.60000e-01 +1.42500 0.00000 0.00000 2.74000e+00 4.00000e-01 4.00000e-01 +1.47500 0.00000 0.00000 3.74000e+00 4.30000e-01 4.30000e-01 +1.52500 0.00000 0.00000 4.14000e+00 4.40000e-01 4.40000e-01 +1.57500 0.00000 0.00000 3.48000e+00 4.00000e-01 4.00000e-01 +1.62500 0.00000 0.00000 2.67000e+00 3.60000e-01 3.60000e-01 +1.67500 0.00000 0.00000 2.52000e+00 3.20000e-01 3.20000e-01 +1.72500 0.00000 0.00000 2.20000e+00 3.00000e-01 3.00000e-01 +1.77500 0.00000 0.00000 2.38000e+00 2.90000e-01 2.90000e-01 +1.82500 0.00000 0.00000 1.39000e+00 2.30000e-01 2.30000e-01 +1.87500 0.00000 0.00000 1.16000e+00 2.10000e-01 2.10000e-01 +1.92500 0.00000 0.00000 1.00000e+00 1.90000e-01 1.90000e-01 +1.97500 0.00000 0.00000 6.50000e-01 1.60000e-01 1.60000e-01 +2.02500 0.00000 0.00000 7.80000e-01 1.60000e-01 1.60000e-01 +2.07500 0.00000 0.00000 5.10000e-01 1.30000e-01 1.30000e-01 +2.12500 0.00000 0.00000 5.00000e-01 1.30000e-01 1.30000e-01 +2.17500 0.00000 0.00000 7.50000e-01 1.30000e-01 1.30000e-01 +2.22500 0.00000 0.00000 2.30000e-01 1.10000e-01 1.10000e-01 +2.27500 0.00000 0.00000 3.50000e-01 1.00000e-01 1.00000e-01 +2.32500 0.00000 0.00000 2.20000e-01 9.00000e-02 9.00000e-02 +2.37500 0.00000 0.00000 3.30000e-01 9.00000e-02 9.00000e-02 +2.42500 0.00000 0.00000 2.20000e-01 7.00000e-02 7.00000e-02 +2.47500 0.00000 0.00000 5.10000e-01 1.00000e-01 1.00000e-01 +2.52500 0.00000 0.00000 2.70000e-01 9.00000e-02 9.00000e-02 +2.57500 0.00000 0.00000 8.00000e-02 5.00000e-02 5.00000e-02 +2.62500 0.00000 0.00000 1.20000e-01 6.00000e-02 6.00000e-02 +2.67500 0.00000 0.00000 2.70000e-01 7.00000e-02 7.00000e-02 +2.72500 0.00000 0.00000 1.10000e-01 5.00000e-02 5.00000e-02 +2.77500 0.00000 0.00000 9.00000e-02 5.00000e-02 5.00000e-02 +2.82500 0.00000 0.00000 3.00000e-02 4.00000e-02 4.00000e-02 +2.87500 0.00000 0.00000 5.00000e-02 4.00000e-02 4.00000e-02 +2.92500 0.00000 0.00000 2.00000e-02 4.00000e-02 4.00000e-02 +2.97500 0.00000 0.00000 9.00000e-02 5.00000e-02 5.00000e-02 +3.02500 0.00000 0.00000 5.00000e-02 5.00000e-02 5.00000e-02 +END YODA_SCATTER2D_V2 +BEGIN YODA_SCATTER2D_V2 /REF/BABAR_2018_I1700745/d03-x01-y01 +IsRef: 1 +Path: /REF/BABAR_2018_I1700745/d03-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.12500 0.00000 0.00000 4.00000e-02 8.00000e-02 8.00000e-02 +1.17500 0.00000 0.00000 3.00000e-02 1.00000e-01 1.00000e-01 +1.22500 0.00000 0.00000 -2.00000e-02 1.00000e-01 1.00000e-01 +1.27500 0.00000 0.00000 1.30000e-01 1.10000e-01 1.10000e-01 +1.32500 0.00000 0.00000 4.10000e-01 1.30000e-01 1.30000e-01 +1.37500 0.00000 0.00000 6.90000e-01 1.80000e-01 1.80000e-01 +1.42500 0.00000 0.00000 2.90000e-01 1.80000e-01 1.80000e-01 +1.47500 0.00000 0.00000 6.80000e-01 1.90000e-01 1.90000e-01 +1.52500 0.00000 0.00000 1.05000e+00 2.10000e-01 2.10000e-01 +1.57500 0.00000 0.00000 1.44000e+00 2.20000e-01 2.20000e-01 +1.62500 0.00000 0.00000 1.40000e+00 2.10000e-01 2.10000e-01 +1.67500 0.00000 0.00000 1.55000e+00 2.00000e-01 2.00000e-01 +1.72500 0.00000 0.00000 9.60000e-01 1.80000e-01 1.80000e-01 +1.77500 0.00000 0.00000 8.80000e-01 1.60000e-01 1.60000e-01 +1.82500 0.00000 0.00000 6.20000e-01 1.40000e-01 1.40000e-01 +1.87500 0.00000 0.00000 9.60000e-01 1.40000e-01 1.40000e-01 +1.92500 0.00000 0.00000 6.10000e-01 1.30000e-01 1.30000e-01 +1.97500 0.00000 0.00000 4.50000e-01 1.10000e-01 1.10000e-01 +2.02500 0.00000 0.00000 4.70000e-01 1.00000e-01 1.00000e-01 +2.07500 0.00000 0.00000 3.30000e-01 9.00000e-02 9.00000e-02 +2.12500 0.00000 0.00000 2.90000e-01 9.00000e-02 9.00000e-02 +2.17500 0.00000 0.00000 2.60000e-01 8.00000e-02 8.00000e-02 +2.22500 0.00000 0.00000 4.00000e-01 8.00000e-02 8.00000e-02 +2.27500 0.00000 0.00000 3.10000e-01 7.00000e-02 7.00000e-02 +2.32500 0.00000 0.00000 2.10000e-01 6.00000e-02 6.00000e-02 +2.37500 0.00000 0.00000 2.30000e-01 6.00000e-02 6.00000e-02 +2.42500 0.00000 0.00000 7.00000e-02 5.00000e-02 5.00000e-02 +2.47500 0.00000 0.00000 1.20000e-01 5.00000e-02 5.00000e-02 +2.52500 0.00000 0.00000 2.10000e-01 5.00000e-02 5.00000e-02 +2.57500 0.00000 0.00000 1.50000e-01 4.00000e-02 4.00000e-02 +2.62500 0.00000 0.00000 1.30000e-01 4.00000e-02 4.00000e-02 +2.67500 0.00000 0.00000 1.20000e-01 4.00000e-02 4.00000e-02 +2.72500 0.00000 0.00000 1.70000e-01 4.00000e-02 4.00000e-02 +2.77500 0.00000 0.00000 1.00000e-01 4.00000e-02 4.00000e-02 +2.82500 0.00000 0.00000 1.10000e-01 4.00000e-02 4.00000e-02 +2.87500 0.00000 0.00000 1.80000e-01 4.00000e-02 4.00000e-02 +2.92500 0.00000 0.00000 1.00000e-01 3.00000e-02 3.00000e-02 +2.97500 0.00000 0.00000 1.40000e-01 6.00000e-02 6.00000e-02 +3.02500 0.00000 0.00000 2.50000e-01 4.00000e-02 4.00000e-02 +3.07500 0.00000 0.00000 8.30000e-01 7.00000e-02 7.00000e-02 +3.12500 0.00000 0.00000 5.20000e-01 5.00000e-02 5.00000e-02 +3.17500 0.00000 0.00000 1.10000e-01 3.00000e-02 3.00000e-02 +3.22500 0.00000 0.00000 8.00000e-02 2.00000e-02 2.00000e-02 +3.27500 0.00000 0.00000 8.00000e-02 2.00000e-02 2.00000e-02 +3.32500 0.00000 0.00000 7.00000e-02 2.00000e-02 2.00000e-02 +3.37500 0.00000 0.00000 6.00000e-02 2.00000e-02 2.00000e-02 +3.42500 0.00000 0.00000 7.00000e-02 2.00000e-02 2.00000e-02 +3.47500 0.00000 0.00000 3.00000e-02 2.00000e-02 2.00000e-02 +3.52500 0.00000 0.00000 7.00000e-02 2.00000e-02 2.00000e-02 +3.57500 0.00000 0.00000 4.00000e-02 2.00000e-02 2.00000e-02 +3.62500 0.00000 0.00000 6.00000e-02 2.00000e-02 2.00000e-02 +3.67500 0.00000 0.00000 1.10000e-01 3.00000e-02 3.00000e-02 +3.72500 0.00000 0.00000 6.00000e-02 2.00000e-02 2.00000e-02 +3.77500 0.00000 0.00000 3.00000e-02 2.00000e-02 2.00000e-02 +3.82500 0.00000 0.00000 3.00000e-02 1.00000e-02 1.00000e-02 +3.87500 0.00000 0.00000 2.00000e-02 1.00000e-02 1.00000e-02 +3.92500 0.00000 0.00000 3.00000e-02 2.00000e-02 2.00000e-02 +3.97500 0.00000 0.00000 4.00000e-02 1.00000e-02 1.00000e-02 +4.02500 0.00000 0.00000 3.00000e-02 1.00000e-02 1.00000e-02 +4.07500 0.00000 0.00000 2.00000e-02 1.00000e-02 1.00000e-02 +4.12500 0.00000 0.00000 3.00000e-02 1.00000e-02 1.00000e-02 +4.17500 0.00000 0.00000 2.00000e-02 1.00000e-02 1.00000e-02 +4.22500 0.00000 0.00000 1.00000e-02 1.00000e-02 1.00000e-02 +4.27500 0.00000 0.00000 1.00000e-02 1.00000e-02 1.00000e-02 +4.32500 0.00000 0.00000 2.00000e-02 1.00000e-02 1.00000e-02 +END YODA_SCATTER2D_V2 +BEGIN YODA_SCATTER2D_V2 /REF/BABAR_2018_I1700745/d04-x01-y01 +IsRef: 1 +Path: /REF/BABAR_2018_I1700745/d04-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.62500 0.00000 0.00000 1.00000e-02 1.00000e-01 1.00000e-01 +1.67500 0.00000 0.00000 -5.00000e-02 8.00000e-02 8.00000e-02 +1.72500 0.00000 0.00000 2.00000e-01 1.00000e-01 1.00000e-01 +1.77500 0.00000 0.00000 5.10000e-01 1.20000e-01 1.20000e-01 +1.82500 0.00000 0.00000 7.10000e-01 1.40000e-01 1.40000e-01 +1.87500 0.00000 0.00000 7.30000e-01 1.40000e-01 1.40000e-01 +1.92500 0.00000 0.00000 1.22000e+00 1.60000e-01 1.60000e-01 +1.97500 0.00000 0.00000 2.22000e+00 2.00000e-01 2.00000e-01 +2.02500 0.00000 0.00000 2.01000e+00 1.90000e-01 1.90000e-01 +2.07500 0.00000 0.00000 1.61000e+00 1.80000e-01 1.80000e-01 +2.12500 0.00000 0.00000 1.90000e+00 1.80000e-01 1.80000e-01 +2.17500 0.00000 0.00000 1.59000e+00 1.60000e-01 1.60000e-01 +2.22500 0.00000 0.00000 1.66000e+00 1.80000e-01 1.80000e-01 +2.27500 0.00000 0.00000 1.29000e+00 1.60000e-01 1.60000e-01 +2.32500 0.00000 0.00000 1.27000e+00 1.50000e-01 1.50000e-01 +2.37500 0.00000 0.00000 1.70000e+00 1.80000e-01 1.80000e-01 +2.42500 0.00000 0.00000 1.30000e+00 1.50000e-01 1.50000e-01 +2.47500 0.00000 0.00000 1.27000e+00 1.60000e-01 1.60000e-01 +2.52500 0.00000 0.00000 1.00000e+00 1.30000e-01 1.30000e-01 +2.57500 0.00000 0.00000 9.50000e-01 1.50000e-01 1.50000e-01 +2.62500 0.00000 0.00000 1.11000e+00 1.60000e-01 1.60000e-01 +2.67500 0.00000 0.00000 6.70000e-01 1.40000e-01 1.40000e-01 +2.72500 0.00000 0.00000 1.07000e+00 1.30000e-01 1.30000e-01 +2.77500 0.00000 0.00000 9.70000e-01 1.40000e-01 1.40000e-01 +2.82500 0.00000 0.00000 6.80000e-01 1.40000e-01 1.40000e-01 +2.87500 0.00000 0.00000 1.00000e+00 1.30000e-01 1.30000e-01 +2.92500 0.00000 0.00000 8.10000e-01 1.30000e-01 1.30000e-01 +2.97500 0.00000 0.00000 9.60000e-01 1.30000e-01 1.30000e-01 +3.02500 0.00000 0.00000 6.10000e-01 1.40000e-01 1.40000e-01 +3.07500 0.00000 0.00000 1.21000e+00 1.60000e-01 1.60000e-01 +3.12500 0.00000 0.00000 1.06000e+00 1.50000e-01 1.50000e-01 +3.17500 0.00000 0.00000 5.00000e-01 1.20000e-01 1.20000e-01 +3.22500 0.00000 0.00000 5.20000e-01 1.10000e-01 1.10000e-01 +3.27500 0.00000 0.00000 2.60000e-01 9.00000e-02 9.00000e-02 +3.32500 0.00000 0.00000 1.50000e-01 1.10000e-01 1.10000e-01 +3.37500 0.00000 0.00000 5.00000e-01 1.00000e-01 1.00000e-01 +3.42500 0.00000 0.00000 1.50000e-01 1.10000e-01 1.10000e-01 +3.47500 0.00000 0.00000 3.40000e-01 1.00000e-01 1.00000e-01 +3.52500 0.00000 0.00000 3.00000e-01 8.00000e-02 8.00000e-02 +3.57500 0.00000 0.00000 1.80000e-01 9.00000e-02 9.00000e-02 +3.62500 0.00000 0.00000 2.00000e-01 1.10000e-01 1.10000e-01 +3.67500 0.00000 0.00000 1.80000e-01 9.00000e-02 9.00000e-02 +3.72500 0.00000 0.00000 2.80000e-01 9.00000e-02 9.00000e-02 +3.77500 0.00000 0.00000 6.00000e-02 9.00000e-02 9.00000e-02 +3.82500 0.00000 0.00000 2.00000e-02 7.00000e-02 7.00000e-02 +3.87500 0.00000 0.00000 8.00000e-02 8.00000e-02 8.00000e-02 +3.92500 0.00000 0.00000 1.20000e-01 7.00000e-02 7.00000e-02 +3.97500 0.00000 0.00000 -2.00000e-02 8.00000e-02 8.00000e-02 +4.02500 0.00000 0.00000 -4.00000e-02 8.00000e-02 8.00000e-02 +4.07500 0.00000 0.00000 1.00000e-01 6.00000e-02 6.00000e-02 +4.12500 0.00000 0.00000 1.40000e-01 7.00000e-02 7.00000e-02 +4.17500 0.00000 0.00000 -6.00000e-02 7.00000e-02 7.00000e-02 +4.22500 0.00000 0.00000 5.00000e-02 6.00000e-02 6.00000e-02 +4.27500 0.00000 0.00000 1.00000e-01 6.00000e-02 6.00000e-02 +4.32500 0.00000 0.00000 4.00000e-02 6.00000e-02 6.00000e-02 +END YODA_SCATTER2D_V2 +BEGIN YODA_SCATTER2D_V2 /REF/BABAR_2018_I1700745/d05-x01-y01 +IsRef: 1 +Path: /REF/BABAR_2018_I1700745/d05-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.52500 0.00000 0.00000 2.00000e-02 1.00000e-01 1.00000e-01 +1.57500 0.00000 0.00000 3.00000e-02 7.00000e-02 7.00000e-02 +1.62500 0.00000 0.00000 2.40000e-01 1.00000e-01 1.00000e-01 +1.67500 0.00000 0.00000 2.00000e-01 1.00000e-01 1.00000e-01 +1.72500 0.00000 0.00000 3.00000e-01 1.10000e-01 1.10000e-01 +1.77500 0.00000 0.00000 7.60000e-01 1.50000e-01 1.50000e-01 +1.82500 0.00000 0.00000 9.60000e-01 1.60000e-01 1.60000e-01 +1.87500 0.00000 0.00000 8.80000e-01 1.60000e-01 1.60000e-01 +1.92500 0.00000 0.00000 1.46000e+00 1.80000e-01 1.80000e-01 +1.97500 0.00000 0.00000 1.62000e+00 2.00000e-01 2.00000e-01 +2.02500 0.00000 0.00000 1.54000e+00 1.90000e-01 1.90000e-01 +2.07500 0.00000 0.00000 1.16000e+00 1.60000e-01 1.60000e-01 +2.12500 0.00000 0.00000 1.26000e+00 1.70000e-01 1.70000e-01 +2.17500 0.00000 0.00000 1.06000e+00 1.40000e-01 1.40000e-01 +2.22500 0.00000 0.00000 8.30000e-01 1.30000e-01 1.30000e-01 +2.27500 0.00000 0.00000 7.40000e-01 1.20000e-01 1.20000e-01 +2.32500 0.00000 0.00000 4.70000e-01 1.00000e-01 1.00000e-01 +2.37500 0.00000 0.00000 6.80000e-01 1.10000e-01 1.10000e-01 +2.42500 0.00000 0.00000 5.80000e-01 1.00000e-01 1.00000e-01 +2.47500 0.00000 0.00000 4.10000e-01 9.00000e-02 9.00000e-02 +2.52500 0.00000 0.00000 4.50000e-01 9.00000e-02 9.00000e-02 +2.57500 0.00000 0.00000 4.80000e-01 9.00000e-02 9.00000e-02 +2.62500 0.00000 0.00000 4.10000e-01 8.00000e-02 8.00000e-02 +2.67500 0.00000 0.00000 3.90000e-01 8.00000e-02 8.00000e-02 +2.72500 0.00000 0.00000 3.50000e-01 7.00000e-02 7.00000e-02 +2.77500 0.00000 0.00000 2.90000e-01 7.00000e-02 7.00000e-02 +2.82500 0.00000 0.00000 2.50000e-01 6.00000e-02 6.00000e-02 +2.87500 0.00000 0.00000 2.20000e-01 6.00000e-02 6.00000e-02 +2.92500 0.00000 0.00000 2.50000e-01 6.00000e-02 6.00000e-02 +2.97500 0.00000 0.00000 1.80000e-01 5.00000e-02 5.00000e-02 +3.02500 0.00000 0.00000 1.50000e-01 5.00000e-02 5.00000e-02 +3.07500 0.00000 0.00000 3.50000e-01 7.00000e-02 7.00000e-02 +3.12500 0.00000 0.00000 2.00000e-01 5.00000e-02 5.00000e-02 +3.17500 0.00000 0.00000 1.40000e-01 4.00000e-02 4.00000e-02 +3.22500 0.00000 0.00000 1.30000e-01 4.00000e-02 4.00000e-02 +3.27500 0.00000 0.00000 9.00000e-02 3.00000e-02 3.00000e-02 +3.32500 0.00000 0.00000 1.30000e-01 4.00000e-02 4.00000e-02 +3.37500 0.00000 0.00000 1.10000e-01 3.00000e-02 3.00000e-02 +3.42500 0.00000 0.00000 1.30000e-01 4.00000e-02 4.00000e-02 +3.47500 0.00000 0.00000 9.00000e-02 3.00000e-02 3.00000e-02 +3.52500 0.00000 0.00000 6.00000e-02 3.00000e-02 3.00000e-02 +3.57500 0.00000 0.00000 1.00000e-01 3.00000e-02 3.00000e-02 +3.62500 0.00000 0.00000 2.00000e-02 2.00000e-02 2.00000e-02 +3.67500 0.00000 0.00000 6.00000e-02 3.00000e-02 3.00000e-02 +3.72500 0.00000 0.00000 5.00000e-02 3.00000e-02 3.00000e-02 +3.77500 0.00000 0.00000 8.00000e-02 2.00000e-02 2.00000e-02 +3.82500 0.00000 0.00000 4.00000e-02 3.00000e-02 3.00000e-02 +3.87500 0.00000 0.00000 7.00000e-02 2.00000e-02 2.00000e-02 +3.92500 0.00000 0.00000 8.00000e-02 3.00000e-02 3.00000e-02 +3.97500 0.00000 0.00000 0.00000e+00 3.00000e-02 3.00000e-02 +4.02500 0.00000 0.00000 5.00000e-02 2.00000e-02 2.00000e-02 +4.07500 0.00000 0.00000 0.00000e+00 3.00000e-02 3.00000e-02 +4.12500 0.00000 0.00000 4.00000e-02 2.00000e-02 2.00000e-02 +4.17500 0.00000 0.00000 3.00000e-02 2.00000e-02 2.00000e-02 +4.22500 0.00000 0.00000 3.00000e-02 2.00000e-02 2.00000e-02 +4.27500 0.00000 0.00000 0.00000e+00 3.00000e-02 3.00000e-02 +4.32500 0.00000 0.00000 2.00000e-02 1.00000e-02 1.00000e-02 +END YODA_SCATTER2D_V2 diff --git a/analyses/pluginBELLE/BELLE_2017_I1613517.cc b/analyses/pluginBELLE/BELLE_2017_I1613517.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginBELLE/BELLE_2017_I1613517.cc @@ -0,0 +1,167 @@ +// -*- 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 BELLE_2017_I1613517 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(BELLE_2017_I1613517); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + + // Initialise and register projections + declare(FinalState(), "FS"); + declare(UnstableParticles(), "UFS"); + + // Book histograms + _c_DpDmS = bookCounter("/TMP/sigma_DpDmS"); + _c_DpSDmS = bookCounter("/TMP/sigma_DpSDmS"); + } + + 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"); + // total hadronic and muonic cross sections + map 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]) + vetoEvent; + // unstable charm analysis + const FinalState& ufs = apply(event, "UFS"); + for(unsigned int ix=0;ix nRes = nCount; + int ncount = ntotal; + findChildren(p1,nRes,ncount); + bool matched=false; + int sign = p1.pdgId()/id1; + // loop over the other fs particles + for(unsigned int iy=ix+1;iy nRes2 = nRes; + int ncount2 = ncount; + findChildren(p2,nRes2,ncount2); + if(ncount2!=0) continue; + matched=true; + for(auto const & val : nRes2) { + if(val.second!=0) { + matched = false; + break; + } + } + if(matched) { + if(id1==413 && id2==413) { + _c_DpSDmS->fill(event.weight()); + } + else if((id1==411 && id2==413) || + (id1==413 && id2==411)) { + _c_DpDmS->fill(event.weight()); + } + break; + } + } + if(matched) break; + } + } + + + /// Normalise histograms etc., after the run + void finalize() { + double fact = crossSection()/ sumOfWeights()/nanobarn; + for(unsigned int iy=1;iy<3;++iy) { + double sigma,error; + if(iy==1) { + sigma = _c_DpDmS->val()*fact; + error = _c_DpDmS->err()*fact; + } + else if(iy==2) { + sigma = _c_DpSDmS->val()*fact; + error = _c_DpSDmS->err()*fact; + } + Scatter2D temphisto(refData(1, 1, iy)); + Scatter2DPtr mult = bookScatter2D(1,1,iy); + 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 _c_DpDmS, _c_DpSDmS; + //@} + + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(BELLE_2017_I1613517); + + +} diff --git a/analyses/pluginBELLE/BELLE_2017_I1613517.info b/analyses/pluginBELLE/BELLE_2017_I1613517.info new file mode 100644 --- /dev/null +++ b/analyses/pluginBELLE/BELLE_2017_I1613517.info @@ -0,0 +1,35 @@ +Name: BELLE_2017_I1613517 +Year: 2017 +Summary: Cross Section for $e^+e^-\to D^\pm D^{*\mp}$ and $D^{*+}D^{*-}$ from threshold to 6 GeV +Experiment: BELLE +Collider: KEKB +InspireID: 1613517 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - Phys.Rev. D97 (2018) no.1, 012002 +RunInfo: e+ e- to hadrons +NeedCrossSection: yes +Beams: [e-, e+] +Description: + 'Measurement of the cross section for $e^+e^-\to D^\pm D^{*\mp}$ and $D^{*+}D^{*-}$ from threshold to 6 GeV.' +BibKey: Zhukova:2017pen +BibTeX: '@article{Zhukova:2017pen, + author = "Zhukova, V. and others", + title = "{Angular analysis of the $e^+ e^- \to D^{(*) \pm} D^{* + \mp}$ process near the open charm threshold using + initial-state radiation}", + collaboration = "Belle", + journal = "Phys. Rev.", + volume = "D97", + year = "2018", + number = "1", + pages = "012002", + doi = "10.1103/PhysRevD.97.012002", + eprint = "1707.09167", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + reportNumber = "BELLE-PREPRINT-2017-15, KEK-PREPRINT-2017-16", + SLACcitation = "%%CITATION = ARXIV:1707.09167;%%" +}' diff --git a/analyses/pluginBELLE/BELLE_2017_I1613517.plot b/analyses/pluginBELLE/BELLE_2017_I1613517.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginBELLE/BELLE_2017_I1613517.plot @@ -0,0 +1,14 @@ +BEGIN PLOT /BELLE_2017_I1613517/d01-x01-y01 +Title=$\sigma(e^+e^-\to D^+D^{*-} +\text{c.c})$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to D^+D^{*-} +\text{c.c} )/nb$ +LogY=0 +ConnectGaps=1 +END PLOT +BEGIN PLOT /BELLE_2017_I1613517/d01-x01-y02 +Title=$\sigma(e^+e^-\to D^{*+}D^{*-})$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to D^{*-}D^{*+})/nb$ +LogY=0 +ConnectGaps=1 +END PLOT diff --git a/analyses/pluginBELLE/BELLE_2017_I1613517.yoda b/analyses/pluginBELLE/BELLE_2017_I1613517.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginBELLE/BELLE_2017_I1613517.yoda @@ -0,0 +1,427 @@ +BEGIN YODA_SCATTER2D_V2 /REF/BELLE_2017_I1613517/d01-x01-y01 +IsRef: 1 +Path: /REF/BELLE_2017_I1613517/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +3.87500 0.00500 0.00500 0.01800 0.02800 0.02800 +3.88500 0.00500 0.00500 0.57900 0.31277 0.31277 +3.89500 0.00500 0.00500 0.89800 0.39763 0.39763 +3.90500 0.00500 0.00500 1.55600 0.45921 0.45921 +3.91500 0.00500 0.00500 1.47000 0.46329 0.46329 +3.92500 0.00500 0.00500 2.89300 0.53573 0.53573 +3.93500 0.00500 0.00500 1.65400 0.50324 0.50324 +3.94500 0.00500 0.00500 2.13500 0.50976 0.50976 +3.95500 0.00500 0.00500 2.61700 0.49578 0.49578 +3.96500 0.00500 0.00500 4.01400 0.56314 0.56314 +3.97500 0.00500 0.00500 2.79600 0.52050 0.52050 +3.98500 0.00500 0.00500 2.83700 0.53782 0.53782 +3.99500 0.00500 0.00500 3.03000 0.53881 0.53881 +4.00500 0.00500 0.00500 3.53200 0.53100 0.53100 +4.01500 0.00500 0.00500 4.42400 0.56820 0.56820 +4.02500 0.00500 0.00500 4.64200 0.56297 0.56297 +4.03500 0.00500 0.00500 3.95300 0.52737 0.52737 +4.04500 0.00500 0.00500 2.71900 0.48067 0.48067 +4.05500 0.00500 0.00500 2.37200 0.44608 0.44608 +4.06500 0.00500 0.00500 1.82500 0.40465 0.40465 +4.07500 0.00500 0.00500 1.46000 0.37964 0.37964 +4.08500 0.00500 0.00500 2.68100 0.39110 0.39110 +4.09500 0.00500 0.00500 1.62900 0.37459 0.37459 +4.10500 0.00500 0.00500 2.36400 0.36911 0.36911 +4.11500 0.00500 0.00500 1.76800 0.36535 0.36535 +4.12500 0.00500 0.00500 2.02200 0.38347 0.38347 +4.13500 0.00500 0.00500 0.82800 0.32450 0.32450 +4.14500 0.00500 0.00500 1.12200 0.32949 0.32949 +4.15500 0.00500 0.00500 1.49200 0.35745 0.35745 +4.16500 0.00500 0.00500 1.32700 0.34147 0.34147 +4.17500 0.00500 0.00500 1.08700 0.33885 0.33885 +4.18500 0.00500 0.00500 0.56200 0.33011 0.33011 +4.19500 0.00500 0.00500 1.03800 0.32150 0.32150 +4.20500 0.00500 0.00500 0.68300 0.29444 0.29444 +4.21500 0.00500 0.00500 0.79300 0.29034 0.29034 +4.22500 0.00500 0.00500 1.08100 0.32162 0.32162 +4.23500 0.00500 0.00500 0.85600 0.29743 0.29743 +4.24500 0.00500 0.00500 0.89400 0.29044 0.29044 +4.25500 0.00500 0.00500 0.46000 0.26748 0.26748 +4.26500 0.00500 0.00500 1.34400 0.29317 0.29317 +4.27500 0.00500 0.00500 0.64700 0.25424 0.25424 +4.28500 0.00500 0.00500 0.92000 0.26112 0.26112 +4.29500 0.00500 0.00500 0.45200 0.25245 0.25245 +4.30500 0.00500 0.00500 0.90300 0.25819 0.25819 +4.31500 0.00500 0.00500 0.96100 0.25612 0.25612 +4.32500 0.00500 0.00500 1.22800 0.27318 0.27318 +4.33500 0.00500 0.00500 0.64600 0.23721 0.23721 +4.34500 0.00500 0.00500 0.68800 0.25039 0.25039 +4.35500 0.00500 0.00500 1.07400 0.24910 0.24910 +4.36500 0.00500 0.00500 1.06400 0.23302 0.23302 +4.37500 0.00500 0.00500 0.82900 0.24529 0.24529 +4.38500 0.00500 0.00500 0.83700 0.24117 0.24117 +4.39500 0.00500 0.00500 0.85200 0.24221 0.24221 +4.40500 0.00500 0.00500 0.29900 0.21340 0.21340 +4.41500 0.00500 0.00500 0.97900 0.22002 0.22002 +4.42500 0.00500 0.00500 0.67000 0.22123 0.22123 +4.43500 0.00500 0.00500 0.52500 0.21119 0.21119 +4.44500 0.00500 0.00500 0.72800 0.21311 0.21311 +4.45500 0.00500 0.00500 0.34000 0.20735 0.20735 +4.46500 0.00500 0.00500 0.61200 0.19509 0.19509 +4.47500 0.00500 0.00500 0.74500 0.21112 0.21112 +4.48500 0.00500 0.00500 0.42300 0.19721 0.19721 +4.49500 0.00500 0.00500 0.24800 0.18527 0.18527 +4.50500 0.00500 0.00500 0.72500 0.20009 0.20009 +4.51500 0.00500 0.00500 0.56400 0.17904 0.17904 +4.52500 0.00500 0.00500 0.29500 0.17623 0.17623 +4.53500 0.00500 0.00500 0.53500 0.19217 0.19217 +4.54500 0.00500 0.00500 0.75700 0.19406 0.19406 +4.55500 0.00500 0.00500 0.43900 0.17410 0.17410 +4.56500 0.00500 0.00500 0.51900 0.18513 0.18513 +4.57500 0.00500 0.00500 0.46000 0.17914 0.17914 +4.58500 0.00500 0.00500 0.46400 0.17314 0.17314 +4.59500 0.00500 0.00500 0.59300 0.17005 0.17005 +4.60500 0.00500 0.00500 0.50600 0.17814 0.17814 +4.61500 0.00500 0.00500 0.35100 0.15711 0.15711 +4.62500 0.00500 0.00500 0.47400 0.15805 0.15805 +4.63500 0.00500 0.00500 0.55800 0.17107 0.17107 +4.64500 0.00500 0.00500 0.33300 0.14506 0.14506 +4.65500 0.00500 0.00500 0.72400 0.16200 0.16200 +4.66500 0.00500 0.00500 0.46000 0.15203 0.15203 +4.67500 0.00500 0.00500 0.46000 0.14705 0.14705 +4.68500 0.00500 0.00500 0.48500 0.15205 0.15205 +4.69500 0.00500 0.00500 0.36500 0.13603 0.13603 +4.70500 0.00500 0.00500 0.31200 0.13203 0.13203 +4.71500 0.00500 0.00500 0.41900 0.13701 0.13701 +4.72500 0.00500 0.00500 0.30300 0.13706 0.13706 +4.73500 0.00500 0.00500 0.32700 0.13809 0.13809 +4.74500 0.00500 0.00500 0.30100 0.12903 0.12903 +4.75500 0.00500 0.00500 0.43000 0.14805 0.14805 +4.76500 0.00500 0.00500 0.15500 0.11204 0.11204 +4.77500 0.00500 0.00500 0.48200 0.14403 0.14403 +4.78500 0.00500 0.00500 0.84600 0.15900 0.15900 +4.79500 0.00500 0.00500 0.42700 0.11900 0.11900 +4.80500 0.00500 0.00500 0.43900 0.14509 0.14509 +4.81500 0.00500 0.00500 0.31000 0.12906 0.12906 +4.82500 0.00500 0.00500 0.39700 0.12704 0.12704 +4.83500 0.00500 0.00500 0.62800 0.14601 0.14601 +4.84500 0.00500 0.00500 0.41600 0.12804 0.12804 +4.85500 0.00500 0.00500 0.44300 0.12500 0.12500 +4.86500 0.00500 0.00500 0.36800 0.11400 0.11400 +4.87500 0.00500 0.00500 0.31100 0.11402 0.11402 +4.88500 0.00500 0.00500 0.22600 0.10602 0.10602 +4.89500 0.00500 0.00500 0.18700 0.11007 0.11007 +4.90500 0.00500 0.00500 0.22500 0.10404 0.10404 +4.91500 0.00500 0.00500 0.33000 0.11804 0.11804 +4.92500 0.00500 0.00500 0.16900 0.08301 0.08301 +4.93500 0.00500 0.00500 0.13500 0.10408 0.10408 +4.94500 0.00500 0.00500 0.34100 0.11002 0.11002 +4.95500 0.00500 0.00500 0.35900 0.09800 0.09800 +4.96500 0.00500 0.00500 0.33400 0.11804 0.11804 +4.97500 0.00500 0.00500 0.33800 0.09800 0.09800 +4.98500 0.00500 0.00500 0.34700 0.10400 0.10400 +4.99500 0.00500 0.00500 0.19700 0.09002 0.09002 +5.00500 0.00500 0.00500 0.21500 0.09302 0.09302 +5.01500 0.00500 0.00500 0.33000 0.10100 0.10100 +5.02500 0.00500 0.00500 0.11300 0.07803 0.07803 +5.03500 0.00500 0.00500 0.28100 0.09902 0.09902 +5.04500 0.00500 0.00500 0.30000 0.09501 0.09501 +5.05500 0.00500 0.00500 0.28900 0.10302 0.10302 +5.06500 0.00500 0.00500 0.26400 0.09302 0.09302 +5.07500 0.00500 0.00500 0.21200 0.09302 0.09302 +5.08500 0.00500 0.00500 0.33400 0.09501 0.09501 +5.09500 0.00500 0.00500 0.33000 0.09801 0.09801 +5.10500 0.00500 0.00500 0.31100 0.09702 0.09702 +5.11500 0.00500 0.00500 0.28300 0.08200 0.08200 +5.12500 0.00500 0.00500 0.43800 0.10400 0.10400 +5.13500 0.00500 0.00500 0.36400 0.09801 0.09801 +5.14500 0.00500 0.00500 0.26200 0.08601 0.08601 +5.15500 0.00500 0.00500 0.16300 0.07301 0.07301 +5.16500 0.00500 0.00500 0.14200 0.07301 0.07301 +5.17500 0.00500 0.00500 0.02800 0.06303 0.06303 +5.18500 0.00500 0.00500 0.20400 0.07100 0.07100 +5.19500 0.00500 0.00500 0.15800 0.07603 0.07603 +5.20500 0.00500 0.00500 0.24000 0.08301 0.08301 +5.21500 0.00500 0.00500 0.15500 0.06400 0.06400 +5.22500 0.00500 0.00500 0.08100 0.06603 0.06603 +5.23500 0.00500 0.00500 0.22300 0.07801 0.07801 +5.24500 0.00500 0.00500 0.15200 0.07503 0.07503 +5.25500 0.00500 0.00500 0.22500 0.07701 0.07701 +5.26500 0.00500 0.00500 0.03800 0.06303 0.06303 +5.27500 0.00500 0.00500 0.06400 0.05401 0.05401 +5.28500 0.00500 0.00500 0.28600 0.08001 0.08001 +5.29500 0.00500 0.00500 0.18400 0.06800 0.06800 +5.30500 0.00500 0.00500 0.34200 0.08000 0.08000 +5.31500 0.00500 0.00500 0.25100 0.07601 0.07601 +5.32500 0.00500 0.00500 0.18500 0.07703 0.07703 +5.33500 0.00500 0.00500 0.26500 0.07801 0.07801 +5.34500 0.00500 0.00500 0.17800 0.05800 0.05800 +5.35500 0.00500 0.00500 0.08400 0.05801 0.05801 +5.36500 0.00500 0.00500 0.14300 0.05600 0.05600 +5.37500 0.00500 0.00500 0.07400 0.06003 0.06003 +5.38500 0.00500 0.00500 0.13500 0.06401 0.06401 +5.39500 0.00500 0.00500 0.13300 0.05701 0.05701 +5.40500 0.00500 0.00500 0.08700 0.05201 0.05201 +5.41500 0.00500 0.00500 0.10900 0.05801 0.05801 +5.42500 0.00500 0.00500 0.16600 0.05400 0.05400 +5.43500 0.00500 0.00500 0.16000 0.06201 0.06201 +5.44500 0.00500 0.00500 0.16900 0.05700 0.05700 +5.45500 0.00500 0.00500 0.16500 0.06001 0.06001 +5.46500 0.00500 0.00500 0.18000 0.05700 0.05700 +5.47500 0.00500 0.00500 0.17000 0.05700 0.05700 +5.48500 0.00500 0.00500 0.16800 0.05400 0.05400 +5.49500 0.00500 0.00500 0.08600 0.05704 0.05704 +5.50500 0.00500 0.00500 0.15800 0.05500 0.05500 +5.51500 0.00500 0.00500 0.21500 0.06701 0.06701 +5.52500 0.00500 0.00500 0.10000 0.04100 0.04100 +5.53500 0.00500 0.00500 0.12000 0.04200 0.04200 +5.54500 0.00500 0.00500 0.14800 0.05801 0.05801 +5.55500 0.00500 0.00500 0.14000 0.05100 0.05100 +5.56500 0.00500 0.00500 0.17300 0.05500 0.05500 +5.57500 0.00500 0.00500 0.11700 0.04900 0.04900 +5.58500 0.00500 0.00500 0.10300 0.05201 0.05201 +5.59500 0.00500 0.00500 0.07500 0.04501 0.04501 +5.60500 0.00500 0.00500 0.10400 0.05001 0.05001 +5.61500 0.00500 0.00500 0.20600 0.05900 0.05900 +5.62500 0.00500 0.00500 0.10000 0.04701 0.04701 +5.63500 0.00500 0.00500 0.19900 0.05600 0.05600 +5.64500 0.00500 0.00500 0.17300 0.05801 0.05801 +5.65500 0.00500 0.00500 0.07800 0.05001 0.05001 +5.66500 0.00500 0.00500 0.07500 0.04000 0.04000 +5.67500 0.00500 0.00500 0.09200 0.04601 0.04601 +5.68500 0.00500 0.00500 0.10400 0.04901 0.04901 +5.69500 0.00500 0.00500 0.17700 0.05401 0.05401 +5.70500 0.00500 0.00500 0.13000 0.04700 0.04700 +5.71500 0.00500 0.00500 0.12700 0.04300 0.04300 +5.72500 0.00500 0.00500 0.15700 0.04900 0.04900 +5.73500 0.00500 0.00500 0.02000 0.03002 0.03002 +5.74500 0.00500 0.00500 0.09700 0.04000 0.04000 +5.75500 0.00500 0.00500 0.18200 0.04700 0.04700 +5.76500 0.00500 0.00500 0.11100 0.04501 0.04501 +5.77500 0.00500 0.00500 0.07600 0.04101 0.04101 +5.78500 0.00500 0.00500 0.12400 0.04400 0.04400 +5.79500 0.00500 0.00500 0.04900 0.03200 0.03200 +5.80500 0.00500 0.00500 0.04300 0.03601 0.03601 +5.81500 0.00500 0.00500 0.11100 0.04301 0.04301 +5.82500 0.00500 0.00500 0.06400 0.03400 0.03400 +5.83500 0.00500 0.00500 0.16700 0.04500 0.04500 +5.84500 0.00500 0.00500 0.04800 0.03200 0.03200 +5.85500 0.00500 0.00500 0.04200 0.03401 0.03401 +5.86500 0.00500 0.00500 0.07200 0.03200 0.03200 +5.87500 0.00500 0.00500 0.04000 0.03202 0.03202 +5.88500 0.00500 0.00500 0.04600 0.03302 0.03302 +5.89500 0.00500 0.00500 0.07000 0.03100 0.03100 +5.90500 0.00500 0.00500 -0.01300 0.02502 0.02502 +5.91500 0.00500 0.00500 0.13500 0.04200 0.04200 +5.92500 0.00500 0.00500 0.07800 0.03601 0.03601 +5.93500 0.00500 0.00500 0.05300 0.03200 0.03200 +5.94500 0.00500 0.00500 0.10600 0.03700 0.03700 +5.95500 0.00500 0.00500 0.05900 0.03000 0.03000 +5.96500 0.00500 0.00500 0.04600 0.02800 0.02800 +5.97500 0.00500 0.00500 0.06000 0.03200 0.03200 +5.98500 0.00500 0.00500 0.07000 0.03300 0.03300 +5.99500 0.00500 0.00500 0.09800 0.03200 0.03200 +END YODA_SCATTER2D_V2 +BEGIN YODA_SCATTER2D_V2 /REF/BELLE_2017_I1613517/d01-x01-y02 +IsRef: 1 +Path: /REF/BELLE_2017_I1613517/d01-x01-y02 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +4.02500 0.00500 0.00500 0.78800 0.15799 0.15799 +4.03500 0.00500 0.00500 2.70200 0.26311 0.26311 +4.04500 0.00500 0.00500 2.95800 0.28778 0.28778 +4.05500 0.00500 0.00500 3.13700 0.28640 0.28640 +4.06500 0.00500 0.00500 2.93200 0.27716 0.27716 +4.07500 0.00500 0.00500 3.62900 0.31385 0.31385 +4.08500 0.00500 0.00500 3.51400 0.30438 0.30438 +4.09500 0.00500 0.00500 2.91900 0.27225 0.27225 +4.10500 0.00500 0.00500 2.57400 0.26565 0.26565 +4.11500 0.00500 0.00500 3.43600 0.29558 0.29558 +4.12500 0.00500 0.00500 3.10000 0.29788 0.29788 +4.13500 0.00500 0.00500 2.94300 0.28898 0.28898 +4.14500 0.00500 0.00500 3.10200 0.28562 0.28562 +4.15500 0.00500 0.00500 2.95600 0.28597 0.28597 +4.16500 0.00500 0.00500 3.17800 0.29954 0.29954 +4.17500 0.00500 0.00500 3.06000 0.29062 0.29062 +4.18500 0.00500 0.00500 2.78000 0.27427 0.27427 +4.19500 0.00500 0.00500 2.57000 0.25775 0.25775 +4.20500 0.00500 0.00500 2.22600 0.24042 0.24042 +4.21500 0.00500 0.00500 1.49500 0.20096 0.20096 +4.22500 0.00500 0.00500 1.10000 0.17293 0.17293 +4.23500 0.00500 0.00500 0.67400 0.13533 0.13533 +4.24500 0.00500 0.00500 0.30300 0.12249 0.12249 +4.25500 0.00500 0.00500 0.37400 0.12678 0.12678 +4.26500 0.00500 0.00500 0.61400 0.14303 0.14303 +4.27500 0.00500 0.00500 0.65000 0.14812 0.14812 +4.28500 0.00500 0.00500 0.74100 0.14983 0.14983 +4.29500 0.00500 0.00500 0.50700 0.12560 0.12560 +4.30500 0.00500 0.00500 0.89700 0.16867 0.16867 +4.31500 0.00500 0.00500 0.73500 0.14591 0.14591 +4.32500 0.00500 0.00500 0.95900 0.15377 0.15377 +4.33500 0.00500 0.00500 1.19100 0.17349 0.17349 +4.34500 0.00500 0.00500 1.07800 0.17051 0.17051 +4.35500 0.00500 0.00500 1.22100 0.17596 0.17596 +4.36500 0.00500 0.00500 1.12700 0.16813 0.16813 +4.37500 0.00500 0.00500 1.24100 0.17624 0.17624 +4.38500 0.00500 0.00500 1.06900 0.15893 0.15893 +4.39500 0.00500 0.00500 1.28500 0.17109 0.17109 +4.40500 0.00500 0.00500 0.91900 0.16130 0.16130 +4.41500 0.00500 0.00500 1.27200 0.17940 0.17940 +4.42500 0.00500 0.00500 0.83100 0.14504 0.14504 +4.43500 0.00500 0.00500 0.95100 0.14822 0.14822 +4.44500 0.00500 0.00500 0.67700 0.13200 0.13200 +4.45500 0.00500 0.00500 0.54500 0.12889 0.12889 +4.46500 0.00500 0.00500 0.72500 0.13147 0.13147 +4.47500 0.00500 0.00500 0.41200 0.10341 0.10341 +4.48500 0.00500 0.00500 0.44000 0.10850 0.10850 +4.49500 0.00500 0.00500 0.75400 0.12491 0.12491 +4.50500 0.00500 0.00500 0.63500 0.11790 0.11790 +4.51500 0.00500 0.00500 0.58400 0.12140 0.12140 +4.52500 0.00500 0.00500 0.57300 0.10651 0.10651 +4.53500 0.00500 0.00500 0.37800 0.10607 0.10607 +4.54500 0.00500 0.00500 0.65200 0.12202 0.12202 +4.55500 0.00500 0.00500 0.27500 0.09166 0.09166 +4.56500 0.00500 0.00500 0.56600 0.12317 0.12317 +4.57500 0.00500 0.00500 0.77600 0.12782 0.12782 +4.58500 0.00500 0.00500 0.45200 0.09767 0.09767 +4.59500 0.00500 0.00500 0.39600 0.09240 0.09240 +4.60500 0.00500 0.00500 0.32300 0.08896 0.08896 +4.61500 0.00500 0.00500 0.38700 0.09535 0.09535 +4.62500 0.00500 0.00500 0.50700 0.10296 0.10296 +4.63500 0.00500 0.00500 0.47900 0.09982 0.09982 +4.64500 0.00500 0.00500 0.43800 0.10439 0.10439 +4.65500 0.00500 0.00500 0.43400 0.10439 0.10439 +4.66500 0.00500 0.00500 0.45600 0.09866 0.09866 +4.67500 0.00500 0.00500 0.55100 0.09849 0.09849 +4.68500 0.00500 0.00500 0.35700 0.08911 0.08911 +4.69500 0.00500 0.00500 0.65300 0.10891 0.10891 +4.70500 0.00500 0.00500 0.39100 0.09420 0.09420 +4.71500 0.00500 0.00500 0.51700 0.09513 0.09513 +4.72500 0.00500 0.00500 0.48100 0.09964 0.09964 +4.73500 0.00500 0.00500 0.56500 0.10435 0.10435 +4.74500 0.00500 0.00500 0.48100 0.10259 0.10259 +4.75500 0.00500 0.00500 0.25400 0.07864 0.07864 +4.76500 0.00500 0.00500 0.25500 0.08161 0.08161 +4.77500 0.00500 0.00500 0.32000 0.07694 0.07694 +4.78500 0.00500 0.00500 0.49800 0.09276 0.09276 +4.79500 0.00500 0.00500 0.46500 0.08963 0.08963 +4.80500 0.00500 0.00500 0.39200 0.09108 0.09108 +4.81500 0.00500 0.00500 0.39900 0.08516 0.08516 +4.82500 0.00500 0.00500 0.23700 0.07344 0.07344 +4.83500 0.00500 0.00500 0.42600 0.08041 0.08041 +4.84500 0.00500 0.00500 0.37400 0.08302 0.08302 +4.85500 0.00500 0.00500 0.34100 0.07497 0.07497 +4.86500 0.00500 0.00500 0.22900 0.07344 0.07344 +4.87500 0.00500 0.00500 0.31900 0.07086 0.07086 +4.88500 0.00500 0.00500 0.15800 0.05423 0.05423 +4.89500 0.00500 0.00500 0.35700 0.07497 0.07497 +4.90500 0.00500 0.00500 0.38400 0.07415 0.07415 +4.91500 0.00500 0.00500 0.22200 0.07234 0.07234 +4.92500 0.00500 0.00500 0.27300 0.06859 0.06859 +4.93500 0.00500 0.00500 0.35100 0.07679 0.07679 +4.94500 0.00500 0.00500 0.32800 0.06972 0.06972 +4.95500 0.00500 0.00500 0.34300 0.08174 0.08174 +4.96500 0.00500 0.00500 0.29300 0.06661 0.06661 +4.97500 0.00500 0.00500 0.22900 0.06339 0.06339 +4.98500 0.00500 0.00500 0.25700 0.06152 0.06152 +4.99500 0.00500 0.00500 0.17000 0.05822 0.05822 +5.00500 0.00500 0.00500 0.23300 0.06140 0.06140 +5.01500 0.00500 0.00500 0.17200 0.05523 0.05523 +5.02500 0.00500 0.00500 0.26000 0.06041 0.06041 +5.03500 0.00500 0.00500 0.24400 0.06140 0.06140 +5.04500 0.00500 0.00500 0.28200 0.06549 0.06549 +5.05500 0.00500 0.00500 0.26100 0.06239 0.06239 +5.06500 0.00500 0.00500 0.18500 0.05423 0.05423 +5.07500 0.00500 0.00500 0.26500 0.06239 0.06239 +5.08500 0.00500 0.00500 0.16000 0.05315 0.05315 +5.09500 0.00500 0.00500 0.18200 0.06313 0.06313 +5.10500 0.00500 0.00500 0.17800 0.05515 0.05515 +5.11500 0.00500 0.00500 0.22100 0.05722 0.05722 +5.12500 0.00500 0.00500 0.25000 0.06129 0.06129 +5.13500 0.00500 0.00500 0.17000 0.05315 0.05315 +5.14500 0.00500 0.00500 0.18300 0.05515 0.05515 +5.15500 0.00500 0.00500 0.04200 0.04201 0.04201 +5.16500 0.00500 0.00500 0.10800 0.04305 0.04305 +5.17500 0.00500 0.00500 0.24500 0.05324 0.05324 +5.18500 0.00500 0.00500 0.22000 0.05415 0.05415 +5.19500 0.00500 0.00500 0.08000 0.03606 0.03606 +5.20500 0.00500 0.00500 0.19600 0.05016 0.05016 +5.21500 0.00500 0.00500 0.07000 0.03801 0.03801 +5.22500 0.00500 0.00500 0.12100 0.04604 0.04604 +5.23500 0.00500 0.00500 0.16600 0.04410 0.04410 +5.24500 0.00500 0.00500 0.14900 0.05104 0.05104 +5.25500 0.00500 0.00500 0.12600 0.04205 0.04205 +5.26500 0.00500 0.00500 0.10800 0.04205 0.04205 +5.27500 0.00500 0.00500 0.19700 0.04909 0.04909 +5.28500 0.00500 0.00500 0.14200 0.04804 0.04804 +5.29500 0.00500 0.00500 0.14500 0.04704 0.04704 +5.30500 0.00500 0.00500 0.07000 0.03202 0.03202 +5.31500 0.00500 0.00500 0.16900 0.04504 0.04504 +5.32500 0.00500 0.00500 0.15300 0.04305 0.04305 +5.33500 0.00500 0.00500 0.19900 0.04704 0.04704 +5.34500 0.00500 0.00500 0.17900 0.04405 0.04405 +5.35500 0.00500 0.00500 0.10400 0.03601 0.03601 +5.36500 0.00500 0.00500 0.24500 0.05104 0.05104 +5.37500 0.00500 0.00500 0.20500 0.05004 0.05004 +5.38500 0.00500 0.00500 0.20100 0.04704 0.04704 +5.39500 0.00500 0.00500 0.11200 0.03701 0.03701 +5.40500 0.00500 0.00500 0.15600 0.04101 0.04101 +5.41500 0.00500 0.00500 0.14700 0.04401 0.04401 +5.42500 0.00500 0.00500 0.13400 0.03901 0.03901 +5.43500 0.00500 0.00500 0.13300 0.03801 0.03801 +5.44500 0.00500 0.00500 0.15200 0.04101 0.04101 +5.45500 0.00500 0.00500 0.16500 0.04101 0.04101 +5.46500 0.00500 0.00500 0.06600 0.02900 0.02900 +5.47500 0.00500 0.00500 0.13800 0.03800 0.03800 +5.48500 0.00500 0.00500 0.19200 0.04200 0.04200 +5.49500 0.00500 0.00500 0.18100 0.04300 0.04300 +5.50500 0.00500 0.00500 0.15800 0.03900 0.03900 +5.51500 0.00500 0.00500 0.11000 0.03700 0.03700 +5.52500 0.00500 0.00500 0.17600 0.04200 0.04200 +5.53500 0.00500 0.00500 0.16300 0.03900 0.03900 +5.54500 0.00500 0.00500 0.07200 0.02800 0.02800 +5.55500 0.00500 0.00500 0.07400 0.03000 0.03000 +5.56500 0.00500 0.00500 0.04200 0.02800 0.02800 +5.57500 0.00500 0.00500 0.08900 0.03200 0.03200 +5.58500 0.00500 0.00500 0.06900 0.03000 0.03000 +5.59500 0.00500 0.00500 0.15600 0.03801 0.03801 +5.60500 0.00500 0.00500 0.03900 0.02100 0.02100 +5.61500 0.00500 0.00500 0.06300 0.02600 0.02600 +5.62500 0.00500 0.00500 0.03500 0.02500 0.02500 +5.63500 0.00500 0.00500 0.10600 0.03701 0.03701 +5.64500 0.00500 0.00500 0.07700 0.02602 0.02602 +5.65500 0.00500 0.00500 0.09400 0.02902 0.02902 +5.66500 0.00500 0.00500 0.05100 0.02100 0.02100 +5.67500 0.00500 0.00500 0.05500 0.02602 0.02602 +5.68500 0.00500 0.00500 0.12300 0.03202 0.03202 +5.69500 0.00500 0.00500 0.09800 0.02802 0.02802 +5.70500 0.00500 0.00500 0.14200 0.03406 0.03406 +5.71500 0.00500 0.00500 0.10200 0.03002 0.03002 +5.72500 0.00500 0.00500 0.08500 0.02902 0.02902 +5.73500 0.00500 0.00500 0.10200 0.02802 0.02802 +5.74500 0.00500 0.00500 0.05900 0.02502 0.02502 +5.75500 0.00500 0.00500 0.09700 0.02802 0.02802 +5.76500 0.00500 0.00500 0.08800 0.02902 0.02902 +5.77500 0.00500 0.00500 0.07100 0.02302 0.02302 +5.78500 0.00500 0.00500 0.05500 0.02502 0.02502 +5.79500 0.00500 0.00500 0.07800 0.02502 0.02502 +5.80500 0.00500 0.00500 0.10300 0.02807 0.02807 +5.81500 0.00500 0.00500 0.08500 0.02508 0.02508 +5.82500 0.00500 0.00500 0.08800 0.02608 0.02608 +5.83500 0.00500 0.00500 0.11400 0.02807 0.02807 +5.84500 0.00500 0.00500 0.05800 0.02102 0.02102 +5.85500 0.00500 0.00500 0.09500 0.02907 0.02907 +5.86500 0.00500 0.00500 0.04100 0.01903 0.01903 +5.87500 0.00500 0.00500 0.04100 0.01803 0.01803 +5.88500 0.00500 0.00500 0.10300 0.02707 0.02707 +5.89500 0.00500 0.00500 0.07000 0.02209 0.02209 +5.90500 0.00500 0.00500 0.08000 0.02309 0.02309 +5.91500 0.00500 0.00500 0.07300 0.02209 0.02209 +5.92500 0.00500 0.00500 0.07900 0.02309 0.02309 +5.93500 0.00500 0.00500 0.06500 0.02010 0.02010 +5.94500 0.00500 0.00500 0.03000 0.01603 0.01603 +5.95500 0.00500 0.00500 0.05600 0.02302 0.02302 +5.96500 0.00500 0.00500 0.10700 0.02518 0.02518 +5.97500 0.00500 0.00500 0.05600 0.01910 0.01910 +5.98500 0.00500 0.00500 0.04200 0.01703 0.01703 +5.99500 0.00500 0.00500 0.03800 0.01503 0.01503 +END YODA_SCATTER2D_V2