*
*=== sgtbcm ===========================================================*
*
*----------------------------------------------------------------------*
*                                                                      *
*     Copyright (C) 2003-2019:  CERN & INFN                            *
*     All Rights Reserved.                                             *
*                                                                      *
*     Include file: sgtbcm                                             *
*                                                                      *
*     Created  on  26  August  1991   by        Alfredo Ferrari        *
*                                                INFN - Milan          *
*                                                                      *
*     Included in the following routines:                              *
*                                                                      *
*                BDTRNS                                                *
*                CKSIGI                                                *
*                FIPROU                                                *
*                FLLEMF                                                *
*                HVBREM                                                *
*                HVPAIR                                                *
*                KASHEA                                                *
*                KASHOP                                                *
*                KASKAD                                                *
*                PHNCVR                                                *
*                PHNVEV                                                *
*                SGTBRE                                                *
*                SGTELS                                                *
*                SGTEMD                                                *
*                SGTINL                                                *
*                SGTPAI                                                *
*                SGTTOT                                                *
*                SIGTAB                                                *
*                ZEROIN                                                *
*                                                                      *
*  !!!! Please note that heavies (d,t,3-He,alphas,HI) must have  !!!!  *
*  !!!! the last indeces in the Ijsig... arrays (i.e. Mpsige-5,  !!!!  *
*  !!!! Mpsige-4, Mpsige-3, Mpsige-2, Mpsige)                    !!!!  *
*                                                                      *
*       Ktsigt(i,m) = Initial address (0 index) for total cross        *
*                     section tabulations for the i_th particle type   *
*                     as defined by Ijsigt and the m_th material       *
*                     Ktsigt does not include EM dissociation if       *
*                     activated                                        *
*       Kisigt(i,m) = Initial address (0 index) for inelastic cross    *
*                     section tabulations for the i_th particle type   *
*                     as defined by Ijsigi and the m_th material       *
*       Kesigt(i,m) = Initial address (0 index) for elastic cross      *
*                     section tabulations for the i_th particle type   *
*                     as defined by Ijsige and the m_th material       *
*       Kbsigt(i,m) = Initial address (0 index) for bremsstr. cross    *
*                     section tabulations for the i_th particle type   *
*                     as defined by Ijsigb and the m_th material       *
*       Kpsigt(i,m) = Initial address (0 index) for pair prod. cross   *
*                     section tabulations for the i_th particle type   *
*                     as defined by Ijsigb and the m_th material       *
*       Kdsigt(i,m) = Initial address (0 index) for delta prod. cross  *
*                     section tabulations for the i_th particle type   *
*                     as defined by Ijsigb and the m_th material       *
*       Kmsigt(i,m) = Initial address (0 index) for EM diss. cross     *
*                     section tabulations for the i_th particle type   *
*                     as defined by Ijsigm and the m_th material       *
*       Kcsigt(i,m) = Initial address (0 index) for Coulomb cross      *
*                     section tabulations for the i_th particle type   *
*                     as defined by Ijsigb and the m_th material       *
*       Lblsgi(i,m) = flag for computing the inelastic cross section   *
*                     when below Eksig0 for the i_th particle type     *
*                     as defined by Ijsigi and the m_th material       *
*                                                                      *
*----------------------------------------------------------------------*
*
*  Minimum energy ratio for which cross section are computed
*  (that is the minimum energy down which the cross sections presently
*   used are reliable)
      PARAMETER ( EKRTMN = 0.6000000000000000000D+00 )
*  Sqr1p5 = sqrt ( 1 / ekrtmn )
      PARAMETER ( SQR1P5 = 1.2909944487358056284D+00 )
*  Alg1p5 = log  ( 1 / ekrtmn )
      PARAMETER ( ALG1P5 = 0.5108256237659906832D+00 )
*  Minimum energy ratio before computing again cross sections
*  (roughly geometric mean between 1 and ekrtmn )
      PARAMETER ( EKRTEF = 0.7750000000000000000D+00 )
      PARAMETER ( EKSIG0 = 1.104746407608999714039030918809046D-03 )
