!  GePartProd.f90
!
!  Calculates the number of each sort of parent particles.
!
!  Detlef Irmscher, Thomas Ullrich, Uni Heidelberg, September 1993
!  Last update: 18.3.95 tu
!
integer function GePartProd (ipart, nch_over_central_over_total_pi0)
    use genesis_m, self=>GePartProd
    implicit none
    integer, intent(in)    :: ipart
    integer, intent(in)    :: nch_over_central_over_total_pi0
    real(kind=double)      :: xpart

    if (particleList(ipart)%relSigmaTotal == 0) then       ! avoid trouble
       GePartProd = 0
       return
    endif

    xpart = pion0_over_nch * nch_over_central_over_total_pi0
    if (is_mass_weight == 0) xpart = xpart*particleList(ipart)%relSigmaTotal
    GePartProd = GeNpois(xpart)

end function GePartProd