yyjet.f

Go to the documentation of this file.

CDECK  ID>, YYJET.
      SUBROUTINE yyjet(NJET,YL,YH,IERR)
C
C  ROUTINE TO RETURN THE VALUES OF YCUT BETWEEN WHICH EVENT IS
C  CLASSIFIED AS N-JET (YL < YH)
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.
C
      IMPLICIT NONE
C     IMPLICIT NONE
      INTEGER nyclmx,njetmx,njet,ierr,imodeo,njeto,nto
      parameter(nyclmx = 500, njetmx = 20)
      INTEGER histor(2,nyclmx)
      REAL yrec(njetmx),yl,yh,pint(10,nyclmx),pjet(10,njetmx,njetmx)
      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(' #### YYJET: YKERN MUST BE CALLED FIRST ! ####')
        yl = -1.
        yh = -1.
        RETURN
      ENDIF
C
CHECK IF INPUT MAKES SENSE
C
      IF(njet.LE.0 .OR. njet.GT.njetmx) THEN
        WRITE(6,1) njet
 1      FORMAT(' #### YYJET: REQUEST FOR NJET=',i12,
     +  ' NOT SUPPORTED ####')
        yl = -1.
        yh = -1.
        RETURN
      ENDIF
      IF(njet.GT.nto) THEN
        WRITE(6,2) njet,nto
 2      FORMAT(' #### YYJET:',i3,' JETS OUT OF',i3,' PARTICLES NOT',
     +  ' POSSIBLE. ####')
        yl = -1.
        yh = -1.
        RETURN
      ENDIF
C
      IF(njet.EQ.1) THEN
        yh = 1.
      ELSE
        yh = yrec(njet-1)
      ENDIF
      yl = yrec(njet)
C
      ierr = 0
      RETURN
      END