Analysis Software
Documentation for sPHENIX simulation software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
pytune.f
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file pytune.f
1 C*********************************************************************
2 
3 C...PYTUNE
4 C...Presets for a few specific underlying-event and min-bias tunes
5 C...Note some tunes require external pdfs to be linked (e.g. 105:QW),
6 C...others require particular versions of pythia (e.g. the SCI and GAL
7 C...models). See below for details.
8  SUBROUTINE pytune(ITUNE)
9 C
10 C ITUNE NAME (detailed descriptions below)
11 C 0 Default : No settings changed => linked Pythia version's defaults.
12 C ====== Old UE, Q2-ordered showers ==========================================
13 C 100 A : Rick Field's CDF Tune A
14 C 101 AW : Rick Field's CDF Tune AW
15 C 102 BW : Rick Field's CDF Tune BW
16 C 103 DW : Rick Field's CDF Tune DW
17 C 104 DWT : Rick Field's CDF Tune DW with slower UE energy scaling
18 C 105 QW : Rick Field's CDF Tune QW (NB: needs CTEQ6.1M pdfs externally)
19 C 106 ATLAS-DC2: Arthur Moraes' (old) ATLAS tune (ATLAS DC2 / Rome)
20 C 107 ACR : Tune A modified with annealing CR
21 C 108 D6 : Rick Field's CDF Tune D6 (NB: needs CTEQ6L pdfs externally)
22 C 109 D6T : Rick Field's CDF Tune D6T (NB: needs CTEQ6L pdfs externally)
23 C ====== Intermediate Models =================================================
24 C 200 IM 1 : Intermediate model: new UE, Q2-ordered showers, annealing CR
25 C 201 APT : Tune A modified to use pT-ordered final-state showers
26 C ====== New UE, interleaved pT-ordered showers, annealing CR ================
27 C 300 S0 : Sandhoff-Skands Tune 0
28 C 301 S1 : Sandhoff-Skands Tune 1
29 C 302 S2 : Sandhoff-Skands Tune 2
30 C 303 S0A : S0 with "Tune A" UE energy scaling
31 C 304 NOCR : New UE "best try" without colour reconnections
32 C 305 Old : New UE, original (primitive) colour reconnections
33 C 306 ATLAS-CSC: Arthur Moraes' (new) ATLAS tune (needs CTEQ6L externally)
34 C ======= The Uppsala models =================================================
35 C ( NB! must be run with special modified Pythia 6.215 version )
36 C ( available from http://www.isv.uu.se/thep/MC/scigal/ )
37 C 400 GAL 0 : Generalized area-law model. Old parameters
38 C 401 SCI 0 : Soft-Colour-Interaction model. Old parameters
39 C 402 GAL 1 : Generalized area-law model. Tevatron MB retuned (Skands)
40 C 403 SCI 1 : Soft-Colour-Interaction model. Tevatron MB retuned (Skands)
41 C
42 C More details;
43 C
44 C Quick Dictionary:
45 C BE : Bose-Einstein
46 C BR : Beam Remnants
47 C CR : Colour Reconnections
48 C HAD: Hadronization
49 C ISR/FSR: Initial-State Radiation / Final-State Radiation
50 C FSI: Final-State Interactions (=CR+BE)
51 C MB : Minimum-bias
52 C MI : Multiple Interactions
53 C UE : Underlying Event
54 C
55 C A (100) and AW (101). Old UE model, Q2-ordered showers.
56 C...*** NB : SHOULD BE RUN WITH PYTHIA 6.2 (e.g. 6.228) ***
57 C...*** CAN ALSO BE RUN WITH PYTHIA 6.406+
58 C...Key feature: extensively compared to CDF data (R.D. Field).
59 C...* Large starting scale for ISR (PARP(67)=4)
60 C...* AW has even more radiation due to smaller mu_R choice in alpha_s.
61 C...* See: http://www.phys.ufl.edu/~rfield/cdf/
62 C
63 C BW (102). Old UE model, Q2-ordered showers.
64 C...*** NB : SHOULD BE RUN WITH PYTHIA 6.2 (e.g. 6.228) ***
65 C...*** CAN ALSO BE RUN WITH PYTHIA 6.406+
66 C...Key feature: extensively compared to CDF data (R.D. Field).
67 C...NB: Can also be run with Pythia 6.2 or 6.312+
68 C...* Small starting scale for ISR (PARP(67)=1)
69 C...* BW has more radiation due to smaller mu_R choice in alpha_s.
70 C...* See: http://www.phys.ufl.edu/~rfield/cdf/
71 C
72 C DW (103) and DWT (104). Old UE model, Q2-ordered showers.
73 C...*** NB : SHOULD BE RUN WITH PYTHIA 6.2 (e.g. 6.228) ***
74 C...*** CAN ALSO BE RUN WITH PYTHIA 6.406+
75 C...Key feature: extensively compared to CDF data (R.D. Field).
76 C...NB: Can also be run with Pythia 6.2 or 6.312+
77 C...* Intermediate starting scale for ISR (PARP(67)=2.5)
78 C...* DWT has a different reference energy, the same as the "S" models
79 C... below, leading to more UE activity at the LHC, but less at RHIC.
80 C...* See: http://www.phys.ufl.edu/~rfield/cdf/
81 C
82 C QW (105). Old UE model, Q2-ordered showers.
83 C...*** NB : SHOULD BE RUN WITH PYTHIA 6.2 (e.g. 6.228) ***
84 C...*** CAN ALSO BE RUN WITH PYTHIA 6.406+
85 C...Key feature: uses CTEQ61 (external pdf library must be linked)
86 C
87 C ATLAS-DC2 (106). Old UE model, Q2-ordered showers.
88 C...*** NB : SHOULD BE RUN WITH PYTHIA 6.2 (e.g. 6.228) ***
89 C...*** CAN ALSO BE RUN WITH PYTHIA 6.406+
90 C...Key feature: tune used by the ATLAS collaboration.
91 C
92 C ACR (107). Old UE model, Q2-ordered showers, annealing CR.
93 C...*** NB : SHOULD BE RUN WITH PYTHIA 6.412+ ***
94 C...Key feature: Tune A modified to use annealing CR.
95 C...NB: PARP(85)=0D0 and amount of CR is regulated by PARP(78).
96 C
97 C D6 (108) and D6T (109). Old UE model, Q2-ordered showers, CTEQ6L PDF.
98 C...Key feature: Like DW and DWT but retuned to use CTEQ6L PDFs.
99 C
100 C...IM1 (200). Intermediate model, Q2-ordered showers.
101 C...Key feature: new UE model with Q2-ordered showers and no interleaving.
102 C...* "Rap" tune of hep-ph/0402078, modified with new annealing CR.
103 C...* See: Sjostrand & Skands: JHEP 03(2004)053, hep-ph/0402078.
104 C
105 C...APT (201). Old UE model, pT-ordered final-state showers
106 C...Key feature: Rick Field's Tune A, but with new final-state showers
107 C
108 C S0 (300) and S0A (303). New UE model, pT-ordered showers.
109 C...Key feature: large amount of multiple interactions
110 C...* Somewhat faster than the other colour annealing scenarios.
111 C...* S0A has a faster energy scaling of the UE IR cutoff, borrowed
112 C... from Tune A, leading to less UE at the LHC, but more at RHIC.
113 C...* Small amount of radiation.
114 C...* Large amount of low-pT MI
115 C...* Low degree of proton lumpiness (broad matter dist.)
116 C...* CR Type S (driven by free triplets), of medium strength.
117 C...* See: Pythia6402 update notes or later.
118 C
119 C S1 (301). New UE model, pT-ordered showers.
120 C...Key feature: large amount of radiation.
121 C...* Large amount of low-pT perturbative ISR
122 C...* Large amount of FSR off ISR partons
123 C...* Small amount of low-pT multiple interactions
124 C...* Moderate degree of proton lumpiness
125 C...* Least aggressive CR type (S+S Type I), but with large strength
126 C...* See: Sandhoff & Skands: FERMILAB-CONF-05-518-T, in hep-ph/0604120.
127 C
128 C S2 (302). New UE model, pT-ordered showers.
129 C...Key feature: very lumpy proton + gg string cluster formation allowed
130 C...* Small amount of radiation
131 C...* Moderate amount of low-pT MI
132 C...* High degree of proton lumpiness (more spiky matter distribution)
133 C...* Most aggressive CR type (S+S Type II), but with small strength
134 C...* See: Sandhoff & Skands: FERMILAB-CONF-05-518-T, in hep-ph/0604120.
135 C
136 C NOCR (304). New UE model, pT-ordered showers.
137 C...Key feature: no colour reconnections (NB: "Best fit" only).
138 C...* NB: <pT>(Nch) problematic in this tune.
139 C...* Small amount of radiation
140 C...* Small amount of low-pT MI
141 C...* Low degree of proton lumpiness
142 C...* Large BR composite x enhancement factor
143 C...* Most clever colour flow without CR ("Lambda ordering")
144 C
145 C ATLAS-CSC (306). New UE mode, pT-ordered showers, CTEQ6L.
146 C...Key feature: 11-parameter ATLAS tune of the new framework.
147 C...* Old (pre-annealing) colour reconnections a la 305.
148 C...* Uses CTEQ6 Leading Order PDFs (must be interfaced externally)
149 C
150 C...The GAL and SCI models (400+) are special and *SHOULD NOT* be run
151 C...with an unmodified Pythia distribution.
152 C...See http://www.isv.uu.se/thep/MC/scigal/ for more information.
153 C
154 C ::: + Future improvements?
155 C Include also QCD K-factor a la M. Heinz / ATLAS TDR ? RDF's QK?
156 C (problem: K-factor affects everything so only works as
157 C intended for min-bias, not for UE ... probably need a
158 C better long-term solution to handle UE as well. Anyway,
159 C Mark uses MSTP(33) and PARP(31)-PARP(33).)
160 
161 C...Global statements
162  IMPLICIT DOUBLE PRECISION(a-h, o-z)
163  INTEGER pyk,pychge,pycomp
164 
165 C...Commonblocks.
166  common/pydat1/mstu(200),paru(200),mstj(200),parj(200)
167  common/pypars/mstp(200),parp(200),msti(200),pari(200)
168 
169 C...SCI and GAL Commonblocks
170  COMMON /scipar/mswi(2),parsci(2)
171 
172 C...Internal parameters
173  parameter(mxtuns=500)
174  CHARACTER*8 chvers, chdoc
175  parameter(chvers='1.012 ',chdoc='Sep 2007')
176  CHARACTER*16 chnams(0:mxtuns), chname
177  CHARACTER*42 chmstj(50), chmstp(51:100), chparp(61:100),
178  & chparj(41:100), ch40
179  CHARACTER*60 ch60
180  CHARACTER*70 ch70
181  DATA (chnams(i),i=0,1)/'Default',' '/
182  DATA (chnams(i),i=100,110)/
183  & 'Tune A','Tune AW','Tune BW','Tune DW','Tune DWT','Tune QW',
184  & 'ATLAS Tune','Tune ACR','Tune D6','Tune D6T',' '/
185  DATA (chnams(i),i=300,310)/
186  & 'Tune S0','Tune S1','Tune S2','Tune S0A','NOCR','Old',
187  5 'ATLAS-CSC Tune','Yale Tune','Yale-K Tune',2*' '/
188  DATA (chnams(i),i=200,210)/
189  & 'IM Tune 1','Tune APT',9*' '/
190  DATA (chnams(i),i=400,410)/
191  & 'GAL Tune 0','SCI Tune 0','GAL Tune 1','SCI Tune 1',7*' '/
192  DATA (chmstj(i),i=11,20)/
193  & 'HAD choice of fragmentation function(s)',4*' ',
194  & 'HAD treatment of small-mass systems',4*' '/
195  DATA (chmstj(i),i=41,50)/
196  & 'FSR type (Q2 or pT) for old framework',9*' '/
197  DATA (chmstp(i),i=51,100)/
198  5 'PDF set','PDF set internal (=1) or pdflib (=2)',8*' ',
199  6 'ISR master switch',6*' ',
200  6 'ISR phase space choice & ME corrections',' ',
201  7 'ISR IR regularization scheme',' ',
202  7 'ISR scheme for FSR off ISR',8*' ',
203  8 'UE model',
204  8 'UE hadron transverse mass distribution',5*' ',
205  8 'BR composite scheme','BR colour scheme',
206  9 'BR primordial kT compensation',
207  9 'BR primordial kT distribution',
208  9 'BR energy partitioning scheme',2*' ',
209  9 'FSI colour (re-)connection model',5*' '/
210  DATA (chparp(i),i=61,100)/
211  6 ' ','ISR IR cutoff',' ','ISR renormalization scale prefactor',
212  6 2*' ','ISR Q2max factor',3*' ',
213  7 'FSR Q2max factor for non-s-channel procs',5*' ',
214  7 'FSI colour reconnection turnoff scale',
215  7 'FSI colour reconnection strength',
216  7 'BR composite x enhancement','BR breakup suppression',
217  8 2*'UE IR cutoff at reference ecm',
218  8 2*'UE mass distribution parameter',
219  8 'UE gg colour correlated fraction','UE total gg fraction',
220  8 2*' ',
221  8 'UE IR cutoff reference ecm','UE IR cutoff ecm scaling power',
222  9 'BR primordial kT width <|kT|>',' ',
223  9 'BR primordial kT UV cutoff',7*' '/
224  DATA (chparj(i),i=41,90)/
225  4 ' ','HAD string parameter b',8*' ',
226  5 3*' ','HAD charm parameter','HAD bottom parameter',5*' ',
227  6 10*' ',10*' ',
228  8 'FSR Lambda_QCD scale','FSR IR cutoff',8*' '/
229  SAVE /pydat1/,/pypars/
230  SAVE /scipar/
231 
232 C...1) Shorthand notation
233  m13=mstu(13)
234  m11=mstu(11)
235  IF (itune.LE.mxtuns.AND.itune.GE.0) THEN
236  chname=chnams(itune)
237  IF (itune.EQ.0) goto 9999
238  ELSE
239  CALL pyerrm(9,'(PYTUNE:) Tune number > max. Using defaults.')
240  goto 9999
241  ENDIF
242 
243 C...2) Hello World
244  IF (m13.GE.1) WRITE(m11,5000) chvers, chdoc
245 
246 C...3) Tune parameters
247 
248 C=============================================================================
249 C...Tunes S0, S1, S2, S0A, NOCR, and RAP (by P. Skands)
250  IF (itune.GE.300.AND.itune.LE.305) THEN
251  IF (m13.GE.1) WRITE(m11,5010) itune, chname
252  IF (mstp(181).LE.5.OR.(mstp(181).EQ.6.AND.mstp(182).LE.405))THEN
253  CALL pyerrm(9,'(PYTUNE:) linked PYTHIA version incompatible'//
254  & ' with tune.')
255  ENDIF
256 
257 C...PDFs
258  mstp(52)=1
259  mstp(51)=7
260 C...ISR
261  parp(64)=1d0
262 C...UE on, new model.
263  mstp(81)=21
264 C...Slow IR cutoff energy scaling by default
265  parp(89)=1800d0
266  parp(90)=0.16d0
267 C...Switch off trial joinings
268  mstp(96)=0
269 C...Primordial kT cutoff
270  parp(93)=5d0
271 
272 C...S0 (300), S0A (303)
273  IF (itune.EQ.300.OR.itune.EQ.303) THEN
274  IF (m13.GE.1) THEN
275  ch60='see P. Skands & D. Wicke, hep-ph/0703081'
276  WRITE(m11,5030) ch60
277  ch60='M. Sandhoff & P. Skands, in hep-ph/0604120'
278  WRITE(m11,5030) ch60
279  ch60='and T. Sjostrand & P. Skands, EPJC39(2005)129'
280  WRITE(m11,5030) ch60
281  ENDIF
282 C...Smooth ISR, low FSR
283  mstp(70)=2
284  mstp(72)=0
285 C...pT0
286  parp(82)=1.85d0
287 C...Transverse density profile.
288  mstp(82)=5
289  parp(83)=1.6d0
290 C...Colour Reconnections
291  mstp(95)=6
292  parp(78)=0.20d0
293  parp(77)=0.0d0
294 C... Reference energy for pT0 and energy scaling pace.
295  IF (itune.EQ.303) parp(90)=0.25d0
296 C...Lambda_FSR scale.
297  parj(81)=0.14d0
298 C...FSR activity.
299  parp(71)=4d0
300 C...Rap order, Valence qq, qq x enhc, BR-g-BR supp
301  mstp(89)=1
302  mstp(88)=0
303  parp(79)=2d0
304  parp(80)=0.01d0
305 
306 C...S1 (301)
307  ELSEIF(itune.EQ.301) THEN
308  IF (m13.GE.1) THEN
309  ch60='see M. Sandhoff & P. Skands, in hep-ph/0604120'
310  WRITE(m11,5030) ch60
311  ch60='and T. Sjostrand & P. Skands, EPJC39(2005)129'
312  WRITE(m11,5030) ch60
313  ENDIF
314 C...Sharp ISR, high FSR
315  mstp(70)=0
316  mstp(72)=1
317 C...pT0
318  parp(82)=2.1d0
319 C...Colour Reconnections
320  mstp(95)=2
321  parp(78)=0.35d0
322 C...Transverse density profile.
323  mstp(82)=5
324  parp(83)=1.4d0
325 C...Lambda_FSR scale.
326  parj(81)=0.14d0
327 C...FSR activity.
328  parp(71)=4d0
329 C...Rap order, Valence qq, qq x enhc, BR-g-BR supp
330  mstp(89)=1
331  mstp(88)=0
332  parp(79)=2d0
333  parp(80)=0.01d0
334 
335 C...S2 (302)
336  ELSEIF(itune.EQ.302) THEN
337  IF (m13.GE.1) THEN
338  ch60='see M. Sandhoff & P. Skands, in hep-ph/0604120'
339  WRITE(m11,5030) ch60
340  ch60='and T. Sjostrand & P. Skands, EPJC39(2005)129'
341  WRITE(m11,5030) ch60
342  ENDIF
343 C...Smooth ISR, low FSR
344  mstp(70)=2
345  mstp(72)=0
346 C...pT0
347  parp(82)=1.9d0
348 C...Transverse density profile.
349  mstp(82)=5
350  parp(83)=1.2d0
351 C...Colour Reconnections
352  mstp(95)=4
353  parp(78)=0.15d0
354 C...Lambda_FSR scale.
355  parj(81)=0.14d0
356 C...FSR activity.
357  parp(71)=4d0
358 C...Rap order, Valence qq, qq x enhc, BR-g-BR supp
359  mstp(89)=1
360  mstp(88)=0
361  parp(79)=2d0
362  parp(80)=0.01d0
363 
364 C...NOCR (304)
365  ELSEIF(itune.EQ.304) THEN
366  IF (m13.GE.1) THEN
367  ch60='"best try" without colour reconnections'
368  WRITE(m11,5030) ch60
369  ch60='see P. Skands & D. Wicke, hep-ph/0703081'
370  WRITE(m11,5030) ch60
371  ch60='and T. Sjostrand & P. Skands, EPJC39(2005)129'
372  WRITE(m11,5030) ch60
373  ENDIF
374 C...Smooth ISR, low FSR
375  mstp(70)=2
376  mstp(72)=0
377 C...pT0
378  parp(82)=2.05d0
379 C...Transverse density profile.
380  mstp(82)=5
381  parp(83)=1.8d0
382 C...Colour Reconnections
383  mstp(95)=0
384 C...Lambda_FSR scale.
385  parj(81)=0.14d0
386 C...FSR activity.
387  parp(71)=4d0
388 C...Lambda order, Valence qq, large qq x enhc, BR-g-BR supp
389  mstp(89)=2
390  mstp(88)=0
391  parp(79)=3d0
392  parp(80)=0.01d0
393 
394 C..."Lo FSR" retune (305)
395  ELSEIF(itune.EQ.305) THEN
396  IF (m13.GE.1) THEN
397  ch60='"Lo FSR retune" with primitive colour reconnections'
398  WRITE(m11,5030) ch60
399  ch60='see T. Sjostrand & P. Skands, EPJC39(2005)129'
400  WRITE(m11,5030) ch60
401  ENDIF
402 C...Smooth ISR, low FSR
403  mstp(70)=2
404  mstp(72)=0
405 C...pT0
406  parp(82)=1.9d0
407 C...Transverse density profile.
408  mstp(82)=5
409  parp(83)=2.0d0
410 C...Colour Reconnections
411  mstp(95)=1
412  parp(78)=1.0d0
413 C...Lambda_FSR scale.
414  parj(81)=0.14d0
415 C...FSR activity.
416  parp(71)=4d0
417 C...Rap order, Valence qq, qq x enhc, BR-g-BR supp
418  mstp(89)=1
419  mstp(88)=0
420  parp(79)=2d0
421  parp(80)=0.01d0
422  ENDIF
423 C...Output
424  IF (m13.GE.1) THEN
425  WRITE(m11,5030) ' '
426  WRITE(m11,5040) 51, mstp(51), chmstp(51)
427  WRITE(m11,5040) 52, mstp(52), chmstp(52)
428  WRITE(m11,5050) 64, parp(64), chparp(64)
429  WRITE(m11,5040) 68, mstp(68), chmstp(68)
430  ch60='(Note: MSTP(68) is not explicitly (re-)set by PYTUNE)'
431  WRITE(m11,5030) ch60
432  WRITE(m11,5040) 70, mstp(70), chmstp(70)
433  WRITE(m11,5040) 72, mstp(72), chmstp(72)
434  WRITE(m11,5050) 71, parp(71), chparp(71)
435  WRITE(m11,5060) 81, parj(81), chparj(81)
436  WRITE(m11,5040) 81, mstp(81), chmstp(81)
437  WRITE(m11,5050) 82, parp(82), chparp(82)
438  WRITE(m11,5050) 89, parp(89), chparp(89)
439  WRITE(m11,5050) 90, parp(90), chparp(90)
440  WRITE(m11,5040) 82, mstp(82), chmstp(82)
441  WRITE(m11,5050) 83, parp(83), chparp(83)
442  WRITE(m11,5040) 88, mstp(88), chmstp(88)
443  WRITE(m11,5040) 89, mstp(89), chmstp(89)
444  WRITE(m11,5050) 79, parp(79), chparp(79)
445  WRITE(m11,5050) 80, parp(80), chparp(80)
446  WRITE(m11,5050) 93, parp(93), chparp(93)
447  WRITE(m11,5040) 95, mstp(95), chmstp(95)
448  WRITE(m11,5050) 78, parp(78), chparp(78)
449  ENDIF
450 
451 C=============================================================================
452 C...ATLAS-CSC 11-parameter tune (By A. Moraes)
453  ELSEIF (itune.EQ.306) THEN
454  IF (m13.GE.1) WRITE(m11,5010) itune, chname
455  IF (mstp(181).LE.5.OR.(mstp(181).EQ.6.AND.mstp(182).LE.405))THEN
456  CALL pyerrm(9,'(PYTUNE:) linked PYTHIA version incompatible'//
457  & ' with tune.')
458  ENDIF
459 
460 C...PDFs
461  mstp(52)=2
462  mstp(54)=2
463  mstp(56)=2
464  mstp(51)=10042
465  mstp(53)=10042
466  mstp(55)=10042
467 C...ISR
468 C PARP(64)=1D0
469 C...UE on, new model.
470  mstp(81)=21
471 C...Energy scaling
472  parp(89)=1800d0
473  parp(90)=0.22d0
474 C...Switch off trial joinings
475  mstp(96)=0
476 C...Primordial kT cutoff
477 
478  IF (m13.GE.1) THEN
479  ch60='see presentations by A. Moraes (ATLAS),'
480  WRITE(m11,5030) ch60
481  ch60='and T. Sjostrand & P. Skands, EPJC39(2005)129'
482  WRITE(m11,5030) ch60
483  WRITE(m11,5030) ' '
484  ch70='NB! This tune requires CTEQ6.1 pdfs to be '//
485  & 'externally linked and'
486  WRITE(m11,5035) ch70
487  ch70='MSTP(51) should be set manually according to '//
488  & 'the library used'
489  WRITE(m11,5035) ch70
490  ENDIF
491 C...Smooth ISR, low FSR
492  mstp(70)=2
493  mstp(72)=0
494 C...pT0
495  parp(82)=1.9d0
496 C...Transverse density profile.
497  mstp(82)=4
498  parp(83)=0.3d0
499  parp(84)=0.5d0
500 C...ISR & FSR in interactions after the first (default)
501  mstp(84)=1
502  mstp(85)=1
503 C...No double-counting (default)
504  mstp(86)=2
505 C...Companion quark parent gluon (1-x) power
506  mstp(87)=4
507 C...Primordial kT compensation along chaings (default = 0 : uniform)
508  mstp(90)=1
509 C...Colour Reconnections
510  mstp(95)=1
511  parp(78)=0.2d0
512 C...Lambda_FSR scale.
513  parj(81)=0.14d0
514 C...Rap order, Valence qq, qq x enhc, BR-g-BR supp
515  mstp(89)=1
516  mstp(88)=0
517 C PARP(79)=2D0
518  parp(80)=0.01d0
519 C...Peterson charm frag, and c and b hadr parameters
520  mstj(11)=3
521  parj(54)=-0.07
522  parj(55)=-0.006
523 C... Output
524  IF (m13.GE.1) THEN
525  WRITE(m11,5030) ' '
526  WRITE(m11,5040) 51, mstp(51), chmstp(51)
527  WRITE(m11,5040) 52, mstp(52), chmstp(52)
528  WRITE(m11,5050) 64, parp(64), chparp(64)
529  WRITE(m11,5040) 68, mstp(68), chmstp(68)
530  ch60='(Note: MSTP(68) is not explicitly (re-)set by PYTUNE)'
531  WRITE(m11,5030) ch60
532  WRITE(m11,5040) 70, mstp(70), chmstp(70)
533  WRITE(m11,5040) 72, mstp(72), chmstp(72)
534  WRITE(m11,5050) 71, parp(71), chparp(71)
535  WRITE(m11,5060) 81, parj(81), chparj(81)
536  WRITE(m11,5040) 81, mstp(81), chmstp(81)
537  WRITE(m11,5050) 82, parp(82), chparp(82)
538  WRITE(m11,5050) 89, parp(89), chparp(89)
539  WRITE(m11,5050) 90, parp(90), chparp(90)
540  WRITE(m11,5040) 82, mstp(82), chmstp(82)
541  WRITE(m11,5050) 83, parp(83), chparp(83)
542  WRITE(m11,5050) 84, parp(84), chparp(84)
543  WRITE(m11,5040) 88, mstp(88), chmstp(88)
544  WRITE(m11,5040) 89, mstp(89), chmstp(89)
545  WRITE(m11,5040) 90, mstp(90), chmstp(90)
546  WRITE(m11,5050) 79, parp(79), chparp(79)
547  WRITE(m11,5050) 80, parp(80), chparp(80)
548  WRITE(m11,5050) 93, parp(93), chparp(93)
549  WRITE(m11,5040) 95, mstp(95), chmstp(95)
550  WRITE(m11,5050) 78, parp(78), chparp(78)
551  WRITE(m11,5070) 11, mstj(11), chmstj(11)
552  WRITE(m11,5060) 54, parj(54), chparj(54)
553  WRITE(m11,5060) 55, parj(55), chparj(55)
554  ENDIF
555 
556 C=============================================================================
557 C...Tunes A, AW, BW, DW, DWT, QW, D6, D6T (by R.D. Field, CDF)
558 C...(100-105,108-109) and ATLAS-DC2 Tune (by A. Moraes, ATLAS) (106)
559  ELSEIF ((itune.GE.100.AND.itune.LE.106).OR.itune.EQ.108.OR.
560  & itune.EQ.109) THEN
561  IF (m13.GE.1.AND.itune.NE.106) THEN
562  WRITE(m11,5010) itune, chname
563  ch60='see R.D. Field (CDF), in hep-ph/0610012'
564  WRITE(m11,5030) ch60
565  ch60='and T. Sjostrand & M. v. Zijl, PRD36(1987)2019'
566  WRITE(m11,5030) ch60
567  ENDIF
568 C...Multiple interactions on, old framework
569  mstp(81)=1
570 C...Fast IR cutoff energy scaling by default
571  parp(89)=1800d0
572  parp(90)=0.25d0
573 C...Default CTEQ5L (internal), except for QW: CTEQ61 (external)
574  mstp(51)=7
575  mstp(52)=1
576  IF (itune.EQ.105) THEN
577  mstp(51)=10150
578  mstp(52)=2
579  ELSEIF(itune.EQ.108.OR.itune.EQ.109) THEN
580  mstp(52)=2
581  mstp(54)=2
582  mstp(56)=2
583  mstp(51)=10042
584  mstp(53)=10042
585  mstp(55)=10042
586  ENDIF
587 C...Double Gaussian matter distribution.
588  mstp(82)=4
589  parp(83)=0.5d0
590  parp(84)=0.4d0
591 C...FSR activity.
592  parp(71)=4d0
593 C...Lambda_FSR scale.
594  parj(81)=0.29d0
595 C...Fragmentation functions and c and b parameters
596  mstj(11)=4
597  parj(54)=-0.05
598  parj(55)=-0.005
599 
600 C...Tune A and AW
601  IF(itune.EQ.100.OR.itune.EQ.101) THEN
602 C...pT0.
603  parp(82)=2.0d0
604 c...String drawing almost completely minimizes string length.
605  parp(85)=0.9d0
606  parp(86)=0.95d0
607 C...ISR cutoff, muR scale factor, and phase space size
608  parp(62)=1d0
609  parp(64)=1d0
610  parp(67)=4d0
611 C...Intrinsic kT, size, and max
612  mstp(91)=1
613  parp(91)=1d0
614  parp(93)=5d0
615 C...AW : higher ISR IR cutoff, but also larger alpha_s and more intrinsic kT.
616  IF (itune.EQ.101) THEN
617  parp(62)=1.25d0
618  parp(64)=0.2d0
619  parp(91)=2.1d0
620  parp(92)=15.0d0
621  ENDIF
622 
623 C...Tune BW (larger alpha_s, more intrinsic kT. Smaller ISR phase space.)
624  ELSEIF (itune.EQ.102) THEN
625 C...pT0.
626  parp(82)=1.9d0
627 c...String drawing completely minimizes string length.
628  parp(85)=1.0d0
629  parp(86)=1.0d0
630 C...ISR cutoff, muR scale factor, and phase space size
631  parp(62)=1.25d0
632  parp(64)=0.2d0
633  parp(67)=1d0
634 C...Intrinsic kT, size, and max
635  mstp(91)=1
636  parp(91)=2.1d0
637  parp(93)=15d0
638 
639 C...Tune DW
640  ELSEIF (itune.EQ.103) THEN
641 C...pT0.
642  parp(82)=1.9d0
643 c...String drawing completely minimizes string length.
644  parp(85)=1.0d0
645  parp(86)=1.0d0
646 C...ISR cutoff, muR scale factor, and phase space size
647  parp(62)=1.25d0
648  parp(64)=0.2d0
649  parp(67)=2.5d0
650 C...Intrinsic kT, size, and max
651  mstp(91)=1
652  parp(91)=2.1d0
653  parp(93)=15d0
654 
655 C...Tune DWT
656  ELSEIF (itune.EQ.104) THEN
657 C...pT0.
658  parp(82)=1.9409d0
659 C...Run II ref scale and slow scaling
660  parp(89)=1960d0
661  parp(90)=0.16d0
662 c...String drawing completely minimizes string length.
663  parp(85)=1.0d0
664  parp(86)=1.0d0
665 C...ISR cutoff, muR scale factor, and phase space size
666  parp(62)=1.25d0
667  parp(64)=0.2d0
668  parp(67)=2.5d0
669 C...Intrinsic kT, size, and max
670  mstp(91)=1
671  parp(91)=2.1d0
672  parp(93)=15d0
673 
674 C...Tune QW
675  ELSEIF(itune.EQ.105) THEN
676  IF (m13.GE.1) THEN
677  WRITE(m11,5030) ' '
678  ch70='NB! This tune requires CTEQ6.1 pdfs to be '//
679  & 'externally linked and'
680  WRITE(m11,5035) ch70
681  ch70='MSTP(51) should be set manually according to '//
682  & 'the library used'
683  WRITE(m11,5035) ch70
684  ENDIF
685 C...pT0.
686  parp(82)=1.1d0
687 c...String drawing completely minimizes string length.
688  parp(85)=1.0d0
689  parp(86)=1.0d0
690 C...ISR cutoff, muR scale factor, and phase space size
691  parp(62)=1.25d0
692  parp(64)=0.2d0
693  parp(67)=2.5d0
694 C...Intrinsic kT, size, and max
695  mstp(91)=1
696  parp(91)=2.1d0
697  parp(93)=15d0
698 
699 C...Tune D6 and D6T
700  ELSEIF(itune.EQ.108.OR.itune.EQ.109) THEN
701  IF (m13.GE.1) THEN
702  WRITE(m11,5030) ' '
703  ch70='NB! This tune requires CTEQ6L pdfs to be '//
704  & 'externally linked and'
705  WRITE(m11,5035) ch70
706  ch70='MSTP(51) should be set manually according to '//
707  & 'the library used'
708  WRITE(m11,5035) ch70
709  ENDIF
710 C...The "Rick" proton, double gauss with 0.5/0.4
711  mstp(82)=4
712  parp(83)=0.5d0
713  parp(84)=0.4d0
714 c...String drawing completely minimizes string length.
715  parp(85)=1.0d0
716  parp(86)=1.0d0
717  IF (itune.EQ.108) THEN
718 C...D6: pT0, Run I ref scale, and fast energy scaling
719  parp(82)=1.8d0
720  parp(89)=1800d0
721  parp(90)=0.25d0
722  ELSE
723 C...D6T: pT0, Run II ref scale, and slow energy scaling
724  parp(82)=1.8387d0
725  parp(89)=1960d0
726  parp(90)=0.16d0
727  ENDIF
728 C...ISR cutoff, muR scale factor, and phase space size
729  parp(62)=1.25d0
730  parp(64)=0.2d0
731  parp(67)=2.5d0
732 C...Intrinsic kT, size, and max
733  mstp(91)=1
734  parp(91)=2.1d0
735  parp(93)=15d0
736 
737 C...Old ATLAS-DC2 5-parameter tune
738  ELSEIF(itune.EQ.106) THEN
739  IF (m13.GE.1) THEN
740  WRITE(m11,5010) itune, chname
741  ch60='see A. Moraes et al., SN-ATLAS-2006-057'
742  WRITE(m11,5030) ch60
743  ch60='and T. Sjostrand & M. v. Zijl, PRD36(1987)2019'
744  WRITE(m11,5030) ch60
745  ENDIF
746 C... pT0.
747  parp(82)=1.8d0
748 C... Different ref and rescaling pacee
749  parp(89)=1000d0
750  parp(90)=0.16d0
751 C... Parameters of mass distribution
752  parp(83)=0.5d0
753  parp(84)=0.5d0
754 C... Old default string drawing
755  parp(85)=0.33d0
756  parp(86)=0.66d0
757 C... ISR, phase space equivalent to Tune B
758  parp(62)=1d0
759  parp(64)=1d0
760  parp(67)=1d0
761 C... FSR
762  parp(71)=4d0
763  parj(81)=0.29d0
764 C... Intrinsic kT
765  mstp(91)=1
766  parp(91)=1d0
767  parp(93)=5d0
768  ENDIF
769 
770 C... Output
771  IF (m13.GE.1) THEN
772  WRITE(m11,5030) ' '
773  WRITE(m11,5040) 51, mstp(51), chmstp(51)
774  WRITE(m11,5040) 52, mstp(52), chmstp(52)
775  WRITE(m11,5050) 62, parp(62), chparp(62)
776  WRITE(m11,5050) 64, parp(64), chparp(64)
777  WRITE(m11,5050) 67, parp(67), chparp(67)
778  WRITE(m11,5040) 68, mstp(68), chmstp(68)
779  ch60='(Note: MSTP(68) is not explicitly (re-)set by PYTUNE)'
780  WRITE(m11,5030) ch60
781  WRITE(m11,5050) 71, parp(71), chparp(71)
782  WRITE(m11,5060) 81, parj(81), chparj(81)
783  WRITE(m11,5040) 81, mstp(81), chmstp(81)
784  WRITE(m11,5050) 82, parp(82), chparp(82)
785  WRITE(m11,5050) 89, parp(89), chparp(89)
786  WRITE(m11,5050) 90, parp(90), chparp(90)
787  WRITE(m11,5040) 82, mstp(82), chmstp(82)
788  WRITE(m11,5050) 83, parp(83), chparp(83)
789  WRITE(m11,5050) 84, parp(84), chparp(84)
790  WRITE(m11,5050) 85, parp(85), chparp(85)
791  WRITE(m11,5050) 86, parp(86), chparp(86)
792  WRITE(m11,5040) 91, mstp(91), chmstp(91)
793  WRITE(m11,5050) 91, parp(91), chparp(91)
794  WRITE(m11,5050) 93, parp(93), chparp(93)
795  WRITE(m11,5070) 11, mstj(11), chmstj(11)
796  WRITE(m11,5060) 54, parj(54), chparj(54)
797  WRITE(m11,5060) 55, parj(55), chparj(55)
798  ENDIF
799 
800 C=============================================================================
801 C... ACR, tune A with new CR (107)
802  ELSEIF(itune.EQ.107) THEN
803  IF (m13.GE.1) THEN
804  WRITE(m11,5010) itune, chname
805  ch60='Tune A modified with new colour reconnections'
806  WRITE(m11,5030) ch60
807  ch60='PARP(85)=0D0 and amount of CR is regulated by PARP(78)'
808  WRITE(m11,5030) ch60
809  ch60='see P. Skands & D. Wicke, hep-ph/0703081,'
810  WRITE(m11,5030) ch60
811  ch60='R.D. Field (CDF), in hep-ph/0610012 (Tune A)'
812  WRITE(m11,5030) ch60
813  ch60='and T. Sjostrand & M. v. Zijl, PRD36(1987)2019'
814  WRITE(m11,5030) ch60
815  ENDIF
816  IF (mstp(181).LE.5.OR.(mstp(181).EQ.6.AND.mstp(182).LE.406))THEN
817  CALL pyerrm(9,'(PYTUNE:) linked PYTHIA version incompatible'//
818  & ' with tune. Using defaults.')
819  goto 9998
820  ENDIF
821  mstp(81)=1
822  parp(89)=1800d0
823  parp(90)=0.25d0
824  mstp(82)=4
825  parp(83)=0.5d0
826  parp(84)=0.4d0
827  mstp(51)=7
828  mstp(52)=1
829  parp(71)=4d0
830  parj(81)=0.29d0
831  parp(82)=2.0d0
832  parp(85)=0.0d0
833  parp(86)=0.66d0
834  parp(62)=1d0
835  parp(64)=1d0
836  parp(67)=4d0
837  mstp(91)=1
838  parp(91)=1d0
839  parp(93)=5d0
840  mstp(95)=6
841  parp(78)=0.25d0
842 C...Fragmentation functions and c and b parameters
843  mstj(11)=4
844  parj(54)=-0.05
845  parj(55)=-0.005
846 C...Output
847  IF (m13.GE.1) THEN
848  WRITE(m11,5030) ' '
849  WRITE(m11,5040) 51, mstp(51), chmstp(51)
850  WRITE(m11,5040) 52, mstp(52), chmstp(52)
851  WRITE(m11,5050) 62, parp(62), chparp(62)
852  WRITE(m11,5050) 64, parp(64), chparp(64)
853  WRITE(m11,5050) 67, parp(67), chparp(67)
854  WRITE(m11,5040) 68, mstp(68), chmstp(68)
855  ch60='(Note: MSTP(68) is not explicitly (re-)set by PYTUNE)'
856  WRITE(m11,5030) ch60
857  WRITE(m11,5050) 71, parp(71), chparp(71)
858  WRITE(m11,5060) 81, parj(81), chparj(81)
859  WRITE(m11,5040) 81, mstp(81), chmstp(81)
860  WRITE(m11,5050) 82, parp(82), chparp(82)
861  WRITE(m11,5050) 89, parp(89), chparp(89)
862  WRITE(m11,5050) 90, parp(90), chparp(90)
863  WRITE(m11,5040) 82, mstp(82), chmstp(82)
864  WRITE(m11,5050) 83, parp(83), chparp(83)
865  WRITE(m11,5050) 84, parp(84), chparp(84)
866  WRITE(m11,5050) 85, parp(85), chparp(85)
867  WRITE(m11,5050) 86, parp(86), chparp(86)
868  WRITE(m11,5040) 91, mstp(91), chmstp(91)
869  WRITE(m11,5050) 91, parp(91), chparp(91)
870  WRITE(m11,5050) 93, parp(93), chparp(93)
871  WRITE(m11,5040) 95, mstp(95), chmstp(95)
872  WRITE(m11,5050) 78, parp(78), chparp(78)
873  WRITE(m11,5070) 11, mstj(11), chmstj(11)
874  WRITE(m11,5060) 54, parj(54), chparj(54)
875  WRITE(m11,5060) 55, parj(55), chparj(55)
876  ENDIF
877 
878 C=============================================================================
879 C... Intermediate model. Rap tune (retuned to post-6.406 IR factorization)
880  ELSEIF(itune.EQ.200) THEN
881  IF (m13.GE.1) THEN
882  WRITE(m11,5010) itune, chname
883  ch60='see T. Sjostrand & P. Skands, JHEP03(2004)053'
884  WRITE(m11,5030) ch60
885  ENDIF
886  IF (mstp(181).LE.5.OR.(mstp(181).EQ.6.AND.mstp(182).LE.405))THEN
887  CALL pyerrm(9,'(PYTUNE:) linked PYTHIA version incompatible'//
888  & ' with tune.')
889  ENDIF
890 C...PDF
891  mstp(51)=7
892  mstp(52)=1
893 C...ISR
894  parp(62)=1d0
895  parp(64)=1d0
896  parp(67)=4d0
897 C...FSR
898  parp(71)=4d0
899  parj(81)=0.29d0
900 C...UE
901  mstp(81)=11
902  parp(82)=2.25d0
903  parp(89)=1800d0
904  parp(90)=0.25d0
905 C... ExpOfPow(1.8) overlap profile
906  mstp(82)=5
907  parp(83)=1.8d0
908 C... Valence qq
909  mstp(88)=0
910 C... Rap Tune
911  mstp(89)=1
912 C... Default diquark, BR-g-BR supp
913  parp(79)=2d0
914  parp(80)=0.01d0
915 C... Final state reconnect.
916  mstp(95)=1
917  parp(78)=0.55d0
918 C...Fragmentation functions and c and b parameters
919  mstj(11)=4
920  parj(54)=-0.05
921  parj(55)=-0.005
922 C... Output
923  IF (m13.GE.1) THEN
924  WRITE(m11,5030) ' '
925  WRITE(m11,5040) 51, mstp(51), chmstp(51)
926  WRITE(m11,5040) 52, mstp(52), chmstp(52)
927  WRITE(m11,5050) 62, parp(62), chparp(62)
928  WRITE(m11,5050) 64, parp(64), chparp(64)
929  WRITE(m11,5050) 67, parp(67), chparp(67)
930  WRITE(m11,5040) 68, mstp(68), chmstp(68)
931  ch60='(Note: MSTP(68) is not explicitly (re-)set by PYTUNE)'
932  WRITE(m11,5030) ch60
933  WRITE(m11,5050) 71, parp(71), chparp(71)
934  WRITE(m11,5060) 81, parj(81), chparj(81)
935  WRITE(m11,5040) 81, mstp(81), chmstp(81)
936  WRITE(m11,5050) 82, parp(82), chparp(82)
937  WRITE(m11,5050) 89, parp(89), chparp(89)
938  WRITE(m11,5050) 90, parp(90), chparp(90)
939  WRITE(m11,5040) 82, mstp(82), chmstp(82)
940  WRITE(m11,5050) 83, parp(83), chparp(83)
941  WRITE(m11,5040) 88, mstp(88), chmstp(88)
942  WRITE(m11,5040) 89, mstp(89), chmstp(89)
943  WRITE(m11,5050) 79, parp(79), chparp(79)
944  WRITE(m11,5050) 80, parp(80), chparp(80)
945  WRITE(m11,5050) 93, parp(93), chparp(93)
946  WRITE(m11,5040) 95, mstp(95), chmstp(95)
947  WRITE(m11,5050) 78, parp(78), chparp(78)
948  WRITE(m11,5070) 11, mstj(11), chmstj(11)
949  WRITE(m11,5060) 54, parj(54), chparj(54)
950  WRITE(m11,5060) 55, parj(55), chparj(55)
951  ENDIF
952 
953 C...APT. Tune A modified to use new pT-ordered FSR.
954  ELSEIF(itune.EQ.201) THEN
955  IF (m13.GE.1) THEN
956  WRITE(m11,5010) itune, chname
957  ch60='see P. Skands & D. Wicke, hep-ph/0703081 (Tune APT),'
958  WRITE(m11,5030) ch60
959  ch60='R.D. Field (CDF), in hep-ph/0610012 (Tune A)'
960  WRITE(m11,5030) ch60
961  ch60='T. Sjostrand & M. v. Zijl, PRD36(1987)2019'
962  WRITE(m11,5030) ch60
963  ch60='and T. Sjostrand & P. Skands, EPJC39(2005)129'
964  WRITE(m11,5030) ch60
965  ENDIF
966  IF (mstp(181).LE.5.OR.(mstp(181).EQ.6.AND.mstp(182).LE.411))THEN
967  CALL pyerrm(9,'(PYTUNE:) linked PYTHIA version incompatible'//
968  & ' with tune.')
969  ENDIF
970 C...First set as if Pythia tune A
971 C...Multiple interactions on, old framework
972  mstp(81)=1
973 C...Fast IR cutoff energy scaling by default
974  parp(89)=1800d0
975  parp(90)=0.25d0
976 C...Default CTEQ5L (internal)
977  mstp(51)=7
978  mstp(52)=1
979 C...Double Gaussian matter distribution.
980  mstp(82)=4
981  parp(83)=0.5d0
982  parp(84)=0.4d0
983 C...FSR activity.
984  parp(71)=4d0
985 c...String drawing almost completely minimizes string length.
986  parp(85)=0.9d0
987  parp(86)=0.95d0
988 C...ISR cutoff, muR scale factor, and phase space size
989  parp(62)=1d0
990  parp(64)=1d0
991  parp(67)=4d0
992 C...Intrinsic kT, size, and max
993  mstp(91)=1
994  parp(91)=1d0
995  parp(93)=5d0
996 C...Use pT-ordered FSR
997  mstj(41)=12
998 C...Lambda_FSR scale for pT-ordering
999  parj(81)=0.14d0
1000 C...Retune pT0
1001  parp(82)=2.1d0
1002 C...Fragmentation functions and c and b parameters
1003  mstj(11)=4
1004  parj(54)=-0.05
1005  parj(55)=-0.005
1006 
1007 C... Output
1008  IF (m13.GE.1) THEN
1009  WRITE(m11,5030) ' '
1010  WRITE(m11,5040) 51, mstp(51), chmstp(51)
1011  WRITE(m11,5040) 52, mstp(52), chmstp(52)
1012  WRITE(m11,5050) 62, parp(62), chparp(62)
1013  WRITE(m11,5050) 64, parp(64), chparp(64)
1014  WRITE(m11,5050) 67, parp(67), chparp(67)
1015  WRITE(m11,5040) 68, mstp(68), chmstp(68)
1016  ch60='(Note: MSTP(68) is not explicitly (re-)set by PYTUNE)'
1017  WRITE(m11,5030) ch60
1018  WRITE(m11,5070) 41, mstj(41), chmstj(41)
1019  WRITE(m11,5050) 71, parp(71), chparp(71)
1020  WRITE(m11,5060) 81, parj(81), chparj(81)
1021  WRITE(m11,5040) 81, mstp(81), chmstp(81)
1022  WRITE(m11,5050) 82, parp(82), chparp(82)
1023  WRITE(m11,5050) 89, parp(89), chparp(89)
1024  WRITE(m11,5050) 90, parp(90), chparp(90)
1025  WRITE(m11,5040) 82, mstp(82), chmstp(82)
1026  WRITE(m11,5050) 83, parp(83), chparp(83)
1027  WRITE(m11,5050) 84, parp(84), chparp(84)
1028  WRITE(m11,5050) 85, parp(85), chparp(85)
1029  WRITE(m11,5050) 86, parp(86), chparp(86)
1030  WRITE(m11,5040) 91, mstp(91), chmstp(91)
1031  WRITE(m11,5050) 91, parp(91), chparp(91)
1032  WRITE(m11,5050) 93, parp(93), chparp(93)
1033  WRITE(m11,5070) 11, mstj(11), chmstj(11)
1034  WRITE(m11,5060) 54, parj(54), chparj(54)
1035  WRITE(m11,5060) 55, parj(55), chparj(55)
1036  ENDIF
1037 
1038 C=============================================================================
1039 C...Uppsala models: Generalized Area Law and Soft Colour Interactions
1040  ELSEIF(chname.EQ.'GAL Tune 0'.OR.chname.EQ.'GAL Tune 1') THEN
1041  IF (m13.GE.1) THEN
1042  WRITE(m11,5010) itune, chname
1043  ch60='see J. Rathsman, PLB452(1999)364'
1044  WRITE(m11,5030) ch60
1045 C ? CH60='A. Edin, G. Ingelman, J. Rathsman, hep-ph/9912539,'
1046 C ? WRITE(M11,5030)
1047  ch60='and T. Sjostrand & M. v. Zijl, PRD36(1987)2019'
1048  WRITE(m11,5030) ch60
1049  WRITE(m11,5030) ' '
1050  ch70='NB! The GAL model must be run with modified '//
1051  & 'Pythia v6.215:'
1052  WRITE(m11,5035) ch70
1053  ch70='available from http://www.isv.uu.se/thep/MC/scigal/'
1054  WRITE(m11,5035) ch70
1055  WRITE(m11,5030) ' '
1056  ENDIF
1057 C...GAL Recommended settings from Uppsala web page (as per 22/08 2006)
1058  mswi(2) = 3
1059  parsci(2) = 0.10
1060  mswi(1) = 2
1061  parsci(1) = 0.44
1062  mstj(16) = 0
1063  parj(42) = 0.45
1064  parj(82) = 2.0
1065  parp(62) = 2.0
1066  mstp(81) = 1
1067  mstp(82) = 1
1068  parp(81) = 1.9
1069  mstp(92) = 1
1070  IF(chname.EQ.'GAL Tune 1') THEN
1071 C...GAL retune (P. Skands) to get better min-bias <Nch> at Tevatron
1072  mstp(82)=4
1073  parp(83)=0.25d0
1074  parp(84)=0.5d0
1075  parp(82) = 1.75
1076  IF (m13.GE.1) THEN
1077  WRITE(m11,5040) 81, mstp(81), chmstp(81)
1078  WRITE(m11,5050) 82, parp(82), chparp(82)
1079  WRITE(m11,5040) 82, mstp(82), chmstp(82)
1080  WRITE(m11,5050) 83, parp(83), chparp(83)
1081  WRITE(m11,5050) 84, parp(84), chparp(84)
1082  ENDIF
1083  ELSE
1084  IF (m13.GE.1) THEN
1085  WRITE(m11,5040) 81, mstp(81), chmstp(81)
1086  WRITE(m11,5050) 81, parp(81), chparp(81)
1087  WRITE(m11,5040) 82, mstp(82), chmstp(82)
1088  ENDIF
1089  ENDIF
1090 C...Output
1091  IF (m13.GE.1) THEN
1092  WRITE(m11,5050) 62, parp(62), chparp(62)
1093  WRITE(m11,5060) 82, parj(82), chparj(82)
1094  WRITE(m11,5040) 92, mstp(92), chmstp(92)
1095  ch40='FSI SCI/GAL selection'
1096  WRITE(m11,6040) 1, mswi(1), ch40
1097  ch40='FSI SCI/GAL sea quark treatment'
1098  WRITE(m11,6040) 2, mswi(2), ch40
1099  ch40='FSI SCI/GAL sea quark treatment parm'
1100  WRITE(m11,6050) 1, parsci(1), ch40
1101  ch40='FSI SCI/GAL string reco probability R_0'
1102  WRITE(m11,6050) 2, parsci(2), ch40
1103  WRITE(m11,5060) 42, parj(42), chparj(42)
1104  WRITE(m11,5070) 16, mstj(16), chmstj(16)
1105  ENDIF
1106  ELSEIF(chname.EQ.'SCI Tune 0'.OR.chname.EQ.'SCI Tune 1') THEN
1107  IF (m13.GE.1) THEN
1108  WRITE(m11,5010) itune, chname
1109  ch60='see A.Edin et al, PLB366(1996)371, Z.Phys.C75(1997)57,'
1110  WRITE(m11,5030) ch60
1111  ch60='and T. Sjostrand & M. v. Zijl, PRD36(1987)2019'
1112  WRITE(m11,5030) ch60
1113  WRITE(m11,5030) ' '
1114  ch70='NB! The SCI model must be run with modified '//
1115  & 'Pythia v6.215:'
1116  WRITE(m11,5035) ch70
1117  ch70='available from http://www.isv.uu.se/thep/MC/scigal/'
1118  WRITE(m11,5035) ch70
1119  WRITE(m11,5030) ' '
1120  ENDIF
1121 C...SCI Recommended settings from Uppsala web page (as per 22/08 2006)
1122  mstp(81)=1
1123  mstp(82)=1
1124  parp(81)=2.2
1125  mstp(92)=1
1126  mswi(2)=2
1127  parsci(2)=0.50
1128  mswi(1)=2
1129  parsci(1)=0.44
1130  mstj(16)=0
1131  IF (chname.EQ.'SCI Tune 1') THEN
1132 C...SCI retune (P. Skands) to get better min-bias <Nch> at Tevatron
1133  mstp(81) = 1
1134  mstp(82) = 3
1135  parp(82) = 2.4
1136  parp(83) = 0.5d0
1137  parp(62) = 1.5
1138  parp(84)=0.25d0
1139  IF (m13.GE.1) THEN
1140  WRITE(m11,5040) 81, mstp(81), chmstp(81)
1141  WRITE(m11,5050) 82, parp(82), chparp(82)
1142  WRITE(m11,5040) 82, mstp(82), chmstp(82)
1143  WRITE(m11,5050) 83, parp(83), chparp(83)
1144  WRITE(m11,5050) 62, parp(62), chparp(62)
1145  ENDIF
1146  ELSE
1147  IF (m13.GE.1) THEN
1148  WRITE(m11,5040) 81, mstp(81), chmstp(81)
1149  WRITE(m11,5050) 81, parp(81), chparp(81)
1150  WRITE(m11,5040) 82, mstp(82), chmstp(82)
1151  ENDIF
1152  ENDIF
1153 C...Output
1154  IF (m13.GE.1) THEN
1155  WRITE(m11,5040) 92, mstp(92), chmstp(92)
1156  ch40='FSI SCI/GAL selection'
1157  WRITE(m11,6040) 1, mswi(1), ch40
1158  ch40='FSI SCI/GAL sea quark treatment'
1159  WRITE(m11,6040) 2, mswi(2), ch40
1160  ch40='FSI SCI/GAL sea quark treatment parm'
1161  WRITE(m11,6050) 1, parsci(1), ch40
1162  ch40='FSI SCI/GAL string reco probability R_0'
1163  WRITE(m11,6050) 2, parsci(2), ch40
1164  WRITE(m11,5070) 16, mstj(16), chmstj(16)
1165  ENDIF
1166 
1167  ELSE
1168  IF (mstu(13).GE.1) WRITE(m11,5020) itune
1169 
1170  ENDIF
1171 
1172  9998 IF (mstu(13).GE.1) WRITE(m11,6000)
1173 
1174  9999 RETURN
1175 
1176  5000 FORMAT(1x,78('*')/' *',76x,'*'/' *',3x,'PYTUNE v',a6,' : ',
1177  & 'Presets for underlying-event (and min-bias)',13x,'*'/' *',
1178  & 20x,'Last Change : ',a8,' - P. Skands',22x,'*'/' *',76x,'*')
1179  5010 FORMAT(' *',3x,i4,1x,a16,52x,'*')
1180  5020 FORMAT(' *',3x,'Tune ',i4, ' not recognized. Using defaults.')
1181  5030 FORMAT(' *',3x,10x,a60,3x,'*')
1182  5035 FORMAT(' *',3x,a70,3x,'*')
1183  5040 FORMAT(' *',5x,'MSTP(',i2,') = ',i12,3x,a42,3x,'*')
1184  5050 FORMAT(' *',5x,'PARP(',i2,') = ',f12.4,3x,a40,5x,'*')
1185  5060 FORMAT(' *',5x,'PARJ(',i2,') = ',f12.4,3x,a40,5x,'*')
1186  5070 FORMAT(' *',5x,'MSTJ(',i2,') = ',i12,3x,a40,5x,'*')
1187  5140 FORMAT(' *',5x,'MSTP(',i3,')= ',i12,3x,a40,5x,'*')
1188  5150 FORMAT(' *',5x,'PARP(',i3,')= ',f12.4,3x,a40,5x,'*')
1189  6000 FORMAT(' *',76x,'*'/1x,32('*'),1x,'END OF PYTUNE',1x,31('*'))
1190  6040 FORMAT(' *',5x,'MSWI(',i1,') = ',i12,3x,a40,5x,'*')
1191  6050 FORMAT(' *',5x,'PARSCI(',i1,')= ',f12.4,3x,a40,5x,'*')
1192 
1193  END