O2DPG/MC/config/PWGHF/pythia8/generator/pythia8_lambdab_probQQtoQ.cfg at 3d977825240dc2d7b8715e04dde2deec5de2cae2 · ercolessi/O2DPG · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
### beams
Beams:idA 2212			# proton
Beams:idB 2212 			# proton
Beams:eCM 13600. 		# GeV

### processes
SoftQCD:inelastic on		# all inelastic processes

### decays
ParticleDecays:limitTau0 on
ParticleDecays:tau0Max 10.

# enable deuteron production by coalescence collisions
HadronLevel:DeuteronProduction = on
DeuteronProduction:channels = {2212 2112 > 22}
DeuteronProduction:models = {0}
DeuteronProduction:norm = 1
DeuteronProduction:parms = {0.5 1} # coalescence momentum p0 in GeV/c

### switching on Pythia Mode2
ColourReconnection:mode 1
ColourReconnection:allowDoubleJunRem off
ColourReconnection:m0 0.3
ColourReconnection:allowJunctions on
ColourReconnection:junctionCorrection 1.20
ColourReconnection:timeDilationMode 2
ColourReconnection:timeDilationPar 0.18
StringPT:sigma 0.335
StringZ:aLund 0.36
StringZ:bLund 0.56
#StringFlav:probQQtoQ 0.078
StringFlav:probQQtoQ 0.17
StringFlav:ProbStoUD 0.2
StringFlav:probQQ1toQQ0join 0.0275,0.0275,0.0275,0.0275
MultiPartonInteractions:pT0Ref 2.15
BeamRemnants:remnantMode 1
BeamRemnants:saturation 5

# Correct decay lengths (wrong in PYTHIA8 decay table)
# Lb
5122:tau0 = 0.4390
# Xic0
4132:tau0 = 0.0455
# OmegaC
4332:tau0 = 0.0803

## Lb decay modes
5122:onMode = off
5122:oneChannel = 1 9.03e-05  0 2212 2212 -2212 -2212 2112
5122:addChannel = 1 0.00181  0 2212 2212 -2212 -2212 2112 111
5122:addChannel = 1 0.0181  0 2212 2212 -2212 -2212 2112 111 111
5122:addChannel = 1 0.0271  0 2212 2212 -2212 -2212 2112 211 -211
5122:addChannel = 1 0.0181  0 2212 2212 -2212 -2212 2112 211 211 -211 -211
5122:addChannel = 1 0.181  0 2212 2212 -2212 -2212 2112 211 -211 111
5122:addChannel = 1 0.199  0 2212 2212 -2212 -2212 2112 211 -211 111 111
5122:addChannel = 1 0.0271  0 2212 2212 -2212 -2212 2112 211 211 -211 -211 111
5122:addChannel = 1 0.00452  0 2212 2212 -2212 -2212 2112 211 111
5122:addChannel = 1 0.0361  0 2212 2212 2112 -2212 -2112 -211 111
5122:addChannel = 1 0.0452  0 2212 2212 2112 -2212 -2112 -211 111 111
5122:addChannel = 1 0.0542  0 2212 2212 2112 -2212 -2112 -211 -211 211
5122:addChannel = 1 0.361  0 2212 2212 2112 -2212 -2112 -211 -211 211 111
5122:addChannel = 1 0.0271  0 2212 2212 2112 -2212 -2112 -211 -211 211 111 111
5122:onIfMatch = 2212 2212 2212 2212 2112
5122:onIfMatch = 2212 2212 2212 2212 2112 111
5122:onIfMatch = 2212 2212 2212 2212 2112 111 111
5122:onIfMatch = 2212 2212 2212 2212 2112 211 211
5122:onIfMatch = 2212 2212 2212 2212 2112 211 211 211 211
5122:onIfMatch = 2212 2212 2212 2212 2112 211 211 111
5122:onIfMatch = 2212 2212 2212 2212 2112 211 211 111 111
5122:onIfMatch = 2212 2212 2212 2212 2112 211 211 211 211 111
5122:onIfMatch = 2212 2212 2212 2212 2112 211 111
5122:onIfMatch = 2212 2212 2112 2212 2112 211 111
5122:onIfMatch = 2212 2212 2112 2212 2112 211 111 111
5122:onIfMatch = 2212 2212 2112 2212 2112 211 211 211
5122:onIfMatch = 2212 2212 2112 2212 2112 211 211 211 111
5122:onIfMatch = 2212 2212 2112 2212 2112 211 211 211 111 111