| e74335a4 |
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 |
git://git.uio.no / u / mrichter / AliRoot.git / blame