bin/BCVEGPY/simparameter.F

0001
0002 c******************************************************

0003 c...     some frequently used parameters.

0004 c******************************************************

0005       subroutine parameters()
0006       implicit double precision(a-h, o-z)
0007         implicit integer(i-n)
0008
0009 #include "inclcon.h"
0010       common/ppp/pp(4,40),guv(4)
0011       common/rconst/pi
0012        double complex colmat, bundamp
0013         common/upcom/ecm,pmbc,pmb,pmc,fbcc,pmomup(5,8),
0014      &  colmat(10,64),bundamp(4),pmomzero(5,8)
0015
0016       hbcvm=pmbc
0017       hbcm=pmbc
0018       fmb=pmb
0019       fmc=pmc
0020
0021       wwpi=pi
0022         hbcp=fbcc
0023
0024       fb1=(3*hbcm**2)/(8*fmb*fmc)
0025       fb2=(3*(fmb-fmc)*hbcm)/(8*fmb*fmc)
0026       fb3=hbcm**2/(8*fmb*fmc)
0027       fb4=((fmb-fmc)*hbcm)/(8*fmb*fmc)
0028
0029 c...the color factor \delta_{ij}/3 of the bound state has

0030 c...been included in the following definiation for

0031 c...the p-wave matrix elements times the factor from bound

0032 c...state, <0|p|0>*c_p.

0033       cbc1p1=dsqrt(hbcp/(72.0d0*hbcm**3))
0034       cbcp0 =dsqrt(hbcp/(216.0d0*hbcm**3))
0035       cbcp1 =dsqrt(hbcp/(144.0d0*hbcm**3))
0036 c...here for simplicity one (hbcm) has been absorbed into

0037 c...matrix element.

0038       cbcp2 =dsqrt(hbcp/(72.0d0*hbcm))
0039
0040       return
0041       end
0042
0043 c************************************************

0044 c...   this is only used for p-wave generation.

0045 c...   some simplified paremters.

0046 c************************************************

0047       subroutine dparameters()
0048       implicit double precision(a-h, o-z)
0049         implicit integer(i-n)
0050
0051 #include "inclcon.h"
0052       common/ppp/pp(4,40),guv(4)
0053
0054       guv(1)= 1.0d0
0055       guv(2)=-1.0d0
0056       guv(3)=-1.0d0
0057       guv(4)=-1.0d0
0058
0059       fmc2=fmc**2
0060       fmb2=fmb**2
0061       fmb3=fmb**3
0062       fmc3=fmc**3
0063
0064       hbcvm2=hbcvm**2
0065       hbcm2 =hbcm**2
0066       hbcm3 =hbcm**3
0067       hbcm4 =hbcm**4
0068       hbcm5 =hbcm**5
0069
0070       dhbcvm2=1.0d0/hbcvm2
0071       dhbcm2 =1.0d0/hbcm2
0072
0073       ffmcfmb=fmc/(fmb+fmc)
0074
0075       return
0076       end
0077
0078 c***********************************************

0079 c******** parameters for s-wave ****************

0080 c***********************************************

0081
0082         subroutine paraswave(ibco)
0083 c...preamble: declarations.

0084       implicit double precision(a-h, o-z)
0085         implicit integer(i-n)
0086
0087         common/rconst/pi
0088         double complex colmat, bundamp
0089        common/upcom/ecm,pmbc,pmb,pmc,fbcc,pmomup(5,8),
0090      &  colmat(10,64),bundamp(4),pmomzero(5,8)
0091         common/wavezero/fbc
0092       common/coloct/ioctet
0093         common/octmatrix/coeoct
0094
0095         ioctet =0
0096         fbc =1.241d0
0097         fbcc=dsqrt(3.0d0*fbc**2/pi/pmbc)
0098
0099 c....for color-octet states.

0100         if(ibco.eq.7.or.ibco.eq.8) then
0101           ioctet =1
0102           coeoct =0.2d0  ! the value of \delta_s(v^2)
0103                          ! defined in prd71,074012(2005)
0104         end if
0105
0106         call parameters()
0107       call dparameters()
0108       call coupling()
0109
0110         return
0111         end
0112
0113 c***********************************************

0114 c******** parameters for p-wave ****************

0115 c***********************************************

0116
0117         subroutine parapwave
0118 c...preamble: declarations.

0119       implicit double precision(a-h, o-z)
0120         implicit integer(i-n)
0121
0122         common/rconst/pi
0123         double complex colmat, bundamp
0124       common/upcom/ecm,pmbc,pmb,pmc,fbcc,pmomup(5,8),
0125      &  colmat(10,64),bundamp(4),pmomzero(5,8)
0126         common/wavezero/fbc
0127       common/coloct/ioctet
0128         common/octmatrix/coeoct
0129
0130         ioctet   =0
0131         fbc =0.44833d0
0132         fbcc=fbc**2*9.0d0/(2.0d0*pi)
0133
0134         call parameters()
0135       call dparameters()
0136       call coupling()
0137
0138         return
0139         end