|
[Rivet-svn] r2246 - trunk/src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgTue Feb 9 11:18:02 GMT 2010
Author: hoeth Date: Tue Feb 9 11:18:01 2010 New Revision: 2246 Log: e+e- multiplicity analyses: Only book the histograms for the correct energies. Modified: trunk/src/Analyses/PDG_Hadron_Multiplicities.cc trunk/src/Analyses/PDG_Hadron_Multiplicities_Ratios.cc Modified: trunk/src/Analyses/PDG_Hadron_Multiplicities.cc ============================================================================== --- trunk/src/Analyses/PDG_Hadron_Multiplicities.cc Mon Feb 8 18:41:49 2010 (r2245) +++ trunk/src/Analyses/PDG_Hadron_Multiplicities.cc Tue Feb 9 11:18:01 2010 (r2246) @@ -2,7 +2,6 @@ #include "Rivet/Analysis.hh" #include "Rivet/RivetAIDA.hh" #include "Rivet/Tools/ParticleIdUtils.hh" -#include "Rivet/Projections/Beam.hh" #include "Rivet/Projections/FinalState.hh" #include "Rivet/Projections/ChargedFinalState.hh" #include "Rivet/Projections/UnstableFinalState.hh" @@ -14,7 +13,7 @@ /// @author Hendrik Hoeth class PDG_HADRON_MULTIPLICITIES : public Analysis { public: - + /// Constructor PDG_HADRON_MULTIPLICITIES() : Analysis("PDG_HADRON_MULTIPLICITIES") { @@ -29,7 +28,7 @@ // First, veto on leptonic events by requiring at least 4 charged FS particles const FinalState& fs = applyProjection<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) { getLog() << Log::DEBUG << "Failed leptonic event cut" << endl; @@ -40,788 +39,727 @@ // Get event weight for histo filling const double weight = e.weight(); - // Get beams and average beam momentum - const ParticlePair& beams = applyProjection<Beam>(e, "Beams").beams(); - const double meanBeamMom = ( beams.first.momentum().vector3().mod() + - beams.second.momentum().vector3().mod() ) / 2.0; - getLog() << Log::DEBUG << "Avg beam momentum = " << meanBeamMom/GeV << " GeV" << endl; + getLog() << Log::DEBUG << "sqrt(S) = " << sqrtS()/GeV << " GeV" << endl; // Final state of unstable particles to get particle spectra const UnstableFinalState& ufs = applyProjection<UnstableFinalState>(e, "UFS"); - - if (2*meanBeamMom/GeV >= 9.5 && 2*meanBeamMom/GeV <= 10.5) { + if (sqrtS()/GeV >= 9.5 && sqrtS()/GeV <= 10.5) { foreach (const Particle& p, ufs.particles()) { const PdgId id = abs(p.pdgId()); switch (id) { case 211: - _hist10MeanMultiPiPlus->fill(_hist10MeanMultiPiPlus->binMean(0), weight); + _histMeanMultiPiPlus->fill(_histMeanMultiPiPlus->binMean(0), weight); break; case 111: - _hist10MeanMultiPi0->fill(_hist10MeanMultiPi0->binMean(0), weight); + _histMeanMultiPi0->fill(_histMeanMultiPi0->binMean(0), weight); break; case 321: - _hist10MeanMultiKPlus->fill(_hist10MeanMultiKPlus->binMean(0), weight); + _histMeanMultiKPlus->fill(_histMeanMultiKPlus->binMean(0), weight); break; case 130: case 310: - _hist10MeanMultiK0->fill(_hist10MeanMultiK0->binMean(0), weight); + _histMeanMultiK0->fill(_histMeanMultiK0->binMean(0), weight); break; case 221: - _hist10MeanMultiEta->fill(_hist10MeanMultiEta->binMean(0), weight); + _histMeanMultiEta->fill(_histMeanMultiEta->binMean(0), weight); break; case 331: - _hist10MeanMultiEtaPrime->fill(_hist10MeanMultiEtaPrime->binMean(0), weight); + _histMeanMultiEtaPrime->fill(_histMeanMultiEtaPrime->binMean(0), weight); break; case 411: - _hist10MeanMultiDPlus->fill(_hist10MeanMultiDPlus->binMean(0), weight); + _histMeanMultiDPlus->fill(_histMeanMultiDPlus->binMean(0), weight); break; case 421: - _hist10MeanMultiD0->fill(_hist10MeanMultiD0->binMean(0), weight); + _histMeanMultiD0->fill(_histMeanMultiD0->binMean(0), weight); break; case 431: - _hist10MeanMultiDPlus_s->fill(_hist10MeanMultiDPlus_s->binMean(0), weight); + _histMeanMultiDPlus_s->fill(_histMeanMultiDPlus_s->binMean(0), weight); break; case 9010221: - _hist10MeanMultiF0_980->fill(_hist10MeanMultiF0_980->binMean(0), weight); + _histMeanMultiF0_980->fill(_histMeanMultiF0_980->binMean(0), weight); break; case 113: - _hist10MeanMultiRho770_0->fill(_hist10MeanMultiRho770_0->binMean(0), weight); + _histMeanMultiRho770_0->fill(_histMeanMultiRho770_0->binMean(0), weight); break; case 223: - _hist10MeanMultiOmega782->fill(_hist10MeanMultiOmega782->binMean(0), weight); + _histMeanMultiOmega782->fill(_histMeanMultiOmega782->binMean(0), weight); break; case 323: - _hist10MeanMultiKStar892Plus->fill(_hist10MeanMultiKStar892Plus->binMean(0), weight); + _histMeanMultiKStar892Plus->fill(_histMeanMultiKStar892Plus->binMean(0), weight); break; case 313: - _hist10MeanMultiKStar892_0->fill(_hist10MeanMultiKStar892_0->binMean(0), weight); + _histMeanMultiKStar892_0->fill(_histMeanMultiKStar892_0->binMean(0), weight); break; case 333: - _hist10MeanMultiPhi1020->fill(_hist10MeanMultiPhi1020->binMean(0), weight); + _histMeanMultiPhi1020->fill(_histMeanMultiPhi1020->binMean(0), weight); break; case 413: - _hist10MeanMultiDStar2010Plus->fill(_hist10MeanMultiDStar2010Plus->binMean(0), weight); + _histMeanMultiDStar2010Plus->fill(_histMeanMultiDStar2010Plus->binMean(0), weight); break; case 423: - _hist10MeanMultiDStar2007_0->fill(_hist10MeanMultiDStar2007_0->binMean(0), weight); + _histMeanMultiDStar2007_0->fill(_histMeanMultiDStar2007_0->binMean(0), weight); break; case 433: - _hist10MeanMultiDStar_s2112Plus->fill(_hist10MeanMultiDStar_s2112Plus->binMean(0), weight); + _histMeanMultiDStar_s2112Plus->fill(_histMeanMultiDStar_s2112Plus->binMean(0), weight); break; case 443: - _hist10MeanMultiJPsi1S->fill(_hist10MeanMultiJPsi1S->binMean(0), weight); + _histMeanMultiJPsi1S->fill(_histMeanMultiJPsi1S->binMean(0), weight); break; case 225: - _hist10MeanMultiF2_1270->fill(_hist10MeanMultiF2_1270->binMean(0), weight); + _histMeanMultiF2_1270->fill(_histMeanMultiF2_1270->binMean(0), weight); break; case 2212: - _hist10MeanMultiP->fill(_hist10MeanMultiP->binMean(0), weight); + _histMeanMultiP->fill(_histMeanMultiP->binMean(0), weight); break; case 3122: - _hist10MeanMultiLambda->fill(_hist10MeanMultiLambda->binMean(0), weight); + _histMeanMultiLambda->fill(_histMeanMultiLambda->binMean(0), weight); break; case 3212: - _hist10MeanMultiSigma0->fill(_hist10MeanMultiSigma0->binMean(0), weight); + _histMeanMultiSigma0->fill(_histMeanMultiSigma0->binMean(0), weight); break; case 3312: - _hist10MeanMultiXiMinus->fill(_hist10MeanMultiXiMinus->binMean(0), weight); + _histMeanMultiXiMinus->fill(_histMeanMultiXiMinus->binMean(0), weight); break; case 2224: - _hist10MeanMultiDelta1232PlusPlus->fill(_hist10MeanMultiDelta1232PlusPlus->binMean(0), weight); + _histMeanMultiDelta1232PlusPlus->fill(_histMeanMultiDelta1232PlusPlus->binMean(0), weight); break; case 3114: - _hist10MeanMultiSigma1385Minus->fill(_hist10MeanMultiSigma1385Minus->binMean(0), weight); - _hist10MeanMultiSigma1385PlusMinus->fill(_hist10MeanMultiSigma1385PlusMinus->binMean(0), weight); + _histMeanMultiSigma1385Minus->fill(_histMeanMultiSigma1385Minus->binMean(0), weight); + _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->binMean(0), weight); break; case 3224: - _hist10MeanMultiSigma1385Plus->fill(_hist10MeanMultiSigma1385Plus->binMean(0), weight); - _hist10MeanMultiSigma1385PlusMinus->fill(_hist10MeanMultiSigma1385PlusMinus->binMean(0), weight); + _histMeanMultiSigma1385Plus->fill(_histMeanMultiSigma1385Plus->binMean(0), weight); + _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->binMean(0), weight); break; case 3324: - _hist10MeanMultiXi1530_0->fill(_hist10MeanMultiXi1530_0->binMean(0), weight); + _histMeanMultiXi1530_0->fill(_histMeanMultiXi1530_0->binMean(0), weight); break; case 3334: - _hist10MeanMultiOmegaMinus->fill(_hist10MeanMultiOmegaMinus->binMean(0), weight); + _histMeanMultiOmegaMinus->fill(_histMeanMultiOmegaMinus->binMean(0), weight); break; case 4122: - _hist10MeanMultiLambda_c_Plus->fill(_hist10MeanMultiLambda_c_Plus->binMean(0), weight); + _histMeanMultiLambda_c_Plus->fill(_histMeanMultiLambda_c_Plus->binMean(0), weight); break; case 4222: case 4112: - _hist10MeanMultiSigma_c_PlusPlus_0->fill(_hist10MeanMultiSigma_c_PlusPlus_0->binMean(0), weight); + _histMeanMultiSigma_c_PlusPlus_0->fill(_histMeanMultiSigma_c_PlusPlus_0->binMean(0), weight); break; case 3124: - _hist10MeanMultiLambda1520->fill(_hist10MeanMultiLambda1520->binMean(0), weight); + _histMeanMultiLambda1520->fill(_histMeanMultiLambda1520->binMean(0), weight); break; } } } - if (2*meanBeamMom/GeV >= 29 && 2*meanBeamMom/GeV <= 35) { + if (sqrtS()/GeV >= 29 && sqrtS()/GeV <= 35) { foreach (const Particle& p, ufs.particles()) { const PdgId id = abs(p.pdgId()); switch (id) { case 211: - _hist32MeanMultiPiPlus->fill(_hist32MeanMultiPiPlus->binMean(0), weight); + _histMeanMultiPiPlus->fill(_histMeanMultiPiPlus->binMean(0), weight); break; case 111: - _hist32MeanMultiPi0->fill(_hist32MeanMultiPi0->binMean(0), weight); + _histMeanMultiPi0->fill(_histMeanMultiPi0->binMean(0), weight); break; case 321: - _hist32MeanMultiKPlus->fill(_hist32MeanMultiKPlus->binMean(0), weight); + _histMeanMultiKPlus->fill(_histMeanMultiKPlus->binMean(0), weight); break; case 130: case 310: - _hist32MeanMultiK0->fill(_hist32MeanMultiK0->binMean(0), weight); + _histMeanMultiK0->fill(_histMeanMultiK0->binMean(0), weight); break; case 221: - _hist32MeanMultiEta->fill(_hist32MeanMultiEta->binMean(0), weight); + _histMeanMultiEta->fill(_histMeanMultiEta->binMean(0), weight); break; case 331: - _hist32MeanMultiEtaPrime->fill(_hist32MeanMultiEtaPrime->binMean(0), weight); + _histMeanMultiEtaPrime->fill(_histMeanMultiEtaPrime->binMean(0), weight); break; case 411: - _hist32MeanMultiDPlus->fill(_hist32MeanMultiDPlus->binMean(0), weight); + _histMeanMultiDPlus->fill(_histMeanMultiDPlus->binMean(0), weight); break; case 421: - _hist32MeanMultiD0->fill(_hist32MeanMultiD0->binMean(0), weight); + _histMeanMultiD0->fill(_histMeanMultiD0->binMean(0), weight); break; case 431: - _hist32MeanMultiDPlus_s->fill(_hist32MeanMultiDPlus_s->binMean(0), weight); + _histMeanMultiDPlus_s->fill(_histMeanMultiDPlus_s->binMean(0), weight); break; case 9010221: - _hist32MeanMultiF0_980->fill(_hist32MeanMultiF0_980->binMean(0), weight); + _histMeanMultiF0_980->fill(_histMeanMultiF0_980->binMean(0), weight); break; case 113: - _hist32MeanMultiRho770_0->fill(_hist32MeanMultiRho770_0->binMean(0), weight); + _histMeanMultiRho770_0->fill(_histMeanMultiRho770_0->binMean(0), weight); break; case 323: - _hist32MeanMultiKStar892Plus->fill(_hist32MeanMultiKStar892Plus->binMean(0), weight); + _histMeanMultiKStar892Plus->fill(_histMeanMultiKStar892Plus->binMean(0), weight); break; case 313: - _hist32MeanMultiKStar892_0->fill(_hist32MeanMultiKStar892_0->binMean(0), weight); + _histMeanMultiKStar892_0->fill(_histMeanMultiKStar892_0->binMean(0), weight); break; case 333: - _hist32MeanMultiPhi1020->fill(_hist32MeanMultiPhi1020->binMean(0), weight); + _histMeanMultiPhi1020->fill(_histMeanMultiPhi1020->binMean(0), weight); break; case 413: - _hist32MeanMultiDStar2010Plus->fill(_hist32MeanMultiDStar2010Plus->binMean(0), weight); + _histMeanMultiDStar2010Plus->fill(_histMeanMultiDStar2010Plus->binMean(0), weight); break; case 423: - _hist32MeanMultiDStar2007_0->fill(_hist32MeanMultiDStar2007_0->binMean(0), weight); + _histMeanMultiDStar2007_0->fill(_histMeanMultiDStar2007_0->binMean(0), weight); break; case 225: - _hist32MeanMultiF2_1270->fill(_hist32MeanMultiF2_1270->binMean(0), weight); + _histMeanMultiF2_1270->fill(_histMeanMultiF2_1270->binMean(0), weight); break; case 325: - _hist32MeanMultiK2Star1430Plus->fill(_hist32MeanMultiK2Star1430Plus->binMean(0), weight); + _histMeanMultiK2Star1430Plus->fill(_histMeanMultiK2Star1430Plus->binMean(0), weight); break; case 315: - _hist32MeanMultiK2Star1430_0->fill(_hist32MeanMultiK2Star1430_0->binMean(0), weight); + _histMeanMultiK2Star1430_0->fill(_histMeanMultiK2Star1430_0->binMean(0), weight); break; case 2212: - _hist32MeanMultiP->fill(_hist32MeanMultiP->binMean(0), weight); + _histMeanMultiP->fill(_histMeanMultiP->binMean(0), weight); break; case 3122: - _hist32MeanMultiLambda->fill(_hist32MeanMultiLambda->binMean(0), weight); + _histMeanMultiLambda->fill(_histMeanMultiLambda->binMean(0), weight); break; case 3312: - _hist32MeanMultiXiMinus->fill(_hist32MeanMultiXiMinus->binMean(0), weight); + _histMeanMultiXiMinus->fill(_histMeanMultiXiMinus->binMean(0), weight); break; case 3114: - _hist32MeanMultiSigma1385Minus->fill(_hist32MeanMultiSigma1385Minus->binMean(0), weight); - _hist32MeanMultiSigma1385PlusMinus->fill(_hist32MeanMultiSigma1385PlusMinus->binMean(0), weight); + _histMeanMultiSigma1385Minus->fill(_histMeanMultiSigma1385Minus->binMean(0), weight); + _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->binMean(0), weight); break; case 3224: - _hist32MeanMultiSigma1385Plus->fill(_hist32MeanMultiSigma1385Plus->binMean(0), weight); - _hist32MeanMultiSigma1385PlusMinus->fill(_hist32MeanMultiSigma1385PlusMinus->binMean(0), weight); + _histMeanMultiSigma1385Plus->fill(_histMeanMultiSigma1385Plus->binMean(0), weight); + _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->binMean(0), weight); break; case 3334: - _hist32MeanMultiOmegaMinus->fill(_hist32MeanMultiOmegaMinus->binMean(0), weight); + _histMeanMultiOmegaMinus->fill(_histMeanMultiOmegaMinus->binMean(0), weight); break; case 4122: - _hist32MeanMultiLambda_c_Plus->fill(_hist32MeanMultiLambda_c_Plus->binMean(0), weight); + _histMeanMultiLambda_c_Plus->fill(_histMeanMultiLambda_c_Plus->binMean(0), weight); break; } } } - - - if (2*meanBeamMom/GeV >= 89.5 && 2*meanBeamMom/GeV <= 91.8) { + if (sqrtS()/GeV >= 89.5 && sqrtS()/GeV <= 91.8) { foreach (const Particle& p, ufs.particles()) { const PdgId id = abs(p.pdgId()); switch (id) { case 211: - _hist91MeanMultiPiPlus->fill(_hist91MeanMultiPiPlus->binMean(0), weight); + _histMeanMultiPiPlus->fill(_histMeanMultiPiPlus->binMean(0), weight); break; case 111: - _hist91MeanMultiPi0->fill(_hist91MeanMultiPi0->binMean(0), weight); + _histMeanMultiPi0->fill(_histMeanMultiPi0->binMean(0), weight); break; case 321: - _hist91MeanMultiKPlus->fill(_hist91MeanMultiKPlus->binMean(0), weight); + _histMeanMultiKPlus->fill(_histMeanMultiKPlus->binMean(0), weight); break; case 130: case 310: - _hist91MeanMultiK0->fill(_hist91MeanMultiK0->binMean(0), weight); + _histMeanMultiK0->fill(_histMeanMultiK0->binMean(0), weight); break; case 221: - _hist91MeanMultiEta->fill(_hist91MeanMultiEta->binMean(0), weight); + _histMeanMultiEta->fill(_histMeanMultiEta->binMean(0), weight); break; case 331: - _hist91MeanMultiEtaPrime->fill(_hist91MeanMultiEtaPrime->binMean(0), weight); + _histMeanMultiEtaPrime->fill(_histMeanMultiEtaPrime->binMean(0), weight); break; case 411: - _hist91MeanMultiDPlus->fill(_hist91MeanMultiDPlus->binMean(0), weight); + _histMeanMultiDPlus->fill(_histMeanMultiDPlus->binMean(0), weight); break; case 421: - _hist91MeanMultiD0->fill(_hist91MeanMultiD0->binMean(0), weight); + _histMeanMultiD0->fill(_histMeanMultiD0->binMean(0), weight); break; case 431: - _hist91MeanMultiDPlus_s->fill(_hist91MeanMultiDPlus_s->binMean(0), weight); + _histMeanMultiDPlus_s->fill(_histMeanMultiDPlus_s->binMean(0), weight); break; case 511: - _hist91MeanMultiBPlus_B0_d->fill(_hist91MeanMultiBPlus_B0_d->binMean(0), weight); + _histMeanMultiBPlus_B0_d->fill(_histMeanMultiBPlus_B0_d->binMean(0), weight); break; case 521: - _hist91MeanMultiBPlus_B0_d->fill(_hist91MeanMultiBPlus_B0_d->binMean(0), weight); - _hist91MeanMultiBPlus_u->fill(_hist91MeanMultiBPlus_u->binMean(0), weight); + _histMeanMultiBPlus_B0_d->fill(_histMeanMultiBPlus_B0_d->binMean(0), weight); + _histMeanMultiBPlus_u->fill(_histMeanMultiBPlus_u->binMean(0), weight); break; case 531: - _hist91MeanMultiB0_s->fill(_hist91MeanMultiB0_s->binMean(0), weight); + _histMeanMultiB0_s->fill(_histMeanMultiB0_s->binMean(0), weight); break; case 9010221: - _hist91MeanMultiF0_980->fill(_hist91MeanMultiF0_980->binMean(0), weight); + _histMeanMultiF0_980->fill(_histMeanMultiF0_980->binMean(0), weight); break; case 9000211: - _hist91MeanMultiA0_980Plus->fill(_hist91MeanMultiA0_980Plus->binMean(0), weight); + _histMeanMultiA0_980Plus->fill(_histMeanMultiA0_980Plus->binMean(0), weight); break; case 113: - _hist91MeanMultiRho770_0->fill(_hist91MeanMultiRho770_0->binMean(0), weight); + _histMeanMultiRho770_0->fill(_histMeanMultiRho770_0->binMean(0), weight); break; case 213: - _hist91MeanMultiRho770Plus->fill(_hist91MeanMultiRho770Plus->binMean(0), weight); + _histMeanMultiRho770Plus->fill(_histMeanMultiRho770Plus->binMean(0), weight); break; case 223: - _hist91MeanMultiOmega782->fill(_hist91MeanMultiOmega782->binMean(0), weight); + _histMeanMultiOmega782->fill(_histMeanMultiOmega782->binMean(0), weight); break; case 323: - _hist91MeanMultiKStar892Plus->fill(_hist91MeanMultiKStar892Plus->binMean(0), weight); + _histMeanMultiKStar892Plus->fill(_histMeanMultiKStar892Plus->binMean(0), weight); break; case 313: - _hist91MeanMultiKStar892_0->fill(_hist91MeanMultiKStar892_0->binMean(0), weight); + _histMeanMultiKStar892_0->fill(_histMeanMultiKStar892_0->binMean(0), weight); break; case 333: - _hist91MeanMultiPhi1020->fill(_hist91MeanMultiPhi1020->binMean(0), weight); + _histMeanMultiPhi1020->fill(_histMeanMultiPhi1020->binMean(0), weight); break; case 413: - _hist91MeanMultiDStar2010Plus->fill(_hist91MeanMultiDStar2010Plus->binMean(0), weight); + _histMeanMultiDStar2010Plus->fill(_histMeanMultiDStar2010Plus->binMean(0), weight); break; case 433: - _hist91MeanMultiDStar_s2112Plus->fill(_hist91MeanMultiDStar_s2112Plus->binMean(0), weight); + _histMeanMultiDStar_s2112Plus->fill(_histMeanMultiDStar_s2112Plus->binMean(0), weight); break; case 513: case 523: case 533: - _hist91MeanMultiBStar->fill(_hist91MeanMultiBStar->binMean(0), weight); + _histMeanMultiBStar->fill(_histMeanMultiBStar->binMean(0), weight); break; case 443: - _hist91MeanMultiJPsi1S->fill(_hist91MeanMultiJPsi1S->binMean(0), weight); + _histMeanMultiJPsi1S->fill(_histMeanMultiJPsi1S->binMean(0), weight); break; case 100443: - _hist91MeanMultiPsi2S->fill(_hist91MeanMultiPsi2S->binMean(0), weight); + _histMeanMultiPsi2S->fill(_histMeanMultiPsi2S->binMean(0), weight); break; case 553: - _hist91MeanMultiUpsilon1S->fill(_hist91MeanMultiUpsilon1S->binMean(0), weight); + _histMeanMultiUpsilon1S->fill(_histMeanMultiUpsilon1S->binMean(0), weight); break; case 20223: - _hist91MeanMultiF1_1285->fill(_hist91MeanMultiF1_1285->binMean(0), weight); + _histMeanMultiF1_1285->fill(_histMeanMultiF1_1285->binMean(0), weight); break; case 20333: - _hist91MeanMultiF1_1420->fill(_hist91MeanMultiF1_1420->binMean(0), weight); + _histMeanMultiF1_1420->fill(_histMeanMultiF1_1420->binMean(0), weight); break; case 445: - _hist91MeanMultiChi_c1_3510->fill(_hist91MeanMultiChi_c1_3510->binMean(0), weight); + _histMeanMultiChi_c1_3510->fill(_histMeanMultiChi_c1_3510->binMean(0), weight); break; case 225: - _hist91MeanMultiF2_1270->fill(_hist91MeanMultiF2_1270->binMean(0), weight); + _histMeanMultiF2_1270->fill(_histMeanMultiF2_1270->binMean(0), weight); break; case 335: - _hist91MeanMultiF2Prime1525->fill(_hist91MeanMultiF2Prime1525->binMean(0), weight); + _histMeanMultiF2Prime1525->fill(_histMeanMultiF2Prime1525->binMean(0), weight); break; case 315: - _hist91MeanMultiK2Star1430_0->fill(_hist91MeanMultiK2Star1430_0->binMean(0), weight); + _histMeanMultiK2Star1430_0->fill(_histMeanMultiK2Star1430_0->binMean(0), weight); break; case 515: case 525: case 535: - _hist91MeanMultiBStarStar->fill(_hist91MeanMultiBStarStar->binMean(0), weight); + _histMeanMultiBStarStar->fill(_histMeanMultiBStarStar->binMean(0), weight); break; case 10433: case 20433: - _hist91MeanMultiDs1Plus->fill(_hist91MeanMultiDs1Plus->binMean(0), weight); + _histMeanMultiDs1Plus->fill(_histMeanMultiDs1Plus->binMean(0), weight); break; case 435: - _hist91MeanMultiDs2Plus->fill(_hist91MeanMultiDs2Plus->binMean(0), weight); + _histMeanMultiDs2Plus->fill(_histMeanMultiDs2Plus->binMean(0), weight); break; case 2212: - _hist91MeanMultiP->fill(_hist91MeanMultiP->binMean(0), weight); + _histMeanMultiP->fill(_histMeanMultiP->binMean(0), weight); break; case 3122: - _hist91MeanMultiLambda->fill(_hist91MeanMultiLambda->binMean(0), weight); + _histMeanMultiLambda->fill(_histMeanMultiLambda->binMean(0), weight); break; case 3212: - _hist91MeanMultiSigma0->fill(_hist91MeanMultiSigma0->binMean(0), weight); + _histMeanMultiSigma0->fill(_histMeanMultiSigma0->binMean(0), weight); break; case 3112: - _hist91MeanMultiSigmaMinus->fill(_hist91MeanMultiSigmaMinus->binMean(0), weight); - _hist91MeanMultiSigmaPlusMinus->fill(_hist91MeanMultiSigmaPlusMinus->binMean(0), weight); + _histMeanMultiSigmaMinus->fill(_histMeanMultiSigmaMinus->binMean(0), weight); + _histMeanMultiSigmaPlusMinus->fill(_histMeanMultiSigmaPlusMinus->binMean(0), weight); break; case 3222: - _hist91MeanMultiSigmaPlus->fill(_hist91MeanMultiSigmaPlus->binMean(0), weight); - _hist91MeanMultiSigmaPlusMinus->fill(_hist91MeanMultiSigmaPlusMinus->binMean(0), weight); + _histMeanMultiSigmaPlus->fill(_histMeanMultiSigmaPlus->binMean(0), weight); + _histMeanMultiSigmaPlusMinus->fill(_histMeanMultiSigmaPlusMinus->binMean(0), weight); break; case 3312: - _hist91MeanMultiXiMinus->fill(_hist91MeanMultiXiMinus->binMean(0), weight); + _histMeanMultiXiMinus->fill(_histMeanMultiXiMinus->binMean(0), weight); break; case 2224: - _hist91MeanMultiDelta1232PlusPlus->fill(_hist91MeanMultiDelta1232PlusPlus->binMean(0), weight); + _histMeanMultiDelta1232PlusPlus->fill(_histMeanMultiDelta1232PlusPlus->binMean(0), weight); break; case 3114: - _hist91MeanMultiSigma1385Minus->fill(_hist91MeanMultiSigma1385Minus->binMean(0), weight); - _hist91MeanMultiSigma1385PlusMinus->fill(_hist91MeanMultiSigma1385PlusMinus->binMean(0), weight); + _histMeanMultiSigma1385Minus->fill(_histMeanMultiSigma1385Minus->binMean(0), weight); + _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->binMean(0), weight); break; case 3224: - _hist91MeanMultiSigma1385Plus->fill(_hist91MeanMultiSigma1385Plus->binMean(0), weight); - _hist91MeanMultiSigma1385PlusMinus->fill(_hist91MeanMultiSigma1385PlusMinus->binMean(0), weight); + _histMeanMultiSigma1385Plus->fill(_histMeanMultiSigma1385Plus->binMean(0), weight); + _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->binMean(0), weight); break; case 3324: - _hist91MeanMultiXi1530_0->fill(_hist91MeanMultiXi1530_0->binMean(0), weight); + _histMeanMultiXi1530_0->fill(_histMeanMultiXi1530_0->binMean(0), weight); break; case 3334: - _hist91MeanMultiOmegaMinus->fill(_hist91MeanMultiOmegaMinus->binMean(0), weight); + _histMeanMultiOmegaMinus->fill(_histMeanMultiOmegaMinus->binMean(0), weight); break; case 4122: - _hist91MeanMultiLambda_c_Plus->fill(_hist91MeanMultiLambda_c_Plus->binMean(0), weight); + _histMeanMultiLambda_c_Plus->fill(_histMeanMultiLambda_c_Plus->binMean(0), weight); break; case 5122: - _hist91MeanMultiLambda_b_0->fill(_hist91MeanMultiLambda_b_0->binMean(0), weight); + _histMeanMultiLambda_b_0->fill(_histMeanMultiLambda_b_0->binMean(0), weight); break; case 3124: - _hist91MeanMultiLambda1520->fill(_hist91MeanMultiLambda1520->binMean(0), weight); + _histMeanMultiLambda1520->fill(_histMeanMultiLambda1520->binMean(0), weight); break; } } } - - - if (2*meanBeamMom/GeV >= 130 && 2*meanBeamMom/GeV <= 200) { + if (sqrtS()/GeV >= 130 && sqrtS()/GeV <= 200) { foreach (const Particle& p, ufs.particles()) { const PdgId id = abs(p.pdgId()); switch (id) { case 211: - _hist165MeanMultiPiPlus->fill(_hist165MeanMultiPiPlus->binMean(0), weight); + _histMeanMultiPiPlus->fill(_histMeanMultiPiPlus->binMean(0), weight); break; case 321: - _hist165MeanMultiKPlus->fill(_hist165MeanMultiKPlus->binMean(0), weight); + _histMeanMultiKPlus->fill(_histMeanMultiKPlus->binMean(0), weight); break; case 130: case 310: - _hist165MeanMultiK0->fill(_hist165MeanMultiK0->binMean(0), weight); + _histMeanMultiK0->fill(_histMeanMultiK0->binMean(0), weight); break; case 2212: - _hist165MeanMultiP->fill(_hist165MeanMultiP->binMean(0), weight); + _histMeanMultiP->fill(_histMeanMultiP->binMean(0), weight); break; case 3122: - _hist165MeanMultiLambda->fill(_hist165MeanMultiLambda->binMean(0), weight); + _histMeanMultiLambda->fill(_histMeanMultiLambda->binMean(0), weight); break; } } } - } void init() { - addProjection(Beam(), "Beams"); addProjection(ChargedFinalState(), "FS"); addProjection(UnstableFinalState(), "UFS"); - _hist10MeanMultiPiPlus = bookHistogram1D( 1, 1, 1); - _hist10MeanMultiPi0 = bookHistogram1D( 2, 1, 1); - _hist10MeanMultiKPlus = bookHistogram1D( 3, 1, 1); - _hist10MeanMultiK0 = bookHistogram1D( 4, 1, 1); - _hist10MeanMultiEta = bookHistogram1D( 5, 1, 1); - _hist10MeanMultiEtaPrime = bookHistogram1D( 6, 1, 1); - _hist10MeanMultiDPlus = bookHistogram1D( 7, 1, 1); - _hist10MeanMultiD0 = bookHistogram1D( 8, 1, 1); - _hist10MeanMultiDPlus_s = bookHistogram1D( 9, 1, 1); - _hist10MeanMultiF0_980 = bookHistogram1D(13, 1, 1); - _hist10MeanMultiRho770_0 = bookHistogram1D(15, 1, 1); - _hist10MeanMultiOmega782 = bookHistogram1D(17, 1, 1); - _hist10MeanMultiKStar892Plus = bookHistogram1D(18, 1, 1); - _hist10MeanMultiKStar892_0 = bookHistogram1D(19, 1, 1); - _hist10MeanMultiPhi1020 = bookHistogram1D(20, 1, 1); - _hist10MeanMultiDStar2010Plus = bookHistogram1D(21, 1, 1); - _hist10MeanMultiDStar2007_0 = bookHistogram1D(22, 1, 1); - _hist10MeanMultiDStar_s2112Plus = bookHistogram1D(23, 1, 1); - _hist10MeanMultiJPsi1S = bookHistogram1D(25, 1, 1); - _hist10MeanMultiF2_1270 = bookHistogram1D(31, 1, 1); - _hist10MeanMultiP = bookHistogram1D(38, 1, 1); - _hist10MeanMultiLambda = bookHistogram1D(39, 1, 1); - _hist10MeanMultiSigma0 = bookHistogram1D(40, 1, 1); - _hist10MeanMultiXiMinus = bookHistogram1D(44, 1, 1); - _hist10MeanMultiDelta1232PlusPlus = bookHistogram1D(45, 1, 1); - _hist10MeanMultiSigma1385Minus = bookHistogram1D(46, 1, 1); - _hist10MeanMultiSigma1385Plus = bookHistogram1D(47, 1, 1); - _hist10MeanMultiSigma1385PlusMinus = bookHistogram1D(48, 1, 1); - _hist10MeanMultiXi1530_0 = bookHistogram1D(49, 1, 1); - _hist10MeanMultiOmegaMinus = bookHistogram1D(50, 1, 1); - _hist10MeanMultiLambda_c_Plus = bookHistogram1D(51, 1, 1); - _hist10MeanMultiSigma_c_PlusPlus_0 = bookHistogram1D(53, 1, 1); - _hist10MeanMultiLambda1520 = bookHistogram1D(54, 1, 1); - - _hist32MeanMultiPiPlus = bookHistogram1D( 1, 1, 2); - _hist32MeanMultiPi0 = bookHistogram1D( 2, 1, 2); - _hist32MeanMultiKPlus = bookHistogram1D( 3, 1, 2); - _hist32MeanMultiK0 = bookHistogram1D( 4, 1, 2); - _hist32MeanMultiEta = bookHistogram1D( 5, 1, 2); - _hist32MeanMultiEtaPrime = bookHistogram1D( 6, 1, 2); - _hist32MeanMultiDPlus = bookHistogram1D( 7, 1, 2); - _hist32MeanMultiD0 = bookHistogram1D( 8, 1, 2); - _hist32MeanMultiDPlus_s = bookHistogram1D( 9, 1, 2); - _hist32MeanMultiF0_980 = bookHistogram1D(13, 1, 2); - _hist32MeanMultiRho770_0 = bookHistogram1D(15, 1, 2); - _hist32MeanMultiKStar892Plus = bookHistogram1D(18, 1, 2); - _hist32MeanMultiKStar892_0 = bookHistogram1D(19, 1, 2); - _hist32MeanMultiPhi1020 = bookHistogram1D(20, 1, 2); - _hist32MeanMultiDStar2010Plus = bookHistogram1D(21, 1, 2); - _hist32MeanMultiDStar2007_0 = bookHistogram1D(22, 1, 2); - _hist32MeanMultiF2_1270 = bookHistogram1D(31, 1, 2); - _hist32MeanMultiK2Star1430Plus = bookHistogram1D(33, 1, 1); - _hist32MeanMultiK2Star1430_0 = bookHistogram1D(34, 1, 1); - _hist32MeanMultiP = bookHistogram1D(38, 1, 2); - _hist32MeanMultiLambda = bookHistogram1D(39, 1, 2); - _hist32MeanMultiXiMinus = bookHistogram1D(44, 1, 2); - _hist32MeanMultiSigma1385Minus = bookHistogram1D(46, 1, 2); - _hist32MeanMultiSigma1385Plus = bookHistogram1D(47, 1, 2); - _hist32MeanMultiSigma1385PlusMinus = bookHistogram1D(48, 1, 2); - _hist32MeanMultiOmegaMinus = bookHistogram1D(50, 1, 2); - _hist32MeanMultiLambda_c_Plus = bookHistogram1D(51, 1, 2); - - _hist91MeanMultiPiPlus = bookHistogram1D( 1, 1, 3); - _hist91MeanMultiPi0 = bookHistogram1D( 2, 1, 3); - _hist91MeanMultiKPlus = bookHistogram1D( 3, 1, 3); - _hist91MeanMultiK0 = bookHistogram1D( 4, 1, 3); - _hist91MeanMultiEta = bookHistogram1D( 5, 1, 3); - _hist91MeanMultiEtaPrime = bookHistogram1D( 6, 1, 3); - _hist91MeanMultiDPlus = bookHistogram1D( 7, 1, 3); - _hist91MeanMultiD0 = bookHistogram1D( 8, 1, 3); - _hist91MeanMultiDPlus_s = bookHistogram1D( 9, 1, 3); - _hist91MeanMultiBPlus_B0_d = bookHistogram1D(10, 1, 1); - _hist91MeanMultiBPlus_u = bookHistogram1D(11, 1, 1); - _hist91MeanMultiB0_s = bookHistogram1D(12, 1, 1); - _hist91MeanMultiF0_980 = bookHistogram1D(13, 1, 3); - _hist91MeanMultiA0_980Plus = bookHistogram1D(14, 1, 1); - _hist91MeanMultiRho770_0 = bookHistogram1D(15, 1, 3); - _hist91MeanMultiRho770Plus = bookHistogram1D(16, 1, 1); - _hist91MeanMultiOmega782 = bookHistogram1D(17, 1, 2); - _hist91MeanMultiKStar892Plus = bookHistogram1D(18, 1, 3); - _hist91MeanMultiKStar892_0 = bookHistogram1D(19, 1, 3); - _hist91MeanMultiPhi1020 = bookHistogram1D(20, 1, 3); - _hist91MeanMultiDStar2010Plus = bookHistogram1D(21, 1, 3); - _hist91MeanMultiDStar_s2112Plus = bookHistogram1D(23, 1, 2); - _hist91MeanMultiBStar = bookHistogram1D(24, 1, 1); - _hist91MeanMultiJPsi1S = bookHistogram1D(25, 1, 2); - _hist91MeanMultiPsi2S = bookHistogram1D(26, 1, 1); - _hist91MeanMultiUpsilon1S = bookHistogram1D(27, 1, 1); - _hist91MeanMultiF1_1285 = bookHistogram1D(28, 1, 1); - _hist91MeanMultiF1_1420 = bookHistogram1D(29, 1, 1); - _hist91MeanMultiChi_c1_3510 = bookHistogram1D(30, 1, 1); - _hist91MeanMultiF2_1270 = bookHistogram1D(31, 1, 3); - _hist91MeanMultiF2Prime1525 = bookHistogram1D(32, 1, 1); - _hist91MeanMultiK2Star1430_0 = bookHistogram1D(34, 1, 2); - _hist91MeanMultiBStarStar = bookHistogram1D(35, 1, 1); - _hist91MeanMultiDs1Plus = bookHistogram1D(36, 1, 1); - _hist91MeanMultiDs2Plus = bookHistogram1D(37, 1, 1); - _hist91MeanMultiP = bookHistogram1D(38, 1, 3); - _hist91MeanMultiLambda = bookHistogram1D(39, 1, 3); - _hist91MeanMultiSigma0 = bookHistogram1D(40, 1, 2); - _hist91MeanMultiSigmaMinus = bookHistogram1D(41, 1, 1); - _hist91MeanMultiSigmaPlus = bookHistogram1D(42, 1, 1); - _hist91MeanMultiSigmaPlusMinus = bookHistogram1D(43, 1, 1); - _hist91MeanMultiXiMinus = bookHistogram1D(44, 1, 3); - _hist91MeanMultiDelta1232PlusPlus = bookHistogram1D(45, 1, 2); - _hist91MeanMultiSigma1385Minus = bookHistogram1D(46, 1, 3); - _hist91MeanMultiSigma1385Plus = bookHistogram1D(47, 1, 3); - _hist91MeanMultiSigma1385PlusMinus = bookHistogram1D(48, 1, 3); - _hist91MeanMultiXi1530_0 = bookHistogram1D(49, 1, 2); - _hist91MeanMultiOmegaMinus = bookHistogram1D(50, 1, 3); - _hist91MeanMultiLambda_c_Plus = bookHistogram1D(51, 1, 3); - _hist91MeanMultiLambda_b_0 = bookHistogram1D(52, 1, 1); - _hist91MeanMultiLambda1520 = bookHistogram1D(54, 1, 2); - - _hist165MeanMultiPiPlus = bookHistogram1D( 1, 1, 4); - _hist165MeanMultiKPlus = bookHistogram1D( 3, 1, 4); - _hist165MeanMultiK0 = bookHistogram1D( 4, 1, 4); - _hist165MeanMultiP = bookHistogram1D(38, 1, 4); - _hist165MeanMultiLambda = bookHistogram1D(39, 1, 4); + if (sqrtS()/GeV >= 9.5 && sqrtS()/GeV <= 10.5) { + _histMeanMultiPiPlus = bookHistogram1D( 1, 1, 1); + _histMeanMultiPi0 = bookHistogram1D( 2, 1, 1); + _histMeanMultiKPlus = bookHistogram1D( 3, 1, 1); + _histMeanMultiK0 = bookHistogram1D( 4, 1, 1); + _histMeanMultiEta = bookHistogram1D( 5, 1, 1); + _histMeanMultiEtaPrime = bookHistogram1D( 6, 1, 1); + _histMeanMultiDPlus = bookHistogram1D( 7, 1, 1); + _histMeanMultiD0 = bookHistogram1D( 8, 1, 1); + _histMeanMultiDPlus_s = bookHistogram1D( 9, 1, 1); + _histMeanMultiF0_980 = bookHistogram1D(13, 1, 1); + _histMeanMultiRho770_0 = bookHistogram1D(15, 1, 1); + _histMeanMultiOmega782 = bookHistogram1D(17, 1, 1); + _histMeanMultiKStar892Plus = bookHistogram1D(18, 1, 1); + _histMeanMultiKStar892_0 = bookHistogram1D(19, 1, 1); + _histMeanMultiPhi1020 = bookHistogram1D(20, 1, 1); + _histMeanMultiDStar2010Plus = bookHistogram1D(21, 1, 1); + _histMeanMultiDStar2007_0 = bookHistogram1D(22, 1, 1); + _histMeanMultiDStar_s2112Plus = bookHistogram1D(23, 1, 1); + _histMeanMultiJPsi1S = bookHistogram1D(25, 1, 1); + _histMeanMultiF2_1270 = bookHistogram1D(31, 1, 1); + _histMeanMultiP = bookHistogram1D(38, 1, 1); + _histMeanMultiLambda = bookHistogram1D(39, 1, 1); + _histMeanMultiSigma0 = bookHistogram1D(40, 1, 1); + _histMeanMultiXiMinus = bookHistogram1D(44, 1, 1); + _histMeanMultiDelta1232PlusPlus = bookHistogram1D(45, 1, 1); + _histMeanMultiSigma1385Minus = bookHistogram1D(46, 1, 1); + _histMeanMultiSigma1385Plus = bookHistogram1D(47, 1, 1); + _histMeanMultiSigma1385PlusMinus = bookHistogram1D(48, 1, 1); + _histMeanMultiXi1530_0 = bookHistogram1D(49, 1, 1); + _histMeanMultiOmegaMinus = bookHistogram1D(50, 1, 1); + _histMeanMultiLambda_c_Plus = bookHistogram1D(51, 1, 1); + _histMeanMultiSigma_c_PlusPlus_0 = bookHistogram1D(53, 1, 1); + _histMeanMultiLambda1520 = bookHistogram1D(54, 1, 1); + } + + if (sqrtS()/GeV >= 29 && sqrtS()/GeV <= 35) { + _histMeanMultiPiPlus = bookHistogram1D( 1, 1, 2); + _histMeanMultiPi0 = bookHistogram1D( 2, 1, 2); + _histMeanMultiKPlus = bookHistogram1D( 3, 1, 2); + _histMeanMultiK0 = bookHistogram1D( 4, 1, 2); + _histMeanMultiEta = bookHistogram1D( 5, 1, 2); + _histMeanMultiEtaPrime = bookHistogram1D( 6, 1, 2); + _histMeanMultiDPlus = bookHistogram1D( 7, 1, 2); + _histMeanMultiD0 = bookHistogram1D( 8, 1, 2); + _histMeanMultiDPlus_s = bookHistogram1D( 9, 1, 2); + _histMeanMultiF0_980 = bookHistogram1D(13, 1, 2); + _histMeanMultiRho770_0 = bookHistogram1D(15, 1, 2); + _histMeanMultiKStar892Plus = bookHistogram1D(18, 1, 2); + _histMeanMultiKStar892_0 = bookHistogram1D(19, 1, 2); + _histMeanMultiPhi1020 = bookHistogram1D(20, 1, 2); + _histMeanMultiDStar2010Plus = bookHistogram1D(21, 1, 2); + _histMeanMultiDStar2007_0 = bookHistogram1D(22, 1, 2); + _histMeanMultiF2_1270 = bookHistogram1D(31, 1, 2); + _histMeanMultiK2Star1430Plus = bookHistogram1D(33, 1, 1); + _histMeanMultiK2Star1430_0 = bookHistogram1D(34, 1, 1); + _histMeanMultiP = bookHistogram1D(38, 1, 2); + _histMeanMultiLambda = bookHistogram1D(39, 1, 2); + _histMeanMultiXiMinus = bookHistogram1D(44, 1, 2); + _histMeanMultiSigma1385Minus = bookHistogram1D(46, 1, 2); + _histMeanMultiSigma1385Plus = bookHistogram1D(47, 1, 2); + _histMeanMultiSigma1385PlusMinus = bookHistogram1D(48, 1, 2); + _histMeanMultiOmegaMinus = bookHistogram1D(50, 1, 2); + _histMeanMultiLambda_c_Plus = bookHistogram1D(51, 1, 2); + } + + if (sqrtS()/GeV >= 89.5 && sqrtS()/GeV <= 91.8) { + _histMeanMultiPiPlus = bookHistogram1D( 1, 1, 3); + _histMeanMultiPi0 = bookHistogram1D( 2, 1, 3); + _histMeanMultiKPlus = bookHistogram1D( 3, 1, 3); + _histMeanMultiK0 = bookHistogram1D( 4, 1, 3); + _histMeanMultiEta = bookHistogram1D( 5, 1, 3); + _histMeanMultiEtaPrime = bookHistogram1D( 6, 1, 3); + _histMeanMultiDPlus = bookHistogram1D( 7, 1, 3); + _histMeanMultiD0 = bookHistogram1D( 8, 1, 3); + _histMeanMultiDPlus_s = bookHistogram1D( 9, 1, 3); + _histMeanMultiBPlus_B0_d = bookHistogram1D(10, 1, 1); + _histMeanMultiBPlus_u = bookHistogram1D(11, 1, 1); + _histMeanMultiB0_s = bookHistogram1D(12, 1, 1); + _histMeanMultiF0_980 = bookHistogram1D(13, 1, 3); + _histMeanMultiA0_980Plus = bookHistogram1D(14, 1, 1); + _histMeanMultiRho770_0 = bookHistogram1D(15, 1, 3); + _histMeanMultiRho770Plus = bookHistogram1D(16, 1, 1); + _histMeanMultiOmega782 = bookHistogram1D(17, 1, 2); + _histMeanMultiKStar892Plus = bookHistogram1D(18, 1, 3); + _histMeanMultiKStar892_0 = bookHistogram1D(19, 1, 3); + _histMeanMultiPhi1020 = bookHistogram1D(20, 1, 3); + _histMeanMultiDStar2010Plus = bookHistogram1D(21, 1, 3); + _histMeanMultiDStar_s2112Plus = bookHistogram1D(23, 1, 2); + _histMeanMultiBStar = bookHistogram1D(24, 1, 1); + _histMeanMultiJPsi1S = bookHistogram1D(25, 1, 2); + _histMeanMultiPsi2S = bookHistogram1D(26, 1, 1); + _histMeanMultiUpsilon1S = bookHistogram1D(27, 1, 1); + _histMeanMultiF1_1285 = bookHistogram1D(28, 1, 1); + _histMeanMultiF1_1420 = bookHistogram1D(29, 1, 1); + _histMeanMultiChi_c1_3510 = bookHistogram1D(30, 1, 1); + _histMeanMultiF2_1270 = bookHistogram1D(31, 1, 3); + _histMeanMultiF2Prime1525 = bookHistogram1D(32, 1, 1); + _histMeanMultiK2Star1430_0 = bookHistogram1D(34, 1, 2); + _histMeanMultiBStarStar = bookHistogram1D(35, 1, 1); + _histMeanMultiDs1Plus = bookHistogram1D(36, 1, 1); + _histMeanMultiDs2Plus = bookHistogram1D(37, 1, 1); + _histMeanMultiP = bookHistogram1D(38, 1, 3); + _histMeanMultiLambda = bookHistogram1D(39, 1, 3); + _histMeanMultiSigma0 = bookHistogram1D(40, 1, 2); + _histMeanMultiSigmaMinus = bookHistogram1D(41, 1, 1); + _histMeanMultiSigmaPlus = bookHistogram1D(42, 1, 1); + _histMeanMultiSigmaPlusMinus = bookHistogram1D(43, 1, 1); + _histMeanMultiXiMinus = bookHistogram1D(44, 1, 3); + _histMeanMultiDelta1232PlusPlus = bookHistogram1D(45, 1, 2); + _histMeanMultiSigma1385Minus = bookHistogram1D(46, 1, 3); + _histMeanMultiSigma1385Plus = bookHistogram1D(47, 1, 3); + _histMeanMultiSigma1385PlusMinus = bookHistogram1D(48, 1, 3); + _histMeanMultiXi1530_0 = bookHistogram1D(49, 1, 2); + _histMeanMultiOmegaMinus = bookHistogram1D(50, 1, 3); + _histMeanMultiLambda_c_Plus = bookHistogram1D(51, 1, 3); + _histMeanMultiLambda_b_0 = bookHistogram1D(52, 1, 1); + _histMeanMultiLambda1520 = bookHistogram1D(54, 1, 2); + } + + if (sqrtS()/GeV >= 130 && sqrtS()/GeV <= 200) { + _histMeanMultiPiPlus = bookHistogram1D( 1, 1, 4); + _histMeanMultiKPlus = bookHistogram1D( 3, 1, 4); + _histMeanMultiK0 = bookHistogram1D( 4, 1, 4); + _histMeanMultiP = bookHistogram1D(38, 1, 4); + _histMeanMultiLambda = bookHistogram1D(39, 1, 4); + } } // Finalize void finalize() { - scale(_hist10MeanMultiPiPlus , 1.0/sumOfWeights()); - scale(_hist10MeanMultiPi0 , 1.0/sumOfWeights()); - scale(_hist10MeanMultiKPlus , 1.0/sumOfWeights()); - scale(_hist10MeanMultiK0 , 1.0/sumOfWeights()); - scale(_hist10MeanMultiEta , 1.0/sumOfWeights()); - scale(_hist10MeanMultiEtaPrime , 1.0/sumOfWeights()); - scale(_hist10MeanMultiDPlus , 1.0/sumOfWeights()); - scale(_hist10MeanMultiD0 , 1.0/sumOfWeights()); - scale(_hist10MeanMultiDPlus_s , 1.0/sumOfWeights()); - scale(_hist10MeanMultiF0_980 , 1.0/sumOfWeights()); - scale(_hist10MeanMultiRho770_0 , 1.0/sumOfWeights()); - scale(_hist10MeanMultiOmega782 , 1.0/sumOfWeights()); - scale(_hist10MeanMultiKStar892Plus , 1.0/sumOfWeights()); - scale(_hist10MeanMultiKStar892_0 , 1.0/sumOfWeights()); - scale(_hist10MeanMultiPhi1020 , 1.0/sumOfWeights()); - scale(_hist10MeanMultiDStar2010Plus , 1.0/sumOfWeights()); - scale(_hist10MeanMultiDStar2007_0 , 1.0/sumOfWeights()); - scale(_hist10MeanMultiDStar_s2112Plus , 1.0/sumOfWeights()); - scale(_hist10MeanMultiJPsi1S , 1.0/sumOfWeights()); - scale(_hist10MeanMultiF2_1270 , 1.0/sumOfWeights()); - scale(_hist10MeanMultiP , 1.0/sumOfWeights()); - scale(_hist10MeanMultiLambda , 1.0/sumOfWeights()); - scale(_hist10MeanMultiSigma0 , 1.0/sumOfWeights()); - scale(_hist10MeanMultiXiMinus , 1.0/sumOfWeights()); - scale(_hist10MeanMultiDelta1232PlusPlus , 1.0/sumOfWeights()); - scale(_hist10MeanMultiSigma1385Minus , 1.0/sumOfWeights()); - scale(_hist10MeanMultiSigma1385Plus , 1.0/sumOfWeights()); - scale(_hist10MeanMultiSigma1385PlusMinus, 1.0/sumOfWeights()); - scale(_hist10MeanMultiXi1530_0 , 1.0/sumOfWeights()); - scale(_hist10MeanMultiOmegaMinus , 1.0/sumOfWeights()); - scale(_hist10MeanMultiLambda_c_Plus , 1.0/sumOfWeights()); - scale(_hist10MeanMultiSigma_c_PlusPlus_0, 1.0/sumOfWeights()); - scale(_hist10MeanMultiLambda1520 , 1.0/sumOfWeights()); - - scale(_hist32MeanMultiPiPlus , 1.0/sumOfWeights()); - scale(_hist32MeanMultiPi0 , 1.0/sumOfWeights()); - scale(_hist32MeanMultiKPlus , 1.0/sumOfWeights()); - scale(_hist32MeanMultiK0 , 1.0/sumOfWeights()); - scale(_hist32MeanMultiEta , 1.0/sumOfWeights()); - scale(_hist32MeanMultiEtaPrime , 1.0/sumOfWeights()); - scale(_hist32MeanMultiDPlus , 1.0/sumOfWeights()); - scale(_hist32MeanMultiD0 , 1.0/sumOfWeights()); - scale(_hist32MeanMultiDPlus_s , 1.0/sumOfWeights()); - scale(_hist32MeanMultiF0_980 , 1.0/sumOfWeights()); - scale(_hist32MeanMultiRho770_0 , 1.0/sumOfWeights()); - scale(_hist32MeanMultiKStar892Plus , 1.0/sumOfWeights()); - scale(_hist32MeanMultiKStar892_0 , 1.0/sumOfWeights()); - scale(_hist32MeanMultiPhi1020 , 1.0/sumOfWeights()); - scale(_hist32MeanMultiDStar2010Plus , 1.0/sumOfWeights()); - scale(_hist32MeanMultiDStar2007_0 , 1.0/sumOfWeights()); - scale(_hist32MeanMultiF2_1270 , 1.0/sumOfWeights()); - scale(_hist32MeanMultiK2Star1430Plus , 1.0/sumOfWeights()); - scale(_hist32MeanMultiK2Star1430_0 , 1.0/sumOfWeights()); - scale(_hist32MeanMultiP , 1.0/sumOfWeights()); - scale(_hist32MeanMultiLambda , 1.0/sumOfWeights()); - scale(_hist32MeanMultiXiMinus , 1.0/sumOfWeights()); - scale(_hist32MeanMultiSigma1385Minus , 1.0/sumOfWeights()); - scale(_hist32MeanMultiSigma1385Plus , 1.0/sumOfWeights()); - scale(_hist32MeanMultiSigma1385PlusMinus, 1.0/sumOfWeights()); - scale(_hist32MeanMultiOmegaMinus , 1.0/sumOfWeights()); - scale(_hist32MeanMultiLambda_c_Plus , 1.0/sumOfWeights()); - - scale(_hist91MeanMultiPiPlus , 1.0/sumOfWeights()); - scale(_hist91MeanMultiPi0 , 1.0/sumOfWeights()); - scale(_hist91MeanMultiKPlus , 1.0/sumOfWeights()); - scale(_hist91MeanMultiK0 , 1.0/sumOfWeights()); - scale(_hist91MeanMultiEta , 1.0/sumOfWeights()); - scale(_hist91MeanMultiEtaPrime , 1.0/sumOfWeights()); - scale(_hist91MeanMultiDPlus , 1.0/sumOfWeights()); - scale(_hist91MeanMultiD0 , 1.0/sumOfWeights()); - scale(_hist91MeanMultiDPlus_s , 1.0/sumOfWeights()); - scale(_hist91MeanMultiBPlus_B0_d , 1.0/sumOfWeights()); - scale(_hist91MeanMultiBPlus_u , 1.0/sumOfWeights()); - scale(_hist91MeanMultiB0_s , 1.0/sumOfWeights()); - scale(_hist91MeanMultiF0_980 , 1.0/sumOfWeights()); - scale(_hist91MeanMultiA0_980Plus , 1.0/sumOfWeights()); - scale(_hist91MeanMultiRho770_0 , 1.0/sumOfWeights()); - scale(_hist91MeanMultiRho770Plus , 1.0/sumOfWeights()); - scale(_hist91MeanMultiOmega782 , 1.0/sumOfWeights()); - scale(_hist91MeanMultiKStar892Plus , 1.0/sumOfWeights()); - scale(_hist91MeanMultiKStar892_0 , 1.0/sumOfWeights()); - scale(_hist91MeanMultiPhi1020 , 1.0/sumOfWeights()); - scale(_hist91MeanMultiDStar2010Plus , 1.0/sumOfWeights()); - scale(_hist91MeanMultiDStar_s2112Plus , 1.0/sumOfWeights()); - scale(_hist91MeanMultiBStar , 1.0/sumOfWeights()); - scale(_hist91MeanMultiJPsi1S , 1.0/sumOfWeights()); - scale(_hist91MeanMultiPsi2S , 1.0/sumOfWeights()); - scale(_hist91MeanMultiUpsilon1S , 1.0/sumOfWeights()); - scale(_hist91MeanMultiF1_1285 , 1.0/sumOfWeights()); - scale(_hist91MeanMultiF1_1420 , 1.0/sumOfWeights()); - scale(_hist91MeanMultiChi_c1_3510 , 1.0/sumOfWeights()); - scale(_hist91MeanMultiF2_1270 , 1.0/sumOfWeights()); - scale(_hist91MeanMultiF2Prime1525 , 1.0/sumOfWeights()); - scale(_hist91MeanMultiK2Star1430_0 , 1.0/sumOfWeights()); - scale(_hist91MeanMultiBStarStar , 1.0/sumOfWeights()); - scale(_hist91MeanMultiDs1Plus , 1.0/sumOfWeights()); - scale(_hist91MeanMultiDs2Plus , 1.0/sumOfWeights()); - scale(_hist91MeanMultiP , 1.0/sumOfWeights()); - scale(_hist91MeanMultiLambda , 1.0/sumOfWeights()); - scale(_hist91MeanMultiSigma0 , 1.0/sumOfWeights()); - scale(_hist91MeanMultiSigmaMinus , 1.0/sumOfWeights()); - scale(_hist91MeanMultiSigmaPlus , 1.0/sumOfWeights()); - scale(_hist91MeanMultiSigmaPlusMinus , 1.0/sumOfWeights()); - scale(_hist91MeanMultiXiMinus , 1.0/sumOfWeights()); - scale(_hist91MeanMultiDelta1232PlusPlus , 1.0/sumOfWeights()); - scale(_hist91MeanMultiSigma1385Minus , 1.0/sumOfWeights()); - scale(_hist91MeanMultiSigma1385Plus , 1.0/sumOfWeights()); - scale(_hist91MeanMultiSigma1385PlusMinus, 1.0/sumOfWeights()); - scale(_hist91MeanMultiXi1530_0 , 1.0/sumOfWeights()); - scale(_hist91MeanMultiOmegaMinus , 1.0/sumOfWeights()); - scale(_hist91MeanMultiLambda_c_Plus , 1.0/sumOfWeights()); - scale(_hist91MeanMultiLambda_b_0 , 1.0/sumOfWeights()); - scale(_hist91MeanMultiLambda1520 , 1.0/sumOfWeights()); - - scale(_hist165MeanMultiPiPlus , 1.0/sumOfWeights()); - scale(_hist165MeanMultiKPlus , 1.0/sumOfWeights()); - scale(_hist165MeanMultiK0 , 1.0/sumOfWeights()); - scale(_hist165MeanMultiP , 1.0/sumOfWeights()); - scale(_hist165MeanMultiLambda , 1.0/sumOfWeights()); - } + if (sqrtS()/GeV >= 9.5 && sqrtS()/GeV <= 10.5) { + scale(_histMeanMultiPiPlus , 1.0/sumOfWeights()); + scale(_histMeanMultiPi0 , 1.0/sumOfWeights()); + scale(_histMeanMultiKPlus , 1.0/sumOfWeights()); + scale(_histMeanMultiK0 , 1.0/sumOfWeights()); + scale(_histMeanMultiEta , 1.0/sumOfWeights()); + scale(_histMeanMultiEtaPrime , 1.0/sumOfWeights()); + scale(_histMeanMultiDPlus , 1.0/sumOfWeights()); + scale(_histMeanMultiD0 , 1.0/sumOfWeights()); + scale(_histMeanMultiDPlus_s , 1.0/sumOfWeights()); + scale(_histMeanMultiF0_980 , 1.0/sumOfWeights()); + scale(_histMeanMultiRho770_0 , 1.0/sumOfWeights()); + scale(_histMeanMultiOmega782 , 1.0/sumOfWeights()); + scale(_histMeanMultiKStar892Plus , 1.0/sumOfWeights()); + scale(_histMeanMultiKStar892_0 , 1.0/sumOfWeights()); + scale(_histMeanMultiPhi1020 , 1.0/sumOfWeights()); + scale(_histMeanMultiDStar2010Plus , 1.0/sumOfWeights()); + scale(_histMeanMultiDStar2007_0 , 1.0/sumOfWeights()); + scale(_histMeanMultiDStar_s2112Plus , 1.0/sumOfWeights()); + scale(_histMeanMultiJPsi1S , 1.0/sumOfWeights()); + scale(_histMeanMultiF2_1270 , 1.0/sumOfWeights()); + scale(_histMeanMultiP , 1.0/sumOfWeights()); + scale(_histMeanMultiLambda , 1.0/sumOfWeights()); + scale(_histMeanMultiSigma0 , 1.0/sumOfWeights()); + scale(_histMeanMultiXiMinus , 1.0/sumOfWeights()); + scale(_histMeanMultiDelta1232PlusPlus , 1.0/sumOfWeights()); + scale(_histMeanMultiSigma1385Minus , 1.0/sumOfWeights()); + scale(_histMeanMultiSigma1385Plus , 1.0/sumOfWeights()); + scale(_histMeanMultiSigma1385PlusMinus, 1.0/sumOfWeights()); + scale(_histMeanMultiXi1530_0 , 1.0/sumOfWeights()); + scale(_histMeanMultiOmegaMinus , 1.0/sumOfWeights()); + scale(_histMeanMultiLambda_c_Plus , 1.0/sumOfWeights()); + scale(_histMeanMultiSigma_c_PlusPlus_0, 1.0/sumOfWeights()); + scale(_histMeanMultiLambda1520 , 1.0/sumOfWeights()); + } + if (sqrtS()/GeV >= 29 && sqrtS()/GeV <= 35) { + scale(_histMeanMultiPiPlus , 1.0/sumOfWeights()); + scale(_histMeanMultiPi0 , 1.0/sumOfWeights()); + scale(_histMeanMultiKPlus , 1.0/sumOfWeights()); + scale(_histMeanMultiK0 , 1.0/sumOfWeights()); + scale(_histMeanMultiEta , 1.0/sumOfWeights()); + scale(_histMeanMultiEtaPrime , 1.0/sumOfWeights()); + scale(_histMeanMultiDPlus , 1.0/sumOfWeights()); + scale(_histMeanMultiD0 , 1.0/sumOfWeights()); + scale(_histMeanMultiDPlus_s , 1.0/sumOfWeights()); + scale(_histMeanMultiF0_980 , 1.0/sumOfWeights()); + scale(_histMeanMultiRho770_0 , 1.0/sumOfWeights()); + scale(_histMeanMultiKStar892Plus , 1.0/sumOfWeights()); + scale(_histMeanMultiKStar892_0 , 1.0/sumOfWeights()); + scale(_histMeanMultiPhi1020 , 1.0/sumOfWeights()); + scale(_histMeanMultiDStar2010Plus , 1.0/sumOfWeights()); + scale(_histMeanMultiDStar2007_0 , 1.0/sumOfWeights()); + scale(_histMeanMultiF2_1270 , 1.0/sumOfWeights()); + scale(_histMeanMultiK2Star1430Plus , 1.0/sumOfWeights()); + scale(_histMeanMultiK2Star1430_0 , 1.0/sumOfWeights()); + scale(_histMeanMultiP , 1.0/sumOfWeights()); + scale(_histMeanMultiLambda , 1.0/sumOfWeights()); + scale(_histMeanMultiXiMinus , 1.0/sumOfWeights()); + scale(_histMeanMultiSigma1385Minus , 1.0/sumOfWeights()); + scale(_histMeanMultiSigma1385Plus , 1.0/sumOfWeights()); + scale(_histMeanMultiSigma1385PlusMinus, 1.0/sumOfWeights()); + scale(_histMeanMultiOmegaMinus , 1.0/sumOfWeights()); + scale(_histMeanMultiLambda_c_Plus , 1.0/sumOfWeights()); + } + + if (sqrtS()/GeV >= 89.5 && sqrtS()/GeV <= 91.8) { + scale(_histMeanMultiPiPlus , 1.0/sumOfWeights()); + scale(_histMeanMultiPi0 , 1.0/sumOfWeights()); + scale(_histMeanMultiKPlus , 1.0/sumOfWeights()); + scale(_histMeanMultiK0 , 1.0/sumOfWeights()); + scale(_histMeanMultiEta , 1.0/sumOfWeights()); + scale(_histMeanMultiEtaPrime , 1.0/sumOfWeights()); + scale(_histMeanMultiDPlus , 1.0/sumOfWeights()); + scale(_histMeanMultiD0 , 1.0/sumOfWeights()); + scale(_histMeanMultiDPlus_s , 1.0/sumOfWeights()); + scale(_histMeanMultiBPlus_B0_d , 1.0/sumOfWeights()); + scale(_histMeanMultiBPlus_u , 1.0/sumOfWeights()); + scale(_histMeanMultiB0_s , 1.0/sumOfWeights()); + scale(_histMeanMultiF0_980 , 1.0/sumOfWeights()); + scale(_histMeanMultiA0_980Plus , 1.0/sumOfWeights()); + scale(_histMeanMultiRho770_0 , 1.0/sumOfWeights()); + scale(_histMeanMultiRho770Plus , 1.0/sumOfWeights()); + scale(_histMeanMultiOmega782 , 1.0/sumOfWeights()); + scale(_histMeanMultiKStar892Plus , 1.0/sumOfWeights()); + scale(_histMeanMultiKStar892_0 , 1.0/sumOfWeights()); + scale(_histMeanMultiPhi1020 , 1.0/sumOfWeights()); + scale(_histMeanMultiDStar2010Plus , 1.0/sumOfWeights()); + scale(_histMeanMultiDStar_s2112Plus , 1.0/sumOfWeights()); + scale(_histMeanMultiBStar , 1.0/sumOfWeights()); + scale(_histMeanMultiJPsi1S , 1.0/sumOfWeights()); + scale(_histMeanMultiPsi2S , 1.0/sumOfWeights()); + scale(_histMeanMultiUpsilon1S , 1.0/sumOfWeights()); + scale(_histMeanMultiF1_1285 , 1.0/sumOfWeights()); + scale(_histMeanMultiF1_1420 , 1.0/sumOfWeights()); + scale(_histMeanMultiChi_c1_3510 , 1.0/sumOfWeights()); + scale(_histMeanMultiF2_1270 , 1.0/sumOfWeights()); + scale(_histMeanMultiF2Prime1525 , 1.0/sumOfWeights()); + scale(_histMeanMultiK2Star1430_0 , 1.0/sumOfWeights()); + scale(_histMeanMultiBStarStar , 1.0/sumOfWeights()); + scale(_histMeanMultiDs1Plus , 1.0/sumOfWeights()); + scale(_histMeanMultiDs2Plus , 1.0/sumOfWeights()); + scale(_histMeanMultiP , 1.0/sumOfWeights()); + scale(_histMeanMultiLambda , 1.0/sumOfWeights()); + scale(_histMeanMultiSigma0 , 1.0/sumOfWeights()); + scale(_histMeanMultiSigmaMinus , 1.0/sumOfWeights()); + scale(_histMeanMultiSigmaPlus , 1.0/sumOfWeights()); + scale(_histMeanMultiSigmaPlusMinus , 1.0/sumOfWeights()); + scale(_histMeanMultiXiMinus , 1.0/sumOfWeights()); + scale(_histMeanMultiDelta1232PlusPlus , 1.0/sumOfWeights()); + scale(_histMeanMultiSigma1385Minus , 1.0/sumOfWeights()); + scale(_histMeanMultiSigma1385Plus , 1.0/sumOfWeights()); + scale(_histMeanMultiSigma1385PlusMinus, 1.0/sumOfWeights()); + scale(_histMeanMultiXi1530_0 , 1.0/sumOfWeights()); + scale(_histMeanMultiOmegaMinus , 1.0/sumOfWeights()); + scale(_histMeanMultiLambda_c_Plus , 1.0/sumOfWeights()); + scale(_histMeanMultiLambda_b_0 , 1.0/sumOfWeights()); + scale(_histMeanMultiLambda1520 , 1.0/sumOfWeights()); + } + + if (sqrtS()/GeV >= 130 && sqrtS()/GeV <= 200) { + scale(_histMeanMultiPiPlus , 1.0/sumOfWeights()); + scale(_histMeanMultiKPlus , 1.0/sumOfWeights()); + scale(_histMeanMultiK0 , 1.0/sumOfWeights()); + scale(_histMeanMultiP , 1.0/sumOfWeights()); + scale(_histMeanMultiLambda , 1.0/sumOfWeights()); + } + } //@} private: - AIDA::IHistogram1D *_hist10MeanMultiPiPlus; - AIDA::IHistogram1D *_hist10MeanMultiPi0; - AIDA::IHistogram1D *_hist10MeanMultiKPlus; - AIDA::IHistogram1D *_hist10MeanMultiK0; - AIDA::IHistogram1D *_hist10MeanMultiEta; - AIDA::IHistogram1D *_hist10MeanMultiEtaPrime; - AIDA::IHistogram1D *_hist10MeanMultiDPlus; - AIDA::IHistogram1D *_hist10MeanMultiD0; - AIDA::IHistogram1D *_hist10MeanMultiDPlus_s; - AIDA::IHistogram1D *_hist10MeanMultiF0_980; - AIDA::IHistogram1D *_hist10MeanMultiRho770_0; - AIDA::IHistogram1D *_hist10MeanMultiOmega782; - AIDA::IHistogram1D *_hist10MeanMultiKStar892Plus; - AIDA::IHistogram1D *_hist10MeanMultiKStar892_0; - AIDA::IHistogram1D *_hist10MeanMultiPhi1020; - AIDA::IHistogram1D *_hist10MeanMultiDStar2010Plus; - AIDA::IHistogram1D *_hist10MeanMultiDStar2007_0; - AIDA::IHistogram1D *_hist10MeanMultiDStar_s2112Plus; - AIDA::IHistogram1D *_hist10MeanMultiJPsi1S; - AIDA::IHistogram1D *_hist10MeanMultiF2_1270; - AIDA::IHistogram1D *_hist10MeanMultiP; - AIDA::IHistogram1D *_hist10MeanMultiLambda; - AIDA::IHistogram1D *_hist10MeanMultiSigma0; - AIDA::IHistogram1D *_hist10MeanMultiXiMinus; - AIDA::IHistogram1D *_hist10MeanMultiDelta1232PlusPlus; - AIDA::IHistogram1D *_hist10MeanMultiSigma1385Minus; - AIDA::IHistogram1D *_hist10MeanMultiSigma1385Plus; - AIDA::IHistogram1D *_hist10MeanMultiSigma1385PlusMinus; - AIDA::IHistogram1D *_hist10MeanMultiXi1530_0; - AIDA::IHistogram1D *_hist10MeanMultiOmegaMinus; - AIDA::IHistogram1D *_hist10MeanMultiLambda_c_Plus; - AIDA::IHistogram1D *_hist10MeanMultiSigma_c_PlusPlus_0; - AIDA::IHistogram1D *_hist10MeanMultiLambda1520; - - AIDA::IHistogram1D *_hist32MeanMultiPiPlus; - AIDA::IHistogram1D *_hist32MeanMultiPi0; - AIDA::IHistogram1D *_hist32MeanMultiKPlus; - AIDA::IHistogram1D *_hist32MeanMultiK0; - AIDA::IHistogram1D *_hist32MeanMultiEta; - AIDA::IHistogram1D *_hist32MeanMultiEtaPrime; - AIDA::IHistogram1D *_hist32MeanMultiDPlus; - AIDA::IHistogram1D *_hist32MeanMultiD0; - AIDA::IHistogram1D *_hist32MeanMultiDPlus_s; - AIDA::IHistogram1D *_hist32MeanMultiF0_980; - AIDA::IHistogram1D *_hist32MeanMultiRho770_0; - AIDA::IHistogram1D *_hist32MeanMultiKStar892Plus; - AIDA::IHistogram1D *_hist32MeanMultiKStar892_0; - AIDA::IHistogram1D *_hist32MeanMultiPhi1020; - AIDA::IHistogram1D *_hist32MeanMultiDStar2010Plus; - AIDA::IHistogram1D *_hist32MeanMultiDStar2007_0; - AIDA::IHistogram1D *_hist32MeanMultiF2_1270; - AIDA::IHistogram1D *_hist32MeanMultiK2Star1430Plus; - AIDA::IHistogram1D *_hist32MeanMultiK2Star1430_0; - AIDA::IHistogram1D *_hist32MeanMultiP; - AIDA::IHistogram1D *_hist32MeanMultiLambda; - AIDA::IHistogram1D *_hist32MeanMultiXiMinus; - AIDA::IHistogram1D *_hist32MeanMultiSigma1385Minus; - AIDA::IHistogram1D *_hist32MeanMultiSigma1385Plus; - AIDA::IHistogram1D *_hist32MeanMultiSigma1385PlusMinus; - AIDA::IHistogram1D *_hist32MeanMultiOmegaMinus; - AIDA::IHistogram1D *_hist32MeanMultiLambda_c_Plus; - - AIDA::IHistogram1D *_hist91MeanMultiPiPlus; - AIDA::IHistogram1D *_hist91MeanMultiPi0; - AIDA::IHistogram1D *_hist91MeanMultiKPlus; - AIDA::IHistogram1D *_hist91MeanMultiK0; - AIDA::IHistogram1D *_hist91MeanMultiEta; - AIDA::IHistogram1D *_hist91MeanMultiEtaPrime; - AIDA::IHistogram1D *_hist91MeanMultiDPlus; - AIDA::IHistogram1D *_hist91MeanMultiD0; - AIDA::IHistogram1D *_hist91MeanMultiDPlus_s; - AIDA::IHistogram1D *_hist91MeanMultiBPlus_B0_d; - AIDA::IHistogram1D *_hist91MeanMultiBPlus_u; - AIDA::IHistogram1D *_hist91MeanMultiB0_s; - AIDA::IHistogram1D *_hist91MeanMultiF0_980; - AIDA::IHistogram1D *_hist91MeanMultiA0_980Plus; - AIDA::IHistogram1D *_hist91MeanMultiRho770_0; - AIDA::IHistogram1D *_hist91MeanMultiRho770Plus; - AIDA::IHistogram1D *_hist91MeanMultiOmega782; - AIDA::IHistogram1D *_hist91MeanMultiKStar892Plus; - AIDA::IHistogram1D *_hist91MeanMultiKStar892_0; - AIDA::IHistogram1D *_hist91MeanMultiPhi1020; - AIDA::IHistogram1D *_hist91MeanMultiDStar2010Plus; - AIDA::IHistogram1D *_hist91MeanMultiDStar_s2112Plus; - AIDA::IHistogram1D *_hist91MeanMultiBStar; - AIDA::IHistogram1D *_hist91MeanMultiJPsi1S; - AIDA::IHistogram1D *_hist91MeanMultiPsi2S; - AIDA::IHistogram1D *_hist91MeanMultiUpsilon1S; - AIDA::IHistogram1D *_hist91MeanMultiF1_1285; - AIDA::IHistogram1D *_hist91MeanMultiF1_1420; - AIDA::IHistogram1D *_hist91MeanMultiChi_c1_3510; - AIDA::IHistogram1D *_hist91MeanMultiF2_1270; - AIDA::IHistogram1D *_hist91MeanMultiF2Prime1525; - AIDA::IHistogram1D *_hist91MeanMultiK2Star1430_0; - AIDA::IHistogram1D *_hist91MeanMultiBStarStar; - AIDA::IHistogram1D *_hist91MeanMultiDs1Plus; - AIDA::IHistogram1D *_hist91MeanMultiDs2Plus; - AIDA::IHistogram1D *_hist91MeanMultiP; - AIDA::IHistogram1D *_hist91MeanMultiLambda; - AIDA::IHistogram1D *_hist91MeanMultiSigma0; - AIDA::IHistogram1D *_hist91MeanMultiSigmaMinus; - AIDA::IHistogram1D *_hist91MeanMultiSigmaPlus; - AIDA::IHistogram1D *_hist91MeanMultiSigmaPlusMinus; - AIDA::IHistogram1D *_hist91MeanMultiXiMinus; - AIDA::IHistogram1D *_hist91MeanMultiDelta1232PlusPlus; - AIDA::IHistogram1D *_hist91MeanMultiSigma1385Minus; - AIDA::IHistogram1D *_hist91MeanMultiSigma1385Plus; - AIDA::IHistogram1D *_hist91MeanMultiSigma1385PlusMinus; - AIDA::IHistogram1D *_hist91MeanMultiXi1530_0; - AIDA::IHistogram1D *_hist91MeanMultiOmegaMinus; - AIDA::IHistogram1D *_hist91MeanMultiLambda_c_Plus; - AIDA::IHistogram1D *_hist91MeanMultiLambda_b_0; - AIDA::IHistogram1D *_hist91MeanMultiLambda1520; - - AIDA::IHistogram1D *_hist165MeanMultiPiPlus; - AIDA::IHistogram1D *_hist165MeanMultiKPlus; - AIDA::IHistogram1D *_hist165MeanMultiK0; - AIDA::IHistogram1D *_hist165MeanMultiP; - AIDA::IHistogram1D *_hist165MeanMultiLambda; + AIDA::IHistogram1D *_histMeanMultiPiPlus; + AIDA::IHistogram1D *_histMeanMultiPi0; + AIDA::IHistogram1D *_histMeanMultiKPlus; + AIDA::IHistogram1D *_histMeanMultiK0; + AIDA::IHistogram1D *_histMeanMultiEta; + AIDA::IHistogram1D *_histMeanMultiEtaPrime; + AIDA::IHistogram1D *_histMeanMultiDPlus; + AIDA::IHistogram1D *_histMeanMultiD0; + AIDA::IHistogram1D *_histMeanMultiDPlus_s; + AIDA::IHistogram1D *_histMeanMultiBPlus_B0_d; + AIDA::IHistogram1D *_histMeanMultiBPlus_u; + AIDA::IHistogram1D *_histMeanMultiB0_s; + AIDA::IHistogram1D *_histMeanMultiF0_980; + AIDA::IHistogram1D *_histMeanMultiA0_980Plus; + AIDA::IHistogram1D *_histMeanMultiRho770_0; + AIDA::IHistogram1D *_histMeanMultiRho770Plus; + AIDA::IHistogram1D *_histMeanMultiOmega782; + AIDA::IHistogram1D *_histMeanMultiKStar892Plus; + AIDA::IHistogram1D *_histMeanMultiKStar892_0; + AIDA::IHistogram1D *_histMeanMultiPhi1020; + AIDA::IHistogram1D *_histMeanMultiDStar2010Plus; + AIDA::IHistogram1D *_histMeanMultiDStar2007_0; + AIDA::IHistogram1D *_histMeanMultiDStar_s2112Plus; + AIDA::IHistogram1D *_histMeanMultiBStar; + AIDA::IHistogram1D *_histMeanMultiJPsi1S; + AIDA::IHistogram1D *_histMeanMultiPsi2S; + AIDA::IHistogram1D *_histMeanMultiUpsilon1S; + AIDA::IHistogram1D *_histMeanMultiF1_1285; + AIDA::IHistogram1D *_histMeanMultiF1_1420; + AIDA::IHistogram1D *_histMeanMultiChi_c1_3510; + AIDA::IHistogram1D *_histMeanMultiF2_1270; + AIDA::IHistogram1D *_histMeanMultiF2Prime1525; + AIDA::IHistogram1D *_histMeanMultiK2Star1430Plus; + AIDA::IHistogram1D *_histMeanMultiK2Star1430_0; + AIDA::IHistogram1D *_histMeanMultiBStarStar; + AIDA::IHistogram1D *_histMeanMultiDs1Plus; + AIDA::IHistogram1D *_histMeanMultiDs2Plus; + AIDA::IHistogram1D *_histMeanMultiP; + AIDA::IHistogram1D *_histMeanMultiLambda; + AIDA::IHistogram1D *_histMeanMultiSigma0; + AIDA::IHistogram1D *_histMeanMultiSigmaMinus; + AIDA::IHistogram1D *_histMeanMultiSigmaPlus; + AIDA::IHistogram1D *_histMeanMultiSigmaPlusMinus; + AIDA::IHistogram1D *_histMeanMultiXiMinus; + AIDA::IHistogram1D *_histMeanMultiDelta1232PlusPlus; + AIDA::IHistogram1D *_histMeanMultiSigma1385Minus; + AIDA::IHistogram1D *_histMeanMultiSigma1385Plus; + AIDA::IHistogram1D *_histMeanMultiSigma1385PlusMinus; + AIDA::IHistogram1D *_histMeanMultiXi1530_0; + AIDA::IHistogram1D *_histMeanMultiOmegaMinus; + AIDA::IHistogram1D *_histMeanMultiLambda_c_Plus; + AIDA::IHistogram1D *_histMeanMultiLambda_b_0; + AIDA::IHistogram1D *_histMeanMultiSigma_c_PlusPlus_0; + AIDA::IHistogram1D *_histMeanMultiLambda1520; //@} Modified: trunk/src/Analyses/PDG_Hadron_Multiplicities_Ratios.cc ============================================================================== --- trunk/src/Analyses/PDG_Hadron_Multiplicities_Ratios.cc Mon Feb 8 18:41:49 2010 (r2245) +++ trunk/src/Analyses/PDG_Hadron_Multiplicities_Ratios.cc Tue Feb 9 11:18:01 2010 (r2246) @@ -2,7 +2,6 @@ #include "Rivet/Analysis.hh" #include "Rivet/RivetAIDA.hh" #include "Rivet/Tools/ParticleIdUtils.hh" -#include "Rivet/Projections/Beam.hh" #include "Rivet/Projections/FinalState.hh" #include "Rivet/Projections/ChargedFinalState.hh" #include "Rivet/Projections/UnstableFinalState.hh" @@ -14,19 +13,15 @@ /// @author Holger Schulz class PDG_HADRON_MULTIPLICITIES_RATIOS : public Analysis { public: - + /// Constructor PDG_HADRON_MULTIPLICITIES_RATIOS() : Analysis("PDG_HADRON_MULTIPLICITIES_RATIOS") { setBeams(ELECTRON, POSITRON); - - _weightedTotalNumPiPlus10 = 0; - _weightedTotalNumPiPlus32 = 0; - _weightedTotalNumPiPlus91 = 0; - _weightedTotalNumPiPlus165 = 0; + _weightedTotalNumPiPlus = 0; } - - + + /// @name Analysis methods //@{ @@ -34,663 +29,669 @@ // First, veto on leptonic events by requiring at least 4 charged FS particles const FinalState& fs = applyProjection<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) { getLog() << Log::DEBUG << "Failed leptonic event cut" << endl; vetoEvent; } getLog() << Log::DEBUG << "Passed leptonic event cut" << endl; - + // Get event weight for histo filling const double weight = e.weight(); - // Get beams and average beam momentum - const ParticlePair& beams = applyProjection<Beam>(e, "Beams").beams(); - const double meanBeamMom = ( beams.first.momentum().vector3().mod() + - beams.second.momentum().vector3().mod() ) / 2.0; - getLog() << Log::DEBUG << "Avg beam momentum = " << meanBeamMom/GeV << " GeV" << endl; + getLog() << Log::DEBUG << "sqrt(S) = " << sqrtS()/GeV << " GeV" << endl; // Final state of unstable particles to get particle spectra const UnstableFinalState& ufs = applyProjection<UnstableFinalState>(e, "UFS"); - if (2*meanBeamMom/GeV >= 9.5 && 2*meanBeamMom/GeV <= 10.5) { + if (sqrtS()/GeV >= 9.5 && sqrtS()/GeV <= 10.5) { foreach (const Particle& p, ufs.particles()) { const PdgId id = abs(p.pdgId()); switch (id) { case 211: - _weightedTotalNumPiPlus10 += weight; + _weightedTotalNumPiPlus += weight; break; case 111: - _hist10MeanMultiPi0->fill(_hist10MeanMultiPi0->binMean(0), weight); + _histMeanMultiPi0->fill(_histMeanMultiPi0->binMean(0), weight); break; case 321: - _hist10MeanMultiKPlus->fill(_hist10MeanMultiKPlus->binMean(0), weight); + _histMeanMultiKPlus->fill(_histMeanMultiKPlus->binMean(0), weight); break; case 130: case 310: - _hist10MeanMultiK0->fill(_hist10MeanMultiK0->binMean(0), weight); + _histMeanMultiK0->fill(_histMeanMultiK0->binMean(0), weight); break; case 221: - _hist10MeanMultiEta->fill(_hist10MeanMultiEta->binMean(0), weight); + _histMeanMultiEta->fill(_histMeanMultiEta->binMean(0), weight); break; case 331: - _hist10MeanMultiEtaPrime->fill(_hist10MeanMultiEtaPrime->binMean(0), weight); + _histMeanMultiEtaPrime->fill(_histMeanMultiEtaPrime->binMean(0), weight); break; case 411: - _hist10MeanMultiDPlus->fill(_hist10MeanMultiDPlus->binMean(0), weight); + _histMeanMultiDPlus->fill(_histMeanMultiDPlus->binMean(0), weight); break; case 421: - _hist10MeanMultiD0->fill(_hist10MeanMultiD0->binMean(0), weight); + _histMeanMultiD0->fill(_histMeanMultiD0->binMean(0), weight); break; case 431: - _hist10MeanMultiDPlus_s->fill(_hist10MeanMultiDPlus_s->binMean(0), weight); + _histMeanMultiDPlus_s->fill(_histMeanMultiDPlus_s->binMean(0), weight); break; case 9010221: - _hist10MeanMultiF0_980->fill(_hist10MeanMultiF0_980->binMean(0), weight); + _histMeanMultiF0_980->fill(_histMeanMultiF0_980->binMean(0), weight); break; case 113: - _hist10MeanMultiRho770_0->fill(_hist10MeanMultiRho770_0->binMean(0), weight); + _histMeanMultiRho770_0->fill(_histMeanMultiRho770_0->binMean(0), weight); break; case 223: - _hist10MeanMultiOmega782->fill(_hist10MeanMultiOmega782->binMean(0), weight); + _histMeanMultiOmega782->fill(_histMeanMultiOmega782->binMean(0), weight); break; case 323: - _hist10MeanMultiKStar892Plus->fill(_hist10MeanMultiKStar892Plus->binMean(0), weight); + _histMeanMultiKStar892Plus->fill(_histMeanMultiKStar892Plus->binMean(0), weight); break; case 313: - _hist10MeanMultiKStar892_0->fill(_hist10MeanMultiKStar892_0->binMean(0), weight); + _histMeanMultiKStar892_0->fill(_histMeanMultiKStar892_0->binMean(0), weight); break; case 333: - _hist10MeanMultiPhi1020->fill(_hist10MeanMultiPhi1020->binMean(0), weight); + _histMeanMultiPhi1020->fill(_histMeanMultiPhi1020->binMean(0), weight); break; case 413: - _hist10MeanMultiDStar2010Plus->fill(_hist10MeanMultiDStar2010Plus->binMean(0), weight); + _histMeanMultiDStar2010Plus->fill(_histMeanMultiDStar2010Plus->binMean(0), weight); break; case 423: - _hist10MeanMultiDStar2007_0->fill(_hist10MeanMultiDStar2007_0->binMean(0), weight); + _histMeanMultiDStar2007_0->fill(_histMeanMultiDStar2007_0->binMean(0), weight); break; case 433: - _hist10MeanMultiDStar_s2112Plus->fill(_hist10MeanMultiDStar_s2112Plus->binMean(0), weight); + _histMeanMultiDStar_s2112Plus->fill(_histMeanMultiDStar_s2112Plus->binMean(0), weight); break; case 443: - _hist10MeanMultiJPsi1S->fill(_hist10MeanMultiJPsi1S->binMean(0), weight); + _histMeanMultiJPsi1S->fill(_histMeanMultiJPsi1S->binMean(0), weight); break; case 225: - _hist10MeanMultiF2_1270->fill(_hist10MeanMultiF2_1270->binMean(0), weight); + _histMeanMultiF2_1270->fill(_histMeanMultiF2_1270->binMean(0), weight); break; case 2212: - _hist10MeanMultiP->fill(_hist10MeanMultiP->binMean(0), weight); + _histMeanMultiP->fill(_histMeanMultiP->binMean(0), weight); break; case 3122: - _hist10MeanMultiLambda->fill(_hist10MeanMultiLambda->binMean(0), weight); + _histMeanMultiLambda->fill(_histMeanMultiLambda->binMean(0), weight); break; case 3212: - _hist10MeanMultiSigma0->fill(_hist10MeanMultiSigma0->binMean(0), weight); + _histMeanMultiSigma0->fill(_histMeanMultiSigma0->binMean(0), weight); break; case 3312: - _hist10MeanMultiXiMinus->fill(_hist10MeanMultiXiMinus->binMean(0), weight); + _histMeanMultiXiMinus->fill(_histMeanMultiXiMinus->binMean(0), weight); break; case 2224: - _hist10MeanMultiDelta1232PlusPlus->fill(_hist10MeanMultiDelta1232PlusPlus->binMean(0), weight); + _histMeanMultiDelta1232PlusPlus->fill(_histMeanMultiDelta1232PlusPlus->binMean(0), weight); break; case 3114: - _hist10MeanMultiSigma1385Minus->fill(_hist10MeanMultiSigma1385Minus->binMean(0), weight); - _hist10MeanMultiSigma1385PlusMinus->fill(_hist10MeanMultiSigma1385PlusMinus->binMean(0), weight); + _histMeanMultiSigma1385Minus->fill(_histMeanMultiSigma1385Minus->binMean(0), weight); + _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->binMean(0), weight); break; case 3224: - _hist10MeanMultiSigma1385Plus->fill(_hist10MeanMultiSigma1385Plus->binMean(0), weight); - _hist10MeanMultiSigma1385PlusMinus->fill(_hist10MeanMultiSigma1385PlusMinus->binMean(0), weight); + _histMeanMultiSigma1385Plus->fill(_histMeanMultiSigma1385Plus->binMean(0), weight); + _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->binMean(0), weight); break; case 3324: - _hist10MeanMultiXi1530_0->fill(_hist10MeanMultiXi1530_0->binMean(0), weight); + _histMeanMultiXi1530_0->fill(_histMeanMultiXi1530_0->binMean(0), weight); break; case 3334: - _hist10MeanMultiOmegaMinus->fill(_hist10MeanMultiOmegaMinus->binMean(0), weight); + _histMeanMultiOmegaMinus->fill(_histMeanMultiOmegaMinus->binMean(0), weight); break; case 4122: - _hist10MeanMultiLambda_c_Plus->fill(_hist10MeanMultiLambda_c_Plus->binMean(0), weight); + _histMeanMultiLambda_c_Plus->fill(_histMeanMultiLambda_c_Plus->binMean(0), weight); break; case 4222: case 4112: - _hist10MeanMultiSigma_c_PlusPlus_0->fill(_hist10MeanMultiSigma_c_PlusPlus_0->binMean(0), weight); + _histMeanMultiSigma_c_PlusPlus_0->fill(_histMeanMultiSigma_c_PlusPlus_0->binMean(0), weight); break; case 3124: - _hist10MeanMultiLambda1520->fill(_hist10MeanMultiLambda1520->binMean(0), weight); + _histMeanMultiLambda1520->fill(_histMeanMultiLambda1520->binMean(0), weight); break; } } } - if (2*meanBeamMom/GeV >= 29 && 2*meanBeamMom/GeV <= 35) { + if (sqrtS()/GeV >= 29 && sqrtS()/GeV <= 35) { foreach (const Particle& p, ufs.particles()) { const PdgId id = abs(p.pdgId()); switch (id) { case 211: - _weightedTotalNumPiPlus32 += weight; + _weightedTotalNumPiPlus += weight; break; case 111: - _hist32MeanMultiPi0->fill(_hist32MeanMultiPi0->binMean(0), weight); + _histMeanMultiPi0->fill(_histMeanMultiPi0->binMean(0), weight); break; case 321: - _hist32MeanMultiKPlus->fill(_hist32MeanMultiKPlus->binMean(0), weight); + _histMeanMultiKPlus->fill(_histMeanMultiKPlus->binMean(0), weight); break; case 130: case 310: - _hist32MeanMultiK0->fill(_hist32MeanMultiK0->binMean(0), weight); + _histMeanMultiK0->fill(_histMeanMultiK0->binMean(0), weight); break; case 221: - _hist32MeanMultiEta->fill(_hist32MeanMultiEta->binMean(0), weight); + _histMeanMultiEta->fill(_histMeanMultiEta->binMean(0), weight); break; case 331: - _hist32MeanMultiEtaPrime->fill(_hist32MeanMultiEtaPrime->binMean(0), weight); + _histMeanMultiEtaPrime->fill(_histMeanMultiEtaPrime->binMean(0), weight); break; case 411: - _hist32MeanMultiDPlus->fill(_hist32MeanMultiDPlus->binMean(0), weight); + _histMeanMultiDPlus->fill(_histMeanMultiDPlus->binMean(0), weight); break; case 421: - _hist32MeanMultiD0->fill(_hist32MeanMultiD0->binMean(0), weight); + _histMeanMultiD0->fill(_histMeanMultiD0->binMean(0), weight); break; case 431: - _hist32MeanMultiDPlus_s->fill(_hist32MeanMultiDPlus_s->binMean(0), weight); + _histMeanMultiDPlus_s->fill(_histMeanMultiDPlus_s->binMean(0), weight); break; case 9010221: - _hist32MeanMultiF0_980->fill(_hist32MeanMultiF0_980->binMean(0), weight); + _histMeanMultiF0_980->fill(_histMeanMultiF0_980->binMean(0), weight); break; case 113: - _hist32MeanMultiRho770_0->fill(_hist32MeanMultiRho770_0->binMean(0), weight); + _histMeanMultiRho770_0->fill(_histMeanMultiRho770_0->binMean(0), weight); break; case 323: - _hist32MeanMultiKStar892Plus->fill(_hist32MeanMultiKStar892Plus->binMean(0), weight); + _histMeanMultiKStar892Plus->fill(_histMeanMultiKStar892Plus->binMean(0), weight); break; case 313: - _hist32MeanMultiKStar892_0->fill(_hist32MeanMultiKStar892_0->binMean(0), weight); + _histMeanMultiKStar892_0->fill(_histMeanMultiKStar892_0->binMean(0), weight); break; case 333: - _hist32MeanMultiPhi1020->fill(_hist32MeanMultiPhi1020->binMean(0), weight); + _histMeanMultiPhi1020->fill(_histMeanMultiPhi1020->binMean(0), weight); break; case 413: - _hist32MeanMultiDStar2010Plus->fill(_hist32MeanMultiDStar2010Plus->binMean(0), weight); + _histMeanMultiDStar2010Plus->fill(_histMeanMultiDStar2010Plus->binMean(0), weight); break; case 423: - _hist32MeanMultiDStar2007_0->fill(_hist32MeanMultiDStar2007_0->binMean(0), weight); + _histMeanMultiDStar2007_0->fill(_histMeanMultiDStar2007_0->binMean(0), weight); break; case 225: - _hist32MeanMultiF2_1270->fill(_hist32MeanMultiF2_1270->binMean(0), weight); + _histMeanMultiF2_1270->fill(_histMeanMultiF2_1270->binMean(0), weight); break; case 325: - _hist32MeanMultiK2Star1430Plus->fill(_hist32MeanMultiK2Star1430Plus->binMean(0), weight); + _histMeanMultiK2Star1430Plus->fill(_histMeanMultiK2Star1430Plus->binMean(0), weight); break; case 315: - _hist32MeanMultiK2Star1430_0->fill(_hist32MeanMultiK2Star1430_0->binMean(0), weight); + _histMeanMultiK2Star1430_0->fill(_histMeanMultiK2Star1430_0->binMean(0), weight); break; case 2212: - _hist32MeanMultiP->fill(_hist32MeanMultiP->binMean(0), weight); + _histMeanMultiP->fill(_histMeanMultiP->binMean(0), weight); break; case 3122: - _hist32MeanMultiLambda->fill(_hist32MeanMultiLambda->binMean(0), weight); + _histMeanMultiLambda->fill(_histMeanMultiLambda->binMean(0), weight); break; case 3312: - _hist32MeanMultiXiMinus->fill(_hist32MeanMultiXiMinus->binMean(0), weight); + _histMeanMultiXiMinus->fill(_histMeanMultiXiMinus->binMean(0), weight); break; case 3114: - _hist32MeanMultiSigma1385Minus->fill(_hist32MeanMultiSigma1385Minus->binMean(0), weight); - _hist32MeanMultiSigma1385PlusMinus->fill(_hist32MeanMultiSigma1385PlusMinus->binMean(0), weight); + _histMeanMultiSigma1385Minus->fill(_histMeanMultiSigma1385Minus->binMean(0), weight); + _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->binMean(0), weight); break; case 3224: - _hist32MeanMultiSigma1385Plus->fill(_hist32MeanMultiSigma1385Plus->binMean(0), weight); - _hist32MeanMultiSigma1385PlusMinus->fill(_hist32MeanMultiSigma1385PlusMinus->binMean(0), weight); + _histMeanMultiSigma1385Plus->fill(_histMeanMultiSigma1385Plus->binMean(0), weight); + _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->binMean(0), weight); break; case 3334: - _hist32MeanMultiOmegaMinus->fill(_hist32MeanMultiOmegaMinus->binMean(0), weight); + _histMeanMultiOmegaMinus->fill(_histMeanMultiOmegaMinus->binMean(0), weight); break; case 4122: - _hist32MeanMultiLambda_c_Plus->fill(_hist32MeanMultiLambda_c_Plus->binMean(0), weight); + _histMeanMultiLambda_c_Plus->fill(_histMeanMultiLambda_c_Plus->binMean(0), weight); break; } } } - - - if (2*meanBeamMom/GeV >= 89.5 && 2*meanBeamMom/GeV <= 91.8) { + if (sqrtS()/GeV >= 89.5 && sqrtS()/GeV <= 91.8) { foreach (const Particle& p, ufs.particles()) { const PdgId id = abs(p.pdgId()); switch (id) { case 211: - _weightedTotalNumPiPlus91 += weight; + _weightedTotalNumPiPlus += weight; break; case 111: - _hist91MeanMultiPi0->fill(_hist91MeanMultiPi0->binMean(0), weight); + _histMeanMultiPi0->fill(_histMeanMultiPi0->binMean(0), weight); break; case 321: - _hist91MeanMultiKPlus->fill(_hist91MeanMultiKPlus->binMean(0), weight); + _histMeanMultiKPlus->fill(_histMeanMultiKPlus->binMean(0), weight); break; case 130: case 310: - _hist91MeanMultiK0->fill(_hist91MeanMultiK0->binMean(0), weight); + _histMeanMultiK0->fill(_histMeanMultiK0->binMean(0), weight); break; case 221: - _hist91MeanMultiEta->fill(_hist91MeanMultiEta->binMean(0), weight); + _histMeanMultiEta->fill(_histMeanMultiEta->binMean(0), weight); break; case 331: - _hist91MeanMultiEtaPrime->fill(_hist91MeanMultiEtaPrime->binMean(0), weight); + _histMeanMultiEtaPrime->fill(_histMeanMultiEtaPrime->binMean(0), weight); break; case 411: - _hist91MeanMultiDPlus->fill(_hist91MeanMultiDPlus->binMean(0), weight); + _histMeanMultiDPlus->fill(_histMeanMultiDPlus->binMean(0), weight); break; case 421: - _hist91MeanMultiD0->fill(_hist91MeanMultiD0->binMean(0), weight); + _histMeanMultiD0->fill(_histMeanMultiD0->binMean(0), weight); break; case 431: - _hist91MeanMultiDPlus_s->fill(_hist91MeanMultiDPlus_s->binMean(0), weight); + _histMeanMultiDPlus_s->fill(_histMeanMultiDPlus_s->binMean(0), weight); break; case 511: - _hist91MeanMultiBPlus_B0_d->fill(_hist91MeanMultiBPlus_B0_d->binMean(0), weight); + _histMeanMultiBPlus_B0_d->fill(_histMeanMultiBPlus_B0_d->binMean(0), weight); break; case 521: - _hist91MeanMultiBPlus_B0_d->fill(_hist91MeanMultiBPlus_B0_d->binMean(0), weight); - _hist91MeanMultiBPlus_u->fill(_hist91MeanMultiBPlus_u->binMean(0), weight); + _histMeanMultiBPlus_B0_d->fill(_histMeanMultiBPlus_B0_d->binMean(0), weight); + _histMeanMultiBPlus_u->fill(_histMeanMultiBPlus_u->binMean(0), weight); break; case 531: - _hist91MeanMultiB0_s->fill(_hist91MeanMultiB0_s->binMean(0), weight); + _histMeanMultiB0_s->fill(_histMeanMultiB0_s->binMean(0), weight); break; case 9010221: - _hist91MeanMultiF0_980->fill(_hist91MeanMultiF0_980->binMean(0), weight); + _histMeanMultiF0_980->fill(_histMeanMultiF0_980->binMean(0), weight); break; case 9000211: - _hist91MeanMultiA0_980Plus->fill(_hist91MeanMultiA0_980Plus->binMean(0), weight); + _histMeanMultiA0_980Plus->fill(_histMeanMultiA0_980Plus->binMean(0), weight); break; case 113: - _hist91MeanMultiRho770_0->fill(_hist91MeanMultiRho770_0->binMean(0), weight); + _histMeanMultiRho770_0->fill(_histMeanMultiRho770_0->binMean(0), weight); break; case 213: - _hist91MeanMultiRho770Plus->fill(_hist91MeanMultiRho770Plus->binMean(0), weight); + _histMeanMultiRho770Plus->fill(_histMeanMultiRho770Plus->binMean(0), weight); break; case 223: - _hist91MeanMultiOmega782->fill(_hist91MeanMultiOmega782->binMean(0), weight); + _histMeanMultiOmega782->fill(_histMeanMultiOmega782->binMean(0), weight); break; case 323: - _hist91MeanMultiKStar892Plus->fill(_hist91MeanMultiKStar892Plus->binMean(0), weight); + _histMeanMultiKStar892Plus->fill(_histMeanMultiKStar892Plus->binMean(0), weight); break; case 313: - _hist91MeanMultiKStar892_0->fill(_hist91MeanMultiKStar892_0->binMean(0), weight); + _histMeanMultiKStar892_0->fill(_histMeanMultiKStar892_0->binMean(0), weight); break; case 333: - _hist91MeanMultiPhi1020->fill(_hist91MeanMultiPhi1020->binMean(0), weight); + _histMeanMultiPhi1020->fill(_histMeanMultiPhi1020->binMean(0), weight); break; case 413: - _hist91MeanMultiDStar2010Plus->fill(_hist91MeanMultiDStar2010Plus->binMean(0), weight); + _histMeanMultiDStar2010Plus->fill(_histMeanMultiDStar2010Plus->binMean(0), weight); break; case 433: - _hist91MeanMultiDStar_s2112Plus->fill(_hist91MeanMultiDStar_s2112Plus->binMean(0), weight); + _histMeanMultiDStar_s2112Plus->fill(_histMeanMultiDStar_s2112Plus->binMean(0), weight); break; case 513: case 523: case 533: - _hist91MeanMultiBStar->fill(_hist91MeanMultiBStar->binMean(0), weight); + _histMeanMultiBStar->fill(_histMeanMultiBStar->binMean(0), weight); break; case 443: - _hist91MeanMultiJPsi1S->fill(_hist91MeanMultiJPsi1S->binMean(0), weight); + _histMeanMultiJPsi1S->fill(_histMeanMultiJPsi1S->binMean(0), weight); break; case 100443: - _hist91MeanMultiPsi2S->fill(_hist91MeanMultiPsi2S->binMean(0), weight); + _histMeanMultiPsi2S->fill(_histMeanMultiPsi2S->binMean(0), weight); break; case 553: - _hist91MeanMultiUpsilon1S->fill(_hist91MeanMultiUpsilon1S->binMean(0), weight); + _histMeanMultiUpsilon1S->fill(_histMeanMultiUpsilon1S->binMean(0), weight); break; case 20223: - _hist91MeanMultiF1_1285->fill(_hist91MeanMultiF1_1285->binMean(0), weight); + _histMeanMultiF1_1285->fill(_histMeanMultiF1_1285->binMean(0), weight); break; case 20333: - _hist91MeanMultiF1_1420->fill(_hist91MeanMultiF1_1420->binMean(0), weight); + _histMeanMultiF1_1420->fill(_histMeanMultiF1_1420->binMean(0), weight); break; case 445: - _hist91MeanMultiChi_c1_3510->fill(_hist91MeanMultiChi_c1_3510->binMean(0), weight); + _histMeanMultiChi_c1_3510->fill(_histMeanMultiChi_c1_3510->binMean(0), weight); break; case 225: - _hist91MeanMultiF2_1270->fill(_hist91MeanMultiF2_1270->binMean(0), weight); + _histMeanMultiF2_1270->fill(_histMeanMultiF2_1270->binMean(0), weight); break; case 335: - _hist91MeanMultiF2Prime1525->fill(_hist91MeanMultiF2Prime1525->binMean(0), weight); + _histMeanMultiF2Prime1525->fill(_histMeanMultiF2Prime1525->binMean(0), weight); break; case 315: - _hist91MeanMultiK2Star1430_0->fill(_hist91MeanMultiK2Star1430_0->binMean(0), weight); + _histMeanMultiK2Star1430_0->fill(_histMeanMultiK2Star1430_0->binMean(0), weight); break; case 515: case 525: case 535: - _hist91MeanMultiBStarStar->fill(_hist91MeanMultiBStarStar->binMean(0), weight); + _histMeanMultiBStarStar->fill(_histMeanMultiBStarStar->binMean(0), weight); break; case 10433: case 20433: - _hist91MeanMultiDs1Plus->fill(_hist91MeanMultiDs1Plus->binMean(0), weight); + _histMeanMultiDs1Plus->fill(_histMeanMultiDs1Plus->binMean(0), weight); break; case 435: - _hist91MeanMultiDs2Plus->fill(_hist91MeanMultiDs2Plus->binMean(0), weight); + _histMeanMultiDs2Plus->fill(_histMeanMultiDs2Plus->binMean(0), weight); break; case 2212: - _hist91MeanMultiP->fill(_hist91MeanMultiP->binMean(0), weight); + _histMeanMultiP->fill(_histMeanMultiP->binMean(0), weight); break; case 3122: - _hist91MeanMultiLambda->fill(_hist91MeanMultiLambda->binMean(0), weight); + _histMeanMultiLambda->fill(_histMeanMultiLambda->binMean(0), weight); break; case 3212: - _hist91MeanMultiSigma0->fill(_hist91MeanMultiSigma0->binMean(0), weight); + _histMeanMultiSigma0->fill(_histMeanMultiSigma0->binMean(0), weight); break; case 3112: - _hist91MeanMultiSigmaMinus->fill(_hist91MeanMultiSigmaMinus->binMean(0), weight); - _hist91MeanMultiSigmaPlusMinus->fill(_hist91MeanMultiSigmaPlusMinus->binMean(0), weight); + _histMeanMultiSigmaMinus->fill(_histMeanMultiSigmaMinus->binMean(0), weight); + _histMeanMultiSigmaPlusMinus->fill(_histMeanMultiSigmaPlusMinus->binMean(0), weight); break; case 3222: - _hist91MeanMultiSigmaPlus->fill(_hist91MeanMultiSigmaPlus->binMean(0), weight); - _hist91MeanMultiSigmaPlusMinus->fill(_hist91MeanMultiSigmaPlusMinus->binMean(0), weight); + _histMeanMultiSigmaPlus->fill(_histMeanMultiSigmaPlus->binMean(0), weight); + _histMeanMultiSigmaPlusMinus->fill(_histMeanMultiSigmaPlusMinus->binMean(0), weight); break; case 3312: - _hist91MeanMultiXiMinus->fill(_hist91MeanMultiXiMinus->binMean(0), weight); + _histMeanMultiXiMinus->fill(_histMeanMultiXiMinus->binMean(0), weight); break; case 2224: - _hist91MeanMultiDelta1232PlusPlus->fill(_hist91MeanMultiDelta1232PlusPlus->binMean(0), weight); + _histMeanMultiDelta1232PlusPlus->fill(_histMeanMultiDelta1232PlusPlus->binMean(0), weight); break; case 3114: - _hist91MeanMultiSigma1385Minus->fill(_hist91MeanMultiSigma1385Minus->binMean(0), weight); - _hist91MeanMultiSigma1385PlusMinus->fill(_hist91MeanMultiSigma1385PlusMinus->binMean(0), weight); + _histMeanMultiSigma1385Minus->fill(_histMeanMultiSigma1385Minus->binMean(0), weight); + _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->binMean(0), weight); break; case 3224: - _hist91MeanMultiSigma1385Plus->fill(_hist91MeanMultiSigma1385Plus->binMean(0), weight); - _hist91MeanMultiSigma1385PlusMinus->fill(_hist91MeanMultiSigma1385PlusMinus->binMean(0), weight); + _histMeanMultiSigma1385Plus->fill(_histMeanMultiSigma1385Plus->binMean(0), weight); + _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->binMean(0), weight); break; case 3324: - _hist91MeanMultiXi1530_0->fill(_hist91MeanMultiXi1530_0->binMean(0), weight); + _histMeanMultiXi1530_0->fill(_histMeanMultiXi1530_0->binMean(0), weight); break; case 3334: - _hist91MeanMultiOmegaMinus->fill(_hist91MeanMultiOmegaMinus->binMean(0), weight); + _histMeanMultiOmegaMinus->fill(_histMeanMultiOmegaMinus->binMean(0), weight); break; case 4122: - _hist91MeanMultiLambda_c_Plus->fill(_hist91MeanMultiLambda_c_Plus->binMean(0), weight); + _histMeanMultiLambda_c_Plus->fill(_histMeanMultiLambda_c_Plus->binMean(0), weight); break; case 5122: - _hist91MeanMultiLambda_b_0->fill(_hist91MeanMultiLambda_b_0->binMean(0), weight); + _histMeanMultiLambda_b_0->fill(_histMeanMultiLambda_b_0->binMean(0), weight); break; case 3124: - _hist91MeanMultiLambda1520->fill(_hist91MeanMultiLambda1520->binMean(0), weight); + _histMeanMultiLambda1520->fill(_histMeanMultiLambda1520->binMean(0), weight); break; } } } - - - if (2*meanBeamMom/GeV >= 130 && 2*meanBeamMom/GeV <= 200) { + if (sqrtS()/GeV >= 130 && sqrtS()/GeV <= 200) { foreach (const Particle& p, ufs.particles()) { const PdgId id = abs(p.pdgId()); switch (id) { case 211: - _weightedTotalNumPiPlus165 += weight; + _weightedTotalNumPiPlus += weight; break; case 321: - _hist165MeanMultiKPlus->fill(_hist165MeanMultiKPlus->binMean(0), weight); + _histMeanMultiKPlus->fill(_histMeanMultiKPlus->binMean(0), weight); break; case 130: case 310: - _hist165MeanMultiK0->fill(_hist165MeanMultiK0->binMean(0), weight); + _histMeanMultiK0->fill(_histMeanMultiK0->binMean(0), weight); break; case 2212: - _hist165MeanMultiP->fill(_hist165MeanMultiP->binMean(0), weight); + _histMeanMultiP->fill(_histMeanMultiP->binMean(0), weight); break; case 3122: - _hist165MeanMultiLambda->fill(_hist165MeanMultiLambda->binMean(0), weight); + _histMeanMultiLambda->fill(_histMeanMultiLambda->binMean(0), weight); break; } } } + } void init() { - // Projections - addProjection(Beam(), "Beams"); addProjection(ChargedFinalState(), "FS"); addProjection(UnstableFinalState(), "UFS"); - // Histos - _hist10MeanMultiPi0 = bookHistogram1D( 2, 1, 1); - _hist10MeanMultiKPlus = bookHistogram1D( 3, 1, 1); - _hist10MeanMultiK0 = bookHistogram1D( 4, 1, 1); - _hist10MeanMultiEta = bookHistogram1D( 5, 1, 1); - _hist10MeanMultiEtaPrime = bookHistogram1D( 6, 1, 1); - _hist10MeanMultiDPlus = bookHistogram1D( 7, 1, 1); - _hist10MeanMultiD0 = bookHistogram1D( 8, 1, 1); - _hist10MeanMultiDPlus_s = bookHistogram1D( 9, 1, 1); - _hist10MeanMultiF0_980 = bookHistogram1D(13, 1, 1); - _hist10MeanMultiRho770_0 = bookHistogram1D(15, 1, 1); - _hist10MeanMultiOmega782 = bookHistogram1D(17, 1, 1); - _hist10MeanMultiKStar892Plus = bookHistogram1D(18, 1, 1); - _hist10MeanMultiKStar892_0 = bookHistogram1D(19, 1, 1); - _hist10MeanMultiPhi1020 = bookHistogram1D(20, 1, 1); - _hist10MeanMultiDStar2010Plus = bookHistogram1D(21, 1, 1); - _hist10MeanMultiDStar2007_0 = bookHistogram1D(22, 1, 1); - _hist10MeanMultiDStar_s2112Plus = bookHistogram1D(23, 1, 1); - _hist10MeanMultiJPsi1S = bookHistogram1D(25, 1, 1); - _hist10MeanMultiF2_1270 = bookHistogram1D(31, 1, 1); - _hist10MeanMultiP = bookHistogram1D(38, 1, 1); - _hist10MeanMultiLambda = bookHistogram1D(39, 1, 1); - _hist10MeanMultiSigma0 = bookHistogram1D(40, 1, 1); - _hist10MeanMultiXiMinus = bookHistogram1D(44, 1, 1); - _hist10MeanMultiDelta1232PlusPlus = bookHistogram1D(45, 1, 1); - _hist10MeanMultiSigma1385Minus = bookHistogram1D(46, 1, 1); - _hist10MeanMultiSigma1385Plus = bookHistogram1D(47, 1, 1); - _hist10MeanMultiSigma1385PlusMinus = bookHistogram1D(48, 1, 1); - _hist10MeanMultiXi1530_0 = bookHistogram1D(49, 1, 1); - _hist10MeanMultiOmegaMinus = bookHistogram1D(50, 1, 1); - _hist10MeanMultiLambda_c_Plus = bookHistogram1D(51, 1, 1); - _hist10MeanMultiSigma_c_PlusPlus_0 = bookHistogram1D(53, 1, 1); - _hist10MeanMultiLambda1520 = bookHistogram1D(54, 1, 1); - - _hist32MeanMultiPi0 = bookHistogram1D( 2, 1, 2); - _hist32MeanMultiKPlus = bookHistogram1D( 3, 1, 2); - _hist32MeanMultiK0 = bookHistogram1D( 4, 1, 2); - _hist32MeanMultiEta = bookHistogram1D( 5, 1, 2); - _hist32MeanMultiEtaPrime = bookHistogram1D( 6, 1, 2); - _hist32MeanMultiDPlus = bookHistogram1D( 7, 1, 2); - _hist32MeanMultiD0 = bookHistogram1D( 8, 1, 2); - _hist32MeanMultiDPlus_s = bookHistogram1D( 9, 1, 2); - _hist32MeanMultiF0_980 = bookHistogram1D(13, 1, 2); - _hist32MeanMultiRho770_0 = bookHistogram1D(15, 1, 2); - _hist32MeanMultiKStar892Plus = bookHistogram1D(18, 1, 2); - _hist32MeanMultiKStar892_0 = bookHistogram1D(19, 1, 2); - _hist32MeanMultiPhi1020 = bookHistogram1D(20, 1, 2); - _hist32MeanMultiDStar2010Plus = bookHistogram1D(21, 1, 2); - _hist32MeanMultiDStar2007_0 = bookHistogram1D(22, 1, 2); - _hist32MeanMultiF2_1270 = bookHistogram1D(31, 1, 2); - _hist32MeanMultiK2Star1430Plus = bookHistogram1D(33, 1, 1); - _hist32MeanMultiK2Star1430_0 = bookHistogram1D(34, 1, 1); - _hist32MeanMultiP = bookHistogram1D(38, 1, 2); - _hist32MeanMultiLambda = bookHistogram1D(39, 1, 2); - _hist32MeanMultiXiMinus = bookHistogram1D(44, 1, 2); - _hist32MeanMultiSigma1385Minus = bookHistogram1D(46, 1, 2); - _hist32MeanMultiSigma1385Plus = bookHistogram1D(47, 1, 2); - _hist32MeanMultiSigma1385PlusMinus = bookHistogram1D(48, 1, 2); - _hist32MeanMultiOmegaMinus = bookHistogram1D(50, 1, 2); - _hist32MeanMultiLambda_c_Plus = bookHistogram1D(51, 1, 2); - - _hist91MeanMultiPi0 = bookHistogram1D( 2, 1, 3); - _hist91MeanMultiKPlus = bookHistogram1D( 3, 1, 3); - _hist91MeanMultiK0 = bookHistogram1D( 4, 1, 3); - _hist91MeanMultiEta = bookHistogram1D( 5, 1, 3); - _hist91MeanMultiEtaPrime = bookHistogram1D( 6, 1, 3); - _hist91MeanMultiDPlus = bookHistogram1D( 7, 1, 3); - _hist91MeanMultiD0 = bookHistogram1D( 8, 1, 3); - _hist91MeanMultiDPlus_s = bookHistogram1D( 9, 1, 3); - _hist91MeanMultiBPlus_B0_d = bookHistogram1D(10, 1, 1); - _hist91MeanMultiBPlus_u = bookHistogram1D(11, 1, 1); - _hist91MeanMultiB0_s = bookHistogram1D(12, 1, 1); - _hist91MeanMultiF0_980 = bookHistogram1D(13, 1, 3); - _hist91MeanMultiA0_980Plus = bookHistogram1D(14, 1, 1); - _hist91MeanMultiRho770_0 = bookHistogram1D(15, 1, 3); - _hist91MeanMultiRho770Plus = bookHistogram1D(16, 1, 1); - _hist91MeanMultiOmega782 = bookHistogram1D(17, 1, 2); - _hist91MeanMultiKStar892Plus = bookHistogram1D(18, 1, 3); - _hist91MeanMultiKStar892_0 = bookHistogram1D(19, 1, 3); - _hist91MeanMultiPhi1020 = bookHistogram1D(20, 1, 3); - _hist91MeanMultiDStar2010Plus = bookHistogram1D(21, 1, 3); - _hist91MeanMultiDStar_s2112Plus = bookHistogram1D(23, 1, 2); - _hist91MeanMultiBStar = bookHistogram1D(24, 1, 1); - _hist91MeanMultiJPsi1S = bookHistogram1D(25, 1, 2); - _hist91MeanMultiPsi2S = bookHistogram1D(26, 1, 1); - _hist91MeanMultiUpsilon1S = bookHistogram1D(27, 1, 1); - _hist91MeanMultiF1_1285 = bookHistogram1D(28, 1, 1); - _hist91MeanMultiF1_1420 = bookHistogram1D(29, 1, 1); - _hist91MeanMultiChi_c1_3510 = bookHistogram1D(30, 1, 1); - _hist91MeanMultiF2_1270 = bookHistogram1D(31, 1, 3); - _hist91MeanMultiF2Prime1525 = bookHistogram1D(32, 1, 1); - _hist91MeanMultiK2Star1430_0 = bookHistogram1D(34, 1, 2); - _hist91MeanMultiBStarStar = bookHistogram1D(35, 1, 1); - _hist91MeanMultiDs1Plus = bookHistogram1D(36, 1, 1); - _hist91MeanMultiDs2Plus = bookHistogram1D(37, 1, 1); - _hist91MeanMultiP = bookHistogram1D(38, 1, 3); - _hist91MeanMultiLambda = bookHistogram1D(39, 1, 3); - _hist91MeanMultiSigma0 = bookHistogram1D(40, 1, 2); - _hist91MeanMultiSigmaMinus = bookHistogram1D(41, 1, 1); - _hist91MeanMultiSigmaPlus = bookHistogram1D(42, 1, 1); - _hist91MeanMultiSigmaPlusMinus = bookHistogram1D(43, 1, 1); - _hist91MeanMultiXiMinus = bookHistogram1D(44, 1, 3); - _hist91MeanMultiDelta1232PlusPlus = bookHistogram1D(45, 1, 2); - _hist91MeanMultiSigma1385Minus = bookHistogram1D(46, 1, 3); - _hist91MeanMultiSigma1385Plus = bookHistogram1D(47, 1, 3); - _hist91MeanMultiSigma1385PlusMinus = bookHistogram1D(48, 1, 3); - _hist91MeanMultiXi1530_0 = bookHistogram1D(49, 1, 2); - _hist91MeanMultiOmegaMinus = bookHistogram1D(50, 1, 3); - _hist91MeanMultiLambda_c_Plus = bookHistogram1D(51, 1, 3); - _hist91MeanMultiLambda_b_0 = bookHistogram1D(52, 1, 1); - _hist91MeanMultiLambda1520 = bookHistogram1D(54, 1, 2); - - _hist165MeanMultiKPlus = bookHistogram1D( 3, 1, 4); - _hist165MeanMultiK0 = bookHistogram1D( 4, 1, 4); - _hist165MeanMultiP = bookHistogram1D(38, 1, 4); - _hist165MeanMultiLambda = bookHistogram1D(39, 1, 4); + if (sqrtS()/GeV >= 9.5 && sqrtS()/GeV <= 10.5) { + _histMeanMultiPi0 = bookHistogram1D( 2, 1, 1); + _histMeanMultiKPlus = bookHistogram1D( 3, 1, 1); + _histMeanMultiK0 = bookHistogram1D( 4, 1, 1); + _histMeanMultiEta = bookHistogram1D( 5, 1, 1); + _histMeanMultiEtaPrime = bookHistogram1D( 6, 1, 1); + _histMeanMultiDPlus = bookHistogram1D( 7, 1, 1); + _histMeanMultiD0 = bookHistogram1D( 8, 1, 1); + _histMeanMultiDPlus_s = bookHistogram1D( 9, 1, 1); + _histMeanMultiF0_980 = bookHistogram1D(13, 1, 1); + _histMeanMultiRho770_0 = bookHistogram1D(15, 1, 1); + _histMeanMultiOmega782 = bookHistogram1D(17, 1, 1); + _histMeanMultiKStar892Plus = bookHistogram1D(18, 1, 1); + _histMeanMultiKStar892_0 = bookHistogram1D(19, 1, 1); + _histMeanMultiPhi1020 = bookHistogram1D(20, 1, 1); + _histMeanMultiDStar2010Plus = bookHistogram1D(21, 1, 1); + _histMeanMultiDStar2007_0 = bookHistogram1D(22, 1, 1); + _histMeanMultiDStar_s2112Plus = bookHistogram1D(23, 1, 1); + _histMeanMultiJPsi1S = bookHistogram1D(25, 1, 1); + _histMeanMultiF2_1270 = bookHistogram1D(31, 1, 1); + _histMeanMultiP = bookHistogram1D(38, 1, 1); + _histMeanMultiLambda = bookHistogram1D(39, 1, 1); + _histMeanMultiSigma0 = bookHistogram1D(40, 1, 1); + _histMeanMultiXiMinus = bookHistogram1D(44, 1, 1); + _histMeanMultiDelta1232PlusPlus = bookHistogram1D(45, 1, 1); + _histMeanMultiSigma1385Minus = bookHistogram1D(46, 1, 1); + _histMeanMultiSigma1385Plus = bookHistogram1D(47, 1, 1); + _histMeanMultiSigma1385PlusMinus = bookHistogram1D(48, 1, 1); + _histMeanMultiXi1530_0 = bookHistogram1D(49, 1, 1); + _histMeanMultiOmegaMinus = bookHistogram1D(50, 1, 1); + _histMeanMultiLambda_c_Plus = bookHistogram1D(51, 1, 1); + _histMeanMultiSigma_c_PlusPlus_0 = bookHistogram1D(53, 1, 1); + _histMeanMultiLambda1520 = bookHistogram1D(54, 1, 1); + } + + if (sqrtS()/GeV >= 29 && sqrtS()/GeV <= 35) { + _histMeanMultiPi0 = bookHistogram1D( 2, 1, 2); + _histMeanMultiKPlus = bookHistogram1D( 3, 1, 2); + _histMeanMultiK0 = bookHistogram1D( 4, 1, 2); + _histMeanMultiEta = bookHistogram1D( 5, 1, 2); + _histMeanMultiEtaPrime = bookHistogram1D( 6, 1, 2); + _histMeanMultiDPlus = bookHistogram1D( 7, 1, 2); + _histMeanMultiD0 = bookHistogram1D( 8, 1, 2); + _histMeanMultiDPlus_s = bookHistogram1D( 9, 1, 2); + _histMeanMultiF0_980 = bookHistogram1D(13, 1, 2); + _histMeanMultiRho770_0 = bookHistogram1D(15, 1, 2); + _histMeanMultiKStar892Plus = bookHistogram1D(18, 1, 2); + _histMeanMultiKStar892_0 = bookHistogram1D(19, 1, 2); + _histMeanMultiPhi1020 = bookHistogram1D(20, 1, 2); + _histMeanMultiDStar2010Plus = bookHistogram1D(21, 1, 2); + _histMeanMultiDStar2007_0 = bookHistogram1D(22, 1, 2); + _histMeanMultiF2_1270 = bookHistogram1D(31, 1, 2); + _histMeanMultiK2Star1430Plus = bookHistogram1D(33, 1, 1); + _histMeanMultiK2Star1430_0 = bookHistogram1D(34, 1, 1); + _histMeanMultiP = bookHistogram1D(38, 1, 2); + _histMeanMultiLambda = bookHistogram1D(39, 1, 2); + _histMeanMultiXiMinus = bookHistogram1D(44, 1, 2); + _histMeanMultiSigma1385Minus = bookHistogram1D(46, 1, 2); + _histMeanMultiSigma1385Plus = bookHistogram1D(47, 1, 2); + _histMeanMultiSigma1385PlusMinus = bookHistogram1D(48, 1, 2); + _histMeanMultiOmegaMinus = bookHistogram1D(50, 1, 2); + _histMeanMultiLambda_c_Plus = bookHistogram1D(51, 1, 2); + } + + if (sqrtS()/GeV >= 89.5 && sqrtS()/GeV <= 91.8) { + _histMeanMultiPi0 = bookHistogram1D( 2, 1, 3); + _histMeanMultiKPlus = bookHistogram1D( 3, 1, 3); + _histMeanMultiK0 = bookHistogram1D( 4, 1, 3); + _histMeanMultiEta = bookHistogram1D( 5, 1, 3); + _histMeanMultiEtaPrime = bookHistogram1D( 6, 1, 3); + _histMeanMultiDPlus = bookHistogram1D( 7, 1, 3); + _histMeanMultiD0 = bookHistogram1D( 8, 1, 3); + _histMeanMultiDPlus_s = bookHistogram1D( 9, 1, 3); + _histMeanMultiBPlus_B0_d = bookHistogram1D(10, 1, 1); + _histMeanMultiBPlus_u = bookHistogram1D(11, 1, 1); + _histMeanMultiB0_s = bookHistogram1D(12, 1, 1); + _histMeanMultiF0_980 = bookHistogram1D(13, 1, 3); + _histMeanMultiA0_980Plus = bookHistogram1D(14, 1, 1); + _histMeanMultiRho770_0 = bookHistogram1D(15, 1, 3); + _histMeanMultiRho770Plus = bookHistogram1D(16, 1, 1); + _histMeanMultiOmega782 = bookHistogram1D(17, 1, 2); + _histMeanMultiKStar892Plus = bookHistogram1D(18, 1, 3); + _histMeanMultiKStar892_0 = bookHistogram1D(19, 1, 3); + _histMeanMultiPhi1020 = bookHistogram1D(20, 1, 3); + _histMeanMultiDStar2010Plus = bookHistogram1D(21, 1, 3); + _histMeanMultiDStar_s2112Plus = bookHistogram1D(23, 1, 2); + _histMeanMultiBStar = bookHistogram1D(24, 1, 1); + _histMeanMultiJPsi1S = bookHistogram1D(25, 1, 2); + _histMeanMultiPsi2S = bookHistogram1D(26, 1, 1); + _histMeanMultiUpsilon1S = bookHistogram1D(27, 1, 1); + _histMeanMultiF1_1285 = bookHistogram1D(28, 1, 1); + _histMeanMultiF1_1420 = bookHistogram1D(29, 1, 1); + _histMeanMultiChi_c1_3510 = bookHistogram1D(30, 1, 1); + _histMeanMultiF2_1270 = bookHistogram1D(31, 1, 3); + _histMeanMultiF2Prime1525 = bookHistogram1D(32, 1, 1); + _histMeanMultiK2Star1430_0 = bookHistogram1D(34, 1, 2); + _histMeanMultiBStarStar = bookHistogram1D(35, 1, 1); + _histMeanMultiDs1Plus = bookHistogram1D(36, 1, 1); + _histMeanMultiDs2Plus = bookHistogram1D(37, 1, 1); + _histMeanMultiP = bookHistogram1D(38, 1, 3); + _histMeanMultiLambda = bookHistogram1D(39, 1, 3); + _histMeanMultiSigma0 = bookHistogram1D(40, 1, 2); + _histMeanMultiSigmaMinus = bookHistogram1D(41, 1, 1); + _histMeanMultiSigmaPlus = bookHistogram1D(42, 1, 1); + _histMeanMultiSigmaPlusMinus = bookHistogram1D(43, 1, 1); + _histMeanMultiXiMinus = bookHistogram1D(44, 1, 3); + _histMeanMultiDelta1232PlusPlus = bookHistogram1D(45, 1, 2); + _histMeanMultiSigma1385Minus = bookHistogram1D(46, 1, 3); + _histMeanMultiSigma1385Plus = bookHistogram1D(47, 1, 3); + _histMeanMultiSigma1385PlusMinus = bookHistogram1D(48, 1, 3); + _histMeanMultiXi1530_0 = bookHistogram1D(49, 1, 2); + _histMeanMultiOmegaMinus = bookHistogram1D(50, 1, 3); + _histMeanMultiLambda_c_Plus = bookHistogram1D(51, 1, 3); + _histMeanMultiLambda_b_0 = bookHistogram1D(52, 1, 1); + _histMeanMultiLambda1520 = bookHistogram1D(54, 1, 2); + } + + if (sqrtS()/GeV >= 130 && sqrtS()/GeV <= 200) { + _histMeanMultiKPlus = bookHistogram1D( 3, 1, 4); + _histMeanMultiK0 = bookHistogram1D( 4, 1, 4); + _histMeanMultiP = bookHistogram1D(38, 1, 4); + _histMeanMultiLambda = bookHistogram1D(39, 1, 4); + } } // Finalize void finalize() { - scale(_hist10MeanMultiPi0 , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiKPlus , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiK0 , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiEta , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiEtaPrime , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiDPlus , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiD0 , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiDPlus_s , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiF0_980 , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiRho770_0 , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiOmega782 , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiKStar892Plus , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiKStar892_0 , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiPhi1020 , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiDStar2010Plus , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiDStar2007_0 , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiDStar_s2112Plus , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiJPsi1S , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiF2_1270 , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiP , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiLambda , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiSigma0 , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiXiMinus , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiDelta1232PlusPlus , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiSigma1385Minus , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiSigma1385Plus , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiSigma1385PlusMinus, 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiXi1530_0 , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiOmegaMinus , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiLambda_c_Plus , 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiSigma_c_PlusPlus_0, 1.0/_weightedTotalNumPiPlus10); - scale(_hist10MeanMultiLambda1520 , 1.0/_weightedTotalNumPiPlus10); - - scale(_hist32MeanMultiPi0 , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiKPlus , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiK0 , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiEta , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiEtaPrime , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiDPlus , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiD0 , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiDPlus_s , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiF0_980 , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiRho770_0 , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiKStar892Plus , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiKStar892_0 , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiPhi1020 , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiDStar2010Plus , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiDStar2007_0 , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiF2_1270 , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiK2Star1430Plus , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiK2Star1430_0 , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiP , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiLambda , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiXiMinus , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiSigma1385Minus , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiSigma1385Plus , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiSigma1385PlusMinus, 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiOmegaMinus , 1.0/_weightedTotalNumPiPlus32); - scale(_hist32MeanMultiLambda_c_Plus , 1.0/_weightedTotalNumPiPlus32); - - scale(_hist91MeanMultiPi0 , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiKPlus , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiK0 , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiEta , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiEtaPrime , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiDPlus , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiD0 , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiDPlus_s , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiBPlus_B0_d , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiBPlus_u , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiB0_s , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiF0_980 , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiA0_980Plus , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiRho770_0 , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiRho770Plus , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiOmega782 , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiKStar892Plus , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiKStar892_0 , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiPhi1020 , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiDStar2010Plus , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiDStar_s2112Plus , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiBStar , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiJPsi1S , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiPsi2S , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiUpsilon1S , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiF1_1285 , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiF1_1420 , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiChi_c1_3510 , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiF2_1270 , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiF2Prime1525 , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiK2Star1430_0 , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiBStarStar , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiDs1Plus , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiDs2Plus , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiP , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiLambda , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiSigma0 , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiSigmaMinus , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiSigmaPlus , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiSigmaPlusMinus , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiXiMinus , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiDelta1232PlusPlus , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiSigma1385Minus , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiSigma1385Plus , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiSigma1385PlusMinus, 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiXi1530_0 , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiOmegaMinus , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiLambda_c_Plus , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiLambda_b_0 , 1.0/_weightedTotalNumPiPlus91); - scale(_hist91MeanMultiLambda1520 , 1.0/_weightedTotalNumPiPlus91); - - scale(_hist165MeanMultiKPlus , 1.0/_weightedTotalNumPiPlus165); - scale(_hist165MeanMultiK0 , 1.0/_weightedTotalNumPiPlus165); - scale(_hist165MeanMultiP , 1.0/_weightedTotalNumPiPlus165); - scale(_hist165MeanMultiLambda , 1.0/_weightedTotalNumPiPlus165); + 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 , 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(_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(_histMeanMultiF2_1270 , 1.0/_weightedTotalNumPiPlus); + scale(_histMeanMultiK2Star1430Plus , 1.0/_weightedTotalNumPiPlus); + scale(_histMeanMultiK2Star1430_0 , 1.0/_weightedTotalNumPiPlus); + scale(_histMeanMultiP , 1.0/_weightedTotalNumPiPlus); + scale(_histMeanMultiLambda , 1.0/_weightedTotalNumPiPlus); + scale(_histMeanMultiXiMinus , 1.0/_weightedTotalNumPiPlus); + scale(_histMeanMultiSigma1385Minus , 1.0/_weightedTotalNumPiPlus); + scale(_histMeanMultiSigma1385Plus , 1.0/_weightedTotalNumPiPlus); + scale(_histMeanMultiSigma1385PlusMinus, 1.0/_weightedTotalNumPiPlus); + scale(_histMeanMultiOmegaMinus , 1.0/_weightedTotalNumPiPlus); + scale(_histMeanMultiLambda_c_Plus , 1.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 , 1.0/_weightedTotalNumPiPlus); + scale(_histMeanMultiK0 , 1.0/_weightedTotalNumPiPlus); + scale(_histMeanMultiP , 1.0/_weightedTotalNumPiPlus); + scale(_histMeanMultiLambda , 1.0/_weightedTotalNumPiPlus); + } } //@} @@ -698,126 +699,61 @@ private: - double _weightedTotalNumPiPlus10; - double _weightedTotalNumPiPlus32; - double _weightedTotalNumPiPlus91; - double _weightedTotalNumPiPlus165; - - AIDA::IHistogram1D *_hist10MeanMultiPi0; - AIDA::IHistogram1D *_hist10MeanMultiKPlus; - AIDA::IHistogram1D *_hist10MeanMultiK0; - AIDA::IHistogram1D *_hist10MeanMultiEta; - AIDA::IHistogram1D *_hist10MeanMultiEtaPrime; - AIDA::IHistogram1D *_hist10MeanMultiDPlus; - AIDA::IHistogram1D *_hist10MeanMultiD0; - AIDA::IHistogram1D *_hist10MeanMultiDPlus_s; - AIDA::IHistogram1D *_hist10MeanMultiF0_980; - AIDA::IHistogram1D *_hist10MeanMultiRho770_0; - AIDA::IHistogram1D *_hist10MeanMultiOmega782; - AIDA::IHistogram1D *_hist10MeanMultiKStar892Plus; - AIDA::IHistogram1D *_hist10MeanMultiKStar892_0; - AIDA::IHistogram1D *_hist10MeanMultiPhi1020; - AIDA::IHistogram1D *_hist10MeanMultiDStar2010Plus; - AIDA::IHistogram1D *_hist10MeanMultiDStar2007_0; - AIDA::IHistogram1D *_hist10MeanMultiDStar_s2112Plus; - AIDA::IHistogram1D *_hist10MeanMultiJPsi1S; - AIDA::IHistogram1D *_hist10MeanMultiF2_1270; - AIDA::IHistogram1D *_hist10MeanMultiP; - AIDA::IHistogram1D *_hist10MeanMultiLambda; - AIDA::IHistogram1D *_hist10MeanMultiSigma0; - AIDA::IHistogram1D *_hist10MeanMultiXiMinus; - AIDA::IHistogram1D *_hist10MeanMultiDelta1232PlusPlus; - AIDA::IHistogram1D *_hist10MeanMultiSigma1385Minus; - AIDA::IHistogram1D *_hist10MeanMultiSigma1385Plus; - AIDA::IHistogram1D *_hist10MeanMultiSigma1385PlusMinus; - AIDA::IHistogram1D *_hist10MeanMultiXi1530_0; - AIDA::IHistogram1D *_hist10MeanMultiOmegaMinus; - AIDA::IHistogram1D *_hist10MeanMultiLambda_c_Plus; - AIDA::IHistogram1D *_hist10MeanMultiSigma_c_PlusPlus_0; - AIDA::IHistogram1D *_hist10MeanMultiLambda1520; - - AIDA::IHistogram1D *_hist32MeanMultiPi0; - AIDA::IHistogram1D *_hist32MeanMultiKPlus; - AIDA::IHistogram1D *_hist32MeanMultiK0; - AIDA::IHistogram1D *_hist32MeanMultiEta; - AIDA::IHistogram1D *_hist32MeanMultiEtaPrime; - AIDA::IHistogram1D *_hist32MeanMultiDPlus; - AIDA::IHistogram1D *_hist32MeanMultiD0; - AIDA::IHistogram1D *_hist32MeanMultiDPlus_s; - AIDA::IHistogram1D *_hist32MeanMultiF0_980; - AIDA::IHistogram1D *_hist32MeanMultiRho770_0; - AIDA::IHistogram1D *_hist32MeanMultiKStar892Plus; - AIDA::IHistogram1D *_hist32MeanMultiKStar892_0; - AIDA::IHistogram1D *_hist32MeanMultiPhi1020; - AIDA::IHistogram1D *_hist32MeanMultiDStar2010Plus; - AIDA::IHistogram1D *_hist32MeanMultiDStar2007_0; - AIDA::IHistogram1D *_hist32MeanMultiF2_1270; - AIDA::IHistogram1D *_hist32MeanMultiK2Star1430Plus; - AIDA::IHistogram1D *_hist32MeanMultiK2Star1430_0; - AIDA::IHistogram1D *_hist32MeanMultiP; - AIDA::IHistogram1D *_hist32MeanMultiLambda; - AIDA::IHistogram1D *_hist32MeanMultiXiMinus; - AIDA::IHistogram1D *_hist32MeanMultiSigma1385Minus; - AIDA::IHistogram1D *_hist32MeanMultiSigma1385Plus; - AIDA::IHistogram1D *_hist32MeanMultiSigma1385PlusMinus; - AIDA::IHistogram1D *_hist32MeanMultiOmegaMinus; - AIDA::IHistogram1D *_hist32MeanMultiLambda_c_Plus; - - AIDA::IHistogram1D *_hist91MeanMultiPi0; - AIDA::IHistogram1D *_hist91MeanMultiKPlus; - AIDA::IHistogram1D *_hist91MeanMultiK0; - AIDA::IHistogram1D *_hist91MeanMultiEta; - AIDA::IHistogram1D *_hist91MeanMultiEtaPrime; - AIDA::IHistogram1D *_hist91MeanMultiDPlus; - AIDA::IHistogram1D *_hist91MeanMultiD0; - AIDA::IHistogram1D *_hist91MeanMultiDPlus_s; - AIDA::IHistogram1D *_hist91MeanMultiBPlus_B0_d; - AIDA::IHistogram1D *_hist91MeanMultiBPlus_u; - AIDA::IHistogram1D *_hist91MeanMultiB0_s; - AIDA::IHistogram1D *_hist91MeanMultiF0_980; - AIDA::IHistogram1D *_hist91MeanMultiA0_980Plus; - AIDA::IHistogram1D *_hist91MeanMultiRho770_0; - AIDA::IHistogram1D *_hist91MeanMultiRho770Plus; - AIDA::IHistogram1D *_hist91MeanMultiOmega782; - AIDA::IHistogram1D *_hist91MeanMultiKStar892Plus; - AIDA::IHistogram1D *_hist91MeanMultiKStar892_0; - AIDA::IHistogram1D *_hist91MeanMultiPhi1020; - AIDA::IHistogram1D *_hist91MeanMultiDStar2010Plus; - AIDA::IHistogram1D *_hist91MeanMultiDStar_s2112Plus; - AIDA::IHistogram1D *_hist91MeanMultiBStar; - AIDA::IHistogram1D *_hist91MeanMultiJPsi1S; - AIDA::IHistogram1D *_hist91MeanMultiPsi2S; - AIDA::IHistogram1D *_hist91MeanMultiUpsilon1S; - AIDA::IHistogram1D *_hist91MeanMultiF1_1285; - AIDA::IHistogram1D *_hist91MeanMultiF1_1420; - AIDA::IHistogram1D *_hist91MeanMultiChi_c1_3510; - AIDA::IHistogram1D *_hist91MeanMultiF2_1270; - AIDA::IHistogram1D *_hist91MeanMultiF2Prime1525; - AIDA::IHistogram1D *_hist91MeanMultiK2Star1430_0; - AIDA::IHistogram1D *_hist91MeanMultiBStarStar; - AIDA::IHistogram1D *_hist91MeanMultiDs1Plus; - AIDA::IHistogram1D *_hist91MeanMultiDs2Plus; - AIDA::IHistogram1D *_hist91MeanMultiP; - AIDA::IHistogram1D *_hist91MeanMultiLambda; - AIDA::IHistogram1D *_hist91MeanMultiSigma0; - AIDA::IHistogram1D *_hist91MeanMultiSigmaMinus; - AIDA::IHistogram1D *_hist91MeanMultiSigmaPlus; - AIDA::IHistogram1D *_hist91MeanMultiSigmaPlusMinus; - AIDA::IHistogram1D *_hist91MeanMultiXiMinus; - AIDA::IHistogram1D *_hist91MeanMultiDelta1232PlusPlus; - AIDA::IHistogram1D *_hist91MeanMultiSigma1385Minus; - AIDA::IHistogram1D *_hist91MeanMultiSigma1385Plus; - AIDA::IHistogram1D *_hist91MeanMultiSigma1385PlusMinus; - AIDA::IHistogram1D *_hist91MeanMultiXi1530_0; - AIDA::IHistogram1D *_hist91MeanMultiOmegaMinus; - AIDA::IHistogram1D *_hist91MeanMultiLambda_c_Plus; - AIDA::IHistogram1D *_hist91MeanMultiLambda_b_0; - AIDA::IHistogram1D *_hist91MeanMultiLambda1520; - - AIDA::IHistogram1D *_hist165MeanMultiKPlus; - AIDA::IHistogram1D *_hist165MeanMultiK0; - AIDA::IHistogram1D *_hist165MeanMultiP; - AIDA::IHistogram1D *_hist165MeanMultiLambda; + double _weightedTotalNumPiPlus; + + AIDA::IHistogram1D *_histMeanMultiPi0; + AIDA::IHistogram1D *_histMeanMultiKPlus; + AIDA::IHistogram1D *_histMeanMultiK0; + AIDA::IHistogram1D *_histMeanMultiEta; + AIDA::IHistogram1D *_histMeanMultiEtaPrime; + AIDA::IHistogram1D *_histMeanMultiDPlus; + AIDA::IHistogram1D *_histMeanMultiD0; + AIDA::IHistogram1D *_histMeanMultiDPlus_s; + AIDA::IHistogram1D *_histMeanMultiBPlus_B0_d; + AIDA::IHistogram1D *_histMeanMultiBPlus_u; + AIDA::IHistogram1D *_histMeanMultiB0_s; + AIDA::IHistogram1D *_histMeanMultiF0_980; + AIDA::IHistogram1D *_histMeanMultiA0_980Plus; + AIDA::IHistogram1D *_histMeanMultiRho770_0; + AIDA::IHistogram1D *_histMeanMultiRho770Plus; + AIDA::IHistogram1D *_histMeanMultiOmega782; + AIDA::IHistogram1D *_histMeanMultiKStar892Plus; + AIDA::IHistogram1D *_histMeanMultiKStar892_0; + AIDA::IHistogram1D *_histMeanMultiPhi1020; + AIDA::IHistogram1D *_histMeanMultiDStar2010Plus; + AIDA::IHistogram1D *_histMeanMultiDStar2007_0; + AIDA::IHistogram1D *_histMeanMultiDStar_s2112Plus; + AIDA::IHistogram1D *_histMeanMultiBStar; + AIDA::IHistogram1D *_histMeanMultiJPsi1S; + AIDA::IHistogram1D *_histMeanMultiPsi2S; + AIDA::IHistogram1D *_histMeanMultiUpsilon1S; + AIDA::IHistogram1D *_histMeanMultiF1_1285; + AIDA::IHistogram1D *_histMeanMultiF1_1420; + AIDA::IHistogram1D *_histMeanMultiChi_c1_3510; + AIDA::IHistogram1D *_histMeanMultiF2_1270; + AIDA::IHistogram1D *_histMeanMultiF2Prime1525; + AIDA::IHistogram1D *_histMeanMultiK2Star1430Plus; + AIDA::IHistogram1D *_histMeanMultiK2Star1430_0; + AIDA::IHistogram1D *_histMeanMultiBStarStar; + AIDA::IHistogram1D *_histMeanMultiDs1Plus; + AIDA::IHistogram1D *_histMeanMultiDs2Plus; + AIDA::IHistogram1D *_histMeanMultiP; + AIDA::IHistogram1D *_histMeanMultiLambda; + AIDA::IHistogram1D *_histMeanMultiSigma0; + AIDA::IHistogram1D *_histMeanMultiSigmaMinus; + AIDA::IHistogram1D *_histMeanMultiSigmaPlus; + AIDA::IHistogram1D *_histMeanMultiSigmaPlusMinus; + AIDA::IHistogram1D *_histMeanMultiXiMinus; + AIDA::IHistogram1D *_histMeanMultiDelta1232PlusPlus; + AIDA::IHistogram1D *_histMeanMultiSigma1385Minus; + AIDA::IHistogram1D *_histMeanMultiSigma1385Plus; + AIDA::IHistogram1D *_histMeanMultiSigma1385PlusMinus; + AIDA::IHistogram1D *_histMeanMultiXi1530_0; + AIDA::IHistogram1D *_histMeanMultiOmegaMinus; + AIDA::IHistogram1D *_histMeanMultiLambda_c_Plus; + AIDA::IHistogram1D *_histMeanMultiLambda_b_0; + AIDA::IHistogram1D *_histMeanMultiSigma_c_PlusPlus_0; + AIDA::IHistogram1D *_histMeanMultiLambda1520; //@}
More information about the Rivet-svn mailing list |