|
JHUGen MELA
JHUGen v7.5.6, MELA v2.4.2
Matrix element calculations as used in JHUGen.
|
Go to the documentation of this file.
11 integer,
public,
parameter ::
dp = selected_real_kind(15)
12 real(8),
public,
parameter ::
tol = 0.0000001d0
13 integer,
public,
parameter ::
invalidmode=-1,
wwmode=00,
zzmode=01,
zgsmode=02,
gszmode=03,
gsgsmode=04,
zgmode=05,
gsgmode=06,
ggmode=07
46 integer,
public,
parameter ::
maxpart = 30
73 integer,
public :: collier_maxnloopprops = -1
74 integer,
public :: collier_maxrank = -1
78 real(8),
public ::
m_v,
ga_v,
m_vprime,
ga_vprime,
m_v_ps,
ga_v_ps,
m_z_ps,
ga_z_ps,
m_w_ps,
ga_w_ps
79 real(
dp),
public,
parameter :: pi =3.141592653589793238462643383279502884197_dp
80 real(
dp),
public,
parameter ::
sqrt2 = 1.4142135623730950488016887242096980786_dp
81 real(
dp),
public,
parameter ::
gamma_0 = 0.5772156649015328606065120900824024310421_dp
82 real(
dp),
public,
parameter ::
pisq = pi**2
83 real(8),
public,
parameter ::
one = 1.0d0,
mone = -1.0d0
84 real(8),
public,
parameter ::
half = 0.5d0,
two = 2.0d0
85 real(8),
public,
parameter ::
zero = 0d0
86 complex(8),
parameter,
public ::
czero = (0d0,0d0)
87 complex(8),
parameter,
public ::
cone = 1.0d0
88 complex(8),
parameter,
public ::
ci=(0d0,1.0d0)
89 complex(8),
parameter,
public ::
ne=(0d0,1.0d0)
90 real(8),
public,
parameter ::
percent=1d0/100d0
93 real(8),
public,
parameter ::
gev=1d0/100d0
94 real(8),
public,
parameter ::
tev=1d+3*
gev
95 real(8),
public,
parameter ::
mev=1d-3*
gev
96 real(8),
public,
parameter ::
kev=1d-6*
gev
97 real(8),
public,
parameter ::
ev=1d-9*
gev
98 real(8),
public,
parameter ::
meter=1d6
105 real(8),
public,
parameter ::
second=1d15
137 character(len=500),
public ::
pmzzfile =
"PMZZdistribution.out"
160 integer,
public,
parameter ::
defaultseeds(1:
nmaxseeds) = (/847362834,470115596,392845769,191039475,372910496,192049687,695820194,218930493,902943834,471748302,123958674,390534012,938576849,386472918,938576483,891928354,593857698,938576432,948576849,192847564 ,218940493,902943854,471748304,123958670,390534013,938576846,386472917,938576484,891928359,593857698,938576434,948576848,192847565/)
180 logical,
public ::
h_dk =.false.
184 character(len=2),
public ::
vh_pc =
"lo"
247 real(8),
public ::
gf = 1.16639d-5/
gev**2
253 real(8),
public ::
xw = 0.23119d0
362 complex(8),
public,
parameter ::
ahg1 = (1.0d0,0d0)
363 complex(8),
public,
parameter ::
ahg2 = (0d0,0d0)
364 complex(8),
public,
parameter ::
ahg3 = (0d0,0d0)
365 complex(8),
public,
parameter ::
ahz1 = (1.0d0,0d0)
366 complex(8),
public,
parameter ::
ahz2 = (0d0,0d0)
367 complex(8),
public,
parameter ::
ahz3 = (0d0,0d0)
371 complex(8),
public ::
ghg2 = (1.0d0,0d0)
372 complex(8),
public ::
ghg3 = (0d0,0d0)
373 complex(8),
public ::
ghg4 = (0d0,0d0)
381 complex(8),
public ::
ghz1 = (2.0d0,0d0)
382 complex(8),
public ::
ghz2 = (0d0,0d0)
383 complex(8),
public ::
ghz3 = (0d0,0d0)
384 complex(8),
public ::
ghz4 = (0d0,0d0)
454 complex(8),
public ::
ghw1 = (0d0,0d0)
455 complex(8),
public ::
ghw2 = (0d0,0d0)
456 complex(8),
public ::
ghw3 = (0d0,0d0)
457 complex(8),
public ::
ghw4 = (0d0,0d0)
721 complex(8),
public ::
gh2g2 = (0d0,0d0)
722 complex(8),
public ::
gh2g3 = (0d0,0d0)
723 complex(8),
public ::
gh2g4 = (0d0,0d0)
731 complex(8),
public ::
gh2z1 = (0.0d0,0d0)
732 complex(8),
public ::
gh2z2 = (0.0d0,0d0)
733 complex(8),
public ::
gh2z3 = (0.0d0,0d0)
734 complex(8),
public ::
gh2z4 = (0.0d0,0d0)
805 complex(8),
public ::
gh2w1 = (0.0d0,0d0)
806 complex(8),
public ::
gh2w2 = (0.0d0,0d0)
807 complex(8),
public ::
gh2w3 = (0.0d0,0d0)
808 complex(8),
public ::
gh2w4 = (0.0d0,0d0)
868 complex(8),
public ::
dv_a= (1.0d0,0d0)
869 complex(8),
public ::
dp_a= (1.0d0,0d0)
870 complex(8),
public ::
dm_a= (1.0d0,0d0)
871 complex(8),
public ::
dv_z= (1.0d0,0d0)
872 complex(8),
public ::
dp_z= (1.0d0,0d0)
873 complex(8),
public ::
dm_z= (1.0d0,0d0)
882 complex(8),
public ::
kappa = (1d0,0d0)
923 complex(8),
public ::
a1 = (0d0,0d0)
924 complex(8),
public ::
a2 = (0d0,0d0)
925 complex(8),
public ::
a3 = (0d0,0d0)
926 complex(8),
public ::
a4 = (0d0,0d0)
927 complex(8),
public ::
a5 = (0d0,0d0)
936 complex(8),
public ::
b1 = (0d0,0d0)
937 complex(8),
public ::
b2 = (0d0,0d0)
938 complex(8),
public ::
b3 = (0d0,0d0)
939 complex(8),
public ::
b4 = (0d0,0d0)
940 complex(8),
public ::
b5 = (0d0,0d0)
941 complex(8),
public ::
b6 = (0d0,0d0)
942 complex(8),
public ::
b7 = (0d0,0d0)
943 complex(8),
public ::
b8 = (0d0,0d0)
944 complex(8),
public ::
b9 = (0d0,0d0)
945 complex(8),
public ::
b10 =(0d0,0d0)
947 complex(8),
public ::
bzzp1 = (0d0,0d0)
948 complex(8),
public ::
bzzp2 = (0d0,0d0)
949 complex(8),
public ::
bzzp3 = (0d0,0d0)
950 complex(8),
public ::
bzzp4 = (0d0,0d0)
951 complex(8),
public ::
bzzp5 = (0d0,0d0)
952 complex(8),
public ::
bzzp6 = (0d0,0d0)
953 complex(8),
public ::
bzzp7 = (0d0,0d0)
954 complex(8),
public ::
bzzp8 = (0d0,0d0)
955 complex(8),
public ::
bzzp9 = (0d0,0d0)
969 complex(8),
public ::
bzgs1 = (0d0,0d0)
970 complex(8),
public ::
bzgs2 = (0d0,0d0)
971 complex(8),
public ::
bzgs3 = (0d0,0d0)
972 complex(8),
public ::
bzgs4 = (0d0,0d0)
973 complex(8),
public ::
bzgs8 = (0d0,0d0)
988 complex(8),
public,
parameter ::
c1 = (1.0d0,0d0)
989 complex(8),
public,
parameter ::
c2 = (0d0,0d0)
990 complex(8),
public,
parameter ::
c3 = (0d0,0d0)
991 complex(8),
public,
parameter ::
c41= (0d0,0d0)
992 complex(8),
public,
parameter ::
c42= (0d0,0d0)
993 complex(8),
public,
parameter ::
c5 = (0d0,0d0)
994 complex(8),
public,
parameter ::
c6 = (0d0,0d0)
995 complex(8),
public,
parameter ::
c7 = (0d0,0d0)
1013 real(
dp),
public,
parameter ::
b0_pdf(0:6) = (/ 0.8753521870054244d0,0.822300539308126d0,0.7692488916108274d0,0.716197243913529d0,0.6631455962162306d0,0.6100939485189321d0,0.5570423008216338d0 /)
1014 real(
dp),
public,
parameter ::
nf = 5.0_dp
1015 real(
dp),
public,
parameter ::
xn = 3.0_dp
1016 real(
dp),
public,
parameter ::
ca = 3.0_dp
1017 real(
dp),
public,
parameter ::
cf = 4.0_dp/3.0_dp
1018 real(
dp),
public,
parameter ::
tr = 0.5_dp
1019 real(
dp),
public,
parameter ::
avegg = 1.0_dp/4.0_dp/64.0_dp
1020 real(
dp),
public,
parameter ::
aveqg = 1.0_dp/4.0_dp/24.0_dp
1021 real(
dp),
public,
parameter ::
aveqq = 1.0_dp/4.0_dp/9.0_dp
1024 real(8),
public,
parameter ::
t3ll= -0.5d0
1025 real(8),
public,
parameter ::
t3lr= 0d0
1026 real(8),
public,
parameter ::
t3nl= 0.5d0
1027 real(8),
public,
parameter ::
t3nr= 0d0
1028 real(8),
public,
parameter ::
t3ul= 0.5d0
1029 real(8),
public,
parameter ::
t3ur= 0d0
1030 real(8),
public,
parameter ::
t3dl= -0.5d0
1031 real(8),
public,
parameter ::
t3dr= 0d0
1032 real(8),
public,
parameter ::
qll = -1d0
1033 real(8),
public,
parameter ::
qlr = -1d0
1034 real(8),
public,
parameter ::
qnl = 0d0
1035 real(8),
public,
parameter ::
qnr = 0d0
1036 real(8),
public,
parameter ::
qul = 2d0/3d0
1037 real(8),
public,
parameter ::
qur = 2d0/3d0
1038 real(8),
public,
parameter ::
qdl = -1d0/3d0
1039 real(8),
public,
parameter ::
qdr = -1d0/3d0
1084 integer,
public,
target ::
up_ = 1
1085 integer,
public,
target ::
dn_ = 2
1086 integer,
public,
target ::
chm_ = 3
1087 integer,
public,
target ::
str_ = 4
1088 integer,
public,
target ::
top_ = 5
1089 integer,
public,
target ::
bot_ = 6
1090 integer,
public,
target ::
elp_ = 7
1091 integer,
public,
target ::
mup_ = 8
1092 integer,
public,
target ::
tap_ = 9
1093 integer,
public,
target ::
glu_ = 10
1094 integer,
public,
target ::
pho_ = 11
1095 integer,
public,
target ::
z0_ = 12
1096 integer,
public,
target ::
wp_ = 13
1097 integer,
public,
target ::
nue_ = 14
1098 integer,
public,
target ::
num_ = 15
1099 integer,
public,
target ::
nut_ = 16
1100 integer,
public,
target ::
hig_ = 25
1101 integer,
public,
target ::
zpr_ = 32
1104 integer,
public,
target ::
gra_ = 39
1106 integer,
public,
target ::
aup_ = -1
1107 integer,
public,
target ::
adn_ = -2
1112 integer,
public,
target ::
elm_ = -7
1113 integer,
public,
target ::
mum_ = -8
1114 integer,
public,
target ::
tam_ = -9
1115 integer,
public,
target ::
wm_ = -13
1161 integer,
intent(in) :: index
1162 real(8),
intent(in) :: swplus, swminus, sww
1163 real(8) :: swplus_signed, swminus_signed, sww_signed, qsqcompoundfactor
1164 logical,
optional :: trywwcoupl
1166 complex(8) :: vvcoupl(1:8)
1168 real(8) :: lambda_v120(1:3)
1169 logical :: forcezzcoupl
1170 logical :: computeqsqcompundcoupl
1172 if(
present(trywwcoupl))
then
1175 forcezzcoupl = .true.
1177 computeqsqcompundcoupl = .false.
1183 if( forcezzcoupl )
then
1185 if(
cz_q2sq.ne.0) swminus_signed=abs(swminus)*dble(sign(1,
cz_q2sq))
1192 elseif(index.eq.2)
then
1196 elseif(index.eq.3)
then
1200 elseif(index.eq.4)
then
1204 elseif(index.eq.5)
then
1208 elseif(index.eq.6)
then
1212 elseif(index.eq.7)
then
1216 elseif(index.eq.8)
then
1220 elseif(index.eq.9)
then
1224 elseif(index.eq.10)
then
1228 elseif(index.eq.11)
then
1232 elseif(index.eq.12)
then
1236 elseif(index.eq.13)
then
1240 elseif(index.eq.14)
then
1244 elseif(index.eq.15)
then
1248 elseif(index.eq.16)
then
1252 elseif(index.eq.17)
then
1256 elseif(index.eq.18)
then
1260 elseif(index.eq.19)
then
1264 elseif(index.eq.20)
then
1268 elseif(index.eq.21)
then
1272 elseif(index.eq.22)
then
1276 elseif(index.eq.23)
then
1283 if(
cw_q2sq.ne.0) swminus_signed=abs(swminus)*dble(sign(1,
cw_q2sq))
1290 elseif(index.eq.2)
then
1294 elseif(index.eq.3)
then
1298 elseif(index.eq.4)
then
1302 elseif(index.eq.12)
then
1306 elseif(index.eq.13)
then
1310 elseif(index.eq.14)
then
1314 elseif(index.eq.15)
then
1318 elseif(index.eq.16)
then
1322 elseif(index.eq.17)
then
1326 elseif(index.eq.18)
then
1330 elseif(index.eq.19)
then
1346 if(computeqsqcompundcoupl)
then
1347 qsqcompoundfactor = (lambda_v120(1)**2 + swplus_signed)*(lambda_v120(2)**2 + swminus_signed)*(lambda_v120(3)**2 + sww_signed)
1348 if(qsqcompoundfactor.ne.0d0)
then
1349 qsqcompoundfactor = (lambda_v120(1)*lambda_v120(2)*lambda_v120(3))**2/qsqcompoundfactor
1356 if(computeqsqcompundcoupl)
then
1357 qsqcompoundfactor = (lambda_v120(1)**2 + swplus_signed)*(lambda_v120(2)**2 + swminus_signed)*(lambda_v120(3)**2 + sww_signed)
1358 if(qsqcompoundfactor.ne.0d0)
then
1359 qsqcompoundfactor = (lambda_v120(1)*lambda_v120(2)*lambda_v120(3))**2/qsqcompoundfactor
1370 integer,
intent(in) :: jhuid
1371 integer,
intent(in) :: hel
1372 logical,
intent(in) :: usewp
1377 if (abs(hel).ne.1)
then
1383 if (absid.eq.abs(
up_) .or. absid.eq.abs(
dn_))
then
1389 elseif (absid.eq.abs(
chm_) .or. absid.eq.abs(
str_))
then
1395 elseif (absid.eq.abs(
top_) .or. absid.eq.abs(
bot_))
then
1401 elseif (absid.eq.abs(
elp_) .or. absid.eq.abs(
nue_))
then
1407 elseif (absid.eq.abs(
mup_) .or. absid.eq.abs(
num_))
then
1413 elseif (absid.eq.abs(
tap_) .or. absid.eq.abs(
nut_))
then
1421 if (absid.eq.abs(
up_))
then
1427 elseif (absid.eq.abs(
dn_))
then
1433 elseif (absid.eq.abs(
chm_))
then
1439 elseif (absid.eq.abs(
str_))
then
1445 elseif (absid.eq.abs(
top_))
then
1451 elseif (absid.eq.abs(
bot_))
then
1457 elseif (absid.eq.abs(
elp_))
then
1463 elseif (absid.eq.abs(
nue_))
then
1469 elseif (absid.eq.abs(
mup_))
then
1475 elseif (absid.eq.abs(
num_))
then
1481 elseif (absid.eq.abs(
tap_))
then
1487 elseif (absid.eq.abs(
nut_))
then
1498 integer,
intent(in) :: pdgid
1499 integer,
intent(in) :: hel
1500 logical,
intent(in) :: usewp
1511 integer :: id1, id2, id1in, id2in
1565 integer :: id1, id2, id1in, id2in
1592 FUNCTION ckm(id1in,id2in)
1595 integer :: id1in, id2in
1600 FUNCTION zffbare(id1in,id2in,h1,h2)
1603 integer :: id1, id2, id1in, id2in
1610 print*,
"Not a valid Z to fermion pair vertex!",id1in,id2in
1614 if((id1in*h1).gt.0d0)
then
1621 if((id1in*h1).gt.0d0)
then
1628 if((id1in*h1).gt.0d0)
then
1635 if((id1in*h1).lt.0d0)
then
1638 print*,
"No right-handed neutrino here!",id1in,h1
1642 print*,
"Not a valid Z to fermion pair vertex!",id1in,id2in
1649 FUNCTION zff(id1in,id2in,h1,h2)
1652 integer :: id1in, id2in
1658 FUNCTION aff(id1in,id2in,h1,h2)
1661 integer :: id1, id2, id1in, id2in
1663 logical,
save :: printedneutrinowarning = .false.
1669 print*,
"Not a valid Z to fermion pair vertex!",id1in,id2in
1673 if((id1in*h1).gt.0d0)
then
1680 if((id1in*h1).gt.0d0)
then
1687 if((id1in*h1).gt.0d0)
then
1695 if (.not. printedneutrinowarning)
then
1696 print*,
"Warning, gamma > nu nu~ gives 0.",id1in,id2in
1697 printedneutrinowarning = .true.
1700 print*,
"Not a valid photon to fermion pair vertex!",id1in,id2in
1715 if( part.eq.0 )
then
1717 elseif( part.eq.1 )
then
1719 elseif( part.eq.2 )
then
1721 elseif( part.eq.3 )
then
1723 elseif( part.eq.4 )
then
1725 elseif( part.eq.5 )
then
1727 elseif( part.eq.6 )
then
1729 elseif( part.eq.-1 )
then
1731 elseif( part.eq.-2 )
then
1733 elseif( part.eq.-3 )
then
1735 elseif( part.eq.-4 )
then
1737 elseif( part.eq.-5 )
then
1739 elseif( part.eq.-6 )
then
1741 elseif( part.eq.21 )
then
1743 elseif( part.eq.11 )
then
1745 elseif( part.eq.22 )
then
1747 elseif( part.eq.23 )
then
1749 elseif( part.eq.24 )
then
1751 elseif( part.eq.-24 )
then
1753 elseif( part.eq.-11 )
then
1755 elseif( part.eq.13 )
then
1757 elseif( part.eq.-13 )
then
1759 elseif( part.eq.15 )
then
1761 elseif( part.eq.-15 )
then
1763 elseif( part.eq.12 )
then
1765 elseif( part.eq.-12)
then
1767 elseif( part.eq.14)
then
1769 elseif( part.eq.-14)
then
1771 elseif( part.eq.16 )
then
1773 elseif( part.eq.-16)
then
1775 elseif( part.eq.+25)
then
1777 elseif( part.eq.-25)
then
1779 elseif( part.eq.32)
then
1781 elseif( part.eq.33)
then
1783 elseif( part.eq.34)
then
1785 elseif( part.eq.-34)
then
1788 print *,
"MYLHE format not implemented for ",part
1803 if( part.eq.0 )
then
1805 elseif( part.eq.
glu_ )
then
1807 elseif( part.eq.
elm_ )
then
1809 elseif( part.eq.
elp_ )
then
1811 elseif( part.eq.
mum_ )
then
1813 elseif( part.eq.
mup_ )
then
1815 elseif( part.eq.
tam_ )
then
1817 elseif( part.eq.
tap_ )
then
1819 elseif( part.eq.
nue_ )
then
1821 elseif( part.eq.
anue_)
then
1823 elseif( part.eq.
num_)
then
1825 elseif( part.eq.
anum_)
then
1827 elseif( part.eq.
nut_ )
then
1829 elseif( part.eq.
anut_)
then
1831 elseif( part.eq.
up_ )
then
1833 elseif( part.eq.
aup_ )
then
1835 elseif( part.eq.
dn_ )
then
1837 elseif( part.eq.
adn_ )
then
1839 elseif( part.eq.
chm_ )
then
1841 elseif( part.eq.
achm_)
then
1843 elseif( part.eq.
str_ )
then
1845 elseif( part.eq.
astr_)
then
1847 elseif( part.eq.
bot_ )
then
1849 elseif( part.eq.
abot_)
then
1851 elseif( part.eq.
top_ )
then
1853 elseif( part.eq.
atop_)
then
1855 elseif( part.eq.
z0_)
then
1857 elseif( part.eq.
wp_)
then
1859 elseif( part.eq.
wm_)
then
1861 elseif( part.eq.
pho_)
then
1863 elseif( part.eq.
hig_)
then
1865 elseif( part.eq.
zpr_)
then
1867 elseif( part.eq.
zpr2_)
then
1869 elseif( part.eq.
wppr_)
then
1871 elseif( part.eq.
wmpr_)
then
1873 elseif( part.eq.
gra_)
then
1878 print *,
"LHE format not implemented for ",part
1891 if( part.eq.
glu_ )
then
1893 elseif( part.eq.
up_ )
then
1895 elseif( part.eq.
aup_ )
then
1897 elseif( part.eq.
dn_ )
then
1899 elseif( part.eq.
adn_ )
then
1901 elseif( part.eq.
chm_ )
then
1903 elseif( part.eq.
achm_)
then
1905 elseif( part.eq.
str_ )
then
1907 elseif( part.eq.
astr_)
then
1909 elseif( part.eq.
bot_ )
then
1911 elseif( part.eq.
abot_)
then
1914 print *,
"Unsuccessful conversion to a parton ME array index from ",part
1929 elseif( part.eq.
pdfup_ )
then
1931 elseif( part.eq.
pdfaup_ )
then
1933 elseif( part.eq.
pdfdn_ )
then
1935 elseif( part.eq.
pdfadn_ )
then
1937 elseif( part.eq.
pdfchm_ )
then
1941 elseif( part.eq.
pdfstr_ )
then
1945 elseif( part.eq.
pdfbot_ )
then
1949 elseif( part.eq.
pdftop_ )
then
1954 print *,
"Unsuccessful conversion to a parton id from the ME array index ",part
1961 subroutine setmass(mass, ipart)
1963 real(8),
intent(in) :: mass
1964 integer,
intent(in) :: ipart
1968 if( part.eq.abs(
elm_) )
then
1970 elseif( part.eq.abs(
mum_) )
then
1972 elseif( part.eq.abs(
tam_) )
then
1974 elseif( part.eq.abs(
chm_) )
then
1976 elseif( part.eq.abs(
bot_) )
then
1978 elseif( part.eq.abs(
top_) )
then
1980 elseif( part.eq.abs(
z0_) )
then
1982 elseif( part.eq.abs(
zpr_) )
then
1984 elseif( part.eq.abs(
wp_) )
then
1986 elseif( part.eq.abs(
wppr_) )
then
1988 elseif( part.eq.abs(
hig_) )
then
1997 real(8),
intent(in) :: width
1998 integer,
intent(in) :: ipart
2002 if( part.eq.abs(
tam_) )
then
2004 elseif( part.eq.abs(
top_) )
then
2006 elseif( part.eq.abs(
z0_) )
then
2008 elseif( part.eq.abs(
zpr_) )
then
2010 elseif( part.eq.abs(
wp_) )
then
2012 elseif( part.eq.abs(
wppr_) )
then
2014 elseif( part.eq.abs(
hig_) )
then
2026 if( part.eq.
glu_ )
then
2028 elseif( abs(part).eq.abs(
elm_) )
then
2030 elseif( abs(part).eq.abs(
mum_) )
then
2032 elseif( abs(part).eq.abs(
tam_) )
then
2034 elseif( abs(part).eq.abs(
nue_) )
then
2036 elseif( abs(part).eq.abs(
num_) )
then
2038 elseif( abs(part).eq.abs(
nut_) )
then
2040 elseif( abs(part).eq.abs(
up_) )
then
2042 elseif( abs(part).eq.abs(
dn_) )
then
2044 elseif( abs(part).eq.abs(
chm_) )
then
2046 elseif( abs(part).eq.abs(
str_) )
then
2048 elseif( abs(part).eq.abs(
bot_) )
then
2050 elseif( abs(part).eq.abs(
top_) )
then
2052 elseif( abs(part).eq.abs(
z0_) )
then
2054 elseif( abs(part).eq.abs(
zpr_) )
then
2056 elseif( abs(part).eq.abs(
wp_) )
then
2058 elseif( abs(part).eq.abs(
wppr_) )
then
2060 elseif( abs(part).eq.abs(
pho_) )
then
2062 elseif( abs(part).eq.abs(
hig_) )
then
2067 print *,
"Error in getMass",part
2079 if( abs(part).eq.abs(
top_) )
then
2081 elseif( abs(part).eq.abs(
z0_) )
then
2083 elseif( abs(part).eq.abs(
zpr_) )
then
2085 elseif( abs(part).eq.abs(
wp_) )
then
2087 elseif( abs(part).eq.abs(
wppr_) )
then
2089 elseif( abs(part).eq.abs(
hig_) )
then
2091 elseif( abs(part).eq.abs(
tam_) )
then
2104 elseif (jh .eq.2)
then
2115 elseif (jh .eq.2)
then
2127 elseif (jh .eq.2)
then
2139 elseif (jh .eq.2)
then
2151 if( part.eq.
glu_ )
then
2153 elseif( part.eq.0 )
then
2155 elseif( part.eq.
elm_ )
then
2157 elseif( part.eq.
elp_ )
then
2159 elseif( part.eq.
mum_ )
then
2161 elseif( part.eq.
mup_ )
then
2163 elseif( part.eq.
tam_ )
then
2165 elseif( part.eq.
tap_ )
then
2167 elseif( part.eq.
nue_ )
then
2169 elseif( part.eq.
anue_ )
then
2171 elseif( part.eq.
num_ )
then
2173 elseif( part.eq.
anum_ )
then
2175 elseif( part.eq.
nut_ )
then
2177 elseif( part.eq.
anut_ )
then
2179 elseif( part.eq.
up_ )
then
2181 elseif( part.eq.
aup_ )
then
2183 elseif( part.eq.
dn_ )
then
2185 elseif( part.eq.
adn_ )
then
2187 elseif( part.eq.
chm_ )
then
2189 elseif( part.eq.
achm_ )
then
2191 elseif( part.eq.
str_ )
then
2193 elseif( part.eq.
astr_ )
then
2195 elseif( part.eq.
bot_ )
then
2197 elseif( part.eq.
abot_ )
then
2199 elseif( part.eq.
top_ )
then
2201 elseif( part.eq.
atop_ )
then
2203 elseif( part.eq.
z0_ )
then
2205 elseif( part.eq.
zpr_ )
then
2207 elseif( part.eq.
wp_ )
then
2209 elseif( part.eq.
wm_ )
then
2211 elseif( part.eq.
wppr_ )
then
2213 elseif( part.eq.
wmpr_ )
then
2215 elseif( part.eq.
pho_ )
then
2217 elseif( part.eq.
hig_ )
then
2220 print *,
"Error in getParticle",part
2231 integer :: partlhe,part
2248 if( dkmode.eq.0 )
then
2250 elseif( dkmode.eq.1 )
then
2252 elseif( dkmode.eq.2 )
then
2254 elseif( dkmode.eq.3 )
then
2256 elseif( dkmode.eq.8 )
then
2258 elseif( dkmode.eq.9 )
then
2260 elseif( dkmode.lt.0 )
then
2262 dkmode.eq.-2*2 .or. &
2263 dkmode.eq.-3*3 .or. &
2264 dkmode.eq.-5*5 .or. &
2265 dkmode.eq.-7*7 .or. &
2266 dkmode.eq.-11*11 .or. &
2267 dkmode.eq.-13*13 .or. &
2283 if( dkmode.eq.4 )
then
2285 elseif( dkmode.eq.5 )
then
2287 elseif( dkmode.eq.6 )
then
2289 elseif( dkmode.eq.10 )
then
2291 elseif( dkmode.eq.11 )
then
2293 elseif( dkmode.lt.0 )
then
2295 dkmode.eq.-2*1 .or. &
2296 dkmode.eq.-3*1 .or. &
2297 dkmode.eq.-5*7 .or. &
2298 dkmode.eq.-5*13 .or. &
2299 dkmode.eq.-11*7 .or. &
2300 dkmode.eq.-11*13 .or. &
2301 dkmode.eq.-17*7 .or. &
2318 if( dkmode.eq.7 )
then
2339 islhedowntypequark = ( abs(parttype).eq.1 .or. abs(parttype).eq.3 .or. abs(parttype).eq.5 )
2352 islheuptypequark = ( abs(parttype).eq.2 .or. abs(parttype).eq.4 .or. abs(parttype).eq.6 )
2405 isalheneutrino = ( abs(parttype).eq.12 .or. abs(parttype).eq.14 .or. abs(parttype).eq.16 )
2412 isalepton = ( abs(parttype).eq.abs(
elp_) .or. abs(parttype).eq.abs(
mup_) .or. abs(parttype).eq.abs(
tap_) )
2418 isalhelepton = ( abs(parttype).eq.11 .or. abs(parttype).eq.13 .or. abs(parttype).eq.15 )
2437 FUNCTION isajet(PartType)
2456 isaboson = ( abs(parttype).eq.abs(
pho_) .or. abs(parttype).eq.abs(
z0_) .or. abs(parttype).eq.abs(
wp_) .or. abs(parttype).eq.abs(
hig_) )
2463 isalheboson = ( abs(parttype).ge.22 .and. abs(parttype).le.25 )
2468 integer,
optional :: useahcoupl
2469 integer,
intent(in) :: id(1:2),hel
2470 integer :: testahcoupl
2474 if(
present(useahcoupl))
then
2475 testahcoupl = useahcoupl
2478 (id(1).eq.
elp_ .and. id(2).eq.
nue_) .or. (id(2).eq.
elp_ .and. id(1).eq.
nue_) .or. &
2479 (id(1).eq.
mup_ .and. id(2).eq.
num_) .or. (id(2).eq.
mup_ .and. id(1).eq.
num_) .or. &
2480 (id(1).eq.
tap_ .and. id(2).eq.
nut_) .or. (id(2).eq.
tap_ .and. id(1).eq.
nut_) .or. &
2481 (id(1).eq.
up_ .and. (id(2).eq.
adn_ .or. id(2).eq.
astr_ .or. id(2).eq.
abot_)) .or. (id(2).eq.
up_ .and. (id(1).eq.
adn_ .or. id(1).eq.
astr_ .or. id(1).eq.
abot_)) .or. &
2484 ) .and. hel.lt.0)
then
2487 (id(1).eq.
elm_ .and. id(2).eq.
anue_) .or. (id(2).eq.
elm_ .and. id(1).eq.
anue_) .or. &
2488 (id(1).eq.
mum_ .and. id(2).eq.
anum_) .or. (id(2).eq.
mum_ .and. id(1).eq.
anum_) .or. &
2489 (id(1).eq.
tam_ .and. id(2).eq.
anut_) .or. (id(2).eq.
tam_ .and. id(1).eq.
anut_) .or. &
2490 (id(1).eq.
aup_ .and. (id(2).eq.
dn_ .or. id(2).eq.
str_ .or. id(2).eq.
bot_)) .or. (id(2).eq.
aup_ .and. (id(1).eq.
dn_ .or. id(1).eq.
str_ .or. id(1).eq.
bot_)) .or. &
2491 (id(1).eq.
achm_ .and. (id(2).eq.
dn_ .or. id(2).eq.
str_ .or. id(2).eq.
bot_)) .or. (id(2).eq.
achm_ .and. (id(1).eq.
dn_ .or. id(1).eq.
str_ .or. id(1).eq.
bot_)) .or. &
2492 (id(1).eq.
atop_ .and. (id(2).eq.
dn_ .or. id(2).eq.
str_ .or. id(2).eq.
bot_)) .or. (id(2).eq.
atop_ .and. (id(1).eq.
dn_ .or. id(1).eq.
str_ .or. id(1).eq.
bot_)) &
2493 ) .and. hel.lt.0)
then
2496 (id(1).eq.
elm_ .and. id(2).eq.
elp_) .or. (id(2).eq.
elm_ .and. id(1).eq.
elp_) .or. &
2497 (id(1).eq.
mum_ .and. id(2).eq.
mup_) .or. (id(2).eq.
mum_ .and. id(1).eq.
mup_) .or. &
2498 (id(1).eq.
tam_ .and. id(2).eq.
tap_) .or. (id(2).eq.
tam_ .and. id(1).eq.
tap_) .or. &
2499 (id(1).eq.
up_ .and. id(2).eq.
aup_) .or. (id(2).eq.
up_ .and. id(1).eq.
aup_) .or. &
2500 (id(1).eq.
dn_ .and. id(2).eq.
adn_) .or. (id(2).eq.
dn_ .and. id(1).eq.
adn_) .or. &
2501 (id(1).eq.
chm_ .and. id(2).eq.
achm_) .or. (id(2).eq.
chm_ .and. id(1).eq.
achm_) .or. &
2502 (id(1).eq.
str_ .and. id(2).eq.
astr_) .or. (id(2).eq.
str_ .and. id(1).eq.
astr_) .or. &
2503 (id(1).eq.
top_ .and. id(2).eq.
atop_) .or. (id(2).eq.
top_ .and. id(1).eq.
atop_) .or. &
2505 ) .and. hel.ne.0)
then
2506 if(testahcoupl.eq.1)
then
2508 elseif(testahcoupl.eq.2)
then
2514 (id(1).eq.
nue_ .and. id(2).eq.
anue_) .or. (id(2).eq.
nue_ .and. id(1).eq.
anue_) .or. &
2515 (id(1).eq.
num_ .and. id(2).eq.
anum_) .or. (id(2).eq.
num_ .and. id(1).eq.
anum_) .or. &
2517 ) .and. hel.lt.0)
then
2528 integer,
intent(in) :: id(1:2),hel
2532 (id(1).eq.
elp_ .and. id(2).eq.
nue_) .or. (id(2).eq.
elp_ .and. id(1).eq.
nue_) .or. &
2533 (id(1).eq.
mup_ .and. id(2).eq.
num_) .or. (id(2).eq.
mup_ .and. id(1).eq.
num_) .or. &
2534 (id(1).eq.
tap_ .and. id(2).eq.
nut_) .or. (id(2).eq.
tap_ .and. id(1).eq.
nut_) .or. &
2535 (id(1).eq.
up_ .and. id(2).eq.
adn_) .or. (id(2).eq.
up_ .and. id(1).eq.
adn_) .or. &
2536 (id(1).eq.
chm_ .and. id(2).eq.
astr_) .or. (id(2).eq.
chm_ .and. id(1).eq.
astr_) .or. &
2537 (id(1).eq.
top_ .and. id(2).eq.
abot_) .or. (id(2).eq.
top_ .and. id(1).eq.
abot_) .or. &
2538 (id(1).eq.
elm_ .and. id(2).eq.
anue_) .or. (id(2).eq.
elm_ .and. id(1).eq.
anue_) .or. &
2539 (id(1).eq.
mum_ .and. id(2).eq.
anum_) .or. (id(2).eq.
mum_ .and. id(1).eq.
anum_) .or. &
2540 (id(1).eq.
tam_ .and. id(2).eq.
anut_) .or. (id(2).eq.
tam_ .and. id(1).eq.
anut_) .or. &
2541 (id(1).eq.
aup_ .and. id(2).eq.
dn_ ) .or. (id(2).eq.
aup_ .and. id(1).eq.
dn_) .or. &
2542 (id(1).eq.
achm_ .and. id(2).eq.
str_) .or. (id(2).eq.
achm_ .and. id(1).eq.
str_) .or. &
2543 (id(1).eq.
atop_ .and. id(2).eq.
bot_) .or. (id(2).eq.
atop_ .and. id(1).eq.
bot_) .or. &
2544 (id(1).eq.
nue_ .and. id(2).eq.
anue_) .or. (id(2).eq.
nue_ .and. id(1).eq.
anue_) .or. &
2545 (id(1).eq.
num_ .and. id(2).eq.
anum_) .or. (id(2).eq.
num_ .and. id(1).eq.
anum_) .or. &
2547 ) .and. hel.lt.0)
then
2550 (id(1).eq.
elm_ .and. id(2).eq.
elp_) .or. (id(2).eq.
elm_ .and. id(1).eq.
elp_) .or. &
2551 (id(1).eq.
mum_ .and. id(2).eq.
mup_) .or. (id(2).eq.
mum_ .and. id(1).eq.
mup_) .or. &
2552 (id(1).eq.
tam_ .and. id(2).eq.
tap_) .or. (id(2).eq.
tam_ .and. id(1).eq.
tap_) .or. &
2553 (id(1).eq.
up_ .and. id(2).eq.
aup_) .or. (id(2).eq.
up_ .and. id(1).eq.
aup_) .or. &
2554 (id(1).eq.
dn_ .and. id(2).eq.
adn_) .or. (id(2).eq.
dn_ .and. id(1).eq.
adn_) .or. &
2555 (id(1).eq.
chm_ .and. id(2).eq.
achm_) .or. (id(2).eq.
chm_ .and. id(1).eq.
achm_) .or. &
2556 (id(1).eq.
str_ .and. id(2).eq.
astr_) .or. (id(2).eq.
str_ .and. id(1).eq.
astr_) .or. &
2557 (id(1).eq.
top_ .and. id(2).eq.
atop_) .or. (id(2).eq.
top_ .and. id(1).eq.
atop_) .or. &
2559 ) .and. hel.ne.0)
then
2577 do i = 1,
size(my_idup)
2592 do i = 1,
size(my_idup)
2607 if( abs(part).eq.
up_ )
then
2609 elseif( abs(part).eq.
dn_ )
then
2611 elseif( abs(part).eq.
chm_ )
then
2613 elseif( abs(part).eq.
str_ )
then
2615 elseif( abs(part).eq.
bot_ )
then
2617 elseif( abs(part).eq.
top_ )
then
2619 elseif( abs(part).eq.
elp_ )
then
2621 elseif( abs(part).eq.
mup_ )
then
2623 elseif( abs(part).eq.
tap_ )
then
2625 elseif( abs(part).eq.
nue_ )
then
2627 elseif( abs(part).eq.
num_ )
then
2629 elseif( abs(part).eq.
nut_ )
then
2632 print *,
"Error: Invalid flavor in SU2flip ",part
2642 integer,
intent(in) :: part(1:2)
2645 ((part(1).eq.2 .and.part(2).eq.-1) &
2646 .or.(part(1).eq.2 .and.part(2).eq.-3) &
2647 .or.(part(1).eq.2 .and.part(2).eq.-5) &
2648 .or.(part(1).eq.4 .and.part(2).eq.-1) &
2649 .or.(part(1).eq.4 .and.part(2).eq.-3) &
2650 .or.(part(1).eq.4 .and.part(2).eq.-5) &
2651 .or.(part(1).eq.-1.and.part(2).eq. 2) &
2652 .or.(part(1).eq.-3.and.part(2).eq. 2) &
2653 .or.(part(1).eq.-5.and.part(2).eq. 2) &
2654 .or.(part(1).eq.-1.and.part(2).eq. 4) &
2655 .or.(part(1).eq.-3.and.part(2).eq. 4) &
2656 .or.(part(1).eq.-5.and.part(2).eq. 4))
2665 integer,
intent(in) :: part(1:2)
2668 ((part(1).eq.-2.and.part(2).eq. 1) &
2669 .or.(part(1).eq.-2.and.part(2).eq. 3) &
2670 .or.(part(1).eq.-2.and.part(2).eq. 5) &
2671 .or.(part(1).eq.-4.and.part(2).eq. 1) &
2672 .or.(part(1).eq.-4.and.part(2).eq. 3) &
2673 .or.(part(1).eq.-4.and.part(2).eq. 5) &
2674 .or.(part(1).eq. 1.and.part(2).eq.-2) &
2675 .or.(part(1).eq. 3.and.part(2).eq.-2) &
2676 .or.(part(1).eq. 5.and.part(2).eq.-2) &
2677 .or.(part(1).eq. 1.and.part(2).eq.-4) &
2678 .or.(part(1).eq. 3.and.part(2).eq.-4) &
2679 .or.(part(1).eq. 5.and.part(2).eq.-4))
2693 if( (abs(part).ge.1 .and. abs(part).le.9) .or. (abs(part).ge.13 .and. abs(part).le.16)) then
2695 elseif( (abs(part).ge.10 .and. abs(part).le.12) ) then
2697 elseif( part.eq.25 .or. part.eq.32 .or. part.eq.39 )
then
2700 print *,
"Error: Invalid flavor in ChargeFlip"
2707 subroutine computeckmelements(inVCKM_ud, inVCKM_us, inVCKM_cd, inVCKM_cs, inVCKM_ts, inVCKM_tb, inVCKM_ub, inVCKM_cb, inVCKM_td)
2709 real(8) :: inVCKM_ud
2710 real(8) :: inVCKM_us
2711 real(8) :: inVCKM_cd
2712 real(8) :: inVCKM_cs
2713 real(8) :: inVCKM_ts
2714 real(8) :: inVCKM_tb
2715 real(8),
optional :: inVCKM_ub
2716 real(8),
optional :: inVCKM_cb
2717 real(8),
optional :: inVCKM_td
2718 real(8) :: sumVsq(1:3),diffVsq
2727 if(
present(invckm_ub))
then
2731 if (diffvsq.ge.0d0)
then
2737 if(
present(invckm_cb))
then
2741 if (diffvsq.ge.0d0)
then
2747 if(
present(invckm_td))
then
2751 if (diffvsq.ge.0d0)
then
2788 vev = 1.0d0/sqrt(
gf*sqrt(2.0d0))
2816 bl = sqrt(2d0*(1d0-
xw))
2854 DOUBLE PRECISION alphasPDF
2865 INTEGER,
PARAMETER :: NF6=6
2866 INTEGER,
PARAMETER :: NF5=5
2867 INTEGER,
PARAMETER :: NF4=4
2868 INTEGER,
PARAMETER :: NF3=3
2869 INTEGER,
PARAMETER :: NF2=2
2870 INTEGER,
PARAMETER :: NF1=1
2871 INTEGER,
PARAMETER :: NF0=0
2873 IF (
mu_ren .LE. 0d0)
THEN
2874 WRITE(6,*) .le.
'ModParameters::EvalAlphaS: Mu_Ren 0, Mu_Ren (GeV) = ',(
mu_ren*
gev)
2878 WRITE(6,*)
'ModParameters::EvalAlphaS: nQflavors_pdf has to be between 0 and 6. nQflavors_pdf = ',
nqflavors_pdf
2882 WRITE(6,*)
'ModParameters::EvalAlphaS: nQflavors_pdf invalid, nQflavors_pdf = ',
nqflavors_pdf
2883 WRITE(6,*)
'ModParameters::EvalAlphaS: Check 28472c5bfee128dde458fd4929b4d3ece9519ab8'
2887 WRITE(6,*)
'ModParameters::EvalAlphaS: nloops_pdf invalid, nloops_pdf = ',
nloops_pdf
2888 WRITE(6,*)
'ModParameters::EvalAlphaS: Check 28472c5bfee128dde458fd4929b4d3ece9519ab8'
2903 real(8) function massfrun(mf,scale)
2927 real(8),
intent(in) :: mf, scale
2928 real(8) scale_temp_ren
2933 real(8) beta0, beta1,gamma0,gamma1
2955 if ( mf.gt.
m_top )
then
2961 beta0 = ( 11d0 - 2d0/3d0 *nfrun )/4d0
2975 massfrun = mf * (as/asmf)**(gamma0/beta0)
2998 complex(dp),
intent(in) :: p(4)
2999 integer,
intent(in) :: i
3000 logical,
intent(in),
optional :: outgoing
3003 real(
dp) :: p0,px,py,pz
3004 real(
dp) :: pv,ct,st,cphi,sphi
3022 st=sqrt(abs(1.0_dp-ct**2))
3034 if ( p0 > 0.0_dp)
then
3041 if (
present(outgoing))
then
3042 if (outgoing) pol = -pol
3054 integer,
intent(in) :: i
3055 complex(dp),
intent(in) :: p(4)
3056 character(len=*),
intent(in):: out
3060 if (out ==
'out')
then
3067 function pol_dk2mom(plepton,antilepton,i,outgoing)
3069 integer,
intent(in) :: i
3071 complex(dp),
intent(in) :: plepton(1:4),antilepton(1:4)
3072 logical,
intent(in),
optional :: outgoing
3073 complex(dp) ::
pol_dk2mom(4),ub(4),v(4),q(4),qsq
3076 q=plepton+antilepton
3077 qsq=q(1)**2-q(2)**2-q(3)**2-q(4)**2
3079 ub(:)=
ubar0(plepton,i)
3080 v(:)=
v0(antilepton,-i)
3083 pol_dk2mom(1)=-(ub(2)*v(4)+v(2)*ub(4)+ub(1)*v(3)+v(1)*ub(3))
3084 pol_dk2mom(2)=-(-ub(1)*v(4)+v(1)*ub(4)-ub(2)*v(3)+v(2)*ub(3))
3085 pol_dk2mom(3)=-
ci*(ub(1)*v(4)+v(1)*ub(4)-ub(2)*v(3)-v(2)*ub(3))
3086 pol_dk2mom(4)=-(ub(2)*v(4)-v(2)*ub(4)-ub(1)*v(3)+v(1)*ub(3))
3094 if (
present(outgoing))
then
3103 integer,
intent(in) :: i
3105 complex(8),
intent(in) :: p(4)
3106 logical,
intent(in),
optional :: outgoing
3109 real(8) :: p0,px,py,pz, pv,pvsq
3110 real(8) :: ct,st,cphi,sphi
3117 pv=px**2 + py**2 + pz**2
3122 if(cdabs(pv/m).lt.1d-8)
then
3131 st= dsqrt(dabs(1.0d0-ct**2))
3135 if (st .lt. 1d-15)
then
3147 if ( p0 .gt. 0.0d0)
then
3154 if (
present(outgoing))
then
3155 if (outgoing) pol = -pol
3158 if(pol.eq.-1 .or. pol.eq.1)
then
3163 else if(pol.eq.0)
then
3169 print *,
"wrong helicity setting in pol_mass"
3177 integer,
intent(in) :: i
3178 complex(dp),
intent(in) :: p(4)
3179 character(len=*),
intent(in):: out
3183 if (out ==
'out')
then
3194 complex(dp),
intent(in) :: p(4)
3195 integer,
intent(in) :: i
3196 complex(dp) ::
ubar0(4)
3197 complex(dp) :: fc, fc2
3198 real(
dp) :: p0,px,py,pz,mass
3205 mass=dsqrt(dabs(p0**2-px**2-py**2-pz**2))
3206 if( mass.lt.1d-4 ) mass=0d0
3212 if (abs(fc2).gt.
tol)
then
3218 elseif (i.eq.-1)
then
3224 stop
'ubar0: i out of range'
3232 elseif (i.eq.-1)
then
3238 stop
'ubar0: i out of range'
3264 complex(dp),
intent(in) :: p(4)
3265 integer,
intent(in) :: i
3266 complex(dp) ::
v0(4)
3267 complex(dp) :: fc2, fc
3268 real(
dp) :: p0,px,py,pz,mass
3274 mass=dsqrt(dabs(p0**2-px**2-py**2-pz**2))
3275 if( mass.lt.1d-4 ) mass=0d0
3281 if (abs(fc2).gt.
tol)
then
3287 elseif (i.eq.-1)
then
3293 stop
'v0: i out of range'
3301 elseif (i.eq.-1)
then
3307 stop
'v0: i out of range'
3337 real(dp) :: p(:,:),two
3338 integer,
parameter :: mxpart=14
3339 complex(dp):: c23(mxpart),f(mxpart),rt(mxpart),za(:,:),zb(:,:),czero,cone,ci
3348 czero=dcmplx(0d0,0d0)
3349 cone=dcmplx(1d0,0d0)
3359 if (p(j,4) .gt. 0d0)
then
3360 rt(j)=dsqrt(p(j,4)+p(j,1))
3361 c23(j)=dcmplx(p(j,3),-p(j,2))
3365 rt(j)=dsqrt(-p(j,4)-p(j,1))
3366 c23(j)=dcmplx(-p(j,3),p(j,2))
3372 s(i,j)=two*(p(i,4)*p(j,4)-p(i,1)*p(j,1)-p(i,2)*p(j,2)-p(i,3)*p(j,3))
3373 za(i,j)=f(i)*f(j)*(c23(i)*dcmplx(rt(j)/rt(i))-c23(j)*dcmplx(rt(i)/rt(j)))
3375 if (abs(s(i,j)).lt.1d-5)
then
3376 zb(i,j)=-(f(i)*f(j))**2*dconjg(za(i,j))
3378 zb(i,j)=-dcmplx(s(i,j))/za(i,j)
3400 integer,
intent(in) :: Nmax, Rmax
3401 integer :: supNmax, supRmax
3402 supnmax = max(nmax, collier_maxnloopprops)
3403 suprmax = max(rmax, collier_maxrank)
3404 if ((supnmax .gt. collier_maxnloopprops) .or. (suprmax .gt. collier_maxrank))
then
3405 call init_cll(supnmax,suprmax,
'')
3410 integer,
intent(in) :: Nmax, Rmax
real(8), public lambda2_w42
integer, public vbfoffsh_run
integer, target, public abot_
real(dp), dimension(0:6), parameter, public b0_pdf
complex(8), public gh2w4_prime
complex(8), public ghzpzp1_prime7
real(dp), parameter, public gamma_0
complex(8), public kappa_tilde_top
complex(8), public ezp_top_right
real(8), parameter, public lambda2_w3
real(dp), public alphas_mz
real(8), public lambda_w31
complex(8), public bzpgs8
integer, target, public zpr2_
complex(8), public ghzpzp4_prime4
real(8), public lambda2_z11
complex(8), public ghwwp4_prime5
complex(8), public gh2w1_prime2
integer, target, public achm_
logical, public includegammastar
complex(8), public ghz4_prime6
integer, dimension(0:10), public debugcounter
integer, dimension(1:2), public ijpartons
complex(8), public ghzzp1_prime
integer, parameter, public nloops_pdf
real(8), public daz_neu_right
integer, target, public glu_
complex(8), public bzpgs3
integer, parameter, public pdfstr_
real(8), parameter, public br_z_nn
logical, public hasleptonfilter
complex(8), public ezp_top_left
complex(8), public gh2z1_prime7
integer, target, public elm_
complex(8), public gh2zgs2
real(8), parameter, public lambda_w1
complex(8), public ghzpzp2_prime7
complex(8), public gh2w2_prime5
complex(8), public ghz2_prime
real(dp), public couplazff
integer, parameter, public krenfacscheme_mj
integer, public widthschemein
real(8), public daz_qdn_left
integer, parameter, public pdfastr_
complex(8), public ghzzp1_prime7
complex(8), public ewp_chm_right
complex(8), public ghz4_prime4
real(8) function zffbare(id1in, id2in, h1, h2)
complex(8), public ghwpwp3
complex(8), public kappa_tilde_4gen_bot
integer, parameter, public zpgsmode
complex(8), public ghzpzp1
complex(8), public ghz3_prime5
integer, parameter, public krenfacscheme_mhstar
complex(8) function hvvspinzerodynamiccoupling(index, sWplus, sWminus, sWW, tryWWcoupl)
complex(8), public ghwwp4_prime3
real(8), parameter, public brlept_z_mm
real(8), public mphotoncutoff
logical, public distinguish_hwwcouplings
integer, public lhapdfmember
integer, dimension(1:nmaxseeds), parameter, public defaultseeds
integer, parameter, public nmaxchannels
real(8), public vprimedecaylengthmm
integer, parameter, public pdfatop_
real(8), public lambda_z10
real(8), public daz_neu_left
complex(8), public bzpzp1
real(8), parameter, public tev
complex(8), public ghzzp3_prime
logical function isagluon(PartType)
complex(8), public ghw2_prime6
integer, parameter, public pdfglu_
integer function coupledvertex(id, hel, useAHcoupl)
integer, parameter, public pdftop_
complex(8), parameter, public c7
complex(8), public bzpzp6
real(8), parameter, public tev_energy
complex(8), public ghzzp4_prime2
complex(8), public gh2z4_prime
complex(8), public ghg3_4gen
integer, parameter, public zgmode
complex(8), public ghwwp2_prime7
complex(8), public ghwwp4_prime4
integer, parameter, public krenfacscheme_mj_mj
real(8), public lambda_z41
real(8), dimension(1:2), public mvh_minmax
real(8), public lambda2_z10
real(8), parameter, public lambda2_w4
complex(8), public ghwpwp1
integer, parameter io_csmaxfile
integer, dimension(0:8), public jetsinevent
complex(8), public kappa2_tilde_4gen_top
integer, public pmzzevals
logical, public doprintpmzz
complex(8), public gh2w3_prime4
logical, public reweightdecay
complex(8), public ezp_up_left
complex(8), public gh2z3_prime2
integer, parameter, public krenfacscheme_custom_scale
complex(8), public ghzgs4
character(len=500), public pmzzfile
complex(8), public zprime_zz_1
real(8), public lambda_z12
character(len=3) function getlheparticle(PartLHE)
real(dp), parameter, public aveqq
integer, dimension(1:5, 1:5), public br_counter
complex(8), public bzpgs2
complex(8), public gh2z2_prime
real(8), public mufacmultiplier
logical, public offshellreson
real(8) function ckm(id1in, id2in)
complex(8), public ghzpzp2_prime6
complex(8), public ghz3_prime2
complex(8), public dzzwpwm
complex(8), public ghw2_prime7
complex(8), public ghzzp2_prime5
real(8), public ga_wprime
complex(8), public ghwwp4_prime7
integer, parameter, public maxpart
real(8), parameter, public brhadr_z_ss
real(8), public lambda_z21
logical, public counttauasany
complex(8), public ghzpzp2
integer, dimension(1:nmaxseeds), public theseeds
complex(8), public ghzzp4_prime3
complex(8), public ghwpwp3_prime6
real(8), public customlowerscalebound
complex(8), public kappa_4gen_bot
complex(8), public bzpzp8
integer, target, public dn_
real(8), public detajetcut
real(8), public murenmultiplier
real(8), parameter, public brlept_w_mn
complex(8), parameter, public ahz2
complex(8), public ghw2_prime3
complex(8), parameter, public ahz3
character(len=500) logfile
complex(8), public ghzzp4_prime4
real(8), dimension(1:2), public m2l_minmax
complex(8), public kappa_bot
complex(8), public ghzpzp3_prime5
real(8), public lambda_z31
complex(8), public ghzpzp3_prime2
complex(dp) function, dimension(4) pol_mless2(p, i, out)
integer, public c2z_q12sq
subroutine spinoru(p, za, zb, s)
real(8), parameter, public qul
complex(8), public ghz3_prime6
subroutine initcollier(Nmax, Rmax)
real(8), parameter, public hbarc2_pbgev2
complex(8), public ghwpwp2_prime2
complex(8), public ghzgs2
integer, parameter, public pdfdn_
real(8), public scale_alpha_w_tn
real(8), parameter, public t3ll
real(8), dimension(-6:+6,-6:+6), public crosssec
integer, parameter, public pmzzsize
real(8), public lambda2_w10
complex(8), public kappa2_4gen_bot
integer, public topdecays
integer, parameter, public krenfacscheme_mjjhstar
integer, parameter, public zpzpmode
complex(8), public ghwpwp3_prime
complex(8), public gh2z2_prime6
complex(8), public gh2w1_prime7
complex(8), public ezp_el_left
complex(8), public ghwwp2_prime5
complex(8), public gh2g3_4gen
complex(8), public ghzpzp1_prime5
complex(8), parameter, public ahg1
complex(8), public gh2z2_prime2
complex(8), public ghz1_prime
real(dp), parameter, public cf
integer function convertlhereverse(Part)
complex(8), public gh2z4_prime6
real(8), dimension(-6:+6,-6:+6), public crosssecmax
complex(8), parameter, public ahg2
integer, parameter, public krenfacscheme_default
real(8), parameter, public brlept_w_en
complex(8), public gh2gsgs2
complex(8), public ghwwp2
real(8), parameter, public one
real(8), parameter, public lambda_z2
complex(8), public gh2zgs4
logical, public fix_channels_ratio
integer, parameter, public zgsmode
integer function convertfrompartindex(Part)
real(dp), parameter, public sqrt2
complex(8), public gh2z1_prime5
complex(8), public ghzpzp3_prime
complex(8), public ezp_bot_right
logical function isajet(PartType)
real(8), public scale_alpha_w_ln
real(8), parameter, public lambda2_z1
real(8), parameter, public br_w_en
complex(8), public ewp_top_right
complex(8), public bzpzp3
real(8), parameter, public two
real(dp), public couplaffsq
complex(8), public ghz1_prime7
real(8) function getmass(Part)
complex(8), public ghwpwp2_prime3
logical, public includevprime
real(8), dimension(1:pmzzsize, 1:2), public pmzzdistribution
real(8), parameter, public hbarc2xsecunit
complex(8), parameter, public c6
complex(8), public ghwpwp3_prime7
complex(8), parameter, public ci
logical, public writeweightedlhe
real(8), dimension(1:nmaxchannels), public crosssec2
real(dp), parameter, public avegg
real(8), parameter, public lambda2_z2
real(8) function gethiggsdecaywidth(jH)
real(8) function zff(id1in, id2in, h1, h2)
real(8), parameter, public lambda
complex(8), public bgsgs1
complex(8), public bgsgs8
complex(8), public ghwpwp1_prime
complex(8), public graviton_qq_left
complex(8), public ghzzp3_prime3
complex(8), public ghzpzp1_prime4
real(8), parameter, public brhadr_w_cs
complex(8), public ghzzp1_prime3
complex(8), public ghzpzp1_prime6
complex(8), public bzzp10
complex(8), public bgsgs4
integer, parameter io_logfile
complex(8), parameter, public c42
complex(8), public ghz2_prime5
complex(8), public gh2z1_prime2
complex(dp) function, dimension(4) pol_mass2(p, i, out)
integer, public ipart_sel
real(8), parameter, public lambda_w2
complex(8), public ghgsgs4
real(8), public m_top_4gen
complex(8), public ghw3_prime5
real(8), parameter, public t3dl
integer, parameter, public wwpmode
integer(8), public rejecounter
complex(8), public gh2w2_prime6
integer, target, public adn_
complex(8), public ghw4_prime4
complex(8), public gh2w4_prime2
real(8), public m_bot_4gen
integer, target, public wm_
complex(8), public kappa2_tilde_bot
real(8), public alpha_dip
integer, public facscheme
real(8), parameter, public zero
complex(8), public ghzpzp4_prime
complex(8), public ghzpzp4_prime5
real(8), parameter, public br_z_ee
integer, target, public chm_
real(8), public etalepcut
complex(8), public kappa2_top
integer, target, public up_
real(8), parameter, public lambda_z3
integer, target, public bot_
integer, parameter, public krenfacscheme_minptj
complex(8), public gh2w4_prime5
complex(8), public ghzzp1
real(8), parameter, public milimeter
integer, parameter, public krenfacscheme_mjj
complex(8), public ghwwp4_prime2
real(8), parameter, public brlept_z_tt
complex(8), public ghzzp1_prime6
real(8) function massfrun(mf, scale)
real(8), parameter, public half
real(8), parameter, public br_w_mn
complex(8), public ghzzp4
integer, parameter io_lheoutfile
character(len=500) lhapdf_data_path
real(8), public lambda_w30
character(len=500) lheprodfile
complex(8), public gh2w2_prime4
complex(8), public ghz3_prime3
real(8), public lambda_w22
real(8), parameter, public br_z_mm
complex(dp) function, dimension(4) v0(p, i)
real(8), parameter, public tol
complex(8), public ghwwp1_prime5
real(8), parameter, public kev
complex(8), public ghw3_prime6
complex(8), public ezp_chm_right
complex(8), public ghwwp1_prime
real(8), parameter, public second
complex(8), public gh2g4_4gen
complex(8), public gh2z4_prime5
logical function isaphoton(DKMode)
integer, parameter io_lheoutfile3
logical function islhedowntypequark(PartType)
integer, public br_z_ll_counter
integer, target, public gra_
real(8), public lambda_z11
complex(8), public zprime_qq_right
real(8), parameter, public br_z_ss
integer, parameter, public zzmode
complex(8), public ghzpgs3
real(dp), parameter, public nf
real(8), parameter, public br_w_ud
real(8), public lambda2_w32
complex(8), public ghz3_prime4
real(8), parameter, public lambda_w4
logical function isalhelightquark(PartType)
logical, public convertlhefile
real(8), parameter, public lambda2_w2
complex(8), public dzawpwm
complex(8), public bzpzp9
complex(8), public gh2gsgs3
complex(8), parameter, public c41
real(8) function ckmbare(id1in, id2in)
complex(8), public ghw3_prime4
complex(8), public ghwwp3_prime
integer(8), dimension(-6:+6,-6:+6), public requevents
logical function coupledvertexisdiagonal(id, hel)
integer, parameter, public krenfacscheme_mjhstar
complex(8), public ghwpwp2
complex(8), public gh2z3_prime7
complex(8), public gh2w3_prime6
complex(8), public gh2zgs1_prime2
logical function isalepton(PartType)
subroutine computeqcdvariables()
integer, parameter, public krenfacscheme_mj_mhstar
complex(8), public ghwpwp3_prime3
real(dp), parameter, public pisq
complex(8), public ghzzp2_prime6
integer, public writefailedevents
subroutine computeewvariables()
complex(8), public ghg2_4gen
real(8), parameter, public t3nr
complex(8), public ghwpwp1_prime5
real(8), parameter, public mone
real(8), parameter, public attobarn
logical function isuptypelightquark(PartType)
real(8), parameter, public brhadr_z_bb
real(8) function aff(id1in, id2in, h1, h2)
real(8), dimension(0:10), public debugvar
complex(8), public ghz1_prime4
integer, parameter, public gsgmode
logical, public unweighted
complex(8), public gh2w2_prime2
complex(8), public ghwwp3_prime5
integer, parameter, public nrenfacschemes
complex(8), public gh2z2_prime7
real(8), public scale_alpha_z_ll
real(8), public crosssectionwithweightserrorsquared
real(8), parameter, public lambda2_zgs1
complex(8), public ghz1_prime2
complex(8), public ghzzp4_prime5
complex(8), public bzpzp5
real(8), public lambda2_z12
complex(8), public ewp_mu_left
complex(8), public ghz4_prime5
logical, public offshellv1
complex(8), public ghw4_prime3
complex(8), public ezp_mu_right
complex(8), public ghzzp2_prime
complex(8), public ezp_chm_left
integer, public vegasnc2_default
character(len=3) function getparticle(Part)
integer, public pmzzmaxindex
real(8), parameter, public lambda_zgs1
complex(8), public gh2z3_prime4
character(len=500) csmaxfile
complex(8), public ghwwp3_prime2
complex(8), public ewp_ta_left
logical, public includeinterference
complex(8), public ghw1_prime6
real(dp), parameter, public ca
complex(8), public gh2z3_prime3
complex(8), public ghzzp2
integer, public pmzzminindex
complex(8), public ghwpwp2_prime
integer, dimension(1:2), public requestnleptons
logical, parameter, public generate_bis
logical, public jetsoppositeeta
complex(8), public ghzzp3_prime7
complex(8), public kappa_4gen_top
integer(8), dimension(1:nmaxchannels), public requevents2
real(8), public lambda2_z40
complex(8), public ghwwp1
complex(8), parameter, public ahz1
logical function isalhequark(PartType)
complex(8), public ghzzp1_prime4
real(8), public lambda2_w12
integer, public printpmzzintervals
real(8), public scale_alpha_w_cs
real(8), public lambda_z32
complex(8), public ghw2_prime5
integer(8), public alertcounter
real(8), parameter, public ev
logical function islheuptypelightquark(PartType)
real(8), parameter, public brlept_w_tn
real(8), public lambda2_w11
real(dp), parameter, public xn
real(8), public lambda_w10
complex(8), public gh2z2_prime5
real(8), public lambda2_w41
real(8), public crosssectionwithweights
complex(8), public ghw2_prime4
complex(8), public ghzpgs4
integer, target, public tam_
integer, public br_w_ll_counter
complex(8), public ghzpzp2_prime3
real(8), public scale_alpha_w_ud
real(dp), public couplwffsq
complex(8), parameter, public c5
logical function isaneutrino(PartType)
complex(8), public zprime_qq_left
real(8), parameter, public brlept_z_ee
logical function isalhelepton(PartType)
complex(8), public gh2z3_prime6
complex(8), public bzpgs4
real(8), parameter, public picobarn
complex(8), public ghwpwp2_prime7
complex(8), public ghw1_prime3
complex(dp) function, dimension(4) ubar0(p, i)
complex(8), public ghwpwp4_prime2
complex(8), public gh2g2_4gen
complex(8), public ghwwp2_prime2
integer, target, public aup_
integer, parameter io_lheoutfile2
real(8), public lambda2_z20
real(8), parameter, public br_w_hadr
complex(8), public gh2w3_prime5
real(8), public lambda_w11
complex(8), public ghzpzp3
complex(8), public daawpwm
real(8), public daz_qdn_right
complex(8), public kappa2_tilde_4gen_bot
complex(8), public ghzpzp4_prime3
subroutine setdefaultckm()
real(8), public ga_zprime
logical function isawdecay(DKMode)
real(8), parameter, public gev
real(8), public daz_qup_right
real(8), parameter, public percent
complex(8), public ghwpwp3_prime2
integer, target, public nue_
complex(8), public ghw1_prime4
complex(8), public ghzpzp4_prime6
real(8), public lambda_z42
integer(8), dimension(1:nmaxchannels), public accepcounter_part2
real(8), public mininputmhstar
real(8), parameter, public milisecond
real(8), public lambda2_w30
complex(8), public ewp_up_right
subroutine sethiggsmass(jH, mass)
subroutine setmass(mass, ipart)
complex(8), public ezp_str_right
integer, parameter, public zzpmode
real(8), parameter, public br_w_ll
complex(8), parameter, public czero
complex(8), public ghzzp2_prime3
real(8), parameter, public t3ul
logical function isuptypequark(PartType)
integer, parameter, public gszmode
complex(8), public gh2w4_prime6
integer, parameter, public pdfup_
complex(8), public ghzpzp2_prime2
complex(8), public ghzpzp1_prime
real(dp), public couplzffsq
complex(8), public ezp_up_right
real(8), public lambda2_z31
complex(8), public ghzzp3
real(dp), parameter, public aveqg
complex(8), parameter, public cone
complex(8), public gh2w2_prime
integer, public renscheme
real(8), public lambda_z30
real(8), parameter, public brlept_z_nn
complex(8), public bzpzp4
complex(8), public ghzgs1_prime2
complex(8), public ghwpwp2_prime4
complex(8), parameter, public c3
integer function su2flip(Part)
complex(8) function vpffcoupling_pdg(pdgid, hel, useWp)
integer, target, public elp_
complex(8), public ghzpzp4_prime2
complex(8), public kappa_tilde
integer, public vegasnc1_default
integer, parameter, public pdfaup_
complex(8), public ghwwp1_prime3
integer, parameter, public nmaxseeds
integer function converttopartindex(Part)
logical, public findcrosssectionwithweights
real(8), public daz_qup_left
character(len=100) lhapdfstring
logical, public readlhefile
complex(8), public gh2z2_prime3
complex(8), public ghzzp3_prime5
complex(8), public kappa2_tilde_top
real(8), parameter, public qdl
complex(8), public kappa_tilde_4gen_top
real(8), parameter, public ctauunit
character(len=500) datafile
complex(8), public ghw2_prime2
integer, public br_w_ud_counter
complex(8), public ghzpgs2
integer function chargeflip(Part)
complex(8), public ghz4_prime7
complex(8), public ghwpwp4_prime5
real(8), public lambda_w20
integer, public c2w_q12sq
real(8), dimension(1:2), public m4l_minmax
complex(8), public ghz3_prime7
real(8), parameter, public br_z_hadr
integer, parameter, public zpzmode
complex(8), public gh2z4_prime7
integer, target, public atop_
real(8), dimension(1:nmaxchannels), public crosssecmax2
real(8), public lambda2_w22
logical, public useunformattedread
complex(8), public bzpzp7
real(8), parameter, public br_z_tt
real(8), parameter, public qnl
real(8), public scale_alpha_z_tt
real(8), public zmass_pdf
complex(8), public ghw4_prime6
complex(8), public ghwpwp2_prime6
complex(8), public ghzzp2_prime7
integer, dimension(1:2), public requestossf
logical, public readcsmax
character(len=2), public vh_pc
real(8), parameter, public spinavg
real(8), parameter, public t3nl
integer, target, public anum_
complex(8), public ewp_chm_left
complex(8), public ghwpwp4_prime7
logical, public foundhiggswidth
real(8), parameter, public lambda_w3
complex(8), public gh2w1_prime5
complex(8), public ghzpzp3_prime7
real(8), public alpha_qed
real(8), parameter, public lambda_q
real(8), public lambda2_w40
logical function isaboson(PartType)
real(8), parameter, public ilc_energy
complex(8), public gh2w3_prime3
real(8), public lambda_w41
complex(8), public ewp_el_right
real(8), parameter, public qnr
subroutine computeckmelements(inVCKM_ud, inVCKM_us, inVCKM_cd, inVCKM_cs, inVCKM_ts, inVCKM_tb, inVCKM_ub, inVCKM_cb, inVCKM_td)
complex(8), public kappa2_bot
complex(8), public ghzpzp2_prime5
real(8), public lambda_w40
complex(8), public ghz1_prime3
complex(8), public ghwwp3
complex(8), public ghzzp4_prime7
complex(8), public gh2w4_prime7
complex(8), public ewp_up_left
integer, target, public anut_
real(8), parameter, public meter
complex(dp) function, dimension(4) pol_dk2mom(plepton, antilepton, i, outgoing)
subroutine sethiggsdecaywidth(jH, width)
integer, target, public anue_
real(8), parameter, public br_z_inv
complex(8), public gh2w1_prime3
integer, parameter, public pdfchm_
complex(8), public gh2w4_prime4
complex(8), public ezp_dn_left
real(8), parameter, public brhadr_z_dd
complex(8), public ghz2_prime6
complex(8), public kappa2_4gen_top
complex(8), public ghz4_prime
logical function coupltolhewp(Part)
complex(8), public gh2w3_prime2
real(8), parameter, public br_w_tn
real(8) function getdecaywidth(Part)
integer, parameter io_tmpfile
logical function islheuptypequark(PartType)
complex(8), public gh2w2_prime3
real(8), public lambda_z40
complex(8), public ghwwp2_prime3
complex(8), public ghg4_4gen
integer, public vegasit1_default
complex(8), public gh2zgs3
integer, public jpart_sel
complex(8), public ezp_nue_right
integer(8), public accepcounter
complex(8), public gh2z1_prime
logical, public reweightinterference
integer function countleptons(MY_IDUP)
complex(8), public ghz3_prime
complex(8), public ghwwp3_prime3
real(8), public scale_alpha_z_nn
complex(8), public gh2z3_prime5
integer, parameter, public wpwmode
real(8), parameter, public lambda2_q
complex(8), public ghwpwp4
complex(8), public dfour_a
complex(8), public ghwpwp4_prime3
complex(8), public ghwpwp4_prime6
real(8), parameter, public t3lr
real(8), public maxinputmhstar
complex(8), public ghz2_prime4
complex(8), public ezp_ta_left
complex(8), public ghzpzp4_prime7
complex(8), public gh2z4_prime4
complex(8), public ghw2_prime
real(8), parameter, public lambda_z1
integer, target, public wppr_
real(8) function scalefactor(id1in, id2in)
integer, parameter, public pdfbot_
complex(8), public ghw1_prime7
real(8), parameter, public femtobarn
complex(8), public ghzzp1_prime2
complex(8), public ghw4_prime
complex(8), public bzpzp10
complex(8), public ezp_mu_left
real(8), parameter, public lhc_energy
complex(8), public ghwwp3_prime4
complex(8), public ghzgs3
complex(8), public ghwwp4
real(8), public breitwignercutoff
real(8), public lambda_z22
complex(8), public ghwwp2_prime
integer function convertlhe(Part)
complex(8), public ezp_dn_right
complex(8), public graviton_qq_right
real(8), parameter, public brhadr_z_uu
integer, parameter, public nqflavors_pdf
complex(8), public gh2z1_prime6
complex(8), public bzpzp2
logical function isalheneutrino(PartType)
real(8), public lambda2_w31
complex(8), public ghwpwp2_prime5
complex(8), public dfour_z
real(8), public scale_alpha_z_dd
real(8), parameter, public symmfac
real(8), parameter, public qur
complex(8), public ghzpzp1_prime2
complex(8), public ezp_nue_left
real(8), public lambda_w42
complex(8), public ghzzp3_prime2
complex(8) function vpffcoupling(jhuid, hel, useWp)
integer, parameter, public krenfacscheme_maxptj
real(8), public mhstarforphasespace
real(8), parameter, public hbarc
integer, parameter, public krenfacscheme_mj_mj_mhstar
integer, parameter, public pdfachm_
complex(8), public ghwpwp1_prime6
complex(8), public gh2w2_prime7
real(8), parameter, public lambda2
logical, public foundhiggsmass
real(8), public lambda_z20
logical, parameter, public generate_as
integer, target, public nut_
complex(8) function, dimension(4) pol_mass(p, i, outgoing)
integer, parameter, public wwmode
integer, parameter, public dp
complex(8), public ghz2_prime7
complex(8), public ghzpzp2_prime4
integer(8), dimension(-6:6,-6:6), public rejecounter_part
complex(8), public ghw4_prime5
complex(8), public ghzpzp1_prime3
complex(8), public ghwwp4_prime6
complex(8), public ghgsgs2
real(8), public lambda2_z32
integer, target, public mum_
integer, target, public zpr_
real(8), parameter, public qlr
integer, parameter io_lheinfile
complex(8), public ghzpgs1_prime2
complex(8), parameter, public c1
complex(8), public ghzpzp3_prime4
complex(8), public ghw3_prime2
logical function isalightquark(PartType)
integer, target, public wp_
real(8), dimension(1:4), parameter, public mom_not_a_particle
real(8), public channels_ratio_fix
real(8), parameter, public t3dr
logical function isalhejet(PartType)
complex(8), public gh2w3_prime
real(8), public lambda2_w20
real(8), parameter, public hbarc2_fbgev2
real(8), parameter, public planckconstant
real(8), parameter, public br_w_cs
real(8), parameter, public lambda2_z3
logical function isalheboson(PartType)
integer, public lenlhapdfstring
integer, dimension(1:2), public requestos
complex(8), public ghzzp2_prime4
complex(8), public gh2w1_prime4
complex(8), public ghzzp3_prime6
integer, public br_z_dd_counter
integer(8), public evalcounter
complex(8), public ghw3_prime3
integer, target, public str_
complex(8), public ghzpzp2_prime
real(8), public lambda2_z41
logical function isaquark(PartType)
complex(8), public ghwwp1_prime6
complex(8), public gh2w3_prime7
complex(8), public ewp_top_left
logical, public randomizevvdecays
integer, public br_z_inv_counter
complex(8), public ghwwp1_prime2
real(8), parameter, public qll
complex(8), public bzpgs1
real(8), parameter, public mev
complex(8), public ghwwp2_prime6
integer, parameter io_stdout
real(8), parameter, public hbarc2
real(8), public scale_alpha_z_uu
complex(8), public ewp_mu_right
real(8), public etajetcut
complex(8), public ghwpwp3_prime4
complex(8), public gh2z1_prime3
complex(8), public ghwwp3_prime7
real(8), public lambda2_z30
complex(8), public ghw1_prime
integer, target, public astr_
real(8), parameter, public t3ur
real(8), parameter, public speedoflight
complex(8), public ghz1_prime6
integer, parameter, public zpgmode
real(8), public lambda2_z22
real(8), public lambda_w12
complex(8), public gh2z3_prime
real(8), parameter, public lambda_z4
real(8), parameter, public brhadr_z_cc
integer, parameter, public not_a_particle_
complex(8), public bgsgs2
integer, dimension(0:8), public leptinevent
integer, parameter, public ggmode
complex(8), public ghzpzp3_prime6
real(8), public lambda_w32
logical function isdowntypequark(PartType)
logical function isazdecay(DKMode)
integer, public decaymode1
complex(8), public ghwwp4_prime
logical, public offshellv2
complex(8), public ezp_str_left
complex(8), public ghwwp1_prime4
complex(8), public ghw1_prime5
integer, target, public num_
complex(8), public ewp_ta_right
integer, parameter, public wpwpmode
complex(8), public ezp_el_right
logical, parameter, public use_dynamic_mg
real(8), public customupperscalebound
real(8), parameter, public barn
complex(8), public gh2z1_prime4
complex(8), public kappa_top
complex(8), public ghwpwp1_prime2
complex(8), public zprime_zz_2
complex(8), public ghz4_prime2
complex(8), public ghzzp1_prime5
real(8), public daz_lep_left
integer, target, public hig_
integer, parameter, public gszpmode
integer, target, public tap_
complex(8), public ghwpwp4_prime
real(8), public collider_energy
complex(8), public gh2w1_prime6
complex(8), public ghzpzp3_prime3
real(dp), parameter, public tr
integer, parameter, public invalidmode
integer function countjets(MY_IDUP)
complex(dp) function, dimension(4) pol_mless(p, i, outgoing)
complex(8), public ghzpzp4
real(8), parameter, public hbar
real(8), public daz_lep_right
complex(8), public ghwpwp4_prime4
complex(8), public ghwpwp1_prime3
integer, public br_z_uu_counter
integer, public decaymode2
integer, parameter, public krenfacscheme_mjj_mhstar
complex(8), public gh2z4_prime3
integer, target, public top_
complex(8), public ghzzp2_prime2
integer, parameter, public gsgsmode
logical function coupltolhewm(Part)
real(8), parameter, public quarkcolavg
complex(8), public ghwwp2_prime4
real(8), public overallcouplvffsq
complex(8), public printpmzz
subroutine setdecaywidth(width, ipart)
integer, target, public wmpr_
complex(8), public ghz2_prime2
complex(8), public gh2gsgs4
real(8), public lambda2_z42
real(8), public ga_vprime
logical function isalhegluon(PartType)
logical, public hbbdecays
integer, public ichann_sel
complex(8), public ghw1_prime2
complex(8), public gh2z2_prime4
complex(8), public ghz1_prime5
real(8), public higgsdecaylengthmm
complex(8), public ghz2_prime3
complex(8), public ghzzp3_prime4
integer, target, public pho_
complex(8), public ghgsgs3
real(8), public lambda_w21
complex(8), public ghw3_prime
integer, target, public z0_
complex(8), public ghzzp4_prime6
real(8), parameter, public qdr
logical, public generateevents
complex(8), public ewp_el_left
real(8), parameter, public br_z_cc
complex(8), public ghwpwp1_prime4
integer, parameter, public pdfadn_
complex(8), public ghwwp3_prime6
complex(8), public ghwpwp3_prime5
complex(8), public ghw3_prime7
complex(8), public gh2w4_prime3
complex(8), public kappa_tilde_bot
complex(8), public ghw4_prime2
complex(8), parameter, public ne
integer, public vegasnc0_default
real(8), public pmzz_mreso
character(len= *), parameter jhugen_version
complex(8), parameter, public ahg3
real(8), parameter, public lambda2_z4
real(8), parameter, public br_z_ll
complex(8), public gh2w1_prime
complex(8), public bgsgs3
real(8), public vprimedecaylengthmasscutofffactor
complex(8), public ghz4_prime3
real(8), parameter, public gluoncolavg
integer, public taudecays
real(8), public lambda2_z21
logical, public filesopened
real(8), parameter, public lambda2_w1
integer, public widthscheme
integer(8), dimension(-6:6,-6:6), public accepcounter_part
integer, parameter, public pdfabot_
complex(8), public ghwpwp1_prime7
complex(8), public ghwwp1_prime7
complex(8), public ezp_ta_right
integer, parameter io_histofile
real(8), parameter, public br_z_dd
complex(8), public ghw4_prime7
real(8), parameter, public br_z_uu
complex(8), public ezp_bot_left
real(8), parameter, public br_z_bb
real(8), parameter, public brhadr_w_ud
complex(8), parameter, public c2
complex(8), public ghzzp4_prime
integer, target, public mup_
integer, dimension(1:2), public requestnjets
real(8), public lambda2_w21
integer, public numhistograms
real(8) function gethiggsmass(jH)
complex(8), public gh2z4_prime2