[cvs] / MG_ME / Models / sm / couplings.f Repository:
ViewVC logotype

Diff of /MG_ME/Models/sm/couplings.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.9, Wed Aug 22 11:58:35 2007 UTC revision 1.10, Tue Dec 4 19:03:58 2007 UTC
# Line 52  Line 52 
52  c     local  c     local
53  c  c
54        double precision  v        double precision  v
55        double precision  ee, ee2, ez, ey, sw, cw, sc2        double precision  ee, ee2, ez, ey, sw, cw, sc2, sin2w, wm
56        double precision  gwne, gwud, lambda, lam4, xt, rew, rqcd        double precision  gwne, gwud, lambda, lam4, xt, rew, rqcd
57        double precision  alphas, alfa, alfaw, mfrun        double precision  alphas, alfa, alfaw, mfrun
58        external          alphas, alfa, alfaw, mfrun        external          alphas, alfa, alfaw, mfrun
59  c  c
60  c     Common to lh_readin and printout  c     Common to lh_readin and printout
61  c  c
62        double precision  alpha, sin2w, gfermi, alfas        double precision  alpha, gfermi, alfas
63        double precision  mtMS,mbMS,mcMS,mtaMS!MSbar masses        double precision  mtMS,mbMS,mcMS,mtaMS!MSbar masses
64        double precision  Vud,Vus             !CKM matrix elements        double precision  Vud,Vus             !CKM matrix elements
65        common/values/    alpha,sin2w,gfermi,alfas,        common/values/    alpha,gfermi,alfas,
66       &                  mtMS,mbMS,mcMS,mtaMS,       &                  mtMS,mbMS,mcMS,mtaMS,
67       &                  Vud       &                  Vud
68  c  c
# Line 99  Line 99 
99  c  c
100  c auxiliary local values  c auxiliary local values
101  c  c
102          wm = sqrt(zmass**2/Two+
103         $     sqrt(zmass**4/Four-Pi/Rt2*alpha/gfermi*zmass**2))
104          sin2w  = One-(wm/zmass)**2
105        cw  = sqrt( One - sin2w )        cw  = sqrt( One - sin2w )
106        ee2 = alpha * Fourpi        ee2 = alpha * Fourpi
107        sw  = sqrt( sin2w )        sw  = sqrt( sin2w )
# Line 106  Line 109 
109        ez  = ee/(sw*cw)        ez  = ee/(sw*cw)
110        ey  = ee*(sw/cw)        ey  = ee*(sw/cw)
111        sc2 = sin2w*( One - sin2w )        sc2 = sin2w*( One - sin2w )
112        v   = Two*wmass*sw/ee   ! the wmass is used to calculate v        v   = Two*wm*sw/ee   ! the wmass is used to calculate v
113        lambda = hmass**2 / (Two * v**2)        lambda = hmass**2 / (Two * v**2)
114  c  c
115  c vector boson couplings  c vector boson couplings
# Line 220  Line 223 
223  c  c
224        integer maxpara        integer maxpara
225        parameter (maxpara=100)        parameter (maxpara=100)
226          double precision  Two, Four, Rt2, Pi
227          parameter( Two = 2.0d0, Four = 4.0d0 )
228          parameter( Rt2   = 1.414213562d0 )
229          parameter( Pi = 3.14159265358979323846d0 )
230  c  c
231  c     local  c     local
232  c  c
# Line 243  Line 250 
250  c  c
251  c     Common to lh_readin and printout  c     Common to lh_readin and printout
252  c  c
253        double precision  alpha, sin2w, gfermi, alfas        double precision  alpha, gfermi, alfas
254        double precision  mtMS,mbMS,mcMS,mtaMS!MSbar masses        double precision  mtMS,mbMS,mcMS,mtaMS!MSbar masses
255        double precision  Vud,Vus             !CKM matrix elements        double precision  Vud,Vus             !CKM matrix elements
256        common/values/    alpha,sin2w,gfermi,alfas,        common/values/    alpha,gfermi,alfas,
257       &                  mtMS,mbMS,mcMS,mtaMS,       &                  mtMS,mbMS,mcMS,mtaMS,
258       &                  Vud       &                  Vud
259  c  c
# Line 328  Line 335 
335         call set_it(n,ivalue,value,name,4,bn,zmass,91.188d0)         call set_it(n,ivalue,value,name,4,bn,zmass,91.188d0)
336         call set_it(n,ivalue,value,name,6,bn,tmass,174.3d0)         call set_it(n,ivalue,value,name,6,bn,tmass,174.3d0)
337         call set_it(n,ivalue,value,name,7,bn,lmass,1.777d0)         call set_it(n,ivalue,value,name,7,bn,lmass,1.777d0)
            bn="mgsmparam"  
        call set_it(n,ivalue,value,name,1,bn,sin2w,0.2312d0)  
        call set_it(n,ivalue,value,name,2,bn,wmass,80.419d0)  
338         bn="mgyukawa"         bn="mgyukawa"
339         call set_it(n,ivalue,value,name,4,bn,mcMS,1.25d0)         call set_it(n,ivalue,value,name,4,bn,mcMS,1.25d0)
340         call set_it(n,ivalue,value,name,5,bn,mbMS,4.2d0)         call set_it(n,ivalue,value,name,5,bn,mbMS,4.2d0)
# Line 345  Line 349 
349         call set_it(n,ivalue,value,name,15,bn,lmass,lmass*1d0)         call set_it(n,ivalue,value,name,15,bn,lmass,lmass*1d0)
350         call set_it(n,ivalue,value,name,25,bn,hmass,120d0)         call set_it(n,ivalue,value,name,25,bn,hmass,120d0)
351         call set_it(n,ivalue,value,name,23,bn,zmass,zmass*1d0)         call set_it(n,ivalue,value,name,23,bn,zmass,zmass*1d0)
352         call set_it(n,ivalue,value,name,24,bn,wmass,wmass*1d0)         call set_it(n,ivalue,value,name,24,bn,wmass,sqrt(zmass**2/Two+
353         $     sqrt(zmass**4/Four-Pi/Rt2*alpha/gfermi*zmass**2)))
354         bn="decay"         bn="decay"
355         call set_it(n,ivalue,value,name,6,bn,twidth,1.5083d0)         call set_it(n,ivalue,value,name,6,bn,twidth,1.5083d0)
356         call set_it(n,ivalue,value,name,25,bn,hwidth,0.0037d0)         call set_it(n,ivalue,value,name,25,bn,hwidth,0.0037d0)

Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10

Mail to administrator
ViewVC Help
Powered by ViewVC 1.0.0