ynjet.f

Go to the documentation of this file.

CDECK  ID>, YNJET.
      SUBROUTINE ynjet(YCUT,NJET,IERR)
C
C  ROUTINE TO RETURN THE NUMBER OF JETS FOR A GIVEN YCUT
C
C  LAST MOD  :  21-Jun-98
C
C  Modification Log.
C  08-Oct-97 D. Chrisman, Increase NYCLMX from 250 to 500.
C  21-Jun-98 D. Chrisman, Want to be able to handle more than 10 jets.
C                         Introduce parameter NJETMX. Use this
C                         in declaration of PJET and YREC in order to remove
C                         hard coded numbers. Also replace the hard
C                         coded number in the "D0 5001" loop with the
C                         parameter NJETMX.
C
      IMPLICIT NONE
C     IMPLICIT NONE
      INTEGER njet,ierr,nyclmx,njetmx,imodeo,njeto,nto,i
      parameter(nyclmx = 500, njetmx = 20)
      REAL ycut,yrec(njetmx),pint(10,nyclmx),pjet(10,njetmx,njetmx)
      INTEGER histor(2,nyclmx)
      COMMON /ycl/yrec,pjet,histor
      COMMON /yint/ imodeo,nto,njeto,pint
      ierr = -1
C
CHECK IF CALL WAS MADE TO YKERN
C
      IF(imodeo.LE.0 .OR. imodeo.GT.7) THEN
        WRITE(6,111)
 111    FORMAT(' #### YNJET: YKERN MUST BE CALLED FIRST ! ####')
        njet = -1
        RETURN
      ENDIF
C
CHECK IF INPUT MAKES SENSE
C
      IF(ycut.LE.0. .OR. ycut.GT.1.) THEN
        WRITE(6,1) ycut
 1      FORMAT(' #### YNJET: INPUT YCUT=',e12.4,' INVALID ####')
        njet = -1
        RETURN
      ENDIF
C
      njet = 1
      DO 5001 i=1,njetmx
        IF(ycut.LT.yrec(i)) njet = i+1
 5001 CONTINUE
      ierr = 0
      RETURN
      END