]>
Commit | Line | Data |
---|---|---|
d30b8254 | 1 | |
2 | +-------------------------------------------------------------+ | |
3 | | | | |
4 | | | | |
5 | | DPMJET 3.0 | | |
6 | | | | |
7 | | | | |
8 | | S. Roesler+), R. Engel#), J. Ranft*) | | |
9 | | | | |
10 | | +) CERN, TIS-RP | | |
11 | | CH-1211 Geneva 23, Switzerland | | |
12 | | Email: Stefan.Roesler@cern.ch | | |
13 | | | | |
14 | | #) University of Delaware, BRI | | |
15 | | Newark, DE 19716, USA | | |
16 | | | | |
17 | | *) University of Siegen, Dept. of Physics | | |
18 | | D-57068 Siegen, Germany | | |
19 | | | | |
20 | | | | |
21 | | http://home.cern.ch/sroesler/dpmjet3.html | | |
22 | | | | |
23 | | | | |
24 | | Monte Carlo models used for event generation: | | |
25 | | PHOJET 1.12, PYTHIA 6.115 and LEPTO 6.5.1 | | |
26 | | | | |
27 | +-------------------------------------------------------------+ | |
28 | ||
29 | ||
30 | ||
31 | INPUT OPTIONS | |
32 | _______________ | |
33 | ||
34 | ||
35 | The input of DPMJET consists of option cards. Option cards have all the | |
36 | same structure and have to be given in fixed format except for the | |
37 | section enclosed by PHOINPUT and ENDINPUT which can be given in free format. | |
38 | ||
39 | CODEWD, (WHAT(I),I=1,6), SDUM (default fixed format is | |
40 | A10,6E10.0,A8) | |
41 | ||
42 | where: | |
43 | ||
44 | - CODEWD is the option keyword | |
45 | - The WHAT-parameters are numerical data | |
46 | - SDUM may contain literal data | |
47 | ||
48 | The order of the input cards is free, with the exception of the START card | |
49 | which initiates event generation and the BEAM card (see below). | |
50 | ||
51 | Most definitions have some default values. If these are acceptable, | |
52 | it is not compulsory that the corresponding option card appear | |
53 | explicitly in the input sequence. Similarly, for most WHAT and/or SDUM | |
54 | parameters a default value is applied if the corresponding field is left | |
55 | blank (or set = 0.0) in the input card. | |
56 | ||
57 | Several option cards may appear more than once in the input sequence. | |
58 | In most cases, each of such additional cards obviously adds more | |
59 | definitions to those already given, provided they are different and not | |
60 | contradictory. In case of conflict, the last given generally overrides | |
61 | the previous one(s). | |
62 | ||
63 | Any line starting with "*" is treated as a comment line. | |
64 | ||
65 | ||
66 | 1) List of input options | |
67 | ------------------------ | |
68 | ||
69 | *** general options *** | |
70 | ||
71 | TITLE title of run | |
72 | START start of event generation | |
73 | STOP stop of the event generation | |
74 | RNDMINIT initialization of random number generator | |
75 | ||
76 | *** interacting particles *** | |
77 | ||
78 | PROJPAR projectile parameters | |
79 | TARPAR target parameters | |
80 | EMULSION definition of nuclear target emulsions or composite targets | |
81 | ||
82 | *** collision energy *** | |
83 | ||
84 | ENERGY interaction energy (per nucleon, lab) | |
85 | MOMENTUM interaction momentum (per nucleon, lab) | |
86 | CMENERGY interaction energy (nucleon-nucleon cms) | |
87 | BEAM definition of beam parameters | |
88 | ||
89 | *** model for hadron / lepton / photon - nucleon interactions *** | |
90 | ||
91 | MODEL model to be used to treat nucleon-nucleon interactions | |
92 | PHOINPUT start of PHOJET-specific input | |
93 | ENDINPUT end of PHOJET-specific input | |
94 | HADRIN HADRIN module | |
95 | LEPTO-CUT parameter CUT in LEPTO-common /LEPTOU/ | |
96 | LEPTO-LST parameter LST in LEPTO-common /LEPTOU/ | |
97 | LEPTO-PARL parameter PARL in LEPTO-common /LEPTOU/ | |
98 | ||
99 | *** Glauber formalism - cross sections *** | |
100 | ||
101 | GLAUB-PAR parameters in Glauber-formalism | |
102 | GLAUB-INI pre-initialization of profile function | |
103 | FLUCTUAT treatment of cross section fluctuations | |
104 | VDM-PAR1 parameters in gamma-nucleus cross section calculation | |
105 | VDM-PAR2 parameters in gamma-nucleus cross section calculation | |
106 | XS-TABLE output of cross section table for requested interaction | |
107 | CENTRAL biasing in impact parameter | |
108 | ||
109 | *** parameters in DPM two-chain approximation *** | |
110 | ||
111 | RECOMBIN chain recombination | |
112 | COMBIJET chain fusion | |
113 | CRONINPT Cronin multiple scattering of partons | |
114 | DIQUARKS sea-diquark/antidiquark-pairs | |
115 | ||
116 | *** hadronization and JETSET-parameters *** | |
117 | ||
118 | LUND-MSTU set parameter MSTU in JETSET-common /LUDAT1/ | |
119 | LUND-PARJ set parameter PARJ in JETSET-common /LUDAT1/ | |
120 | LUND-PARU set parameter PARJ in JETSET-common /LUDAT1/ | |
121 | POPCORN Popcorn-effect in fragmentation | |
122 | PARDECAY decay of Sigma0, Asigma0, pion^0 | |
123 | ||
124 | *** nuclear fragmentation *** | |
125 | ||
126 | FERMI Fermi momentum of nucleons | |
127 | TAUFOR formation time suppressed intranuclear cascade | |
128 | PAULI treatment of Pauli's principle | |
129 | COULOMB treatment of Coulomb force | |
130 | EVAP evaporation / fragmentation module | |
131 | ||
132 | *** output and checks *** | |
133 | ||
134 | FRAME Lorentz-frame in which final state is given in DTEVT1 | |
135 | HISTOGRAM activate different classes of histograms | |
136 | EMCCHECK extended energy-momentum / quantum-number conservation check | |
137 | ||
138 | *** lepton tagger *** | |
139 | ||
140 | L-TAG lepton tagger (lepton-nucleus interactions only) | |
141 | L-ETAG lepton tagger (lepton-nucleus interactions only) | |
142 | ECMS-CUT lepton tagger (lepton-nucleus interactions only) | |
143 | ||
144 | *** for code development only *** | |
145 | ||
146 | INTPT intrinsic transverse momenta of partons | |
147 | OUTLEVEL output control switches | |
148 | RESONANC treatment of low mass chains | |
149 | SEASU3 treatment of strange-quarks at chain ends | |
150 | XCUTS thresholds for x-sampling | |
151 | ||
152 | ||
153 | 2) Description of input options | |
154 | ------------------------------- | |
155 | ||
156 | ||
157 | *** general options *** | |
158 | ||
159 | ______________________________________________________________________ | |
160 | ||
161 | control card: codewd = TITLE | |
162 | ||
163 | what (1..6), sdum no meaning | |
164 | ||
165 | Note: The control-card following this card must consist of | |
166 | a string of characters usually giving the title of | |
167 | the run. | |
168 | ______________________________________________________________________ | |
169 | ||
170 | control card: codewd = START | |
171 | ||
172 | what (1) = number of events default: 100. | |
173 | ______________________________________________________________________ | |
174 | ||
175 | control card: codewd = STOP | |
176 | ||
177 | stop of the event generation | |
178 | ||
179 | what (1..6) no meaning | |
180 | ______________________________________________________________________ | |
181 | ||
182 | control card: RNDMINIT | |
183 | ||
184 | initialization of random number generator | |
185 | ||
186 | what (1..4) values for initialization (= 1..168) | |
187 | what (5..6), sdum no meaning | |
188 | ______________________________________________________________________ | |
189 | ||
190 | ||
191 | *** interacting particles *** | |
192 | ||
193 | ______________________________________________________________________ | |
194 | ||
195 | control card: codewd = PROJPAR | |
196 | ||
197 | definition of projectile properties | |
198 | ||
199 | what (1) = (nucleus) mass number | |
200 | (photon) virtuality Q^2 | |
201 | (lepton, PHOJET) | |
202 | maximum virtuality Q^2 of emitted photon | |
203 | (otherwise) no meaning | |
204 | what (2) = (nucleus) charge number | |
205 | (otherwise) no meaning | |
206 | what (3..6) no meaning | |
207 | sdum (hadrons,photons,leptons) particle code word | |
208 | ||
209 | Note: In general, projectile nuclei are defined by what (1) and | |
210 | what (2). All other projectiles are defined by sdum. | |
211 | ______________________________________________________________________ | |
212 | ||
213 | control card: codewd = TARPAR | |
214 | ||
215 | definition of target properties | |
216 | ||
217 | what (1) = (nucleus) mass number | |
218 | (otherwise) no meaning | |
219 | what (2) = (nucleus) charge number | |
220 | (otherwise) no meaning | |
221 | what (3..6) no meaning | |
222 | sdum (hadrons) particle code word | |
223 | ||
224 | Note: In general, target nuclei are defined by what (1) and | |
225 | what (2). Target hadrons are defined by sdum. | |
226 | ______________________________________________________________________ | |
227 | ||
228 | ||
229 | control card: codewd = EMULSION | |
230 | ||
231 | definition of nuclear target emulsions or composite targets | |
232 | ||
233 | what(1) mass number of emulsion component | |
234 | what(2) charge of emulsion component | |
235 | what(3) fraction of events with this target | |
236 | what(4,5,6) as what(1,2,3) but for a further component | |
237 | default: no emulsion | |
238 | sdum no meaning | |
239 | ||
240 | Note: If this input-card is once used with valid parameters | |
241 | TARPAR is obsolete. Not the absolute values of the fractions | |
242 | are important but only relative values. | |
243 | This control card can be repeatedly used to define | |
244 | emulsions / composite targets consisting of up to 10 | |
245 | elements. | |
246 | ______________________________________________________________________ | |
247 | ||
248 | ||
249 | *** collision energy *** | |
250 | ||
251 | ______________________________________________________________________ | |
252 | ||
253 | control card: codewd = ENERGY | |
254 | ||
255 | definition of laboratory energy | |
256 | ||
257 | what (1) > 0: what (1) = total energy per nucleon (GeV) | |
258 | < 0: |what(1)| = kinetic energy per nucleon (GeV) | |
259 | default: 200 GeV | |
260 | if |what(2)| > 0: min. total/kinetic energy per nucleon | |
261 | for variable energy runs | |
262 | what (2) max. energy per nucleon for variable energy runs | |
263 | > 0: what (2) = total energy per nucleon (GeV) | |
264 | < 0: |what(1)| = kinetic energy per nucleon (GeV) | |
265 | ______________________________________________________________________ | |
266 | ||
267 | control card: codewd = MOMENTUM | |
268 | ||
269 | definition of laboratory momentum | |
270 | ||
271 | what (1) = momentum per nucleon (GeV/c) of projectile in Lab. | |
272 | default: 200 GeV/c | |
273 | what (2..6), sdum no meaning | |
274 | ______________________________________________________________________ | |
275 | ||
276 | control card: codewd = CMENERGY | |
277 | ||
278 | what (1) = nucleon-nucleon c.m. energy default: none | |
279 | what (2..6), sdum no meaning | |
280 | ______________________________________________________________________ | |
281 | ||
282 | control card: codewd = BEAM | |
283 | ||
284 | definition of beam parameters | |
285 | ||
286 | what (1/2) > 0 : energy per nucleon of beam 1/2 (GeV) | |
287 | < 0 : |what(1/2)| = energy per charge of | |
288 | beam 1/2 (GeV) | |
289 | (beam 1 is directed into positive z-direction) | |
290 | what (3) beam crossing angle, defined as 2x angle between | |
291 | one beam and the z-axis (micro rad) | |
292 | what (4) angle with x-axis defining the collision plane | |
293 | what (5..6), sdum no meaning | |
294 | ||
295 | Note: This card requires previously defined projectile and | |
296 | target identities (PROJPAR, TARPAR) ! | |
297 | ______________________________________________________________________ | |
298 | ||
299 | ||
300 | *** model for hadron / lepton / photon - nucleon interactions *** | |
301 | ||
302 | ______________________________________________________________________ | |
303 | ||
304 | control card: codewd = MODEL | |
305 | ||
306 | Model used to describe nucleon(hadron,photon,lepton)-nucleon | |
307 | interactions | |
308 | ||
309 | what (1) (only if sdum = LEPTO) | |
310 | variable INTER (see LEPTO-manual) | |
311 | = 1 gamma exchange | |
312 | = 2 W+- exchange | |
313 | = 3 Z0 exchange | |
314 | = 4 gamma/Z0 exchange | |
315 | ||
316 | sdum = DTUNUC two-chain model as for versions 1.xx | |
317 | (nucleon/hadron-nucleon interactions only) | |
318 | = PHOJET multiple chains including minijets | |
319 | = LEPTO DIS | |
320 | default: PHOJET | |
321 | what (2..6) no meaning | |
322 | ______________________________________________________________________ | |
323 | ||
324 | control card: codewd = PHOINPUT | |
325 | ||
326 | Start of PHOJET-specific input. | |
327 | For details and a list of PHOJET input cards see the PHOJET-manual | |
328 | available at | |
329 | http://lepton.bartol.udel.edu/~eng/phojet.html | |
330 | Note: This part of the input has to be closed by the ENDINPUT-card | |
331 | ||
332 | what (1..6), sdum no meaning | |
333 | ______________________________________________________________________ | |
334 | ||
335 | control card: codewd = ENDINPUT | |
336 | ||
337 | End of PHOJET-specific input. | |
338 | what (1..6), sdum no meaning | |
339 | ______________________________________________________________________ | |
340 | ||
341 | control card: codewd = HADRIN | |
342 | ||
343 | HADRIN module | |
344 | ||
345 | what (1) = 0. elastic/inelastic interactions with probab. | |
346 | as defined by cross-sections | |
347 | = 1. inelastic interactions forced | |
348 | = 2. elastic interactions forced default: 1 | |
349 | what (2) upper threshold in total energy (GeV) below | |
350 | which interactions are sampled by HADRIN | |
351 | default: steady transition btw. HADRIN and | |
352 | DPM in the range 4-7 GeV | |
353 | what (3..6), sdum no meaning | |
354 | ______________________________________________________________________ | |
355 | ||
356 | control card: codewd = LEPTO-CUT | |
357 | ||
358 | set parameter CUT in LEPTO-common /LEPTOU/ | |
359 | ||
360 | what (1) = index in CUT-array | |
361 | what (2) = new value of CUT( int(what(1)) ) | |
362 | what (3), what(4) and what (5), what(6) further | |
363 | parameter in the same way as what (1) and what (2) | |
364 | default: default-LEPTO parameters | |
365 | ||
366 | Note: see LEPTO-manual. | |
367 | ______________________________________________________________________ | |
368 | ||
369 | control card: codewd = LEPTO-LST | |
370 | ||
371 | set parameter LST in LEPTO-common /LEPTOU/ | |
372 | ||
373 | what (1) = index in LST-array | |
374 | what (2) = new value of LST( int(what(1)) ) | |
375 | what (3), what(4) and what (5), what(6) further | |
376 | parameter in the same way as what (1) and what (2) | |
377 | default: default-LEPTO parameters | |
378 | ||
379 | Note: see LEPTO-manual. | |
380 | ______________________________________________________________________ | |
381 | ||
382 | control card: codewd = LEPTO-PARL | |
383 | ||
384 | set parameter PARL in LEPTO-common /LEPTOU/ | |
385 | ||
386 | what (1) = index in PARL-array | |
387 | what (2) = new value of PARL( int(what(1)) ) | |
388 | what (3), what(4) and what (5), what(6) further | |
389 | parameter in the same way as what (1) and what (2) | |
390 | default: default-LEPTO parameters | |
391 | ||
392 | Note: see LEPTO-manual. | |
393 | ______________________________________________________________________ | |
394 | ||
395 | ||
396 | *** Glauber formalism - cross sections *** | |
397 | ||
398 | ______________________________________________________________________ | |
399 | ||
400 | control card: codewd = GLAUB-PAR | |
401 | ||
402 | parameters in Glauber-formalism | |
403 | ||
404 | what (1) number of nucleon configurations sampled in integration | |
405 | over nuclear density default: 1000 | |
406 | what (2) number of bins for integration over impact-parameter and | |
407 | for profile-function calculation default: 49 | |
408 | what (3) = 1 calculation of tot., el. and qel. cross sections | |
409 | otherwise calculation of production cross sections only | |
410 | default: 0 | |
411 | what (4) = 1 read pre-calculated impact-parameter distribution | |
412 | from "sdum".glb for fixed projectile/target/energy | |
413 | runs | |
414 | = -1 dump calculated impact-parameter distribution | |
415 | into "sdum".glb for fixed or variable projectile/ | |
416 | target/energy runs | |
417 | = 100 read pre-calculated impact-parameter distribution | |
418 | from "sdum".glb for variable projectile/target/ | |
419 | energy runs | |
420 | default: 0 | |
421 | what (5..6) no meaning | |
422 | sdum if |what (4)| = 1 name of in/output-file (sdum.glb) | |
423 | ______________________________________________________________________ | |
424 | ||
425 | control card: codewd = GLAUB-INI | |
426 | ||
427 | pre-initialization of profile function | |
428 | ||
429 | what (1) lower energy limit for initialization | |
430 | > 0 Lab. frame | |
431 | < 0 nucleon-nucleon cms | |
432 | what (2) upper energy limit for initialization | |
433 | > 0 Lab. frame | |
434 | < 0 nucleon-nucleon cms | |
435 | what (3) > 0 # of equidistant lin. bins in E | |
436 | < 0 # of equidistant log. bins in E | |
437 | what (4) maximum projectile mass number for which the Glauber | |
438 | data are initialized for each projectile mass number | |
439 | (if <= mass given with the PROJPAR-card) | |
440 | default: 18 | |
441 | what (5) steps in mass number starting from what (4) | |
442 | up to mass number defined with PROJPAR-card | |
443 | for which Glauber data are initialized | |
444 | default: 5 | |
445 | what (6) no meaning | |
446 | sdum no meaning | |
447 | ______________________________________________________________________ | |
448 | ||
449 | control card: codewd = FLUCTUAT | |
450 | ||
451 | Treatment of cross section fluctuations | |
452 | ||
453 | what (1) = 1 cross section fluctuations treated default: 0. | |
454 | what (2..6), sdum no meaning | |
455 | ______________________________________________________________________ | |
456 | ||
457 | ||
458 | control card: codewd = VDM-PAR1 | |
459 | ||
460 | parameters in gamma-nucleus cross section calculation | |
461 | ||
462 | what (1) = Lambda^2 default: 2. | |
463 | what (2) lower limit in M^2 integration | |
464 | = 1 (3m_pi)^2 | |
465 | = 2 (m_rho0)^2 | |
466 | = 3 (m_phi)^2 default: 1 | |
467 | what (3) upper limit in M^2 integration | |
468 | = 1 s/2 | |
469 | = 2 s/4 | |
470 | = 3 s default: 3 | |
471 | what (4) CKMT F_2 structure function | |
472 | = 2212 proton | |
473 | = 100 deuteron default: 2212 | |
474 | what (5) calculation of gamma-nucleon xsections | |
475 | = 1 according to CKMT-parametrization of F_2 | |
476 | = 2 integrating SIGVP over M^2 | |
477 | = 3 using SIGGA | |
478 | = 4 PHOJET cross sections default: 4 | |
479 | ||
480 | what (6), sdum no meaning | |
481 | ______________________________________________________________________ | |
482 | ||
483 | control card: codewd = VDM-PAR2 | |
484 | ||
485 | parameters in gamma-nucleus cross section calculation | |
486 | ||
487 | what (1) = 0 no suppression of shadowing by direct photon | |
488 | processes | |
489 | = 1 suppression .. default: 1 | |
490 | what (2) = 0 no suppression of shadowing by anomalous | |
491 | component if photon-F_2 | |
492 | = 1 suppression .. default: 1 | |
493 | what (3) = 0 no suppression of shadowing by coherence | |
494 | length of the photon | |
495 | = 1 suppression .. default: 1 | |
496 | what (4) = 1 longitudinal polarized photons are taken into | |
497 | account | |
498 | eps*R*Q^2/M^2 = what(4)*Q^2/M^2 default: 0 | |
499 | what (5..6), sdum no meaning | |
500 | ______________________________________________________________________ | |
501 | ||
502 | control card: codewd = XS-TABLE | |
503 | ||
504 | output of cross section table for requested interaction | |
505 | - particle production deactivated ! - | |
506 | ||
507 | what (1) lower energy limit for tabulation | |
508 | > 0 Lab. frame | |
509 | < 0 nucleon-nucleon cms | |
510 | what (2) upper energy limit for tabulation | |
511 | > 0 Lab. frame | |
512 | < 0 nucleon-nucleon cms | |
513 | what (3) > 0 # of equidistant lin. bins in E | |
514 | < 0 # of equidistant log. bins in E | |
515 | what (4) lower limit of particle virtuality (photons) | |
516 | what (5) upper limit of particle virtuality (photons) | |
517 | what (6) > 0 # of equidistant lin. bins in Q^2 | |
518 | < 0 # of equidistant log. bins in Q^2 | |
519 | ______________________________________________________________________ | |
520 | ||
521 | control card: codewd = CENTRAL | |
522 | ||
523 | Biasing in impact parameter | |
524 | ||
525 | what (1) = 1. central production | |
526 | (not recommended, has to be updated) | |
527 | what (1) < 0 and > -100 | |
528 | what (2) = min. impact parameter | |
529 | what (3) = max. impact parameter | |
530 | what (1) < -99 | |
531 | what (2) = fraction of cross section | |
532 | what (4..6), sdum no meaning | |
533 | ||
534 | Note: if what (1) = -1 : evaporation is suppressed | |
535 | if what (1) < -1 : evaporation is allowed | |
536 | ______________________________________________________________________ | |
537 | ||
538 | ||
539 | *** parameters in DPM two-chain approximation *** | |
540 | ||
541 | ______________________________________________________________________ | |
542 | ||
543 | control card: codewd = RECOMBIN | |
544 | ||
545 | Chain recombination | |
546 | (recombine S-S and V-V chains to V-S chains) | |
547 | ||
548 | what (1) = -1. recombination switched off default: 1 | |
549 | what (2..6), sdum no meaning | |
550 | ||
551 | Note: Limited applicability for MODEL = PHOJET. | |
552 | ______________________________________________________________________ | |
553 | ||
554 | control card: codewd = COMBIJET | |
555 | ||
556 | chain fusion (2 q-aq --> qq-aqaq) | |
557 | ||
558 | what (1) = 1 fusion treated default: 0. | |
559 | what (2) minimum number of uncombined chains from | |
560 | single projectile or target nucleons | |
561 | default: 0. | |
562 | what (3..6), sdum no meaning | |
563 | ||
564 | Note: Limited applicability for MODEL = PHOJET. | |
565 | ______________________________________________________________________ | |
566 | ||
567 | control card: codewd = CRONINPT | |
568 | ||
569 | Cronin effect (multiple scattering of partons at chain ends) | |
570 | ||
571 | what (1) = -1 Cronin effect not treated | |
572 | default (h+A): 1 (A+A): 0 | |
573 | what (2) scattering parameter default: 0.64 | |
574 | what (3..6), sdum no meaning | |
575 | ||
576 | Note: The Cronin-treatment should not be invoked for A+A int. | |
577 | ______________________________________________________________________ | |
578 | ||
579 | control card: codewd = DIQUARKS | |
580 | ||
581 | what (1) = -1. sea-diquark/antidiquark-pairs not treated | |
582 | default: -1. | |
583 | what (2..6), sdum no meaning | |
584 | ||
585 | Note: Limited applicability for MODEL = PHOJET. | |
586 | ______________________________________________________________________ | |
587 | ||
588 | ||
589 | *** hadronization and JETSET-parameters *** | |
590 | ||
591 | ______________________________________________________________________ | |
592 | ||
593 | control card: codewd = LUND-MSTU | |
594 | ||
595 | parameter MSTU in JETSET-common /LUDAT1/ | |
596 | ||
597 | what (1) = index according to LUND-common block | |
598 | what (2) = new value of MSTU( int(what(1)) ) | |
599 | what (3), what(4) and what (5), what(6) further | |
600 | parameter used in the same way as what (1) and what (2) | |
601 | ||
602 | Note: The use of this card is not recommended. Some parameters | |
603 | can presently not be changed with this card anyway. | |
604 | Limited applicability for MODEL = PHOJET. | |
605 | ______________________________________________________________________ | |
606 | ||
607 | control card: codewd = LUND-PARJ | |
608 | ||
609 | parameter PARJ in JETSET-common /LUDAT1/ | |
610 | ||
611 | what (1) = index according to LUND-common block | |
612 | what (2) = new value of PARJ( int(what(1)) ) | |
613 | what (3), what(4) and what (5), what(6) further | |
614 | parameter used in the same way as what (1) and what (2) | |
615 | ||
616 | Note: The use of this card is not recommended. Some parameters | |
617 | can presently not be changed with this card anyway. | |
618 | Limited applicability for MODEL = PHOJET. | |
619 | ______________________________________________________________________ | |
620 | ||
621 | control card: codewd = LUND-PARU | |
622 | ||
623 | parameter PARJ in JETSET-common /LUDAT1/ | |
624 | ||
625 | what (1) = index according to LUND-common block | |
626 | what (2) = new value of PARU( int(what(1)) ) | |
627 | what (3), what(4) and what (5), what(6) further | |
628 | parameter used in the same way as what (1) and what (2) | |
629 | ||
630 | Note: The use of this card is not recommended. Some parameters | |
631 | can presently not be changed with this card anyway. | |
632 | Limited applicability for MODEL = PHOJET. | |
633 | ______________________________________________________________________ | |
634 | ||
635 | control card: codewd = POPCORN | |
636 | ||
637 | "Popcorn-effect" in fragmentation | |
638 | ||
639 | what (1) < 0 Popcorn-effect switched off (MSTJ(12) = 1) | |
640 | >=0 Popcorn-effect treated (PARJ(5) = what (1)) | |
641 | default: 0.15 | |
642 | what (2..6), sdum no meaning | |
643 | ______________________________________________________________________ | |
644 | ||
645 | control card: codewd = PARDECAY | |
646 | ||
647 | what (1) = 1. Sigma0/Asigma0 decay treated by JETSET | |
648 | = 2. pion^0 decay after intranucl. cascade | |
649 | default: 0 (no such decays) | |
650 | what (2..6), sdum no meaning | |
651 | ______________________________________________________________________ | |
652 | ||
653 | ||
654 | *** nuclear fragmentation *** | |
655 | ||
656 | ______________________________________________________________________ | |
657 | ||
658 | control card: codewd = FERMI | |
659 | ||
660 | what (1) = -1 Fermi-motion of nucleons not treated default: 1 | |
661 | what (2) = scale factor for Fermi-momentum default: 0.68 | |
662 | what (3..6), sdum no meaning | |
663 | ______________________________________________________________________ | |
664 | ||
665 | control card: codewd = TAUFOR | |
666 | ||
667 | formation time suppressed intranuclear cascade | |
668 | ||
669 | what (1) formation time (in fm/c) default: 3.1 fm/c | |
670 | what (2) number of generations followed default: 25 | |
671 | what (3) = 1. p_t-dependent formation zone | |
672 | = 2. constant formation zone default: 1 | |
673 | what (4) modus of selection of nucleus where the | |
674 | cascade if followed first | |
675 | = 1. proj./target-nucleus with probab. 1/2 | |
676 | = 2. nucleus with highest mass | |
677 | = 3. proj. nucleus if particle is moving in pos. z | |
678 | targ. nucleus if particle is moving in neg. z | |
679 | default: 1 | |
680 | what (5..6), sdum no meaning | |
681 | ______________________________________________________________________ | |
682 | ||
683 | control card: codewd = PAULI | |
684 | ||
685 | what (1) = -1 Pauli's principle for secondary | |
686 | interactions not treated default: 1 | |
687 | what (2..6), sdum no meaning | |
688 | ______________________________________________________________________ | |
689 | ||
690 | control card: codewd = COULOMB | |
691 | ||
692 | what (1) = -1. Coulomb-energy treatment switched off default: 1 | |
693 | what (2..6), sdum no meaning | |
694 | ______________________________________________________________________ | |
695 | ||
696 | control card: codewd = EVAP | |
697 | ||
698 | evaporation module of FLUKA | |
699 | ||
700 | The following options and defaults apply only if the code is linked | |
701 | to the FLUKA-library (see README file). | |
702 | ||
703 | what (1) =< -1 ==> evaporation is switched off | |
704 | >= 1 ==> evaporation is performed | |
705 | ||
706 | what (1) = i1 + i2*10 + i3*100 + i4*10000 | |
707 | (i1, i2, i3, i4 >= 0 ) | |
708 | ||
709 | i1 is the flag for selecting the T=0 level density option used | |
710 | = 1: standard EVAP level densities with Cook pairing | |
711 | energies | |
712 | = 2: Z,N-dependent Gilbert & Cameron level densities | |
713 | (default) | |
714 | = 3: Julich A-dependent level densities | |
715 | = 4: Z,N-dependent Brancazio & Cameron level densities | |
716 | ||
717 | i2 >= 1: high energy fission activated | |
718 | (default high energy fission is activated) | |
719 | ||
720 | i3 = 0: No energy dependence for level densities | |
721 | = 1: Standard Ignyatuk (1975, 1st) energy dependence | |
722 | for level densities (default) | |
723 | = 2: Standard Ignyatuk (1975, 1st) energy dependence | |
724 | for level densities with NOT used set of parameters | |
725 | = 3: Standard Ignyatuk (1975, 1st) energy dependence | |
726 | for level densities with NOT used set of parameters | |
727 | = 4: Second Ignyatuk (1975, 2nd) energy dependence | |
728 | for level densities | |
729 | = 5: Second Ignyatuk (1975, 2nd) energy dependence | |
730 | for level densities with fit 1 Iljinov & Mebel set of | |
731 | parameters | |
732 | = 6: Second Ignyatuk (1975, 2nd) energy dependence | |
733 | for level densities with fit 2 Iljinov & Mebel set of | |
734 | parameters | |
735 | = 7: Second Ignyatuk (1975, 2nd) energy dependence | |
736 | for level densities with fit 3 Iljinov & Mebel set of | |
737 | parameters | |
738 | = 8: Second Ignyatuk (1975, 2nd) energy dependence | |
739 | for level densities with fit 4 Iljinov & Mebel set of | |
740 | parameters | |
741 | ||
742 | i4 >= 1: Original Gilbert and Cameron pairing energies used | |
743 | (default Cook's modified pairing energies) | |
744 | ||
745 | what (2) = ig + 10 * if (ig and if must have the same sign) | |
746 | ||
747 | ig =< -1 ==> deexcitation gammas are not produced | |
748 | (if the evaporation step is not performed | |
749 | they are never produced) | |
750 | if =< -1 ==> Fermi Break Up is not invoked | |
751 | (if the evaporation step is not performed | |
752 | it is never invoked) | |
753 | The default is: deexcitation gamma are produced and | |
754 | Fermi break up is activated | |
755 | what (3..6), sdum no meaning | |
756 | ______________________________________________________________________ | |
757 | ||
758 | ||
759 | *** output and checks *** | |
760 | ||
761 | ______________________________________________________________________ | |
762 | ||
763 | control card: codewd = FRAME | |
764 | ||
765 | frame in which final state is given in DTEVT1 | |
766 | ||
767 | what (1) = 1 target rest frame (laboratory) | |
768 | = 2 nucleon-nucleon cms default: 1 | |
769 | ______________________________________________________________________ | |
770 | ||
771 | ||
772 | control card: codewd = HISTOGRAM | |
773 | ||
774 | activate different classes of histograms | |
775 | ||
776 | default: no histograms | |
777 | ______________________________________________________________________ | |
778 | ||
779 | control card: codewd = EMCCHECK | |
780 | ||
781 | extended energy-momentum / quantum-number conservation check | |
782 | ||
783 | what (1) = -1 extended check not performed default: -1. | |
784 | what (2..6), sdum no meaning | |
785 | ______________________________________________________________________ | |
786 | ||
787 | ||
788 | *** lepton tagger *** | |
789 | ||
790 | ______________________________________________________________________ | |
791 | ||
792 | control card: codewd = L-TAG | |
793 | (lepton-nucleus interactions with MODEL=PHOJET only) | |
794 | ||
795 | lepton tagger: | |
796 | definition of kinematic cuts for radiated photon and | |
797 | outgoing lepton detection in lepton-nucleus interactions | |
798 | ||
799 | what (1) = y_min | |
800 | what (2) = y_max | |
801 | what (3) = Q^2_min | |
802 | what (4) = Q^2_max | |
803 | what (5) = theta_min (Lab) | |
804 | what (6) = theta_max (Lab) default: no cuts | |
805 | sdum no meaning | |
806 | ______________________________________________________________________ | |
807 | ||
808 | ||
809 | control card: codewd = L-ETAG | |
810 | (lepton-nucleus interactions with MODEL=PHOJET only) | |
811 | ||
812 | lepton tagger: | |
813 | what (1) = min. outgoing lepton energy (in Lab) | |
814 | what (2) = min. photon energy (in Lab) | |
815 | what (3) = max. photon energy (in Lab) default: no cuts | |
816 | what (2..6), sdum no meaning | |
817 | ______________________________________________________________________ | |
818 | ||
819 | control card: codewd = ECMS-CUT | |
820 | (lepton-nucleus interactions with MODEL=PHOJET only) | |
821 | ||
822 | what (1) = min. c.m. energy to be sampled | |
823 | what (2) = max. c.m. energy to be sampled | |
824 | what (3) = min x_Bj to be sampled default: no cuts | |
825 | what (3..6), sdum no meaning | |
826 | ______________________________________________________________________ | |
827 | ||
828 | ||
829 | *** for code development only *** | |
830 | ||
831 | ______________________________________________________________________ | |
832 | ||
833 | control card: codewd = INTPT | |
834 | ||
835 | what (1) = -1 intrinsic transverse momenta of partons | |
836 | not treated default: 1 | |
837 | what (2..6), sdum no meaning | |
838 | ______________________________________________________________________ | |
839 | ||
840 | control card: codewd = OUTLEVEL | |
841 | ||
842 | output control switches | |
843 | ||
844 | what (1) = internal rejection informations default: 0 | |
845 | what (2) = energy-momentum conservation check output default: 0 | |
846 | what (3..6) internal warning messages default: 0 | |
847 | ______________________________________________________________________ | |
848 | ||
849 | control card: codewd = RESONANC | |
850 | ||
851 | treatment of low mass chains | |
852 | ||
853 | what (1) = -1 low chain masses are not corrected for resonance | |
854 | masses default: 1. | |
855 | what (2) = -1 massless partons default: 1. (massive) | |
856 | what (3) = -1 chain-system containing chain of too small | |
857 | mass is rejected (note: this does not fully | |
858 | apply to S-S chains) default: 0. | |
859 | what (4..6), sdum no meaning | |
860 | ______________________________________________________________________ | |
861 | ||
862 | control card: codewd = SEASU3 | |
863 | ||
864 | Treatment of strange-quarks at chain ends | |
865 | ||
866 | what (1) (SEASQ) strange-quark suppression factor | |
867 | iflav = 1.+rndm*(2.+SEASQ) default: 1. | |
868 | what (2..6), sdum no meaning | |
869 | ______________________________________________________________________ | |
870 | ||
871 | control card: codewd = XCUTS | |
872 | ||
873 | thresholds for x-sampling | |
874 | ||
875 | what (1) defines lower threshold for val.-q x-value (CVQ) | |
876 | default: 1. | |
877 | what (2) defines lower threshold for val.-qq x-value (CDQ) | |
878 | default: 2. | |
879 | what (3) defines lower threshold for sea-q x-value (CSEA) | |
880 | default: 0.2 | |
881 | what (4) sea-q x-values in S-S chains (SSMIMA) default: 0.14 | |
882 | what (5) not used default: 2. | |
883 | what (6), sdum no meaning | |
884 | ||
885 | Note: Lower thresholds (what(1..3)) are def. as x_thr=CXXX/ECM | |
886 | ______________________________________________________________________ | |
887 | ||
888 | ||
889 | ||
890 | EVENT HISTORY - COMMON /DTEVT1/ | |
891 | _________________________________ | |
892 | ||
893 | ||
894 | 1) Common block DTEVT1 | |
895 | ---------------------- | |
896 | ||
897 | The complete history of each event can be found in COMMON /DTEVT1/ | |
898 | and COMMON /EXTEVT/. | |
899 | ||
900 | COMMON /DTEVT1/ : | |
901 | NHKK number of entries in common block | |
902 | NEVHKK number of the event | |
903 | ISTHKK(i) status code for entry i | |
904 | IDHKK(i) identifier for the entry | |
905 | (for particles: identifier according to the PDG scheme) | |
906 | JMOHKK(1,i) pointer to the entry of the first mother of entry i | |
907 | JMOHKK(2,i) pointer to the entry of the second mother of entry i | |
908 | JDAHKK(1,i) pointer to the entry of the first daughter of entry i | |
909 | JDAHKK(2,i) pointer to the entry of the second daughter of entry i | |
910 | PHKK(1..3,i) 3-momentum | |
911 | PHKK(4,i) energy | |
912 | PHKK(5,i) mass | |
913 | VHKK / WHKK spatial position of particle in target / projectile | |
914 | rest frame | |
915 | ||
916 | 2) Final state particles | |
917 | ------------------------ | |
918 | ||
919 | The final state particles from the actual event (number NEVHKK) | |
920 | can be found in DTEVT1 and identified by their status: | |
921 | ||
922 | ISTHKK(i) = 1 final state particle produced in | |
923 | photon-/hadron-/nucleon-nucleon collisions or | |
924 | in intranuclear cascade processes | |
925 | -1 nucleons, deuterons, H-3, He-3, He-4 evaporated | |
926 | from excited nucleus, fragmentation and fission | |
927 | products (A > 4) and photons produced in nuclear | |
928 | deexcitation processes | |
929 | 1001 residual nucleus (ground state) | |
930 | ||
931 | The types of these particles/nuclei are given in IDHKK as follows | |
932 | ||
933 | all final state particles except nuclei : | |
934 | IDHKK(i)=particle identifier according to PDG numbering scheme | |
935 | nuclei (A > 1: evaporation / fragmentation / fission products, and | |
936 | residual nuclei) : | |
937 | IDHKK(i)=80000, IDRES(i)=mass number, IDXRES(i)=charge number | |
938 | ||
939 | The 4-momenta and masses can be found in PHKK (target nucleus rest frame | |
940 | unless defined by the FRAME-card): | |
941 | PHKK(1..3,i) 3-momentum (p_x,p_y,p_z) | |
942 | PHKK(4,i) energy | |
943 | PHKK(5,i) mass |