Rivet analyses referencePDG_HADRON_MULTIPLICITIESHadron multiplicities in hadronic $e^+e^-$ eventsExperiment: PDG (Various) Spires ID: 7857373 Status: VALIDATED Authors:
Beam energies: (5.0, 5.0); (17.5, 17.5); (45.6, 45.6); (88.5, 88.5) GeV Run details:
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}
|