Back to home page

sPhenix code displayed by LXR
 
 

    

File indexing completed on 2025-08-05 08:21:18 

0001  
0002  
0003 C***********************************************************************
0004  
0005 C...PYSTAT
0006 C...Prints out information about cross-sections, decay widths, branching
0007 C...ratios, kinematical limits, status codes and parameter values.
0008  
0009       SUBROUTINE PYSTAT(MSTAT)
0010  
0011 C...Double precision and integer declarations.
0012       IMPLICIT DOUBLE PRECISION(A-H, O-Z)
0013       IMPLICIT INTEGER(I-N)
0014       INTEGER PYK,PYCHGE,PYCOMP
0015 C...Parameter statement to help give large particle numbers.
0016       PARAMETER (KSUSY1=1000000,KSUSY2=2000000,KTECHN=3000000,
0017      &KEXCIT=4000000,KDIMEN=5000000)
0018       PARAMETER (EPS=1D-3)
0019 C...Commonblocks.
0020       COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
0021       COMMON/PYDAT2/KCHG(500,4),PMAS(500,4),PARF(2000),VCKM(4,4)
0022       COMMON/PYDAT3/MDCY(500,3),MDME(8000,2),BRAT(8000),KFDP(8000,5)
0023       COMMON/PYSUBS/MSEL,MSELPD,MSUB(500),KFIN(2,-40:40),CKIN(200)
0024       COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
0025       COMMON/PYINT1/MINT(400),VINT(400)
0026       COMMON/PYINT2/ISET(500),KFPR(500,2),COEF(500,20),ICOL(40,4,2)
0027       COMMON/PYINT4/MWID(500),WIDS(500,5)
0028       COMMON/PYINT5/NGENPD,NGEN(0:500,3),XSEC(0:500,3)
0029       COMMON/PYINT6/PROC(0:500)
0030       CHARACTER PROC*28, CHTMP*16
0031       COMMON/PYMSSM/IMSS(0:99),RMSS(0:99)
0032       COMMON/PYMSRV/RVLAM(3,3,3), RVLAMP(3,3,3), RVLAMB(3,3,3)
0033       SAVE /PYDAT1/,/PYDAT2/,/PYDAT3/,/PYSUBS/,/PYPARS/,/PYINT1/,
0034      &/PYINT2/,/PYINT4/,/PYINT5/,/PYINT6/,/PYMSSM/,/PYMSRV/
0035 C...Local arrays, character variables and data.
0036       DIMENSION WDTP(0:400),WDTE(0:400,0:5),NMODES(0:20),PBRAT(10)
0037       CHARACTER PROGA(6)*28,CHAU*16,CHKF*16,CHD1*16,CHD2*16,CHD3*16,
0038      &CHIN(2)*12,STATE(-1:5)*4,CHKIN(21)*18,DISGA(2)*28,
0039      &PROGG9(13)*28,PROGG4(4)*28,PROGG2(2)*28,PROGP4(4)*28
0040       CHARACTER*24 CHD0, CHDC(10)
0041       CHARACTER*6 DNAME(3)
0042       DATA PROGA/
0043      &'VMD/hadron * VMD            ','VMD/hadron * direct         ',
0044      &'VMD/hadron * anomalous      ','direct * direct             ',
0045      &'direct * anomalous          ','anomalous * anomalous       '/
0046       DATA DISGA/'e * VMD','e * anomalous'/
0047       DATA PROGG9/
0048      &'direct * direct             ','direct * VMD                ',
0049      &'direct * anomalous          ','VMD * direct                ',
0050      &'VMD * VMD                   ','VMD * anomalous             ',
0051      &'anomalous * direct          ','anomalous * VMD             ',
0052      &'anomalous * anomalous       ','DIS * VMD                   ',
0053      &'DIS * anomalous             ','VMD * DIS                   ',
0054      &'anomalous * DIS             '/
0055       DATA PROGG4/
0056      &'direct * direct             ','direct * resolved           ',
0057      &'resolved * direct           ','resolved * resolved         '/
0058       DATA PROGG2/
0059      &'direct * hadron             ','resolved * hadron           '/
0060       DATA PROGP4/
0061      &'VMD * hadron                ','direct * hadron             ',
0062      &'anomalous * hadron          ','DIS * hadron                '/
0063       DATA STATE/'----','off ','on  ','on/+','on/-','on/1','on/2'/,
0064      &CHKIN/' m_hard (GeV/c^2) ',' p_T_hard (GeV/c) ',
0065      &'m_finite (GeV/c^2)','   y*_subsystem   ','     y*_large     ',
0066      &'     y*_small     ','    eta*_large    ','    eta*_small    ',
0067      &'cos(theta*)_large ','cos(theta*)_small ','       x_1        ',
0068      &'       x_2        ','       x_F        ',' cos(theta_hard)  ',
0069      &'m''_hard (GeV/c^2) ','       tau        ','        y*        ',
0070      &'cos(theta_hard^-) ','cos(theta_hard^+) ','      x_T^2       ',
0071      &'       tau''       '/
0072       DATA DNAME /'q     ','lepton','nu    '/
0073  
0074 C...Cross-sections.
0075       IF(MSTAT.LE.1) THEN
0076         IF(MINT(121).GT.1) CALL PYSAVE(5,0)
0077         WRITE(MSTU(11),5000)
0078         WRITE(MSTU(11),5100)
0079         WRITE(MSTU(11),5200) 0,PROC(0),NGEN(0,3),NGEN(0,1),XSEC(0,3)
0080         DO 100 I=1,500
0081           IF(MSUB(I).NE.1) GOTO 100
0082           WRITE(MSTU(11),5200) I,PROC(I),NGEN(I,3),NGEN(I,1),XSEC(I,3)
0083   100   CONTINUE
0084         IF(MINT(121).GT.1) THEN
0085           WRITE(MSTU(11),5300)
0086           DO 110 IGA=1,MINT(121)
0087             CALL PYSAVE(3,IGA)
0088             IF(MINT(121).EQ.2.AND.MSTP(14).EQ.10) THEN
0089               WRITE(MSTU(11),5200) IGA,DISGA(IGA),NGEN(0,3),NGEN(0,1),
0090      &        XSEC(0,3)
0091             ELSEIF(MINT(121).EQ.9.OR.MINT(121).EQ.13) THEN
0092               WRITE(MSTU(11),5200) IGA,PROGG9(IGA),NGEN(0,3),NGEN(0,1),
0093      &        XSEC(0,3)
0094             ELSEIF(MINT(121).EQ.4.AND.MSTP(14).EQ.30) THEN
0095               WRITE(MSTU(11),5200) IGA,PROGP4(IGA),NGEN(0,3),NGEN(0,1),
0096      &        XSEC(0,3)
0097             ELSEIF(MINT(121).EQ.4) THEN
0098               WRITE(MSTU(11),5200) IGA,PROGG4(IGA),NGEN(0,3),NGEN(0,1),
0099      &        XSEC(0,3)
0100             ELSEIF(MINT(121).EQ.2) THEN
0101               WRITE(MSTU(11),5200) IGA,PROGG2(IGA),NGEN(0,3),NGEN(0,1),
0102      &        XSEC(0,3)
0103             ELSE
0104               WRITE(MSTU(11),5200) IGA,PROGA(IGA),NGEN(0,3),NGEN(0,1),
0105      &        XSEC(0,3)
0106             ENDIF
0107   110     CONTINUE
0108           CALL PYSAVE(5,0)
0109         ENDIF
0110         WRITE(MSTU(11),5400) MSTU(23),MSTU(30),MSTU(27),
0111      &  1D0-DBLE(NGEN(0,3))/MAX(1D0,DBLE(NGEN(0,2)))
0112  
0113 C...Decay widths and branching ratios.
0114       ELSEIF(MSTAT.EQ.2) THEN
0115         WRITE(MSTU(11),5500)
0116         WRITE(MSTU(11),5600)
0117         DO 140 KC=1,500
0118           KF=KCHG(KC,4)
0119           CALL PYNAME(KF,CHKF)
0120           IOFF=0
0121           IF(KC.LE.22) THEN
0122             IF(KC.GT.2*MSTP(1).AND.KC.LE.10) GOTO 140
0123             IF(KC.GT.10+2*MSTP(1).AND.KC.LE.20) GOTO 140
0124             IF(KC.LE.5.OR.(KC.GE.11.AND.KC.LE.16)) IOFF=1
0125             IF(KC.EQ.18.AND.PMAS(18,1).LT.1D0) IOFF=1
0126             IF(KC.EQ.21.OR.KC.EQ.22) IOFF=1
0127           ELSE
0128             IF(MWID(KC).LE.0) GOTO 140
0129             IF(IMSS(1).LE.0.AND.(KF/KSUSY1.EQ.1.OR.
0130      &      KF/KSUSY1.EQ.2)) GOTO 140
0131           ENDIF
0132 C...Off-shell branchings.
0133           IF(IOFF.EQ.1) THEN
0134             NGP=0
0135             IF(KC.LE.20) NGP=(MOD(KC,10)+1)/2
0136             IF(NGP.LE.MSTP(1)) WRITE(MSTU(11),5700) KF,CHKF(1:10),
0137      &      PMAS(KC,1),0D0,0D0,STATE(MDCY(KC,1)),0D0
0138             DO 120 J=1,MDCY(KC,3)
0139               IDC=J+MDCY(KC,2)-1
0140               NGP1=0
0141               IF(IABS(KFDP(IDC,1)).LE.20) NGP1=
0142      &        (MOD(IABS(KFDP(IDC,1)),10)+1)/2
0143               NGP2=0
0144               IF(IABS(KFDP(IDC,2)).LE.20) NGP2=
0145      &        (MOD(IABS(KFDP(IDC,2)),10)+1)/2
0146               CALL PYNAME(KFDP(IDC,1),CHD1)
0147               CALL PYNAME(KFDP(IDC,2),CHD2)
0148               IF(KFDP(IDC,3).EQ.0) THEN
0149                 IF(MDME(IDC,2).EQ.102.AND.NGP1.LE.MSTP(1).AND.
0150      &          NGP2.LE.MSTP(1)) WRITE(MSTU(11),5800) IDC,CHD1(1:10),
0151      &          CHD2(1:10),0D0,0D0,STATE(MDME(IDC,1)),0D0
0152               ELSE
0153                 CALL PYNAME(KFDP(IDC,3),CHD3)
0154                 IF(MDME(IDC,2).EQ.102.AND.NGP1.LE.MSTP(1).AND.
0155      &          NGP2.LE.MSTP(1)) WRITE(MSTU(11),5900) IDC,CHD1(1:10),
0156      &          CHD2(1:10),CHD3(1:10),0D0,0D0,STATE(MDME(IDC,1)),0D0
0157               ENDIF
0158   120       CONTINUE
0159 C...On-shell decays.
0160           ELSE
0161             CALL PYWIDT(KF,PMAS(KC,1)**2,WDTP,WDTE)
0162             BRFIN=1D0
0163             IF(WDTE(0,0).LE.0D0) BRFIN=0D0
0164             WRITE(MSTU(11),5700) KF,CHKF(1:10),PMAS(KC,1),WDTP(0),1D0,
0165      &      STATE(MDCY(KC,1)),BRFIN
0166             DO 130 J=1,MDCY(KC,3)
0167               IDC=J+MDCY(KC,2)-1
0168               NGP1=0
0169               IF(IABS(KFDP(IDC,1)).LE.20) NGP1=
0170      &        (MOD(IABS(KFDP(IDC,1)),10)+1)/2
0171               NGP2=0
0172               IF(IABS(KFDP(IDC,2)).LE.20) NGP2=
0173      &        (MOD(IABS(KFDP(IDC,2)),10)+1)/2
0174               BRPRI=0D0
0175               IF(WDTP(0).GT.0D0) BRPRI=WDTP(J)/WDTP(0)
0176               BRFIN=0D0
0177               IF(WDTE(0,0).GT.0D0) BRFIN=WDTE(J,0)/WDTE(0,0)
0178               CALL PYNAME(KFDP(IDC,1),CHD1)
0179               CALL PYNAME(KFDP(IDC,2),CHD2)
0180               IF(KFDP(IDC,3).EQ.0) THEN
0181                 IF(NGP1.LE.MSTP(1).AND.NGP2.LE.MSTP(1))
0182      &          WRITE(MSTU(11),5800) IDC,CHD1(1:10),
0183      &          CHD2(1:10),WDTP(J),BRPRI,
0184      &          STATE(MDME(IDC,1)),BRFIN
0185               ELSE
0186                 CALL PYNAME(KFDP(IDC,3),CHD3)
0187                 IF(NGP1.LE.MSTP(1).AND.NGP2.LE.MSTP(1))
0188      &          WRITE(MSTU(11),5900) IDC,CHD1(1:10),
0189      &          CHD2(1:10),CHD3(1:10),WDTP(J),BRPRI,
0190      &          STATE(MDME(IDC,1)),BRFIN
0191               ENDIF
0192   130       CONTINUE
0193           ENDIF
0194   140   CONTINUE
0195         WRITE(MSTU(11),6000)
0196  
0197 C...Allowed incoming partons/particles at hard interaction.
0198       ELSEIF(MSTAT.EQ.3) THEN
0199         WRITE(MSTU(11),6100)
0200         CALL PYNAME(MINT(11),CHAU)
0201         CHIN(1)=CHAU(1:12)
0202         CALL PYNAME(MINT(12),CHAU)
0203         CHIN(2)=CHAU(1:12)
0204         WRITE(MSTU(11),6200) CHIN(1),CHIN(2)
0205         DO 150 I=-20,22
0206           IF(I.EQ.0) GOTO 150
0207           IA=IABS(I)
0208           IF(IA.GT.MSTP(58).AND.IA.LE.10) GOTO 150
0209           IF(IA.GT.10+2*MSTP(1).AND.IA.LE.20) GOTO 150
0210           CALL PYNAME(I,CHAU)
0211           WRITE(MSTU(11),6300) CHAU,STATE(KFIN(1,I)),CHAU,
0212      &    STATE(KFIN(2,I))
0213   150   CONTINUE
0214         WRITE(MSTU(11),6400)
0215  
0216 C...User-defined limits on kinematical variables.
0217       ELSEIF(MSTAT.EQ.4) THEN
0218         WRITE(MSTU(11),6500)
0219         WRITE(MSTU(11),6600)
0220         SHRMAX=CKIN(2)
0221         IF(SHRMAX.LT.0D0) SHRMAX=VINT(1)
0222         WRITE(MSTU(11),6700) CKIN(1),CHKIN(1),SHRMAX
0223         PTHMIN=MAX(CKIN(3),CKIN(5))
0224         PTHMAX=CKIN(4)
0225         IF(PTHMAX.LT.0D0) PTHMAX=0.5D0*SHRMAX
0226         WRITE(MSTU(11),6800) CKIN(3),PTHMIN,CHKIN(2),PTHMAX
0227         WRITE(MSTU(11),6900) CHKIN(3),CKIN(6)
0228         DO 160 I=4,14
0229           WRITE(MSTU(11),6700) CKIN(2*I-1),CHKIN(I),CKIN(2*I)
0230   160   CONTINUE
0231         SPRMAX=CKIN(32)
0232         IF(SPRMAX.LT.0D0) SPRMAX=VINT(1)
0233         WRITE(MSTU(11),6700) CKIN(31),CHKIN(15),SPRMAX
0234         WRITE(MSTU(11),7000)
0235  
0236 C...Status codes and parameter values.
0237       ELSEIF(MSTAT.EQ.5) THEN
0238         WRITE(MSTU(11),7100)
0239         WRITE(MSTU(11),7200)
0240         DO 170 I=1,100
0241           WRITE(MSTU(11),7300) I,MSTP(I),PARP(I),100+I,MSTP(100+I),
0242      &    PARP(100+I)
0243   170   CONTINUE
0244  
0245 C...List of all processes implemented in the program.
0246       ELSEIF(MSTAT.EQ.6) THEN
0247         WRITE(MSTU(11),7400)
0248         WRITE(MSTU(11),7500)
0249         DO 180 I=1,500
0250           IF(ISET(I).LT.0) GOTO 180
0251           WRITE(MSTU(11),7600) I,PROC(I),ISET(I),KFPR(I,1),KFPR(I,2)
0252   180   CONTINUE
0253         WRITE(MSTU(11),7700)
0254  
0255       ELSEIF(MSTAT.EQ.7) THEN
0256       WRITE (MSTU(11),8000)
0257       NMODES(0)=0
0258       NMODES(10)=0
0259       NMODES(9)=0
0260       DO 290 ILR=1,2
0261         DO 280 KFSM=1,16
0262           KFSUSY=ILR*KSUSY1+KFSM
0263           NRVDC=0
0264 C...SDOWN DECAYS
0265           IF (KFSM.EQ.1.OR.KFSM.EQ.3.OR.KFSM.EQ.5) THEN
0266             NRVDC=3
0267             DO 190 I=1,NRVDC
0268               PBRAT(I)=0D0
0269               NMODES(I)=0
0270   190       CONTINUE
0271             CALL PYNAME(KFSUSY,CHTMP)
0272             CHD0=CHTMP//' '
0273             CHDC(1)=DNAME(3) // ' + ' // DNAME(1)
0274             CHDC(2)=DNAME(2) // ' + ' // DNAME(1)
0275             CHDC(3)=DNAME(1) // ' + ' // DNAME(1)
0276             KC=PYCOMP(KFSUSY)
0277             DO 200 J=1,MDCY(KC,3)
0278               IDC=J+MDCY(KC,2)-1
0279               ID1=IABS(KFDP(IDC,1))
0280               ID2=IABS(KFDP(IDC,2))
0281               IF (KFDP(IDC,3).EQ.0) THEN
0282                 IF ((ID1.EQ.12.OR.ID1.EQ.14.OR.ID1.EQ.16).AND.(ID2
0283      &               .EQ.1.OR.ID2.EQ.3.OR.ID2.EQ.5)) THEN
0284                   PBRAT(1)=PBRAT(1)+BRAT(IDC)
0285                   NMODES(1)=NMODES(1)+1
0286                   IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0287                   IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0288                 ELSE IF ((ID1.EQ.11.OR.ID1.EQ.13.OR.ID1.EQ.15).AND
0289      &                 .(ID2.EQ.2.OR.ID2.EQ.4.OR.ID2.EQ.6)) THEN
0290                   PBRAT(2)=PBRAT(2)+BRAT(IDC)
0291                   NMODES(2)=NMODES(2)+1
0292                   IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0293                   IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0294                 ELSE IF ((ID1.EQ.2.OR.ID1.EQ.4.OR.ID1.EQ.6).AND
0295      &                 .(ID2.EQ.1.OR.ID2.EQ.3.OR.ID2.EQ.5)) THEN
0296                   PBRAT(3)=PBRAT(3)+BRAT(IDC)
0297                   NMODES(3)=NMODES(3)+1
0298                   IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0299                   IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0300                 ENDIF
0301               ENDIF
0302   200       CONTINUE
0303           ENDIF
0304 C...SUP DECAYS
0305           IF (KFSM.EQ.2.OR.KFSM.EQ.4.OR.KFSM.EQ.6) THEN
0306             NRVDC=2
0307             DO 210 I=1,NRVDC
0308               NMODES(I)=0
0309               PBRAT(I)=0D0
0310   210       CONTINUE
0311             CALL PYNAME(KFSUSY,CHTMP)
0312             CHD0=CHTMP//' '
0313             CHDC(1)=DNAME(2) // ' + ' // DNAME(1)
0314             CHDC(2)=DNAME(1) // ' + ' // DNAME(1)
0315             KC=PYCOMP(KFSUSY)
0316             DO 220 J=1,MDCY(KC,3)
0317               IDC=J+MDCY(KC,2)-1
0318               ID1=IABS(KFDP(IDC,1))
0319               ID2=IABS(KFDP(IDC,2))
0320               IF (KFDP(IDC,3).EQ.0) THEN
0321                 IF ((ID1.EQ.11.OR.ID1.EQ.13.OR.ID1.EQ.15).AND.(ID2
0322      &               .EQ.1.OR.ID2.EQ.3.OR.ID2.EQ.5)) THEN
0323                   PBRAT(1)=PBRAT(1)+BRAT(IDC)
0324                   NMODES(1)=NMODES(1)+1
0325                   IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0326                   IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0327                 ELSE IF ((ID1.EQ.1.OR.ID1.EQ.3.OR.ID1.EQ.5).AND.(ID2
0328      &               .EQ.1.OR.ID2.EQ.3.OR.ID2.EQ.5)) THEN
0329                   PBRAT(2)=PBRAT(2)+BRAT(IDC)
0330                   NMODES(2)=NMODES(2)+1
0331                   IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0332                   IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0333                 ENDIF
0334               ENDIF
0335   220       CONTINUE
0336           ENDIF
0337 C...SLEPTON DECAYS
0338           IF (KFSM.EQ.11.OR.KFSM.EQ.13.OR.KFSM.EQ.15) THEN
0339             NRVDC=2
0340             DO 230 I=1,NRVDC
0341               PBRAT(I)=0D0
0342               NMODES(I)=0
0343   230       CONTINUE
0344             CALL PYNAME(KFSUSY,CHTMP)
0345             CHD0=CHTMP//' '
0346             CHDC(1)=DNAME(3) // ' + ' // DNAME(2)
0347             CHDC(2)=DNAME(1) // ' + ' // DNAME(1)
0348             KC=PYCOMP(KFSUSY)
0349             DO 240 J=1,MDCY(KC,3)
0350               IDC=J+MDCY(KC,2)-1
0351               ID1=IABS(KFDP(IDC,1))
0352               ID2=IABS(KFDP(IDC,2))
0353               IF (KFDP(IDC,3).EQ.0) THEN
0354                 IF ((ID1.EQ.12.OR.ID1.EQ.14.OR.ID1.EQ.16).AND.(ID2
0355      &               .EQ.11.OR.ID2.EQ.13.OR.ID2.EQ.15)) THEN
0356                   PBRAT(1)=PBRAT(1)+BRAT(IDC)
0357                   NMODES(1)=NMODES(1)+1
0358                   IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0359                   IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0360                 ENDIF
0361                 IF ((ID1.EQ.2.OR.ID1.EQ.4.OR.ID1.EQ.6).AND.(ID2
0362      &               .EQ.1.OR.ID2.EQ.3.OR.ID2.EQ.5)) THEN
0363                   PBRAT(2)=PBRAT(2)+BRAT(IDC)
0364                   NMODES(2)=NMODES(2)+1
0365                   IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0366                   IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0367                 ENDIF
0368               ENDIF
0369   240       CONTINUE
0370           ENDIF
0371 C...SNEUTRINO DECAYS
0372           IF ((KFSM.EQ.12.OR.KFSM.EQ.14.OR.KFSM.EQ.16).AND.ILR.EQ.1)
0373      &         THEN
0374             NRVDC=2
0375             DO 250 I=1,NRVDC
0376               PBRAT(I)=0D0
0377               NMODES(I)=0
0378   250       CONTINUE
0379             CALL PYNAME(KFSUSY,CHTMP)
0380             CHD0=CHTMP//' '
0381             CHDC(1)=DNAME(2) // ' + ' // DNAME(2)
0382             CHDC(2)=DNAME(1) // ' + ' // DNAME(1)
0383             KC=PYCOMP(KFSUSY)
0384             DO 260 J=1,MDCY(KC,3)
0385               IDC=J+MDCY(KC,2)-1
0386               ID1=IABS(KFDP(IDC,1))
0387               ID2=IABS(KFDP(IDC,2))
0388               IF (KFDP(IDC,3).EQ.0) THEN
0389                 IF ((ID1.EQ.11.OR.ID1.EQ.13.OR.ID1.EQ.15).AND.(ID2
0390      &               .EQ.11.OR.ID2.EQ.13.OR.ID2.EQ.15)) THEN
0391                   PBRAT(1)=PBRAT(1)+BRAT(IDC)
0392                   NMODES(1)=NMODES(1)+1
0393                   IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0394                   IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0395                 ENDIF
0396                 IF ((ID1.EQ.1.OR.ID1.EQ.3.OR.ID1.EQ.5).AND.(ID2
0397      &               .EQ.1.OR.ID2.EQ.3.OR.ID2.EQ.5)) THEN
0398                   NMODES(2)=NMODES(2)+1
0399                   PBRAT(2)=PBRAT(2)+BRAT(IDC)
0400                   IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0401                   IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0402                 ENDIF
0403               ENDIF
0404   260       CONTINUE
0405           ENDIF
0406           IF (NRVDC.NE.0) THEN
0407             DO 270 I=1,NRVDC
0408               WRITE (MSTU(11),8200) CHD0, CHDC(I), PBRAT(I), NMODES(I)
0409               NMODES(0)=NMODES(0)+NMODES(I)
0410   270       CONTINUE
0411           ENDIF
0412   280   CONTINUE
0413   290 CONTINUE
0414       DO 370 KFSM=21,37
0415         KFSUSY=KSUSY1+KFSM
0416         NRVDC=0
0417 C...NEUTRALINO DECAYS
0418         IF (KFSM.EQ.22.OR.KFSM.EQ.23.OR.KFSM.EQ.25.OR.KFSM.EQ.35) THEN
0419           NRVDC=4
0420           DO 300 I=1,NRVDC
0421             PBRAT(I)=0D0
0422             NMODES(I)=0
0423   300     CONTINUE
0424           CALL PYNAME(KFSUSY,CHTMP)
0425           CHD0=CHTMP//' '
0426           CHDC(1)=DNAME(3) // ' + ' // DNAME(2) // ' + ' // DNAME(2)
0427           CHDC(2)=DNAME(3) // ' + ' // DNAME(1) // ' + ' // DNAME(1)
0428           CHDC(3)=DNAME(2) // ' + ' // DNAME(1) // ' + ' // DNAME(1)
0429           CHDC(4)=DNAME(1) // ' + ' // DNAME(1) // ' + ' // DNAME(1)
0430           KC=PYCOMP(KFSUSY)
0431           DO 310 J=1,MDCY(KC,3)
0432             IDC=J+MDCY(KC,2)-1
0433             ID1=IABS(KFDP(IDC,1))
0434             ID2=IABS(KFDP(IDC,2))
0435             ID3=IABS(KFDP(IDC,3))
0436             IF ((ID1.EQ.12.OR.ID1.EQ.14.OR.ID1.EQ.16).AND.(ID2
0437      &           .EQ.11.OR.ID2.EQ.13.OR.ID2.EQ.15).AND.(ID3.EQ.11.OR
0438      &           .ID3.EQ.13.OR.ID3.EQ.15)) THEN
0439               PBRAT(1)=PBRAT(1)+BRAT(IDC)
0440               NMODES(1)=NMODES(1)+1
0441               IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0442               IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0443             ELSE IF ((ID1.EQ.12.OR.ID1.EQ.14.OR.ID1.EQ.16).AND
0444      &             .(ID2.EQ.1.OR.ID2.EQ.3.OR.ID2.EQ.5).AND.(ID3.EQ.1
0445      &             .OR.ID3.EQ.3.OR.ID3.EQ.5)) THEN
0446               PBRAT(2)=PBRAT(2)+BRAT(IDC)
0447               NMODES(2)=NMODES(2)+1
0448               IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0449               IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0450             ELSE IF ((ID1.EQ.11.OR.ID1.EQ.13.OR.ID1.EQ.15).AND
0451      &             .(ID2.EQ.2.OR.ID2.EQ.4.OR.ID2.EQ.6).AND.(ID3.EQ.1
0452      &             .OR.ID3.EQ.3.OR.ID3.EQ.5)) THEN
0453               PBRAT(3)=PBRAT(3)+BRAT(IDC)
0454               NMODES(3)=NMODES(3)+1
0455               IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0456               IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0457             ELSE IF ((ID1.EQ.2.OR.ID1.EQ.4.OR.ID1.EQ.6).AND
0458      &             .(ID2.EQ.1.OR.ID2.EQ.3.OR.ID2.EQ.5).AND.(ID3.EQ.1
0459      &             .OR.ID3.EQ.3.OR.ID3.EQ.5)) THEN
0460               PBRAT(4)=PBRAT(4)+BRAT(IDC)
0461               NMODES(4)=NMODES(4)+1
0462               IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0463               IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0464             ENDIF
0465   310     CONTINUE
0466         ENDIF
0467 C...CHARGINO DECAYS
0468         IF (KFSM.EQ.24.OR.KFSM.EQ.37) THEN
0469           NRVDC=5
0470           DO 320 I=1,NRVDC
0471             PBRAT(I)=0D0
0472             NMODES(I)=0
0473   320     CONTINUE
0474           CALL PYNAME(KFSUSY,CHTMP)
0475           CHD0=CHTMP//' '
0476           CHDC(1)=DNAME(3) // ' + ' // DNAME(3) // ' + ' // DNAME(2)
0477           CHDC(2)=DNAME(2) // ' + ' // DNAME(2) // ' + ' // DNAME(2)
0478           CHDC(3)=DNAME(3) // ' + ' // DNAME(1) // ' + ' // DNAME(1)
0479           CHDC(4)=DNAME(2) // ' + ' // DNAME(1) // ' + ' // DNAME(1)
0480           CHDC(5)=DNAME(1) // ' + ' // DNAME(1) // ' + ' // DNAME(1)
0481           KC=PYCOMP(KFSUSY)
0482           DO 330 J=1,MDCY(KC,3)
0483             IDC=J+MDCY(KC,2)-1
0484             ID1=IABS(KFDP(IDC,1))
0485             ID2=IABS(KFDP(IDC,2))
0486             ID3=IABS(KFDP(IDC,3))
0487             IF ((ID1.EQ.12.OR.ID1.EQ.14.OR.ID1.EQ.16).AND.(ID2
0488      &           .EQ.11.OR.ID2.EQ.13.OR.ID2.EQ.15).AND.(ID3.EQ.12.OR
0489      &           .ID3.EQ.14.OR.ID3.EQ.16)) THEN
0490               PBRAT(1)=PBRAT(1)+BRAT(IDC)
0491               NMODES(1)=NMODES(1)+1
0492               IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0493               IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0494             ELSE IF ((ID1.EQ.12.OR.ID1.EQ.14.OR.ID1.EQ.16).AND
0495      &             .(ID2.EQ.12.OR.ID2.EQ.14.OR.ID2.EQ.16).AND.(ID3.EQ
0496      &             .11.OR.ID3.EQ.13.OR.ID3.EQ.15)) THEN
0497               PBRAT(1)=PBRAT(1)+BRAT(IDC)
0498               NMODES(1)=NMODES(1)+1
0499               IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0500               IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0501             ELSE IF ((ID1.EQ.11.OR.ID1.EQ.13.OR.ID1.EQ.15).AND
0502      &             .(ID2.EQ.11.OR.ID2.EQ.13.OR.ID2.EQ.15).AND.(ID3.EQ
0503      &             .11.OR.ID3.EQ.13.OR.ID3.EQ.15)) THEN
0504               PBRAT(2)=PBRAT(2)+BRAT(IDC)
0505               NMODES(2)=NMODES(2)+1
0506               IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0507               IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0508             ELSE IF ((ID1.EQ.12.OR.ID1.EQ.14.OR.ID1.EQ.16).AND
0509      &             .(ID2.EQ.1.OR.ID2.EQ.3.OR.ID2.EQ.5).AND.(ID3.EQ
0510      &             .2.OR.ID3.EQ.4.OR.ID3.EQ.6)) THEN
0511               PBRAT(3)=PBRAT(3)+BRAT(IDC)
0512               NMODES(3)=NMODES(3)+1
0513               IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0514               IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0515             ELSE IF ((ID1.EQ.12.OR.ID1.EQ.14.OR.ID1.EQ.16).AND
0516      &             .(ID2.EQ.2.OR.ID2.EQ.4.OR.ID2.EQ.6).AND.(ID3.EQ
0517      &             .1.OR.ID3.EQ.3.OR.ID3.EQ.5)) THEN
0518               PBRAT(3)=PBRAT(3)+BRAT(IDC)
0519               NMODES(3)=NMODES(3)+1
0520               IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0521               IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0522             ELSE IF ((ID1.EQ.11.OR.ID1.EQ.13.OR.ID1.EQ.15).AND
0523      &             .(ID2.EQ.2.OR.ID2.EQ.4.OR.ID2.EQ.6).AND.(ID3.EQ
0524      &             .2.OR.ID3.EQ.4.OR.ID3.EQ.6)) THEN
0525               PBRAT(4)=PBRAT(4)+BRAT(IDC)
0526               NMODES(4)=NMODES(4)+1
0527               IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0528               IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0529             ELSE IF ((ID1.EQ.11.OR.ID1.EQ.13.OR.ID1.EQ.15).AND
0530      &             .(ID2.EQ.1.OR.ID2.EQ.3.OR.ID2.EQ.5).AND.(ID3.EQ
0531      &             .1.OR.ID3.EQ.3.OR.ID3.EQ.5)) THEN
0532               PBRAT(4)=PBRAT(4)+BRAT(IDC)
0533               NMODES(4)=NMODES(4)+1
0534               IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0535               IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0536             ELSE IF ((ID1.EQ.2.OR.ID1.EQ.4.OR.ID1.EQ.6).AND
0537      &             .(ID2.EQ.2.OR.ID2.EQ.4.OR.ID2.EQ.6).AND.(ID3.EQ
0538      &             .1.OR.ID3.EQ.3.OR.ID3.EQ.5)) THEN
0539               PBRAT(5)=PBRAT(5)+BRAT(IDC)
0540               NMODES(5)=NMODES(5)+1
0541               IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0542               IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0543             ELSE IF ((ID1.EQ.1.OR.ID1.EQ.3.OR.ID1.EQ.5).AND
0544      &             .(ID2.EQ.1.OR.ID2.EQ.3.OR.ID2.EQ.5).AND.(ID3.EQ
0545      &             .1.OR.ID3.EQ.3.OR.ID3.EQ.5)) THEN
0546               PBRAT(5)=PBRAT(5)+BRAT(IDC)
0547               NMODES(5)=NMODES(5)+1
0548               IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0549               IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0550             ENDIF
0551   330     CONTINUE
0552         ENDIF
0553 C...GLUINO DECAYS
0554         IF (KFSM.EQ.21) THEN
0555           NRVDC=3
0556           DO 340 I=1,NRVDC
0557             PBRAT(I)=0D0
0558             NMODES(I)=0
0559   340     CONTINUE
0560           CALL PYNAME(KFSUSY,CHTMP)
0561           CHD0=CHTMP//' '
0562           CHDC(1)=DNAME(3) // ' + ' // DNAME(1) // ' + ' // DNAME(1)
0563           CHDC(2)=DNAME(2) // ' + ' // DNAME(1) // ' + ' // DNAME(1)
0564           CHDC(3)=DNAME(1) // ' + ' // DNAME(1) // ' + ' // DNAME(1)
0565           KC=PYCOMP(KFSUSY)
0566           DO 350 J=1,MDCY(KC,3)
0567             IDC=J+MDCY(KC,2)-1
0568             ID1=IABS(KFDP(IDC,1))
0569             ID2=IABS(KFDP(IDC,2))
0570             ID3=IABS(KFDP(IDC,3))
0571             IF ((ID1.EQ.12.OR.ID1.EQ.14.OR.ID1.EQ.16).AND.(ID2
0572      &           .EQ.1.OR.ID2.EQ.3.OR.ID2.EQ.5).AND.(ID3.EQ.1.OR
0573      &           .ID3.EQ.3.OR.ID3.EQ.5)) THEN
0574               PBRAT(1)=PBRAT(1)+BRAT(IDC)
0575               NMODES(1)=NMODES(1)+1
0576               IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0577               IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0578             ELSE IF ((ID1.EQ.11.OR.ID1.EQ.13.OR.ID1.EQ.15).AND
0579      &             .(ID2.EQ.2.OR.ID2.EQ.4.OR.ID2.EQ.6).AND.(ID3.EQ.1
0580      &             .OR.ID3.EQ.3.OR.ID3.EQ.5)) THEN
0581               PBRAT(2)=PBRAT(2)+BRAT(IDC)
0582               NMODES(2)=NMODES(2)+1
0583               IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0584               IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0585             ELSE IF ((ID1.EQ.2.OR.ID1.EQ.4.OR.ID1.EQ.6).AND
0586      &             .(ID2.EQ.1.OR.ID2.EQ.3.OR.ID2.EQ.5).AND.(ID3.EQ.1
0587      &             .OR.ID3.EQ.3.OR.ID3.EQ.5)) THEN
0588               PBRAT(3)=PBRAT(3)+BRAT(IDC)
0589               NMODES(3)=NMODES(3)+1
0590               IF (BRAT(IDC).GT.0D0) NMODES(10)=NMODES(10)+1
0591               IF (BRAT(IDC).GT.EPS) NMODES(9)=NMODES(9)+1
0592             ENDIF
0593   350     CONTINUE
0594         ENDIF
0595  
0596         IF (NRVDC.NE.0) THEN
0597           DO 360 I=1,NRVDC
0598             WRITE (MSTU(11),8200) CHD0, CHDC(I), PBRAT(I), NMODES(I)
0599             NMODES(0)=NMODES(0)+NMODES(I)
0600   360     CONTINUE
0601         ENDIF
0602   370 CONTINUE
0603       WRITE (MSTU(11),8100) NMODES(0), NMODES(10), NMODES(9)
0604  
0605       IF (IMSS(51).GE.1.OR.IMSS(52).GE.1.OR.IMSS(53).GE.1) THEN
0606         WRITE (MSTU(11),8500)
0607         DO 400 IRV=1,3
0608           DO 390 JRV=1,3
0609             DO 380 KRV=1,3
0610               WRITE (MSTU(11),8700) IRV,JRV,KRV,RVLAM(IRV,JRV,KRV)
0611      &             ,RVLAMP(IRV,JRV,KRV),RVLAMB(IRV,JRV,KRV)
0612   380       CONTINUE
0613   390     CONTINUE
0614   400   CONTINUE
0615         WRITE (MSTU(11),8600)
0616       ENDIF
0617       ENDIF
0618  
0619 C...Formats for printouts.
0620  5000 FORMAT('1',9('*'),1X,'PYSTAT:  Statistics on Number of ',
0621      &'Events and Cross-sections',1X,9('*'))
0622  5100 FORMAT(/1X,78('=')/1X,'I',34X,'I',28X,'I',12X,'I'/1X,'I',12X,
0623      &'Subprocess',12X,'I',6X,'Number of points',6X,'I',4X,'Sigma',3X,
0624      &'I'/1X,'I',34X,'I',28X,'I',12X,'I'/1X,'I',34('-'),'I',28('-'),
0625      &'I',4X,'(mb)',4X,'I'/1X,'I',34X,'I',28X,'I',12X,'I'/1X,'I',1X,
0626      &'N:o',1X,'Type',25X,'I',4X,'Generated',9X,'Tried',1X,'I',12X,
0627      &'I'/1X,'I',34X,'I',28X,'I',12X,'I'/1X,78('=')/1X,'I',34X,'I',28X,
0628      &'I',12X,'I')
0629  5200 FORMAT(1X,'I',1X,I3,1X,A28,1X,'I',1X,I12,1X,I13,1X,'I',1X,1P,
0630      &D10.3,1X,'I')
0631  5300 FORMAT(1X,'I',34X,'I',28X,'I',12X,'I'/1X,78('=')/
0632      &1X,'I',34X,'I',28X,'I',12X,'I')
0633  5400 FORMAT(1X,'I',34X,'I',28X,'I',12X,'I'/1X,78('=')//
0634      &1X,'********* Total number of errors, excluding junctions =',
0635      &1X,I8,' *************'/
0636      &1X,'********* Total number of errors, including junctions =',
0637      &1X,I8,' *************'/
0638      &1X,'********* Total number of warnings =                   ',
0639      &1X,I8,' *************'/
0640      &1X,'********* Fraction of events that fail fragmentation ',
0641      &'cuts =',1X,F8.5,' *********'/)
0642  5500 FORMAT('1',27('*'),1X,'PYSTAT:  Decay Widths and Branching ',
0643      &'Ratios',1X,27('*'))
0644  5600 FORMAT(/1X,98('=')/1X,'I',49X,'I',13X,'I',12X,'I',6X,'I',12X,'I'/
0645      &1X,'I',5X,'Mother  -->  Branching/Decay Channel',8X,'I',1X,
0646      &'Width (GeV)',1X,'I',7X,'B.R.',1X,'I',1X,'Stat',1X,'I',2X,
0647      &'Eff. B.R.',1X,'I'/1X,'I',49X,'I',13X,'I',12X,'I',6X,'I',12X,'I'/
0648      &1X,98('='))
0649  5700 FORMAT(1X,'I',49X,'I',13X,'I',12X,'I',6X,'I',12X,'I'/1X,'I',1X,
0650      &I8,2X,A10,3X,'(m =',F10.3,')',2X,'-->',5X,'I',2X,1P,D10.3,0P,1X,
0651      &'I',1X,1P,D10.3,0P,1X,'I',1X,A4,1X,'I',1X,1P,D10.3,0P,1X,'I')
0652  5800 FORMAT(1X,'I',1X,I8,2X,A10,1X,'+',1X,A10,15X,'I',2X,
0653      &1P,D10.3,0P,1X,'I',1X,1P,D10.3,0P,1X,'I',1X,A4,1X,'I',1X,
0654      &1P,D10.3,0P,1X,'I')
0655  5900 FORMAT(1X,'I',1X,I8,2X,A10,1X,'+',1X,A10,1X,'+',1X,A10,2X,'I',2X,
0656      &1P,D10.3,0P,1X,'I',1X,1P,D10.3,0P,1X,'I',1X,A4,1X,'I',1X,
0657      &1P,D10.3,0P,1X,'I')
0658  6000 FORMAT(1X,'I',49X,'I',13X,'I',12X,'I',6X,'I',12X,'I'/1X,98('='))
0659  6100 FORMAT('1',7('*'),1X,'PYSTAT: Allowed Incoming Partons/',
0660      &'Particles at Hard Interaction',1X,7('*'))
0661  6200 FORMAT(/1X,78('=')/1X,'I',38X,'I',37X,'I'/1X,'I',1X,
0662      &'Beam particle:',1X,A12,10X,'I',1X,'Target particle:',1X,A12,7X,
0663      &'I'/1X,'I',38X,'I',37X,'I'/1X,'I',1X,'Content',6X,'State',19X,
0664      &'I',1X,'Content',6X,'State',18X,'I'/1X,'I',38X,'I',37X,'I'/1X,
0665      &78('=')/1X,'I',38X,'I',37X,'I')
0666  6300 FORMAT(1X,'I',1X,A9,5X,A4,19X,'I',1X,A9,5X,A4,18X,'I')
0667  6400 FORMAT(1X,'I',38X,'I',37X,'I'/1X,78('='))
0668  6500 FORMAT('1',12('*'),1X,'PYSTAT: User-Defined Limits on ',
0669      &'Kinematical Variables',1X,12('*'))
0670  6600 FORMAT(/1X,78('=')/1X,'I',76X,'I')
0671  6700 FORMAT(1X,'I',16X,1P,D10.3,0P,1X,'<',1X,A,1X,'<',1X,1P,D10.3,0P,
0672      &16X,'I')
0673  6800 FORMAT(1X,'I',3X,1P,D10.3,0P,1X,'(',1P,D10.3,0P,')',1X,'<',1X,A,
0674      &1X,'<',1X,1P,D10.3,0P,16X,'I')
0675  6900 FORMAT(1X,'I',29X,A,1X,'=',1X,1P,D10.3,0P,16X,'I')
0676  7000 FORMAT(1X,'I',76X,'I'/1X,78('='))
0677  7100 FORMAT('1',12('*'),1X,'PYSTAT: Summary of Status Codes and ',
0678      &'Parameter Values',1X,12('*'))
0679  7200 FORMAT(/3X,'I',4X,'MSTP(I)',9X,'PARP(I)',20X,'I',4X,'MSTP(I)',9X,
0680      &'PARP(I)'/)
0681  7300 FORMAT(1X,I3,5X,I6,6X,1P,D10.3,0P,18X,I3,5X,I6,6X,1P,D10.3)
0682  7400 FORMAT('1',13('*'),1X,'PYSTAT: List of implemented processes',
0683      &1X,13('*'))
0684  7500 FORMAT(/1X,65('=')/1X,'I',34X,'I',28X,'I'/1X,'I',12X,
0685      &'Subprocess',12X,'I',1X,'ISET',2X,'KFPR(I,1)',2X,'KFPR(I,2)',1X,
0686      &'I'/1X,'I',34X,'I',28X,'I'/1X,65('=')/1X,'I',34X,'I',28X,'I')
0687  7600 FORMAT(1X,'I',1X,I3,1X,A28,1X,'I',1X,I4,1X,I10,1X,I10,1X,'I')
0688  7700 FORMAT(1X,'I',34X,'I',28X,'I'/1X,65('='))
0689  8000 FORMAT(1X/ 1X/
0690      &     17X,'Sums over R-Violating branching ratios',1X/ 1X
0691      &     /1X,70('=')/1X,'I',50X,'I',11X,'I',5X,'I'/1X,'I',4X
0692      &     ,'Mother  -->  Sum over final state flavours',4X,'I',2X
0693      &     ,'BR(sum)',2X,'I',2X,'N',2X,'I'/1X,'I',50X,'I',11X,'I',5X,'I'
0694      &     /1X,70('=')/1X,'I',50X,'I',11X,'I',5X,'I')
0695  8100 FORMAT(1X,'I',50X,'I',11X,'I',5X,'I'/1X,70('=')/1X,'I',1X
0696      &     ,'Total number of R-Violating modes :',3X,I5,24X,'I'/
0697      &     1X,'I',1X,'Total number with non-vanishing BR :',2X,I5,24X
0698      &     ,'I'/1X,'I',1X,'Total number with BR > 0.001 :',8X,I5,24X,'I'
0699      &     /1X,70('='))
0700  8200 FORMAT(1X,'I',1X,A9,1X,'-->',1X,A24,11X,
0701      &     'I',2X,1P,D8.2,0P,1X,'I',2X,I2,1X,'I')
0702  8300 FORMAT(1X,'I',50X,'I',11X,'I',5X,'I')
0703  8500 FORMAT(1X/ 1X/
0704      &     1X,'R-Violating couplings',1X/ 1X /
0705      &     1X,55('=')/
0706      &     1X,'I',1X,'IJK',1X,'I',2X,'LAMBDA(IJK)',2X,'I',2X
0707      &     ,'LAMBDA''(IJK)',1X,'I',1X,"LAMBDA''(IJK)",1X,'I'/1X,'I',5X
0708      &     ,'I',15X,'I',15X,'I',15X,'I')
0709  8600 FORMAT(1X,55('='))
0710  8700 FORMAT(1X,'I',1X,I1,I1,I1,1X,'I',1X,1P,D13.3,0P,1X,'I',1X,1P
0711      &     ,D13.3,0P,1X,'I',1X,1P,D13.3,0P,1X,'I')
0712  
0713       RETURN
0714       END
[ Source navigation ] [ Diff markup ] [ Identifier search ] [ general search ]

This page was automatically generated by the 2.3.7 LXR engine.

The LXR team
Valid CSS 2.1! Valid HTML 4.01!