diff --git a/analyses/pluginRHIC/BRAHMS_2004_I647076.cc b/analyses/pluginRHIC/BRAHMS_2004_I647076.cc --- a/analyses/pluginRHIC/BRAHMS_2004_I647076.cc +++ b/analyses/pluginRHIC/BRAHMS_2004_I647076.cc @@ -1,164 +1,217 @@ // -*- C++ -*- #include "Rivet/Analysis.hh" #include "Rivet/Projections/SingleValueProjection.hh" +#include "Rivet/Projections/ImpactParameterProjection.hh" #include "Rivet/Projections/FinalState.hh" +#include "Rivet/Projections/UnstableFinalState.hh" #include "Rivet/Projections/ChargedFinalState.hh" namespace Rivet { /// @brief BRAHMS Centrality projection. class BRAHMSCentrality : public SingleValueProjection { public: // Constructor BRAHMSCentrality() : SingleValueProjection() { - setName("BRAHMSCentrality"); // Using here the BRAHMS reaction centrality from eg. 1602.01183, which // might not be correct. - const ChargedFinalState cfs(Cuts::pT > 0.1*GeV && Cuts::abseta < 2.2); - this->declare(cfs,"ChargedFinalState"); + declare(ChargedFinalState(Cuts::pT > 0.1*GeV && Cuts::abseta < 2.2), + "ChargedFinalState"); } // Destructor virtual ~BRAHMSCentrality() {} + // Clone on the heap. + DEFAULT_RIVET_PROJ_CLONE(BRAHMSCentrality); + + protected: // Do the projection. Count the number of charged particles in // the specified range. virtual void project(const Event& e) { clear(); set(apply (e, "ChargedFinalState").particles().size()); } // Compare to another projection. virtual int compare(const Projection& p) const { // This projection is only used for the analysis below. return UNDEFINED; } - // Clone this projection - virtual unique_ptr clone() const { - return unique_ptr(new BRAHMSCentrality(*this)); - } }; - /// @brief Brahms centrality analysis + /// @brief Brahms centrality calibration analysis based on the + // BrahmsCentrality projection. No data is given for this + // analysis, so one MUST do a calibration run. class BRAHMS_2004_CENTRALITY : public Analysis { public: // Constructor BRAHMS_2004_CENTRALITY() : Analysis("BRAHMS_2004_CENTRALITY") {} // Initialize the analysis void init() { - declare(BRAHMSCentrality(),"BCEN"); + declare(BRAHMSCentrality(),"Centrality"); + declare(ImpactParameterProjection(), "IMP"); + + // The central multiplicity. mult = bookHisto1D("mult",450,0,4500); - imp = bookHisto1D("imp",100,0,20); + + // Safeguard against filling preloaded histograms. + done = (mult->numEntries() > 0); + + // The impact parameter. + imp = bookHisto1D("mult_IMP",100,0,20); } // Analyse a single event void analyze(const Event& event) { - // Get and fill in the impact parameter value if the information - // is valid. - const HepMC::GenEvent* ge = event.genEvent(); - const HepMC::HeavyIon* hi = ge->heavy_ion(); - if (hi && hi->is_valid()) - imp->fill(hi->impact_parameter(), event.weight()); - - mult->fill(apply(event,"BCEN")(), event.weight()); + if (done) return; + // Fill impact parameter. + imp->fill(apply(event,"IMP")(), event.weight()); + // Fill multiplicity. + mult->fill(apply(event,"Centrality")(), event.weight()); } // Finalize the analysis void finalize() { - if(mult->numEntries()) mult->normalize(); - if(imp->numEntries()) imp->normalize(); + // Normalize the distributions, safeguarding against + // yoda normalization error. + if(mult->numEntries() > 0) mult->normalize(); + if(imp->numEntries() > 0) imp->normalize(); } private: // Histograms. Histo1DPtr mult; Histo1DPtr imp; + // Flag to test if we have preloaded histograms. + bool done; }; // The hook for the plugin system DECLARE_RIVET_PLUGIN(BRAHMS_2004_CENTRALITY); - /// @brief Brahms pT spectra for id particles AuAu @ 200GeV/nn + /// @brief Brahms pT spectra for id particles (pi+, pi-, K+, K-) + // in small bins of rapidity, 5% central collisions. + // System: AuAu @ 200GeV/nn. class BRAHMS_2004_I647076 : public Analysis { public: /// Constructor DEFAULT_RIVET_ANALYSIS_CTOR(BRAHMS_2004_I647076); /// @name Analysis methods //@{ /// Book histograms and initialise projections before the run void init() { // Initialise and register projections - // Centrality Projection - declareCentrality(BRAHMSCentrality(), "BRAHMS_2004_Centrality","mult","mult"); - // TODO: Change to PrimaryParticles depending on feed down. - declare(FinalState(Cuts::abseta < 5 && Cuts::pT > 100*MeV), "FS"); - - rapIntervals = {{-0.1,0.},{0.,0.1}}; + // Centrality Projection. + declareCentrality(BRAHMSCentrality(), "BRAHMS_2004_CENTRALITY","mult","BCEN"); + // TODO: Feed down correction is unclear. + declare(UnstableFinalState(Cuts::abseta < 5 && Cuts::pT > 100*MeV), "FS"); + // The measured rapidity intervals for pions. + rapIntervalsPi = {{-0.1,0.},{0.,0.1},{0.4,0.6},{0.6,0.8},{0.8,1.0}, + {1.0,1.2},{1.2,1.4},{2.1,2.3},{2.4,2.6},{3.0,3.1},{3.1,3.2},{3.2,3.3}, + {3.3,3.4},{3.4,3.66}}; + // The measured rapidity intervals for kaons. + rapIntervalsK = {{-0.1,0.},{0.,0.1},{0.4,0.6},{0.6,0.8},{0.8,1.0}, + {1.0,1.2},{2.0,2.2},{2.3,2.5},{2.9,3.0},{3.0,3.1},{3.1,3.2},{3.2,3.4}}; // Book histograms - for (int i = 1, N = rapIntervals.size(); i <= N; ++i) { + for (int i = 1, N = rapIntervalsPi.size(); i <= N; ++i) { piPlus.push_back(bookHisto1D(1, 1, i)); + piMinus.push_back(bookHisto1D(1, 1, 14 + i)); } + for (int i = 1, N = rapIntervalsK.size(); i <= N; ++i) { + kPlus.push_back(bookHisto1D(2, 1, i)); + kMinus.push_back(bookHisto1D(2, 1, 12 + i)); + } + // Counter for accepted sum of weights (centrality cut). centSow = bookCounter("centSow"); } /// Perform the per-event analysis void analyze(const Event& event) { const double w = event.weight(); - const CentralityProjection& cent = apply(event,"mult"); - if (cent() > 5.0) return; + // Reject all non-central events. The paper does not speak of + // any other event trigger, which in any case should matter + // little for central events. + if(apply(event,"BCEN")() > 5.0) return; + // Keep track of sum of weights. centSow->fill(w); - const FinalState& fs = apply(event,"FS"); + const UnstableFinalState& fs = apply(event,"FS"); + // Loop over particles. for (const auto& p : fs.particles()) { - double y = p.rapidity(); - for (int i = 0, N = rapIntervals.size(); i < N; ++i) { - if (y > rapIntervals[i].first && y <= rapIntervals[i].second) { - double dy = rapIntervals[i].second - rapIntervals[i].first; - double pT = p.pT(); - double nWeight = w / ( 2*M_PI*pT*dy); - if (p.pid() == 211) piPlus[i]->fill(pT, nWeight); - } - - } + const double y = p.rapidity(); + const double pT = p.pT(); + const int id = p.pid(); + // First pions. + if (abs(id) == 211) { + // Protect against decaying K0S and Lambda + if (p.hasAncestor(310) || p.hasAncestor(-310) || + p.hasAncestor(3122) || p.hasAncestor(3122)) continue; + for (int i = 0, N = rapIntervalsPi.size(); i < N; ++i) { + if (y > rapIntervalsPi[i].first && y <= rapIntervalsPi[i].second) { + const double dy = rapIntervalsPi[i].second - rapIntervalsPi[i].first; + const double nWeight = w / ( 2*M_PI*pT*dy); + if (id == 211) piPlus[i]->fill(pT, nWeight); + else piMinus[i]->fill(pT, nWeight); + break; + } + } + } + // Then kaons. + else if (abs(id) == 321) { + for (int i = 0, N = rapIntervalsK.size(); i < N; ++i) { + if (y > rapIntervalsK[i].first && y <= rapIntervalsK[i].second) { + const double dy = rapIntervalsK[i].second - rapIntervalsK[i].first; + const double nWeight = w / ( 2*M_PI*pT*dy); + if (id == 321) kPlus[i]->fill(pT, nWeight); + else kMinus[i]->fill(pT, nWeight); + break; + } + } + } } - } /// Normalise histograms etc., after the run void finalize() { - for (int i = 0, N = rapIntervals.size(); i < N; ++i) { + // Normalize all histograms to per-event yields. + for (int i = 0, N = rapIntervalsPi.size(); i < N; ++i) { piPlus[i]->scaleW(1./centSow->sumW()); + piMinus[i]->scaleW(1./centSow->sumW()); + } + for (int i = 0, N = rapIntervalsK.size(); i < N; ++i) { + kPlus[i]->scaleW(1./centSow->sumW()); + kMinus[i]->scaleW(1./centSow->sumW()); } } //@} // The rapidity intervals. - vector > rapIntervals; + vector > rapIntervalsPi; + vector > rapIntervalsK; /// @name Histograms //@{ vector piPlus; + vector piMinus; + vector kPlus; + vector kMinus; CounterPtr centSow; //@} - - }; - // The hook for the plugin system DECLARE_RIVET_PLUGIN(BRAHMS_2004_I647076); - - } diff --git a/analyses/pluginRHIC/BRAHMS_2004_I647076.info b/analyses/pluginRHIC/BRAHMS_2004_I647076.info --- a/analyses/pluginRHIC/BRAHMS_2004_I647076.info +++ b/analyses/pluginRHIC/BRAHMS_2004_I647076.info @@ -1,45 +1,45 @@ Name: BRAHMS_2004_I647076 Year: 2004 -Summary: Charged meson spectra for various rapiditues in central AuAu collisions +Summary: Charged meson pT spectra for various rapidities in central Au--Au collisions at 200 GeV. Experiment: BRAHMS Collider: RHIC InspireID: 647076 Status: UNVALIDATED Authors: - - Christian Bierlich References: - 'Phys.Rev.Lett. 94 (2005) 162301' - '10.1103/PhysRevLett.94.162301' - 'arXiv:043050' RunInfo: Au Au minimum bias events, only 5% most central events used in analysis. NeedCrossSection: no Options: - cent=REF,GEN,IMP,USR Beams: [1000791970, 1000791970] # This is _total_ energy of beams, so this becomes 197*200=39400 Energies: [39400] Description: - 'Measurements of invariant pT spectra for pions and kaons in AuAu collisions, over a - broad rapidity range $-0.1 < y < 3.5$. The paper quotes several interesting derived - quantities, like total dN/dy and , which can be fitted from the spectra produced in - this analysis' + 'Measurements of invariant $p_\perp$ spectra for pions and kaons in the 5% most central Au--Au collisions + at $\sqrt{s_{NN}} = 200$ GeV, over a broad rapidity range $-0.1 < y < 3.5$. The paper quotes several derived + quantities, like total $dN/dy$ and $\langle p_\perp \rangle$, which can be fitted from the spectra produced in + this analysis.' Keywords: [] BibKey: Bearden:2004yx BibTeX: '@article{Bearden:2004yx, author = "Bearden, I. G. and others", title = "{Charged meson rapidity distributions in central Au+Au collisions at s(NN)**(1/2) = 200-GeV}", collaboration = "BRAHMS", journal = "Phys. Rev. Lett.", volume = "94", year = "2005", pages = "162301", doi = "10.1103/PhysRevLett.94.162301", eprint = "nucl-ex/0403050", archivePrefix = "arXiv", primaryClass = "nucl-ex", SLACcitation = "%%CITATION = NUCL-EX/0403050;%%" }' ToDo: - Implement the analysis, test it, remove this ToDo, and mark as VALIDATED :-) diff --git a/analyses/pluginRHIC/BRAHMS_2004_I647076.plot b/analyses/pluginRHIC/BRAHMS_2004_I647076.plot --- a/analyses/pluginRHIC/BRAHMS_2004_I647076.plot +++ b/analyses/pluginRHIC/BRAHMS_2004_I647076.plot @@ -0,0 +1,260 @@ +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y01 +Title=Invariant $p_\perp$ of $\pi^+$, $-0.1 < y \leq 0$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y02 +Title=Invariant $p_\perp$ of $\pi^+$, $0 < y \leq 0.1$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y03 +Title=Invariant $p_\perp$ of $\pi^+$, $0.4 < y \leq 0.6$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y04 +Title=Invariant $p_\perp$ of $\pi^+$, $0.6 < y \leq 0.8$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y05 +Title=Invariant $p_\perp$ of $\pi^+$, $0.8 < y \leq 1.0$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y06 +Title=Invariant $p_\perp$ of $\pi^+$, $1.0 < y \leq 1.2$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y07 +Title=Invariant $p_\perp$ of $\pi^+$, $1.2 < y \leq 1.4$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y08 +Title=Invariant $p_\perp$ of $\pi^+$, $2.1 < y \leq 2.3$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y09 +Title=Invariant $p_\perp$ of $\pi^+$, $2.4 < y \leq 2.6$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y10 +Title=Invariant $p_\perp$ of $\pi^+$, $3.0 < y \leq 3.1$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y11 +Title=Invariant $p_\perp$ of $\pi^+$, $3.1 < y \leq 3.2$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y12 +Title=Invariant $p_\perp$ of $\pi^+$, $3.2 < y \leq 3.4$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y13 +Title=Invariant $p_\perp$ of $\pi^+$, $3.3 < y \leq 3.4$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y14 +Title=Invariant $p_\perp$ of $\pi^+$, $3.4 < y \leq 3.66$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y15 +Title=Invariant $p_\perp$ of $\pi^-$, $-0.1 < y \leq 0$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y16 +Title=Invariant $p_\perp$ of $\pi^-$, $0 < y \leq 0.1$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y17 +Title=Invariant $p_\perp$ of $\pi^-$, $0.4 < y \leq 0.6$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y18 +Title=Invariant $p_\perp$ of $\pi^-$, $0.6 < y \leq 0.8$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y19 +Title=Invariant $p_\perp$ of $\pi^-$, $0.8 < y \leq 1.0$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y20 +Title=Invariant $p_\perp$ of $\pi^-$, $1.0 < y \leq 1.2$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y21 +Title=Invariant $p_\perp$ of $\pi^-$, $1.2 < y \leq 1.4$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y22 +Title=Invariant $p_\perp$ of $\pi^-$, $2.1 < y \leq 2.3$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y23 +Title=Invariant $p_\perp$ of $\pi^-$, $2.4 < y \leq 2.6$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y24 +Title=Invariant $p_\perp$ of $\pi^-$, $3.0 < y \leq 3.1$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y25 +Title=Invariant $p_\perp$ of $\pi^-$, $3.1 < y \leq 3.2$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y26 +Title=Invariant $p_\perp$ of $\pi^-$, $3.2 < y \leq 3.4$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y27 +Title=Invariant $p_\perp$ of $\pi^-$, $3.3 < y \leq 3.4$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d01-x01-y28 +Title=Invariant $p_\perp$ of $\pi^-$, $3.4 < y \leq 3.66$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y01 +Title=Invariant $p_\perp$ of $K^+$, $-0.1 < y \leq 0$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y02 +Title=Invariant $p_\perp$ of $K^+$, $0 < y \leq 0.1$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y03 +Title=Invariant $p_\perp$ of $K^+$, $0.4 < y \leq 0.6$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y04 +Title=Invariant $p_\perp$ of $K^+$, $0.6 < y \leq 0.8$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y05 +Title=Invariant $p_\perp$ of $K^+$, $0.8 < y \leq 1.0$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y06 +Title=Invariant $p_\perp$ of $K^+$, $1.0 < y \leq 1.2$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y07 +Title=Invariant $p_\perp$ of $K^+$, $2.0 < y \leq 2.2$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y08 +Title=Invariant $p_\perp$ of $K^+$, $2.3 < y \leq 2.5$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y09 +Title=Invariant $p_\perp$ of $K^+$, $2.9 < y \leq 3.0$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y10 +Title=Invariant $p_\perp$ of $K^+$, $3.0 < y \leq 3.1$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y11 +Title=Invariant $p_\perp$ of $K^+$, $3.1 < y \leq 3.2$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y12 +Title=Invariant $p_\perp$ of $K^+$, $3.2 < y \leq 3.4$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y013 +Title=Invariant $p_\perp$ of $K^-$, $-0.1 < y \leq 0$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y14 +Title=Invariant $p_\perp$ of $K^-$, $0 < y \leq 0.1$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y15 +Title=Invariant $p_\perp$ of $K^-$, $0.4 < y \leq 0.6$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y16 +Title=Invariant $p_\perp$ of $K^-$, $0.6 < y \leq 0.8$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y17 +Title=Invariant $p_\perp$ of $K^-$, $0.8 < y \leq 1.0$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y18 +Title=Invariant $p_\perp$ of $K^-$, $1.0 < y \leq 1.2$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y19 +Title=Invariant $p_\perp$ of $K^-$, $2.0 < y \leq 2.2$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y20 +Title=Invariant $p_\perp$ of $K^-$, $2.3 < y \leq 2.5$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y21 +Title=Invariant $p_\perp$ of $K^-$, $2.9 < y \leq 3.0$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y22 +Title=Invariant $p_\perp$ of $K^-$, $3.0 < y \leq 3.1$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y23 +Title=Invariant $p_\perp$ of $K^-$, $3.1 < y \leq 3.2$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT +# BEGIN PLOT /BRAHMS_2004_I647076/d02-x01-y24 +Title=Invariant $p_\perp$ of $K^-$, $3.2 < y \leq 3.4$ +XLabel=$p_\perp$ [GeV] +YLabel=$\frac{1}{N}\frac{d^2N}{2\pi p_\perp dp_\perp dy}$ [GeV$^{-2}$] +# END PLOT diff --git a/analyses/pluginRHIC/BRAHMS_2004_I647076.yoda b/analyses/pluginRHIC/BRAHMS_2004_I647076.yoda --- a/analyses/pluginRHIC/BRAHMS_2004_I647076.yoda +++ b/analyses/pluginRHIC/BRAHMS_2004_I647076.yoda @@ -1,1136 +1,1226 @@ BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y01 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y01 -Title: ~ +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 5.000000e-02 5.000000e-02 3.198000e+02 1.740200e+00 1.740200e+00 3.500000e-01 5.000000e-02 5.000000e-02 2.034640e+02 5.246430e+00 5.246430e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.169420e+02 2.422670e+00 2.422670e+00 5.500000e-01 5.000000e-02 5.000000e-02 7.712500e+01 1.234830e+00 1.234830e+00 6.500000e-01 5.000000e-02 5.000000e-02 4.784130e+01 7.363690e-01 7.363690e-01 7.500000e-01 5.000000e-02 5.000000e-02 3.128740e+01 5.076160e-01 5.076160e-01 8.500000e-01 5.000000e-02 5.000000e-02 2.001260e+01 3.590320e-01 3.590320e-01 9.500000e-01 5.000000e-02 5.000000e-02 1.310690e+01 2.651790e-01 2.651790e-01 1.050000e+00 5.000000e-02 5.000000e-02 9.408240e+00 2.075910e-01 2.075910e-01 1.150000e+00 5.000000e-02 5.000000e-02 6.103240e+00 1.588070e-01 1.588070e-01 1.250000e+00 5.000000e-02 5.000000e-02 4.272470e+00 1.265950e-01 1.265950e-01 1.350000e+00 5.000000e-02 5.000000e-02 3.030110e+00 1.000630e-01 1.000630e-01 1.450000e+00 5.000000e-02 5.000000e-02 1.893260e+00 7.634400e-02 7.634400e-02 1.550000e+00 5.000000e-02 5.000000e-02 1.336390e+00 6.112500e-02 6.112500e-02 1.650000e+00 5.000000e-02 5.000000e-02 9.495600e-01 4.997700e-02 4.997700e-02 1.750000e+00 5.000000e-02 5.000000e-02 6.206000e-01 3.909400e-02 3.909400e-02 1.850000e+00 5.000000e-02 5.000000e-02 4.620100e-01 3.300100e-02 3.300100e-02 1.950000e+00 5.000000e-02 5.000000e-02 2.901900e-01 2.535400e-02 2.535400e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y02 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y02 Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 5.000000e-02 5.000000e-02 3.184450e+02 1.130830e+01 1.130830e+01 3.500000e-01 5.000000e-02 5.000000e-02 2.032250e+02 5.402560e+00 5.402560e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.223860e+02 2.870350e+00 2.870350e+00 5.500000e-01 5.000000e-02 5.000000e-02 7.289920e+01 1.705040e+00 1.705040e+00 6.500000e-01 5.000000e-02 5.000000e-02 4.914630e+01 1.148850e+00 1.148850e+00 7.500000e-01 5.000000e-02 5.000000e-02 3.088610e+01 7.645460e-01 7.645460e-01 8.500000e-01 5.000000e-02 5.000000e-02 2.046810e+01 5.320430e-01 5.320430e-01 9.500000e-01 5.000000e-02 5.000000e-02 1.290520e+01 3.726960e-01 3.726960e-01 1.050000e+00 5.000000e-02 5.000000e-02 9.260470e+00 2.908130e-01 2.908130e-01 1.150000e+00 5.000000e-02 5.000000e-02 6.070850e+00 2.219720e-01 2.219720e-01 1.250000e+00 5.000000e-02 5.000000e-02 4.046470e+00 1.646490e-01 1.646490e-01 1.350000e+00 5.000000e-02 5.000000e-02 2.724590e+00 1.282960e-01 1.282960e-01 1.450000e+00 5.000000e-02 5.000000e-02 1.891660e+00 1.022890e-01 1.022890e-01 1.550000e+00 5.000000e-02 5.000000e-02 1.313200e+00 8.021640e-02 8.021640e-02 1.650000e+00 5.000000e-02 5.000000e-02 9.134310e-01 6.458930e-02 6.458930e-02 1.750000e+00 5.000000e-02 5.000000e-02 7.018450e-01 5.463860e-02 5.463860e-02 1.850000e+00 5.000000e-02 5.000000e-02 5.209110e-01 4.466770e-02 4.466770e-02 1.950000e+00 5.000000e-02 5.000000e-02 2.741660e-01 3.144900e-02 3.144900e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y03 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y03 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 1.500000e-01 5.000000e-02 5.000000e-02 4.694910e+02 9.207470e+01 9.207470e+01 2.500000e-01 5.000000e-02 5.000000e-02 3.008980e+02 7.253180e+00 7.253180e+00 3.500000e-01 5.000000e-02 5.000000e-02 1.979610e+02 3.417690e+00 3.417690e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.197250e+02 2.092720e+00 2.092720e+00 5.500000e-01 5.000000e-02 5.000000e-02 7.794100e+01 1.442360e+00 1.442360e+00 6.500000e-01 5.000000e-02 5.000000e-02 4.648950e+01 9.993700e-01 9.993700e-01 7.500000e-01 5.000000e-02 5.000000e-02 2.988480e+01 7.370540e-01 7.370540e-01 8.500000e-01 5.000000e-02 5.000000e-02 1.981500e+01 5.540640e-01 5.540640e-01 9.500000e-01 5.000000e-02 5.000000e-02 1.300980e+01 4.227630e-01 4.227630e-01 1.050000e+00 5.000000e-02 5.000000e-02 8.489510e+00 3.243670e-01 3.243670e-01 1.150000e+00 5.000000e-02 5.000000e-02 6.116240e+00 2.605610e-01 2.605610e-01 1.250000e+00 5.000000e-02 5.000000e-02 3.890150e+00 1.987780e-01 1.987780e-01 1.350000e+00 5.000000e-02 5.000000e-02 2.621570e+00 1.563900e-01 1.563900e-01 1.450000e+00 5.000000e-02 5.000000e-02 1.759500e+00 1.234930e-01 1.234930e-01 1.550000e+00 5.000000e-02 5.000000e-02 1.201790e+00 9.812600e-02 9.812600e-02 1.650000e+00 5.000000e-02 5.000000e-02 8.262660e-01 7.914190e-02 7.914190e-02 1.750000e+00 5.000000e-02 5.000000e-02 5.794280e-01 9.794120e-02 9.794120e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y04 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y04 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 1.500000e-01 5.000000e-02 5.000000e-02 4.341380e+02 3.487080e+01 3.487080e+01 2.500000e-01 5.000000e-02 5.000000e-02 2.990630e+02 6.489110e+00 6.489110e+00 3.500000e-01 5.000000e-02 5.000000e-02 1.908930e+02 2.616450e+00 2.616450e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.192630e+02 1.755950e+00 1.755950e+00 5.500000e-01 5.000000e-02 5.000000e-02 7.551800e+01 1.240500e+00 1.240500e+00 6.500000e-01 5.000000e-02 5.000000e-02 4.588700e+01 8.756680e-01 8.756680e-01 7.500000e-01 5.000000e-02 5.000000e-02 2.951330e+01 6.477450e-01 6.477450e-01 8.500000e-01 5.000000e-02 5.000000e-02 1.920820e+01 4.875740e-01 4.875740e-01 9.500000e-01 5.000000e-02 5.000000e-02 1.247770e+01 3.669900e-01 3.669900e-01 1.050000e+00 5.000000e-02 5.000000e-02 8.061120e+00 2.796370e-01 2.796370e-01 1.150000e+00 5.000000e-02 5.000000e-02 5.692850e+00 2.246790e-01 2.246790e-01 1.250000e+00 5.000000e-02 5.000000e-02 3.701040e+00 1.731270e-01 1.731270e-01 1.350000e+00 5.000000e-02 5.000000e-02 2.599270e+00 1.391360e-01 1.391360e-01 1.450000e+00 5.000000e-02 5.000000e-02 1.638780e+00 1.062260e-01 1.062260e-01 1.550000e+00 5.000000e-02 5.000000e-02 1.174510e+00 1.146210e-01 1.146210e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y05 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y05 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 5.000000e-02 5.000000e-02 3.041250e+02 3.963740e+00 3.963740e+00 3.500000e-01 5.000000e-02 5.000000e-02 1.914750e+02 1.904210e+00 1.904210e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.183370e+02 1.085200e+00 1.085200e+00 5.500000e-01 5.000000e-02 5.000000e-02 7.271430e+01 7.075670e-01 7.075670e-01 6.500000e-01 5.000000e-02 5.000000e-02 4.580660e+01 4.995240e-01 4.995240e-01 7.500000e-01 5.000000e-02 5.000000e-02 2.914740e+01 3.634070e-01 3.634070e-01 8.500000e-01 5.000000e-02 5.000000e-02 1.888770e+01 2.711840e-01 2.711840e-01 9.500000e-01 5.000000e-02 5.000000e-02 1.246390e+01 2.057410e-01 2.057410e-01 1.050000e+00 5.000000e-02 5.000000e-02 7.959150e+00 1.553470e-01 1.553470e-01 1.150000e+00 5.000000e-02 5.000000e-02 5.218900e+00 1.195410e-01 1.195410e-01 1.250000e+00 5.000000e-02 5.000000e-02 3.263770e+00 9.031200e-02 9.031200e-02 1.350000e+00 5.000000e-02 5.000000e-02 2.255490e+00 8.860400e-02 8.860400e-02 1.450000e+00 5.000000e-02 5.000000e-02 1.742260e+00 1.528060e-01 1.528060e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y06 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y06 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 5.000000e-02 5.000000e-02 2.972910e+02 5.527230e+00 5.527230e+00 3.500000e-01 5.000000e-02 5.000000e-02 1.795720e+02 2.610730e+00 2.610730e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.143310e+02 1.603640e+00 1.603640e+00 5.500000e-01 5.000000e-02 5.000000e-02 7.115950e+01 1.061020e+00 1.061020e+00 6.500000e-01 5.000000e-02 5.000000e-02 4.318620e+01 7.254350e-01 7.254350e-01 7.500000e-01 5.000000e-02 5.000000e-02 2.722040e+01 5.230830e-01 5.230830e-01 8.500000e-01 5.000000e-02 5.000000e-02 1.786750e+01 3.900860e-01 3.900860e-01 9.500000e-01 5.000000e-02 5.000000e-02 1.206900e+01 2.992120e-01 2.992120e-01 1.050000e+00 5.000000e-02 5.000000e-02 7.144770e+00 2.176100e-01 2.176100e-01 1.150000e+00 5.000000e-02 5.000000e-02 4.275400e+00 2.054310e-01 2.054310e-01 1.250000e+00 5.000000e-02 5.000000e-02 2.976970e+00 1.686050e-01 1.686050e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y07 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y07 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 1.500000e-01 5.000000e-02 5.000000e-02 4.112930e+02 1.837520e+01 1.837520e+01 2.500000e-01 5.000000e-02 5.000000e-02 2.936310e+02 4.883680e+00 4.883680e+00 3.500000e-01 5.000000e-02 5.000000e-02 1.738920e+02 2.393560e+00 2.393560e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.063340e+02 1.494250e+00 1.494250e+00 5.500000e-01 5.000000e-02 5.000000e-02 6.700870e+01 1.028350e+00 1.028350e+00 6.500000e-01 5.000000e-02 5.000000e-02 4.196770e+01 7.329000e-01 7.329000e-01 7.500000e-01 5.000000e-02 5.000000e-02 2.621790e+01 5.337270e-01 5.337270e-01 8.500000e-01 5.000000e-02 5.000000e-02 1.515470e+01 3.802960e-01 3.802960e-01 9.500000e-01 5.000000e-02 5.000000e-02 9.978050e+00 3.061850e-01 3.061850e-01 1.050000e+00 5.000000e-02 5.000000e-02 6.633520e+00 3.792130e-01 3.792130e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y08 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y08 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 5.000000e-02 5.000000e-02 2.129720e+02 1.146600e+01 1.146600e+01 3.500000e-01 5.000000e-02 5.000000e-02 1.270710e+02 3.891910e+00 3.891910e+00 4.500000e-01 5.000000e-02 5.000000e-02 8.429280e+01 2.421250e+00 2.421250e+00 5.500000e-01 5.000000e-02 5.000000e-02 5.371060e+01 1.592870e+00 1.592870e+00 6.500000e-01 5.000000e-02 5.000000e-02 3.262240e+01 7.525800e-01 7.525800e-01 7.500000e-01 5.000000e-02 5.000000e-02 2.116870e+01 5.262600e-01 5.262600e-01 8.500000e-01 5.000000e-02 5.000000e-02 1.363870e+01 3.930600e-01 3.930600e-01 9.500000e-01 5.000000e-02 5.000000e-02 8.591067e+00 2.708600e-01 2.708600e-01 1.050000e+00 5.000000e-02 5.000000e-02 5.589277e+00 2.028700e-01 2.028700e-01 1.150000e+00 5.000000e-02 5.000000e-02 3.266547e+00 1.733700e-01 1.733700e-01 1.250000e+00 5.000000e-02 5.000000e-02 2.206237e+00 1.412300e-01 1.412300e-01 1.350000e+00 5.000000e-02 5.000000e-02 1.430700e+00 1.110400e-01 1.110400e-01 1.450000e+00 5.000000e-02 5.000000e-02 8.327570e-01 8.369500e-02 8.369500e-02 1.550000e+00 5.000000e-02 5.000000e-02 8.739630e-01 8.569900e-02 8.569900e-02 1.650000e+00 5.000000e-02 5.000000e-02 4.267110e-01 5.861300e-02 5.861300e-02 1.750000e+00 5.000000e-02 5.000000e-02 4.215370e-01 5.902700e-02 5.902700e-02 1.850000e+00 5.000000e-02 5.000000e-02 1.441420e-01 3.603500e-02 3.603500e-02 1.950000e+00 5.000000e-02 5.000000e-02 7.014290e-02 2.479900e-02 2.479900e-02 2.050000e+00 5.000000e-02 5.000000e-02 2.891980e-02 1.669600e-02 1.669600e-02 2.150000e+00 5.000000e-02 5.000000e-02 3.285590e-02 1.896900e-02 1.896900e-02 2.250000e+00 5.000000e-02 5.000000e-02 5.341070e-02 2.670500e-02 2.670500e-02 2.350000e+00 5.000000e-02 5.000000e-02 4.775460e-02 2.757100e-02 2.757100e-02 2.450000e+00 5.000000e-02 5.000000e-02 1.629680e-02 1.029600e-02 1.029600e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y09 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y09 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 4.500000e-01 5.000000e-02 5.000000e-02 6.796160e+01 1.862140e+00 1.862140e+00 5.500000e-01 5.000000e-02 5.000000e-02 4.324420e+01 1.114330e+00 1.114330e+00 6.500000e-01 5.000000e-02 5.000000e-02 2.884020e+01 6.590420e-01 6.590420e-01 7.500000e-01 5.000000e-02 5.000000e-02 1.763720e+01 4.484190e-01 4.484190e-01 8.500000e-01 5.000000e-02 5.000000e-02 9.167380e+00 3.908980e-01 3.908980e-01 9.500000e-01 5.000000e-02 5.000000e-02 5.750220e+00 4.066020e-01 4.066020e-01 1.050000e+00 5.000000e-02 5.000000e-02 4.066340e+00 3.388620e-01 3.388620e-01 1.150000e+00 5.000000e-02 5.000000e-02 2.636040e+00 2.704520e-01 2.704520e-01 1.250000e+00 5.000000e-02 5.000000e-02 1.578260e+00 2.147750e-01 2.147750e-01 1.350000e+00 5.000000e-02 5.000000e-02 1.185000e+00 1.975000e-01 1.975000e-01 1.450000e+00 5.000000e-02 5.000000e-02 7.104750e-01 1.629940e-01 1.629940e-01 1.550000e+00 5.000000e-02 5.000000e-02 4.567960e-01 1.377290e-01 1.377290e-01 1.650000e+00 5.000000e-02 5.000000e-02 2.411020e-01 1.078240e-01 1.078240e-01 1.750000e+00 5.000000e-02 5.000000e-02 1.753900e-01 1.012610e-01 1.012610e-01 1.850000e+00 5.000000e-02 5.000000e-02 7.465000e-02 3.465000e-02 3.465000e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y10 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y10 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 3.500000e-01 5.000000e-02 5.000000e-02 8.737660e+01 2.440350e+00 2.440350e+00 4.500000e-01 5.000000e-02 5.000000e-02 5.318970e+01 8.770450e-01 8.770450e-01 5.500000e-01 5.000000e-02 5.000000e-02 3.225710e+01 7.641370e-01 7.641370e-01 6.500000e-01 5.000000e-02 5.000000e-02 1.971550e+01 6.256540e-01 6.256540e-01 7.500000e-01 5.000000e-02 5.000000e-02 1.190100e+01 4.696960e-01 4.696960e-01 8.500000e-01 5.000000e-02 5.000000e-02 8.097480e+00 3.715380e-01 3.715380e-01 9.500000e-01 5.000000e-02 5.000000e-02 5.288130e+00 2.942390e-01 2.942390e-01 1.050000e+00 5.000000e-02 5.000000e-02 3.075290e+00 1.972800e-01 1.972800e-01 1.150000e+00 5.000000e-02 5.000000e-02 2.059260e+00 1.470900e-01 1.470900e-01 1.250000e+00 5.000000e-02 5.000000e-02 1.547590e+00 1.169870e-01 1.169870e-01 1.350000e+00 5.000000e-02 5.000000e-02 8.503300e-01 8.298400e-02 8.298400e-02 1.450000e+00 5.000000e-02 5.000000e-02 5.246490e-01 6.558100e-02 6.558100e-02 1.550000e+00 5.000000e-02 5.000000e-02 3.939510e-01 5.939000e-02 5.939000e-02 1.650000e+00 5.000000e-02 5.000000e-02 1.654710e-01 4.136800e-02 4.136800e-02 1.750000e+00 5.000000e-02 5.000000e-02 1.960040e-01 4.753800e-02 4.753800e-02 1.850000e+00 5.000000e-02 5.000000e-02 4.266100e-02 2.463000e-02 2.463000e-02 1.950000e+00 5.000000e-02 5.000000e-02 7.216200e-02 3.608200e-02 3.608200e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y11 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y11 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 4.500000e-01 5.000000e-02 5.000000e-02 4.990830e+01 9.458850e-01 9.458850e-01 5.500000e-01 5.000000e-02 5.000000e-02 2.840650e+01 9.641810e-01 9.641810e-01 6.500000e-01 5.000000e-02 5.000000e-02 2.051750e+01 7.049900e-01 7.049900e-01 7.500000e-01 5.000000e-02 5.000000e-02 1.207450e+01 4.908970e-01 4.908970e-01 8.500000e-01 5.000000e-02 5.000000e-02 6.997800e+00 3.324750e-01 3.324750e-01 9.500000e-01 5.000000e-02 5.000000e-02 4.849640e+00 2.424820e-01 2.424820e-01 1.050000e+00 5.000000e-02 5.000000e-02 3.298230e+00 1.910610e-01 1.910610e-01 1.150000e+00 5.000000e-02 5.000000e-02 2.015910e+00 1.486150e-01 1.486150e-01 1.250000e+00 5.000000e-02 5.000000e-02 1.260900e+00 1.321790e-01 1.321790e-01 1.350000e+00 5.000000e-02 5.000000e-02 6.567850e-01 1.110170e-01 1.110170e-01 1.450000e+00 5.000000e-02 5.000000e-02 6.462480e-01 1.377810e-01 1.377810e-01 1.550000e+00 5.000000e-02 5.000000e-02 2.944760e-01 1.472380e-01 1.472380e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y12 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y12 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 5.000000e-02 5.000000e-02 1.233910e+02 2.784280e+00 2.784280e+00 3.500000e-01 5.000000e-02 5.000000e-02 6.925240e+01 1.391460e+00 1.391460e+00 4.500000e-01 5.000000e-02 5.000000e-02 5.367010e+01 1.889280e+00 1.889280e+00 5.500000e-01 5.000000e-02 5.000000e-02 3.164830e+01 9.585980e-01 9.585980e-01 6.500000e-01 5.000000e-02 5.000000e-02 1.757810e+01 7.236800e-01 7.236800e-01 7.500000e-01 5.000000e-02 5.000000e-02 1.208510e+01 4.925520e-01 4.925520e-01 8.500000e-01 5.000000e-02 5.000000e-02 7.463730e+00 3.361500e-01 3.361500e-01 9.500000e-01 5.000000e-02 5.000000e-02 4.400600e+00 2.620520e-01 2.620520e-01 1.050000e+00 5.000000e-02 5.000000e-02 2.632150e+00 2.556570e-01 2.556570e-01 1.150000e+00 5.000000e-02 5.000000e-02 2.490760e+00 4.210160e-01 4.210160e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y13 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y13 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 5.000000e-02 5.000000e-02 1.074270e+02 7.693030e+00 7.693030e+00 3.500000e-01 5.000000e-02 5.000000e-02 5.918960e+01 3.706590e+00 3.706590e+00 4.500000e-01 5.000000e-02 5.000000e-02 4.433330e+01 1.261020e+00 1.261020e+00 5.500000e-01 5.000000e-02 5.000000e-02 2.328380e+01 6.485250e-01 6.485250e-01 6.500000e-01 5.000000e-02 5.000000e-02 1.528760e+01 4.235120e-01 4.235120e-01 7.500000e-01 5.000000e-02 5.000000e-02 9.882980e+00 2.906760e-01 2.906760e-01 8.500000e-01 5.000000e-02 5.000000e-02 5.497380e+00 2.016790e-01 2.016790e-01 9.500000e-01 5.000000e-02 5.000000e-02 3.354030e+00 1.597160e-01 1.597160e-01 1.050000e+00 5.000000e-02 5.000000e-02 2.234180e+00 1.541730e-01 1.541730e-01 1.150000e+00 5.000000e-02 5.000000e-02 1.448840e+00 1.407240e-01 1.407240e-01 1.250000e+00 5.000000e-02 5.000000e-02 8.366100e-01 1.220320e-01 1.220320e-01 1.350000e+00 5.000000e-02 5.000000e-02 5.039600e-01 1.156160e-01 1.156160e-01 1.450000e+00 5.000000e-02 5.000000e-02 2.860400e-01 1.279210e-01 1.279210e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y14 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y14 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 3.500000e-01 5.000000e-02 5.000000e-02 6.583180e+01 1.956650e+00 1.956650e+00 4.500000e-01 5.000000e-02 5.000000e-02 3.554520e+01 8.608330e-01 8.608330e-01 5.500000e-01 5.000000e-02 5.000000e-02 2.416650e+01 5.745780e-01 5.745780e-01 6.500000e-01 5.000000e-02 5.000000e-02 1.405870e+01 3.957450e-01 3.957450e-01 7.500000e-01 5.000000e-02 5.000000e-02 8.835690e+00 3.388330e-01 3.388330e-01 8.500000e-01 5.000000e-02 5.000000e-02 5.544910e+00 3.267370e-01 3.267370e-01 9.500000e-01 5.000000e-02 5.000000e-02 3.201930e+00 3.579860e-01 3.579860e-01 1.050000e+00 5.000000e-02 5.000000e-02 2.029600e+00 1.014800e-01 1.014800e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y15 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y15 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 5.000000e-02 5.000000e-02 2.976810e+02 1.121930e+01 1.121930e+01 3.500000e-01 5.000000e-02 5.000000e-02 2.033150e+02 5.270690e+00 5.270690e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.193320e+02 2.762490e+00 2.762490e+00 5.500000e-01 5.000000e-02 5.000000e-02 7.792150e+01 1.754700e+00 1.754700e+00 6.500000e-01 5.000000e-02 5.000000e-02 4.675840e+01 1.113290e+00 1.113290e+00 7.500000e-01 5.000000e-02 5.000000e-02 3.150530e+01 7.693380e-01 7.693380e-01 8.500000e-01 5.000000e-02 5.000000e-02 2.152140e+01 5.219710e-01 5.219710e-01 9.500000e-01 5.000000e-02 5.000000e-02 1.341770e+01 3.697300e-01 3.697300e-01 1.050000e+00 5.000000e-02 5.000000e-02 9.027010e+00 2.713130e-01 2.713130e-01 1.150000e+00 5.000000e-02 5.000000e-02 5.969500e+00 1.979960e-01 1.979960e-01 1.250000e+00 5.000000e-02 5.000000e-02 4.209080e+00 1.524790e-01 1.524790e-01 1.350000e+00 5.000000e-02 5.000000e-02 2.568980e+00 1.121200e-01 1.121200e-01 1.450000e+00 5.000000e-02 5.000000e-02 1.861030e+00 8.763240e-02 8.763240e-02 1.550000e+00 5.000000e-02 5.000000e-02 1.256980e+00 6.898600e-02 6.898600e-02 1.650000e+00 5.000000e-02 5.000000e-02 8.146510e-01 5.336960e-02 5.336960e-02 1.750000e+00 5.000000e-02 5.000000e-02 6.682830e-01 4.600650e-02 4.600650e-02 1.850000e+00 5.000000e-02 5.000000e-02 4.303710e-01 3.598950e-02 3.598950e-02 1.950000e+00 5.000000e-02 5.000000e-02 3.091530e-01 2.895480e-02 2.895480e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y16 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y16 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 5.000000e-02 5.000000e-02 3.272500e+02 1.160640e+01 1.160640e+01 3.500000e-01 5.000000e-02 5.000000e-02 1.946980e+02 5.300970e+00 5.300970e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.146480e+02 2.553410e+00 2.553410e+00 5.500000e-01 5.000000e-02 5.000000e-02 7.171600e+01 1.258180e+00 1.258180e+00 6.500000e-01 5.000000e-02 5.000000e-02 4.664460e+01 7.411380e-01 7.411380e-01 7.500000e-01 5.000000e-02 5.000000e-02 3.039480e+01 4.962130e-01 4.962130e-01 8.500000e-01 5.000000e-02 5.000000e-02 1.980710e+01 3.472790e-01 3.472790e-01 9.500000e-01 5.000000e-02 5.000000e-02 1.262820e+01 2.538890e-01 2.538890e-01 1.050000e+00 5.000000e-02 5.000000e-02 9.313920e+00 2.030530e-01 2.030530e-01 1.150000e+00 5.000000e-02 5.000000e-02 6.039680e+00 1.556330e-01 1.556330e-01 1.250000e+00 5.000000e-02 5.000000e-02 4.081380e+00 1.223930e-01 1.223930e-01 1.350000e+00 5.000000e-02 5.000000e-02 2.676000e+00 9.379360e-02 9.379360e-02 1.450000e+00 5.000000e-02 5.000000e-02 1.796050e+00 7.483530e-02 7.483530e-02 1.550000e+00 5.000000e-02 5.000000e-02 1.323520e+00 6.253030e-02 6.253030e-02 1.650000e+00 5.000000e-02 5.000000e-02 8.808570e-01 4.931850e-02 4.931850e-02 1.750000e+00 5.000000e-02 5.000000e-02 6.166410e-01 4.066010e-02 4.066010e-02 1.850000e+00 5.000000e-02 5.000000e-02 4.853180e-01 3.493400e-02 3.493400e-02 1.950000e+00 5.000000e-02 5.000000e-02 3.275290e-01 2.808540e-02 2.808540e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y17 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y17 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 1.500000e-01 5.000000e-02 5.000000e-02 4.564020e+02 8.625190e+01 8.625190e+01 2.500000e-01 5.000000e-02 5.000000e-02 2.964270e+02 7.306390e+00 7.306390e+00 3.500000e-01 5.000000e-02 5.000000e-02 2.000640e+02 3.501280e+00 3.501280e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.241590e+02 2.171550e+00 2.171550e+00 5.500000e-01 5.000000e-02 5.000000e-02 7.670470e+01 1.454790e+00 1.454790e+00 6.500000e-01 5.000000e-02 5.000000e-02 4.694460e+01 1.019570e+00 1.019570e+00 7.500000e-01 5.000000e-02 5.000000e-02 3.091340e+01 7.585120e-01 7.585120e-01 8.500000e-01 5.000000e-02 5.000000e-02 2.112470e+01 5.788140e-01 5.788140e-01 9.500000e-01 5.000000e-02 5.000000e-02 1.325200e+01 4.310880e-01 4.310880e-01 1.050000e+00 5.000000e-02 5.000000e-02 8.208010e+00 3.221930e-01 3.221930e-01 1.150000e+00 5.000000e-02 5.000000e-02 6.119380e+00 2.630930e-01 2.630930e-01 1.250000e+00 5.000000e-02 5.000000e-02 4.203710e+00 2.086260e-01 2.086260e-01 1.350000e+00 5.000000e-02 5.000000e-02 2.613730e+00 1.573280e-01 1.573280e-01 1.450000e+00 5.000000e-02 5.000000e-02 1.869980e+00 1.281290e-01 1.281290e-01 1.550000e+00 5.000000e-02 5.000000e-02 1.167970e+00 9.733100e-02 9.733100e-02 1.650000e+00 5.000000e-02 5.000000e-02 9.035700e-01 8.318100e-02 8.318100e-02 1.750000e+00 5.000000e-02 5.000000e-02 7.941900e-01 1.158450e-01 1.158450e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y18 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y18 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 1.500000e-01 5.000000e-02 5.000000e-02 3.401470e+02 2.624290e+01 2.624290e+01 2.500000e-01 5.000000e-02 5.000000e-02 3.008720e+02 5.827090e+00 5.827090e+00 3.500000e-01 5.000000e-02 5.000000e-02 1.937510e+02 2.439480e+00 2.439480e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.184350e+02 1.638470e+00 1.638470e+00 5.500000e-01 5.000000e-02 5.000000e-02 7.508010e+01 1.169850e+00 1.169850e+00 6.500000e-01 5.000000e-02 5.000000e-02 4.627210e+01 8.384060e-01 8.384060e-01 7.500000e-01 5.000000e-02 5.000000e-02 2.974240e+01 6.232970e-01 6.232970e-01 8.500000e-01 5.000000e-02 5.000000e-02 1.982100e+01 4.764050e-01 4.764050e-01 9.500000e-01 5.000000e-02 5.000000e-02 1.275630e+01 3.589410e-01 3.589410e-01 1.050000e+00 5.000000e-02 5.000000e-02 8.329670e+00 2.756720e-01 2.756720e-01 1.150000e+00 5.000000e-02 5.000000e-02 5.833590e+00 2.211220e-01 2.211220e-01 1.250000e+00 5.000000e-02 5.000000e-02 3.633700e+00 1.670780e-01 1.670780e-01 1.350000e+00 5.000000e-02 5.000000e-02 2.492940e+00 1.330630e-01 1.330630e-01 1.450000e+00 5.000000e-02 5.000000e-02 1.688880e+00 1.055550e-01 1.055550e-01 1.550000e+00 5.000000e-02 5.000000e-02 1.149800e+00 1.122090e-01 1.122090e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y19 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y19 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 5.000000e-02 5.000000e-02 3.138530e+02 4.092270e+00 4.092270e+00 3.500000e-01 5.000000e-02 5.000000e-02 1.921930e+02 1.939360e+00 1.939360e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.190280e+02 1.116270e+00 1.116270e+00 5.500000e-01 5.000000e-02 5.000000e-02 7.324270e+01 7.290450e-01 7.290450e-01 6.500000e-01 5.000000e-02 5.000000e-02 4.547210e+01 5.120230e-01 5.120230e-01 7.500000e-01 5.000000e-02 5.000000e-02 2.982340e+01 3.785440e-01 3.785440e-01 8.500000e-01 5.000000e-02 5.000000e-02 1.912040e+01 2.809700e-01 2.809700e-01 9.500000e-01 5.000000e-02 5.000000e-02 1.217020e+01 2.095200e-01 2.095200e-01 1.050000e+00 5.000000e-02 5.000000e-02 8.078060e+00 1.612070e-01 1.612070e-01 1.150000e+00 5.000000e-02 5.000000e-02 5.082490e+00 1.214950e-01 1.214950e-01 1.250000e+00 5.000000e-02 5.000000e-02 3.312310e+00 9.383700e-02 9.383700e-02 1.350000e+00 5.000000e-02 5.000000e-02 2.445530e+00 9.371300e-02 9.371300e-02 1.450000e+00 5.000000e-02 5.000000e-02 1.530490e+00 1.408930e-01 1.408930e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y20 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y20 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 5.000000e-02 5.000000e-02 2.937950e+02 5.274160e+00 5.274160e+00 3.500000e-01 5.000000e-02 5.000000e-02 1.804090e+02 2.210160e+00 2.210160e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.114330e+02 1.296440e+00 1.296440e+00 5.500000e-01 5.000000e-02 5.000000e-02 6.891490e+01 8.623800e-01 8.623800e-01 6.500000e-01 5.000000e-02 5.000000e-02 4.303390e+01 6.025950e-01 6.025950e-01 7.500000e-01 5.000000e-02 5.000000e-02 2.795390e+01 4.424890e-01 4.424890e-01 8.500000e-01 5.000000e-02 5.000000e-02 1.736510e+01 3.222400e-01 3.222400e-01 9.500000e-01 5.000000e-02 5.000000e-02 1.137730e+01 2.444610e-01 2.444610e-01 1.050000e+00 5.000000e-02 5.000000e-02 6.994810e+00 1.818830e-01 1.818830e-01 1.150000e+00 5.000000e-02 5.000000e-02 4.595630e+00 1.593250e-01 1.593250e-01 1.250000e+00 5.000000e-02 5.000000e-02 3.161950e+00 2.141540e-01 2.141540e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y21 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y21 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 1.500000e-01 5.000000e-02 5.000000e-02 3.958560e+02 1.962180e+01 1.962180e+01 2.500000e-01 5.000000e-02 5.000000e-02 2.866740e+02 6.100840e+00 6.100840e+00 3.500000e-01 5.000000e-02 5.000000e-02 1.822020e+02 3.178480e+00 3.178480e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.104590e+02 1.949310e+00 1.949310e+00 5.500000e-01 5.000000e-02 5.000000e-02 6.719420e+01 1.272580e+00 1.272580e+00 6.500000e-01 5.000000e-02 5.000000e-02 4.329200e+01 8.930450e-01 8.930450e-01 7.500000e-01 5.000000e-02 5.000000e-02 2.582070e+01 6.227740e-01 6.227740e-01 8.500000e-01 5.000000e-02 5.000000e-02 1.634950e+01 4.562700e-01 4.562700e-01 9.500000e-01 5.000000e-02 5.000000e-02 1.034970e+01 3.596780e-01 3.596780e-01 1.050000e+00 5.000000e-02 5.000000e-02 6.064570e+00 4.052060e-01 4.052060e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y22 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y22 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 5.000000e-02 5.000000e-02 2.055600e+02 7.655460e+00 7.655460e+00 3.500000e-01 5.000000e-02 5.000000e-02 1.369450e+02 2.784950e+00 2.784950e+00 4.500000e-01 5.000000e-02 5.000000e-02 8.591320e+01 1.686520e+00 1.686520e+00 5.500000e-01 5.000000e-02 5.000000e-02 5.286110e+01 1.091370e+00 1.091370e+00 6.500000e-01 5.000000e-02 5.000000e-02 3.167450e+01 5.203740e-01 5.203740e-01 7.500000e-01 5.000000e-02 5.000000e-02 1.901100e+01 3.394810e-01 3.394810e-01 8.500000e-01 5.000000e-02 5.000000e-02 1.166970e+01 2.790390e-01 2.790390e-01 9.500000e-01 5.000000e-02 5.000000e-02 7.901210e+00 2.092350e-01 2.092350e-01 1.050000e+00 5.000000e-02 5.000000e-02 5.541160e+00 1.661680e-01 1.661680e-01 1.150000e+00 5.000000e-02 5.000000e-02 3.084650e+00 1.419830e-01 1.419830e-01 1.250000e+00 5.000000e-02 5.000000e-02 1.958180e+00 1.106830e-01 1.106830e-01 1.350000e+00 5.000000e-02 5.000000e-02 1.327370e+00 8.908700e-02 8.908700e-02 1.450000e+00 5.000000e-02 5.000000e-02 9.436730e-01 8.340900e-02 8.340900e-02 1.550000e+00 5.000000e-02 5.000000e-02 8.183760e-01 9.449700e-02 9.449700e-02 1.650000e+00 5.000000e-02 5.000000e-02 3.225220e-01 5.888400e-02 5.888400e-02 1.750000e+00 5.000000e-02 5.000000e-02 2.302380e-01 5.024200e-02 5.024200e-02 1.850000e+00 5.000000e-02 5.000000e-02 2.380370e-01 5.194300e-02 5.194300e-02 1.950000e+00 5.000000e-02 5.000000e-02 1.520670e-01 4.217500e-02 4.217500e-02 2.050000e+00 5.000000e-02 5.000000e-02 5.279100e-02 2.639500e-02 2.639500e-02 2.150000e+00 5.000000e-02 5.000000e-02 5.323000e-02 2.661500e-02 2.661500e-02 2.250000e+00 5.000000e-02 5.000000e-02 6.762400e-02 3.024200e-02 3.024200e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y23 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y23 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 4.500000e-01 5.000000e-02 5.000000e-02 7.181910e+01 9.632560e-01 9.632560e-01 5.500000e-01 5.000000e-02 5.000000e-02 4.398850e+01 6.124900e-01 6.124900e-01 6.500000e-01 5.000000e-02 5.000000e-02 2.694120e+01 4.098500e-01 4.098500e-01 7.500000e-01 5.000000e-02 5.000000e-02 1.826360e+01 3.338910e-01 3.338910e-01 8.500000e-01 5.000000e-02 5.000000e-02 1.039880e+01 3.415420e-01 3.415420e-01 9.500000e-01 5.000000e-02 5.000000e-02 6.446800e+00 3.866530e-01 3.866530e-01 1.050000e+00 5.000000e-02 5.000000e-02 3.964640e+00 3.560350e-01 3.560350e-01 1.150000e+00 5.000000e-02 5.000000e-02 2.846050e+00 3.426240e-01 3.426240e-01 1.250000e+00 5.000000e-02 5.000000e-02 2.196580e+00 3.311460e-01 3.311460e-01 1.350000e+00 5.000000e-02 5.000000e-02 1.317930e+00 2.946970e-01 2.946970e-01 1.450000e+00 5.000000e-02 5.000000e-02 6.857450e-01 2.424470e-01 2.424470e-01 1.550000e+00 5.000000e-02 5.000000e-02 8.816070e-01 3.332160e-01 3.332160e-01 1.650000e+00 5.000000e-02 5.000000e-02 5.874150e-01 2.937070e-01 2.937070e-01 1.750000e+00 5.000000e-02 5.000000e-02 7.053660e-01 4.301540e-01 4.301540e-01 1.850000e+00 5.000000e-02 5.000000e-02 4.568660e-01 3.230530e-01 3.230530e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y24 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y24 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 3.500000e-01 5.000000e-02 5.000000e-02 8.650690e+01 3.672010e+00 3.672010e+00 4.500000e-01 5.000000e-02 5.000000e-02 5.243730e+01 1.297620e+00 1.297620e+00 5.500000e-01 5.000000e-02 5.000000e-02 3.106830e+01 1.118170e+00 1.118170e+00 6.500000e-01 5.000000e-02 5.000000e-02 1.960520e+01 8.397940e-01 8.397940e-01 7.500000e-01 5.000000e-02 5.000000e-02 1.150940e+01 4.985270e-01 4.985270e-01 8.500000e-01 5.000000e-02 5.000000e-02 6.987870e+00 3.268790e-01 3.268790e-01 9.500000e-01 5.000000e-02 5.000000e-02 4.439180e+00 2.425380e-01 2.425380e-01 1.050000e+00 5.000000e-02 5.000000e-02 3.305780e+00 1.847990e-01 1.847990e-01 1.150000e+00 5.000000e-02 5.000000e-02 1.868160e+00 1.283060e-01 1.283060e-01 1.250000e+00 5.000000e-02 5.000000e-02 1.252160e+00 9.631900e-02 9.631900e-02 1.350000e+00 5.000000e-02 5.000000e-02 8.186510e-01 7.351700e-02 7.351700e-02 1.450000e+00 5.000000e-02 5.000000e-02 5.613590e-01 6.237300e-02 6.237300e-02 1.550000e+00 5.000000e-02 5.000000e-02 3.560040e-01 5.192800e-02 5.192800e-02 1.650000e+00 5.000000e-02 5.000000e-02 1.843260e-01 4.022300e-02 4.022300e-02 1.750000e+00 5.000000e-02 5.000000e-02 1.224180e-01 3.533900e-02 3.533900e-02 1.850000e+00 5.000000e-02 5.000000e-02 1.303250e-01 4.121200e-02 4.121200e-02 1.950000e+00 5.000000e-02 5.000000e-02 6.277800e-02 3.138900e-02 3.138900e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y25 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y25 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 4.500000e-01 5.000000e-02 5.000000e-02 5.113720e+01 1.420480e+00 1.420480e+00 5.500000e-01 5.000000e-02 5.000000e-02 3.204600e+01 1.094670e+00 1.094670e+00 6.500000e-01 5.000000e-02 5.000000e-02 2.061720e+01 6.329530e-01 6.329530e-01 7.500000e-01 5.000000e-02 5.000000e-02 1.216790e+01 4.416650e-01 4.416650e-01 8.500000e-01 5.000000e-02 5.000000e-02 7.001200e+00 2.985320e-01 2.985320e-01 9.500000e-01 5.000000e-02 5.000000e-02 4.525580e+00 2.126300e-01 2.126300e-01 1.050000e+00 5.000000e-02 5.000000e-02 3.106700e+00 1.707600e-01 1.707600e-01 1.150000e+00 5.000000e-02 5.000000e-02 1.991660e+00 1.355150e-01 1.355150e-01 1.250000e+00 5.000000e-02 5.000000e-02 1.040010e+00 1.090230e-01 1.090230e-01 1.350000e+00 5.000000e-02 5.000000e-02 6.507900e-01 1.016360e-01 1.016360e-01 1.450000e+00 5.000000e-02 5.000000e-02 4.306700e-01 1.044530e-01 1.044530e-01 1.550000e+00 5.000000e-02 5.000000e-02 1.280000e-01 9.051400e-02 9.051400e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y26 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y26 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 5.000000e-02 5.000000e-02 1.249130e+02 4.170720e+00 4.170720e+00 3.500000e-01 5.000000e-02 5.000000e-02 7.474840e+01 2.163220e+00 2.163220e+00 4.500000e-01 5.000000e-02 5.000000e-02 5.532390e+01 1.777260e+00 1.777260e+00 5.500000e-01 5.000000e-02 5.000000e-02 3.278720e+01 8.845270e-01 8.845270e-01 6.500000e-01 5.000000e-02 5.000000e-02 1.734440e+01 6.415050e-01 6.415050e-01 7.500000e-01 5.000000e-02 5.000000e-02 1.392150e+01 4.758270e-01 4.758270e-01 8.500000e-01 5.000000e-02 5.000000e-02 8.396030e+00 3.226860e-01 3.226860e-01 9.500000e-01 5.000000e-02 5.000000e-02 4.477570e+00 2.439090e-01 2.439090e-01 1.050000e+00 5.000000e-02 5.000000e-02 2.537130e+00 2.316080e-01 2.316080e-01 1.150000e+00 5.000000e-02 5.000000e-02 2.646420e+00 3.945050e-01 3.945050e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y27 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d01-x01-y27 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 5.000000e-02 5.000000e-02 1.504620e+02 1.373520e+01 1.373520e+01 3.500000e-01 5.000000e-02 5.000000e-02 5.292920e+01 4.224200e+00 4.224200e+00 4.500000e-01 5.000000e-02 5.000000e-02 4.462350e+01 1.211800e+00 1.211800e+00 5.500000e-01 5.000000e-02 5.000000e-02 2.345630e+01 6.513100e-01 6.513100e-01 6.500000e-01 5.000000e-02 5.000000e-02 1.548500e+01 4.151900e-01 4.151900e-01 7.500000e-01 5.000000e-02 5.000000e-02 1.023350e+01 2.721400e-01 2.721400e-01 8.500000e-01 5.000000e-02 5.000000e-02 5.897710e+00 1.876300e-01 1.876300e-01 9.500000e-01 5.000000e-02 5.000000e-02 3.515710e+00 1.429300e-01 1.429300e-01 1.050000e+00 5.000000e-02 5.000000e-02 2.156900e+00 1.286700e-01 1.286700e-01 1.150000e+00 5.000000e-02 5.000000e-02 1.477210e+00 1.182700e-01 1.182700e-01 1.250000e+00 5.000000e-02 5.000000e-02 9.290110e-01 1.065600e-01 1.065600e-01 1.350000e+00 5.000000e-02 5.000000e-02 4.019260e-01 8.569100e-02 8.569100e-02 1.450000e+00 5.000000e-02 5.000000e-02 2.003180e-01 8.958500e-02 8.958500e-02 END YODA_SCATTER2D_V2 + +BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d01-x01-y28 +Variations: [""] +IsRef: 1 +Path: /REF/BRAHMS_2004_I647076/d01-x01-y28 +Title: ~ +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +3.500000e-01 5.000000e-02 5.000000e-02 7.007270e+01 2.231570e+00 2.231570e+00 +4.500000e-01 5.000000e-02 5.000000e-02 3.611770e+01 8.553500e-01 8.553500e-01 +5.500000e-01 5.000000e-02 5.000000e-02 2.569250e+01 5.356000e-01 5.356000e-01 +6.500000e-01 5.000000e-02 5.000000e-02 1.451510e+01 3.463800e-01 3.463800e-01 +7.500000e-01 5.000000e-02 5.000000e-02 9.398100e+00 2.979400e-01 2.979400e-01 +8.500000e-01 5.000000e-02 5.000000e-02 5.123400e+00 2.692800e-01 2.692800e-01 +9.500000e-01 5.000000e-02 5.000000e-02 2.483800e+00 2.694100e-01 2.694100e-01 +1.050000e+00 5.000000e-02 5.000000e-02 2.795700e+00 9.884400e-01 9.884400e-01 +END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y01 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y01 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 5.000000e-02 5.000000e-02 1.708910e+01 8.544560e+00 8.544560e+00 3.500000e-01 5.000000e-02 5.000000e-02 1.920920e+01 3.567070e+00 3.567070e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.727850e+01 1.763480e+00 1.763480e+00 5.500000e-01 5.000000e-02 5.000000e-02 1.496710e+01 9.763400e-01 9.763400e-01 6.500000e-01 5.000000e-02 5.000000e-02 1.135820e+01 5.729400e-01 5.729400e-01 7.500000e-01 5.000000e-02 5.000000e-02 8.156420e+00 3.901700e-01 3.901700e-01 8.500000e-01 5.000000e-02 5.000000e-02 6.372640e+00 2.861400e-01 2.861400e-01 9.500000e-01 5.000000e-02 5.000000e-02 4.274170e+00 2.008100e-01 2.008100e-01 1.050000e+00 5.000000e-02 5.000000e-02 3.460700e+00 1.649800e-01 1.649800e-01 1.150000e+00 5.000000e-02 5.000000e-02 2.641520e+00 1.307700e-01 1.307700e-01 1.250000e+00 5.000000e-02 5.000000e-02 1.814340e+00 1.012600e-01 1.012600e-01 1.350000e+00 5.000000e-02 5.000000e-02 1.393010e+00 8.056000e-02 8.056000e-02 1.450000e+00 5.000000e-02 5.000000e-02 8.167500e-01 5.863000e-02 5.863000e-02 1.550000e+00 5.000000e-02 5.000000e-02 7.230900e-01 5.125000e-02 5.125000e-02 1.650000e+00 5.000000e-02 5.000000e-02 5.337800e-01 4.206000e-02 4.206000e-02 1.750000e+00 5.000000e-02 5.000000e-02 3.648200e-01 3.330000e-02 3.330000e-02 1.850000e+00 5.000000e-02 5.000000e-02 2.471600e-01 2.649000e-02 2.649000e-02 1.950000e+00 5.000000e-02 5.000000e-02 2.202000e-01 2.417000e-02 2.417000e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y02 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y02 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 5.000000e-02 5.000000e-02 1.848240e+01 1.067080e+01 1.067080e+01 3.500000e-01 5.000000e-02 5.000000e-02 2.320910e+01 4.466600e+00 4.466600e+00 4.500000e-01 5.000000e-02 5.000000e-02 2.030190e+01 2.360050e+00 2.360050e+00 5.500000e-01 5.000000e-02 5.000000e-02 1.290070e+01 1.258980e+00 1.258980e+00 6.500000e-01 5.000000e-02 5.000000e-02 1.066070e+01 8.481200e-01 8.481200e-01 7.500000e-01 5.000000e-02 5.000000e-02 7.557390e+00 5.664500e-01 5.664500e-01 8.500000e-01 5.000000e-02 5.000000e-02 6.173100e+00 4.219800e-01 4.219800e-01 9.500000e-01 5.000000e-02 5.000000e-02 4.163380e+00 2.879800e-01 2.879800e-01 1.050000e+00 5.000000e-02 5.000000e-02 3.839380e+00 2.499200e-01 2.499200e-01 1.150000e+00 5.000000e-02 5.000000e-02 2.681650e+00 1.877500e-01 1.877500e-01 1.250000e+00 5.000000e-02 5.000000e-02 2.062530e+00 1.496300e-01 1.496300e-01 1.350000e+00 5.000000e-02 5.000000e-02 1.362170e+00 1.108500e-01 1.108500e-01 1.450000e+00 5.000000e-02 5.000000e-02 7.693400e-01 7.852000e-02 7.852000e-02 1.550000e+00 5.000000e-02 5.000000e-02 6.328800e-01 6.598000e-02 6.598000e-02 1.650000e+00 5.000000e-02 5.000000e-02 5.874800e-01 6.091000e-02 6.091000e-02 1.750000e+00 5.000000e-02 5.000000e-02 3.584100e-01 4.480000e-02 4.480000e-02 1.850000e+00 5.000000e-02 5.000000e-02 2.393200e-01 3.490000e-02 3.490000e-02 1.950000e+00 5.000000e-02 5.000000e-02 2.037800e-01 3.107000e-02 3.107000e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y03 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y03 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 1.500000e-01 5.000000e-02 5.000000e-02 2.892530e+01 6.817750e+00 6.817750e+00 2.500000e-01 5.000000e-02 5.000000e-02 2.208150e+01 2.395070e+00 2.395070e+00 3.500000e-01 5.000000e-02 5.000000e-02 2.427990e+01 1.508680e+00 1.508680e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.729780e+01 1.041210e+00 1.041210e+00 5.500000e-01 5.000000e-02 5.000000e-02 1.324420e+01 7.831480e-01 7.831480e-01 6.500000e-01 5.000000e-02 5.000000e-02 1.112650e+01 6.134210e-01 6.134210e-01 7.500000e-01 5.000000e-02 5.000000e-02 8.749370e+00 4.838410e-01 4.838410e-01 8.500000e-01 5.000000e-02 5.000000e-02 6.335110e+00 3.713710e-01 3.713710e-01 9.500000e-01 5.000000e-02 5.000000e-02 4.385490e+00 2.824940e-01 2.824940e-01 1.050000e+00 5.000000e-02 5.000000e-02 3.428480e+00 2.322060e-01 2.322060e-01 1.150000e+00 5.000000e-02 5.000000e-02 2.432830e+00 1.833810e-01 1.833810e-01 1.250000e+00 5.000000e-02 5.000000e-02 1.699170e+00 1.446430e-01 1.446430e-01 1.350000e+00 5.000000e-02 5.000000e-02 1.230300e+00 1.178410e-01 1.178410e-01 1.450000e+00 5.000000e-02 5.000000e-02 8.523300e-01 9.355600e-02 9.355600e-02 1.550000e+00 5.000000e-02 5.000000e-02 5.860800e-01 7.443200e-02 7.443200e-02 1.650000e+00 5.000000e-02 5.000000e-02 4.628900e-01 6.612700e-02 6.612700e-02 1.750000e+00 5.000000e-02 5.000000e-02 3.416600e-01 8.821700e-02 8.821700e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y04 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y04 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 3.500000e-01 5.000000e-02 5.000000e-02 1.961940e+01 7.695400e-01 7.695400e-01 4.500000e-01 5.000000e-02 5.000000e-02 1.659680e+01 4.870900e-01 4.870900e-01 5.500000e-01 5.000000e-02 5.000000e-02 1.341470e+01 3.900200e-01 3.900200e-01 6.500000e-01 5.000000e-02 5.000000e-02 1.134120e+01 3.490000e-01 3.490000e-01 7.500000e-01 5.000000e-02 5.000000e-02 8.262360e+00 2.941480e-01 2.941480e-01 8.500000e-01 5.000000e-02 5.000000e-02 6.115460e+00 2.472030e-01 2.472030e-01 9.500000e-01 5.000000e-02 5.000000e-02 4.369030e+00 2.045990e-01 2.045990e-01 1.050000e+00 5.000000e-02 5.000000e-02 3.041940e+00 1.664480e-01 1.664480e-01 1.150000e+00 5.000000e-02 5.000000e-02 2.352020e+00 1.426120e-01 1.426120e-01 1.250000e+00 5.000000e-02 5.000000e-02 2.119960e+00 1.314740e-01 1.314740e-01 1.350000e+00 5.000000e-02 5.000000e-02 1.305270e+00 1.001090e-01 1.001090e-01 1.450000e+00 5.000000e-02 5.000000e-02 8.659600e-01 8.371500e-02 8.371500e-02 1.550000e+00 5.000000e-02 5.000000e-02 5.985700e-01 1.058140e-01 1.058140e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y05 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y05 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 5.000000e-02 5.000000e-02 2.345250e+01 4.281830e+00 4.281830e+00 3.500000e-01 5.000000e-02 5.000000e-02 2.113240e+01 1.602040e+00 1.602040e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.471280e+01 9.231630e-01 9.231630e-01 5.500000e-01 5.000000e-02 5.000000e-02 1.225640e+01 5.758580e-01 5.758580e-01 6.500000e-01 5.000000e-02 5.000000e-02 1.111370e+01 3.909770e-01 3.909770e-01 7.500000e-01 5.000000e-02 5.000000e-02 7.735160e+00 2.591380e-01 2.591380e-01 8.500000e-01 5.000000e-02 5.000000e-02 5.988510e+00 1.938850e-01 1.938850e-01 9.500000e-01 5.000000e-02 5.000000e-02 4.296010e+00 1.453150e-01 1.453150e-01 1.050000e+00 5.000000e-02 5.000000e-02 3.206490e+00 1.151800e-01 1.151800e-01 1.150000e+00 5.000000e-02 5.000000e-02 2.281980e+00 9.020300e-02 9.020300e-02 1.250000e+00 5.000000e-02 5.000000e-02 1.859230e+00 7.964100e-02 7.964100e-02 1.350000e+00 5.000000e-02 5.000000e-02 1.368230e+00 8.925400e-02 8.925400e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y06 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y06 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 3.500000e-01 5.000000e-02 5.000000e-02 1.786330e+01 3.994360e+00 3.994360e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.569950e+01 1.135980e+00 1.135980e+00 5.500000e-01 5.000000e-02 5.000000e-02 1.327680e+01 6.085410e-01 6.085410e-01 6.500000e-01 5.000000e-02 5.000000e-02 1.021950e+01 4.165180e-01 4.165180e-01 7.500000e-01 5.000000e-02 5.000000e-02 7.782710e+00 3.038640e-01 3.038640e-01 8.500000e-01 5.000000e-02 5.000000e-02 5.175930e+00 2.107800e-01 2.107800e-01 9.500000e-01 5.000000e-02 5.000000e-02 4.070060e+00 1.622840e-01 1.622840e-01 1.050000e+00 5.000000e-02 5.000000e-02 2.880720e+00 1.298730e-01 1.298730e-01 1.150000e+00 5.000000e-02 5.000000e-02 2.163200e+00 1.242720e-01 1.242720e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y07 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y07 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 5.500000e-01 5.000000e-02 5.000000e-02 9.187640e+00 1.958810e+00 1.958810e+00 6.500000e-01 5.000000e-02 5.000000e-02 7.758600e+00 2.073570e+00 2.073570e+00 7.500000e-01 5.000000e-02 5.000000e-02 6.309230e+00 1.214210e+00 1.214210e+00 8.500000e-01 5.000000e-02 5.000000e-02 5.011930e+00 8.025510e-01 8.025510e-01 9.500000e-01 5.000000e-02 5.000000e-02 3.048820e+00 4.447170e-01 4.447170e-01 1.050000e+00 5.000000e-02 5.000000e-02 2.355430e+00 3.235430e-01 3.235430e-01 1.150000e+00 5.000000e-02 5.000000e-02 1.139200e+00 1.872840e-01 1.872840e-01 1.250000e+00 5.000000e-02 5.000000e-02 1.306570e+00 1.829560e-01 1.829560e-01 1.350000e+00 5.000000e-02 5.000000e-02 1.018810e+00 1.518750e-01 1.518750e-01 1.450000e+00 5.000000e-02 5.000000e-02 4.999100e-01 1.042380e-01 1.042380e-01 1.550000e+00 5.000000e-02 5.000000e-02 3.901300e-01 1.042660e-01 1.042660e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y08 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y08 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 6.500000e-01 5.000000e-02 5.000000e-02 6.622050e+00 3.908870e-01 3.908870e-01 7.500000e-01 5.000000e-02 5.000000e-02 4.627760e+00 3.120040e-01 3.120040e-01 8.500000e-01 5.000000e-02 5.000000e-02 3.467990e+00 3.558080e-01 3.558080e-01 9.500000e-01 5.000000e-02 5.000000e-02 2.419260e+00 3.176640e-01 3.176640e-01 1.050000e+00 5.000000e-02 5.000000e-02 1.585410e+00 2.506760e-01 2.506760e-01 1.150000e+00 5.000000e-02 5.000000e-02 1.621610e+00 3.180250e-01 3.180250e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y09 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y09 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 8.500000e-01 5.000000e-02 5.000000e-02 2.954070e+00 3.530790e-01 3.530790e-01 9.500000e-01 5.000000e-02 5.000000e-02 2.052690e+00 1.966120e-01 1.966120e-01 1.050000e+00 5.000000e-02 5.000000e-02 1.269430e+00 1.376890e-01 1.376890e-01 1.150000e+00 5.000000e-02 5.000000e-02 7.875800e-01 1.000230e-01 1.000230e-01 1.250000e+00 5.000000e-02 5.000000e-02 6.314000e-01 8.220100e-02 8.220100e-02 1.350000e+00 5.000000e-02 5.000000e-02 4.941100e-01 6.852100e-02 6.852100e-02 1.450000e+00 5.000000e-02 5.000000e-02 3.266800e-01 5.299400e-02 5.299400e-02 1.550000e+00 5.000000e-02 5.000000e-02 1.943900e-01 3.887800e-02 3.887800e-02 1.650000e+00 5.000000e-02 5.000000e-02 1.594100e-01 3.478600e-02 3.478600e-02 1.750000e+00 5.000000e-02 5.000000e-02 8.125000e-02 2.449600e-02 2.449600e-02 1.850000e+00 5.000000e-02 5.000000e-02 6.533000e-02 2.177700e-02 2.177700e-02 1.950000e+00 5.000000e-02 5.000000e-02 2.217000e-02 1.280200e-02 1.280200e-02 2.050000e+00 5.000000e-02 5.000000e-02 4.120000e-02 1.842400e-02 1.842400e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y10 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y10 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 7.500000e-01 5.000000e-02 5.000000e-02 4.030230e+00 4.450650e-01 4.450650e-01 8.500000e-01 5.000000e-02 5.000000e-02 2.820810e+00 2.427770e-01 2.427770e-01 9.500000e-01 5.000000e-02 5.000000e-02 2.112680e+00 1.867360e-01 1.867360e-01 1.050000e+00 5.000000e-02 5.000000e-02 1.115160e+00 1.279180e-01 1.279180e-01 1.150000e+00 5.000000e-02 5.000000e-02 5.011800e-01 8.353100e-02 8.353100e-02 1.250000e+00 5.000000e-02 5.000000e-02 6.490500e-01 8.673400e-02 8.673400e-02 1.350000e+00 5.000000e-02 5.000000e-02 5.403700e-01 7.353500e-02 7.353500e-02 1.450000e+00 5.000000e-02 5.000000e-02 3.689900e-01 6.066200e-02 6.066200e-02 1.550000e+00 5.000000e-02 5.000000e-02 1.924500e-01 4.536200e-02 4.536200e-02 1.650000e+00 5.000000e-02 5.000000e-02 1.408100e-01 4.064700e-02 4.064700e-02 1.750000e+00 5.000000e-02 5.000000e-02 1.187600e-01 3.958700e-02 3.958700e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y11 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y11 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 6.500000e-01 5.000000e-02 5.000000e-02 3.592700e+00 5.185620e-01 5.185620e-01 7.500000e-01 5.000000e-02 5.000000e-02 3.436590e+00 2.376430e-01 2.376430e-01 8.500000e-01 5.000000e-02 5.000000e-02 2.472080e+00 1.633590e-01 1.633590e-01 9.500000e-01 5.000000e-02 5.000000e-02 1.589760e+00 1.219290e-01 1.219290e-01 1.050000e+00 5.000000e-02 5.000000e-02 1.249680e+00 1.045030e-01 1.045030e-01 1.150000e+00 5.000000e-02 5.000000e-02 6.273800e-01 7.293200e-02 7.293200e-02 1.250000e+00 5.000000e-02 5.000000e-02 4.522300e-01 6.154000e-02 6.154000e-02 1.350000e+00 5.000000e-02 5.000000e-02 3.822000e-01 5.697500e-02 5.697500e-02 1.450000e+00 5.000000e-02 5.000000e-02 2.748700e-01 5.729200e-02 5.729200e-02 1.550000e+00 5.000000e-02 5.000000e-02 1.312200e-01 3.507000e-02 3.507000e-02 1.650000e+00 5.000000e-02 5.000000e-02 1.095400e-01 3.302600e-02 3.302600e-02 1.750000e+00 5.000000e-02 5.000000e-02 5.292000e-02 2.366700e-02 2.366700e-02 1.850000e+00 5.000000e-02 5.000000e-02 2.295000e-02 1.622700e-02 1.622700e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y12 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y12 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 5.500000e-01 5.000000e-02 5.000000e-02 6.445410e+00 4.186740e-01 4.186740e-01 6.500000e-01 5.000000e-02 5.000000e-02 4.351030e+00 2.387940e-01 2.387940e-01 7.500000e-01 5.000000e-02 5.000000e-02 3.067660e+00 1.745130e-01 1.745130e-01 8.500000e-01 5.000000e-02 5.000000e-02 1.880630e+00 1.270810e-01 1.270810e-01 9.500000e-01 5.000000e-02 5.000000e-02 1.331580e+00 1.098270e-01 1.098270e-01 1.050000e+00 5.000000e-02 5.000000e-02 1.423460e+00 1.258170e-01 1.258170e-01 1.150000e+00 5.000000e-02 5.000000e-02 8.283400e-01 1.019610e-01 1.019610e-01 1.250000e+00 5.000000e-02 5.000000e-02 5.077300e-01 8.347000e-02 8.347000e-02 1.350000e+00 5.000000e-02 5.000000e-02 3.626300e-01 7.402200e-02 7.402200e-02 1.450000e+00 5.000000e-02 5.000000e-02 2.474200e-01 6.612500e-02 6.612500e-02 1.550000e+00 5.000000e-02 5.000000e-02 1.842300e-01 6.513500e-02 6.513500e-02 1.650000e+00 5.000000e-02 5.000000e-02 8.424000e-02 4.863600e-02 4.863600e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y13 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y13 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 5.000000e-02 5.000000e-02 2.129900e+01 9.525190e+00 9.525190e+00 3.500000e-01 5.000000e-02 5.000000e-02 1.338850e+01 2.993750e+00 2.993750e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.604920e+01 1.878420e+00 1.878420e+00 5.500000e-01 5.000000e-02 5.000000e-02 1.530300e+01 1.253670e+00 1.253670e+00 6.500000e-01 5.000000e-02 5.000000e-02 9.904920e+00 6.969100e-01 6.969100e-01 7.500000e-01 5.000000e-02 5.000000e-02 6.788490e+00 4.576800e-01 4.576800e-01 8.500000e-01 5.000000e-02 5.000000e-02 5.904070e+00 3.358710e-01 3.358710e-01 9.500000e-01 5.000000e-02 5.000000e-02 4.128640e+00 2.290160e-01 2.290160e-01 1.050000e+00 5.000000e-02 5.000000e-02 3.420700e+00 1.807890e-01 1.807890e-01 1.150000e+00 5.000000e-02 5.000000e-02 2.145380e+00 1.286710e-01 1.286710e-01 1.250000e+00 5.000000e-02 5.000000e-02 1.702360e+00 1.037950e-01 1.037950e-01 1.350000e+00 5.000000e-02 5.000000e-02 1.171820e+00 7.847100e-02 7.847100e-02 1.450000e+00 5.000000e-02 5.000000e-02 8.489900e-01 6.241900e-02 6.241900e-02 1.550000e+00 5.000000e-02 5.000000e-02 6.069700e-01 4.955900e-02 4.955900e-02 1.650000e+00 5.000000e-02 5.000000e-02 4.631900e-01 4.159600e-02 4.159600e-02 1.750000e+00 5.000000e-02 5.000000e-02 3.675100e-01 3.488300e-02 3.488300e-02 1.850000e+00 5.000000e-02 5.000000e-02 2.219800e-01 2.616100e-02 2.616100e-02 1.950000e+00 5.000000e-02 5.000000e-02 1.714900e-01 2.213900e-02 2.213900e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y14 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y14 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 5.000000e-02 5.000000e-02 2.549520e+01 1.274760e+01 1.274760e+01 3.500000e-01 5.000000e-02 5.000000e-02 1.669880e+01 3.830960e+00 3.830960e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.687840e+01 1.923470e+00 1.923470e+00 5.500000e-01 5.000000e-02 5.000000e-02 1.393790e+01 9.855570e-01 9.855570e-01 6.500000e-01 5.000000e-02 5.000000e-02 9.918000e+00 5.459680e-01 5.459680e-01 7.500000e-01 5.000000e-02 5.000000e-02 7.314350e+00 3.747260e-01 3.747260e-01 8.500000e-01 5.000000e-02 5.000000e-02 5.812430e+00 2.764690e-01 2.764690e-01 9.500000e-01 5.000000e-02 5.000000e-02 4.380020e+00 2.104910e-01 2.104910e-01 1.050000e+00 5.000000e-02 5.000000e-02 3.065660e+00 1.606840e-01 1.606840e-01 1.150000e+00 5.000000e-02 5.000000e-02 2.258630e+00 1.258680e-01 1.258680e-01 1.250000e+00 5.000000e-02 5.000000e-02 1.514000e+00 9.575400e-02 9.575400e-02 1.350000e+00 5.000000e-02 5.000000e-02 1.354670e+00 8.401300e-02 8.401300e-02 1.450000e+00 5.000000e-02 5.000000e-02 1.033480e+00 6.999600e-02 6.999600e-02 1.550000e+00 5.000000e-02 5.000000e-02 6.601100e-01 5.371900e-02 5.371900e-02 1.650000e+00 5.000000e-02 5.000000e-02 4.324000e-01 4.180200e-02 4.180200e-02 1.750000e+00 5.000000e-02 5.000000e-02 3.584700e-01 3.658600e-02 3.658600e-02 1.850000e+00 5.000000e-02 5.000000e-02 1.908100e-01 2.572900e-02 2.572900e-02 1.950000e+00 5.000000e-02 5.000000e-02 1.555500e-01 2.268900e-02 2.268900e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y15 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y15 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 1.500000e-01 5.000000e-02 5.000000e-02 3.267370e+01 1.033230e+01 1.033230e+01 2.500000e-01 5.000000e-02 5.000000e-02 2.487560e+01 2.416130e+00 2.416130e+00 3.500000e-01 5.000000e-02 5.000000e-02 2.196010e+01 1.380620e+00 1.380620e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.712990e+01 9.760660e-01 9.760660e-01 5.500000e-01 5.000000e-02 5.000000e-02 1.265970e+01 7.459790e-01 7.459790e-01 6.500000e-01 5.000000e-02 5.000000e-02 1.157610e+01 6.223380e-01 6.223380e-01 7.500000e-01 5.000000e-02 5.000000e-02 6.911770e+00 4.311440e-01 4.311440e-01 8.500000e-01 5.000000e-02 5.000000e-02 5.775310e+00 3.561210e-01 3.561210e-01 9.500000e-01 5.000000e-02 5.000000e-02 4.623810e+00 2.912730e-01 2.912730e-01 1.050000e+00 5.000000e-02 5.000000e-02 3.400000e+00 2.324190e-01 2.324190e-01 1.150000e+00 5.000000e-02 5.000000e-02 2.236470e+00 1.768080e-01 1.768080e-01 1.250000e+00 5.000000e-02 5.000000e-02 1.730480e+00 1.467770e-01 1.467770e-01 1.350000e+00 5.000000e-02 5.000000e-02 9.003400e-01 1.012960e-01 1.012960e-01 1.450000e+00 5.000000e-02 5.000000e-02 7.687900e-01 8.937000e-02 8.937000e-02 1.550000e+00 5.000000e-02 5.000000e-02 5.547800e-01 7.284600e-02 7.284600e-02 1.650000e+00 5.000000e-02 5.000000e-02 3.723900e-01 5.963100e-02 5.963100e-02 1.750000e+00 5.000000e-02 5.000000e-02 3.043900e-01 8.442300e-02 8.442300e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y16 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y16 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 3.500000e-01 5.000000e-02 5.000000e-02 2.150310e+01 1.004770e+00 1.004770e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.700410e+01 6.038330e-01 6.038330e-01 5.500000e-01 5.000000e-02 5.000000e-02 1.377530e+01 4.474000e-01 4.474000e-01 6.500000e-01 5.000000e-02 5.000000e-02 9.846580e+00 3.407550e-01 3.407550e-01 7.500000e-01 5.000000e-02 5.000000e-02 8.054080e+00 2.919600e-01 2.919600e-01 8.500000e-01 5.000000e-02 5.000000e-02 5.473630e+00 2.306860e-01 2.306860e-01 9.500000e-01 5.000000e-02 5.000000e-02 4.227650e+00 1.973300e-01 1.973300e-01 1.050000e+00 5.000000e-02 5.000000e-02 2.900950e+00 1.587330e-01 1.587330e-01 1.150000e+00 5.000000e-02 5.000000e-02 2.262850e+00 1.364550e-01 1.364550e-01 1.250000e+00 5.000000e-02 5.000000e-02 1.671930e+00 1.137600e-01 1.137600e-01 1.350000e+00 5.000000e-02 5.000000e-02 1.122710e+00 9.047000e-02 9.047000e-02 1.450000e+00 5.000000e-02 5.000000e-02 7.798800e-01 7.760100e-02 7.760100e-02 1.550000e+00 5.000000e-02 5.000000e-02 5.820200e-01 1.028870e-01 1.028870e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y17 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y17 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 5.000000e-02 5.000000e-02 2.209050e+01 4.509210e+00 4.509210e+00 3.500000e-01 5.000000e-02 5.000000e-02 1.879550e+01 1.509690e+00 1.509690e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.623950e+01 7.750710e-01 7.750710e-01 5.500000e-01 5.000000e-02 5.000000e-02 1.292460e+01 4.682090e-01 4.682090e-01 6.500000e-01 5.000000e-02 5.000000e-02 9.693250e+00 3.139950e-01 3.139950e-01 7.500000e-01 5.000000e-02 5.000000e-02 7.465610e+00 2.323940e-01 2.323940e-01 8.500000e-01 5.000000e-02 5.000000e-02 5.426950e+00 1.734460e-01 1.734460e-01 9.500000e-01 5.000000e-02 5.000000e-02 3.877390e+00 1.321410e-01 1.321410e-01 1.050000e+00 5.000000e-02 5.000000e-02 3.076000e+00 1.093700e-01 1.093700e-01 1.150000e+00 5.000000e-02 5.000000e-02 2.109460e+00 8.485500e-02 8.485500e-02 1.250000e+00 5.000000e-02 5.000000e-02 1.521670e+00 7.141600e-02 7.141600e-02 1.350000e+00 5.000000e-02 5.000000e-02 1.156260e+00 8.259000e-02 8.259000e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y18 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y18 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 3.500000e-01 5.000000e-02 5.000000e-02 1.433280e+01 3.830600e+00 3.830600e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.624540e+01 1.335370e+00 1.335370e+00 5.500000e-01 5.000000e-02 5.000000e-02 1.170410e+01 7.122860e-01 7.122860e-01 6.500000e-01 5.000000e-02 5.000000e-02 9.441070e+00 4.868860e-01 4.868860e-01 7.500000e-01 5.000000e-02 5.000000e-02 6.842260e+00 3.404150e-01 3.404150e-01 8.500000e-01 5.000000e-02 5.000000e-02 5.280700e+00 2.570600e-01 2.570600e-01 9.500000e-01 5.000000e-02 5.000000e-02 4.061900e+00 1.965690e-01 1.965690e-01 1.050000e+00 5.000000e-02 5.000000e-02 2.995440e+00 1.612690e-01 1.612690e-01 1.150000e+00 5.000000e-02 5.000000e-02 2.213420e+00 1.631750e-01 1.631750e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y19 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y19 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 5.000000e-02 5.000000e-02 1.657120e+01 5.856120e+00 5.856120e+00 3.500000e-01 5.000000e-02 5.000000e-02 1.752920e+01 1.855870e+00 1.855870e+00 4.500000e-01 5.000000e-02 5.000000e-02 1.378560e+01 1.092960e+00 1.092960e+00 5.500000e-01 5.000000e-02 5.000000e-02 7.515910e+00 1.082920e+00 1.082920e+00 6.500000e-01 5.000000e-02 5.000000e-02 5.845030e+00 1.036250e+00 1.036250e+00 7.500000e-01 5.000000e-02 5.000000e-02 4.314750e+00 6.909120e-01 6.909120e-01 8.500000e-01 5.000000e-02 5.000000e-02 3.198840e+00 4.523850e-01 4.523850e-01 9.500000e-01 5.000000e-02 5.000000e-02 2.679550e+00 3.180040e-01 3.180040e-01 1.050000e+00 5.000000e-02 5.000000e-02 1.736700e+00 2.223620e-01 2.223620e-01 1.150000e+00 5.000000e-02 5.000000e-02 1.277900e+00 1.663680e-01 1.663680e-01 1.250000e+00 5.000000e-02 5.000000e-02 1.005360e+00 1.331630e-01 1.331630e-01 1.350000e+00 5.000000e-02 5.000000e-02 7.066700e-01 1.065340e-01 1.065340e-01 1.450000e+00 5.000000e-02 5.000000e-02 4.666700e-01 1.099950e-01 1.099950e-01 1.550000e+00 5.000000e-02 5.000000e-02 4.026000e-01 2.012990e-01 2.012990e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y20 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y20 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 4.500000e-01 5.000000e-02 5.000000e-02 8.300650e+00 9.921180e-01 9.921180e-01 5.500000e-01 5.000000e-02 5.000000e-02 8.730900e+00 7.486690e-01 7.486690e-01 6.500000e-01 5.000000e-02 5.000000e-02 6.172290e+00 5.832270e-01 5.832270e-01 7.500000e-01 5.000000e-02 5.000000e-02 4.724140e+00 5.064810e-01 5.064810e-01 8.500000e-01 5.000000e-02 5.000000e-02 2.869700e+00 4.099570e-01 4.099570e-01 9.500000e-01 5.000000e-02 5.000000e-02 2.443320e+00 3.912440e-01 3.912440e-01 1.050000e+00 5.000000e-02 5.000000e-02 1.256390e+00 3.047180e-01 3.047180e-01 1.150000e+00 5.000000e-02 5.000000e-02 6.874000e-01 3.437010e-01 3.437010e-01 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y21 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y21 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 8.500000e-01 5.000000e-02 5.000000e-02 2.030980e+00 1.650290e-01 1.650290e-01 9.500000e-01 5.000000e-02 5.000000e-02 1.412310e+00 1.022580e-01 1.022580e-01 1.050000e+00 5.000000e-02 5.000000e-02 9.211600e-01 8.589900e-02 8.589900e-02 1.150000e+00 5.000000e-02 5.000000e-02 6.063400e-01 6.615700e-02 6.615700e-02 1.250000e+00 5.000000e-02 5.000000e-02 4.435600e-01 5.339800e-02 5.339800e-02 1.350000e+00 5.000000e-02 5.000000e-02 2.707800e-01 3.908300e-02 3.908300e-02 1.450000e+00 5.000000e-02 5.000000e-02 1.787700e-01 3.021800e-02 3.021800e-02 1.550000e+00 5.000000e-02 5.000000e-02 1.337200e-01 2.483100e-02 2.483100e-02 1.650000e+00 5.000000e-02 5.000000e-02 7.015000e-02 1.811300e-02 1.811300e-02 1.750000e+00 5.000000e-02 5.000000e-02 4.385000e-02 1.461700e-02 1.461700e-02 1.850000e+00 5.000000e-02 5.000000e-02 3.860000e-02 1.459000e-02 1.459000e-02 1.950000e+00 5.000000e-02 5.000000e-02 2.320000e-02 1.160300e-02 1.160300e-02 2.050000e+00 5.000000e-02 5.000000e-02 1.257000e-02 8.893000e-03 8.893000e-03 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y22 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y22 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 7.500000e-01 5.000000e-02 5.000000e-02 2.876370e+00 1.948120e-01 1.948120e-01 8.500000e-01 5.000000e-02 5.000000e-02 1.966780e+00 1.259100e-01 1.259100e-01 9.500000e-01 5.000000e-02 5.000000e-02 1.049810e+00 8.841000e-02 8.841000e-02 1.050000e+00 5.000000e-02 5.000000e-02 1.120460e+00 9.304900e-02 9.304900e-02 1.150000e+00 5.000000e-02 5.000000e-02 5.913300e-01 7.118800e-02 7.118800e-02 1.250000e+00 5.000000e-02 5.000000e-02 4.347400e-01 6.409900e-02 6.409900e-02 1.350000e+00 5.000000e-02 5.000000e-02 3.051400e-01 5.666400e-02 5.666400e-02 1.450000e+00 5.000000e-02 5.000000e-02 1.050100e-01 3.500200e-02 3.500200e-02 1.550000e+00 5.000000e-02 5.000000e-02 1.751500e-01 4.857700e-02 4.857700e-02 1.650000e+00 5.000000e-02 5.000000e-02 6.577000e-02 3.288500e-02 3.288500e-02 1.750000e+00 5.000000e-02 5.000000e-02 4.074000e-02 2.881000e-02 2.881000e-02 END YODA_SCATTER2D_V2 BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y23 +Variations: [""] IsRef: 1 Path: /REF/BRAHMS_2004_I647076/d02-x01-y23 -Title: +Title: ~ Type: Scatter2D --- # xval xerr- xerr+ yval yerr- yerr+ 6.500000e-01 5.000000e-02 5.000000e-02 3.320390e+00 3.461740e-01 3.461740e-01 7.500000e-01 5.000000e-02 5.000000e-02 2.527710e+00 1.529840e-01 1.529840e-01 8.500000e-01 5.000000e-02 5.000000e-02 1.737180e+00 1.107580e-01 1.107580e-01 9.500000e-01 5.000000e-02 5.000000e-02 9.933800e-01 7.953400e-02 7.953400e-02 1.050000e+00 5.000000e-02 5.000000e-02 9.144000e-01 7.392400e-02 7.392400e-02 1.150000e+00 5.000000e-02 5.000000e-02 4.793500e-01 5.230100e-02 5.230100e-02 1.250000e+00 5.000000e-02 5.000000e-02 3.068100e-01 4.099900e-02 4.099900e-02 1.350000e+00 5.000000e-02 5.000000e-02 2.125400e-01 3.447900e-02 3.447900e-02 1.450000e+00 5.000000e-02 5.000000e-02 1.663600e-01 3.144000e-02 3.144000e-02 1.550000e+00 5.000000e-02 5.000000e-02 6.205000e-02 1.962200e-02 1.962200e-02 1.650000e+00 5.000000e-02 5.000000e-02 7.496000e-02 2.260000e-02 2.260000e-02 1.750000e+00 5.000000e-02 5.000000e-02 2.987000e-02 1.493600e-02 1.493600e-02 END YODA_SCATTER2D_V2 + +BEGIN YODA_SCATTER2D_V2 /REF/BRAHMS_2004_I647076/d02-x01-y24 +Variations: [""] +IsRef: 1 +Path: /REF/BRAHMS_2004_I647076/d02-x01-y24 +Title: ~ +Type: Scatter2D +--- +# xval xerr- xerr+ yval yerr- yerr+ +5.500000e-01 5.000000e-02 5.000000e-02 3.835580e+00 2.739700e-01 2.739700e-01 +6.500000e-01 5.000000e-02 5.000000e-02 2.606860e+00 1.528170e-01 1.528170e-01 +7.500000e-01 5.000000e-02 5.000000e-02 1.840450e+00 1.126330e-01 1.126330e-01 +8.500000e-01 5.000000e-02 5.000000e-02 1.343360e+00 9.098300e-02 9.098300e-02 +9.500000e-01 5.000000e-02 5.000000e-02 8.608900e-01 7.550500e-02 7.550500e-02 +1.050000e+00 5.000000e-02 5.000000e-02 7.215200e-01 7.605400e-02 7.605400e-02 +1.150000e+00 5.000000e-02 5.000000e-02 3.779500e-01 5.763700e-02 5.763700e-02 +1.250000e+00 5.000000e-02 5.000000e-02 3.031300e-01 5.749500e-02 5.749500e-02 +1.350000e+00 5.000000e-02 5.000000e-02 2.117500e-01 5.227600e-02 5.227600e-02 +1.450000e+00 5.000000e-02 5.000000e-02 1.631200e-01 4.799600e-02 4.799600e-02 +1.550000e+00 5.000000e-02 5.000000e-02 4.723000e-02 1.559600e-02 1.559600e-02 +1.650000e+00 5.000000e-02 5.000000e-02 6.040000e-02 3.486700e-02 3.486700e-02 +END YODA_SCATTER2D_V2 diff --git a/include/Rivet/Tools/Correlators.hh b/include/Rivet/Tools/Correlators.hh --- a/include/Rivet/Tools/Correlators.hh +++ b/include/Rivet/Tools/Correlators.hh @@ -1,1136 +1,1136 @@ // -*- C++ -*- #ifndef RIVET_Correlators_HH #define RIVET_Correlators_HH #include "Rivet/Projection.hh" #include "Rivet/Projections/ParticleFinder.hh" #include #include "YODA/Scatter2D.h" #include "Rivet/Analysis.hh" /* File contains tools for calculating flow coefficients using * correlators. * Classes: * Correlators: Calculates single event correlators of a given harmonic. * Cumulants: An additional base class for flow analyses * (use as: class MyAnalysis : public Analysis, Cumulants {};) * Includes a framework for calculating cumulants and flow coefficients * from single event correlators, including automatic handling of statistical * errors. Contains multiple internal sub-classes: * CorBinBase: Base class for correlators binned in event or particle observables. * CorSingleBin: A simple bin for correlators. * CorBin: Has the interface of a simple bin, but does automatic calculation * of statistical errors by a bootstrap method. * ECorrelator: Data type for event averaged correlators. * @author Vytautas Vislavicius, Christine O. Rasmussen, Christian Bierlich. */ namespace Rivet { /* @brief Projection for calculating correlators for flow measurements. * * A projection which calculates Q-vectors and P-vectors, and projects * them out as correlators. Implementation follows the description of * the ''Generic Framework'' * Phys. Rev. C 83 (2011) 044913, arXiv: 1010.0233 * Phys. Rev. C 89 (2014) 064904, arXiv: 1312.3572 * */ class Correlators : public Projection { public: // Constructor. Takes parameters @parm fsp, the Final State // projection correlators should be constructed from, @parm nMaxIn, // the maximal sum of harmonics, eg. for // c_2{2} = {2,-2} = 2 + 2 = 4 // c_2{4} = {2,2,-2,-2} = 2 + 2 + 2 + 2 = 8 // c_4{2} = {4,-4} = 4 + 4 = 8 // c_4{4} = {4,4,-4,-4} = 4 + 4 + 4 + 4 = 16. // @parm pMaxIn is the maximal number of particles // you want to correlate and @parm pTbinEdgesIn is the (lower) // edges of pT bins, the last one the upper edge of the final bin. Correlators(const ParticleFinder& fsp, int nMaxIn = 2, int pMaxIn = 0, vector pTbinEdgesIn = {}); // Constructor which takes a Scatter2D to estimate bin edges. Correlators(const ParticleFinder& fsp, int nMaxIn, int pMaxIn, const Scatter2DPtr hIn); /// @brief Integrated correlator of @parm n harmonic, with the /// number of powers being the size of @parm n. /// Eg. @parm n should be: /// <<2>>_2 => n = {2, -2} /// <<4>>_2 => n = {2, 2, -2, -2} /// <<2>>_4 => n = {4, -4} /// <<4>>_4 => n = {4, 4, -4, 4} and so on. const pair intCorrelator(vector n) const; /// @brief pT differential correlator of @parm n harmonic, with the /// number of powers being the size of @parm n. /// The method can include overflow/underflow bins in the /// beginning/end of the returned vector, by toggling /// @parm overflow = true. const vector> pTBinnedCorrelators(vector n, bool overflow = false) const; /// @brief Integrated correlator of @parm n1 harmonic, with the /// number of powers being the size of @parm n1. This method /// imposes an eta gap, correlating with another phase space, /// where another Correlators projection @parm other should be /// defined. The harmonics of the other phase space is given /// as @parm n2. /// To get eg. integrated <<4>>_2, @parm n1 should be: /// n1 = {2, 2} and n2 = {-2, -2} const pair intCorrelatorGap(const Correlators& other, vector n1, vector n2) const; /// @brief pT differential correlators of @parm n1 harmonic, with the /// number of powers being the size of @parm n1. This method /// imposes an eta gap, correlating with another phase space, /// where another Correlators projection @parm other should be /// defined. The harmonics of the other phase space is given /// as @parm n2. /// To get eg. differential <<4'>_2, @parm n1 should be: /// n1 = {2, 2} and @parm n2: n2 = {-2, -2}. /// To get eg. differential <<2'>>_4, @parm n1 should be: /// n1 = {4} and @parm n2: n2 = {-4}. /// The method can include overflow/underflow /// bins in the beginning/end of the returned vector, by toggling /// @parm overflow = true. const vector> pTBinnedCorrelatorsGap( const Correlators& other, vector n1, vector n2, bool oveflow = false) const; /// @brief Construct a harmonic vectors from @parm n harmonics /// and @parm m number of particles. /// TODO: In C++14 this can be done much nicer with TMP. static vector hVec(int n, int m) { if (m%2 != 0) { cout << "Harmonic Vector: Number of particles " "must be an even number." << endl; return {}; } vector ret = {}; for (int i = 0; i < m; ++i) { if (i < m/2) ret.push_back(n); else ret.push_back(-n); } return ret; } /// @brief Return the maximal values for n, p to be used in the /// constructor of Correlators(xxx, nMax, pMax, xxxx) static pair getMaxValues(vector< vector >& hList){ int nMax = 0, pMax = 0; for (vector h : hList) { int tmpN = 0, tmpP = 0; for ( int i = 0; i < int(h.size()); ++i) { tmpN += abs(h[i]); ++tmpP; } if (tmpN > nMax) nMax = tmpN; if (tmpP > pMax) pMax = tmpP; } return make_pair(nMax,pMax); } // Clone on the heap. DEFAULT_RIVET_PROJ_CLONE(Correlators); protected: // @brief Project function. Loops over array and calculates Q vectors // and P vectors if needed void project(const Event& e); // @brief Compare against other projection. Test if harmonics, pT-bins // and underlying final state are similar. int compare(const Projection& p) const { const Correlators* other = dynamic_cast(&p); if (nMax != other->nMax) return UNDEFINED; if (pMax != other->pMax) return UNDEFINED; if (pTbinEdges != other->pTbinEdges) return UNDEFINED; return mkPCmp(p, "FS"); }; // @brief Calculate correlators from one particle. void fillCorrelators(const Particle& p, const double& weight); // @brief Return a Q-vector. const complex getQ(int n, int p) const { bool isNeg = (n < 0); if (isNeg) return conj( qVec[abs(n)][p] ); else return qVec[n][p]; }; // @brief Return a P-vector. const complex getP(int n, int p, double pT = 0.) const { bool isNeg = (n < 0); map::const_iterator pTitr = pVec.lower_bound(pT); if (pTitr == pVec.end()) return DBL_NAN; if (isNeg) return conj( pTitr->second[abs(n)][p] ); else return pTitr->second[n][p]; }; private: // Find correlators by recursion. Order = M (# of particles), // n's are harmonics, p's are the powers of the weights const complex recCorr(int order, vector n, vector p, bool useP, double pT = 0.) const; // Two-particle correlator Eq. (19) p. 6 // Flag if p-vectors or q-vectors should be used to // calculate the correlator. const complex twoPartCorr(int n1, int n2, int p1 = 1, int p2 = 1, double pT = 0., bool useP = false) const; // Set elements in vectors to zero. void setToZero(); // Shorthands for setting and comparing to zero. const complex _ZERO = {0., 0.}; const double _TINY = 1e-10; // Shorthand typedefs for vec>. typedef vector< vector> > Vec2D; // Define Q-vectors and p-vectors Vec2D qVec; // Q[n][p] map pVec; // p[pT][n][p] // The max values of n and p to be calculated. int nMax, pMax; // pT bin-edges vector pTbinEdges; bool isPtDiff; /// End class Correlators }; /// @brief Tools for flow analyses. /// The following are helper classes to construct event averaged correlators /// as well as cummulants and flow coefficents from the basic event // correlators defined above. They are all encapsulated in a Cumulants // class, which can be used as a(nother) base class for flow analyses, // to ensure access. class CumulantAnalysis : public Analysis { private: // Number of bins used for bootstrap calculation of statistical // uncertainties. It is hard coded, and shout NOT be changed unless there // are good reasons to do so. static const int BOOT_BINS = 9; // Enum for choosing the method of error calculation. enum Error { VARIANCE, ENVELOPE }; // The desired error method. Can be changed in the analysis constructor // by setting it appropriately. Error errorMethod; /// @brief Base class for correlator bins. class CorBinBase { public: CorBinBase() {} virtual ~CorBinBase() {}; // Derived class should have fill and mean defined. virtual void fill(const pair& cor, const double& weight) = 0; virtual const double mean() const = 0; }; /// @brief The CorSingleBin is the basic quantity filled in an ECorrelator. /// It is a simple counter with an even simpler structure than normal /// YODA type DBNs, but added functionality to test for out of /// bounds correlators. class CorSingleBin : public CorBinBase { public: /// @brief The default constructor. CorSingleBin() : _sumWX(0.), _sumW(0.), _sumW2(0.), _numEntries(0.) {} ~CorSingleBin() {} /// @brief Fill a correlator bin with the return type from a /// Correlator (a pair giving numerator and denominator of _event). void fill(const pair& cor, const double& weight) { // Test if denominator for the single event average is zero. if (cor.second < 1e-10) return; // The single event average is then cor.first / cor.second. // With weights this becomes just: _sumWX += cor.first * weight; _sumW += weight * cor.second; _sumW2 += weight * weight * cor.second * cor.second; _numEntries += 1.; } const double mean() const { if (_sumW < 1e-10) return 0; return _sumWX / _sumW; } // @brief Sum of weights. const double sumW() const { return _sumW; } const double sumW2() const { return _sumW2; } // @brief Sum of weight * X. const double sumWX() const { return _sumWX; } // @brief Number of entries. const double numEntries() const { return _numEntries; } void addContent(double ne, double sw, double sw2, double swx) { _numEntries += ne; _sumW += sw; _sumW2 += sw2; _sumWX += swx; } private: double _sumWX, _sumW, _sumW2, _numEntries; }; // End of CorSingleBin sub-class. /// @brief The CorBin is the basic bin quantity in ECorrelators. /// It consists of several CorSingleBins, to facilitate /// bootstrapping calculation of statistical uncertainties. class CorBin : public CorBinBase { public: // @brief The constructor. nBins signifies the period of the bootstrap // calculation, and should never be changed here, but in its definition // above -- and only if there are good reasons to do so. CorBin() : binIndex(0), nBins(BOOT_BINS) { for(size_t i = 0; i < nBins; ++i) bins.push_back(CorSingleBin()); } // Destructor must be implemented. ~CorBin() {} // @brief Fill the correct underlying bin and take a step. void fill(const pair& cor, const double& weight) { // Test if denominator for the single event average is zero. if (cor.second < 1e-10) return; // Fill the correct bin. bins[binIndex].fill(cor, weight); if (binIndex == nBins - 1) binIndex = 0; else ++binIndex; } // @brief Calculate the total sample mean with all // available statistics. const double mean() const { double sow = 0; double sowx = 0; for(auto b : bins) { if (b.sumW() < 1e-10) continue; sow += b.sumW(); sowx += b.sumWX(); } return sowx / sow; } // @brief Return a copy of the bins. const vector getBins() const { return bins; } // @brief Return the bins as pointers to the base class. template const vector getBinPtrs() { vector ret(bins.size()); transform(bins.begin(), bins.end(), ret.begin(), [](CorSingleBin& b) {return &b;}); return ret; } private: vector bins; size_t binIndex; size_t nBins; }; // End of CorBin sub-class. public: /// @brief The ECorrelator is a helper class to calculate all event /// averages of correlators, in order to construct cumulants. /// It can be binned in any variable. class ECorrelator { public: /// @brief Constructor. Takes as argument the desired harmonic and number /// of correlated particles as a generic framework style vector, eg, /// {2, -2} for <<2>>_2, no binning. /// TODO: Implement functionality for this if needed. //ECorrelator(vector h) : h1(h), h2({}), // binX(0), binContent(0), reference() { //}; /// @brief Constructor. Takes as argument the desired harmonic and number /// of correlated particles as a generic framework style vector, eg, /// {2, -2} for <<2>>_2 and binning. ECorrelator(vector h, vector binIn) : h1(h), h2({}), binX(binIn), binContent(binIn.size() - 1), reference() {}; /// @brief Constructor for gapped correlator. Takes as argument the /// desired harmonics for the two final states, and binning. ECorrelator(vector h1In, vector h2In, vector binIn) : h1(h1In), h2(h2In), binX(binIn), binContent(binIn.size() - 1), reference() {}; /// @brief Fill the appropriate bin given an input (per event) /// observable, eg. centrality. void fill(const double& obs, const Correlators& c, const double weight = 1.0) { int index = getBinIndex(obs); if (index < 0) return; binContent[index].fill(c.intCorrelator(h1), weight); } /// @brief Fill the appropriate bin given an input (per event) /// observable, eg. centrality. Using a rapidity gap between /// two Correlators. void fill(const double& obs, const Correlators& c1, const Correlators& c2, const double weight = 1.0) { if (!h2.size()) { cout << "Trying to fill gapped correlator, but harmonics behind " "the gap (h2) are not given!" << endl; return; } int index = getBinIndex(obs); if (index < 0) return; binContent[index].fill(c1.intCorrelatorGap(c2, h1, h2), weight); } /// @brief Fill the bins with the appropriate correlator, taking the /// binning directly from the Correlators object, and filling also the /// reference flow. void fill(const Correlators& c, const double& weight = 1.0) { vector< pair > diffCorr = c.pTBinnedCorrelators(h1); // We always skip overflow when calculating the all event average. if (diffCorr.size() != binX.size() - 1) cout << "Tried to fill event with wrong binning (ungapped)" << endl; for (size_t i = 0; i < diffCorr.size(); ++i) { int index = getBinIndex(binX[i]); if (index < 0) return; binContent[index].fill(diffCorr[i], weight); } reference.fill(c.intCorrelator(h1), weight); } /// @brief Fill bins with the appropriate correlator, taking the binning /// directly from the Correlators object, and also the reference flow. /// Using a rapidity gap between two Correlators. void fill(const Correlators& c1, const Correlators& c2, const double& weight = 1.0) { if (!h2.size()) { cout << "Trying to fill gapped correlator, but harmonics behind " "the gap (h2) are not given!" << endl; return; } vector< pair > diffCorr = c1.pTBinnedCorrelatorsGap(c2, h1, h2); // We always skip overflow when calculating the all event average. if (diffCorr.size() != binX.size() - 1) cout << "Tried to fill event with wrong binning (gapped)" << endl; for (size_t i = 0; i < diffCorr.size(); ++i) { int index = getBinIndex(binX[i]); if (index < 0) return; binContent[index].fill(diffCorr[i], weight); } reference.fill(c1.intCorrelatorGap(c2, h1, h2), weight); } /// @brief Get a copy of the bin contents. const vector getBins() const { return binContent; } // @brief Return the bins as pointers to the base class. const vector getBinPtrs() { vector ret(binContent.size()); transform(binContent.begin(), binContent.end(), ret.begin(), [](CorBin& b) {return &b;}); return ret; } /// @brief Get a copy of the bin x-values. const vector getBinX() const { return binX; } /// @brief Get a copy of the @ret h1 harmonic vector. const vector getH1() const { return h1; } /// @brief Get a copy of the @ret h2 harmonic vector. const vector getH2() const { return h2; } /// @brief Replace reference flow bin with another reference /// flow bin, eg. calculated in another phase space or with /// other pid. void setReference(CorBin refIn) { reference = refIn; } /// @brief Extract the reference flow from a differential event /// averaged correlator. const CorBin getReference() const { if (reference.mean() < 1e-10) cout << "Warning: ECorrelator, reference bin is zero." << endl; return reference; } /// @brief set the @parm prIn list of profile histograms associated /// with the internal bins. Is called automatically when booking, no /// need to call it yourself. void setProfs(list prIn) { profs = prIn; } /// @brief Fill bins with content from preloaded histograms. void fillFromProfs() { list::iterator hItr = profs.begin(); auto refs = reference.getBinPtrs(); for (size_t i = 0; i < profs.size(); ++i, ++hItr) { for (size_t j = 0; j < binX.size() - 1; ++j) { const YODA::ProfileBin1D& pBin = (*hItr)->binAt(binX[j]); auto tmp = binContent[j].getBinPtrs(); tmp[i]->addContent(pBin.numEntries(), pBin.sumW(), pBin.sumW2(), pBin.sumWY()); } // Get the reference flow from the underflow bin of the histogram. const YODA::Dbn2D& uBin = (*hItr)->underflow(); refs[i]->addContent(uBin.numEntries(), uBin.sumW(), uBin.sumW2(), uBin.sumWY()); } // End loop of bootstrapped correlators. } /// @brief begin() iterator for the list of associated profile histograms. list::iterator profBegin() { return profs.begin(); } /// @brief end() iterator for the list of associated profile histograms. list::iterator profEnd() { return profs.end(); } private: /// @brief Get correct bin index for a given @parm obs value. const int getBinIndex(const double& obs) const { // Find the correct index of binContent. // If we are in overflow, just skip. if (obs >= binX.back()) return -1; // If we are in underflow, ditto. if (obs < binX[0]) return -1; int index = 0; for (int i = 0, N = binX.size() - 1; i < N; ++i, ++index) if (obs >= binX[i] && obs < binX[i + 1]) break; return index; } // The harmonics vectors. vector h1; vector h2; // The bins. vector binX; vector binContent; // The reference flow. CorBin reference; // The profile histograms associated with the CorBins for streaming. list profs; }; // End of ECorrelator sub-class. // @brief Get the correct max N and max P for the set of // booked correlators. const pair getMaxValues() const { vector< vector> harmVecs; for ( auto eItr = eCorrPtrs.begin(); eItr != eCorrPtrs.end(); ++eItr) { vector h1 = (*eItr)->getH1(); vector h2 = (*eItr)->getH2(); if (h1.size() > 0) harmVecs.push_back(h1); if (h2.size() > 0) harmVecs.push_back(h2); } if (harmVecs.size() == 0) { cout << "Warning: You tried to extract max values from harmonic " "vectors, but have not booked any." << endl; return pair(); } return Correlators::getMaxValues(harmVecs); } // Typedeffing shared pointer to ECorrelator. typedef shared_ptr ECorrPtr; // @brief Book an ECorrelator in the same way as a histogram. ECorrPtr bookECorrelator(const string name, const vector& h, const Scatter2DPtr hIn) { vector binIn; for (auto b : hIn->points()) binIn.push_back(b.xMin()); binIn.push_back(hIn->points().back().xMax()); ECorrPtr ecPtr = ECorrPtr(new ECorrelator(h, binIn)); list eCorrProfs; for (int i = 0; i < BOOT_BINS; ++i) { Profile1DPtr tmp = bookProfile1D(name+"-"+to_string(i),*hIn); - tmp->setPath(this->name()+"/FINAL/" + name+"-"+to_string(i)); + tmp->setPath("/TMP/"+this->name()+"/FINAL/" + name+"-"+to_string(i)); //tmp->setPath(tmp->path()+"FINAL"); eCorrProfs.push_back(tmp); } ecPtr->setProfs(eCorrProfs); eCorrPtrs.push_back(ecPtr); return ecPtr; } // @brief Book a gapped ECorrelator with two harmonic vectors. ECorrPtr bookECorrelator(const string name, const vector& h1, const vector& h2, const Scatter2DPtr hIn ) { vector binIn; for (auto b : hIn->points()) binIn.push_back(b.xMin()); binIn.push_back(hIn->points().back().xMax()); ECorrPtr ecPtr = ECorrPtr(new ECorrelator(h1, h2, binIn)); list eCorrProfs; for (int i = 0; i < BOOT_BINS; ++i) { Profile1DPtr tmp = bookProfile1D(name+"-"+to_string(i),*hIn); - tmp->setPath(this->name()+"/FINAL/" + name+"-"+to_string(i)); + tmp->setPath("/TMP/"+this->name()+"/FINAL/" + name+"-"+to_string(i)); //tmp->setPath(tmp->path()+"FINAL"); eCorrProfs.push_back(tmp); } ecPtr->setProfs(eCorrProfs); eCorrPtrs.push_back(ecPtr); return ecPtr; } // @brief Short hand for gapped correlators which splits the harmonic // vector into negative and positive components. ECorrPtr bookECorrelatorGap (const string name, const vector& h, const Scatter2DPtr hIn) { const vector h1(h.begin(), h.begin() + h.size() / 2); const vector h2(h.begin() + h.size() / 2, h.end()); return bookECorrelator(name, h1, h2, hIn); } // @brief Templated version of correlator booking which takes // @parm N desired harmonic and @parm M number of particles. template ECorrPtr bookECorrelator(const string name, const Scatter2DPtr hIn) { return bookECorrelator(name, Correlators::hVec(N, M), hIn); } // @brief Templated version of gapped correlator booking which takes // @parm N desired harmonic and @parm M number of particles. template ECorrPtr bookECorrelatorGap(const string name, const Scatter2DPtr hIn) { return bookECorrelatorGap(name, Correlators::hVec(N, M), hIn); } // @brief The stream method MUST be called in finalize() if one wants to stream // correlators to the yoda file, in order to do reentrant finalize // (ie. multi-histogram merging) for the analysis. void stream() { for (auto ecItr = eCorrPtrs.begin(); ecItr != eCorrPtrs.end(); ++ecItr){ (*ecItr)->fillFromProfs(); corrPlot(list((*ecItr)->profBegin(), (*ecItr)->profEnd()), *ecItr); } } private: /// Bookkeeping of the event averaged correlators. list eCorrPtrs; public: // @brief Constructor. Use CumulantAnalysis as base class for the // analysis to have access to functionality. CumulantAnalysis (string n) : Analysis(n), errorMethod(VARIANCE) {}; // @brief Helper method for turning correlators into Scatter2Ds. // Takes @parm h a pointer to the resulting Scatter2D, @parm binx // the x-bins and a function @parm func defining the transformation. // Makes no attempt at statistical errors. // See usage in the methods below. // Can also be used directly in the analysis if a user wants to // perform an unforseen transformation from correlators to // Scatter2D. template static void fillScatter(Scatter2DPtr h, vector& binx, T func) { vector points; for (int i = 0, N = binx.size() - 1; i < N; ++i) { double xMid = (binx[i] + binx[i + 1]) / 2.0; double xeMin = fabs(xMid - binx[i]); double xeMax = fabs(xMid - binx[i + 1]); double yVal = func(i); if (std::isnan(yVal)) yVal = 0.; double yErr = 0; points.push_back(YODA::Point2D(xMid, yVal, xeMin, xeMax, yErr, yErr)); } h->reset(); h->points().clear(); for (int i = 0, N = points.size(); i < N; ++i) h->addPoint(points[i]); } // @brief Helper method for turning correlators into Scatter2Ds // with error estimates. // Takes @parm h a pointer to the resulting Scatter2D, @parm binx // the x-bins, a function @parm func defining the transformation // and a vector of errors @parm err. // See usage in the methods below. // Can also be used directly in the analysis if a user wants to // perform an unforseen transformation from correlators to // Scatter2D. template const void fillScatter(Scatter2DPtr h, vector& binx, F func, vector >& yErr) const { vector points; for (int i = 0, N = binx.size() - 1; i < N; ++i) { double xMid = (binx[i] + binx[i + 1]) / 2.0; double xeMin = fabs(xMid - binx[i]); double xeMax = fabs(xMid - binx[i + 1]); double yVal = func(i); if (std::isnan(yVal)) points.push_back(YODA::Point2D(xMid, 0., xeMin, xeMax,0., 0.)); else points.push_back(YODA::Point2D(xMid, yVal, xeMin, xeMax, yErr[i].first, yErr[i].second)); } h->reset(); h->points().clear(); for (int i = 0, N = points.size(); i < N; ++i) h->addPoint(points[i]); } /// @brief Take the @parm n th power of all points in @parm hIn, /// and put the result in @parm hOut. Optionally put a /// @parm k constant below the root. static void nthPow(Scatter2DPtr hOut, const Scatter2DPtr hIn, const double& n, const double& k = 1.0) { if (n == 0 || n == 1) { cout << "Error: Do not take the 0th or 1st power of a Scatter2D," " use scale instead." << endl; return; } if (hIn->points().size() != hOut->points().size()) { cout << "nthRoot: Scatterplots: " << hIn->name() << " and " << hOut->name() << " not initialized with same length." << endl; return; } vector points; // The error pre-factor is k^(1/n) / n by Taylors formula. double eFac = pow(k,1./n)/n; for (auto b : hIn->points()) { double yVal = pow(k * b.y(),n); if (std::isnan(yVal)) points.push_back(YODA::Point2D(b.x(), 0., b.xErrMinus(), b.xErrPlus(), 0, 0 )); else { double yemin = abs(eFac * pow(yVal,1./(n - 1.))) * b.yErrMinus(); if (std::isnan(yemin)) yemin = b.yErrMinus(); double yemax = abs(eFac * pow(yVal,1./(n - 1.))) * b.yErrPlus(); if (std::isnan(yemax)) yemax = b.yErrPlus(); points.push_back(YODA::Point2D(b.x(), yVal, b.xErrMinus(), b.xErrPlus(), yemin, yemax )); } } hOut->reset(); hOut->points().clear(); for (int i = 0, N = points.size(); i < N; ++i) hOut->addPoint(points[i]); } /// @brief Take the @parm n th power of all points in @parm h, /// and put the result back in the same Scatter2D. Optionally put a /// @parm k constant below the root. static void nthPow(Scatter2DPtr h, const double& n, const double& k = 1.0) { if (n == 0 || n == 1) { cout << "Error: Do not take the 0th or 1st power of a Scatter2D," " use scale instead." << endl; return; } vector points; vector pIn = h->points(); // The error pre-factor is k^(1/n) / n by Taylors formula. double eFac = pow(k,1./n)/n; for (auto b : pIn) { double yVal = pow(k * b.y(),n); if (std::isnan(yVal)) points.push_back(YODA::Point2D(b.x(), 0., b.xErrMinus(), b.xErrPlus(), 0, 0 )); else { double yemin = abs(eFac * pow(yVal,1./(n - 1.))) * b.yErrMinus(); if (std::isnan(yemin)) yemin = b.yErrMinus(); double yemax = abs(eFac * pow(yVal,1./(n - 1.))) * b.yErrPlus(); if (std::isnan(yemax)) yemax = b.yErrPlus(); points.push_back(YODA::Point2D(b.x(), yVal, b.xErrMinus(), b.xErrPlus(), yemin, yemax )); } } h->reset(); h->points().clear(); for (int i = 0, N = points.size(); i < N; ++i) h->addPoint(points[i]); } // @brief Calculate the bootstrapped sample variance for the observable // given by correlators and the transformation @parm func. template static pair sampleVariance(T func) { // First we calculate the mean (two pass calculation). double avg = 0.; for (int i = 0; i < BOOT_BINS; ++i) avg += func(i); avg /= double(BOOT_BINS); // Then we find the variance. double var = 0.; for (int i = 0; i < BOOT_BINS; ++i) var += pow(func(i) - avg, 2.); var /= (double(BOOT_BINS) - 1); return pair(var, var); } // @brief Calculate the bootstrapped sample envelope for the observable // given by correlators and the transformation @parm func. template static pair sampleEnvelope(T func) { // First we calculate the mean. double avg = 0.; for (int i = 0; i < BOOT_BINS; ++i) avg += func(i); avg /= double(BOOT_BINS); double yMax = avg; double yMin = avg; // Then we find the envelope using the mean as initial value. for (int i = 0; i < BOOT_BINS; ++i) { double yVal = func(i); if (yMin > yVal) yMin = yVal; else if (yMax < yVal) yMax = yVal; } return pair(fabs(avg - yMin), fabs(yMax - avg)); } // @brief Selection method for which sample error to use, given // in the constructor. template const pair sampleError(T func) const { if (errorMethod == VARIANCE) return sampleVariance(func); else if (errorMethod == ENVELOPE) return sampleEnvelope(func); else cout << "Error: Error method not found!" << endl; return pair(0.,0.); } // @brief Two particle integrated cn. const void cnTwoInt(Scatter2DPtr h, ECorrPtr e2) const { vector bins = e2->getBins(); vector binx = e2->getBinX(); // Assert bin size. if (binx.size() - 1 != bins.size()){ cout << "cnTwoInt: Bin size (x,y) differs!" << endl; return; } vector binPtrs; // The mean value of the cumulant. auto cn = [&] (int i) { return binPtrs[i]->mean(); }; // Error calculation. vector > yErr; for (int j = 0, N = bins.size(); j < N; ++j) { binPtrs = bins[j].getBinPtrs(); yErr.push_back(sampleError(cn)); } binPtrs = e2->getBinPtrs(); fillScatter(h, binx, cn, yErr); } // @brief Two particle integrated vn. const void vnTwoInt(Scatter2DPtr h, ECorrPtr e2) const { cnTwoInt(h, e2); nthPow(h, 0.5); } // @brief Put an event averaged correlator into a Scatter2D. // Reduces to cnTwoInt, but better with a proper name. const void corrPlot(Scatter2DPtr h, ECorrPtr e) const { cnTwoInt(h, e); } // @brief Put an event averaged correlator into Profile1Ds, one // for each bootstrapping bin. const void corrPlot(list profs, ECorrPtr e) const { vector corBins = e->getBins(); vector binx = e->getBinX(); auto ref = e->getReference(); auto refBins = ref.getBinPtrs(); // Assert bin size. if (binx.size() - 1 != corBins.size()){ cout << "corrPlot: Bin size (x,y) differs!" << endl; return; } list::iterator hItr = profs.begin(); // Loop over the boostrapped correlators. for (size_t i = 0; i < profs.size(); ++i, ++hItr) { vector profBins; // Numbers for the summary distribution double ne = 0., sow = 0., sow2 = 0.; for (size_t j = 0, N = binx.size() - 1; j < N; ++j) { vector binPtrs = corBins[j].getBinPtrs(); // Construct bin of the profiled quantities. We have no information // (and no desire to add it) of sumWX of the profile, so really // we should use a Dbn1D - but that does not work for Profile1D's. profBins.push_back( YODA::ProfileBin1D((*hItr)->bin(j).xEdges(), YODA::Dbn2D( binPtrs[i]->numEntries(), binPtrs[i]->sumW(), binPtrs[i]->sumW2(), 0., 0., binPtrs[i]->sumWX(), 0, 0))); ne += binPtrs[i]->numEntries(); sow += binPtrs[i]->sumW(); sow2 += binPtrs[i]->sumW2(); } // Reset the bins of the profiles. (*hItr)->reset(); (*hItr)->bins().clear(); // Add our new bins. // The total distribution (*hItr)->setTotalDbn(YODA::Dbn2D(ne,sow,sow2,0.,0.,0.,0.,0.)); // The bins. for (int j = 0, N = profBins.size(); j < N; ++j) (*hItr)->addBin(profBins[j]); // The reference flow in the underflow bin. (*hItr)->setUnderflow(YODA::Dbn2D(refBins[i]->numEntries(), refBins[i]->sumW(), refBins[i]->sumW2(), 0., 0., refBins[i]->sumWX(), 0., 0.)); } // End loop of bootstrapped correlators. } // @brief Four particle integrated cn. const void cnFourInt(Scatter2DPtr h, ECorrPtr e2, ECorrPtr e4) const { auto e2bins = e2->getBins(); auto e4bins = e4->getBins(); auto binx = e2->getBinX(); if (binx.size() - 1 != e2bins.size()){ cout << "cnFourInt: Bin size (x,y) differs!" << endl; return; } if (binx != e4->getBinX()) { cout << "Error in cnFourInt: Correlator x-binning differs!" << endl; return; } vector e2binPtrs; vector e4binPtrs; auto cn = [&] (int i) { double e22 = e2binPtrs[i]->mean() * e2binPtrs[i]->mean(); return e4binPtrs[i]->mean() - 2. * e22; }; // Error calculation. vector > yErr; for (int j = 0, N = e2bins.size(); j < N; ++j) { e2binPtrs = e2bins[j].getBinPtrs(); e4binPtrs = e4bins[j].getBinPtrs(); yErr.push_back(sampleError(cn)); } // Put the bin ptrs back in place. e2binPtrs = e2->getBinPtrs(); e4binPtrs = e4->getBinPtrs(); fillScatter(h, binx, cn, yErr); } // @brief Four particle integrated vn const void vnFourInt(Scatter2DPtr h, ECorrPtr e2, ECorrPtr e4) const { cnFourInt(h, e2, e4); nthPow(h, 0.25, -1.0); } // @brief Six particle integrated cn. const void cnSixInt(Scatter2DPtr h, ECorrPtr e2, ECorrPtr e4, ECorrPtr e6) const { auto e2bins = e2->getBins(); auto e4bins = e4->getBins(); auto e6bins = e6->getBins(); auto binx = e2->getBinX(); if (binx.size() - 1 != e2bins.size()){ cout << "cnSixInt: Bin size (x,y) differs!" << endl; return; } if (binx != e4->getBinX() || binx != e6->getBinX()) { cout << "Error in cnSixInt: Correlator x-binning differs!" << endl; return; } vector e2binPtrs; vector e4binPtrs; vector e6binPtrs; auto cn = [&] (int i) { double e23 = pow(e2binPtrs[i]->mean(), 3.0); return e6binPtrs[i]->mean() - 9.*e2binPtrs[i]->mean()*e4binPtrs[i]->mean() + 12.*e23; }; // Error calculation. vector > yErr; for (int j = 0, N = e2bins.size(); j < N; ++j) { e2binPtrs = e2bins[j].getBinPtrs(); e4binPtrs = e4bins[j].getBinPtrs(); e6binPtrs = e6bins[j].getBinPtrs(); yErr.push_back(sampleError(cn)); } // Put the bin ptrs back in place. e2binPtrs = e2->getBinPtrs(); e4binPtrs = e4->getBinPtrs(); e6binPtrs = e6->getBinPtrs(); fillScatter(h, binx, cn, yErr); } // @brief Six particle integrated vn const void vnSixInt(Scatter2DPtr h, ECorrPtr e2, ECorrPtr e4, ECorrPtr e6) const { cnSixInt(h, e2, e4, e6); nthPow(h, 1./6., 0.25); } // @brief Eight particle integrated cn. const void cnEightInt(Scatter2DPtr h, ECorrPtr e2, ECorrPtr e4, ECorrPtr e6, ECorrPtr e8) const { auto e2bins = e2->getBins(); auto e4bins = e4->getBins(); auto e6bins = e6->getBins(); auto e8bins = e8->getBins(); auto binx = e2->getBinX(); if (binx.size() - 1 != e2bins.size()){ cout << "cnEightInt: Bin size (x,y) differs!" << endl; return; } if (binx != e4->getBinX() || binx != e6->getBinX() || binx != e8->getBinX()) { cout << "Error in cnEightInt: Correlator x-binning differs!" << endl; return; } vector e2binPtrs; vector e4binPtrs; vector e6binPtrs; vector e8binPtrs; auto cn = [&] (int i ) { double e22 = e2binPtrs[i]->mean() * e2binPtrs[i]->mean(); double e24 = e22 * e22; double e42 = e4binPtrs[i]->mean() * e4binPtrs[i]->mean(); return e8binPtrs[i]->mean() - 16. * e6binPtrs[i]->mean() * e2binPtrs[i]->mean() - 18. * e42 + 144. * e4binPtrs[i]->mean()*e22 - 144. * e24; }; // Error calculation. vector > yErr; for (int j = 0, N = e2bins.size(); j < N; ++j) { e2binPtrs = e2bins[j].getBinPtrs(); e4binPtrs = e4bins[j].getBinPtrs(); e6binPtrs = e6bins[j].getBinPtrs(); e8binPtrs = e8bins[j].getBinPtrs(); yErr.push_back(sampleError(cn)); } // Put the bin ptrs back in place. e2binPtrs = e2->getBinPtrs(); e4binPtrs = e4->getBinPtrs(); e6binPtrs = e6->getBinPtrs(); e8binPtrs = e8->getBinPtrs(); fillScatter(h, binx, cn, yErr); } // @brief Eight particle integrated vn const void vnEightInt(Scatter2DPtr h, ECorrPtr e2, ECorrPtr e4, ECorrPtr e6, ECorrPtr e8) const { cnEightInt(h, e2, e4, e6, e8); nthPow(h, 1./8., -1./33.); } // @brief Two particle differential vn. const void vnTwoDiff(Scatter2DPtr h, ECorrPtr e2Dif) const { auto e2bins = e2Dif->getBins(); auto ref = e2Dif->getReference(); auto binx = e2Dif->getBinX(); if (binx.size() -1 != e2bins.size()) { cout << "vnTwoDif: Bin size (x,y) differs!" << endl; return; } vector e2binPtrs; vector refPtrs; auto vn = [&] (int i) { // Test reference flow. if (ref.mean() <= 0) return 0.; return e2binPtrs[i]->mean() / sqrt(ref.mean()); }; // We need here a seperate error function, as we don't // iterate over the reference flow. auto vnerr = [&] (int i) { // Test reference flow. if (refPtrs[i]->mean() <=0) return 0.; return e2binPtrs[i]->mean() / sqrt(refPtrs[i]->mean()); }; // Error calculation. vector > yErr; refPtrs = ref.getBinPtrs(); for (int j = 0, N = e2bins.size(); j < N; ++j) { e2binPtrs = e2bins[j].getBinPtrs(); yErr.push_back(sampleError(vnerr)); } // Put the e2binPtrs back in place. e2binPtrs = e2Dif->getBinPtrs(); fillScatter(h, binx, vn); } // @brief Four particle differential vn. const void vnFourDiff(Scatter2DPtr h, ECorrPtr e2Dif, ECorrPtr e4Dif) const { auto e2bins = e2Dif->getBins(); auto e4bins = e4Dif->getBins(); auto ref2 = e2Dif->getReference(); auto ref4 = e4Dif->getReference(); auto binx = e2Dif->getBinX(); if (binx.size() - 1 != e2bins.size()){ cout << "vnFourDif: Bin size (x,y) differs!" << endl; return; } if (binx != e4Dif->getBinX()) { cout << "Error in vnFourDif: Correlator x-binning differs!" << endl; return; } vector e2binPtrs; vector e4binPtrs; vector ref2Ptrs; vector ref4Ptrs; double denom = 2 * ref2.mean() * ref2.mean() - ref4.mean(); auto vn = [&] (int i) { // Test denominator. if (denom <= 0 ) return 0.; return ((2 * ref2.mean() * e2bins[i].mean() - e4bins[i].mean()) / pow(denom, 0.75)); }; // We need here a seperate error function, as we don't // iterate over the reference flow. auto vnerr = [&] (int i) { double denom2 = 2 * ref2Ptrs[i]->mean() * ref2Ptrs[i]->mean() - ref4Ptrs[i]->mean(); // Test denominator. if (denom2 <= 0) return 0.; return ((2 * ref2Ptrs[i]->mean() * e2binPtrs[i]->mean() - e4binPtrs[i]->mean()) / pow(denom2, 0.75)); }; // Error calculation. vector > yErr; ref2Ptrs = ref2.getBinPtrs(); ref4Ptrs = ref4.getBinPtrs(); for (int j = 0, N = e2bins.size(); j < N; ++j) { e2binPtrs = e2bins[j].getBinPtrs(); e4binPtrs = e4bins[j].getBinPtrs(); yErr.push_back(sampleError(vnerr)); } // Put the binPtrs back in place. e2binPtrs = e2Dif->getBinPtrs(); e4binPtrs = e4Dif->getBinPtrs(); fillScatter(h, binx, vn, yErr); } }; // End class CumulantAnalysis. } // End Rivet namespace. #endif