00001 REAL FUNCTION PKORB(IF1,IF2)
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015 INTEGER IF1
00016 INTEGER IF2
00017
00018
00019
00020
00021
00022
00023 INTEGER JAK1,JAK2,JAKP,JAKM,KTOM
00024 COMMON / JAKI / JAK1,JAK2,JAKP,JAKM,KTOM
00025 REAL*4 RRR(1)
00026 REAL PARM(4,100)
00027 integer imixpp(300)
00028 INTEGER INIT,I,J
00029 REAL C1270,C1402,A1270_KSPI,A1270_KRHO,A1402_KSPI,A1402_KRHO
00030 REAL CG1,CG2,R,BRA1,BRKS
00031 SAVE INIT,PARM
00032 DATA INIT/0/
00033
00034
00035
00036 PKORB = 0.
00037
00038
00039
00040 IF (INIT.EQ.0) THEN
00041 INIT = 1
00042
00043
00044 DO I=1,4
00045 DO J=1,100
00046 PARM(I,J) = 0
00047 END DO
00048 END DO
00049
00050
00051
00052 PARM(1, 1) = 1.777000
00053 PARM(1, 2) = 0.
00054 PARM(1, 3) = 0.000511
00055 PARM(1, 4) = 0.
00056 PARM(1, 5) = 0.105658
00057 PARM(1, 6) = 0.
00058 PARM(1, 7) = 0.134976
00059 PARM(1, 8) = 0.139570
00060 PARM(1, 9) = 0.769900
00061 PARM(1,10) = 1.275000
00062 PARM(1,11) = 0.493677
00063 PARM(1,12) = 0.497670
00064 PARM(1,13) = 0.891590
00065 PARM(1,14) = 0.781940
00066 PARM(1,15) = 1.370000
00067 PARM(1,16) = 1.700000
00068 PARM(1,17) = 1.461000
00069 PARM(1,18) = 1.300000
00070 PARM(1,19) = 1.270000
00071 PARM(1,20) = 1.402000
00072 PARM(1,21) = 1.465000
00073 PARM(1,22) = 1.700000
00074
00075
00076 PARM(2, 1) = 0.
00077 PARM(2, 2) = 0.
00078 PARM(2, 3) = 0.
00079 PARM(2, 4) = 0.
00080 PARM(2, 5) = 0.
00081 PARM(2, 6) = 0.
00082 PARM(2, 7) = 0.
00083 PARM(2, 8) = 0.
00084 PARM(2, 9) = 0.1512
00085 PARM(2,10) = 0.700
00086 PARM(2,11) = 0.
00087 PARM(2,12) = 0.
00088 PARM(2,13) = 0.0498
00089 PARM(2,14) = 0.00843
00090 PARM(2,15) = 0.510
00091 PARM(2,16) = 0.235
00092 PARM(2,17) = 0.250
00093 PARM(2,18) = 0.400
00094 PARM(2,19) = 0.090
00095 PARM(2,20) = 0.174
00096 PARM(2,21) = 0.310
00097 PARM(2,22) = 0.235
00098
00099
00100
00101
00102 PARM(3,15) = -0.145
00103
00104 IMIXPP(205)=1
00105 IMIXPP(207)=1
00106 IMIXPP(209)=1
00107 IMIXPP(211)=1
00108 IMIXPP(201)=1
00109 IMIXPP(203)=1
00110 IMIXPP(213)=1
00111 IMIXPP(215)=1
00112
00113
00114 IF (IMIXPP(205).NE.0) PARM(3,15) = -0.110
00115 IF (IMIXPP(207).NE.0) PARM(3,16) = -0.038
00116 IF (IMIXPP(209).NE.0) PARM(3,17) = 0.00
00117 IF (IMIXPP(211).NE.0) PARM(3,18) = 0.00
00118 IF (IMIXPP(201).NE.0) PARM(3,19) = 1.0
00119 IF (IMIXPP(203).NE.0) PARM(3,20) = 0.8
00120 IF (IMIXPP(213).NE.0) PARM(3,21) = -0.110
00121 IF (IMIXPP(215).NE.0) PARM(3,22) = -0.110
00122
00123 PRINT *,' KORB: rho/rhop -> pi-pi0 mixing:'
00124 PRINT *,' KORB: rho =',PARM(1,9) ,PARM(2,9)
00125 PRINT *,' KORB: rhop =',PARM(1,15),PARM(2,15),PARM(3,15)
00126 PRINT *,' KORB: K*/K*prime -> Kpi mixing:'
00127 PRINT *,' KORB: kstp =',PARM(1,16),PARM(2,16),PARM(3,16)
00128 PRINT *,' KORB: a1/a1prime -> 3pi, KKpi mixing:'
00129 PRINT *,' KORB: a1 =',PARM(1,10),PARM(2,10)
00130 PRINT *,' KORB: a1prim=',PARM(1,17),PARM(2,17),PARM(3,17)
00131 PRINT *,' KORB: K1A/K1B -> Kpipi mixing:'
00132 PRINT *,' KORB: K1A =',PARM(1,19),PARM(2,19),PARM(3,19)
00133 PRINT *,' KORB: K1B =',PARM(1,20),PARM(2,20),PARM(3,20)
00134 PRINT *,' KORB: rho/rhop/rhopp -> 4pi mixing:'
00135 PRINT *,' KORB: rho =',PARM(1,9) ,PARM(2,9)
00136 PRINT *,' KORB: rhopp =',PARM(1,21),PARM(2,21),PARM(3,21)
00137 PRINT *,' KORB: rhoppp=',PARM(1,22),PARM(2,22),PARM(3,22)
00138
00139
00140
00141 PARM(3,31) = 0.
00142 PARM(3,32) = 0.1242
00143 PARM(3,33) = 0.1604
00144 PARM(3,34) = 0.2711
00145 PARM(3,35) = 0.4443
00146
00147 PARM(3,36) = 0.
00148 PARM(3,37) = 1.0
00149
00150
00151
00152 IPLIST=0
00153 IF (IPLIST.EQ.5) THEN
00154 PARM(3,31) = 0.0000
00155 PARM(3,32) = 0.1242
00156 PARM(3,33) = 0.1604
00157 PARM(3,34) = 0.2711
00158 PARM(3,35) = 0.4443
00159 PARM(3,36) = 0.0000
00160 PARM(3,37) = 1.0000
00161 END IF
00162
00163 PRINT *,' KORB: 3PI-PI0 PARAMS:',(PARM(3,I),I=31,35)
00164 PRINT *,' KORB: PI-3PI0 PARAMS:',(PARM(3,I),I=36,37)
00165
00166
00167
00168
00169 IF (ABS(PARM(3,31)-0.0000).GT.0.0001 .OR.
00170 1 ABS(PARM(3,32)-0.1242).GT.0.0001 .OR.
00171 1 ABS(PARM(3,33)-0.1604).GT.0.0001 .OR.
00172 1 ABS(PARM(3,34)-0.2711).GT.0.0001 .OR.
00173 1 ABS(PARM(3,35)-0.4443).GT.0.0001 ) THEN
00174 PARM(3,38) = -1.
00175 ELSE
00176 PARM(3,38) = 6.9673671E-14
00177 END IF
00178
00179 IF (ABS(PARM(3,36)-0.0000).GT.0.0001 .OR.
00180 1 ABS(PARM(3,37)-1.0000).GT.0.0001 ) THEN
00181 PARM(3,39) = -1.
00182 ELSE
00183 PARM(3,39) = 3.5374880E-13
00184 END IF
00185
00186
00187
00188 PARM(3,42) = -0.40
00189 PARM(3,43) = 0.00
00190 PARM(3,44) = -0.20+3.1416
00191 PARM(3,45) = -1.50
00192
00193
00194 PARM(3,51) = -0.10
00195 PARM(3,52) = 1.00
00196 PARM(3,53) = -0.10
00197 PARM(3,54) = -0.04
00198
00199
00200 PARM(3,55) = 1.00
00201 PARM(3,56) = 0.14
00202 PARM(3,57) = -0.05
00203 PARM(3,58) = -0.05
00204
00205
00206 PARM(3,59) = 1.000
00207 PARM(3,60) = -0.145
00208 PARM(3,61) = 0.000
00209
00210
00211
00212 PARM(4,1) = 0.4920
00213 PARM(4,2) = 0.4920
00214 PARM(4,3) = 0.6660
00215 PARM(4,4) = 0.6660
00216 PARM(4,5) = 0.5
00217 PARM(4,6) = 0.5
00218
00219
00220 C1270 = PARM(3,19)
00221 C1402 = PARM(3,20)
00222 IF (C1270.EQ.0.AND.C1402.EQ.0.) THEN
00223 C1270 = 1.
00224 C1402 = 0.6
00225 END IF
00226
00227 A1270_KSPI = SQRT(0.16)
00228 A1270_KRHO = SQRT(0.42)
00229 A1402_KSPI = SQRT(0.94)
00230 A1402_KRHO = SQRT(0.03)
00231
00232 CG1 = -SQRT(2./3.)
00233 CG2 = SQRT(1./3.)
00234
00235 PARM(3,81) = C1270*A1270_KSPI*CG1
00236 PARM(3,82) = C1402*A1402_KSPI*CG1
00237 PARM(3,83) = C1270*A1270_KRHO*CG1
00238 PARM(3,84) = C1402*A1402_KRHO*CG1
00239 PARM(3,85) = C1270*A1270_KSPI*CG2
00240 PARM(3,86) = C1402*A1402_KSPI*CG2
00241 PARM(3,87) = C1270*A1270_KRHO*CG2
00242 PARM(3,88) = C1402*A1402_KRHO*CG2
00243
00244 END IF
00245
00246
00247 R = 0.
00248 IF (IF1.GE.1 .AND. IF1.LE.4 .AND. IF2.GE.1 .AND. IF2.LE.100) THEN
00249 R = PARM(IF1,IF2)
00250
00251
00252 IF (IF1.EQ.4.AND.JAK1.EQ.5) THEN
00253 IF (IF2.EQ.11) THEN
00254
00255 R = BRA1
00256 ELSE IF (IF2.EQ.1) THEN
00257 BRA1 = R
00258 CALL RANMAR(RRR,1)
00259 IF (RRR(1).LT.BRA1) THEN
00260 R = 1.
00261 ELSE
00262 R = 0.
00263 END IF
00264 BRA1 = R
00265 END IF
00266 ELSEIF (IF1.EQ.4.AND.JAK1.EQ.7) THEN
00267 IF (IF2.EQ.13) THEN
00268
00269 R = BRKS
00270 ELSE IF (IF2.EQ.3) THEN
00271 BRKS = R
00272 CALL RANMAR(RRR,1)
00273 IF (RRR(1).LT.BRKS) THEN
00274 R = 1.
00275 ELSE
00276 R = 0.
00277 END IF
00278 BRKS = R
00279 END IF
00280 END IF
00281
00282 END IF
00283
00284 PKORB = R
00285 RETURN
00286 END