diff --git a/analyses/pluginNovosibirsk/SND_2016_I1471515.cc b/analyses/pluginNovosibirsk/SND_2016_I1471515.cc --- a/analyses/pluginNovosibirsk/SND_2016_I1471515.cc +++ b/analyses/pluginNovosibirsk/SND_2016_I1471515.cc @@ -1,122 +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 SND_2016_I1471515 : public Analysis { public: /// Constructor DEFAULT_RIVET_ANALYSIS_CTOR(SND_2016_I1471515); /// @name Analysis methods //@{ /// Book histograms and initialise projections before the run void init() { // Initialise and register projections declare(FinalState(), "FS"); declare(UnstableParticles(), "UFS"); _nOmegaEtaPi= bookCounter("/TMP/nOmegaEtaPi"); } void findChildren(const Particle & p,map & nRes, int &ncount) { foreach(const Particle &child, p.children()) { if(child.children().empty()) { nRes[child.pdgId()]-=1; --ncount; } else findChildren(child,nRes,ncount); } } /// Perform the per-event analysis void analyze(const Event& event) { const FinalState& fs = apply(event, "FS"); map nCount; int ntotal(0); foreach (const Particle& p, fs.particles()) { nCount[p.pdgId()] += 1; ++ntotal; } const FinalState& ufs = apply(event, "UFS"); foreach (const Particle& p, ufs.particles()) { if(p.children().empty()) continue; if(p.pdgId()!=223) continue; map nRes = nCount; int ncount = ntotal; findChildren(p,nRes,ncount); foreach (const Particle& p2, ufs.particles()) { if(p2.pdgId()!=221) continue; map nResB = nRes; int ncountB = ncount; findChildren(p2,nResB,ncountB); if(ncountB!=1) continue; bool matched2 = true; for(auto const & val : nResB) { if(val.first==111) { if(val.second!=1) { matched2 = false; break; } } else if(val.second!=0) { matched2 = false; break; } } if(matched2) { _nOmegaEtaPi->fill(event.weight()); break; } } } } /// Normalise histograms etc., after the run void finalize() { double sigma = _nOmegaEtaPi->val(); double error = _nOmegaEtaPi->err(); sigma *= crossSection()/ sumOfWeights() /nanobarn; error *= crossSection()/ sumOfWeights() /nanobarn; - Scatter2D temphisto(refData(1, 1, 1)); - Scatter2DPtr mult = bookScatter2D(1, 1, 1); + Scatter2D temphisto(refData(1, 1, 6)); + Scatter2DPtr mult = bookScatter2D(1, 1, 6); for (size_t b = 0; b < temphisto.numPoints(); b++) { const double x = temphisto.point(b).x(); pair ex = temphisto.point(b).xErrs(); pair ex2 = ex; if(ex2.first ==0.) ex2. first=0.0001; if(ex2.second==0.) ex2.second=0.0001; if (inRange(sqrtS()/MeV, x-ex2.first, x+ex2.second)) { mult->addPoint(x, sigma, ex, make_pair(error,error)); } else { mult->addPoint(x, 0., ex, make_pair(0.,.0)); } } } //@} /// @name Histograms //@{ CounterPtr _nOmegaEtaPi; //@} }; // The hook for the plugin system DECLARE_RIVET_PLUGIN(SND_2016_I1471515); } diff --git a/analyses/pluginNovosibirsk/SND_2016_I1471515.yoda b/analyses/pluginNovosibirsk/SND_2016_I1471515.yoda --- a/analyses/pluginNovosibirsk/SND_2016_I1471515.yoda +++ b/analyses/pluginNovosibirsk/SND_2016_I1471515.yoda @@ -1,18 +1,136 @@ BEGIN YODA_SCATTER2D_V2 /REF/SND_2016_I1471515/d01-x01-y01 +Variations: [""] +ErrorBreakdown: {0: {}, 1: {}, 2: {}, 3: {}, 4: {}, 5: {}, 6: {}, 7: {}, 8: {}, 9: {}, 10: {}, 11: {}} + IsRef: 1 Path: /REF/SND_2016_I1471515/d01-x01-y01 -Title: +Title: doi:10.17182/hepdata.82577.v1/t1 Type: Scatter2D --- -# xval xerr- xerr+ yval yerr- yerr+ -1725.00000 25.00000 25.00000 4.00000e-01 9.00000e-01 3.00000e-01 -1775.00000 25.00000 25.00000 1.10000e+00 1.10000e+00 6.00000e-01 -1812.50000 12.50000 12.50000 2.10000e+00 1.30000e+00 8.00000e-01 -1837.50000 12.50000 12.50000 2.80000e+00 1.50000e+00 1.00000e+00 -1862.50000 12.50000 12.50000 2.20000e+00 1.10000e+00 8.00000e-01 -1887.50000 12.50000 12.50000 3.00000e+00 2.90000e+00 1.60000e+00 -1912.50000 12.50000 12.50000 1.90000e+00 9.00000e-01 6.00000e-01 -1937.50000 12.50000 12.50000 2.80000e+00 1.10000e+00 8.00000e-01 -1962.50000 12.50000 12.50000 1.50000e+00 1.00000e+00 7.00000e-01 -1987.50000 12.50000 12.50000 1.90000e+00 9.00000e-01 6.00000e-01 +# xval xerr- xerr+ yval yerr- yerr+ +1.550000e+03 5.000000e+01 5.000000e+01 5.888000e+03 0.000000e+00 0.000000e+00 +1.650000e+03 5.000000e+01 5.000000e+01 5.004000e+03 0.000000e+00 0.000000e+00 +1.725000e+03 2.500000e+01 2.500000e+01 2.261000e+03 0.000000e+00 0.000000e+00 +1.775000e+03 2.500000e+01 2.500000e+01 2.392000e+03 0.000000e+00 0.000000e+00 +1.812500e+03 1.250000e+01 1.250000e+01 2.373000e+03 0.000000e+00 0.000000e+00 +1.837500e+03 1.250000e+01 1.250000e+01 1.897000e+03 0.000000e+00 0.000000e+00 +1.862500e+03 1.250000e+01 1.250000e+01 2.527000e+03 0.000000e+00 0.000000e+00 +1.887500e+03 1.250000e+01 1.250000e+01 6.620000e+02 0.000000e+00 0.000000e+00 +1.912500e+03 1.250000e+01 1.250000e+01 3.459000e+03 0.000000e+00 0.000000e+00 +1.937500e+03 1.250000e+01 1.250000e+01 2.361000e+03 0.000000e+00 0.000000e+00 +1.962500e+03 1.250000e+01 1.250000e+01 2.077000e+03 0.000000e+00 0.000000e+00 +1.987500e+03 1.250000e+01 1.250000e+01 2.682000e+03 0.000000e+00 0.000000e+00 END YODA_SCATTER2D_V2 +BEGIN YODA_SCATTER2D_V2 /REF/SND_2016_I1471515/d01-x01-y02 +Variations: [""] +ErrorBreakdown: {0: {}, 1: {}, 2: {}, 3: {}, 4: {}, 5: {}, 6: {}, 7: {}, 8: {}, 9: {}, 10: {}, 11: {}} + +IsRef: 1 +Path: /REF/SND_2016_I1471515/d01-x01-y02 +Title: doi:10.17182/hepdata.82577.v1/t1 +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.550000e+03 5.000000e+01 5.000000e+01 0.000000e+00 0.000000e+00 0.000000e+00 +1.650000e+03 5.000000e+01 5.000000e+01 0.000000e+00 0.000000e+00 0.000000e+00 +1.725000e+03 2.500000e+01 2.500000e+01 1.000000e+00 0.000000e+00 0.000000e+00 +1.775000e+03 2.500000e+01 2.500000e+01 3.000000e+00 0.000000e+00 0.000000e+00 +1.812500e+03 1.250000e+01 1.250000e+01 6.000000e+00 0.000000e+00 0.000000e+00 +1.837500e+03 1.250000e+01 1.250000e+01 7.000000e+00 0.000000e+00 0.000000e+00 +1.862500e+03 1.250000e+01 1.250000e+01 8.000000e+00 0.000000e+00 0.000000e+00 +1.887500e+03 1.250000e+01 1.250000e+01 3.000000e+00 0.000000e+00 0.000000e+00 +1.912500e+03 1.250000e+01 1.250000e+01 9.000000e+00 0.000000e+00 0.000000e+00 +1.937500e+03 1.250000e+01 1.250000e+01 1.100000e+01 0.000000e+00 0.000000e+00 +1.962500e+03 1.250000e+01 1.250000e+01 5.000000e+00 0.000000e+00 0.000000e+00 +1.987500e+03 1.250000e+01 1.250000e+01 9.000000e+00 0.000000e+00 0.000000e+00 +END YODA_SCATTER2D_V2 +BEGIN YODA_SCATTER2D_V2 /REF/SND_2016_I1471515/d01-x01-y03 +Variations: [""] +ErrorBreakdown: {0: {}, 1: {}, 2: {}, 3: {}, 4: {}, 5: {}, 6: {}, 7: {}, 8: {}, 9: {}, 10: {}, 11: {}} + +IsRef: 1 +Path: /REF/SND_2016_I1471515/d01-x01-y03 +Title: doi:10.17182/hepdata.82577.v1/t1 +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.550000e+03 5.000000e+01 5.000000e+01 2.000000e-02 0.000000e+00 0.000000e+00 +1.650000e+03 5.000000e+01 5.000000e+01 2.100000e-01 0.000000e+00 0.000000e+00 +1.725000e+03 2.500000e+01 2.500000e+01 1.300000e-01 0.000000e+00 0.000000e+00 +1.775000e+03 2.500000e+01 2.500000e+01 9.000000e-02 0.000000e+00 0.000000e+00 +1.812500e+03 1.250000e+01 1.250000e+01 9.000000e-02 0.000000e+00 0.000000e+00 +1.837500e+03 1.250000e+01 1.250000e+01 5.000000e-02 0.000000e+00 0.000000e+00 +1.862500e+03 1.250000e+01 1.250000e+01 6.000000e-02 0.000000e+00 0.000000e+00 +1.887500e+03 1.250000e+01 1.250000e+01 2.000000e-02 0.000000e+00 0.000000e+00 +1.912500e+03 1.250000e+01 1.250000e+01 7.000000e-02 0.000000e+00 0.000000e+00 +1.937500e+03 1.250000e+01 1.250000e+01 5.000000e-02 0.000000e+00 0.000000e+00 +1.962500e+03 1.250000e+01 1.250000e+01 6.000000e-02 0.000000e+00 0.000000e+00 +1.987500e+03 1.250000e+01 1.250000e+01 6.000000e-02 0.000000e+00 0.000000e+00 +END YODA_SCATTER2D_V2 +BEGIN YODA_SCATTER2D_V2 /REF/SND_2016_I1471515/d01-x01-y04 +Variations: [""] +ErrorBreakdown: {0: {}, 1: {}, 2: {}, 3: {}, 4: {}, 5: {}, 6: {}, 7: {}, 8: {}, 9: {}, 10: {}, 11: {}} + +IsRef: 1 +Path: /REF/SND_2016_I1471515/d01-x01-y04 +Title: doi:10.17182/hepdata.82577.v1/t1 +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.550000e+03 5.000000e+01 5.000000e+01 2.850000e+00 0.000000e+00 0.000000e+00 +1.650000e+03 5.000000e+01 5.000000e+01 4.030000e+00 0.000000e+00 0.000000e+00 +1.725000e+03 2.500000e+01 2.500000e+01 4.250000e+00 0.000000e+00 0.000000e+00 +1.775000e+03 2.500000e+01 2.500000e+01 4.290000e+00 0.000000e+00 0.000000e+00 +1.812500e+03 1.250000e+01 1.250000e+01 4.610000e+00 0.000000e+00 0.000000e+00 +1.837500e+03 1.250000e+01 1.250000e+01 5.020000e+00 0.000000e+00 0.000000e+00 +1.862500e+03 1.250000e+01 1.250000e+01 5.340000e+00 0.000000e+00 0.000000e+00 +1.887500e+03 1.250000e+01 1.250000e+01 5.580000e+00 0.000000e+00 0.000000e+00 +1.912500e+03 1.250000e+01 1.250000e+01 5.120000e+00 0.000000e+00 0.000000e+00 +1.937500e+03 1.250000e+01 1.250000e+01 5.980000e+00 0.000000e+00 0.000000e+00 +1.962500e+03 1.250000e+01 1.250000e+01 5.640000e+00 0.000000e+00 0.000000e+00 +1.987500e+03 1.250000e+01 1.250000e+01 6.140000e+00 0.000000e+00 0.000000e+00 +END YODA_SCATTER2D_V2 +BEGIN YODA_SCATTER2D_V2 /REF/SND_2016_I1471515/d01-x01-y05 +Variations: [""] +ErrorBreakdown: {0: {}, 1: {}, 2: {}, 3: {}, 4: {}, 5: {}, 6: {}, 7: {}, 8: {}, 9: {}, 10: {}, 11: {}} + +IsRef: 1 +Path: /REF/SND_2016_I1471515/d01-x01-y05 +Title: doi:10.17182/hepdata.82577.v1/t1 +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.550000e+03 5.000000e+01 5.000000e+01 7.790000e-01 0.000000e+00 0.000000e+00 +1.650000e+03 5.000000e+01 5.000000e+01 8.150000e-01 0.000000e+00 0.000000e+00 +1.725000e+03 2.500000e+01 2.500000e+01 8.130000e-01 0.000000e+00 0.000000e+00 +1.775000e+03 2.500000e+01 2.500000e+01 8.100000e-01 0.000000e+00 0.000000e+00 +1.812500e+03 1.250000e+01 1.250000e+01 8.230000e-01 0.000000e+00 0.000000e+00 +1.837500e+03 1.250000e+01 1.250000e+01 8.410000e-01 0.000000e+00 0.000000e+00 +1.862500e+03 1.250000e+01 1.250000e+01 8.550000e-01 0.000000e+00 0.000000e+00 +1.887500e+03 1.250000e+01 1.250000e+01 8.650000e-01 0.000000e+00 0.000000e+00 +1.912500e+03 1.250000e+01 1.250000e+01 8.710000e-01 0.000000e+00 0.000000e+00 +1.937500e+03 1.250000e+01 1.250000e+01 8.800000e-01 0.000000e+00 0.000000e+00 +1.962500e+03 1.250000e+01 1.250000e+01 8.870000e-01 0.000000e+00 0.000000e+00 +1.987500e+03 1.250000e+01 1.250000e+01 8.930000e-01 0.000000e+00 0.000000e+00 +END YODA_SCATTER2D_V2 +BEGIN YODA_SCATTER2D_V2 /REF/SND_2016_I1471515/d01-x01-y06 +Variations: [""] +ErrorBreakdown: {2: {error: {dn: -0.3, up: 0.9}}, 3: {error: {dn: -0.6, up: 1.1}}, 4: {error: {dn: -0.8, up: 1.3}}, 5: {error: {dn: -1.0, up: 1.5}}, 6: {error: {dn: -0.8, up: 1.1}}, 7: {error: {dn: -1.6, up: 2.9}}, 8: {error: {dn: -0.6, up: 0.9}}, 9: {error: {dn: -0.8, up: 1.1}}, 10: {error: {dn: -0.7, up: 1.0}}, 11: {error: {dn: -0.6, up: 0.9}}} + +IsRef: 1 +Path: /REF/SND_2016_I1471515/d01-x01-y06 +Title: doi:10.17182/hepdata.82577.v1/t1 +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +1.725000e+03 2.500000e+01 2.500000e+01 4.000000e-01 3.000000e-01 9.000000e-01 +1.775000e+03 2.500000e+01 2.500000e+01 1.100000e+00 6.000000e-01 1.100000e+00 +1.812500e+03 1.250000e+01 1.250000e+01 2.100000e+00 8.000000e-01 1.300000e+00 +1.837500e+03 1.250000e+01 1.250000e+01 2.800000e+00 1.000000e+00 1.500000e+00 +1.862500e+03 1.250000e+01 1.250000e+01 2.200000e+00 8.000000e-01 1.100000e+00 +1.887500e+03 1.250000e+01 1.250000e+01 3.000000e+00 1.600000e+00 2.900000e+00 +1.912500e+03 1.250000e+01 1.250000e+01 1.900000e+00 6.000000e-01 9.000000e-01 +1.937500e+03 1.250000e+01 1.250000e+01 2.800000e+00 8.000000e-01 1.100000e+00 +1.962500e+03 1.250000e+01 1.250000e+01 1.500000e+00 7.000000e-01 1.000000e+00 +1.987500e+03 1.250000e+01 1.250000e+01 1.900000e+00 6.000000e-01 9.000000e-01 +END YODA_SCATTER2D_V2