HIJING/hipyset1_35/pystat_hijing.F

   1 * $Id$
   2
   3 C***********************************************************************
   4
   5       SUBROUTINE PYSTAT_HIJING(MSTAT)
   6
   7 C...Prints out information about cross-sections, decay widths, branching
   8 C...ratios, kinematical limits, status codes and parameter values.
   9 #include "ludat1_hijing.inc"
  10 #include "ludat2_hijing.inc"
  11 #include "ludat3_hijing.inc"
  12 #include "pysubs_hijing.inc"
  13 #include "pypars_hijing.inc"
  14 #include "pyint1_hijing.inc"
  15 #include "pyint4_hijing.inc"
  16 #include "pyint5_hijing.inc"
  17 #include "pyint6_hijing.inc"
  18       CHARACTER CHAU*16,CHPA(-40:40)*12,CHIN(2)*12,
  19      &STATE(-1:5)*4,CHKIN(21)*18
  20       DATA STATE/'----','off ','on  ','on/+','on/-','on/1','on/2'/,
  21      &CHKIN/' m_hard (GeV/c^2) ',' p_T_hard (GeV/c) ',
  22      &'m_finite (GeV/c^2)','   y*_subsystem   ','     y*_large     ',
  23      &'     y*_small     ','    eta*_large    ','    eta*_small    ',
  24      &'cos(theta*)_large ','cos(theta*)_small ','       x_1        ',
  25      &'       x_2        ','       x_F        ',' cos(theta_hard)  ',
  26      &'m''_hard (GeV/c^2) ','       tau        ','        y*        ',
  27      &'cos(theta_hard^-) ','cos(theta_hard^+) ','      x_T^2       ',
  28      &'       tau''       '/
  29
  30 C...Cross-sections.
  31       IF(MSTAT.LE.1) THEN
  32         WRITE(MSTU(11),1000)
  33         WRITE(MSTU(11),1100)
  34         WRITE(MSTU(11),1200) 0,PROC(0),NGEN(0,3),NGEN(0,1),XSEC(0,3)
  35         DO 100 I=1,200
  36         IF(MSUB(I).NE.1) GOTO 100
  37         WRITE(MSTU(11),1200) I,PROC(I),NGEN(I,3),NGEN(I,1),XSEC(I,3)
  38   100   CONTINUE
  39         WRITE(MSTU(11),1300) 1.-FLOAT(NGEN(0,3))/
  40      &  MAX(1.,FLOAT(NGEN(0,2)))
  41
  42 C...Decay widths and branching ratios.
  43       ELSEIF(MSTAT.EQ.2) THEN
  44         DO 110 KF=-40,40
  45         CALL LUNAME_HIJING(KF,CHAU)
  46   110   CHPA(KF)=CHAU(1:12)
  47         WRITE(MSTU(11),1400)
  48         WRITE(MSTU(11),1500)
  49 C...Off-shell branchings.
  50         DO 130 I=1,17
  51         KC=I
  52         IF(I.GE.9) KC=I+2
  53         IF(I.EQ.17) KC=21
  54         WRITE(MSTU(11),1600) CHPA(KC),0.,0.,STATE(MDCY(KC,1)),0.
  55         DO 120 J=1,MDCY(KC,3)
  56         IDC=J+MDCY(KC,2)-1
  57   120   IF(MDME(IDC,2).EQ.102) WRITE(MSTU(11),1700) CHPA(KFDP(IDC,1)),
  58      &  CHPA(KFDP(IDC,2)),0.,0.,STATE(MDME(IDC,1)),0.
  59   130   CONTINUE
  60 C...On-shell decays.
  61         DO 150 I=1,6
  62         KC=I+22
  63         IF(I.EQ.4) KC=32
  64         IF(I.EQ.5) KC=37
  65         IF(I.EQ.6) KC=40
  66         IF(WIDE(KC,0).GT.0.) THEN
  67           WRITE(MSTU(11),1600) CHPA(KC),WIDP(KC,0),1.,
  68      &    STATE(MDCY(KC,1)),1.
  69           DO 140 J=1,MDCY(KC,3)
  70           IDC=J+MDCY(KC,2)-1
  71   140     WRITE(MSTU(11),1700) CHPA(KFDP(IDC,1)),CHPA(KFDP(IDC,2)),
  72      &    WIDP(KC,J),WIDP(KC,J)/WIDP(KC,0),STATE(MDME(IDC,1)),
  73      &    WIDE(KC,J)/WIDE(KC,0)
  74         ELSE
  75           WRITE(MSTU(11),1600) CHPA(KC),WIDP(KC,0),1.,
  76      &    STATE(MDCY(KC,1)),0.
  77         ENDIF
  78   150   CONTINUE
  79         WRITE(MSTU(11),1800)
  80
  81 C...Allowed incoming partons/particles at hard interaction.
  82       ELSEIF(MSTAT.EQ.3) THEN
  83         WRITE(MSTU(11),1900)
  84         CALL LUNAME_HIJING(MINT(11),CHAU)
  85         CHIN(1)=CHAU(1:12)
  86         CALL LUNAME_HIJING(MINT(12),CHAU)
  87         CHIN(2)=CHAU(1:12)
  88         WRITE(MSTU(11),2000) CHIN(1),CHIN(2)
  89         DO 160 KF=-40,40
  90         CALL LUNAME_HIJING(KF,CHAU)
  91   160   CHPA(KF)=CHAU(1:12)
  92         IF(MINT(43).EQ.1) THEN
  93           WRITE(MSTU(11),2100) CHPA(MINT(11)),STATE(KFIN(1,MINT(11))),
  94      &    CHPA(MINT(12)),STATE(KFIN(2,MINT(12)))
  95         ELSEIF(MINT(43).EQ.2) THEN
  96           WRITE(MSTU(11),2100) CHPA(MINT(11)),STATE(KFIN(1,MINT(11))),
  97      &    CHPA(-MSTP(54)),STATE(KFIN(2,-MSTP(54)))
  98           DO 170 I=-MSTP(54)+1,-1
  99   170     WRITE(MSTU(11),2200) CHPA(I),STATE(KFIN(2,I))
 100           DO 180 I=1,MSTP(54)
 101   180     WRITE(MSTU(11),2200) CHPA(I),STATE(KFIN(2,I))
 102           WRITE(MSTU(11),2200) CHPA(21),STATE(KFIN(2,21))
 103         ELSEIF(MINT(43).EQ.3) THEN
 104           WRITE(MSTU(11),2100) CHPA(-MSTP(54)),STATE(KFIN(1,-MSTP(54))),
 105      &    CHPA(MINT(12)),STATE(KFIN(2,MINT(12)))
 106           DO 190 I=-MSTP(54)+1,-1
 107   190     WRITE(MSTU(11),2300) CHPA(I),STATE(KFIN(1,I))
 108           DO 200 I=1,MSTP(54)
 109   200     WRITE(MSTU(11),2300) CHPA(I),STATE(KFIN(1,I))
 110           WRITE(MSTU(11),2300) CHPA(21),STATE(KFIN(1,21))
 111         ELSEIF(MINT(43).EQ.4) THEN
 112           DO 210 I=-MSTP(54),-1
 113   210     WRITE(MSTU(11),2100) CHPA(I),STATE(KFIN(1,I)),CHPA(I),
 114      &    STATE(KFIN(2,I))
 115           DO 220 I=1,MSTP(54)
 116   220     WRITE(MSTU(11),2100) CHPA(I),STATE(KFIN(1,I)),CHPA(I),
 117      &    STATE(KFIN(2,I))
 118           WRITE(MSTU(11),2100) CHPA(21),STATE(KFIN(1,21)),CHPA(21),
 119      &    STATE(KFIN(2,21))
 120         ENDIF
 121         WRITE(MSTU(11),2400)
 122
 123 C...User-defined and derived limits on kinematical variables.
 124       ELSEIF(MSTAT.EQ.4) THEN
 125         WRITE(MSTU(11),2500)
 126         WRITE(MSTU(11),2600)
 127         SHRMAX=CKIN(2)
 128         IF(SHRMAX.LT.0.) SHRMAX=VINT(1)
 129         WRITE(MSTU(11),2700) CKIN(1),CHKIN(1),SHRMAX
 130         PTHMIN=MAX(CKIN(3),CKIN(5))
 131         PTHMAX=CKIN(4)
 132         IF(PTHMAX.LT.0.) PTHMAX=0.5*SHRMAX
 133         WRITE(MSTU(11),2800) CKIN(3),PTHMIN,CHKIN(2),PTHMAX
 134         WRITE(MSTU(11),2900) CHKIN(3),CKIN(6)
 135         DO 230 I=4,14
 136   230   WRITE(MSTU(11),2700) CKIN(2*I-1),CHKIN(I),CKIN(2*I)
 137         SPRMAX=CKIN(32)
 138         IF(SPRMAX.LT.0.) SPRMAX=VINT(1)
 139         WRITE(MSTU(11),2700) CKIN(31),CHKIN(13),SPRMAX
 140         WRITE(MSTU(11),3000)
 141         WRITE(MSTU(11),3100)
 142         WRITE(MSTU(11),2600)
 143         DO 240 I=16,21
 144   240   WRITE(MSTU(11),2700) VINT(I-5),CHKIN(I),VINT(I+15)
 145         WRITE(MSTU(11),3000)
 146
 147 C...Status codes and parameter values.
 148       ELSEIF(MSTAT.EQ.5) THEN
 149         WRITE(MSTU(11),3200)
 150         WRITE(MSTU(11),3300)
 151         DO 250 I=1,100
 152   250   WRITE(MSTU(11),3400) I,MSTP(I),PARP(I),100+I,MSTP(100+I),
 153      &  PARP(100+I)
 154       ENDIF
 155
 156 C...Formats for printouts.
 157  1000 FORMAT('1',9('*'),1X,'PYSTAT_HIJING:  Statistics on Number of ',
 158      &'Events and Cross-sections',1X,9('*'))
 159  1100 FORMAT(/1X,78('=')/1X,'I',34X,'I',28X,'I',12X,'I'/1X,'I',12X,
 160      &'Subprocess',12X,'I',6X,'Number of points',6X,'I',4X,'Sigma',3X,
 161      &'I'/1X,'I',34X,'I',28X,'I',12X,'I'/1X,'I',34('-'),'I',28('-'),
 162      &'I',4X,'(mb)',4X,'I'/1X,'I',34X,'I',28X,'I',12X,'I'/1X,'I',1X,
 163      &'N:o',1X,'Type',25X,'I',4X,'Generated',9X,'Tried',1X,'I',12X,
 164      &'I'/1X,'I',34X,'I',28X,'I',12X,'I'/1X,78('=')/1X,'I',34X,'I',28X,
 165      &'I',12X,'I')
 166  1200 FORMAT(1X,'I',1X,I3,1X,A28,1X,'I',1X,I12,1X,I13,1X,'I',1X,1P,
 167      &E10.3,1X,'I')
 168  1300 FORMAT(1X,'I',34X,'I',28X,'I',12X,'I'/1X,78('=')//
 169      &1X,'********* Fraction of events that fail fragmentation ',
 170      &'cuts =',1X,F8.5,' *********'/)
 171  1400 FORMAT('1',17('*'),1X
 172      $     ,'PYSTAT_HIJING:  Decay Widths and Branching ','Ratios',1X
 173      $     ,17('*'))
 174  1500 FORMAT(/1X,78('=')/1X,'I',29X,'I',13X,'I',12X,'I',6X,'I',12X,'I'/
 175      &1X,'I',1X,'Branching/Decay Channel',5X,'I',1X,'Width (GeV)',1X,
 176      &'I',7X,'B.R.',1X,'I',1X,'Stat',1X,'I',2X,'Eff. B.R.',1X,'I'/1X,
 177      &'I',29X,'I',13X,'I',12X,'I',6X,'I',12X,'I'/1X,78('='))
 178  1600 FORMAT(1X,'I',29X,'I',13X,'I',12X,'I',6X,'I',12X,'I'/1X,'I',1X,
 179      &A12,1X,'->',13X,'I',2X,1P,E10.3,0P,1X,'I',1X,1P,E10.3,0P,1X,'I',
 180      &1X,A4,1X,'I',1X,1P,E10.3,0P,1X,'I')
 181  1700 FORMAT(1X,'I',1X,A12,1X,'+',1X,A12,1X,'I',2X,1P,E10.3,0P,1X,'I',
 182      &1X,1P,E10.3,0P,1X,'I',1X,A4,1X,'I',1X,1P,E10.3,0P,1X,'I')
 183  1800 FORMAT(1X,'I',29X,'I',13X,'I',12X,'I',6X,'I',12X,'I'/1X,78('='))
 184  1900 FORMAT('1',7('*'),1X,'PYSTAT_HIJING: Allowed Incoming Partons/',
 185      &'Particles at Hard Interaction',1X,7('*'))
 186  2000 FORMAT(/1X,78('=')/1X,'I',38X,'I',37X,'I'/1X,'I',1X,
 187      &'Beam particle:',1X,A,10X,'I',1X,'Target particle:',1X,A,7X,
 188      &'I'/1X,'I',38X,'I',37X,'I'/1X,'I',1X,'Content',9X,'State',16X,
 189      &'I',1X,'Content',9X,'State',15X,'I'/1X,'I',38X,'I',37X,'I'/1X,
 190      &78('=')/1X,'I',38X,'I',37X,'I')
 191  2100 FORMAT(1X,'I',1X,A,5X,A,16X,'I',1X,A,5X,A,15X,'I')
 192  2200 FORMAT(1X,'I',38X,'I',1X,A,5X,A,15X,'I')
 193  2300 FORMAT(1X,'I',1X,A,5X,A,16X,'I',37X,'I')
 194  2400 FORMAT(1X,'I',38X,'I',37X,'I'/1X,78('='))
 195  2500 FORMAT('1',12('*'),1X,'PYSTAT_HIJING: User-Defined Limits on ',
 196      &'Kinematical Variables',1X,12('*'))
 197  2600 FORMAT(/1X,78('=')/1X,'I',76X,'I')
 198  2700 FORMAT(1X,'I',16X,1P,E10.3,0P,1X,'<',1X,A,1X,'<',1X,1P,E10.3,0P,
 199      &16X,'I')
 200  2800 FORMAT(1X,'I',3X,1P,E10.3,0P,1X,'(',1P,E10.3,0P,')',1X,'<',1X,A,
 201      &1X,'<',1X,1P,E10.3,0P,16X,'I')
 202  2900 FORMAT(1X,'I',29X,A,1X,'=',1X,1P,E10.3,0P,16X,'I')
 203  3000 FORMAT(1X,'I',76X,'I'/1X,78('='))
 204  3100 FORMAT(////1X,5('*'),1X
 205      $     ,'PYSTAT_HIJING: Derived Limits on Kinematical '
 206      $     ,'Variables Used in Generation',1X,5('*'))
 207  3200 FORMAT('1',12('*'),1X
 208      $     ,'PYSTAT_HIJING: Summary of Status Codes and '
 209      $     ,'Parameter Values',1X,12('*'))
 210  3300 FORMAT(/3X,'I',4X,'MSTP(I)',9X,'PARP(I)',20X,'I',4X,'MSTP(I)',9X,
 211      &'PARP(I)'/)
 212  3400 FORMAT(1X,I3,5X,I6,6X,1P,E10.3,0P,18X,I3,5X,I6,6X,1P,E10.3)
 213
 214       RETURN
 215       END