*
*=== Cphspc ===========================================================*
*
*----------------------------------------------------------------------*
*                                                                      *
*     Copyright (C) 2003-2019:  CERN & INFN                            *
*     All Rights Reserved.                                             *
*                                                                      *
*     Chain PHase-SPace break-up Common:                               *
*                                                                      *
*     Created on  10 January 2011  by    Alfredo Ferrari & Paola Sala  *
*                                           CERN - EN     Infn - Milan *
*                                                                      *
*          Sdphsp (i) = (global) spin, degeneracy, weight, and mass    *
*                       factor of the i_th break-up                    *
*          Afphsp (i) = (minimum) mass of the i_th break-up            *
*          Bfphsp (i) = mass factor in the phase space expression      *
*        Amphsp (j,k) = mass(es) (GeV) of the mass spectrum of the     *
*                       k_th particle                                  *
*        Wephsp (j,k) = weight of the j_th component of the mass       *
*                       spectrum of the k_th particle                  *
*          Gaphsp (k) = width (GeV) of the k_th particle               *
*          Kpphsp (k) = (Part) id of the k_th particle                 *
*          Kcphsp (k) = (Part) id of the k_th particle charge conjugate*
*          Icphsp (k) = electric charge of the k_th particle           *
*          Ibphsp (k) = baryonic charge of the k_th particle           *
*          Jsphsp (k) = spin (in hbar/2 units) of the k_th particle    *
*          Ipphsp (k) = parity of the k_th particle                    *
*          Lmphsp (k) = number of masses in the mass spectrum of the   *
*                       k_th particle                                  *
*       Iqphsp(m,n,k) = quark composition (m=1,2,3) of the n_th quark  *
*                       conbination of the k_th particle               *
*        Ifphsp (1,i) = net u quark content of the i_th break up comb- *
*                       ination                                        *
*        Ifphsp (2,i) = net d quark content of the i_th break up comb- *
*                       ination                                        *
*        Ifphsp (3,i) = net s quark content of the i_th break up comb- *
*                       ination                                        *
*        Ifphsp (4,i) = net c quark content of the i_th break up comb- *
*                       ination                                        *
*        Ifphsp (5,i) = first  particle to be emitted                  *
*        Ifphsp (6,i) = second particle to be emitted (if negative it  *
*                       is a combination of particles of index=| |)    *
*        Ifphsp (7,i) = (global) multiplicity of the i_th break-up     *
*        Ifphsp (8,i) = number of resonances  of the i_th break-up     *
*      Idphsp (l,m,n) = id of the break-up set corresponding to quarks *
*                       l,m,n                                          *
*          K0phsp (l) = starting break-up index for the l_th break-up  *
*                       set (eg l=1 <-> uuu)                           *
*          K1phsp (l) = last     break-up index for the l_th break-up  *
*                       set (eg l=1 <-> uuu)                           *
*          Mcphsp (k) = combination id corresponding to the k_th       *
*                       break-up index. Given a break-up set l         *
*                       (eg l=1 <-> uuu) all possible k's are those    *
*                       in the range [K0phsp(l),K1phsp(l)]. Then the   *
*                       corresponding break-up combinations i are given*
*                       by i = Mcphsp(k)                               *
*                                                                      *
*----------------------------------------------------------------------*
*  1 proton rms (fm)
      PARAMETER ( R0PHSP = 0.8 D+00 )
*  Maximum number of particles  emitted:
      PARAMETER ( MXPHSP =  6 )
*  Maximum number of resonances emitted:
      PARAMETER ( MXRPHS =  2 )
*  Maximum number of combinations:
      PARAMETER ( MXFPHS = 80000 )
*  Total number of hadrons considered
      PARAMETER ( NPPHSP = 35 )
*  Maximum number of mass intervals for resonances
      PARAMETER ( NMPHSP =  6 )
*  Final index of unflavoured baryons:
      PARAMETER ( MBPHSP =  8 )
*  Final index of unflavoured mesons:
      PARAMETER ( MMPHSP = 17 )
*  Final index of strange (+1) baryons:
      PARAMETER ( MBSPHS = 24 )
*  Final index of strange (+2) baryons:
      PARAMETER ( M2SPHS = 26 )
*  Final index of strange (+3) baryons:
      PARAMETER ( M3SPHS = 27 )
*  Number of different chains supported
      PARAMETER ( MXCPHS = 14 )
*  Invariant mass cut-off for baryon (unflavoured) chains:
      PARAMETER ( CTBPHS = 2.3 D+00 )
*  Mass cut-off smearing parameter for baryon (unflavoured) chains:
      PARAMETER ( SMBPHS = 0.3 D+00 )
*  Invariant mass cut-off for meson (unflavoured) chains:
      PARAMETER ( CTMPHS = 1.3 D+00 )
*  Mass cut-off smearing parameter for meson (unflavoured) chains:
      PARAMETER ( SMMPHS = 0.3 D+00 )
*  Mass cut-off smearing parameter for hyperon chains:
      PARAMETER ( SMHPHS = 0.3 D+00 )
*  Invariant mass cut-off for hyperon chains:
      PARAMETER ( CTHPHS = 2.6 D+00 )
*  Invariant mass cut-off for meson strangeness +/-1 chains:
      PARAMETER ( CTKPHS = 1.6 D+00 )
*  Mass cut-off smearing parameter meson strangeness +/-1 chains:
      PARAMETER ( SMKPHS = 0.3 D+00 )
*
      COMMON / CPHSPS /                          AMPHSP (NMPHSP,NPPHSP),
     &                WEPHSP (NMPHSP,NPPHSP),   DMPHSP, GAPHSP (NPPHSP),
     &                SDPHSP (MXFPHS), AFPHSP (MXFPHS), BFPHSP (MXFPHS),
     &                IQPHSP (3,2,NPPHSP), IFPHSP (8,MXFPHS),
     &                LMPHSP (NPPHSP), IBPHSP (NPPHSP), ICPHSP (NPPHSP),
     &                JSPHSP (NPPHSP), IPPHSP (NPPHSP), KPPHSP (NPPHSP),
     &                KCPHSP (NPPHSP), K0PHSP (MXCPHS), K1PHSP (MXCPHS),
     &                MCPHSP (MXFPHS), IDPHSP (-4:4,-4:4,-4:4,-4:4),
     &                IQCPHS (4,MXCPHS)