9
double precision Zero, One, Two, Three, Four, Half, Rt2
10
parameter( Zero = 0.0d0, One = 1.0d0, Two = 2.0d0 )
11
logical print_couplings
12
data print_couplings/.false./
18
c Common to lh_readin and printout
20
double precision alpha, sin2w, gfermi, alfas
21
double precision mtMS,mbMS,mcMS,mtaMS!MSbar masses
22
double precision Vud,Vus !CKM matrix elements
23
common/values/ alpha,sin2w,gfermi,alfas,
24
& mtMS,mbMS,mcMS,mtaMS,
29
9 format( 1x,a,f7.3,' GeV ',a,e12.4,' GeV' )
30
10 format( 1x,a,f7.3,' GeV ',a,f7.4,' GeV' )
31
11 format( 1x,a,f10.7,2x,f10.7,a,f10.7,2x,f10.7 )
32
12 format( 1x,a,f6.2,a )
33
13 format( 1x,a,f6.4,a )
34
14 format( 1x,2(a,f10.7,', ',f10.7) )
35
15 format( 1x,a,f9.5,a )
36
16 format( 1x,a,f7.5 )
37
17 format( 1x,a,f8.4 )
38
18 format( 1x,a,f8.4,' GeV' )
39
19 format( 1x,a,f6.4,a,f6.4 )
40
20 format( 1x,a,f11.5,1x,f11.5 )
41
21 format( 1x,a,f8.4,' GeV',1x,a,f8.4,' GeV' )
42
22 format( 1x,a,f10.8,a,f6.4 )
44
24 format( 1x,a,f7.3,' GeV ',a,f7.4,' GeV ')
48
write(*,*) '*****************************************************'
49
write(*,*) '* SUMMARY OF THE SM PARAMETERS *'
50
write(*,*) '*****************************************************'
52
write(6,*) ' EW Params '
53
write(6,*) '----------------------'
54
write(6,23) ' GF (10^-5*GeV^-2) ',gfermi*1d5
55
write(6,23) ' 1/alpha ',1d0/alpha
56
write(6,23) ' sin^2(W) ',sin2w
57
write(6,23) ' M_Z (GeV) ',zmass
58
write(6,23) ' M_W (GeV) ',wmass
61
write(6,*) 'Boson masses and widths:'
62
write(6,*) '------------------------'
64
write(6,24) 'Z mass = ',zmass, 'Z width = ',zwidth
65
write(6,24) 'W mass = ',wmass, 'W width = ',wwidth
66
write(6,24) 'H mass = ',hmass, 'H width = ',hwidth
68
write(6,*) 'Fermion masses and widths:'
69
write(6,*) '--------------------------'
71
write(6,24) 'top mass = ', tmass, 'top width = ', twidth
72
write(6,10) 'bottom mass = ', bmass, 'bottom width = ', Zero
73
write(6,10) 'charm mass = ', cmass, 'charm width = ', Zero
74
write(6,9) 'tau mass = ', lmass, 'tau width = ', lwidth
75
write(6,*) 'all other quark and lepton masses set to zero'
77
if(print_couplings) then
79
write(6,*) 'Boson couplings:'
80
write(6,*) '----------------'
82
write(6,20) 'gwwa = ', gwwa
83
write(6,20) 'gwwz = ', gwwz
84
write(6,20) 'gwwh = ', gwwh
85
write(6,20) 'gzzh = ', gzzh
86
write(6,20) 'ghhh = ', ghhh
88
write(6,20) 'gwwhh = ', gwwhh
89
write(6,20) 'gzzhh = ', gzzhh
90
write(6,20) 'ghhhh = ', ghhhh
92
write(6,*) 'FFV couplings:'
93
write(6,*) '--------------'
95
write(6,11) 'gal(L) = ',gal(1), ' gal(R) = ',gal(2)
96
write(6,11) 'gau(L) = ',gau(1), ' gau(R) = ',gau(2)
97
write(6,11) 'gad(L) = ',gad(1), ' gad(R) = ',gad(2)
99
write(6,11) 'gwf(L) = ',gwf(1), ' gwf(R) = ',gwf(2)
101
write(6,11) 'gzn(L) = ',gzn(1), ' gzn(R) = ',gzn(2)
102
write(6,11) 'gzl(L) = ',gzl(1), ' gzl(R) = ',gzl(2)
103
write(6,11) 'gzu(L) = ',gzu(1), ' gzu(R) = ',gzu(2)
104
write(6,11) 'gzd(L) = ',gzd(1), ' gzd(R) = ',gzd(2)
106
write(6,*) 'FFH couplings:'
107
write(6,*) '--------------'
109
write(6,14) 'gHtop(L) = ',ghtop(1), ' gHtop(R) = ',ghtop(2)
110
write(6,14) 'gHbot(L) = ',ghbot(1), ' gHbot(R) = ',ghbot(2)
111
write(6,14) 'gHcha(L) = ',ghcha(1), ' gHcha(R) = ',ghcha(2)
112
write(6,14) 'gHtau(L) = ',ghtau(1), ' gHtau(R) = ',ghtau(2)
114
write(6,*) 'Strong couplings:'
115
write(6,*) '-----------------'
117
write(6,14) 'gg(1) = ',gg(1) , ' gg(2) = ',gg(2)
added
removed
MG4_DECAY/decay_printout.f
