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/FinalState.hh"
#include "Rivet/Projections/ChargedFinalState.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);
}