diff --git a/analyses/pluginBES/BESIII_2016_I1495838.cc b/analyses/pluginBES/BESIII_2016_I1495838.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2016_I1495838.cc @@ -0,0 +1,119 @@ +// -*- 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 BESIII_2016_I1495838 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(BESIII_2016_I1495838); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + declare(FinalState(), "FS"); + declare(UnstableParticles(), "UFS"); + _nPsi = bookCounter("TMP/psi"); + } + + 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 psi + if(p.pdgId()==443) { + map nRes = nCount; + int ncount = ntotal; + findChildren(p,nRes,ncount); + // psi 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) { + _nPsi->fill(event.weight()); + break; + } + } + } + } + + + /// Normalise histograms etc., after the run + void finalize() { + double fact = crossSection()/ sumOfWeights() /picobarn; + double sigma = _nPsi->val()*fact; + double error = _nPsi->err()*fact; + for(unsigned int ix=1;ix<3;++ix) { + Scatter2D temphisto(refData(ix, 1, 1)); + Scatter2DPtr mult = bookScatter2D(ix, 1, 1); + for (size_t b = 0; b < temphisto.numPoints(); b++) { + const double x = temphisto.point(b).x(); + pair ex = temphisto.point(b).xErrs(); + pair ex2 = ex; + if(ex2.first ==0.) ex2. first=0.0001; + if(ex2.second==0.) ex2.second=0.0001; + if (inRange(sqrtS()/GeV, x-ex2.first, x+ex2.second)) { + mult->addPoint(x, sigma, ex, make_pair(error,error)); + } + else { + mult->addPoint(x, 0., ex, make_pair(0.,.0)); + } + } + } + } + + //@} + + + /// @name Histograms + //@{ + CounterPtr _nPsi; + //@} + + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(BESIII_2016_I1495838); + + +} diff --git a/analyses/pluginBES/BESIII_2016_I1495838.info b/analyses/pluginBES/BESIII_2016_I1495838.info new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2016_I1495838.info @@ -0,0 +1,36 @@ +Name: BESIII_2016_I1495838 +Year: 2016 +Summary: Cross section for $e^+e^-\to J/\psi\pi^+\pi^-$ at energies between 3.77 and 4.6 GeV +Experiment: BESIII +Collider: BEPC +InspireID: 1495838 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - Phys.Rev.Lett. 118 (2017) no.9, 092001 +RunInfo: e+e- to hadrons +NeedCrossSection: yes +Beams: [e+,e-] +Luminosity_fb: +Description: + 'Measurement of the cross section for $e^+e^-\to J/\psi\pi^+\pi^-$ at energies between 3.77 and 4.6 GeV.' +Keywords: [] +BibKey: Ablikim:2016qzw +BibTeX: '@article{Ablikim:2016qzw, + author = "Ablikim, Medina and others", + title = "{Precise measurement of the $e^+e^-\to \pi^+\pi^-J/\psi$ + cross section at center-of-mass energies from 3.77 to 4.60 + GeV}", + collaboration = "BESIII", + journal = "Phys. Rev. Lett.", + volume = "118", + year = "2017", + number = "9", + pages = "092001", + doi = "10.1103/PhysRevLett.118.092001", + eprint = "1611.01317", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + SLACcitation = "%%CITATION = ARXIV:1611.01317;%%" +}' diff --git a/analyses/pluginBES/BESIII_2016_I1495838.plot b/analyses/pluginBES/BESIII_2016_I1495838.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2016_I1495838.plot @@ -0,0 +1,12 @@ +BEGIN PLOT /BESIII_2016_I1495838/d01-x01-y01 +Title=$\sigma(e^+e^-\to \pi^+\pi^-J/\psi)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to \pi^+\pi^-J/\psi)$/pb +ConnectGaps=1 +END PLOT +BEGIN PLOT /BESIII_2016_I1495838/d02-x01-y01 +Title=$\sigma(e^+e^-\to \pi^+\pi^-J/\psi)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to \pi^+\pi^-J/\psi)$/pb +ConnectGaps=1 +END PLOT \ No newline at end of file diff --git a/analyses/pluginBES/BESIII_2016_I1495838.yoda b/analyses/pluginBES/BESIII_2016_I1495838.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2016_I1495838.yoda @@ -0,0 +1,138 @@ +BEGIN YODA_SCATTER2D_V2 /REF/BESIII_2016_I1495838/d01-x01-y01 +IsRef: 1 +Path: /REF/BESIII_2016_I1495838/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +3.77300 0.00000 0.00000 2.85000e+01 1.80278e+00 1.80278e+00 +3.80770 0.00000 0.00000 1.67000e+01 3.44819e+00 3.44819e+00 +3.89620 0.00000 0.00000 1.71000e+01 3.54401e+00 3.54401e+00 +4.00760 0.00000 0.00000 1.60000e+01 1.48661e+00 1.48661e+00 +4.08550 0.00000 0.00000 1.50000e+01 3.22800e+00 3.22800e+00 +4.18860 0.00000 0.00000 1.55000e+01 3.90512e+00 3.90512e+00 +4.20770 0.00000 0.00000 5.34000e+01 6.22656e+00 6.22656e+00 +4.21710 0.00000 0.00000 6.03000e+01 6.68880e+00 6.68880e+00 +4.22630 0.00000 0.00000 8.51000e+01 5.12445e+00 5.12445e+00 +4.24170 0.00000 0.00000 8.44000e+01 7.98123e+00 7.98123e+00 +4.25800 0.00000 0.00000 5.95000e+01 3.67696e+00 3.67696e+00 +4.30790 0.00000 0.00000 5.20000e+01 6.44127e+00 6.44127e+00 +4.35830 0.00000 0.00000 2.54000e+01 1.92094e+00 1.92094e+00 +4.38740 0.00000 0.00000 2.00000e+01 3.41760e+00 3.41760e+00 +4.41560 0.00000 0.00000 1.21000e+01 9.21954e-01 9.21954e-01 +4.46710 0.00000 0.00000 1.33000e+01 2.24722e+00 2.24722e+00 +4.52710 0.00000 0.00000 1.06000e+01 1.99249e+00 1.99249e+00 +4.57450 0.00000 0.00000 1.34000e+01 3.29848e+00 3.29848e+00 +4.59950 0.00000 0.00000 6.40000e+00 8.06226e-01 8.06226e-01 +END YODA_SCATTER2D_V2 +BEGIN YODA_SCATTER2D_V2 /REF/BESIII_2016_I1495838/d02-x01-y01 +IsRef: 1 +Path: /REF/BESIII_2016_I1495838/d02-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +3.88740 0.00000 0.00000 3.88740e+00 1.31137e+01 9.11976e+00 +3.89240 0.00000 0.00000 3.89240e+00 1.38232e+01 9.83260e+00 +3.89740 0.00000 0.00000 3.89740e+00 1.34536e+01 9.37710e+00 +3.90240 0.00000 0.00000 3.90240e+00 1.42425e+01 1.02591e+01 +3.90740 0.00000 0.00000 3.90740e+00 1.28316e+01 8.54751e+00 +3.91240 0.00000 0.00000 3.91240e+00 1.34183e+01 9.22659e+00 +3.91740 0.00000 0.00000 3.91740e+00 1.14018e+01 6.60303e+00 +3.92240 0.00000 0.00000 3.92240e+00 1.71747e+01 1.27012e+01 +3.92740 0.00000 0.00000 3.92740e+00 1.56627e+01 1.11879e+01 +3.93240 0.00000 0.00000 3.93240e+00 1.24064e+01 8.10987e+00 +3.93740 0.00000 0.00000 3.93740e+00 1.27252e+01 8.53756e+00 +3.94240 0.00000 0.00000 3.94240e+00 1.26396e+01 8.65794e+00 +3.94740 0.00000 0.00000 3.94740e+00 1.48570e+01 1.10766e+01 +3.95240 0.00000 0.00000 3.95240e+00 1.30384e+01 9.35361e+00 +3.95740 0.00000 0.00000 3.95740e+00 1.39517e+01 1.04690e+01 +3.96240 0.00000 0.00000 3.96240e+00 1.23329e+01 8.44808e+00 +3.96740 0.00000 0.00000 3.96740e+00 1.32242e+01 9.13510e+00 +3.97240 0.00000 0.00000 3.97240e+00 1.27142e+01 9.01998e+00 +3.97740 0.00000 0.00000 3.97740e+00 1.17107e+01 7.91581e+00 +3.98240 0.00000 0.00000 3.98240e+00 1.45774e+01 1.09037e+01 +3.98740 0.00000 0.00000 3.98740e+00 8.42140e+00 1.21149e+01 +3.99240 0.00000 0.00000 3.99240e+00 6.60076e+00 1.05005e+01 +3.99740 0.00000 0.00000 3.99740e+00 8.96772e+00 1.27475e+01 +4.00240 0.00000 0.00000 4.00240e+00 1.10653e+01 1.49482e+01 +4.00740 0.00000 0.00000 4.00740e+00 1.43503e+01 1.02703e+01 +4.00940 0.00000 0.00000 4.00940e+00 1.33135e+01 8.92020e+00 +4.01140 0.00000 0.00000 4.01140e+00 1.53639e+01 1.09895e+01 +4.01340 0.00000 0.00000 4.01340e+00 1.38232e+01 9.23472e+00 +4.01540 0.00000 0.00000 4.01540e+00 1.36297e+01 9.34345e+00 +4.01740 0.00000 0.00000 4.01740e+00 1.73277e+01 1.31685e+01 +4.02240 0.00000 0.00000 4.02240e+00 1.67084e+01 1.23471e+01 +4.02740 0.00000 0.00000 4.02740e+00 7.71622e+00 6.02080e+00 +4.03240 0.00000 0.00000 4.03240e+00 1.52555e+01 1.09772e+01 +4.03740 0.00000 0.00000 4.03740e+00 1.17004e+01 7.00071e+00 +4.04740 0.00000 0.00000 4.04740e+00 1.26004e+01 8.00062e+00 +4.05240 0.00000 0.00000 4.05240e+00 1.44679e+01 1.02956e+01 +4.05740 0.00000 0.00000 4.05740e+00 1.39291e+01 9.24392e+00 +4.06240 0.00000 0.00000 4.06240e+00 1.34239e+01 9.23472e+00 +4.06740 0.00000 0.00000 4.06740e+00 1.23146e+01 8.22192e+00 +4.07740 0.00000 0.00000 4.07740e+00 1.36444e+01 9.96092e+00 +4.08740 0.00000 0.00000 4.08740e+00 1.29190e+01 9.12688e+00 +4.09740 0.00000 0.00000 4.09740e+00 1.17068e+01 7.61052e+00 +4.10740 0.00000 0.00000 4.10740e+00 1.26143e+01 8.52115e+00 +4.11740 0.00000 0.00000 4.11740e+00 1.12071e+01 7.31095e+00 +4.12740 0.00000 0.00000 4.12740e+00 1.27142e+01 8.52115e+00 +4.13740 0.00000 0.00000 4.13740e+00 1.51954e+01 1.12294e+01 +4.14240 0.00000 0.00000 4.14240e+00 1.25196e+01 8.62844e+00 +4.14740 0.00000 0.00000 4.14740e+00 1.14158e+01 7.32462e+00 +4.15740 0.00000 0.00000 4.15740e+00 1.55929e+01 1.19218e+01 +4.16740 0.00000 0.00000 4.16740e+00 6.51230e+00 4.81664e+00 +4.17740 0.00000 0.00000 4.17740e+00 1.44779e+01 1.03097e+01 +4.18740 0.00000 0.00000 4.18740e+00 1.12040e+01 7.00643e+00 +4.19240 0.00000 0.00000 4.19240e+00 1.46874e+01 1.07201e+01 +4.19740 0.00000 0.00000 4.19740e+00 1.56285e+01 1.16726e+01 +4.20040 0.00000 0.00000 4.20040e+00 2.00960e+01 1.58493e+01 +4.20340 0.00000 0.00000 4.20340e+00 1.59706e+01 1.19942e+01 +4.20740 0.00000 0.00000 4.20740e+00 1.51954e+01 1.12294e+01 +4.21240 0.00000 0.00000 4.21240e+00 2.02000e+01 1.65895e+01 +4.21740 0.00000 0.00000 4.21740e+00 1.98479e+01 1.63248e+01 +4.22240 0.00000 0.00000 4.22240e+00 2.05915e+01 1.73081e+01 +4.22740 0.00000 0.00000 4.22740e+00 2.40052e+01 2.09745e+01 +4.23240 0.00000 0.00000 4.23240e+00 1.86344e+01 1.55242e+01 +4.23740 0.00000 0.00000 4.23740e+00 2.21721e+01 1.89114e+01 +4.24040 0.00000 0.00000 4.24040e+00 1.87984e+01 1.57965e+01 +4.24240 0.00000 0.00000 4.24240e+00 1.57334e+01 1.25929e+01 +4.24540 0.00000 0.00000 4.24540e+00 1.68671e+01 1.39463e+01 +4.24740 0.00000 0.00000 4.24740e+00 1.70294e+01 1.40246e+01 +4.25240 0.00000 0.00000 4.25240e+00 1.49459e+01 1.20075e+01 +4.25740 0.00000 0.00000 4.25740e+00 1.76570e+01 1.49616e+01 +4.26240 0.00000 0.00000 4.26240e+00 1.62595e+01 1.33405e+01 +4.26740 0.00000 0.00000 4.26740e+00 1.77539e+01 1.49616e+01 +4.27240 0.00000 0.00000 4.27240e+00 1.62975e+01 1.33645e+01 +4.27740 0.00000 0.00000 4.27740e+00 1.42706e+01 1.12178e+01 +4.28240 0.00000 0.00000 4.28240e+00 1.69650e+01 1.40431e+01 +4.28740 0.00000 0.00000 4.28740e+00 1.60854e+01 1.32699e+01 +4.29740 0.00000 0.00000 4.29740e+00 1.44257e+01 1.12614e+01 +4.30740 0.00000 0.00000 4.30740e+00 1.68122e+01 1.38051e+01 +4.31740 0.00000 0.00000 4.31740e+00 1.34974e+01 1.06513e+01 +4.32740 0.00000 0.00000 4.32740e+00 1.32966e+01 1.01272e+01 +4.33740 0.00000 0.00000 4.33740e+00 1.34213e+01 1.03576e+01 +4.34740 0.00000 0.00000 4.34740e+00 1.14280e+01 8.23893e+00 +4.35740 0.00000 0.00000 4.35740e+00 1.49626e+01 1.18068e+01 +4.36740 0.00000 0.00000 4.36740e+00 1.38452e+01 1.07870e+01 +4.37740 0.00000 0.00000 4.37740e+00 1.22409e+01 8.95600e+00 +4.38740 0.00000 0.00000 4.38740e+00 1.32544e+01 9.67471e+00 +4.39240 0.00000 0.00000 4.39240e+00 1.32740e+01 9.60260e+00 +4.39740 0.00000 0.00000 4.39740e+00 1.21413e+01 8.26075e+00 +4.40740 0.00000 0.00000 4.40740e+00 1.10041e+01 6.20725e+00 +4.41740 0.00000 0.00000 4.41740e+00 1.23406e+01 8.65794e+00 +4.42240 0.00000 0.00000 4.42240e+00 1.24487e+01 8.67006e+00 +4.42740 0.00000 0.00000 4.42740e+00 1.19151e+01 7.62365e+00 +4.43740 0.00000 0.00000 4.43740e+00 1.21499e+01 8.47172e+00 +4.44740 0.00000 0.00000 4.44740e+00 1.02020e+01 6.40312e+00 +4.45740 0.00000 0.00000 4.45740e+00 9.40851e+00 6.11310e+00 +4.47740 0.00000 0.00000 4.47740e+00 1.12219e+01 7.73175e+00 +4.49740 0.00000 0.00000 4.49740e+00 8.30060e+00 4.30116e+00 +4.51740 0.00000 0.00000 4.51740e+00 1.02240e+01 6.73647e+00 +4.53740 0.00000 0.00000 4.53740e+00 1.06301e+01 7.54255e+00 +4.54740 0.00000 0.00000 4.54740e+00 1.08374e+01 7.45453e+00 +4.55740 0.00000 0.00000 4.55740e+00 1.00045e+01 6.20725e+00 +4.56740 0.00000 0.00000 4.56740e+00 1.06118e+01 6.81836e+00 +4.57740 0.00000 0.00000 4.57740e+00 1.11113e+01 7.51665e+00 +4.58740 0.00000 0.00000 4.58740e+00 8.70057e+00 4.40114e+00 +END YODA_SCATTER2D_V2 diff --git a/analyses/pluginBES/BESIII_2018_I1627871.cc b/analyses/pluginBES/BESIII_2018_I1627871.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2018_I1627871.cc @@ -0,0 +1,136 @@ +// -*- 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 BESIII_2018_I1627871 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(BESIII_2018_I1627871); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + // Initialise and register projections + declare(FinalState(), "FS"); + declare(UnstableParticles(), "UFS"); + _nLambda= bookCounter( "/TMP/nLambda" ); + } + + 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 + 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) { + _nLambda->fill(event.weight()); + break; + } + } + if(matched) break; + } + } + + + /// Normalise histograms etc., after the run + void finalize() { + double sigma = _nLambda->val(); + double error = _nLambda->err(); + sigma *= crossSection()/ sumOfWeights() /picobarn; + error *= crossSection()/ sumOfWeights() /picobarn; + Scatter2D temphisto(refData(1, 1, 1)); + Scatter2DPtr mult = bookScatter2D(1, 1, 1); + for (size_t b = 0; b < temphisto.numPoints(); b++) { + const double x = temphisto.point(b).x(); + pair ex = temphisto.point(b).xErrs(); + pair ex2 = ex; + if(ex2.first ==0.) ex2. first=0.0001; + if(ex2.second==0.) ex2.second=0.0001; + if (inRange(sqrtS()/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 _nLambda; + //@} + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(BESIII_2018_I1627871); + + +} diff --git a/analyses/pluginBES/BESIII_2018_I1627871.info b/analyses/pluginBES/BESIII_2018_I1627871.info new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2018_I1627871.info @@ -0,0 +1,34 @@ +Name: BESIII_2018_I1627871 +Year: 2018 +Summary: Cross section for $e^+e^-\to\Lambda\bar{\Lambda}$ between threshold and 3.08 GeV +Experiment: BESIII +Collider: BEPC II +InspireID: 1627871 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - Phys.Rev. D97 (2018) no.3, 032013 +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.08 GeV.' +Keywords: +BibKey: Ablikim:2017pyl +BibTeX: '@article{Ablikim:2017pyl, + author = "Ablikim, Medina and others", + title = "{Observation of a cross-section enhancement near mass + threshold in $e^{+}e^{-}\rightarrow\Lambda\bar{\Lambda}$}", + collaboration = "BESIII", + journal = "Phys. Rev.", + volume = "D97", + year = "2018", + number = "3", + pages = "032013", + doi = "10.1103/PhysRevD.97.032013", + eprint = "1709.10236", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + SLACcitation = "%%CITATION = ARXIV:1709.10236;%%" +}' diff --git a/analyses/pluginBES/BESIII_2018_I1627871.plot b/analyses/pluginBES/BESIII_2018_I1627871.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2018_I1627871.plot @@ -0,0 +1,6 @@ +BEGIN PLOT /BESIII_2018_I1627871/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 diff --git a/analyses/pluginBES/BESIII_2018_I1627871.yoda b/analyses/pluginBES/BESIII_2018_I1627871.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2018_I1627871.yoda @@ -0,0 +1,12 @@ +BEGIN YODA_SCATTER2D_V2 /REF/BESIII_2018_I1627871/d01-x01-y01 +IsRef: 1 +Path: /REF/BESIII_2018_I1627871/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +2.23240 0.00000 0.00000 3.05000e+02 5.76281e+01 7.98812e+01 +2.40000 0.00000 0.00000 1.28000e+02 2.61725e+01 2.61725e+01 +2.80000 0.00000 0.00000 1.48000e+01 5.53624e+00 5.53624e+00 +3.08000 0.00000 0.00000 4.20000e+00 1.30000e+00 1.30000e+00 +END YODA_SCATTER2D_V2 diff --git a/analyses/pluginBES/BESIII_2018_I1691798.cc b/analyses/pluginBES/BESIII_2018_I1691798.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2018_I1691798.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 BESIII_2018_I1691798 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(BESIII_2018_I1691798); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + // Initialise and register projections + declare(FinalState(), "FS"); + _nKKpi= bookCounter( "/TMP/nKKpi" ); + } + + + /// Perform the per-event analysis + void analyze(const Event& event) { + const FinalState& fs = apply(event, "FS"); + + map nCount; + int ntotal(0); + foreach (const Particle& p, fs.particles()) { + nCount[p.pdgId()] += 1; + ++ntotal; + } + + if(ntotal==3 && nCount[310]==1 && + ((nCount[ 321]=1 && nCount[-211] ==1) || + (nCount[-321]=1 && nCount[ 211] ==1))) + _nKKpi->fill(event.weight()); + } + + + /// Normalise histograms etc., after the run + void finalize() { + double sigma = _nKKpi->val(); + double error = _nKKpi->err(); + sigma *= crossSection()/ sumOfWeights() /picobarn; + error *= crossSection()/ sumOfWeights() /picobarn; + Scatter2D temphisto(refData(1, 1, 1)); + Scatter2DPtr mult = bookScatter2D(1, 1, 1); + for (size_t b = 0; b < temphisto.numPoints(); b++) { + const double x = temphisto.point(b).x(); + pair ex = temphisto.point(b).xErrs(); + pair ex2 = ex; + if(ex2.first ==0.) ex2. first=0.0001; + if(ex2.second==0.) ex2.second=0.0001; + if (inRange(sqrtS()/GeV, x-ex2.first, x+ex2.second)) { + mult->addPoint(x, sigma, ex, make_pair(error,error)); + } + else { + mult->addPoint(x, 0., ex, make_pair(0.,.0)); + } + } + } + + //@} + + + /// @name Histograms + //@{ + CounterPtr _nKKpi; + //@} + + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(BESIII_2018_I1691798); + + +} diff --git a/analyses/pluginBES/BESIII_2018_I1691798.info b/analyses/pluginBES/BESIII_2018_I1691798.info new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2018_I1691798.info @@ -0,0 +1,31 @@ +Name: BESIII_2018_I1691798 +Year: 2018 +Summary: Cross Section for $K^0_SK^\pm\pi^\mp$ between 3.8 and 4.6 GeV +Experiment: BESIII +Collider: BEPC II +InspireID: 1691798 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - arXiv:1808.08733 [hep-ex] +RunInfo: e+e- to hadrons +NeedCrossSection: yes +Beams: [e+,e-] +Luminosity_fb: +Description: + 'Measurement of the cross section for $K^0_SK^\pm\pi^\mp$ between 3.8 and 4.6 GeV.' +Keywords: +BibKey: Ablikim:2018lmx +BibTeX: '@article{Ablikim:2018lmx, + author = "Ablikim, Medina and others", + title = "{Precision measurements of $\sigma_{B}(e^+e^- \to + K_{S}^{0}K^{\pm}\pi^{\mp})$ at center-of-mass energies + between 3.8 and 4.6 GeV}", + collaboration = "BESIII", + year = "2018", + eprint = "1808.08733", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + SLACcitation = "%%CITATION = ARXIV:1808.08733;%%" +}' diff --git a/analyses/pluginBES/BESIII_2018_I1691798.plot b/analyses/pluginBES/BESIII_2018_I1691798.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2018_I1691798.plot @@ -0,0 +1,6 @@ +BEGIN PLOT /BESIII_2018_I1691798/d01-x01-y01 +Title=$\sigma(e^+e^-\to K^0_SK^\pm \pi^\mp)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to K^0_SK^\pm \pi^\mp)$/pb +ConnectGaps=1 +END PLOT diff --git a/analyses/pluginBES/BESIII_2018_I1691798.yoda b/analyses/pluginBES/BESIII_2018_I1691798.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2018_I1691798.yoda @@ -0,0 +1,23 @@ +BEGIN YODA_SCATTER2D_V2 /REF/BESIII_2018_I1691798/d01-x01-y01 +IsRef: 1 +Path: /REF/BESIII_2018_I1691798/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +3.80800 0.00000 0.00000 1.74000e+01 1.89568e+00 1.88053e+00 +3.89600 0.00000 0.00000 1.01100e+01 1.44568e+00 1.42299e+00 +4.00800 0.00000 0.00000 9.08000e+00 9.53362e-01 9.53362e-01 +4.08600 0.00000 0.00000 8.91000e+00 1.37971e+00 1.35680e+00 +4.18900 0.00000 0.00000 9.32000e+00 1.48741e+00 1.45588e+00 +4.20800 0.00000 0.00000 7.92000e+00 1.27342e+00 1.25841e+00 +4.21700 0.00000 0.00000 8.59000e+00 1.31913e+00 1.30419e+00 +4.22600 0.00000 0.00000 7.69000e+00 8.53464e-01 8.53464e-01 +4.24200 0.00000 0.00000 7.28000e+00 1.26507e+00 1.23434e+00 +4.25800 0.00000 0.00000 6.79000e+00 8.08022e-01 8.08022e-01 +4.30800 0.00000 0.00000 5.22000e+00 1.11104e+00 1.08743e+00 +4.35800 0.00000 0.00000 5.90000e+00 7.70000e-01 7.70000e-01 +4.38700 0.00000 0.00000 5.39000e+00 1.09659e+00 1.06588e+00 +4.41600 0.00000 0.00000 5.59000e+00 7.15821e-01 7.15821e-01 +4.60000 0.00000 0.00000 4.18000e+00 6.39062e-01 6.39062e-01 +END YODA_SCATTER2D_V2 diff --git a/analyses/pluginBES/BESIII_2018_I1699641.cc b/analyses/pluginBES/BESIII_2018_I1699641.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2018_I1699641.cc @@ -0,0 +1,138 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/Projections/FinalState.hh" +#include "Rivet/Projections/FastJets.hh" + +namespace Rivet { + + + /// @brief Add a short analysis description here + class BESIII_2018_I1699641 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(BESIII_2018_I1699641); + + + /// @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 + _cKKpipi = bookCounter("TMP/2Kpipi" ); + _cKKpieta = bookCounter("TMP/2Kpieta"); + } + + + 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; + } + // K K pi pi + if(ntotal==4 && nCount[310]==1 && nCount[111]==1 && + ((nCount[ 321]==1 &&nCount[-211]==1) || + (nCount[-321]==1 &&nCount[ 211]==1) )) + _cKKpipi->fill(event.weight()); + // eta resonance + const FinalState& ufs = apply(event, "UFS"); + foreach (const Particle& p, ufs.particles()) { + if(p.children().empty()) continue; + if(p.pdgId()!=221) continue; + map nRes=nCount; + int ncount = ntotal; + findChildren(p,nRes,ncount); + if(ncount!=3) continue; + bool matched=true; + for(auto const & val : nRes) { + if(abs(val.first)==321 || abs(val.first)==211) { + continue; + } + else if(abs(val.first)==310) { + if(val.second!=1) { + matched = false; + break; + } + } + else if(val.second!=0) { + matched = false; + break; + } + } + if(matched==false) continue; + if((nCount[ 321] == 1 && nCount[-211] ==1) || + (nCount[-321] == 1 && nCount[ 211] ==1)) + _cKKpieta->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 = _cKKpipi->val(); + error = _cKKpipi->err(); + } + else if(ix==2) { + sigma = _cKKpieta->val(); + error = _cKKpieta->err(); + } + sigma *= crossSection()/ sumOfWeights() /picobarn; + error *= crossSection()/ sumOfWeights() /picobarn; + 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 _cKKpipi,_cKKpieta; + //@} + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(BESIII_2018_I1699641); + + +} diff --git a/analyses/pluginBES/BESIII_2018_I1699641.info b/analyses/pluginBES/BESIII_2018_I1699641.info new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2018_I1699641.info @@ -0,0 +1,31 @@ +Name: BESIII_2018_I1699641 +Year: 2018 +Summary: Cross section for $e^+e^-\to K^0_SK^\pm\pi^\mp\pi^0$ and $K^0_SK^\pm\pi^\mp\eta$ between 3.90 to 4.60 GeV +Experiment: BESIII +Collider: BEPC II +InspireID: 1699641 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - arXiv:1810.09395 +RunInfo: e+e- to hadrons +NeedCrossSection: yes +Beams: [e+,e-] +Description: + 'Measurement of the cross section for $e^+e^-\to K^0_SK^\pm\pi^\mp\pi^0$ and $K^0_SK^\pm\pi^\mp\eta$ between 3.90 to 4.60 GeV.' +Keywords: +BibKey: Ablikim:2018ddb +BibTeX: '@article{Ablikim:2018ddb, + author = "Ablikim, Medina and others", + title = "{Measurements of $e^+e^- \to + K_{S}^{0}K^{\pm}\pi^{\mp}\pi^0$ and + $K_{S}^{0}K^{\pm}\pi^{\mp}\eta$ at center-of-mass energies + from $3.90$ to $4.60~\mathrm{GeV}$}", + collaboration = "BESIII", + year = "2018", + eprint = "1810.09395", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + SLACcitation = "%%CITATION = ARXIV:1810.09395;%%" +}' diff --git a/analyses/pluginBES/BESIII_2018_I1699641.plot b/analyses/pluginBES/BESIII_2018_I1699641.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2018_I1699641.plot @@ -0,0 +1,12 @@ +BEGIN PLOT /BESIII_2018_I1699641/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)$/pb +ConnectGaps=1 +END PLOT +BEGIN PLOT /BESIII_2018_I1699641/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)$/pb +ConnectGaps=1 +END PLOT diff --git a/analyses/pluginBES/BESIII_2018_I1699641.yoda b/analyses/pluginBES/BESIII_2018_I1699641.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2018_I1699641.yoda @@ -0,0 +1,50 @@ +BEGIN YODA_SCATTER2D_V2 /REF/BESIII_2018_I1699641/d01-x01-y01 +IsRef: 1 +Path: /REF/BESIII_2018_I1699641/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +3.89600 0.00000 0.00000 7.44100e+01 4.32908e+00 4.32908e+00 +4.00800 0.00000 0.00000 5.80200e+01 2.49802e+00 2.49802e+00 +4.08600 0.00000 0.00000 4.75200e+01 3.42533e+00 3.42533e+00 +4.18900 0.00000 0.00000 4.53800e+01 3.56702e+00 3.56702e+00 +4.20800 0.00000 0.00000 4.09500e+01 3.15119e+00 3.15119e+00 +4.21700 0.00000 0.00000 3.82800e+01 3.03315e+00 3.03315e+00 +4.22600 0.00000 0.00000 3.92300e+01 1.96023e+00 1.96023e+00 +4.24200 0.00000 0.00000 3.74600e+01 2.98371e+00 2.98371e+00 +4.25800 0.00000 0.00000 3.86500e+01 1.96898e+00 1.96898e+00 +4.30800 0.00000 0.00000 3.68600e+01 3.19130e+00 3.19130e+00 +4.35800 0.00000 0.00000 3.35300e+01 1.88106e+00 1.88106e+00 +4.38700 0.00000 0.00000 3.36800e+01 2.80713e+00 2.80713e+00 +4.41600 0.00000 0.00000 3.03800e+01 1.72630e+00 1.72630e+00 +4.46700 0.00000 0.00000 2.67000e+01 2.07470e+00 2.07470e+00 +4.52700 0.00000 0.00000 2.65100e+01 2.08327e+00 2.08327e+00 +4.57500 0.00000 0.00000 2.48700e+01 2.54629e+00 2.54629e+00 +4.60000 0.00000 0.00000 2.27100e+01 1.53782e+00 1.53782e+00 +END YODA_SCATTER2D_V2 +BEGIN YODA_SCATTER2D_V2 /REF/BESIII_2018_I1699641/d02-x01-y01 +IsRef: 1 +Path: /REF/BESIII_2018_I1699641/d02-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +3.89600 0.00000 0.00000 2.72300e+01 3.59004e+00 3.59004e+00 +4.00800 0.00000 0.00000 1.98800e+01 1.65118e+00 1.65118e+00 +4.08500 0.00000 0.00000 1.47100e+01 2.72075e+00 2.72075e+00 +4.18900 0.00000 0.00000 1.77500e+01 3.16735e+00 3.16735e+00 +4.20800 0.00000 0.00000 1.42000e+01 2.56049e+00 2.56049e+00 +4.21700 0.00000 0.00000 1.00500e+01 2.15000e+00 2.15000e+00 +4.22600 0.00000 0.00000 1.56100e+01 1.24101e+00 1.24101e+00 +4.24200 0.00000 0.00000 1.46300e+01 2.52951e+00 2.52951e+00 +4.25800 0.00000 0.00000 1.43500e+01 1.28254e+00 1.28254e+00 +4.30800 0.00000 0.00000 1.26700e+01 2.60085e+00 2.60085e+00 +4.35800 0.00000 0.00000 1.13800e+01 1.18507e+00 1.18507e+00 +4.38700 0.00000 0.00000 1.30500e+01 2.25424e+00 2.25424e+00 +4.41600 0.00000 0.00000 1.06200e+01 1.07931e+00 1.07931e+00 +4.46700 0.00000 0.00000 1.06200e+01 1.64441e+00 1.64441e+00 +4.52700 0.00000 0.00000 7.27000e+00 1.56719e+00 1.56719e+00 +4.57500 0.00000 0.00000 9.23000e+00 2.08940e+00 2.08940e+00 +4.60000 0.00000 0.00000 8.67000e+00 1.07313e+00 1.07313e+00 +END YODA_SCATTER2D_V2 diff --git a/analyses/pluginBES/BESIII_2018_I1704558.cc b/analyses/pluginBES/BESIII_2018_I1704558.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2018_I1704558.cc @@ -0,0 +1,81 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/Projections/FinalState.hh" + +namespace Rivet { + + + /// @brief Add a short analysis description here + class BESIII_2018_I1704558 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(BESIII_2018_I1704558); + + + /// @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() /picobarn; + error *= crossSection()/ sumOfWeights() /picobarn; + Scatter2D temphisto(refData(1, 1, 1)); + Scatter2DPtr mult = bookScatter2D(1, 1, 1); + for (size_t b = 0; b < temphisto.numPoints(); b++) { + const double x = temphisto.point(b).x(); + pair ex = temphisto.point(b).xErrs(); + pair ex2 = ex; + if(ex2.first ==0.) ex2. first=0.0001; + if(ex2.second==0.) ex2.second=0.0001; + if (inRange(sqrtS()/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(BESIII_2018_I1704558); + + +} diff --git a/analyses/pluginBES/BESIII_2018_I1704558.info b/analyses/pluginBES/BESIII_2018_I1704558.info new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2018_I1704558.info @@ -0,0 +1,29 @@ +Name: BESIII_2018_I1704558 +Year: 2018 +Summary: Cross section for $e^+e^-\to K^+K^-$ between 2.00 and 3.08 GeV +Experiment: BESIII +Collider: BEPC II +InspireID: 1704558 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - arXiv:1811.08742 +RunInfo: e+ e- to hadrons between 2.00 and 3.08 GeV +NeedCrossSection: yes +Beams: [e+, e-] +Description: + 'Measurement of the cross section for $e^+e^-\to K^+K^-$ at energies between 2.00 and 3.08 GeV. Useful for comparing models of the kaon form factor.' +Keywords: [] +BibKey: Ablikim:2018iyx +BibTeX: '@article{Ablikim:2018iyx, + author = "Ablikim, M. and others", + title = "{Measurement of $e^{+} e^{-} \rightarrow K^{+} K^{-}$ + cross section at $\sqrt{s} = 2.00 - 3.08$ GeV}", + collaboration = "BESIII", + year = "2018", + eprint = "1811.08742", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + SLACcitation = "%%CITATION = ARXIV:1811.08742;%%" +}' diff --git a/analyses/pluginBES/BESIII_2018_I1704558.plot b/analyses/pluginBES/BESIII_2018_I1704558.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2018_I1704558.plot @@ -0,0 +1,6 @@ +BEGIN PLOT /BESIII_2018_I1704558/d01-x01-y01 +Title=$\sigma(e^+e^-\to K^+K^-)$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to K^+K^-)$/pb +ConnectGaps=1 +END PLOT diff --git a/analyses/pluginBES/BESIII_2018_I1704558.yoda b/analyses/pluginBES/BESIII_2018_I1704558.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2018_I1704558.yoda @@ -0,0 +1,30 @@ +BEGIN YODA_SCATTER2D_V2 /REF/BESIII_2018_I1704558/d01-x01-y01 +IsRef: 1 +Path: /REF/BESIII_2018_I1704558/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +2.00000 0.00000 0.00000 3.51500e+02 9.23255e+00 9.23255e+00 +2.05000 0.00000 0.00000 2.96100e+02 1.36605e+01 1.36605e+01 +2.10000 0.00000 0.00000 2.20600e+02 6.76831e+00 6.76831e+00 +2.12500 0.00000 0.00000 1.92000e+02 3.55528e+00 3.55528e+00 +2.15000 0.00000 0.00000 1.71700e+02 1.10856e+01 1.10856e+01 +2.17500 0.00000 0.00000 1.84200e+02 6.45678e+00 6.45678e+00 +2.20000 0.00000 0.00000 2.31400e+02 6.67008e+00 6.67008e+00 +2.23240 0.00000 0.00000 2.53200e+02 7.24500e+00 7.24500e+00 +2.30940 0.00000 0.00000 1.84000e+02 5.05964e+00 5.05964e+00 +2.38640 0.00000 0.00000 1.01500e+02 3.74299e+00 3.74299e+00 +2.39600 0.00000 0.00000 1.02600e+02 2.94788e+00 2.94788e+00 +2.50000 0.00000 0.00000 8.81000e+01 1.24274e+01 1.24274e+01 +2.64440 0.00000 0.00000 5.66000e+01 2.80000e+00 2.80000e+00 +2.64640 0.00000 0.00000 5.67000e+01 2.61534e+00 2.61534e+00 +2.70000 0.00000 0.00000 3.63000e+01 8.55336e+00 8.55336e+00 +2.80000 0.00000 0.00000 3.79000e+01 8.99111e+00 8.99111e+00 +2.90000 0.00000 0.00000 3.04000e+01 1.90788e+00 1.90788e+00 +2.95000 0.00000 0.00000 2.53000e+01 2.53180e+00 2.53180e+00 +2.98100 0.00000 0.00000 2.80000e+01 2.40000e+00 2.40000e+00 +3.00000 0.00000 0.00000 2.44000e+01 2.49800e+00 2.49800e+00 +3.02000 0.00000 0.00000 2.48000e+01 2.49800e+00 2.49800e+00 +3.08000 0.00000 0.00000 2.53000e+01 2.21133e+00 2.21133e+00 +END YODA_SCATTER2D_V2 diff --git a/analyses/pluginBES/BESIII_2019_I1718337.cc b/analyses/pluginBES/BESIII_2019_I1718337.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2019_I1718337.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 BESIII_2019_I1718337 : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(BESIII_2019_I1718337); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + + // Initialise and register projections + declare(FinalState(), "FS"); + + // Book histograms + _nproton = bookCounter("TMP/proton"); + } + + + /// 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::PROTON) vetoEvent; + } + _nproton->fill(event.weight()); + } + + + /// Normalise histograms etc., after the run + void finalize() { + double sigma = _nproton->val(); + double error = _nproton->err(); + sigma *= crossSection()/ sumOfWeights() /picobarn; + error *= crossSection()/ sumOfWeights() /picobarn; + Scatter2D temphisto(refData(1, 1, 1)); + Scatter2DPtr mult = bookScatter2D(1, 1, 1); + for (size_t b = 0; b < temphisto.numPoints(); b++) { + const double x = temphisto.point(b).x(); + pair ex = temphisto.point(b).xErrs(); + pair ex2 = ex; + if(ex2.first ==0.) ex2. first=0.0001; + if(ex2.second==0.) ex2.second=0.0001; + if (inRange(sqrtS()/GeV, x-ex2.first, x+ex2.second)) { + mult->addPoint(x, sigma, ex, make_pair(error,error)); + } + else { + mult->addPoint(x, 0., ex, make_pair(0.,.0)); + } + } + } + + //@} + + + /// @name Histograms + //@{ + CounterPtr _nproton; + //@} + + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(BESIII_2019_I1718337); + + +} diff --git a/analyses/pluginBES/BESIII_2019_I1718337.info b/analyses/pluginBES/BESIII_2019_I1718337.info new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2019_I1718337.info @@ -0,0 +1,29 @@ +Name: BESIII_2019_I1718337 +Year: 2019 +Summary: Cross section for $e^+e^-\to$ proton and antiproton between 2. and 3.8GeV using ISR +Experiment: BESIII +Collider: BEPC II +InspireID: 1718337 +Status: VALIDATED +Authors: + - Peter Richardson +References: + - arXiv:1902.00665 +RunInfo: e+ e- to hadrons +NeedCrossSection: yes +Beams: [e+, e-] +Description: + 'Measurement of the cross section for $e^+e^-\to p\bar{p}$ for energies between 2. and 3.8 GeV using initial-state radiation.' +Keywords: [] +BibKey: Ablikim:2019njl +BibTeX: '@article{Ablikim:2019njl, + author = "Ablikim, M. and others", + title = "{Study of $e^+ e^- \to p \bar p$ via initial state + radiation at BESIII}", + collaboration = "BESIII", + year = "2019", + eprint = "1902.00665", + archivePrefix = "arXiv", + primaryClass = "hep-ex", + SLACcitation = "%%CITATION = ARXIV:1902.00665;%%" +}' diff --git a/analyses/pluginBES/BESIII_2019_I1718337.plot b/analyses/pluginBES/BESIII_2019_I1718337.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2019_I1718337.plot @@ -0,0 +1,5 @@ +BEGIN PLOT /BESIII_2019_I1718337/d01-x01-y01 +Title=$\sigma(e^+e^-\to p\bar{p})$ +XLabel=$\sqrt{s}$/GeV +YLabel=$\sigma(e^+e^-\to p\bar{p})$/pb +END PLOT diff --git a/analyses/pluginBES/BESIII_2019_I1718337.yoda b/analyses/pluginBES/BESIII_2019_I1718337.yoda new file mode 100644 --- /dev/null +++ b/analyses/pluginBES/BESIII_2019_I1718337.yoda @@ -0,0 +1,38 @@ +BEGIN YODA_SCATTER2D_V2 /REF/BESIII_2019_I1718337/d01-x01-y01 +IsRef: 1 +Path: /REF/BESIII_2019_I1718337/d01-x01-y01 +Title: +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +2.01250 0.01250 0.01250 7.97000e+02 9.36002e+01 9.36002e+01 +2.03750 0.01250 0.01250 8.33000e+02 8.29277e+01 8.29277e+01 +2.06250 0.01250 0.01250 7.23000e+02 6.76757e+01 6.76757e+01 +2.08750 0.01250 0.01250 7.49000e+02 5.78014e+01 5.78014e+01 +2.11250 0.01250 0.01250 6.54000e+02 5.63028e+01 5.63028e+01 +2.13750 0.01250 0.01250 6.37000e+02 4.94065e+01 4.94065e+01 +2.16250 0.01250 0.01250 5.57000e+02 4.74342e+01 4.74342e+01 +2.18750 0.01250 0.01250 4.67000e+02 3.92046e+01 3.92046e+01 +2.21250 0.01250 0.01250 3.71000e+02 3.42053e+01 3.42053e+01 +2.23750 0.01250 0.01250 3.10000e+02 2.90689e+01 2.90689e+01 +2.26250 0.01250 0.01250 2.25000e+02 2.26274e+01 2.26274e+01 +2.28750 0.01250 0.01250 1.92000e+02 1.97990e+01 1.97990e+01 +2.32500 0.02500 0.02500 1.36100e+02 1.13146e+01 1.13146e+01 +2.37500 0.02500 0.02500 1.16300e+02 1.21037e+01 1.21037e+01 +2.42500 0.02500 0.02500 1.26100e+02 9.56713e+00 9.56713e+00 +2.47500 0.02500 0.02500 1.00100e+02 9.12853e+00 9.12853e+00 +2.52500 0.02500 0.02500 6.74000e+01 6.86222e+00 6.86222e+00 +2.57500 0.02500 0.02500 6.11000e+01 5.90339e+00 5.90339e+00 +2.62500 0.02500 0.02500 4.10000e+01 4.70106e+00 4.70106e+00 +2.67500 0.02500 0.02500 3.36000e+01 3.94081e+00 3.94081e+00 +2.72500 0.02500 0.02500 3.07000e+01 4.24264e+00 4.24264e+00 +2.77500 0.02500 0.02500 2.68000e+01 3.61248e+00 3.61248e+00 +2.82500 0.02500 0.02500 2.16000e+01 3.25269e+00 3.25269e+00 +2.87500 0.02500 0.02500 2.04000e+01 2.84253e+00 2.84253e+00 +2.92500 0.02500 0.02500 1.02000e+01 2.72029e+00 2.72029e+00 +2.97500 0.02500 0.02500 1.41000e+01 2.64008e+00 2.64008e+00 +3.10000 0.10000 0.10000 1.11000e+01 1.69706e+00 1.69706e+00 +3.30000 0.10000 0.10000 3.59000e+00 6.51153e-01 6.51153e-01 +3.50000 0.10000 0.10000 2.18000e+00 3.92046e-01 3.92046e-01 +3.70000 0.10000 0.10000 6.40000e-01 2.62488e-01 2.62488e-01 +END YODA_SCATTER2D_V2