DPMJET/flukapro/(HADFLG)

   1 *$ CREATE HADFLG.ADD
   2 *COPY HADFLG
   3 *
   4 *=== hadflg ===========================================================*
   5 *
   6 *
   7 *----------------------------------------------------------------------*
   8 *                                                                      *
   9 *     Created on 13 september 1991 by    Alfredo Ferrari & Paola Sala  *
  10 *                                                   Infn - Milan       *
  11 *                                                                      *
  12 *     Last change on 08-sep-00     by    Alfredo Ferrari               *
  13 *                                                                      *
  14 *     Included in the following routines:                              *
  15 *                                                                      *
  16 *                                         Blkdt4                       *
  17 *                                         Eventv                       *
  18 *                                         Ferhav                       *
  19 *                                         Hadriv                       *
  20 *                                         Hdncin                       *
  21 *                                         Hhflev                       *
  22 *                                         Hncnev                       *
  23 *                                         Nucdcy                       *
  24 *                                         Nuncln                       *
  25 *                                         Phycrd                       *
  26 *                                         Rchanv                       *
  27 *                                         Twopar                       *
  28 *                                                                      *
  29 *                    Lnwhdr = flag for using the (still on development)*
  30 *                             new Hadriv                               *
  31 *                    Ielflg = Flag to control the elastic channel      *
  32 *                             in Hadriv:                               *
  33 *                             -1: elastic channel reduced according    *
  34 *                                 to a rough Pauli blocking factor     *
  35 *                              0: elastic channel fully accounted for  *
  36 *                             +1: elastic channel suppressed           *
  37 *                    Icxflg = Flag to control the charge exchange      *
  38 *                             channel in Hadriv:                       *
  39 *                             -1: ch. ex. channel reduced according    *
  40 *                                 to a rough Pauli blocking factor     *
  41 *                              0: ch. ex. channel fully accounted for  *
  42 *                             +1: ch. ex. channel suppressed           *
  43 *                    Istflg = Flag to control the strangeness exchange *
  44 *                             channel in Hadriv:                       *
  45 *                             -1: st. ex. channel reduced according    *
  46 *                                 to a rough Pauli blocking factor     *
  47 *                              0: st. ex. channel fully accounted for  *
  48 *                             +1: st. ex. channel suppressed           *
  49 *                    Ianflg = Flag to control the annhiliation channel *
  50 *                             in Hadriv:                               *
  51 *                             -1: annih.  channel only one accounted   *
  52 *                                 for (!!! note please !!!)            *
  53 *                              0: annih.  channel fully accounted for  *
  54 *                             +1: annih.  channel suppressed           *
  55 *             Innure(1,1,k) = Projectile index for the Hadriv reaction *
  56 *                             k (first entrance channel)               *
  57 *             Innure(2,1,k) = Target index for the Hadriv reaction     *
  58 *                             k (first entrance channel)               *
  59 *             Innure(1,2,k) = Projectile index for the Hadriv reaction *
  60 *                             k (second entrance channel, if any)      *
  61 *             Innure(2,2,k) = Target index for the Hadriv reaction     *
  62 *                             k (second entrance channel, if any)      *
  63 *                             It holds:                                *
  64 *                        Nure (Innure(1,j,k),(Innure(2,j,k)/8+1)) = k  *
  65 *                  Areso1/2 = mass   of the 1st/2nd resonance produced *
  66 *                             in Hadriv                                *
  67 *                  Greso1/2 = width  of the 1st/2nd resonance produced *
  68 *                             in Hadriv                                *
  69 *                  Ereso1/2 = energy of the 1st/2nd resonance produced *
  70 *                             in Hadriv                                *
  71 *                  Preso1/2 = momen. of the 1st/2nd resonance produced *
  72 *                             in Hadriv                                *
  73 *              Px,y,zres1/2 = mom. comp. of the 1st/2nd resonance pro- *
  74 *                             duced in Hadrin                          *
  75 *                  Jreso1/2 = index of the 1st/2nd resonance produced  *
  76 *                             in Hadriv                                *
  77 *               Kreso1/2(m) = pointer (inside ...j arrays) of 1st decay*
  78 *                             particles originating from the 1st/2nd   *
  79 *                             resonance produced in Hadriv             *
  80 *                             Kreso.. = k1 + 1000 * k2 + 1000000 * k3  *
  81 *                  Nreso1/2 = number of final particles originating    *
  82 *                             from the 1st/2nd resonance produced in   *
  83 *                             Hadriv                                   *
  84 *                 Aresoj(i) = mass   of the i_th intermediate resonance*
  85 *                             produced in Hadriv                       *
  86 *                 Gresoj(i) = width  of the i_th intermediate resonance*
  87 *                             produced in Hadriv                       *
  88 *                 Eresoj(i) = energy of the i_th intermediate resonance*
  89 *                             produced in Hadriv                       *
  90 *                 Presoj(i) = momen. of the i_th intermediate resonance*
  91 *                             produced in Hadriv                       *
  92 *             Px,y,zresj(i) = mom. comp. of the i_th intermediate reso-*
  93 *                             nance produced in Hadrin                 *
  94 *                 Jresoj(i) = index of the i_th intermediate resonance *
  95 *                             produced in Hadriv                       *
  96 *                 Kresoj(i) = pointer (inside ...j arrays) of the decay*
  97 *                             particles originating from the i_th int- *
  98 *                             ermediate resonance produced in Hadriv   *
  99 *                             Kresoj = k1 + 1000 * k2 + 1000000 * k3   *
 100 *                    Qqreso = Q^2=-t of the resonance creation process *
 101 *                             in Hadriv (projectile and resonance 1)   *
 102 *                    Uureso = u of the resonance creation process in   *
 103 *                             Hadriv (target and resonance 1)          *
 104 *                    Rhresp = local proton  nuclear density (fm^-3)    *
 105 *                    Rhresn = local neutron nuclear density (fm^-3)    *
 106 *                                                                      *
 107 *----------------------------------------------------------------------*
 108 *
 109       PARAMETER ( MXRSDC = 25 )
 110       PARAMETER ( PPAMXB = 0.6       D+00 )
 111       PARAMETER ( PAUMXB = 2.D+00 * 0.8 D+00 )
 112       PARAMETER ( PPAMXM = 0.6       D+00 )
 113       PARAMETER ( PAUMXM = 2.D+00 * 0.4 D+00 )
 114 *
 115       COMMON / HADFLG / QQRESO, UURESO, RHRESP, RHRESN,
 116      &               ARESOJ (MXRSDC), GRESOJ (MXRSDC), ERESOJ (MXRSDC),
 117      &               PRESOJ (MXRSDC), PXRESJ (MXRSDC), PYRESJ (MXRSDC),
 118      &               PZRESJ (MXRSDC), JRESOJ (MXRSDC), KDCRSJ (MXRSDC),
 119      &                  IELFLG, ICXFLG, ISTFLG, IANFLG, IOLDHD, LNWHDR,
 120      &                  INNURE (2,2,16), IKCHXG (16),
 121      &                  NRESO1, NRESO2
 122       LOGICAL LNWHDR
 123       EQUIVALENCE ( ARESO1, ARESOJ (1) ), ( ARESO2, ARESOJ (2) )
 124       EQUIVALENCE ( GRESO1, GRESOJ (1) ), ( GRESO2, GRESOJ (2) )
 125       EQUIVALENCE ( ERESO1, ERESOJ (1) ), ( ERESO2, ERESOJ (2) )
 126       EQUIVALENCE ( PRESO1, PRESOJ (1) ), ( PRESO2, PRESOJ (2) )
 127       EQUIVALENCE ( PXRES1, PXRESJ (1) ), ( PXRES2, PXRESJ (2) )
 128       EQUIVALENCE ( PYRES1, PYRESJ (1) ), ( PYRES2, PYRESJ (2) )
 129       EQUIVALENCE ( PZRES1, PZRESJ (1) ), ( PZRES2, PZRESJ (2) )
 130       EQUIVALENCE ( KRESO1, KDCRSJ (1) ), ( KRESO2, KDCRSJ (2) )
 131       EQUIVALENCE ( JRESO1, JRESOJ (1) ), ( JRESO2, JRESOJ (2) )
 132