00001
00002
00003
00004
00005 XM2 = 1.402
00006 GAM2 = 0.174
00007 FORM1 = FPIKM(SQRT(S1),AMPI,AMPI)
00008 FORM1 = FORM1*WIGFOR(QQ,XM2,GAM2)
00009
00010
00011
00012
00013 XM2 = 1.270
00014 GAM2 = 0.090
00015 AMRES = AMRO
00016 GAMRES = GAMRO
00017 AM3 = AMPI
00018 AM2 = AMPI
00019 AM1 = AMK
00020 L1 = 0
00021 L2 = 1
00022 FORM1 = BWIGK1(QQ,XM2,GAM2,AMRES,GAMRES,AM3,AM2,AM1,L1,L2)
00023 FORM1 = FORM1*BWIGKST(S1,AMRES,GAMRES,AM2,AM3,L2)
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048 +DECK,BWIGK1.
00049 COMPLEX FUNCTION BWIGK1(S,M,G,XMRES0,XGAM0,XMRES1,XMRES2,XMLEFT
00050 * ,L1,L2)
00051 IMPLICIT NONE
00052 REAL
00053 * S,M,G
00054 * ,XMRES1,XMRES2,XMLEFT
00055 * ,MRES1,MRES2,MLEFT
00056 * ,GAM
00057 * ,EPS,W
00058 * ,XMRES0,XGAM0
00059 * ,MRES0,GAMRES
00060 * ,RHOK1,RHOS
00061 * ,A,B,C
00062
00063 INTEGER L1,L2,N1,N2
00064 REAL GAUSS
00065 EXTERNAL GAUSS
00066 REAL FUNC
00067 EXTERNAL FUNC
00068
00069 COMMON /K1MASS/ W,MRES0,GAM0,MRES1,MRES2,MLEFT,N1,N2
00070
00071 COMPLEX WIGNER
00072 WIGNER(A,B,C)= CMPLX(1.0,0.0)/CMPLX(A-B**2,B*C)
00073
00074 N1 = L1
00075 N2 = L2
00076 MRES0 = XMRES0
00077 GAM0 = XGAM0
00078 MRES1 = XMRES1
00079 MRES2 = XMRES2
00080 MLEFT = XMLEFT
00081
00082 EPS = 0.00001
00083 W = M
00084 RHOK1 = GAUSS(FUNC,XMRES1+XMRES2,M-XMLEFT,EPS)
00085 W = SQRT(S)
00086 RHOS = GAUSS(FUNC,XMRES1+XMRES2,SQRT(S)-XMLEFT,EPS)
00087 GAM = G*(RHOS/RHOK1)
00088 BWIGK1 = (M**2)*WIGNER(S,M,GAM)
00089
00090 RETURN
00091 END
00092
00093
00094
00095
00096
00097
00098 +DECK,FUNC.
00099 REAL FUNCTION FUNC(X)
00100 IMPLICIT NONE
00101 REAL
00102 * X,W
00103 * ,MRES0,MRES1,MRES2,MLEFT
00104 * ,GAM,GAM0
00105 * ,Q,Q0
00106 * ,PIVAL
00107 * ,BWIG
00108 INTEGER L1,L2
00109 COMMON /K1MASS/ W,MRES0,GAM0,MRES1,MRES2,MLEFT,L1,L2
00110
00111 PIVAL = ACOS(-1.0)
00112
00113 L2 = 1
00114 IF (X.GT.(MRES1+MRES2)) THEN
00115 Q0 = (MRES0**2-(MRES1+MRES2)**2)*(MRES0**2-(MRES1-MRES2)**2)
00116 Q0 = SQRT(Q0)/(2.0*MRES0)
00117
00118 Q = (X**2-(MRES1+MRES2)**2)*(X**2-(MRES1-MRES2)**2)
00119 Q = SQRT(Q)/(2.0*X)
00120
00121 GAM = GAM0*(Q/Q0)**(2*L2+1)
00122 BWIG = (X-MRES0)**2 + (GAM/2.0)**2
00123 BWIG = (GAM/(2.0*PIVAL))/BWIG
00124 ELSE
00125 BWIG = 0.0
00126 END IF
00127
00128 L1 = 0
00129 Q = (W**2-(MLEFT+X)**2)*(W**2-(MLEFT-X)**2)
00130 Q = SQRT(Q)/(2.0*W)
00131
00132 FUNC = BWIG*(Q)**(2*L1+1)/W
00133
00134 RETURN
00135 END
00136
00137
00138
00139
00140
00141
00142 +DECK,BWIGKST.
00143 COMPLEX FUNCTION BWIGKST(S,M,G,XM1,XM2,L)
00144 IMPLICIT NONE
00145 REAL
00146 * S,M,G
00147 * ,XM1,XM2
00148 * ,GAM,Q,Q0
00149 * ,A,B,C
00150 INTEGER L
00151 COMPLEX WIGNER
00152 WIGNER(A,B,C)= CMPLX(1.0,0.0)/CMPLX(A-B**2,B*C)
00153
00154 IF (S.GT.(XM1+XM2)**2.AND.M.GT.(XM1+XM2)) THEN
00155 Q0 = (M**2-(XM1+XM2)**2)*(M**2-(XM1-XM2)**2)
00156 Q0 = SQRT(Q0)/(2.0*M)
00157
00158 Q = (S-(XM1+XM2)**2)*(S-(XM1-XM2)**2)
00159 Q = SQRT(Q)/(2.0*SQRT(S))
00160
00161 GAM = G*(Q/Q0)**(2*L+1)
00162 BWIGKST = (M**2)*WIGNER(S,M,GAM)
00163 ELSE
00164 BWIGKST = CMPLX(0.0,0.0)
00165 END IF
00166
00167 RETURN
00168 END
00169
00170