*     PARAMETER ( EKSIG0 = 6.14229867868362340504446476139873D-03 )
      PARAMETER ( EKSIG1 = 50.0000000000000000000000000000000D-03 )
      PARAMETER ( EKSIG2 = 641.959232274432000000000000000000D-03 )
      PARAMETER ( EKSIG3 = 11677.1704101562500000000000000000D-03 )
      PARAMETER ( ELKEG0 =-6.80813946572698539519330850737315D+00 )
*     PARAMETER ( ELKEG0 =-5.09255622924913791440217028831937D+00 )
      PARAMETER ( ELKEG1 =-2.99573227355399099343522357614254D+00 )
      PARAMETER ( ELKEG2 =-4.43230478438626226471171223979336D-01 )
      PARAMETER ( ELKEG3 = 2.45763568864610059849066495004851D+00 )
      PARAMETER ( SIGRT0 = 1.1000000000000000000D+00 )
      PARAMETER ( SIGRT1 = 1.2000000000000000000D+00 )
      PARAMETER ( SIGRT2 = 1.2500000000000000000D+00 )
      PARAMETER ( SIGRT3 = 1.4000000000000000000D+00 )
      PARAMETER ( ASGRT0 = 9.531017980432486004395212328076517D-02 )
      PARAMETER ( ASGRT1 = 0.182321556793954626211718025154515D+00 )
      PARAMETER ( ASGRT2 = 0.223143551314209755766295090309834D+00 )
      PARAMETER ( ASGRT3 = 0.336472236621212930504593410216992D+00 )
      PARAMETER ( NPSIG0 = 40 )
*     PARAMETER ( NPSIG0 = 22 )
      PARAMETER ( NPSIG1 = 14 )
      PARAMETER ( NPSIG2 = 13 )
      PARAMETER ( MPSIGT = 36 )
      PARAMETER ( MPSIGE = 35 )
      PARAMETER ( MPSIGI = 39 )
      PARAMETER ( MPSIGB = 14 )
      PARAMETER ( MPSIGM = 15 )
      PARAMETER ( MPSIGH =  5 )
*
      LOGICAL LKSHRT, LK0BAR, LKAON0, LKAOSL, LEXDCY, LBLSGI
      REAL * 4 SGDLIX, SGDLJX, SGDLME, SGDLKE, SGDLIE, SGDLJE
      COMMON / SGTBCM / EKSIG4, ELKE0, ELKE1, ESSGTB, DESGTB, EKIOLD,
     &                  EKSDTH (MPSIGB,MXXMDF), EKSDMX (MPSIGB,MXXMDF),
     &                  SGDLMX (MPSIGB,MXXMDF), SGDLIX (MPSIGB,MXXMDF),
     &                  SGDLJX (MPSIGB,MXXMDF), SGDLKE (MPSIGB,MXXMDF),
     &                  SGDLIE (MPSIGB,MXXMDF), SGDLJE (MPSIGB,MXXMDF),
     &                  SGDLME (MPSIGB,MXXMDF),
     &                  KDLTIE (MPSIGB,MXXMDF), KDLTIX (MPSIGB,MXXMDF),
     &                  KTSIGT (MPSIGT,MXXMDF), KESIGT (MPSIGE,MXXMDF),
     &                  KISIGT (MPSIGI,MXXMDF), KBSIGT (MPSIGB,MXXMDF),
     &                  KPSIGT (MPSIGB,MXXMDF), KDSIGT (MPSIGB,MXXMDF),
     &                  KMSIGT (MPSIGM,MXXMDF), KCSIGT (MPSIGB,MXXMDF),
     &                  LBLSGI (MPSIGI,MXXMDF),
     &                  IJSIGT (-6:NALLWP), IJSIGE (-6:NALLWP),
     &                  IJSIGI (-6:NALLWP), IJSIGB (-6:NALLWP),
     &                  IJSIGM (-6:NALLWP),
     &                  NPSIG3, IESGTB, NPSGMN, NPSGMX, LKSHRT, LK0BAR,
     &                  LKAON0, LKAOSL, LEXDCY
      SAVE / SGTBCM /