SUBROUTINE SUSCHK(ok) C ----------------------------------------------------------------------------- C Purpose: Check the starting parameters for SUSANA. C Author: JMB 19/03/92 C ----------------------------------------------------------------------------- IMPLICIT NONE C INCLUDE 'SUSANA.INC' C INTEGER i, check LOGICAL ok C ok = .TRUE. C -- Check that the photino is the LSP. DO 1 i = 1, 7 IF ( msus(i).LE.msus(8) ) THEN WRITE (isustr(1),*) ' Photino must be the LSP!' ok = .FALSE. RETURN ENDIF 1 CONTINUE C IF (mp2.GE.epcms) THEN WRITE(isustr(1),9004) ok = .FALSE. RETURN ENDIF C C -- Check that photino decay is kinematically allowed. IF ( msus(8)**2 .LT.(mdd2+msin2 +2.0*sqrt(mdd2*msin2)))THEN WRITE (isustr(1),9001) WRITE (isustr(1),9006) isustr(6) = MOD(isustr(6),100) ENDIF C IF (msus(8)**2 .LT.(msin2+mdu2 +2.0*sqrt(msin2*mdu2)))THEN WRITE (isustr(1),9003) WRITE (isustr(1),9006) isustr(6) = 100*(isustr(6)/100) ENDIF C IF (isustr(6).EQ.0) THEN ok = .FALSE. WRITE (isustr(1),9002) ENDIF C RETURN 9001 FORMAT(8X, & 'Photino -> real(q q) nu decay channel is forbidden!') 9002 FORMAT(8X,' NO ALLOWED PHOTINO DECAYS!') 9003 FORMAT(8X, & 'Photino -> real(q Q) e+ decay channel is forbidden!') 9004 FORMAT(8X, & 'Real photino production forbidden over whole kinematic' &, ' range!') 9006 FORMAT & (8X,'Branching fraction for this channel' &, ' will be assumed zero.') END