ISAJET/code/logerr.F

   1 #include "isajet/pilot.h"
   2       SUBROUTINE LOGERR(IMSG,I,IERR)
   3 C
   4 C          ERROR MESSAGES
   5 C
   6 #include "isajet/itapes.inc"
   7 #include "isajet/jetlim.inc"
   8 #include "isajet/primar.inc"
   9 #include "isajet/jetpar.inc"
  10 #include "isajet/const.inc"
  11 #include "isajet/dylim.inc"
  12 #include "isajet/keys.inc"
  13 #include "isajet/q1q2.inc"
  14
  15 C
  16 C        ERRORS IN JET PARAMETERS
  17 C
  18       IERR=IERR+1
  19       IF(IMSG.EQ.0) WRITE(ITLIS,81)
  20 81    FORMAT(//5X,'DEFAULT LIMITS HAVE BEEN SET')
  21       IF(IMSG.EQ.1) WRITE(ITLIS,1001) I,PMIN(I),PMAX(I)
  22  1001 FORMAT(//10X,'BAD LIMITS FOR P(',I2,')=',2E12.4)
  23       IF(IMSG.EQ.2) WRITE(ITLIS,1002) I,PTMIN(I),PTMAX(I)
  24  1002 FORMAT(//10X,'BAD LIMITS FOR PT(',I2,')=',2E12.4)
  25       IF(IMSG.EQ.3) WRITE(ITLIS,1003) I,THMIN(I),THMAX(I)
  26  1003 FORMAT(//10X,'BAD LIMITS FOR THETA(',I2,')=',2E12.4)
  27       IF(IMSG.EQ.4) WRITE(ITLIS,1004) I,XJMIN(I),XJMAX(I)
  28  1004 FORMAT(//10X,'BAD LIMITS FOR X(',I2,')=',2E12.4)
  29       IF(IMSG.EQ.5) WRITE(ITLIS,1005) I,XJ(I),P(I)
  30  1005 FORMAT(//5X,'X AND P FOR JET',I2,' ARE INCOMPATIBLE',2E12.4)
  31       IF(IMSG.EQ.6) WRITE(ITLIS,1006) I,THMIN(I),THMAX(I)
  32  1006 FORMAT(//10X,'LIMITS FOR THETA MUST BE .GT.0 AND .LT.PI. PRESENT'
  33      C  ,' LIMITS FOR JET NO.',I3,' ARE',2E12.4)
  34       IF(IMSG.EQ.7) WRITE(ITLIS,1007) I,XJ(I),X1,X2
  35  1007 FORMAT(//5X,'FIXED X VALUE FOR JET NO.',I3,' IS',E12.4,2X,
  36      C  'THIS IS INCOMPATIBLE WITH ALLOWED X LIMITS',2E12.4)
  37 C
  38 C           ERRORS IN W(Z0) PARAMETERS
  39 C
  40       IF(IMSG.EQ.101) WRITE(ITLIS,901) XW,XWMIN,XWMAX
  41   901 FORMAT(//5X,'CHOICE OF PARAMETERS GIVES A FIXED XW',E12.4,
  42      C  ' ,THIS VALUE IS INCOMPATIBLE WITH THE LIMITS',2E12.4)
  43       IF(IMSG.EQ.102) WRITE(ITLIS,902) YW,YWMIN,YWMAX
  44   902 FORMAT(//5X,'CHOICE OF PARAMETERS GIVES A FIXED YW',
  45      C  E12.4,' ,THIS VALUE IS INCOMPATIBLE WITH THE LIMITS ')
  46       IF(IMSG.EQ.103) WRITE(ITLIS,903) QMW,QMIN,QMAX
  47   903 FORMAT(//5X,'CHOICE OF PARAMETERS GIVES A FIXED QMW',
  48      C  E12.4,' ,THIS VALUE IS INCOMPATIBLE WITH THE LIMITS',
  49      C  E12.4)
  50       IF(IMSG.EQ.104) WRITE(ITLIS,904) XW,YW,QTW
  51   904 FORMAT(//5X,'FIXED VALUES FOR XW,YW,AND QTW',3E12.4,
  52      C  ' ARE UNPHYSICAL')
  53       IF(IMSG.EQ.105) WRITE(ITLIS,905) QTW,QTMIN,QTMAX
  54   905 FORMAT(//5X,'CHOICE OF PARAMETERS GIVES A FIXED QTW',E12.4
  55      C  ,' ,THIS VALUE IS INCOMPATIBLE WITH THE LIMITS',2E12.4)
  56       IF(IMSG.EQ.106) WRITE(ITLIS,906) XW,YW,QMW
  57   906 FORMAT(//5X,'FIXED VALUS FOR XW,YW,AND QMW',3E12.4,
  58      C  ' ARE UNPHYSICAL')
  59       IF(IMSG.EQ.107) WRITE(ITLIS,907) QTMIN,QTMAX
  60   907 FORMAT(//5X,'BAD LIMITS FOR QTW',2E12.4)
  61       IF(IMSG.EQ.108) WRITE(ITLIS,908) QMIN,QMAX
  62   908 FORMAT(//5X,'BAD LIMITS FOR QMW',2E12.4)
  63       IF(IMSG.EQ.109) WRITE(ITLIS,909) THWMIN,THWMAX
  64   909 FORMAT(//5X,'BAD LIMITS FOR THW',2E12.4,2X,' REMEMBER TH MUST',
  65      C  ' BE IN RADIANS AND LIE BETWEEN 0 AND PI')
  66       IF(IMSG.EQ.110) WRITE(ITLIS,910) PHWMIN,PHWMAX
  67   910 FORMAT(//5X,'BAD LIMITS FOR PHW',2E12.4,' ,REMEMBER PHW MUST',
  68      C  ' BE IN RADIANS AND PHMAX-PHMIN MUST BE LESS THAN 2PI')
  69       IF(IMSG.EQ.111) WRITE(ITLIS,911) XWMIN,XWMAX
  70   911 FORMAT(//5X,'BAD LIMITS FOR XW',2E12.4)
  71       IF(IMSG.EQ.112) WRITE(ITLIS,912) YWMIN,YWMAX
  72   912 FORMAT(//5X,'BAD LIMITS FOR YW',2E12.4)
  73       IF(IMSG.EQ.113) WRITE(ITLIS,913)
  74   913 FORMAT(//5X,'SORRY, BUT YOU CANNOT FIX THETA FOR DRELLYAN EVENTS.'
  75      C,'  THINK OF SOMETHING ELSE.')
  76       IF(IMSG.EQ.114) WRITE(ITLIS,914)
  77   914 FORMAT(//5X,'YOU CANNOT FIX PARAMETERS FOR THE DECAY OF A',
  78      C  ' DRELL YAN JET')
  79       IF(IMSG.EQ.115) WRITE(ITLIS,915)
  80   915 FORMAT(//5X,'YOU CANNOT FIX QTW,QMW,YW AND XW SIMULTANEUOSLY')
  81 C
  82 C       ERRORS IN E+E- PARAMETERS
  83 C
  84       IF(IMSG.EQ.116)
  85      1WRITE(ITLIS,631) THMIN(1),THMAX(1),THMIN(2),THMAX(2)
  86 631   FORMAT(//10X,'THETA LIMITS',2E12.4,' FOR JET 1 AND',2E12.4
  87      C  ,' FOR JET 2 ARE INCOMPATIBLE')
  88 C
  89       RETURN
  90       END