Rivet analyses
Ratios (w.r.t. π+/π−) of hadron multiplicities in hadronic e+e− events
Experiment: PDG (Various)
Spires ID: 7857373
Status: VALIDATED
Authors: - Holger Schulz
References: - Phys. Lett. B, 667, 1 (2008)
Beams: e+ e-
Beam energies: (5.0, 5.0); (17.5, 17.5); (45.6, 45.6); (88.5, 88.5)GeV
Run details: - Hadronic events in e+e− collisions.
Ratios (w.r.t. π+/π−) of hadron multiplicities in hadronic e+e− events, taken from Review of Particle Properties 2008, table 40.1, page 355. Average hadron multiplicities per hadronic e+e− annihilation event at $\sqrt{s} \approx$ 10, 29–35, 91, and 130–200 GeV, normalised to the pion multiplicity. The numbers are averages from various experiments. Correlations of the systematic uncertainties were considered for the calculation of the averages.
Source
code:PDG_HADRON_MULTIPLICITIES_RATIOS.cc
// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/ChargedFinalState.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/UnstableParticles.hh"
namespace Rivet {
/// @brief Implementation of PDG hadron multiplicities as ratios to \f$ \pi^\pm \f$ multiplicity
/// @author Holger Schulz
class PDG_HADRON_MULTIPLICITIES_RATIOS : public Analysis {
public:
/// Constructor
RIVET_DEFAULT_ANALYSIS_CTOR(PDG_HADRON_MULTIPLICITIES_RATIOS);
/// @name Analysis methods
/// @{
void analyze(const Event& e) {
// First, veto on leptonic events by requiring at least 4 charged FS particles
const FinalState& fs = apply<FinalState>(e, "FS");
const size_t numParticles = fs.particles().size();
// Even if we only generate hadronic events, we still need a cut on numCharged >= 2.
if (numParticles < 2) {
MSG_DEBUG("Failed leptonic event cut");
vetoEvent;
}
MSG_DEBUG("Passed leptonic event cut");
MSG_DEBUG("sqrt(S) = " << sqrtS() / GeV << " GeV");
// Final state of unstable particles to get particle spectra
const UnstableParticles& ufs = apply<UnstableParticles>(e, "UFS");
if (sqrtS() / GeV >= 9.5 && sqrtS() / GeV <= 10.5) {
for (const Particle& p : ufs.particles()) {
const PdgId id = p.abspid();
switch (id) {
case 211: _weightedTotalNumPiPlus->fill(); break;
case 111: _histMeanMultiPi0->fill(_histMeanMultiPi0->bin(1).xMid()); break;
case 321: _histMeanMultiKPlus->fill(_histMeanMultiKPlus->bin(1).xMid()); break;
case 130:
case 310: _histMeanMultiK0->fill(_histMeanMultiK0->bin(1).xMid()); break;
case 221: _histMeanMultiEta->fill(_histMeanMultiEta->bin(1).xMid()); break;
case 331: _histMeanMultiEtaPrime->fill(_histMeanMultiEtaPrime->bin(1).xMid()); break;
case 411: _histMeanMultiDPlus->fill(_histMeanMultiDPlus->bin(1).xMid()); break;
case 421: _histMeanMultiD0->fill(_histMeanMultiD0->bin(1).xMid()); break;
case 431: _histMeanMultiDPlus_s->fill(_histMeanMultiDPlus_s->bin(1).xMid()); break;
case 9010221: _histMeanMultiF0_980->fill(_histMeanMultiF0_980->bin(1).xMid()); break;
case 113: _histMeanMultiRho770_0->fill(_histMeanMultiRho770_0->bin(1).xMid()); break;
case 223: _histMeanMultiOmega782->fill(_histMeanMultiOmega782->bin(1).xMid()); break;
case 323: _histMeanMultiKStar892Plus->fill(_histMeanMultiKStar892Plus->bin(1).xMid()); break;
case 313: _histMeanMultiKStar892_0->fill(_histMeanMultiKStar892_0->bin(1).xMid()); break;
case 333: _histMeanMultiPhi1020->fill(_histMeanMultiPhi1020->bin(1).xMid()); break;
case 413: _histMeanMultiDStar2010Plus->fill(_histMeanMultiDStar2010Plus->bin(1).xMid()); break;
case 423: _histMeanMultiDStar2007_0->fill(_histMeanMultiDStar2007_0->bin(1).xMid()); break;
case 433:
_histMeanMultiDStar_s2112Plus->fill(_histMeanMultiDStar_s2112Plus->bin(1).xMid());
break;
case 443: _histMeanMultiJPsi1S->fill(_histMeanMultiJPsi1S->bin(1).xMid()); break;
case 225: _histMeanMultiF2_1270->fill(_histMeanMultiF2_1270->bin(1).xMid()); break;
case 2212: _histMeanMultiP->fill(_histMeanMultiP->bin(1).xMid()); break;
case 3122: _histMeanMultiLambda->fill(_histMeanMultiLambda->bin(1).xMid()); break;
case 3212: _histMeanMultiSigma0->fill(_histMeanMultiSigma0->bin(1).xMid()); break;
case 3312: _histMeanMultiXiMinus->fill(_histMeanMultiXiMinus->bin(1).xMid()); break;
case 2224:
_histMeanMultiDelta1232PlusPlus->fill(_histMeanMultiDelta1232PlusPlus->bin(1).xMid());
break;
case 3114:
_histMeanMultiSigma1385Minus->fill(_histMeanMultiSigma1385Minus->bin(1).xMid());
_histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->bin(1).xMid());
break;
case 3224:
_histMeanMultiSigma1385Plus->fill(_histMeanMultiSigma1385Plus->bin(1).xMid());
_histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->bin(1).xMid());
break;
case 3324: _histMeanMultiXi1530_0->fill(_histMeanMultiXi1530_0->bin(1).xMid()); break;
case 3334: _histMeanMultiOmegaMinus->fill(_histMeanMultiOmegaMinus->bin(1).xMid()); break;
case 4122: _histMeanMultiLambda_c_Plus->fill(_histMeanMultiLambda_c_Plus->bin(1).xMid()); break;
case 4222:
case 4112:
_histMeanMultiSigma_c_PlusPlus_0->fill(_histMeanMultiSigma_c_PlusPlus_0->bin(1).xMid());
break;
case 102134: _histMeanMultiLambda1520->fill(_histMeanMultiLambda1520->bin(1).xMid()); break;
}
}
}
if (sqrtS() / GeV >= 29 && sqrtS() / GeV <= 35) {
for (const Particle& p : ufs.particles()) {
const PdgId id = p.abspid();
switch (id) {
case 211: _weightedTotalNumPiPlus->fill(); break;
case 111: _histMeanMultiPi0->fill(_histMeanMultiPi0->bin(1).xMid()); break;
case 321: _histMeanMultiKPlus->fill(_histMeanMultiKPlus->bin(1).xMid()); break;
case 130:
case 310: _histMeanMultiK0->fill(_histMeanMultiK0->bin(1).xMid()); break;
case 221: _histMeanMultiEta->fill(_histMeanMultiEta->bin(1).xMid()); break;
case 331: _histMeanMultiEtaPrime->fill(_histMeanMultiEtaPrime->bin(1).xMid()); break;
case 411: _histMeanMultiDPlus->fill(_histMeanMultiDPlus->bin(1).xMid()); break;
case 421: _histMeanMultiD0->fill(_histMeanMultiD0->bin(1).xMid()); break;
case 431: _histMeanMultiDPlus_s->fill(_histMeanMultiDPlus_s->bin(1).xMid()); break;
case 9010221: _histMeanMultiF0_980->fill(_histMeanMultiF0_980->bin(1).xMid()); break;
case 113: _histMeanMultiRho770_0->fill(_histMeanMultiRho770_0->bin(1).xMid()); break;
case 323: _histMeanMultiKStar892Plus->fill(_histMeanMultiKStar892Plus->bin(1).xMid()); break;
case 313: _histMeanMultiKStar892_0->fill(_histMeanMultiKStar892_0->bin(1).xMid()); break;
case 333: _histMeanMultiPhi1020->fill(_histMeanMultiPhi1020->bin(1).xMid()); break;
case 413: _histMeanMultiDStar2010Plus->fill(_histMeanMultiDStar2010Plus->bin(1).xMid()); break;
case 423: _histMeanMultiDStar2007_0->fill(_histMeanMultiDStar2007_0->bin(1).xMid()); break;
case 225: _histMeanMultiF2_1270->fill(_histMeanMultiF2_1270->bin(1).xMid()); break;
case 325: _histMeanMultiK2Star1430Plus->fill(_histMeanMultiK2Star1430Plus->bin(1).xMid()); break;
case 315: _histMeanMultiK2Star1430_0->fill(_histMeanMultiK2Star1430_0->bin(1).xMid()); break;
case 2212: _histMeanMultiP->fill(_histMeanMultiP->bin(1).xMid()); break;
case 3122: _histMeanMultiLambda->fill(_histMeanMultiLambda->bin(1).xMid()); break;
case 3312: _histMeanMultiXiMinus->fill(_histMeanMultiXiMinus->bin(1).xMid()); break;
case 3114:
_histMeanMultiSigma1385Minus->fill(_histMeanMultiSigma1385Minus->bin(1).xMid());
_histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->bin(1).xMid());
break;
case 3224:
_histMeanMultiSigma1385Plus->fill(_histMeanMultiSigma1385Plus->bin(1).xMid());
_histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->bin(1).xMid());
break;
case 3334: _histMeanMultiOmegaMinus->fill(_histMeanMultiOmegaMinus->bin(1).xMid()); break;
case 4122: _histMeanMultiLambda_c_Plus->fill(_histMeanMultiLambda_c_Plus->bin(1).xMid()); break;
}
}
}
if (sqrtS() / GeV >= 89.5 && sqrtS() / GeV <= 91.8) {
for (const Particle& p : ufs.particles()) {
const PdgId id = p.abspid();
switch (id) {
case 211: _weightedTotalNumPiPlus->fill(); break;
case 111: _histMeanMultiPi0->fill(_histMeanMultiPi0->bin(1).xMid()); break;
case 321: _histMeanMultiKPlus->fill(_histMeanMultiKPlus->bin(1).xMid()); break;
case 130:
case 310: _histMeanMultiK0->fill(_histMeanMultiK0->bin(1).xMid()); break;
case 221: _histMeanMultiEta->fill(_histMeanMultiEta->bin(1).xMid()); break;
case 331: _histMeanMultiEtaPrime->fill(_histMeanMultiEtaPrime->bin(1).xMid()); break;
case 411: _histMeanMultiDPlus->fill(_histMeanMultiDPlus->bin(1).xMid()); break;
case 421: _histMeanMultiD0->fill(_histMeanMultiD0->bin(1).xMid()); break;
case 431: _histMeanMultiDPlus_s->fill(_histMeanMultiDPlus_s->bin(1).xMid()); break;
case 511: _histMeanMultiBPlus_B0_d->fill(_histMeanMultiBPlus_B0_d->bin(1).xMid()); break;
case 521:
_histMeanMultiBPlus_B0_d->fill(_histMeanMultiBPlus_B0_d->bin(1).xMid());
_histMeanMultiBPlus_u->fill(_histMeanMultiBPlus_u->bin(1).xMid());
break;
case 531: _histMeanMultiB0_s->fill(_histMeanMultiB0_s->bin(1).xMid()); break;
case 9010221: _histMeanMultiF0_980->fill(_histMeanMultiF0_980->bin(1).xMid()); break;
case 9000211: _histMeanMultiA0_980Plus->fill(_histMeanMultiA0_980Plus->bin(1).xMid()); break;
case 113: _histMeanMultiRho770_0->fill(_histMeanMultiRho770_0->bin(1).xMid()); break;
case 213: _histMeanMultiRho770Plus->fill(_histMeanMultiRho770Plus->bin(1).xMid()); break;
case 223: _histMeanMultiOmega782->fill(_histMeanMultiOmega782->bin(1).xMid()); break;
case 323: _histMeanMultiKStar892Plus->fill(_histMeanMultiKStar892Plus->bin(1).xMid()); break;
case 313: _histMeanMultiKStar892_0->fill(_histMeanMultiKStar892_0->bin(1).xMid()); break;
case 333: _histMeanMultiPhi1020->fill(_histMeanMultiPhi1020->bin(1).xMid()); break;
case 413: _histMeanMultiDStar2010Plus->fill(_histMeanMultiDStar2010Plus->bin(1).xMid()); break;
case 433:
_histMeanMultiDStar_s2112Plus->fill(_histMeanMultiDStar_s2112Plus->bin(1).xMid());
break;
case 513:
case 523:
case 533: _histMeanMultiBStar->fill(_histMeanMultiBStar->bin(1).xMid()); break;
case 443: _histMeanMultiJPsi1S->fill(_histMeanMultiJPsi1S->bin(1).xMid()); break;
case 100443: _histMeanMultiPsi2S->fill(_histMeanMultiPsi2S->bin(1).xMid()); break;
case 553: _histMeanMultiUpsilon1S->fill(_histMeanMultiUpsilon1S->bin(1).xMid()); break;
case 20223: _histMeanMultiF1_1285->fill(_histMeanMultiF1_1285->bin(1).xMid()); break;
case 20333: _histMeanMultiF1_1420->fill(_histMeanMultiF1_1420->bin(1).xMid()); break;
case 445: _histMeanMultiChi_c1_3510->fill(_histMeanMultiChi_c1_3510->bin(1).xMid()); break;
case 225: _histMeanMultiF2_1270->fill(_histMeanMultiF2_1270->bin(1).xMid()); break;
case 335: _histMeanMultiF2Prime1525->fill(_histMeanMultiF2Prime1525->bin(1).xMid()); break;
case 315: _histMeanMultiK2Star1430_0->fill(_histMeanMultiK2Star1430_0->bin(1).xMid()); break;
case 515:
case 525:
case 535: _histMeanMultiBStarStar->fill(_histMeanMultiBStarStar->bin(1).xMid()); break;
case 10433:
case 20433: _histMeanMultiDs1Plus->fill(_histMeanMultiDs1Plus->bin(1).xMid()); break;
case 435: _histMeanMultiDs2Plus->fill(_histMeanMultiDs2Plus->bin(1).xMid()); break;
case 2212: _histMeanMultiP->fill(_histMeanMultiP->bin(1).xMid()); break;
case 3122: _histMeanMultiLambda->fill(_histMeanMultiLambda->bin(1).xMid()); break;
case 3212: _histMeanMultiSigma0->fill(_histMeanMultiSigma0->bin(1).xMid()); break;
case 3112:
_histMeanMultiSigmaMinus->fill(_histMeanMultiSigmaMinus->bin(1).xMid());
_histMeanMultiSigmaPlusMinus->fill(_histMeanMultiSigmaPlusMinus->bin(1).xMid());
break;
case 3222:
_histMeanMultiSigmaPlus->fill(_histMeanMultiSigmaPlus->bin(1).xMid());
_histMeanMultiSigmaPlusMinus->fill(_histMeanMultiSigmaPlusMinus->bin(1).xMid());
break;
case 3312: _histMeanMultiXiMinus->fill(_histMeanMultiXiMinus->bin(1).xMid()); break;
case 2224:
_histMeanMultiDelta1232PlusPlus->fill(_histMeanMultiDelta1232PlusPlus->bin(1).xMid());
break;
case 3114:
_histMeanMultiSigma1385Minus->fill(_histMeanMultiSigma1385Minus->bin(1).xMid());
_histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->bin(1).xMid());
break;
case 3224:
_histMeanMultiSigma1385Plus->fill(_histMeanMultiSigma1385Plus->bin(1).xMid());
_histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->bin(1).xMid());
break;
case 3324: _histMeanMultiXi1530_0->fill(_histMeanMultiXi1530_0->bin(1).xMid()); break;
case 3334: _histMeanMultiOmegaMinus->fill(_histMeanMultiOmegaMinus->bin(1).xMid()); break;
case 4122: _histMeanMultiLambda_c_Plus->fill(_histMeanMultiLambda_c_Plus->bin(1).xMid()); break;
case 5122: _histMeanMultiLambda_b_0->fill(_histMeanMultiLambda_b_0->bin(1).xMid()); break;
case 102134: _histMeanMultiLambda1520->fill(_histMeanMultiLambda1520->bin(1).xMid()); break;
}
}
}
if (sqrtS() / GeV >= 130 && sqrtS() / GeV <= 200) {
for (const Particle& p : ufs.particles()) {
const PdgId id = p.abspid();
switch (id) {
case 211: _weightedTotalNumPiPlus->fill(); break;
case 321: _histMeanMultiKPlus->fill(_histMeanMultiKPlus->bin(1).xMid()); break;
case 130:
case 310: _histMeanMultiK0->fill(_histMeanMultiK0->bin(1).xMid()); break;
case 2212: _histMeanMultiP->fill(_histMeanMultiP->bin(1).xMid()); break;
case 3122: _histMeanMultiLambda->fill(_histMeanMultiLambda->bin(1).xMid()); break;
}
}
}
}
void init() {
declare(ChargedFinalState(), "FS");
declare(UnstableParticles(), "UFS");
book(_weightedTotalNumPiPlus, "TMP/PiPlus");
if (sqrtS() / GeV >= 9.5 && sqrtS() / GeV <= 10.5) {
book(_histMeanMultiPi0, 2, 1, 1);
book(_histMeanMultiKPlus, 3, 1, 1);
book(_histMeanMultiK0, 4, 1, 1);
book(_histMeanMultiEta, 5, 1, 1);
book(_histMeanMultiEtaPrime, 6, 1, 1);
book(_histMeanMultiDPlus, 7, 1, 1);
book(_histMeanMultiD0, 8, 1, 1);
book(_histMeanMultiDPlus_s, 9, 1, 1);
book(_histMeanMultiF0_980, 13, 1, 1);
book(_histMeanMultiRho770_0, 15, 1, 1);
book(_histMeanMultiOmega782, 17, 1, 1);
book(_histMeanMultiKStar892Plus, 18, 1, 1);
book(_histMeanMultiKStar892_0, 19, 1, 1);
book(_histMeanMultiPhi1020, 20, 1, 1);
book(_histMeanMultiDStar2010Plus, 21, 1, 1);
book(_histMeanMultiDStar2007_0, 22, 1, 1);
book(_histMeanMultiDStar_s2112Plus, 23, 1, 1);
book(_histMeanMultiJPsi1S, 25, 1, 1);
book(_histMeanMultiF2_1270, 31, 1, 1);
book(_histMeanMultiP, 38, 1, 1);
book(_histMeanMultiLambda, 39, 1, 1);
book(_histMeanMultiSigma0, 40, 1, 1);
book(_histMeanMultiXiMinus, 44, 1, 1);
book(_histMeanMultiDelta1232PlusPlus, 45, 1, 1);
book(_histMeanMultiSigma1385Minus, 46, 1, 1);
book(_histMeanMultiSigma1385Plus, 47, 1, 1);
book(_histMeanMultiSigma1385PlusMinus, 48, 1, 1);
book(_histMeanMultiXi1530_0, 49, 1, 1);
book(_histMeanMultiOmegaMinus, 50, 1, 1);
book(_histMeanMultiLambda_c_Plus, 51, 1, 1);
book(_histMeanMultiSigma_c_PlusPlus_0, 53, 1, 1);
book(_histMeanMultiLambda1520, 54, 1, 1);
}
if (sqrtS() / GeV >= 29 && sqrtS() / GeV <= 35) {
book(_histMeanMultiPi0, 2, 1, 2);
book(_histMeanMultiKPlus, 3, 1, 2);
book(_histMeanMultiK0, 4, 1, 2);
book(_histMeanMultiEta, 5, 1, 2);
book(_histMeanMultiEtaPrime, 6, 1, 2);
book(_histMeanMultiDPlus, 7, 1, 2);
book(_histMeanMultiD0, 8, 1, 2);
book(_histMeanMultiDPlus_s, 9, 1, 2);
book(_histMeanMultiF0_980, 13, 1, 2);
book(_histMeanMultiRho770_0, 15, 1, 2);
book(_histMeanMultiKStar892Plus, 18, 1, 2);
book(_histMeanMultiKStar892_0, 19, 1, 2);
book(_histMeanMultiPhi1020, 20, 1, 2);
book(_histMeanMultiDStar2010Plus, 21, 1, 2);
book(_histMeanMultiDStar2007_0, 22, 1, 2);
book(_histMeanMultiF2_1270, 31, 1, 2);
book(_histMeanMultiK2Star1430Plus, 33, 1, 1);
book(_histMeanMultiK2Star1430_0, 34, 1, 1);
book(_histMeanMultiP, 38, 1, 2);
book(_histMeanMultiLambda, 39, 1, 2);
book(_histMeanMultiXiMinus, 44, 1, 2);
book(_histMeanMultiSigma1385Minus, 46, 1, 2);
book(_histMeanMultiSigma1385Plus, 47, 1, 2);
book(_histMeanMultiSigma1385PlusMinus, 48, 1, 2);
book(_histMeanMultiOmegaMinus, 50, 1, 2);
book(_histMeanMultiLambda_c_Plus, 51, 1, 2);
}
if (sqrtS() / GeV >= 89.5 && sqrtS() / GeV <= 91.8) {
book(_histMeanMultiPi0, 2, 1, 3);
book(_histMeanMultiKPlus, 3, 1, 3);
book(_histMeanMultiK0, 4, 1, 3);
book(_histMeanMultiEta, 5, 1, 3);
book(_histMeanMultiEtaPrime, 6, 1, 3);
book(_histMeanMultiDPlus, 7, 1, 3);
book(_histMeanMultiD0, 8, 1, 3);
book(_histMeanMultiDPlus_s, 9, 1, 3);
book(_histMeanMultiBPlus_B0_d, 10, 1, 1);
book(_histMeanMultiBPlus_u, 11, 1, 1);
book(_histMeanMultiB0_s, 12, 1, 1);
book(_histMeanMultiF0_980, 13, 1, 3);
book(_histMeanMultiA0_980Plus, 14, 1, 1);
book(_histMeanMultiRho770_0, 15, 1, 3);
book(_histMeanMultiRho770Plus, 16, 1, 1);
book(_histMeanMultiOmega782, 17, 1, 2);
book(_histMeanMultiKStar892Plus, 18, 1, 3);
book(_histMeanMultiKStar892_0, 19, 1, 3);
book(_histMeanMultiPhi1020, 20, 1, 3);
book(_histMeanMultiDStar2010Plus, 21, 1, 3);
book(_histMeanMultiDStar_s2112Plus, 23, 1, 2);
book(_histMeanMultiBStar, 24, 1, 1);
book(_histMeanMultiJPsi1S, 25, 1, 2);
book(_histMeanMultiPsi2S, 26, 1, 1);
book(_histMeanMultiUpsilon1S, 27, 1, 1);
book(_histMeanMultiF1_1285, 28, 1, 1);
book(_histMeanMultiF1_1420, 29, 1, 1);
book(_histMeanMultiChi_c1_3510, 30, 1, 1);
book(_histMeanMultiF2_1270, 31, 1, 3);
book(_histMeanMultiF2Prime1525, 32, 1, 1);
book(_histMeanMultiK2Star1430_0, 34, 1, 2);
book(_histMeanMultiBStarStar, 35, 1, 1);
book(_histMeanMultiDs1Plus, 36, 1, 1);
book(_histMeanMultiDs2Plus, 37, 1, 1);
book(_histMeanMultiP, 38, 1, 3);
book(_histMeanMultiLambda, 39, 1, 3);
book(_histMeanMultiSigma0, 40, 1, 2);
book(_histMeanMultiSigmaMinus, 41, 1, 1);
book(_histMeanMultiSigmaPlus, 42, 1, 1);
book(_histMeanMultiSigmaPlusMinus, 43, 1, 1);
book(_histMeanMultiXiMinus, 44, 1, 3);
book(_histMeanMultiDelta1232PlusPlus, 45, 1, 2);
book(_histMeanMultiSigma1385Minus, 46, 1, 3);
book(_histMeanMultiSigma1385Plus, 47, 1, 3);
book(_histMeanMultiSigma1385PlusMinus, 48, 1, 3);
book(_histMeanMultiXi1530_0, 49, 1, 2);
book(_histMeanMultiOmegaMinus, 50, 1, 3);
book(_histMeanMultiLambda_c_Plus, 51, 1, 3);
book(_histMeanMultiLambda_b_0, 52, 1, 1);
book(_histMeanMultiLambda1520, 54, 1, 2);
}
if (sqrtS() / GeV >= 130 && sqrtS() / GeV <= 200) {
book(_histMeanMultiKPlus, 3, 1, 4);
book(_histMeanMultiK0, 4, 1, 4);
book(_histMeanMultiP, 38, 1, 4);
book(_histMeanMultiLambda, 39, 1, 4);
}
}
// Finalize
void finalize() {
if (sqrtS() / GeV >= 9.5 && sqrtS() / GeV <= 10.5) {
scale(_histMeanMultiPi0, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiKPlus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiK0, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiEta, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiEtaPrime, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiDPlus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiD0, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiDPlus_s, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiF0_980, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiRho770_0, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiOmega782, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiKStar892Plus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiKStar892_0, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiPhi1020, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiDStar2010Plus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiDStar2007_0, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiDStar_s2112Plus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiJPsi1S, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiF2_1270, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiP, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiLambda, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiSigma0, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiXiMinus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiDelta1232PlusPlus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiSigma1385Minus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiSigma1385Plus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiSigma1385PlusMinus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiXi1530_0, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiOmegaMinus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiLambda_c_Plus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiSigma_c_PlusPlus_0, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiLambda1520, 1.0 / *_weightedTotalNumPiPlus);
}
if (sqrtS() / GeV >= 29 && sqrtS() / GeV <= 35) {
scale(_histMeanMultiPi0, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiKPlus, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiK0, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiEta, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiEtaPrime, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiDPlus, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiD0, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiDPlus_s, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiF0_980, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiRho770_0, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiKStar892Plus, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiKStar892_0, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiPhi1020, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiDStar2010Plus, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiDStar2007_0, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiF2_1270, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiK2Star1430Plus, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiK2Star1430_0, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiP, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiLambda, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiXiMinus, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiSigma1385Minus, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiSigma1385Plus, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiSigma1385PlusMinus, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiOmegaMinus, 5.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiLambda_c_Plus, 5.0 / *_weightedTotalNumPiPlus);
}
if (sqrtS() / GeV >= 89.5 && sqrtS() / GeV <= 91.8) {
scale(_histMeanMultiPi0, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiKPlus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiK0, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiEta, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiEtaPrime, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiDPlus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiD0, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiDPlus_s, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiBPlus_B0_d, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiBPlus_u, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiB0_s, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiF0_980, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiA0_980Plus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiRho770_0, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiRho770Plus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiOmega782, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiKStar892Plus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiKStar892_0, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiPhi1020, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiDStar2010Plus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiDStar_s2112Plus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiBStar, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiJPsi1S, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiPsi2S, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiUpsilon1S, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiF1_1285, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiF1_1420, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiChi_c1_3510, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiF2_1270, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiF2Prime1525, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiK2Star1430_0, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiBStarStar, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiDs1Plus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiDs2Plus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiP, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiLambda, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiSigma0, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiSigmaMinus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiSigmaPlus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiSigmaPlusMinus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiXiMinus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiDelta1232PlusPlus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiSigma1385Minus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiSigma1385Plus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiSigma1385PlusMinus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiXi1530_0, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiOmegaMinus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiLambda_c_Plus, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiLambda_b_0, 1.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiLambda1520, 1.0 / *_weightedTotalNumPiPlus);
}
if (sqrtS() / GeV >= 130 && sqrtS() / GeV <= 200) {
scale(_histMeanMultiKPlus, 70.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiK0, 70.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiP, 70.0 / *_weightedTotalNumPiPlus);
scale(_histMeanMultiLambda, 70.0 / *_weightedTotalNumPiPlus);
}
}
/// @}
private:
CounterPtr _weightedTotalNumPiPlus;
Histo1DPtr _histMeanMultiPi0;
Histo1DPtr _histMeanMultiKPlus;
Histo1DPtr _histMeanMultiK0;
Histo1DPtr _histMeanMultiEta;
Histo1DPtr _histMeanMultiEtaPrime;
Histo1DPtr _histMeanMultiDPlus;
Histo1DPtr _histMeanMultiD0;
Histo1DPtr _histMeanMultiDPlus_s;
Histo1DPtr _histMeanMultiBPlus_B0_d;
Histo1DPtr _histMeanMultiBPlus_u;
Histo1DPtr _histMeanMultiB0_s;
Histo1DPtr _histMeanMultiF0_980;
Histo1DPtr _histMeanMultiA0_980Plus;
Histo1DPtr _histMeanMultiRho770_0;
Histo1DPtr _histMeanMultiRho770Plus;
Histo1DPtr _histMeanMultiOmega782;
Histo1DPtr _histMeanMultiKStar892Plus;
Histo1DPtr _histMeanMultiKStar892_0;
Histo1DPtr _histMeanMultiPhi1020;
Histo1DPtr _histMeanMultiDStar2010Plus;
Histo1DPtr _histMeanMultiDStar2007_0;
Histo1DPtr _histMeanMultiDStar_s2112Plus;
Histo1DPtr _histMeanMultiBStar;
Histo1DPtr _histMeanMultiJPsi1S;
Histo1DPtr _histMeanMultiPsi2S;
Histo1DPtr _histMeanMultiUpsilon1S;
Histo1DPtr _histMeanMultiF1_1285;
Histo1DPtr _histMeanMultiF1_1420;
Histo1DPtr _histMeanMultiChi_c1_3510;
Histo1DPtr _histMeanMultiF2_1270;
Histo1DPtr _histMeanMultiF2Prime1525;
Histo1DPtr _histMeanMultiK2Star1430Plus;
Histo1DPtr _histMeanMultiK2Star1430_0;
Histo1DPtr _histMeanMultiBStarStar;
Histo1DPtr _histMeanMultiDs1Plus;
Histo1DPtr _histMeanMultiDs2Plus;
Histo1DPtr _histMeanMultiP;
Histo1DPtr _histMeanMultiLambda;
Histo1DPtr _histMeanMultiSigma0;
Histo1DPtr _histMeanMultiSigmaMinus;
Histo1DPtr _histMeanMultiSigmaPlus;
Histo1DPtr _histMeanMultiSigmaPlusMinus;
Histo1DPtr _histMeanMultiXiMinus;
Histo1DPtr _histMeanMultiDelta1232PlusPlus;
Histo1DPtr _histMeanMultiSigma1385Minus;
Histo1DPtr _histMeanMultiSigma1385Plus;
Histo1DPtr _histMeanMultiSigma1385PlusMinus;
Histo1DPtr _histMeanMultiXi1530_0;
Histo1DPtr _histMeanMultiOmegaMinus;
Histo1DPtr _histMeanMultiLambda_c_Plus;
Histo1DPtr _histMeanMultiLambda_b_0;
Histo1DPtr _histMeanMultiSigma_c_PlusPlus_0;
Histo1DPtr _histMeanMultiLambda1520;
/// @}
};
RIVET_DECLARE_PLUGIN(PDG_HADRON_MULTIPLICITIES_RATIOS);
}