!  GeIsDalitzDecay.f90
!
!  Returns true if argument is a 3-body decay with two 
!  electrons or muons, false otherwise.
!
!  Detlef Irmscher, Thomas Ullrich, Uni Heidelberg, September 1993
!  Last update: 24.3.95 tu  - added muons
!
logical function GeIsDalitzDecay (thisDecay)
    use genesis_m
    implicit none
    type(decay)             :: thisDecay
    integer                 :: ibody, nleptons
    
    GeIsDalitzDecay = .false.

    ! check if 3 body decay
    if (thisDecay%nbody /= 3)   return

    ! check if 2 leptons are present
    nleptons = 0
    do ibody = 1, 3
        if (abs(thisDecay%child(ibody)%id) == 12 .or. &
	    abs(thisDecay%child(ibody)%id) == 14 )  nleptons = nleptons + 1
    enddo

    if (nleptons == 2)  GeIsDalitzDecay = .true.   

end function GeIsDalitzDecay