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