Newer
Older
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
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
Pythia6-EIC
===========
.. code:: sql
%%writefile pythia6.eic
ep_Q2-10-1000.18x275GeV.100ev.txt ! output file name
11 ! lepton beam type
275, 18 ! proton and electron beam energy
10000,1 ! Number of events <nub to generate>, <num to display>
1e-09, 0.99 ! xmin and xmax
1e-09,1.00 ! ymin and ymax
10,1000 ! Q2min and Q2max
F2PY,1998 ! F2-Model, R-Parametrisation
0 ! switch for rad corrections; 0:no, 1:yes, 2:gen.lookup table
1 ! Pythia-Model = 0 standard GVMD generation in Pythia-x and Q2; = 1 GVMD model with generation in y and Q2 as for radgen
1,1 ! A-Tar and Z-Tar
1,1 ! nuclear pdf parameter1: nucleon mass number A, charge number Z
201 ! nuclear pdf parameter2: correction order x*100+y x= 1:LO, 2:NLO y:error set
! PMAS(4,1)=1.27 ! charm mass
MSEL=2
MSTP(14)=30
MSTP(15)=0
MSTP(16)=1
MSTP(17)=4 ! MSTP 17=6 is the R-rho measured as by hermes, =4 Default
MSTP(18)=3
MSTP(19)=1 ! Hermes MSTP-19=1 different Q2 suppression, default = 4
MSTP(20)=0 ! Hermes MSTP(20)=0 , default MSTP(20)=3
MSTP(32)=8
MSTP(38)=4
!MSTP(51)=10150 ! if pdflib is linked than non pythia-pdfs are available, like MSTP(51)=4046
MSTP(51)=7
MSTP(52)=1 ! ---> pdflib used MSTP 52=2
MSTP(53)=3
MSTP(54)=1
MSTP(55)=5
MSTP(56)=1
MSTP(57)=1
MSTP(58)=5
MSTP(59)=1
MSTP(60)=7
MSTP(61)=2
MSTP(71)=1
MSTP(81)=0
MSTP(82)=1
MSTP(91)=1
MSTP(92)=3 ! hermes MSTP(92)=4
MSTP(93)=1
MSTP(101)=3
MSTP(102)=1
MSTP(111)=1
MSTP(121)=0
! ----------- Now all the PARPs -----------
PARP(13)=1
PARP(18)=0.40 ! hermes PARP(18)=0.17
PARP(81)=1.9
PARP(89)=1800
PARP(90)=0.16
PARP(91)=0.40
PARP(93)=5.
PARP(99)=0.40
PARP(100)=5
PARP(102)=0.28
PARP(103)=1.0
PARP(104)=0.8
PARP(111)=2.
PARP(161)=3.00
PARP(162)=24.6
PARP(163)=18.8
PARP(164)=11.5
PARP(165)=0.47679
PARP(166)=0.67597 ! PARP165/166 are linked to MSTP17 as R_rho of HERMES is used
! PARP(166)=0.5
! ----------- Now come all the switches for Jetset -----------
PARJ(1)=0.100
PARJ(2)=0.300
PARJ(11)=0.25
PARJ(12)=0.25
PARJ(21)= 0.40
PARJ(32)=1.0
PARJ(33)= 0.80
PARJ(41)= 0.30
PARJ(42)= 0.58
PARJ(45)= 0.5
!----------------------------------------------------------------------
MSTJ(1)=1
MSTJ(12)=1
MSTJ(45)=5
MSTU(16)=2
MSTU(112)=5
MSTU(113)=5
MSTU(114)=5
! ----------- Now all the CKINs for pythia ----------
CKIN(1)=1.
CKIN(2)=-1.
CKIN(3)=0.
CKIN(4)=-1.
CKIN(5)=1.00
CKIN(6)=1.00
CKIN(7)=-10.
CKIN(8)=10.
CKIN(9)=-40.
CKIN(10)=40.
CKIN(11)=-40.
CKIN(12)=40.
CKIN(13)=-40.
CKIN(14)=40.
CKIN(15)=-40.
CKIN(16)=40.
CKIN(17)=-1.
CKIN(18)=1.
CKIN(19)=-1.
CKIN(20)=1.
CKIN(21)=0.
CKIN(22)=1.
CKIN(23)=0.
CKIN(24)=1.
CKIN(25)=-1.
CKIN(26)=1.
CKIN(27)=-1.
CKIN(28)=1.
CKIN(31)=2.
CKIN(32)=-1.
CKIN(35)=0.
CKIN(36)=-1
CKIN(37)=0.
CKIN(38)=-1.
CKIN(39)=4.
CKIN(40)=-1.
CKIN(65)=1.e-09 ! Min for Q^2
CKIN(66)=-1. ! Max for Q^2
CKIN(67)=0.
CKIN(68)=-1.
CKIN(77)=2.0
CKIN(78)=-1.