1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
18 // Implementation of AliDecayer using Pythia8
19 // Author: andreas.morsch@cern.ch
23 #include "AliDecayerPythia8.h"
24 #include "ParticleData.h"
26 ClassImp(AliDecayerPythia8)
28 Bool_t AliDecayerPythia8::fgInit = kFALSE;
30 AliDecayerPythia8::AliDecayerPythia8():
38 void AliDecayerPythia8::Init()
44 // fPythia->SetDecayTable();
47 // Switch on heavy flavor decays
50 Int_t heavy[14] = {411, 421, 431, 4122, 4132, 4232, 4332, 511, 521, 531, 5122, 5132, 5232, 5332};
51 // fPythia->ResetDecayTable();
52 for (j=0; j < 14; j++) {
53 if (fDecay == kNoDecayHeavy) {
54 TPythia8::Instance()->ReadString(Form("%d:onMode = off", heavy[j]));
56 TPythia8::Instance()->ReadString(Form("%d:onMode = on", heavy[j]));
61 //...Switch off decay of pi0, K0S, Lambda, Sigma+-, Xi0-, Omega-.
63 if (fDecay != kNeutralPion) {
64 TPythia8::Instance()->ReadString("111:onMode = off");
66 TPythia8::Instance()->ReadString("111:onMode = on");
69 TPythia8::Instance()->ReadString("310:onMode = off");
70 TPythia8::Instance()->ReadString("3112:onMode = off");
71 TPythia8::Instance()->ReadString("3212:onMode = off");
72 TPythia8::Instance()->ReadString("3222:onMode = off");
73 TPythia8::Instance()->ReadString("3312:onMode = off");
74 TPythia8::Instance()->ReadString("3322:onMode = off");
75 TPythia8::Instance()->ReadString("3334:onMode = off");
76 // .. Force decay channels
80 void AliDecayerPythia8::ForceDecay()
83 // Force a particle decay mode
84 // Switch heavy flavour production off if requested
85 if (!fHeavyFlavour) SwitchOffHeavyFlavour();
87 Decay_t decay = fDecay;
88 TPythia8::Instance()->ReadString("HadronLevel:Decay = on");
90 if (decay == kNoDecayHeavy) return;
98 TPythia8::Instance()->ReadString("511:onMode = off");
99 TPythia8::Instance()->ReadString("511:onIfAny = 13 443 100443");
101 TPythia8::Instance()->ReadString("521:onMode = off");
102 TPythia8::Instance()->ReadString("521:onIfAny = 13 443 100443");
104 TPythia8::Instance()->ReadString("531:onMode = off");
105 TPythia8::Instance()->ReadString("531:onIfAny = 13 443 100443");
107 TPythia8::Instance()->ReadString("5122:onMode = off");
108 TPythia8::Instance()->ReadString("5122:onIfAny = 13 443 100443");
110 TPythia8::Instance()->ReadString("5132:onMode = off");
111 TPythia8::Instance()->ReadString("5132:onIfAny = 13 443 100443");
113 TPythia8::Instance()->ReadString("5232:onMode = off");
114 TPythia8::Instance()->ReadString("5232:onIfAny = 13 443 100443");
116 TPythia8::Instance()->ReadString("5332:onMode = off");
117 TPythia8::Instance()->ReadString("5332:onIfAny = 13 443 100443");
119 TPythia8::Instance()->ReadString("100443:onMode = off");
120 TPythia8::Instance()->ReadString("100443:onIfAny = 443");
122 TPythia8::Instance()->ReadString("443:onMode = off");
123 TPythia8::Instance()->ReadString("443:onIfAll = 13 13");
125 TPythia8::Instance()->ReadString("411:onMode = off");
126 TPythia8::Instance()->ReadString("411:onIfAll = 13");
128 TPythia8::Instance()->ReadString("421:onMode = off");
129 TPythia8::Instance()->ReadString("421:onIfAll = 13");
131 TPythia8::Instance()->ReadString("431:onMode = off");
132 TPythia8::Instance()->ReadString("431:onIfAll = 13");
134 TPythia8::Instance()->ReadString("4122:onMode = off");
135 TPythia8::Instance()->ReadString("4122:onIfAll = 13");
137 TPythia8::Instance()->ReadString("4132:onMode = off");
138 TPythia8::Instance()->ReadString("4132:onIfAll = 13");
140 TPythia8::Instance()->ReadString("4232:onMode = off");
141 TPythia8::Instance()->ReadString("4232:onIfAll = 13");
143 TPythia8::Instance()->ReadString("4332:onMode = off");
144 TPythia8::Instance()->ReadString("4332:onIfAll = 13");
147 case kChiToJpsiGammaToMuonMuon:
148 // Chi_1c -> J/Psi Gamma
149 TPythia8::Instance()->ReadString("20443:onMode = off");
150 TPythia8::Instance()->ReadString("20443:onIfAll = 443 22");
151 // Chi_2c -> J/Psi Gamma
152 TPythia8::Instance()->ReadString("445:onMode = off");
153 TPythia8::Instance()->ReadString("445:onIfAll = 443 22");
155 TPythia8::Instance()->ReadString("443:onMode = off");
156 TPythia8::Instance()->ReadString("443:onIfAll = 13 13");
158 case kChiToJpsiGammaToElectronElectron:
159 // Chi_1c -> J/Psi Gamma
160 TPythia8::Instance()->ReadString("20443:onMode = off");
161 TPythia8::Instance()->ReadString("20443:onIfAll = 443 22");
162 // Chi_2c -> J/Psi Gamma
163 TPythia8::Instance()->ReadString("445:onMode = off");
164 TPythia8::Instance()->ReadString("445:onIfAll = 443 22");
166 TPythia8::Instance()->ReadString("443:onMode = off");
167 TPythia8::Instance()->ReadString("443:onIfAll = 11 11");
172 TPythia8::Instance()->ReadString("511:onMode = off");
173 TPythia8::Instance()->ReadString("511:onIfAny = 13");
175 TPythia8::Instance()->ReadString("521:onMode = off");
176 TPythia8::Instance()->ReadString("521:onIfAny = 13");
178 TPythia8::Instance()->ReadString("531:onMode = off");
179 TPythia8::Instance()->ReadString("531:onIfAny = 13");
181 TPythia8::Instance()->ReadString("5122:onMode = off");
182 TPythia8::Instance()->ReadString("5122:onIfAny = 13");
184 TPythia8::Instance()->ReadString("5132:onMode = off");
185 TPythia8::Instance()->ReadString("5132:onIfAny = 13");
187 TPythia8::Instance()->ReadString("5232:onMode = off");
188 TPythia8::Instance()->ReadString("5232:onIfAny = 13");
190 TPythia8::Instance()->ReadString("5332:onMode = off");
191 TPythia8::Instance()->ReadString("5332:onIfAny = 13");
195 TPythia8::Instance()->ReadString("411:onMode = off");
196 TPythia8::Instance()->ReadString("411:onIfAll = 13");
198 TPythia8::Instance()->ReadString("421:onMode = off");
199 TPythia8::Instance()->ReadString("421:onIfAll = 13");
201 TPythia8::Instance()->ReadString("431:onMode = off");
202 TPythia8::Instance()->ReadString("431:onIfAll = 13");
204 TPythia8::Instance()->ReadString("4122:onMode = off");
205 TPythia8::Instance()->ReadString("4122:onIfAll = 13");
207 TPythia8::Instance()->ReadString("4132:onMode = off");
208 TPythia8::Instance()->ReadString("4132:onIfAll = 13");
210 TPythia8::Instance()->ReadString("4232:onMode = off");
211 TPythia8::Instance()->ReadString("4232:onIfAll = 13");
213 TPythia8::Instance()->ReadString("4332:onMode = off");
214 TPythia8::Instance()->ReadString("4332:onIfAll = 13");
216 TPythia8::Instance()->ReadString("511:onMode = off");
217 TPythia8::Instance()->ReadString("511:onIfAny = 13");
219 TPythia8::Instance()->ReadString("521:onMode = off");
220 TPythia8::Instance()->ReadString("521:onIfAny = 13");
222 TPythia8::Instance()->ReadString("531:onMode = off");
223 TPythia8::Instance()->ReadString("531:onIfAny = 13");
225 TPythia8::Instance()->ReadString("5122:onMode = off");
226 TPythia8::Instance()->ReadString("5122:onIfAny = 13");
228 TPythia8::Instance()->ReadString("5132:onMode = off");
229 TPythia8::Instance()->ReadString("5132:onIfAny = 13");
231 TPythia8::Instance()->ReadString("5232:onMode = off");
232 TPythia8::Instance()->ReadString("5232:onIfAny = 13");
234 TPythia8::Instance()->ReadString("5332:onMode = off");
235 TPythia8::Instance()->ReadString("5332:onIfAny = 13");
240 TPythia8::Instance()->ReadString("113:onMode = off");
241 TPythia8::Instance()->ReadString("113:onIfAll = 13 13");
243 TPythia8::Instance()->ReadString("221:onMode = off");
244 TPythia8::Instance()->ReadString("221:onIfAll = 13 13");
246 TPythia8::Instance()->ReadString("223:onMode = off");
247 TPythia8::Instance()->ReadString("223:onIfAll = 13 13");
249 TPythia8::Instance()->ReadString("333:onMode = off");
250 TPythia8::Instance()->ReadString("333:onIfAll = 13 13");
252 TPythia8::Instance()->ReadString("443:onMode = off");
253 TPythia8::Instance()->ReadString("443:onIfAll = 13 13");
255 TPythia8::Instance()->ReadString("100443:onMode = off");
256 TPythia8::Instance()->ReadString("100443:onIfAll = 13 13");
258 TPythia8::Instance()->ReadString("553:onMode = off");
259 TPythia8::Instance()->ReadString("553:onIfAll = 13 13");
261 TPythia8::Instance()->ReadString("100553:onMode = off");
262 TPythia8::Instance()->ReadString("100553:onIfAll = 13 13");
264 TPythia8::Instance()->ReadString("200553:onMode = off");
265 TPythia8::Instance()->ReadString("200553:onIfAll = 13 13");
267 case kBSemiElectronic:
269 TPythia8::Instance()->ReadString("511:onMode = off");
270 TPythia8::Instance()->ReadString("511:onIfAny = 11");
272 TPythia8::Instance()->ReadString("521:onMode = off");
273 TPythia8::Instance()->ReadString("521:onIfAny = 11");
275 TPythia8::Instance()->ReadString("531:onMode = off");
276 TPythia8::Instance()->ReadString("531:onIfAny = 11");
278 TPythia8::Instance()->ReadString("5122:onMode = off");
279 TPythia8::Instance()->ReadString("5122:onIfAny = 11");
281 TPythia8::Instance()->ReadString("5132:onMode = off");
282 TPythia8::Instance()->ReadString("5132:onIfAny = 11");
284 TPythia8::Instance()->ReadString("5232:onMode = off");
285 TPythia8::Instance()->ReadString("5232:onIfAny = 11");
287 TPythia8::Instance()->ReadString("5332:onMode = off");
288 TPythia8::Instance()->ReadString("5332:onIfAny = 11");
290 case kSemiElectronic:
292 TPythia8::Instance()->ReadString("411:onMode = off");
293 TPythia8::Instance()->ReadString("411:onIfAll = 11");
295 TPythia8::Instance()->ReadString("421:onMode = off");
296 TPythia8::Instance()->ReadString("421:onIfAll = 11");
298 TPythia8::Instance()->ReadString("431:onMode = off");
299 TPythia8::Instance()->ReadString("431:onIfAll = 11");
301 TPythia8::Instance()->ReadString("4122:onMode = off");
302 TPythia8::Instance()->ReadString("4122:onIfAll = 11");
304 TPythia8::Instance()->ReadString("4132:onMode = off");
305 TPythia8::Instance()->ReadString("4132:onIfAll = 11");
307 TPythia8::Instance()->ReadString("4232:onMode = off");
308 TPythia8::Instance()->ReadString("4232:onIfAll = 11");
310 TPythia8::Instance()->ReadString("4332:onMode = off");
311 TPythia8::Instance()->ReadString("4332:onIfAll = 11");
313 TPythia8::Instance()->ReadString("511:onMode = off");
314 TPythia8::Instance()->ReadString("511:onIfAny = 11");
316 TPythia8::Instance()->ReadString("521:onMode = off");
317 TPythia8::Instance()->ReadString("521:onIfAny = 11");
319 TPythia8::Instance()->ReadString("531:onMode = off");
320 TPythia8::Instance()->ReadString("531:onIfAny = 11");
322 TPythia8::Instance()->ReadString("5122:onMode = off");
323 TPythia8::Instance()->ReadString("5122:onIfAny = 11");
325 TPythia8::Instance()->ReadString("5132:onMode = off");
326 TPythia8::Instance()->ReadString("5132:onIfAny = 11");
328 TPythia8::Instance()->ReadString("5232:onMode = off");
329 TPythia8::Instance()->ReadString("5232:onIfAny = 11");
331 TPythia8::Instance()->ReadString("5332:onMode = off");
332 TPythia8::Instance()->ReadString("5332:onIfAny = 11");
336 TPythia8::Instance()->ReadString("113:onMode = off");
337 TPythia8::Instance()->ReadString("113:onIfAll = 11 11");
339 TPythia8::Instance()->ReadString("221:onMode = off");
340 TPythia8::Instance()->ReadString("221:onIfAll = 11 11");
342 TPythia8::Instance()->ReadString("223:onMode = off");
343 TPythia8::Instance()->ReadString("223:onIfAll = 11 11");
345 TPythia8::Instance()->ReadString("333:onMode = off");
346 TPythia8::Instance()->ReadString("333:onIfAll = 11 11");
348 TPythia8::Instance()->ReadString("443:onMode = off");
349 TPythia8::Instance()->ReadString("443:onIfAll = 11 11");
351 TPythia8::Instance()->ReadString("100443:onMode = off");
352 TPythia8::Instance()->ReadString("100443:onIfAll = 11 11");
354 TPythia8::Instance()->ReadString("553:onMode = off");
355 TPythia8::Instance()->ReadString("553:onIfAll = 11 11");
357 TPythia8::Instance()->ReadString("100553:onMode = off");
358 TPythia8::Instance()->ReadString("100553:onIfAll = 11 11");
360 TPythia8::Instance()->ReadString("200553:onMode = off");
361 TPythia8::Instance()->ReadString("200553:onIfAll = 11 11");
364 // B0 -> J/Psi (Psi') X
365 TPythia8::Instance()->ReadString("511:onMode = off");
366 TPythia8::Instance()->ReadString("511:onIfAny = 443 100443");
367 // B+/- -> J/Psi (Psi') X
368 TPythia8::Instance()->ReadString("521:onMode = off");
369 TPythia8::Instance()->ReadString("521:onIfAny = 443 100443");
370 // B_s -> J/Psi (Psi') X
371 TPythia8::Instance()->ReadString("531:onMode = off");
372 TPythia8::Instance()->ReadString("531:onIfAny = 443 100443");
373 // Lambda_b -> J/Psi (Psi') X
374 TPythia8::Instance()->ReadString("5122:onMode = off");
375 TPythia8::Instance()->ReadString("5122:onIfAny = 443 100443");
378 TPythia8::Instance()->ReadString("443:onMode = off");
379 TPythia8::Instance()->ReadString("443:onIfAll = 13 13");
381 TPythia8::Instance()->ReadString("100443:onMode = off");
382 TPythia8::Instance()->ReadString("100443:onIfAll = 13 13");
384 case kBPsiPrimeDiMuon:
386 TPythia8::Instance()->ReadString("511:onMode = off");
387 TPythia8::Instance()->ReadString("511:onIfAny = 100443");
389 TPythia8::Instance()->ReadString("521:onMode = off");
390 TPythia8::Instance()->ReadString("521:onIfAny = 100443");
392 TPythia8::Instance()->ReadString("531:onMode = off");
393 TPythia8::Instance()->ReadString("531:onIfAny = 100443");
394 // Lambda_b -> Psi' X
395 TPythia8::Instance()->ReadString("5122:onMode = off");
396 TPythia8::Instance()->ReadString("5122:onIfAny = 100443");
399 TPythia8::Instance()->ReadString("100443:onMode = off");
400 TPythia8::Instance()->ReadString("100443:onIfAll = 13 13");
402 case kBJpsiDiElectron:
404 TPythia8::Instance()->ReadString("511:onMode = off");
405 TPythia8::Instance()->ReadString("511:onIfAny = 443");
407 TPythia8::Instance()->ReadString("521:onMode = off");
408 TPythia8::Instance()->ReadString("521:onIfAny = 443");
410 TPythia8::Instance()->ReadString("531:onMode = off");
411 TPythia8::Instance()->ReadString("531:onIfAny = 443");
413 TPythia8::Instance()->ReadString("5122:onMode = off");
414 TPythia8::Instance()->ReadString("5122:onIfAny = 443");
417 TPythia8::Instance()->ReadString("443:onMode = off");
418 TPythia8::Instance()->ReadString("443:onIfAll = 11 11");
423 TPythia8::Instance()->ReadString("511:onMode = off");
424 TPythia8::Instance()->ReadString("511:onIfAny = 443");
426 TPythia8::Instance()->ReadString("521:onMode = off");
427 TPythia8::Instance()->ReadString("521:onIfAny = 443");
429 TPythia8::Instance()->ReadString("531:onMode = off");
430 TPythia8::Instance()->ReadString("531:onIfAny = 443");
432 TPythia8::Instance()->ReadString("5122:onMode = off");
433 TPythia8::Instance()->ReadString("5122:onIfAny = 443");
435 case kBPsiPrimeDiElectron:
437 TPythia8::Instance()->ReadString("511:onMode = off");
438 TPythia8::Instance()->ReadString("511:onIfAny = 100443");
440 TPythia8::Instance()->ReadString("521:onMode = off");
441 TPythia8::Instance()->ReadString("521:onIfAny = 100443");
443 TPythia8::Instance()->ReadString("531:onMode = off");
444 TPythia8::Instance()->ReadString("531:onIfAny = 100443");
445 // Lambda_b -> Psi' X
446 TPythia8::Instance()->ReadString("5122:onMode = off");
447 TPythia8::Instance()->ReadString("5122:onIfAny = 100443");
450 TPythia8::Instance()->ReadString("100443:onMode = off");
451 TPythia8::Instance()->ReadString("100443:onIfAll = 11 11");
455 TPythia8::Instance()->ReadString("211:onMode = off");
456 TPythia8::Instance()->ReadString("211:onIfAny = 13");
460 TPythia8::Instance()->ReadString("321:onMode = off");
461 TPythia8::Instance()->ReadString("321:onIfAny = 13");
465 TPythia8::Instance()->ReadString("211:onMode = off");
466 TPythia8::Instance()->ReadString("211:onIfAny = 13");
467 TPythia8::Instance()->ReadString("321:onMode = off");
468 TPythia8::Instance()->ReadString("321:onIfAny = 13");
472 TPythia8::Instance()->ReadString("24:onMode = off");
473 TPythia8::Instance()->ReadString("24:onIfAny = 13");
477 TPythia8::Instance()->ReadString("24:onMode = off");
478 TPythia8::Instance()->ReadString("24:onIfAny = 4");
480 case kWToCharmToMuon:
482 TPythia8::Instance()->ReadString("24:onMode = off");
483 TPythia8::Instance()->ReadString("24:onIfAny = 4");
485 TPythia8::Instance()->ReadString("411:onMode = off");
486 TPythia8::Instance()->ReadString("411:onIfAll = 13");
488 TPythia8::Instance()->ReadString("421:onMode = off");
489 TPythia8::Instance()->ReadString("421:onIfAll = 13");
491 TPythia8::Instance()->ReadString("431:onMode = off");
492 TPythia8::Instance()->ReadString("431:onIfAll = 13");
494 TPythia8::Instance()->ReadString("4122:onMode = off");
495 TPythia8::Instance()->ReadString("4122:onIfAll = 13");
497 TPythia8::Instance()->ReadString("4132:onMode = off");
498 TPythia8::Instance()->ReadString("4132:onIfAll = 13");
500 TPythia8::Instance()->ReadString("4232:onMode = off");
501 TPythia8::Instance()->ReadString("4232:onIfAll = 13");
503 TPythia8::Instance()->ReadString("4332:onMode = off");
504 TPythia8::Instance()->ReadString("4332:onIfAll = 13");
508 TPythia8::Instance()->ReadString("23:onMode = off");
509 TPythia8::Instance()->ReadString("23:onIfAll = 13 13");
513 TPythia8::Instance()->ReadString("23:onMode = off");
514 TPythia8::Instance()->ReadString("23:onIfAll = 11 11");
519 case kHadronicDWithout4Bodies:
524 TPythia8::Instance()->ReadString("333:onMode = off");
525 TPythia8::Instance()->ReadString("333:onIfAll = 321 321");
529 TPythia8::Instance()->ReadString("3334:onMode = off");
530 TPythia8::Instance()->ReadString("3334:onIfAll = 3122 321 ");
533 TPythia8::Instance()->ReadString("3122:onMode = off");
534 TPythia8::Instance()->ReadString("3122:onIfAll = 2212 211 ");
538 TPythia8::Instance()->ReadString("HadronLevel:Decay = off");
547 Float_t AliDecayerPythia8::GetPartialBranchingRatio(Int_t ipart)
549 // Get the partial branching ration for the forced decay channels
551 Pythia8::Pythia* thePythia = TPythia8::Instance()->Pythia8();
552 Pythia8::ParticleDataTable table = thePythia->particleData;
553 Pythia8::ParticleDataEntry* pd = table.particleDataPtr(ipart);
554 Pythia8::DecayTable decays = pd->decay;
557 Int_t nc = decays.size();
560 // Loop over decay channels
561 for (Int_t ic = 0; ic < nc; ic++) {
562 Pythia8::DecayChannel& decCh = decays[ic];
563 for (Int_t i = 0; i < decCh.multiplicity(); i++) {
564 br += decCh.bRatio();
571 Float_t AliDecayerPythia8::GetLifetime(Int_t kf)
573 // Return lifetime of particle
574 Pythia8::Pythia* thePythia = TPythia8::Instance()->Pythia8();
575 Pythia8::ParticleDataTable table = thePythia->particleData;
576 Float_t tau = table.tau0(kf);
580 void AliDecayerPythia8::SwitchOffHeavyFlavour()
582 // Switch off heavy flavour production
584 // Maximum number of quark flavours used in pdf
585 TPythia8::Instance()->ReadString("PDFinProcess:nQuarkIn = 3");
586 // Maximum number of flavors that can be used in showers
587 TPythia8::Instance()->ReadString("SpaceShower:nQuarkIn = 3");
588 TPythia8::Instance()->ReadString("TimeShower:nGammaToQuark = 3");
589 TPythia8::Instance()->ReadString("TimeShower:nGluonToQuark = 3");
593 void AliDecayerPythia8::ForceHadronicD(Int_t optUse4Bodies)
596 // Force golden D decay modes
599 TPythia8::Instance()->ReadString("313:onMode = off");
600 TPythia8::Instance()->ReadString("313:onIfAll = 321 211");
602 TPythia8::Instance()->ReadString("333:onMode = off");
603 TPythia8::Instance()->ReadString("333:onIfAll = 321 321");
604 // for D0 -> rho0 pi+ k-
605 TPythia8::Instance()->ReadString("113:onMode = off");
606 TPythia8::Instance()->ReadString("113:onIfAll = 211 211");
607 // for Lambda_c -> Delta++ K-
608 TPythia8::Instance()->ReadString("2224:onMode = off");
609 TPythia8::Instance()->ReadString("2224:onIfAll = 2212 211");
610 // for Lambda_c -> Lambda(1520) K-
611 TPythia8::Instance()->ReadString("3124:onMode = off");
612 TPythia8::Instance()->ReadString("3124:onIfAll = 2212 321");
615 TPythia8::Instance()->ReadString("411:onMode = off");
616 TPythia8::Instance()->ReadString("421:onMode = off");
617 TPythia8::Instance()->ReadString("431:onMode = off");
618 TPythia8::Instance()->ReadString("4112:onMode = off");
619 TPythia8::Instance()->ReadString("4122:onMode = off");
622 TPythia8::Instance()->ReadString("411:onIfMatch = 321 211 211");
624 TPythia8::Instance()->ReadString("411:onIfMatch = 313 211");
626 TPythia8::Instance()->ReadString("421:onIfMatch = 321 211");
630 TPythia8::Instance()->ReadString("421:onIfMatch = 321 211 211 211");
632 TPythia8::Instance()->ReadString("421:onIfMatch = 321 211 113");
634 TPythia8::Instance()->ReadString("421:onIfMatch = 313 211 211");
638 TPythia8::Instance()->ReadString("431:onIfMatch = 321 313");
640 TPythia8::Instance()->ReadString("431:onIfMatch = 333 211");
643 TPythia8::Instance()->ReadString("4122:onIfMatch = 2212 313");
644 // Lambda_c -> Delta K
645 TPythia8::Instance()->ReadString("4122:onIfMatch = 2224 321");
646 // Lambda_c -> Lambda(1520) pi
647 TPythia8::Instance()->ReadString("4122:onIfMatch = 3124 211");
648 // Lambda_c -> p K pi
649 TPythia8::Instance()->ReadString("4122:onIfMatch = 2212 321 211");
650 // Lambda_c -> Lambda pi
651 TPythia8::Instance()->ReadString("4122:onIfMatch = 3122 211");