00001 ************************************************************************************
00002 * *
00003 * Pseudo Class RRes *
00004 * *
00005 ************************************************************************************
00006 *
00007 * Parameters and physical constants
00008 *
00009 * Maarten Boonekamp, sept. 2001
00010 ************************************************************************
00011 DOUBLE PRECISION alQED,pi,xsmu
00012 PARAMETER ( alQED = 1d0/137d0 )
00013 PARAMETER ( pi = 3.14159265354 )
00014 PARAMETER ( xsmu = 4d0*pi*alQED**2/3d0 )
00015
00016 ************************************************************************
00017 *
00018 * PARTICLE DATA
00019 *
00020 * Rho(770) and Omega(782) are taken from CMD-2 F_pi fit
00021 * (hep-ex/9904027), all others are from PDG'2000.
00022 * Only 1-- resonances with known widths into electron pairs are
00023 * included (others are taken out by setting zero width into electrons).
00024 *
00025 * Maarten Boonekamp, sept. 2001
00026 ************************************************************************
00027
00028 C Declarations
00029 C ============
00030 DOUBLE PRECISION xmd, xmu, xms, xmc, xmb
00031 DOUBLE PRECISION xcd, xcu, xcs, xcc, xcb
00032 DOUBLE PRECISION xmpic, xmpi0, xmkac, xmka0, xmDcd, xmDcu, xmDcs
00033 DOUBLE PRECISION xmBbu, xmBbd, xmBbs, xmBbc
00034 DOUBLE PRECISION xmrho, xwrho, rhoee, xmrh2, xwrh2, rh2ee
00035 DOUBLE PRECISION xmrh3, xwrh3, rh3ee, xmome, xwome, omeee, omepi
00036 DOUBLE PRECISION xmom2, xwom2, om2ee, xmom3, xwom3, om3ee
00037 DOUBLE PRECISION xmphi, xwphi, phiee, xmph2, xwph2, ph2ee
00038 DOUBLE PRECISION xmpsi, xwpsi, psiee, xmps2, xwps2, ps2ee
00039 DOUBLE PRECISION xmps3, xwps3, ps3ee, xmps4, xwps4, ps4ee
00040 DOUBLE PRECISION xmps5, xwps5, ps5ee, xmps6, xwps6, ps6ee
00041 DOUBLE PRECISION xmps7, xwps7, ps7ee, xmps8, xwps8, ps8ee
00042 DOUBLE PRECISION xmups, xwups, upsee, xmup2, xwup2, up2ee
00043 DOUBLE PRECISION xmup3, xwup3, up3ee, xmup4, xwup4, up4ee
00044 DOUBLE PRECISION xmup5, xwup5, up5ee, xmup6, xwup6, up6ee
00045 DOUBLE PRECISION xmz, xwz, zwee, xmw, xww
00046 common/rresparam/xmrho, xwrho, rhoee, xmrh2, xwrh2, rh2ee,
00047 & xmrh3, xwrh3, rh3ee, xmome, xwome, omeee, xmom2, xwom2, om2ee,
00048 & xmom3, xwom3, om3ee, xmphi, xwphi, phiee, xmph2, xwph2, ph2ee,
00049 & xmpsi, xwpsi, psiee, xmps2, xwps2, ps2ee, xmps3, xwps3, ps3ee,
00050 & xmps4, xwps4, ps4ee, xmps5, xwps5, ps5ee, xmps6, xwps6, ps6ee,
00051 & xmps7, xwps7, ps7ee, xmps8, xwps8, ps8ee, xmups, xwups, upsee,
00052 & xmup2, xwup2, up2ee, xmup3, xwup3, up3ee, xmup4, xwup4, up4ee,
00053 & xmup5, xwup5, up5ee, xmup6, xwup6, up6ee, xmz, xwz, zwee,
00054 & xmw, xww
00055
00056 DOUBLE PRECISION ps3tononDDb, ps3D0toDp
00057 double precision ps4toD0D0b, ps4toDpDm, ps4toDsDs
00058 double precision ps4toD0D0Star, ps4toDpDmStar
00059 double precision ps4toD0StarD0Star, ps4toDpStarDmStar
00060 double precision ps5toD0D0b, ps5toDpDm, ps5toDsDs
00061 double precision ps5toD0D0Star, ps5toDpDmStar, ps5toDsDsStar
00062 double precision ps5toD0StarD0Star, ps5toDpStarDmStar
00063
00064 common/psidecaybrat/ps3tononddb,ps3D0toDp,
00065 & ps4toD0D0b, ps4toDpDm, ps4toDsDs,
00066 & ps4toD0D0Star, ps4toDpDmStar,
00067 & ps4toD0StarD0Star, ps4toDpStarDmStar,
00068 & ps5toD0D0b, ps5toDpDm, ps5toDsDs,
00069 & ps5toD0D0Star, ps5toDpDmStar,ps5toDsDsStar,
00070 & ps5toD0StarD0Star, ps5toDpStarDmStar
00071
00072 C QUARKS
00073 C ======
00074 c ... d
00075 PARAMETER ( xmd = 0.006d0, xcd = -1d0/3d0 )
00076 c ... u
00077 PARAMETER ( xmu = 0.003d0, xcu = 2d0/3d0 )
00078 c ... s
00079 PARAMETER ( xms = 0.122d0, xcs = -1d0/3d0 )
00080 c ... c
00081 PARAMETER ( xmc = 1.250d0, xcc = 2d0/3d0 )
00082 c ... b
00083 PARAMETER ( xmb = 4.200d0, xcb = -1d0/3d0 )
00084
00085 C SCALAR MESONS
00086 C =============
00087
00088 c ... pion
00089 PARAMETER ( xmpic= 0.13957d0, xmpi0= 0.13497d0 )
00090 c ... kaon
00091 PARAMETER ( xmkac= 0.49368d0, xmka0= 0.49767d0 )
00092 c ... D mesons
00093 PARAMETER ( xmDcd= 1.8697 d0 ) ! D+/- pingrg:2009-10-10,PDG08
00094 PARAMETER ( xmDcu= 1.8648 d0 ) ! D0
00095 PARAMETER ( xmDcs= 1.9685 d0 ) ! Ds
00096 c ... B mesons
00097 PARAMETER ( xmBbu= 5.2790 d0 ) ! B+/-
00098 PARAMETER ( xmBbd= 5.2794 d0 ) ! B0
00099 PARAMETER ( xmBbs= 5.3696 d0 ) ! Bs
00100 PARAMETER ( xmBbc= 6.4 d0 ) ! Bc
00101
00102 C VECTOR MESONS
00103 C =============
00104
00105 cC ... Rho family
00106 cc ...... rho(770)
00107 c PARAMETER ( xmrho= 0.77457d0, xwrho=147.65d-3, rhoee=6.89 d-6 )
00108 cc ...... rho(1450) <-- NOT INCLUDED
00109 c PARAMETER ( xmrh2= 1.465 d0, xwrh2=310 d-3, rh2ee=0.0000d-6 )
00110 cc ...... rho(1700) <-- NOT INCLUDED
00111 c PARAMETER ( xmrh3= 1.700 d0, xwrh3=240 d-3, rh3ee=0.0000d-6 )
00112 c
00113 cC ... Omega family
00114 cc ...... omega(782)
00115 c PARAMETER ( xmome= 0.78194d0, xwome=8.41 d-3, omeee=0.60 d-6 )
00116 PARAMETER ( omepi= 0.00011d0 )
00117 cc ...... omega(1420) <-- NOT INCLUDED
00118 c PARAMETER ( xmom2= 1.419 d0, xwom2=174 d-3, om2ee=0.0000d-6 )
00119 cc ...... omega(1650) <-- NOT INCLUDED
00120 c PARAMETER ( xmom3= 1.649 d0, xwom3=220 d-3, om3ee=0.0000d-6 )
00121
00122 cC ... Phi family
00123 cc ...... phi(1020)
00124 c PARAMETER ( xmphi= 1.01942d0, xwphi=4.46 d-3, phiee=1.33 d-6 )
00125 cc ...... phi(1680) <-- NOT INCLUDED
00126 c PARAMETER ( xmph2= 1.680 d0, xwph2=150 d-3, ph2ee=0.0000d-6 )
00127 c
00128 cC ... Psi family
00129 cc ...... J/Psi(1S)
00130 c PARAMETER ( xmpsi= 3.096916d0, xwpsi=0.091 d-3, psiee=5.40 d-6 )
00131 cc ...... Psi(2S)
00132 c PARAMETER ( xmps2= 3.686093d0, xwps2=0.281 d-3, ps2ee=2.12 d-6 )
00133 cc ...... Psi(3770)
00134 cC PARAMETER ( xmps3= 3.7700 d0, xwps3=23.6 d-3, ps3ee=0.26 d-6 )
00135 c PARAMETER ( xmps3= 3.7729 d0, xwps3=26.2 d-3, ps3ee=0.24 d-6 )
00136 cc ...... Psi(4040)
00137 cc PARAMETER ( xmps4= 4.040 d0, xwps4=52 d-3, ps4ee=0.75 d-6 )
00138 c PARAMETER ( xmps4= 4.040 d0, xwps4=74 d-3, ps4ee=0.75 d-6 )
00139 cc ...... Psi(4160)
00140 c PARAMETER ( xmps5= 4.159 d0, xwps5=78 d-3, ps5ee=0.77 d-6 )
00141 cc ...... Psi(4415)
00142 cc PARAMETER ( xmps6= 4.415 d0, xwps6=43 d-3, ps6ee=0.47 d-6 )
00143 c PARAMETER ( xmps6= 4.415 d0, xwps6=77 d-3, ps6ee=0.47 d-6 )
00144 cc ...... Y(4260)
00145 c PARAMETER ( xmps7= 4.259 d0, xwps7=88 d-3, ps7ee=0.24 d-6 )
00146 c
00147 c The cross section of Y(4260) is about 1/3 of psi(4030)
00148 c
00149
00150 cC ... Upsilon family
00151 cc ...... Upsilon(1S)
00152 c PARAMETER ( xmups= 9.46030d0, xwups=0.0525d-3, upsee=1.32 d-6 )
00153 cc ...... Upsilon(2S)
00154 c PARAMETER ( xmup2=10.02326d0, xwup2=0.044 d-3, up2ee=0.52 d-6 )
00155 cc ...... Upsilon(3S) <-- NOT INCLUDED
00156 c PARAMETER ( xmup3=10.3552 d0, xwup3=0.026 d-3, up3ee=0.0000d-6 )
00157 cc ...... Upsilon(4S)
00158 c PARAMETER ( xmup4=10.5800 d0, xwup4= 14 d-3, up4ee=0.248 d-6 )
00159 cc ...... Upsilon(10860)
00160 c PARAMETER ( xmup5=10.865 d0, xwup5=110 d-3, up5ee=0.31 d-6 )
00161 cc ...... Upsilon(11020)
00162 c PARAMETER ( xmup6=11.019 d0, xwup6= 79 d-3, up6ee=0.13 d-6 )
00163
00164 C VECTOR BOSONS
00165 C =============
00166
00167 cc ... SU(2)
00168 cc ...... Z
00169 c PARAMETER ( xmz = 91.1876d0, xwz = 2.4952d0, zwee = 0.08391d0)
00170 cc ...... W
00171 c PARAMETER ( xmw = 80.43 d0, xww = 2.11 d0 )
00172
00173
00174 ************************************************************************
00175 * Some declarations
00176 ************************************************************************
00177
00178 DOUBLE PRECISION p_pi ! decaying pion momentum
00179 DOUBLE PRECISION sdepw ! s-dependent width
00180 DOUBLE COMPLEX BrWig ! vector meson Breit-Wigner
00181
00182 ************************************************************************
00183 *
00184 * Inline function definitions
00185 *
00186 * Maarten Boonekamp, sept. 2001
00187 ************************************************************************
00188 DOUBLE PRECISION beta
00189 DOUBLE PRECISION X,Y,Z,BW_sq
00190 C Variables : X == nominal width, Y == nominal mass, Z == c.o.m. energy squared
00191 c ... kinematic suppression
00192 beta(Y,Z) = DSQRT(1d0 - MIN(Z,4d0*Y**2)/Z)
00193 c ... pion momentum in 2-pi decay
00194 p_pi(Z) = 0.5d0 * DSQRT(MAX(Z,4d0*xmpic**2) - 4d0*xmpic**2)
00195 c ... s-dependent width
00196 sdepw(X,Y,Z) = X * (p_pi(Z) / p_pi(Y**2))**3 * DSQRT(Y**2 / Z)
00197 c ... Breit-Wigner
00198 BrWig(X,Y,Z) = DCMPLX(Y**2,0d0)/DCMPLX(Y**2 - Z,Y * sdepw(X,Y,Z))
00199 BW_sq(X,Y,Z) = X / ((DSQRT(Z)-Y)**2 + X**2/4d0)
00200
00201 ************************************************************************************
00202 * *
00203 * End of Pseudo Class RRes *
00204 * *
00205 ************************************************************************************
