← Back to branch summary

~useakat/cfmc/gg_4g_cfmc

« back to all changes in this revision

Viewing changes to hml/hml/gmh2ff.f

  • Committer: useakat at gmail
  • Date: 2012-10-01 07:45:50 UTC
  • mfrom: (4.1.13 gg_3g_cfmc_dev)
  • Revision ID: useakat@gmail.com-20121001074550-zs09eu30khm4yvyx
independent ver.

Show diffs side-by-side

added added

removed removed

Lines of Context:
hml/hml/gmh2ff.f
 
1
C* (Inputs)
 
2
C*    AMH    : (R*4) : Higss mass.
 
3
C*    FM    : (R*4) : fermion mass.
 
4
C*    ZM    : (R*4) : Z mass.
 
5
C* (Output)
 
6
C*    GM     : (R*4) : (1/1)*( Higss partical width ) for F = L.
 
7
C*                   : (1/3)*( Higss partical width ) for F = Q.
 
8
C*
 
9
 
 
10
      subroutine gmh2ff(hm,fm,zm,gm)
 
11
 
 
12
      implicit none
 
13
 
 
14
      include 'hmparm.inc'
 
15
      include 'coupl.inc'
 
16
 
 
17
      double precision hm,fm,zm,gm
 
18
      double precision gf1,fact
 
19
      double precision ef,bt
 
20
 
 
21
c--
 
22
c  Initialize constants.
 
23
c--
 
24
      gf1  = pi*alpha/sin2w/(1.d0-sin2w)/sqrt2/zmass**2
 
25
      fact = gf1/4.d0/sqrt2/pi
 
26
c--
 
27
c  Calculate width.
 
28
c--
 
29
      ef = 0.5d0*hmass
 
30
      bt = (ef-fm)*(ef+fm)
 
31
      if(bt.le.0.d0) then
 
32
         gm = 0.d0
 
33
      else
 
34
         bt   = sqrt(bt)/ef
 
35
         gm = fact*fm**2*hmass*bt**3
 
36
      endif
 
37
 
 
38
 
 
39
      return
 
40
      end