← Back to branch summary

~maddevelopers/mg5amcnlo/3.0.1

« back to all changes in this revision

Viewing changes to DECAY/decay_printout.f

move ./decay to ./mg5decay; resolve unit tests (n.b. __init__ does not check keys of input dictionaries, followed last revision)

Show diffs side-by-side

added added

removed removed

Lines of Context:
DECAY/decay_printout.f
 
1
      subroutine printout
 
2
      implicit none
 
3
c
 
4
c     local
 
5
c
 
6
      integer i,iformat
 
7
      character*2 ab(2)
 
8
      real*8 ene
 
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./
 
13
c
 
14
c     include
 
15
c
 
16
      include 'coupl.inc'
 
17
c
 
18
c     Common to lh_readin and printout
 
19
c
 
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,
 
25
     &                  Vud
 
26
c
 
27
c output all info
 
28
c
 
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 )
 
43
 23   format( 1x,a,f8.4)
 
44
 24   format( 1x,a,f7.3,' GeV        ',a,f7.4,' GeV ')
 
45
 
 
46
 
 
47
      write(6,*)
 
48
      write(*,*) '*****************************************************'
 
49
      write(*,*) '*          SUMMARY OF THE SM PARAMETERS             *'
 
50
      write(*,*) '*****************************************************'
 
51
      write(6,*)
 
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
 
59
      write(6,*)
 
60
      write(6,*)
 
61
      write(6,*)  'Boson masses and widths:'
 
62
      write(6,*)  '------------------------'
 
63
      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
 
67
      write(6,*)
 
68
      write(6,*)  'Fermion masses and widths:'
 
69
      write(6,*)  '--------------------------'
 
70
      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'
 
76
      write(6,*)
 
77
      if(print_couplings) then
 
78
      write(6,*)
 
79
      write(6,*) 'Boson couplings:'
 
80
      write(6,*) '----------------'
 
81
      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
 
87
      write(6,*)
 
88
      write(6,20) 'gwwhh = ', gwwhh
 
89
      write(6,20) 'gzzhh = ', gzzhh
 
90
      write(6,20) 'ghhhh = ', ghhhh
 
91
      write(6,*)
 
92
      write(6,*) 'FFV couplings:'
 
93
      write(6,*) '--------------'
 
94
      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)
 
98
      write(6,*)
 
99
      write(6,11) 'gwf(L)   =  ',gwf(1), '    gwf(R)   =  ',gwf(2)
 
100
      write(6,*)
 
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)
 
105
      write(6,*)
 
106
      write(6,*) 'FFH couplings:'
 
107
      write(6,*) '--------------'
 
108
      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)
 
113
      write(6,*)
 
114
      write(6,*) 'Strong couplings:'
 
115
      write(6,*) '-----------------'
 
116
      write(6,*)
 
117
      write(6,14) 'gg(1)    =  ',gg(1)   , '    gg(2)    =  ',gg(2)
 
118
      write(6,*)
 
119
      endif
 
120
 
 
121
      return
 
122
      end
 
123