- AliPhysicsSelection: protected writing of fHistBunchCrossing and fHistTriggerPattern
[u/mrichter/AliRoot.git] / DPMJET / doc / dpmjet3.manual
CommitLineData
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