rivet is hosted by Hepforge, IPPP Durham

Rivet analyses reference

PDG_HADRON_MULTIPLICITIES

Hadron multiplicities in hadronic $e^+e^-$ events
Experiment: PDG (Various)
Spires ID: 7857373
Status: VALIDATED
Authors:
  • Hendrik Hoeth
References:
  • Phys. Lett. B, 667, 1 (2008)
Beams: e+ e-
Beam energies: (5.0, 5.0); (17.5, 17.5); (45.6, 45.6); (88.5, 88.5) GeV
Run details:
  • Hadronic events in $e^+ e^-$ collisions. Beam energy must be specified as analysis option "ENERGY" when rivet-merging samples.

Hadron multiplicities in hadronic $e^+e^-$ events, taken from Review of Particle Properties 2008, table 40.1, page 355. Average hadron multiplicities per hadronic $e^+e^-$ annihilation event at $\sqrt{s} \approx {}$ 10, 29--35, 91, and 130--200 GeV. The numbers are averages from various experiments. Correlations of the systematic uncertainties were considered for the calculation of the averages. Beam energy must be specified as analysis option "ENERGY" when rivet-merging samples.

Source code: PDG_HADRON_MULTIPLICITIES.cc
  1// -*- C++ -*-
  2#include "Rivet/Analysis.hh"
  3#include "Rivet/Projections/FinalState.hh"
  4#include "Rivet/Projections/ChargedFinalState.hh"
  5#include "Rivet/Projections/UnstableParticles.hh"
  6
  7namespace Rivet {
  8
  9
 10  /// @brief Implementation of PDG hadron multiplicities
 11  /// @author Hendrik Hoeth
 12  class PDG_HADRON_MULTIPLICITIES : public Analysis {
 13  public:
 14
 15    /// Constructor
 16    PDG_HADRON_MULTIPLICITIES() : Analysis("PDG_HADRON_MULTIPLICITIES")
 17    {
 18    }
 19
 20
 21    /// @name Analysis methods
 22    /// @{
 23
 24    void analyze(const Event& e) {
 25      // First, veto on leptonic events by requiring at least 4 charged FS particles
 26      const FinalState& fs = apply<FinalState>(e, "FS");
 27      const size_t numParticles = fs.particles().size();
 28
 29      // Even if we only generate hadronic events, we still need a cut on numCharged >= 2.
 30      if (numParticles < 2) {
 31        MSG_DEBUG("Failed leptonic event cut");
 32        vetoEvent;
 33      }
 34      MSG_DEBUG("Passed leptonic event cut");
 35
 36      MSG_DEBUG("sqrt(s) = " << sqrtS()/GeV << " GeV");
 37
 38      // Final state of unstable particles to get particle spectra
 39      const UnstableParticles& ufs = apply<UnstableParticles>(e, "UFS");
 40
 41      if (sqrtS()/GeV >= 9.5 && sqrtS()/GeV <= 10.5) {
 42        for (const Particle& p : ufs.particles()) {
 43          const PdgId id = p.abspid();
 44          switch (id) {
 45             case 211:
 46                _histMeanMultiPiPlus->fill(_histMeanMultiPiPlus->bin(1).xMid());
 47                break;
 48             case 111:
 49                _histMeanMultiPi0->fill(_histMeanMultiPi0->bin(1).xMid());
 50                break;
 51             case 321:
 52                _histMeanMultiKPlus->fill(_histMeanMultiKPlus->bin(1).xMid());
 53                break;
 54             case 130:
 55             case 310:
 56                _histMeanMultiK0->fill(_histMeanMultiK0->bin(1).xMid());
 57                break;
 58             case 221:
 59                _histMeanMultiEta->fill(_histMeanMultiEta->bin(1).xMid());
 60                break;
 61             case 331:
 62                _histMeanMultiEtaPrime->fill(_histMeanMultiEtaPrime->bin(1).xMid());
 63                break;
 64             case 411:
 65                _histMeanMultiDPlus->fill(_histMeanMultiDPlus->bin(1).xMid());
 66                break;
 67             case 421:
 68                _histMeanMultiD0->fill(_histMeanMultiD0->bin(1).xMid());
 69                break;
 70             case 431:
 71                _histMeanMultiDPlus_s->fill(_histMeanMultiDPlus_s->bin(1).xMid());
 72                break;
 73             case 9010221:
 74                _histMeanMultiF0_980->fill(_histMeanMultiF0_980->bin(1).xMid());
 75                break;
 76             case 113:
 77                _histMeanMultiRho770_0->fill(_histMeanMultiRho770_0->bin(1).xMid());
 78                break;
 79             case 223:
 80                _histMeanMultiOmega782->fill(_histMeanMultiOmega782->bin(1).xMid());
 81                break;
 82             case 323:
 83                _histMeanMultiKStar892Plus->fill(_histMeanMultiKStar892Plus->bin(1).xMid());
 84                break;
 85             case 313:
 86                _histMeanMultiKStar892_0->fill(_histMeanMultiKStar892_0->bin(1).xMid());
 87                break;
 88             case 333:
 89                _histMeanMultiPhi1020->fill(_histMeanMultiPhi1020->bin(1).xMid());
 90                break;
 91             case 413:
 92                _histMeanMultiDStar2010Plus->fill(_histMeanMultiDStar2010Plus->bin(1).xMid());
 93                break;
 94             case 423:
 95                _histMeanMultiDStar2007_0->fill(_histMeanMultiDStar2007_0->bin(1).xMid());
 96                break;
 97             case 433:
 98                _histMeanMultiDStar_s2112Plus->fill(_histMeanMultiDStar_s2112Plus->bin(1).xMid());
 99                break;
100             case 443:
101                _histMeanMultiJPsi1S->fill(_histMeanMultiJPsi1S->bin(1).xMid());
102                break;
103             case 225:
104                _histMeanMultiF2_1270->fill(_histMeanMultiF2_1270->bin(1).xMid());
105                break;
106             case 2212:
107                _histMeanMultiP->fill(_histMeanMultiP->bin(1).xMid());
108                break;
109             case 3122:
110                _histMeanMultiLambda->fill(_histMeanMultiLambda->bin(1).xMid());
111                break;
112             case 3212:
113                _histMeanMultiSigma0->fill(_histMeanMultiSigma0->bin(1).xMid());
114                break;
115             case 3312:
116                _histMeanMultiXiMinus->fill(_histMeanMultiXiMinus->bin(1).xMid());
117                break;
118             case 2224:
119                _histMeanMultiDelta1232PlusPlus->fill(_histMeanMultiDelta1232PlusPlus->bin(1).xMid());
120                break;
121             case 3114:
122                _histMeanMultiSigma1385Minus->fill(_histMeanMultiSigma1385Minus->bin(1).xMid());
123                _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->bin(1).xMid());
124                break;
125             case 3224:
126                _histMeanMultiSigma1385Plus->fill(_histMeanMultiSigma1385Plus->bin(1).xMid());
127                _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->bin(1).xMid());
128                break;
129             case 3324:
130                _histMeanMultiXi1530_0->fill(_histMeanMultiXi1530_0->bin(1).xMid());
131                break;
132             case 3334:
133                _histMeanMultiOmegaMinus->fill(_histMeanMultiOmegaMinus->bin(1).xMid());
134                break;
135             case 4122:
136                _histMeanMultiLambda_c_Plus->fill(_histMeanMultiLambda_c_Plus->bin(1).xMid());
137                break;
138             case 4222:
139             case 4112:
140                _histMeanMultiSigma_c_PlusPlus_0->fill(_histMeanMultiSigma_c_PlusPlus_0->bin(1).xMid());
141                break;
142             case 3124:
143                _histMeanMultiLambda1520->fill(_histMeanMultiLambda1520->bin(1).xMid());
144                break;
145          }
146        }
147      }
148
149      if (sqrtS()/GeV >= 29 && sqrtS()/GeV <= 35) {
150        for (const Particle& p : ufs.particles()) {
151          const PdgId id = p.abspid();
152          switch (id) {
153             case 211:
154                _histMeanMultiPiPlus->fill(_histMeanMultiPiPlus->bin(1).xMid());
155                break;
156             case 111:
157                _histMeanMultiPi0->fill(_histMeanMultiPi0->bin(1).xMid());
158                break;
159             case 321:
160                _histMeanMultiKPlus->fill(_histMeanMultiKPlus->bin(1).xMid());
161                break;
162             case 130:
163             case 310:
164                _histMeanMultiK0->fill(_histMeanMultiK0->bin(1).xMid());
165                break;
166             case 221:
167                _histMeanMultiEta->fill(_histMeanMultiEta->bin(1).xMid());
168                break;
169             case 331:
170                _histMeanMultiEtaPrime->fill(_histMeanMultiEtaPrime->bin(1).xMid());
171                break;
172             case 411:
173                _histMeanMultiDPlus->fill(_histMeanMultiDPlus->bin(1).xMid());
174                break;
175             case 421:
176                _histMeanMultiD0->fill(_histMeanMultiD0->bin(1).xMid());
177                break;
178             case 431:
179                _histMeanMultiDPlus_s->fill(_histMeanMultiDPlus_s->bin(1).xMid());
180                break;
181             case 9010221:
182                _histMeanMultiF0_980->fill(_histMeanMultiF0_980->bin(1).xMid());
183                break;
184             case 113:
185                _histMeanMultiRho770_0->fill(_histMeanMultiRho770_0->bin(1).xMid());
186                break;
187             case 323:
188                _histMeanMultiKStar892Plus->fill(_histMeanMultiKStar892Plus->bin(1).xMid());
189                break;
190             case 313:
191                _histMeanMultiKStar892_0->fill(_histMeanMultiKStar892_0->bin(1).xMid());
192                break;
193             case 333:
194                _histMeanMultiPhi1020->fill(_histMeanMultiPhi1020->bin(1).xMid());
195                break;
196             case 413:
197                _histMeanMultiDStar2010Plus->fill(_histMeanMultiDStar2010Plus->bin(1).xMid());
198                break;
199             case 423:
200                _histMeanMultiDStar2007_0->fill(_histMeanMultiDStar2007_0->bin(1).xMid());
201                break;
202             case 225:
203                _histMeanMultiF2_1270->fill(_histMeanMultiF2_1270->bin(1).xMid());
204                break;
205             case 325:
206                _histMeanMultiK2Star1430Plus->fill(_histMeanMultiK2Star1430Plus->bin(1).xMid());
207                break;
208             case 315:
209                _histMeanMultiK2Star1430_0->fill(_histMeanMultiK2Star1430_0->bin(1).xMid());
210                break;
211             case 2212:
212                _histMeanMultiP->fill(_histMeanMultiP->bin(1).xMid());
213                break;
214             case 3122:
215                _histMeanMultiLambda->fill(_histMeanMultiLambda->bin(1).xMid());
216                break;
217             case 3312:
218                _histMeanMultiXiMinus->fill(_histMeanMultiXiMinus->bin(1).xMid());
219                break;
220             case 3114:
221                _histMeanMultiSigma1385Minus->fill(_histMeanMultiSigma1385Minus->bin(1).xMid());
222                _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->bin(1).xMid());
223                break;
224             case 3224:
225                _histMeanMultiSigma1385Plus->fill(_histMeanMultiSigma1385Plus->bin(1).xMid());
226                _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->bin(1).xMid());
227                break;
228             case 3334:
229                _histMeanMultiOmegaMinus->fill(_histMeanMultiOmegaMinus->bin(1).xMid());
230                break;
231             case 4122:
232                _histMeanMultiLambda_c_Plus->fill(_histMeanMultiLambda_c_Plus->bin(1).xMid());
233                break;
234          }
235        }
236      }
237
238      if (sqrtS()/GeV >= 89.5 && sqrtS()/GeV <= 91.8) {
239        for (const Particle& p : ufs.particles()) {
240          const PdgId id = p.abspid();
241          switch (id) {
242             case 211:
243                _histMeanMultiPiPlus->fill(_histMeanMultiPiPlus->bin(1).xMid());
244                break;
245             case 111:
246                _histMeanMultiPi0->fill(_histMeanMultiPi0->bin(1).xMid());
247                break;
248             case 321:
249                _histMeanMultiKPlus->fill(_histMeanMultiKPlus->bin(1).xMid());
250                break;
251             case 130:
252             case 310:
253                _histMeanMultiK0->fill(_histMeanMultiK0->bin(1).xMid());
254                break;
255             case 221:
256                _histMeanMultiEta->fill(_histMeanMultiEta->bin(1).xMid());
257                break;
258             case 331:
259                _histMeanMultiEtaPrime->fill(_histMeanMultiEtaPrime->bin(1).xMid());
260                break;
261             case 411:
262                _histMeanMultiDPlus->fill(_histMeanMultiDPlus->bin(1).xMid());
263                break;
264             case 421:
265                _histMeanMultiD0->fill(_histMeanMultiD0->bin(1).xMid());
266                break;
267             case 431:
268                _histMeanMultiDPlus_s->fill(_histMeanMultiDPlus_s->bin(1).xMid());
269                break;
270             case 511:
271                _histMeanMultiBPlus_B0_d->fill(_histMeanMultiBPlus_B0_d->bin(1).xMid());
272                break;
273             case 521:
274                _histMeanMultiBPlus_B0_d->fill(_histMeanMultiBPlus_B0_d->bin(1).xMid());
275                _histMeanMultiBPlus_u->fill(_histMeanMultiBPlus_u->bin(1).xMid());
276                break;
277             case 531:
278                _histMeanMultiB0_s->fill(_histMeanMultiB0_s->bin(1).xMid());
279                break;
280             case 9010221:
281                _histMeanMultiF0_980->fill(_histMeanMultiF0_980->bin(1).xMid());
282                break;
283             case 9000211:
284                _histMeanMultiA0_980Plus->fill(_histMeanMultiA0_980Plus->bin(1).xMid());
285                break;
286             case 113:
287                _histMeanMultiRho770_0->fill(_histMeanMultiRho770_0->bin(1).xMid());
288                break;
289             case 213:
290                _histMeanMultiRho770Plus->fill(_histMeanMultiRho770Plus->bin(1).xMid());
291                break;
292             case 223:
293                _histMeanMultiOmega782->fill(_histMeanMultiOmega782->bin(1).xMid());
294                break;
295             case 323:
296                _histMeanMultiKStar892Plus->fill(_histMeanMultiKStar892Plus->bin(1).xMid());
297                break;
298             case 313:
299                _histMeanMultiKStar892_0->fill(_histMeanMultiKStar892_0->bin(1).xMid());
300                break;
301             case 333:
302                _histMeanMultiPhi1020->fill(_histMeanMultiPhi1020->bin(1).xMid());
303                break;
304             case 413:
305                _histMeanMultiDStar2010Plus->fill(_histMeanMultiDStar2010Plus->bin(1).xMid());
306                break;
307             case 433:
308                _histMeanMultiDStar_s2112Plus->fill(_histMeanMultiDStar_s2112Plus->bin(1).xMid());
309                break;
310             case 513:
311             case 523:
312             case 533:
313                _histMeanMultiBStar->fill(_histMeanMultiBStar->bin(1).xMid());
314                break;
315             case 443:
316                _histMeanMultiJPsi1S->fill(_histMeanMultiJPsi1S->bin(1).xMid());
317                break;
318             case 100443:
319                _histMeanMultiPsi2S->fill(_histMeanMultiPsi2S->bin(1).xMid());
320                break;
321             case 553:
322                _histMeanMultiUpsilon1S->fill(_histMeanMultiUpsilon1S->bin(1).xMid());
323                break;
324             case 20223:
325                _histMeanMultiF1_1285->fill(_histMeanMultiF1_1285->bin(1).xMid());
326                break;
327             case 20333:
328                _histMeanMultiF1_1420->fill(_histMeanMultiF1_1420->bin(1).xMid());
329                break;
330             case 445:
331                _histMeanMultiChi_c1_3510->fill(_histMeanMultiChi_c1_3510->bin(1).xMid());
332                break;
333             case 225:
334                _histMeanMultiF2_1270->fill(_histMeanMultiF2_1270->bin(1).xMid());
335                break;
336             case 335:
337                _histMeanMultiF2Prime1525->fill(_histMeanMultiF2Prime1525->bin(1).xMid());
338                break;
339             case 315:
340                _histMeanMultiK2Star1430_0->fill(_histMeanMultiK2Star1430_0->bin(1).xMid());
341                break;
342             case 515:
343             case 525:
344             case 535:
345                _histMeanMultiBStarStar->fill(_histMeanMultiBStarStar->bin(1).xMid());
346                break;
347             case 10433:
348             case 20433:
349                _histMeanMultiDs1Plus->fill(_histMeanMultiDs1Plus->bin(1).xMid());
350                break;
351             case 435:
352                _histMeanMultiDs2Plus->fill(_histMeanMultiDs2Plus->bin(1).xMid());
353                break;
354             case 2212:
355                _histMeanMultiP->fill(_histMeanMultiP->bin(1).xMid());
356                break;
357             case 3122:
358                _histMeanMultiLambda->fill(_histMeanMultiLambda->bin(1).xMid());
359                break;
360             case 3212:
361                _histMeanMultiSigma0->fill(_histMeanMultiSigma0->bin(1).xMid());
362                break;
363             case 3112:
364                _histMeanMultiSigmaMinus->fill(_histMeanMultiSigmaMinus->bin(1).xMid());
365                _histMeanMultiSigmaPlusMinus->fill(_histMeanMultiSigmaPlusMinus->bin(1).xMid());
366                break;
367             case 3222:
368                _histMeanMultiSigmaPlus->fill(_histMeanMultiSigmaPlus->bin(1).xMid());
369                _histMeanMultiSigmaPlusMinus->fill(_histMeanMultiSigmaPlusMinus->bin(1).xMid());
370                break;
371             case 3312:
372                _histMeanMultiXiMinus->fill(_histMeanMultiXiMinus->bin(1).xMid());
373                break;
374             case 2224:
375                _histMeanMultiDelta1232PlusPlus->fill(_histMeanMultiDelta1232PlusPlus->bin(1).xMid());
376                break;
377             case 3114:
378                _histMeanMultiSigma1385Minus->fill(_histMeanMultiSigma1385Minus->bin(1).xMid());
379                _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->bin(1).xMid());
380                break;
381             case 3224:
382                _histMeanMultiSigma1385Plus->fill(_histMeanMultiSigma1385Plus->bin(1).xMid());
383                _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->bin(1).xMid());
384                break;
385             case 3324:
386                _histMeanMultiXi1530_0->fill(_histMeanMultiXi1530_0->bin(1).xMid());
387                break;
388             case 3334:
389                _histMeanMultiOmegaMinus->fill(_histMeanMultiOmegaMinus->bin(1).xMid());
390                break;
391             case 4122:
392                _histMeanMultiLambda_c_Plus->fill(_histMeanMultiLambda_c_Plus->bin(1).xMid());
393                break;
394             case 5122:
395                _histMeanMultiLambda_b_0->fill(_histMeanMultiLambda_b_0->bin(1).xMid());
396                break;
397             case 3124:
398                _histMeanMultiLambda1520->fill(_histMeanMultiLambda1520->bin(1).xMid());
399                break;
400          }
401        }
402      }
403
404      if (sqrtS()/GeV >= 130 && sqrtS()/GeV <= 200) {
405        for (const Particle& p : ufs.particles()) {
406          const PdgId id = p.abspid();
407          switch (id) {
408             case 211:
409                _histMeanMultiPiPlus->fill(_histMeanMultiPiPlus->bin(1).xMid());
410                break;
411             case 321:
412                _histMeanMultiKPlus->fill(_histMeanMultiKPlus->bin(1).xMid());
413                break;
414             case 130:
415             case 310:
416                _histMeanMultiK0->fill(_histMeanMultiK0->bin(1).xMid());
417                break;
418             case 2212:
419                _histMeanMultiP->fill(_histMeanMultiP->bin(1).xMid());
420                break;
421             case 3122:
422                _histMeanMultiLambda->fill(_histMeanMultiLambda->bin(1).xMid());
423                break;
424          }
425        }
426      }
427
428    }
429
430
431
432    void init() {
433      declare(ChargedFinalState(), "FS");
434      declare(UnstableParticles(), "UFS");
435
436      if (sqrtS()/GeV >= 9.5 && sqrtS()/GeV <= 10.5) {
437        book(_histMeanMultiPiPlus             , 1, 1, 1);
438        book(_histMeanMultiPi0                , 2, 1, 1);
439        book(_histMeanMultiKPlus              , 3, 1, 1);
440        book(_histMeanMultiK0                 , 4, 1, 1);
441        book(_histMeanMultiEta                , 5, 1, 1);
442        book(_histMeanMultiEtaPrime           , 6, 1, 1);
443        book(_histMeanMultiDPlus              , 7, 1, 1);
444        book(_histMeanMultiD0                 , 8, 1, 1);
445        book(_histMeanMultiDPlus_s            , 9, 1, 1);
446        book(_histMeanMultiF0_980             ,13, 1, 1);
447        book(_histMeanMultiRho770_0           ,15, 1, 1);
448        book(_histMeanMultiOmega782           ,17, 1, 1);
449        book(_histMeanMultiKStar892Plus       ,18, 1, 1);
450        book(_histMeanMultiKStar892_0         ,19, 1, 1);
451        book(_histMeanMultiPhi1020            ,20, 1, 1);
452        book(_histMeanMultiDStar2010Plus      ,21, 1, 1);
453        book(_histMeanMultiDStar2007_0        ,22, 1, 1);
454        book(_histMeanMultiDStar_s2112Plus    ,23, 1, 1);
455        book(_histMeanMultiJPsi1S             ,25, 1, 1);
456        book(_histMeanMultiF2_1270            ,31, 1, 1);
457        book(_histMeanMultiP                  ,38, 1, 1);
458        book(_histMeanMultiLambda             ,39, 1, 1);
459        book(_histMeanMultiSigma0             ,40, 1, 1);
460        book(_histMeanMultiXiMinus            ,44, 1, 1);
461        book(_histMeanMultiDelta1232PlusPlus  ,45, 1, 1);
462        book(_histMeanMultiSigma1385Minus     ,46, 1, 1);
463        book(_histMeanMultiSigma1385Plus      ,47, 1, 1);
464        book(_histMeanMultiSigma1385PlusMinus ,48, 1, 1);
465        book(_histMeanMultiXi1530_0           ,49, 1, 1);
466        book(_histMeanMultiOmegaMinus         ,50, 1, 1);
467        book(_histMeanMultiLambda_c_Plus      ,51, 1, 1);
468        book(_histMeanMultiSigma_c_PlusPlus_0 ,53, 1, 1);
469        book(_histMeanMultiLambda1520         ,54, 1, 1);
470      }
471
472      if (sqrtS()/GeV >= 29 && sqrtS()/GeV <= 35) {
473        book(_histMeanMultiPiPlus             , 1, 1, 2);
474        book(_histMeanMultiPi0                , 2, 1, 2);
475        book(_histMeanMultiKPlus              , 3, 1, 2);
476        book(_histMeanMultiK0                 , 4, 1, 2);
477        book(_histMeanMultiEta                , 5, 1, 2);
478        book(_histMeanMultiEtaPrime           , 6, 1, 2);
479        book(_histMeanMultiDPlus              , 7, 1, 2);
480        book(_histMeanMultiD0                 , 8, 1, 2);
481        book(_histMeanMultiDPlus_s            , 9, 1, 2);
482        book(_histMeanMultiF0_980             ,13, 1, 2);
483        book(_histMeanMultiRho770_0           ,15, 1, 2);
484        book(_histMeanMultiKStar892Plus       ,18, 1, 2);
485        book(_histMeanMultiKStar892_0         ,19, 1, 2);
486        book(_histMeanMultiPhi1020            ,20, 1, 2);
487        book(_histMeanMultiDStar2010Plus      ,21, 1, 2);
488        book(_histMeanMultiDStar2007_0        ,22, 1, 2);
489        book(_histMeanMultiF2_1270            ,31, 1, 2);
490        book(_histMeanMultiK2Star1430Plus     ,33, 1, 1);
491        book(_histMeanMultiK2Star1430_0       ,34, 1, 1);
492        book(_histMeanMultiP                  ,38, 1, 2);
493        book(_histMeanMultiLambda             ,39, 1, 2);
494        book(_histMeanMultiXiMinus            ,44, 1, 2);
495        book(_histMeanMultiSigma1385Minus     ,46, 1, 2);
496        book(_histMeanMultiSigma1385Plus      ,47, 1, 2);
497        book(_histMeanMultiSigma1385PlusMinus ,48, 1, 2);
498        book(_histMeanMultiOmegaMinus         ,50, 1, 2);
499        book(_histMeanMultiLambda_c_Plus      ,51, 1, 2);
500      }
501
502      if (sqrtS()/GeV >= 89.5 && sqrtS()/GeV <= 91.8) {
503        book(_histMeanMultiPiPlus             , 1, 1, 3);
504        book(_histMeanMultiPi0                , 2, 1, 3);
505        book(_histMeanMultiKPlus              , 3, 1, 3);
506        book(_histMeanMultiK0                 , 4, 1, 3);
507        book(_histMeanMultiEta                , 5, 1, 3);
508        book(_histMeanMultiEtaPrime           , 6, 1, 3);
509        book(_histMeanMultiDPlus              , 7, 1, 3);
510        book(_histMeanMultiD0                 , 8, 1, 3);
511        book(_histMeanMultiDPlus_s            , 9, 1, 3);
512        book(_histMeanMultiBPlus_B0_d         ,10, 1, 1);
513        book(_histMeanMultiBPlus_u            ,11, 1, 1);
514        book(_histMeanMultiB0_s               ,12, 1, 1);
515        book(_histMeanMultiF0_980             ,13, 1, 3);
516        book(_histMeanMultiA0_980Plus         ,14, 1, 1);
517        book(_histMeanMultiRho770_0           ,15, 1, 3);
518        book(_histMeanMultiRho770Plus         ,16, 1, 1);
519        book(_histMeanMultiOmega782           ,17, 1, 2);
520        book(_histMeanMultiKStar892Plus       ,18, 1, 3);
521        book(_histMeanMultiKStar892_0         ,19, 1, 3);
522        book(_histMeanMultiPhi1020            ,20, 1, 3);
523        book(_histMeanMultiDStar2010Plus      ,21, 1, 3);
524        book(_histMeanMultiDStar_s2112Plus    ,23, 1, 2);
525        book(_histMeanMultiBStar              ,24, 1, 1);
526        book(_histMeanMultiJPsi1S             ,25, 1, 2);
527        book(_histMeanMultiPsi2S              ,26, 1, 1);
528        book(_histMeanMultiUpsilon1S          ,27, 1, 1);
529        book(_histMeanMultiF1_1285            ,28, 1, 1);
530        book(_histMeanMultiF1_1420            ,29, 1, 1);
531        book(_histMeanMultiChi_c1_3510        ,30, 1, 1);
532        book(_histMeanMultiF2_1270            ,31, 1, 3);
533        book(_histMeanMultiF2Prime1525        ,32, 1, 1);
534        book(_histMeanMultiK2Star1430_0       ,34, 1, 2);
535        book(_histMeanMultiBStarStar          ,35, 1, 1);
536        book(_histMeanMultiDs1Plus            ,36, 1, 1);
537        book(_histMeanMultiDs2Plus            ,37, 1, 1);
538        book(_histMeanMultiP                  ,38, 1, 3);
539        book(_histMeanMultiLambda             ,39, 1, 3);
540        book(_histMeanMultiSigma0             ,40, 1, 2);
541        book(_histMeanMultiSigmaMinus         ,41, 1, 1);
542        book(_histMeanMultiSigmaPlus          ,42, 1, 1);
543        book(_histMeanMultiSigmaPlusMinus     ,43, 1, 1);
544        book(_histMeanMultiXiMinus            ,44, 1, 3);
545        book(_histMeanMultiDelta1232PlusPlus  ,45, 1, 2);
546        book(_histMeanMultiSigma1385Minus     ,46, 1, 3);
547        book(_histMeanMultiSigma1385Plus      ,47, 1, 3);
548        book(_histMeanMultiSigma1385PlusMinus ,48, 1, 3);
549        book(_histMeanMultiXi1530_0           ,49, 1, 2);
550        book(_histMeanMultiOmegaMinus         ,50, 1, 3);
551        book(_histMeanMultiLambda_c_Plus      ,51, 1, 3);
552        book(_histMeanMultiLambda_b_0         ,52, 1, 1);
553        book(_histMeanMultiLambda1520         ,54, 1, 2);
554      }
555
556      if (sqrtS()/GeV >= 130 && sqrtS()/GeV <= 200) {
557        book(_histMeanMultiPiPlus            , 1, 1, 4);
558        book(_histMeanMultiKPlus             , 3, 1, 4);
559        book(_histMeanMultiK0                , 4, 1, 4);
560        book(_histMeanMultiP                 ,38, 1, 4);
561        book(_histMeanMultiLambda            ,39, 1, 4);
562      }
563    }
564
565
566
567    // Finalize
568    void finalize() {
569      if (sqrtS()/GeV >= 9.5 && sqrtS()/GeV <= 10.5) {
570        scale(_histMeanMultiPiPlus            , 1.0/sumOfWeights());
571        scale(_histMeanMultiPi0               , 1.0/sumOfWeights());
572        scale(_histMeanMultiKPlus             , 1.0/sumOfWeights());
573        scale(_histMeanMultiK0                , 1.0/sumOfWeights());
574        scale(_histMeanMultiEta               , 1.0/sumOfWeights());
575        scale(_histMeanMultiEtaPrime          , 1.0/sumOfWeights());
576        scale(_histMeanMultiDPlus             , 1.0/sumOfWeights());
577        scale(_histMeanMultiD0                , 1.0/sumOfWeights());
578        scale(_histMeanMultiDPlus_s           , 1.0/sumOfWeights());
579        scale(_histMeanMultiF0_980            , 1.0/sumOfWeights());
580        scale(_histMeanMultiRho770_0          , 1.0/sumOfWeights());
581        scale(_histMeanMultiOmega782          , 1.0/sumOfWeights());
582        scale(_histMeanMultiKStar892Plus      , 1.0/sumOfWeights());
583        scale(_histMeanMultiKStar892_0        , 1.0/sumOfWeights());
584        scale(_histMeanMultiPhi1020           , 1.0/sumOfWeights());
585        scale(_histMeanMultiDStar2010Plus     , 1.0/sumOfWeights());
586        scale(_histMeanMultiDStar2007_0       , 1.0/sumOfWeights());
587        scale(_histMeanMultiDStar_s2112Plus   , 1.0/sumOfWeights());
588        scale(_histMeanMultiJPsi1S            , 1.0/sumOfWeights());
589        scale(_histMeanMultiF2_1270           , 1.0/sumOfWeights());
590        scale(_histMeanMultiP                 , 1.0/sumOfWeights());
591        scale(_histMeanMultiLambda            , 1.0/sumOfWeights());
592        scale(_histMeanMultiSigma0            , 1.0/sumOfWeights());
593        scale(_histMeanMultiXiMinus           , 1.0/sumOfWeights());
594        scale(_histMeanMultiDelta1232PlusPlus , 1.0/sumOfWeights());
595        scale(_histMeanMultiSigma1385Minus    , 1.0/sumOfWeights());
596        scale(_histMeanMultiSigma1385Plus     , 1.0/sumOfWeights());
597        scale(_histMeanMultiSigma1385PlusMinus, 1.0/sumOfWeights());
598        scale(_histMeanMultiXi1530_0          , 1.0/sumOfWeights());
599        scale(_histMeanMultiOmegaMinus        , 1.0/sumOfWeights());
600        scale(_histMeanMultiLambda_c_Plus     , 1.0/sumOfWeights());
601        scale(_histMeanMultiSigma_c_PlusPlus_0, 1.0/sumOfWeights());
602        scale(_histMeanMultiLambda1520        , 1.0/sumOfWeights());
603      }
604
605      if (sqrtS()/GeV >= 29 && sqrtS()/GeV <= 35) {
606        scale(_histMeanMultiPiPlus            , 5.0/sumOfWeights());
607        scale(_histMeanMultiPi0               , 5.0/sumOfWeights());
608        scale(_histMeanMultiKPlus             , 5.0/sumOfWeights());
609        scale(_histMeanMultiK0                , 5.0/sumOfWeights());
610        scale(_histMeanMultiEta               , 5.0/sumOfWeights());
611        scale(_histMeanMultiEtaPrime          , 5.0/sumOfWeights());
612        scale(_histMeanMultiDPlus             , 5.0/sumOfWeights());
613        scale(_histMeanMultiD0                , 5.0/sumOfWeights());
614        scale(_histMeanMultiDPlus_s           , 5.0/sumOfWeights());
615        scale(_histMeanMultiF0_980            , 5.0/sumOfWeights());
616        scale(_histMeanMultiRho770_0          , 5.0/sumOfWeights());
617        scale(_histMeanMultiKStar892Plus      , 5.0/sumOfWeights());
618        scale(_histMeanMultiKStar892_0        , 5.0/sumOfWeights());
619        scale(_histMeanMultiPhi1020           , 5.0/sumOfWeights());
620        scale(_histMeanMultiDStar2010Plus     , 5.0/sumOfWeights());
621        scale(_histMeanMultiDStar2007_0       , 5.0/sumOfWeights());
622        scale(_histMeanMultiF2_1270           , 5.0/sumOfWeights());
623        scale(_histMeanMultiK2Star1430Plus    , 5.0/sumOfWeights());
624        scale(_histMeanMultiK2Star1430_0      , 5.0/sumOfWeights());
625        scale(_histMeanMultiP                 , 5.0/sumOfWeights());
626        scale(_histMeanMultiLambda            , 5.0/sumOfWeights());
627        scale(_histMeanMultiXiMinus           , 5.0/sumOfWeights());
628        scale(_histMeanMultiSigma1385Minus    , 5.0/sumOfWeights());
629        scale(_histMeanMultiSigma1385Plus     , 5.0/sumOfWeights());
630        scale(_histMeanMultiSigma1385PlusMinus, 5.0/sumOfWeights());
631        scale(_histMeanMultiOmegaMinus        , 5.0/sumOfWeights());
632        scale(_histMeanMultiLambda_c_Plus     , 5.0/sumOfWeights());
633      }
634
635      if (sqrtS()/GeV >= 89.5 && sqrtS()/GeV <= 91.8) {
636        scale(_histMeanMultiPiPlus            , 1.0/sumOfWeights());
637        scale(_histMeanMultiPi0               , 1.0/sumOfWeights());
638        scale(_histMeanMultiKPlus             , 1.0/sumOfWeights());
639        scale(_histMeanMultiK0                , 1.0/sumOfWeights());
640        scale(_histMeanMultiEta               , 1.0/sumOfWeights());
641        scale(_histMeanMultiEtaPrime          , 1.0/sumOfWeights());
642        scale(_histMeanMultiDPlus             , 1.0/sumOfWeights());
643        scale(_histMeanMultiD0                , 1.0/sumOfWeights());
644        scale(_histMeanMultiDPlus_s           , 1.0/sumOfWeights());
645        scale(_histMeanMultiBPlus_B0_d        , 1.0/sumOfWeights());
646        scale(_histMeanMultiBPlus_u           , 1.0/sumOfWeights());
647        scale(_histMeanMultiB0_s              , 1.0/sumOfWeights());
648        scale(_histMeanMultiF0_980            , 1.0/sumOfWeights());
649        scale(_histMeanMultiA0_980Plus        , 1.0/sumOfWeights());
650        scale(_histMeanMultiRho770_0          , 1.0/sumOfWeights());
651        scale(_histMeanMultiRho770Plus        , 1.0/sumOfWeights());
652        scale(_histMeanMultiOmega782          , 1.0/sumOfWeights());
653        scale(_histMeanMultiKStar892Plus      , 1.0/sumOfWeights());
654        scale(_histMeanMultiKStar892_0        , 1.0/sumOfWeights());
655        scale(_histMeanMultiPhi1020           , 1.0/sumOfWeights());
656        scale(_histMeanMultiDStar2010Plus     , 1.0/sumOfWeights());
657        scale(_histMeanMultiDStar_s2112Plus   , 1.0/sumOfWeights());
658        scale(_histMeanMultiBStar             , 1.0/sumOfWeights());
659        scale(_histMeanMultiJPsi1S            , 1.0/sumOfWeights());
660        scale(_histMeanMultiPsi2S             , 1.0/sumOfWeights());
661        scale(_histMeanMultiUpsilon1S         , 1.0/sumOfWeights());
662        scale(_histMeanMultiF1_1285           , 1.0/sumOfWeights());
663        scale(_histMeanMultiF1_1420           , 1.0/sumOfWeights());
664        scale(_histMeanMultiChi_c1_3510       , 1.0/sumOfWeights());
665        scale(_histMeanMultiF2_1270           , 1.0/sumOfWeights());
666        scale(_histMeanMultiF2Prime1525       , 1.0/sumOfWeights());
667        scale(_histMeanMultiK2Star1430_0      , 1.0/sumOfWeights());
668        scale(_histMeanMultiBStarStar         , 1.0/sumOfWeights());
669        scale(_histMeanMultiDs1Plus           , 1.0/sumOfWeights());
670        scale(_histMeanMultiDs2Plus           , 1.0/sumOfWeights());
671        scale(_histMeanMultiP                 , 1.0/sumOfWeights());
672        scale(_histMeanMultiLambda            , 1.0/sumOfWeights());
673        scale(_histMeanMultiSigma0            , 1.0/sumOfWeights());
674        scale(_histMeanMultiSigmaMinus        , 1.0/sumOfWeights());
675        scale(_histMeanMultiSigmaPlus         , 1.0/sumOfWeights());
676        scale(_histMeanMultiSigmaPlusMinus    , 1.0/sumOfWeights());
677        scale(_histMeanMultiXiMinus           , 1.0/sumOfWeights());
678        scale(_histMeanMultiDelta1232PlusPlus , 1.0/sumOfWeights());
679        scale(_histMeanMultiSigma1385Minus    , 1.0/sumOfWeights());
680        scale(_histMeanMultiSigma1385Plus     , 1.0/sumOfWeights());
681        scale(_histMeanMultiSigma1385PlusMinus, 1.0/sumOfWeights());
682        scale(_histMeanMultiXi1530_0          , 1.0/sumOfWeights());
683        scale(_histMeanMultiOmegaMinus        , 1.0/sumOfWeights());
684        scale(_histMeanMultiLambda_c_Plus     , 1.0/sumOfWeights());
685        scale(_histMeanMultiLambda_b_0        , 1.0/sumOfWeights());
686        scale(_histMeanMultiLambda1520        , 1.0/sumOfWeights());
687      }
688
689      if (sqrtS()/GeV >= 130 && sqrtS()/GeV <= 200) {
690        scale(_histMeanMultiPiPlus           , 70.0/sumOfWeights());
691        scale(_histMeanMultiKPlus            , 70.0/sumOfWeights());
692        scale(_histMeanMultiK0               , 70.0/sumOfWeights());
693        scale(_histMeanMultiP                , 70.0/sumOfWeights());
694        scale(_histMeanMultiLambda           , 70.0/sumOfWeights());
695      }
696    }
697
698    /// @}
699
700
701  private:
702
703    Histo1DPtr _histMeanMultiPiPlus;
704    Histo1DPtr _histMeanMultiPi0;
705    Histo1DPtr _histMeanMultiKPlus;
706    Histo1DPtr _histMeanMultiK0;
707    Histo1DPtr _histMeanMultiEta;
708    Histo1DPtr _histMeanMultiEtaPrime;
709    Histo1DPtr _histMeanMultiDPlus;
710    Histo1DPtr _histMeanMultiD0;
711    Histo1DPtr _histMeanMultiDPlus_s;
712    Histo1DPtr _histMeanMultiBPlus_B0_d;
713    Histo1DPtr _histMeanMultiBPlus_u;
714    Histo1DPtr _histMeanMultiB0_s;
715    Histo1DPtr _histMeanMultiF0_980;
716    Histo1DPtr _histMeanMultiA0_980Plus;
717    Histo1DPtr _histMeanMultiRho770_0;
718    Histo1DPtr _histMeanMultiRho770Plus;
719    Histo1DPtr _histMeanMultiOmega782;
720    Histo1DPtr _histMeanMultiKStar892Plus;
721    Histo1DPtr _histMeanMultiKStar892_0;
722    Histo1DPtr _histMeanMultiPhi1020;
723    Histo1DPtr _histMeanMultiDStar2010Plus;
724    Histo1DPtr _histMeanMultiDStar2007_0;
725    Histo1DPtr _histMeanMultiDStar_s2112Plus;
726    Histo1DPtr _histMeanMultiBStar;
727    Histo1DPtr _histMeanMultiJPsi1S;
728    Histo1DPtr _histMeanMultiPsi2S;
729    Histo1DPtr _histMeanMultiUpsilon1S;
730    Histo1DPtr _histMeanMultiF1_1285;
731    Histo1DPtr _histMeanMultiF1_1420;
732    Histo1DPtr _histMeanMultiChi_c1_3510;
733    Histo1DPtr _histMeanMultiF2_1270;
734    Histo1DPtr _histMeanMultiF2Prime1525;
735    Histo1DPtr _histMeanMultiK2Star1430Plus;
736    Histo1DPtr _histMeanMultiK2Star1430_0;
737    Histo1DPtr _histMeanMultiBStarStar;
738    Histo1DPtr _histMeanMultiDs1Plus;
739    Histo1DPtr _histMeanMultiDs2Plus;
740    Histo1DPtr _histMeanMultiP;
741    Histo1DPtr _histMeanMultiLambda;
742    Histo1DPtr _histMeanMultiSigma0;
743    Histo1DPtr _histMeanMultiSigmaMinus;
744    Histo1DPtr _histMeanMultiSigmaPlus;
745    Histo1DPtr _histMeanMultiSigmaPlusMinus;
746    Histo1DPtr _histMeanMultiXiMinus;
747    Histo1DPtr _histMeanMultiDelta1232PlusPlus;
748    Histo1DPtr _histMeanMultiSigma1385Minus;
749    Histo1DPtr _histMeanMultiSigma1385Plus;
750    Histo1DPtr _histMeanMultiSigma1385PlusMinus;
751    Histo1DPtr _histMeanMultiXi1530_0;
752    Histo1DPtr _histMeanMultiOmegaMinus;
753    Histo1DPtr _histMeanMultiLambda_c_Plus;
754    Histo1DPtr _histMeanMultiLambda_b_0;
755    Histo1DPtr _histMeanMultiSigma_c_PlusPlus_0;
756    Histo1DPtr _histMeanMultiLambda1520;
757
758    /// @}
759
760  };
761
762
763
764  RIVET_DECLARE_PLUGIN(PDG_HADRON_MULTIPLICITIES);
765
766}