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