Page MenuHomeHEPForge
Files F10881587

parameters.2HDM.f90
No OneTemporary
Actions

parameters.2HDM.f90
View Options

! $Id: parameters.SM.f90,v 1.4 2006/06/16 13:31:48 kilian Exp $
!
! Copyright (C) 1999-2015 by
! Wolfgang Kilian <kilian@physik.uni-siegen.de>
! Thorsten Ohl <ohl@physik.uni-wuerzburg.de>
! Juergen Reuter <juergen.reuter@desy.de>
! with contributions from
! Christian Speckner <cnspeckn@googlemail.com>
! David Gordo Gomez (only this file)
!
! WHIZARD is free software; you can redistribute it and/or modify it
! under the terms of the GNU General Public License as published by
! the Free Software Foundation; either version 2, or (at your option)
! any later version.
!
! WHIZARD is distributed in the hope that it will be useful, but
! WITHOUT ANY WARRANTY; without even the implied warranty of
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
! GNU General Public License for more details.
!
! You should have received a copy of the GNU General Public License
! along with this program; if not, write to the Free Software
! Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
module parameters_2hdm
use kinds
use constants
use sm_physics !NODEP!
implicit none
private
real(default), dimension(37), public :: mass, width
real(default), public :: as
complex(default), public :: gs, igs
real(default), public :: e, g, e_em
real(default), public :: sinthw, costhw, sin2thw, tanthw
real(default), public :: qelep, qeup, qedwn
complex(default), public :: qlep, qup, qdwn, gcc, qw, &
gzww, gwww, &
ghbb, ghtt, ghcc, ghtautau, ghmm, &
iqw, igzww, igwww, gw4, gzzww, gazww, gaaww
real(default), public :: vev
complex(default), dimension(2), public :: &
gncneu, gnclep, gncup, gncdwn, &
gh1uu, gh2uu, gh3uu, gh1cc, gh2cc, &
gh3cc, gh1tt, gh2tt, gh3tt, gh1dd, &
gh2dd, gh3dd, gh1ss, gh2ss, gh3ss, &
gh1bb, gh2bb, gh3bb, ghud, ghcs, ghtb, &
ghdu, ghsc, ghbt, &
gh1e1e1, gh2e1e1, gh3e1e1, gh1e2e2, gh2e2e2, &
gh3e2e2, gh1e3e3, gh2e3e3, gh3e3e3
real(default), public :: cotanthw2, sinthw2, &
R11, R12, R13, R21, R22, R23, R31, R32, R33, &
lamb1, lamb2, lamb3, lamb4, lamb5R, lamb5I, lamb6R, &
lamb6I, lamb7R, lamb7I
complex(default), public :: gAHpHm, gZHpHm, gZh1h2, gZh1h3, gZh2h3, &
gWpHmh1, gWpHmh2, gWpHmh3, gWmHph1, gWmHph2, gWmHph3, gh1ZZ, &
gh2ZZ, gh3ZZ, gh1WpWm, gh2WpWm, gh3WpWm, ghhWpWm, ghhZZ, gHpHmAA, &
gHpHmZZ, gHpHmAZ, gHpHmWpWm, gh1HpAWm, gh2HpAWm, gh3HpAWm, &
gh1HpZWm, gh2HpZWm, gh3HpZWm, gh1HpAWmC, gh2HpAWmC, gh3HpAWmC, &
gh1HpZWmC, gh2HpZWmC, gh3HpZWmC, gh1HpHm, gh2HpHm, gh3HpHm, &
gh111, gh112, gh113, gh221, gh222, gh223, gh331, gh332, gh333, &
gh123, gHpHmHpHm, gHpHm11, gHpHm12, gHpHm13, gHpHm22, gHpHm23, &
gHpHm33, gh1111, gh1112, gh1113, gh1122, gh1123, gh1133, gh1222, &
gh1223, gh1233, gh1333, gh2222, gh2223, gh2233, gh2333, gh3333, &
lamb5,lamb6,lamb7, &
zetau, zetad, zetal, &
ghe1n1, ghe2n2, ghe3n3, ghn1e1, ghn2e2, ghn3e3
public :: import_from_whizard, model_update_alpha_s
contains
subroutine import_from_whizard (par_array)
real(default), dimension(57), intent(in) :: par_array
type :: parameter_set
real(default) :: gf
real(default) :: mZ
real(default) :: mW
real(default) :: alphas
real(default) :: me
real(default) :: mmu
real(default) :: mtau
real(default) :: ms
real(default) :: mc
real(default) :: mb
real(default) :: mtop
real(default) :: wtop
real(default) :: wZ
real(default) :: wW
real(default) :: khgaz
real(default) :: khgaga
real(default) :: khgg
real(default) :: xi0
real(default) :: xipm
real(default) :: mh1
real(default) :: wh1
real(default) :: mh2
real(default) :: wh2
real(default) :: mh3
real(default) :: wh3
real(default) :: mHp
real(default) :: wHp
real(default) :: R_11
real(default) :: R_12
real(default) :: R_13
real(default) :: R_21
real(default) :: R_22
real(default) :: R_23
real(default) :: R_31
real(default) :: R_32
real(default) :: R_33
real(default) :: lamb_1
real(default) :: lamb_2
real(default) :: lamb_3
real(default) :: lamb_4
real(default) :: lamb_5R
real(default) :: lamb_5I
real(default) :: lamb_6R
real(default) :: lamb_6I
real(default) :: lamb_7R
real(default) :: lamb_7I
real(default) :: zetauR
real(default) :: zetauI
real(default) :: zetadR
real(default) :: zetadI
real(default) :: zetalR
real(default) :: zetalI
real(default) :: v
real(default) :: cw
real(default) :: sw
real(default) :: ee
end type parameter_set
type(parameter_set) :: par
!!! This corresponds to 1/alpha = 137.03598949333
real(default), parameter :: &
alpha = 1.0_default/137.03598949333_default
e_em = sqrt(4.0_default * PI * alpha)
par%gf = par_array(1)
par%mZ = par_array(2)
par%mW = par_array(3)
par%alphas = par_array(4)
par%me = par_array(5)
par%mmu = par_array(6)
par%mtau = par_array(7)
par%ms = par_array(8)
par%mc = par_array(9)
par%mb = par_array(10)
par%mtop = par_array(11)
par%wtop = par_array(12)
par%wZ = par_array(13)
par%wW = par_array(14)
par%khgaz = par_array(15)
par%khgaga = par_array(16)
par%khgg = par_array(17)
par%xi0 = par_array(18)
par%xipm = par_array(19)
par%mh1 = par_array(20)
par%wh1 = par_array(21)
par%mh2 = par_array(22)
par%wh2 = par_array(23)
par%mh3 = par_array(24)
par%wh3 = par_array(25)
par%mHp = par_array(26)
par%wHp = par_array(27)
par%R_11 = par_array(28)
par%R_12 = par_array(29)
par%R_13 = par_array(30)
par%R_21 = par_array(31)
par%R_22 = par_array(32)
par%R_23 = par_array(33)
par%R_31 = par_array(34)
par%R_32 = par_array(35)
par%R_33 = par_array(36)
par%lamb_1 = par_array(37)
par%lamb_2 = par_array(38)
par%lamb_3 = par_array(39)
par%lamb_4 = par_array(40)
par%lamb_5R = par_array(41)
par%lamb_5I = par_array(42)
par%lamb_6R = par_array(43)
par%lamb_6I = par_array(44)
par%lamb_7R = par_array(45)
par%lamb_7I = par_array(46)
par%zetauR = par_array(47)
par%zetauI = par_array(48)
par%zetadR = par_array(49)
par%zetadI = par_array(50)
par%zetalR = par_array(51)
par%zetalI = par_array(52)
par%v = par_array(53)
par%cw = par_array(54)
par%sw = par_array(55)
par%ee = par_array(56)
mass(1:37) = 0
width(1:37) = 0
mass(3) = par%ms
mass(4) = par%mc
mass(5) = par%mb
mass(6) = par%mtop
width(6) = par%wtop
mass(11) = par%me
mass(13) = par%mmu
mass(15) = par%mtau
mass(23) = par%mZ
width(23) = par%wZ
mass(24) = par%mW
width(24) = par%wW
mass(25) = par%mh1
width(25) = par%wh1
mass(35) = par%mh2
width(35) = par%wh2
mass(36) = par%mh3
width(36) = par%wh3
mass(27) = par%mHp
width(27) = par%wHp
mass(26) = par%xi0 * mass(23)
width(26) = 0
mass(37) = par%xipm * mass(24)
width(37) = 0
R11 = par%R_11
R12 = par%R_12
R13 = par%R_13
R21 = par%R_21
R22 = par%R_22
R23 = par%R_23
R31 = par%R_31
R32 = par%R_32
R33 = par%R_33
lamb1 = par%lamb_1
lamb2 = par%lamb_2
lamb3 = par%lamb_3
lamb4 = par%lamb_4
lamb5R = par%lamb_5R
lamb5I = par%lamb_5I
lamb6R = par%lamb_6R
lamb6I = par%lamb_6I
lamb7R = par%lamb_7R
lamb7I = par%lamb_7I
lamb5 = par%lamb_5R + imago*par%lamb_5I
lamb6 = par%lamb_6R + imago*par%lamb_6I
lamb7 = par%lamb_7R + imago*par%lamb_7I
zetau= par%zetauR + imago*par%zetauI
zetad= par%zetadR + imago*par%zetadI
zetal= par%zetalR + imago*par%zetalI
vev = par%v
e = par%ee
sinthw = par%sw
sin2thw = par%sw**2
costhw = par%cw
tanthw = sinthw/costhw
qelep = - 1
qeup = 2.0_default / 3.0_default
qedwn = - 1.0_default / 3.0_default
g = e / sinthw
gcc = - g / 2 / sqrt (2.0_default)
gncneu(1) = - g / 2 / costhw * ( + 0.5_default)
gnclep(1) = - g / 2 / costhw * ( - 0.5_default - 2 * qelep * sin2thw)
gncup(1) = - g / 2 / costhw * ( + 0.5_default - 2 * qeup * sin2thw)
gncdwn(1) = - g / 2 / costhw * ( - 0.5_default - 2 * qedwn * sin2thw)
gncneu(2) = - g / 2 / costhw * ( + 0.5_default)
gnclep(2) = - g / 2 / costhw * ( - 0.5_default)
gncup(2) = - g / 2 / costhw * ( + 0.5_default)
gncdwn(2) = - g / 2 / costhw * ( - 0.5_default)
qlep = - e * qelep
qup = - e * qeup
qdwn = - e * qedwn
qw = e
iqw = (0,1)*qw
gzww = g * costhw
igzww = (0,1)*gzww
gwww = g
igwww = (0,1)*gwww
gw4 = gwww**2
gzzww = gzww**2
gazww = gzww * qw
gaaww = qw**2
ghtt = - mass(6) / vev
ghbb = - mass(5) / vev
ghcc = - mass(4) / vev
ghtautau = - mass(15) / vev
ghmm = - mass(13) / vev
!!!
cotanthw2 = (1.0-tanthw**2)/(2.0*tanthw)
sinthw2 = 2.0*sinthw*costhw
!Gauge boson with 2 Higgs coupligns
gAHpHm = e
gZHpHm = e * cotanthw2
gZh1h2 = - (0,1)*(e/sinthw2)*(R13*R22-R12*R23)/2.0
gZh1h3 = -(0,1)*(e/sinthw2)*(R13*R32-R12*R33)/2.0
gZh2h3 = -(0,1)*(e/sinthw2)*(R23*R32-R22*R33)/2.0
gWpHmh1 = -(g/2.0)*(R12+(0,1)*R13)
gWpHmh2 = -(g/2.0)*(R22+(0,1)*R23)
gWpHmh3 = -(g/2.0)*(R32+(0,1)*R33)
gWmHph1 = (g/2.0)*(R12-(0,1)*R13)
gWmHph2 = (g/2.0)*(R22-(0,1)*R23)
gWmHph3 = (g/2.0)*(R32-(0,1)*R33)
!Higgs + Gauge boson + Gauge boson couplings
gh1ZZ = (mass(23)**2) * R11 / vev
gh2ZZ = (mass(23)**2) * R21 / vev
gh3ZZ = (mass(23)**2) * R31 / vev
gh1WpWm = 2.0*(mass(24)**2) * R11 / vev
gh2WpWm = 2.0*(mass(24)**2) * R21 / vev
gh3WpWm = 2.0*(mass(24)**2) * R31 / vev
!Higgs+Higgs+gauge+gauge
ghhWpWm = (mass(24)/vev)**2
ghhZZ = ((mass(23)/vev)**2) / 2.0
gHpHmAA = e**2
gHpHmZZ = (e*cotanthw2)**2
gHpHmAZ = 2.0*cotanthw2*e**2
gHpHmWpWm = (g**2)/2.0
gh1HpAWm = e*g*(R12-(0,1)*R13) / 2.0
gh2HpAWm = e*g*(R22-(0,1)*R23) / 2.0
gh3HpAWm = e*g*(R32-(0,1)*R33) / 2.0
gh1HpZWm = - tanthw*gh1HpAWm
gh2HpZWm = - tanthw*gh2HpAWm
gh3HpZWm = - tanthw*gh3HpAWm
gh1HpAWmC = conjg (gh1HpAWm)
gh2HpAWmC = conjg (gh2HpAWm)
gh3HpAWmC = conjg (gh3HpAWm)
gh1HpZWmC = conjg (gh1HpZWm)
gh2HpZWmC = conjg (gh2HpZWm)
gh3HpZWmC = conjg (gh3HpZWm)
!!! Cubic Higgs
gh1HpHm = vev * (lamb3*R11 - R13*aimag(lamb7) &
+ R12*real(lamb7))
gh2HpHm = vev * (lamb3*R21 - R23*aimag(lamb7) &
+ R22*real(lamb7))
gh3HpHm = vev * (lamb3*R31 - R33*aimag(lamb7) + R32*real(lamb7))
gh111 = (vev*(2*lamb1*R11**3 + R11*(-4*R12*R13*aimag(lamb5) &
+ R13**2*(lamb3 + lamb4 - 2*real(lamb5)) + &
R12**2*(lamb3 + lamb4 + 2*real(lamb5))) + &
R11**2*(-3*R13*aimag(lamb6) + 3*R12*real(lamb6)) + &
(R12**2 + R13**2)*(-(R13*aimag(lamb7)) + R12*real(lamb7))))/2.0
gh112 = (vev*(3*R11**2*(2*lamb1*R21 - R23*aimag(lamb6) &
+ R22*real(lamb6)) + 2*R11*(R13*(-2*R22*aimag(lamb5) &
- 3*R21*aimag(lamb6) + R23*(lamb3 + lamb4 - &
2*real(lamb5))) + R12*(-2*R23*aimag(lamb5) + R22*(lamb3 + lamb4 + &
2*real(lamb5)) + 3*R21*real(lamb6))) + R13**2*(-3*R23*aimag(lamb7) + &
R21*(lamb3 + lamb4 - 2*real(lamb5)) + R22*real(lamb7)) + &
R12**2*(-(R23*aimag(lamb7)) + R21*(lamb3 + lamb4 + 2*real(lamb5)) + &
3*R22*real(lamb7)) - 2*R12*R13*(2*R21*aimag(lamb5) + &
R22*aimag(lamb7) - R23*real(lamb7))))/2.0
gh113 = (vev*(3*R11**2*(2*lamb1*R31 - R33*aimag(lamb6) + &
R32*real(lamb6)) + 2*R11*(R13*(-2*R32*aimag(lamb5) - &
3*R31*aimag(lamb6) + R33*(lamb3 + lamb4 - &
2*real(lamb5))) + R12*(-2*R33*aimag(lamb5) + R32*(lamb3 + lamb4 + &
2*real(lamb5)) + 3*R31*real(lamb6))) + R13**2*(-3*R33*aimag(lamb7) + &
R31*(lamb3 + lamb4 - 2*real(lamb5)) + R32*real(lamb7)) + &
R12**2*(-(R33*aimag(lamb7)) + R31*(lamb3 + lamb4 + 2*real(lamb5)) + &
3*R32*real(lamb7)) - 2*R12*R13*(2*R31*aimag(lamb5) + &
R32*aimag(lamb7) - R33*real(lamb7))))/2.0
gh221 = (vev*(R11*(6*lamb1*R21**2 - 4*R22*R23*aimag(lamb5) + &
R23**2*(lamb3 + lamb4 - 2*real(lamb5)) + R22**2*(lamb3 + &
lamb4 + 2*real(lamb5)) + R21*(-6*R23*aimag(lamb6) + &
6*R22*real(lamb6))) - R13*(3*R21**2*aimag(lamb6) + &
R22**2*aimag(lamb7) + 3*R23**2*aimag(lamb7) - &
2*R21*(-2*R22*aimag(lamb5) + R23*(lamb3 + lamb4 - 2*real(lamb5))) &
- 2*R22*R23*real(lamb7)) + R12*(-2*R22*R23*aimag(lamb7) + &
2*R21*(-2*R23*aimag(lamb5) + R22*(lamb3 + &
lamb4 + 2*real(lamb5))) + 3*R21**2*real(lamb6) + &
3*R22**2*real(lamb7) + R23**2*real(lamb7))))/2.0
gh222 = (vev*(2*lamb1*R21**3 + R21*(-4*R22*R23*aimag(lamb5) &
+ R23**2*(lamb3 + lamb4 - 2*real(lamb5)) + R22**2*(lamb3 + &
lamb4 + 2*real(lamb5))) + R21**2*(-3*R23*aimag(lamb6) + &
3*R22*real(lamb6)) + (R22**2 + R23**2)*(-(R23*aimag(lamb7)) &
+ R22*real(lamb7))))/2.0
gh223 = (vev*(3*R21**2*(2*lamb1*R31 - R33*aimag(lamb6) + &
R32*real(lamb6)) + 2*R21*(R23*(-2*R32*aimag(lamb5) - &
3*R31*aimag(lamb6) + R33*(lamb3 + lamb4 - &
2*real(lamb5))) + R22*(-2*R33*aimag(lamb5) + R32*(lamb3 + lamb4 + &
2*real(lamb5)) + 3*R31*real(lamb6))) + R23**2*(-3*R33*aimag(lamb7) + &
R31*(lamb3 + lamb4 - 2*real(lamb5)) + R32*real(lamb7)) + &
R22**2*(-(R33*aimag(lamb7)) + R31*(lamb3 + lamb4 + 2*real(lamb5)) + &
3*R32*real(lamb7)) - 2*R22*R23*(2*R31*aimag(lamb5) + &
R32*aimag(lamb7) - R33*real(lamb7))))/2.0
gh331 = (vev*(R11*(6*lamb1*R31**2 - 4*R32*R33*aimag(lamb5) + &
R33**2*(lamb3 + lamb4 - 2*real(lamb5)) + R32**2*(lamb3 + lamb4 &
+ 2*real(lamb5)) + R31*(-6*R33*aimag(lamb6) + 6*R32*real(lamb6))) &
- R13*(3*R31**2*aimag(lamb6) + R32**2*aimag(lamb7) + &
3*R33**2*aimag(lamb7) - 2*R31*(-2*R32*aimag(lamb5) + &
R33*(lamb3 + lamb4 - 2*real(lamb5))) - 2*R32*R33*real(lamb7)) + &
R12*(-2*R32*R33*aimag(lamb7) + 2*R31*(-2*R33*aimag(lamb5) + &
R32*(lamb3 + lamb4 + 2*real(lamb5))) + 3*R31**2*real(lamb6) + &
3*R32**2*real(lamb7) + R33**2*real(lamb7))))/2.0
gh332 = (vev*(R21*(6*lamb1*R31**2 - 4*R32*R33*aimag(lamb5) + &
R33**2*(lamb3 + lamb4 - 2*real(lamb5)) + R32**2*(lamb3 + lamb4 &
+ 2*real(lamb5)) + R31*(-6*R33*aimag(lamb6) + 6*R32*real(lamb6))) &
- R23*(3*R31**2*aimag(lamb6) + R32**2*aimag(lamb7) + &
3*R33**2*aimag(lamb7) - 2*R31*(-2*R32*aimag(lamb5) + &
R33*(lamb3 + lamb4 - 2*real(lamb5))) - 2*R32*R33*real(lamb7)) + &
R22*(-2*R32*R33*aimag(lamb7) + 2*R31*(-2*R33*aimag(lamb5) + &
R32*(lamb3 + lamb4 + 2*real(lamb5))) + 3*R31**2*real(lamb6) + &
3*R32**2*real(lamb7) + R33**2*real(lamb7))))/2.0
gh333 = (vev*(2*lamb1*R31**3 + R31*(-4*R32*R33*aimag(lamb5) + &
R33**2*(lamb3 + lamb4 - 2*real(lamb5)) + R32**2*(lamb3 + &
lamb4 + 2*real(lamb5))) + R31**2*(-3*R33*aimag(lamb6) + &
3*R32*real(lamb6)) + (R32**2 + R33**2)*(-(R33*aimag(lamb7)) + &
R32*real(lamb7))))/2.0
gh123 = vev*(R11*(R23*(-2*R32*aimag(lamb5) - 3*R31*aimag(lamb6) + &
R33*(lamb3 + lamb4 - 2*real(lamb5))) + R22*(-2*R33*aimag(lamb5) &
+ R32*(lamb3 + lamb4 + 2*real(lamb5)) + 3*R31*real(lamb6)) + &
3*R21*(2*lamb1*R31 - R33*aimag(lamb6) + R32*real(lamb6))) + &
R13*(R21*(-2*R32*aimag(lamb5) - 3*R31*aimag(lamb6) + R33*(lamb3 &
+ lamb4 - 2*real(lamb5))) + R23*(-3*R33*aimag(lamb7) + &
R31*(lamb3 + lamb4 - 2*real(lamb5)) + &
R32*real(lamb7)) + R22*(-2*R31*aimag(lamb5) - R32*aimag(lamb7) + &
R33*real(lamb7))) + R12*(R21*(-2*R33*aimag(lamb5) + R32*(lamb3 + &
lamb4 + 2*real(lamb5)) + 3*R31*real(lamb6)) + &
R22*(-(R33*aimag(lamb7)) + R31*(lamb3 + lamb4 + 2*real(lamb5)) &
+ 3*R32*real(lamb7)) + R23*(-2*R31*aimag(lamb5) - R32*aimag(lamb7) &
+ R33*real(lamb7))))
!!! Quartic charged Higgs
gHpHmHpHm = lamb2
gHpHm11 = (lamb3*R11**2)/2. + lamb7R*R11*R12 + lamb2*R12**2 - &
lamb7I*R11*R13 + lamb2*R13**2
gHpHm12 = lamb3*R11*R21 + lamb7R*R12*R21 - lamb7I*R13*R21 + &
lamb7R*R11*R22 + 2*lamb2*R12*R22 - lamb7I*R11*R23 + 2*lamb2*R13*R23
gHpHm13 = lamb3*R11*R31 + lamb7R*R12*R31 - lamb7I*R13*R31 + &
lamb7R*R11*R32 + 2*lamb2*R12*R32 - lamb7I*R11*R33 + 2*lamb2*R13*R33
gHpHm22 = (lamb3*R21**2)/2. + lamb7R*R21*R22 + lamb2*R22**2 - &
lamb7I*R21*R23 + lamb2*R23**2
gHpHm23 = lamb3*R21*R31 + lamb7R*R22*R31 - lamb7I*R23*R31 + &
lamb7R*R21*R32 + 2*lamb2*R22*R32 - lamb7I*R21*R33 + 2*lamb2*R23*R33
gHpHm33 = (lamb3*R31**2)/2. + lamb7R*R31*R32 + lamb2*R32**2 - &
lamb7I*R31*R33 + lamb2*R33**2
!!! Quartic neutral higgs
gh1111 = (lamb1*R11**4 + R11**3*(2*lamb6R*R12 - 2*lamb6I*R13) + &
2*R11*(lamb7R*R12 - lamb7I*R13)*(R12**2 + R13**2) + lamb2*(R12**2 + &
R13**2)**2 + R11**2*((lamb3 + lamb4 + 2*lamb5R)*R12**2 - &
4*lamb5I*R12*R13 + (lamb3 + lamb4 - 2*lamb5R)*R13**2))/4.0
gh1112 = (R11**3*(2*lamb1*R21 + lamb6R*R22 - lamb6I*R23) + &
(R12**2 + R13**2)*(lamb7R*R12*R21 - lamb7I*R13*R21 + &
2*lamb2*R12*R22 + 2*lamb2*R13*R23) + R11**2*(R12*(3*lamb6R*R21 + &
(lamb3 + lamb4 + 2*lamb5R)*R22 - 2*lamb5I*R23) + &
R13*(-3*lamb6I*R21 - 2*lamb5I*R22 + (lamb3 + lamb4 - &
2*lamb5R)*R23)) + R11*(R13**2*((lamb3 + lamb4 - &
2*lamb5R)*R21 + lamb7R*R22 - 3*lamb7I*R23) + R12**2*((lamb3 + &
lamb4 + 2*lamb5R)*R21 + 3*lamb7R*R22 - lamb7I*R23) - &
2*R12*R13*(2*lamb5I*R21 + lamb7I*R22 - lamb7R*R23)))/2.0
gh1113 = (R11**3*(2*lamb1*R31 + lamb6R*R32 - lamb6I*R33) + (R12**2 + &
R13**2)*(lamb7R*R12*R31 - lamb7I*R13*R31 + 2*lamb2*R12*R32 + &
2*lamb2*R13*R33) + R11**2*(R12*(3*lamb6R*R31 + (lamb3 + lamb4 + &
2*lamb5R)*R32 - 2*lamb5I*R33) + R13*(-3*lamb6I*R31 - 2*lamb5I*R32 + &
(lamb3 + lamb4 - 2*lamb5R)*R33)) + R11*(R13**2*((lamb3 + lamb4 - &
2*lamb5R)*R31 + lamb7R*R32 - 3*lamb7I*R33) + R12**2*((lamb3 + &
lamb4 + 2*lamb5R)*R31 + 3*lamb7R*R32 - lamb7I*R33) - &
2*R12*R13*(2*lamb5I*R31 + lamb7I*R32 - lamb7R*R33)))/2.0
gh1122 = (4*R12*R13*(-(lamb5I*R21**2) - lamb7I*R21*R22 + lamb7R*R21*R23 + &
2*lamb2*R22*R23) + R12**2*((lamb3 + lamb4 + 2*lamb5R)*R21**2 + &
6*lamb7R*R21*R22 + 6*lamb2*R22**2 - 2*lamb7I*R21*R23 + &
2*lamb2*R23**2) + R13**2*((lamb3 + lamb4 - 2*lamb5R)*R21**2 + &
2*lamb7R*R21*R22 + 2*lamb2*R22**2 - 6*lamb7I*R21*R23 + &
6*lamb2*R23**2) + R11**2*(6*lamb1*R21**2 + (lamb3 + lamb4 + &
2*lamb5R)*R22**2 - 4*lamb5I*R22*R23 + (lamb3 + lamb4 - &
2*lamb5R)*R23**2 + R21*(6*lamb6R*R22 - 6*lamb6I*R23)) + &
2*R11*(R12*(3*lamb6R*R21**2 + 3*lamb7R*R22**2 - 2*lamb7I*R22*R23 &
+ lamb7R*R23**2 + 2*R21*((lamb3 + lamb4 + 2*lamb5R)*R22 &
- 2*lamb5I*R23)) - R13*(3*lamb6I*R21**2 + &
lamb7I*R22**2 - 2*lamb7R*R22*R23 + 3*lamb7I*R23**2 - &
2*R21*(-2*lamb5I*R22 + (lamb3 + lamb4 - 2*lamb5R)*R23))))/4.0
gh1123 = (2*R12*R13*(-2*lamb5I*R21*R31 - lamb7I*R22*R31 + &
lamb7R*R23*R31 - lamb7I*R21*R32 + 2*lamb2*R23*R32 + lamb7R*R21*R33 &
+ 2*lamb2*R22*R33) + R11**2*(R22*(3*lamb6R*R31 + (lamb3 + lamb4 + &
2*lamb5R)*R32 - 2*lamb5I*R33) + R23*(-3*lamb6I*R31 - &
2*lamb5I*R32 + (lamb3 + lamb4 - 2*lamb5R)*R33) + &
3*R21*(2*lamb1*R31 + lamb6R*R32 - lamb6I*R33)) + &
R13**2*(R22*(lamb7R*R31 + 2*lamb2*R32) + R23*(-3*lamb7I*R31 &
+ 6*lamb2*R33) + R21*((lamb3 + lamb4 - 2*lamb5R)*R31 + &
lamb7R*R32 - 3*lamb7I*R33)) + R12**2*(3*R22*(lamb7R*R31 + &
2*lamb2*R32) + R23*(-(lamb7I*R31) + 2*lamb2*R33) + R21*((lamb3 &
+ lamb4 + 2*lamb5R)*R31 + 3*lamb7R*R32 - lamb7I*R33)) + &
2*R11*(R13*(R21*(-3*lamb6I*R31 - 2*lamb5I*R32 + (lamb3 + lamb4 - &
2*lamb5R)*R33) + R23*((lamb3 + lamb4 - 2*lamb5R)*R31 + lamb7R*R32 - &
3*lamb7I*R33) + R22*(-2*lamb5I*R31 - lamb7I*R32 + lamb7R*R33)) + &
R12*(R21*(3*lamb6R*R31 + (lamb3 + lamb4 + 2*lamb5R)*R32 - &
2*lamb5I*R33) + R22*((lamb3 + lamb4 + 2*lamb5R)*R31 + &
3*lamb7R*R32 - lamb7I*R33) + R23*(-2*lamb5I*R31 - &
lamb7I*R32 + lamb7R*R33))))/2.0
gh1133 = (4*R12*R13*(-(lamb5I*R31**2) - lamb7I*R31*R32 + lamb7R*R31*R33 + &
2*lamb2*R32*R33) + R12**2*((lamb3 + lamb4 + 2*lamb5R)*R31**2 + &
6*lamb7R*R31*R32 + 6*lamb2*R32**2 - 2*lamb7I*R31*R33 + &
2*lamb2*R33**2) + R13**2*((lamb3 + lamb4 - 2*lamb5R)*R31**2 + &
2*lamb7R*R31*R32 + 2*lamb2*R32**2 - 6*lamb7I*R31*R33 + &
6*lamb2*R33**2) + R11**2*(6*lamb1*R31**2 + (lamb3 + lamb4 + &
2*lamb5R)*R32**2 - 4*lamb5I*R32*R33 + (lamb3 + lamb4 - &
2*lamb5R)*R33**2 + R31*(6*lamb6R*R32 - 6*lamb6I*R33)) + &
2*R11*(R12*(3*lamb6R*R31**2 + 3*lamb7R*R32**2 - 2*lamb7I*R32*R33 &
+ lamb7R*R33**2 + 2*R31*((lamb3 + lamb4 + 2*lamb5R)*R32 - &
2*lamb5I*R33)) - R13*(3*lamb6I*R31**2 + lamb7I*R32**2 - &
2*lamb7R*R32*R33 + 3*lamb7I*R33**2 - 2*R31*(-2*lamb5I*R32 + &
(lamb3 + lamb4 - 2*lamb5R)*R33))))/4.0
gh1222 = (R13*(-(lamb6I*R21**3) + (lamb3 + lamb4 - &
2*lamb5R)*R21**2*R23 - 3*lamb7I*R21*R23**2 + 2*lamb2*R23**3 + &
R22**2*(-(lamb7I*R21) + 2*lamb2*R23) + 2*R21*R22*(-(lamb5I*R21) &
+ lamb7R*R23)) + R11*(2*lamb1*R21**3 + R21**2*(3*lamb6R*R22 - &
3*lamb6I*R23) + (lamb7R*R22 - lamb7I*R23)*(R22**2 + R23**2) + &
R21*((lamb3 + lamb4 + 2*lamb5R)*R22**2 - 4*lamb5I*R22*R23 + &
(lamb3 + lamb4 - 2*lamb5R)*R23**2)) + R12*(3*lamb7R*R21*R22**2 &
+ 2*lamb2*R22**3 + R22*((lamb3 + lamb4 + 2*lamb5R)*R21**2 &
- 2*lamb7I*R21*R23 + 2*lamb2*R23**2) + R21*(lamb6R*R21**2 - &
2*lamb5I*R21*R23 + lamb7R*R23**2)))/2.0
gh1223 = (R13*(R22**2*(-(lamb7I*R31) + 2*lamb2*R33) + &
R23**2*(-3*lamb7I*R31 + 6*lamb2*R33) + R21**2*(-3*lamb6I*R31 - &
2*lamb5I*R32 + (lamb3 + lamb4 - 2*lamb5R)*R33) + 2*R21*R23*((lamb3 &
+ lamb4 - 2*lamb5R)*R31 + lamb7R*R32 - 3*lamb7I*R33) + &
2*R22*(-2*lamb5I*R21*R31 + lamb7R*R23*R31 - lamb7I*R21*R32 + &
2*lamb2*R23*R32 + lamb7R*R21*R33)) + R11*(3*R21**2*(2*lamb1*R31 &
+ lamb6R*R32 - lamb6I*R33) + R23**2*((lamb3 + &
lamb4 - 2*lamb5R)*R31 + lamb7R*R32 - 3*lamb7I*R33) + &
R22**2*((lamb3 + lamb4 + 2*lamb5R)*R31 + 3*lamb7R*R32 - &
lamb7I*R33) - 2*R22*R23*(2*lamb5I*R31 + lamb7I*R32 - &
lamb7R*R33) + 2*R21*(R22*(3*lamb6R*R31 + (lamb3 + lamb4 &
+ 2*lamb5R)*R32 - 2*lamb5I*R33) + R23*(-3*lamb6I*R31 - &
2*lamb5I*R32 + (lamb3 + lamb4 - 2*lamb5R)*R33))) + &
R12*(3*R22**2*(lamb7R*R31 + 2*lamb2*R32) + R23**2*(lamb7R*R31 &
+ 2*lamb2*R32) + R21**2*(3*lamb6R*R31 + (lamb3 + lamb4 + &
2*lamb5R)*R32 - 2*lamb5I*R33) - 2*R21*R23*(2*lamb5I*R31 + &
lamb7I*R32 - lamb7R*R33) + 2*R22*(R23*(-(lamb7I*R31) + &
2*lamb2*R33) + R21*((lamb3 + lamb4 + 2*lamb5R)*R31 + &
3*lamb7R*R32 - lamb7I*R33))))/2.0
gh1233 = (R12*(2*R23*(-(lamb5I*R31**2) - lamb7I*R31*R32 + &
lamb7R*R31*R33 + 2*lamb2*R32*R33) + R22*((lamb3 + lamb4 + &
2*lamb5R)*R31**2 + 6*lamb7R*R31*R32 + 6*lamb2*R32**2 - &
2*lamb7I*R31*R33 + 2*lamb2*R33**2) + R21*(3*lamb6R*R31**2 + &
3*lamb7R*R32**2 - 2*lamb7I*R32*R33 + lamb7R*R33**2 + &
2*R31*((lamb3 + lamb4 + 2*lamb5R)*R32 - 2*lamb5I*R33))) + &
R13*(2*R22*(-(lamb5I*R31**2) - lamb7I*R31*R32 + lamb7R*R31*R33 + &
2*lamb2*R32*R33) + R23*((lamb3 + lamb4 - 2*lamb5R)*R31**2 + &
2*lamb7R*R31*R32 + 2*lamb2*R32**2 - 6*lamb7I*R31*R33 + &
6*lamb2*R33**2) - R21*(3*lamb6I*R31**2 + lamb7I*R32**2 - &
2*lamb7R*R32*R33 + 3*lamb7I*R33**2 - &
2*R31*(-2*lamb5I*R32 + (lamb3 + lamb4 - 2*lamb5R)*R33))) + &
R11*(R22*(3*lamb6R*R31**2 + 3*lamb7R*R32**2 - 2*lamb7I*R32*R33 + &
lamb7R*R33**2 + 2*R31*((lamb3 + lamb4 + 2*lamb5R)*R32 - &
2*lamb5I*R33)) - R23*(3*lamb6I*R31**2 + lamb7I*R32**2 - &
2*lamb7R*R32*R33 + 3*lamb7I*R33**2 - 2*R31*(-2*lamb5I*R32 + &
(lamb3 + lamb4 - 2*lamb5R)*R33)) + R21*(6*lamb1*R31**2 + &
(lamb3 + lamb4 + 2*lamb5R)*R32**2 - 4*lamb5I*R32*R33 + &
(lamb3 + lamb4 - 2*lamb5R)*R33**2 + R31*(6*lamb6R*R32 - &
6*lamb6I*R33))))/2.
gh1333 = R13*(-(lamb6I*R31**3) + (lamb3 + lamb4 - 2*lamb5R)*R31**2*R33 &
- 3*lamb7I*R31*R33**2 + 2*lamb2*R33**3 + R32**2*(-(lamb7I*R31) &
+ 2*lamb2*R33) + 2*R31*R32*(-(lamb5I*R31) + lamb7R*R33)) + &
R11*(2*lamb1*R31**3 + R31**2*(3*lamb6R*R32 - 3*lamb6I*R33) + &
(lamb7R*R32 - lamb7I*R33)*(R32**2 + R33**2) + R31*((lamb3 + &
lamb4 + 2*lamb5R)*R32**2 - 4*lamb5I*R32*R33 + (lamb3 + lamb4 &
- 2*lamb5R)*R33**2)) + R12*(3*lamb7R*R31*R32**2 + &
2*lamb2*R32**3 + R32*((lamb3 + lamb4 + 2*lamb5R)*R31**2 - &
2*lamb7I*R31*R33 + 2*lamb2*R33**2) + R31*(lamb6R*R31**2 - &
2*lamb5I*R31*R33 + lamb7R*R33**2))
gh2222 = (lamb1*R21**4 + R21**3*(2*lamb6R*R22 - 2*lamb6I*R23) + &
2*R21*(lamb7R*R22 - lamb7I*R23)*(R22**2 + R23**2) + lamb2*(R22**2 + &
R23**2)**2 + R21**2*((lamb3 + lamb4 + 2*lamb5R)*R22**2 - &
4*lamb5I*R22*R23 + (lamb3 + lamb4 - 2*lamb5R)*R23**2))/4.0
gh2223 = (R21**3*(2*lamb1*R31 + lamb6R*R32 - lamb6I*R33) + &
(R22**2 + R23**2)*(lamb7R*R22*R31 - lamb7I*R23*R31 + &
2*lamb2*R22*R32 + 2*lamb2*R23*R33) + R21**2*(R22*(3*lamb6R*R31 &
+ (lamb3 + lamb4 + 2*lamb5R)*R32 - 2*lamb5I*R33) + &
R23*(-3*lamb6I*R31 - 2*lamb5I*R32 + &
(lamb3 + lamb4 - 2*lamb5R)*R33)) + R21*(R23**2*((lamb3 + lamb4 - &
2*lamb5R)*R31 + lamb7R*R32 - 3*lamb7I*R33) + R22**2*((lamb3 &
+ lamb4 + 2*lamb5R)*R31 + 3*lamb7R*R32 - lamb7I*R33) - &
2*R22*R23*(2*lamb5I*R31 + lamb7I*R32 - lamb7R*R33)))/2.0
gh2233 = (4*R22*R23*(-(lamb5I*R31**2) - lamb7I*R31*R32 + &
lamb7R*R31*R33 + 2*lamb2*R32*R33) + R22**2*((lamb3 + &
lamb4 + 2*lamb5R)*R31**2 + 6*lamb7R*R31*R32 + 6*lamb2*R32**2 &
- 2*lamb7I*R31*R33 + 2*lamb2*R33**2) + R23**2*((lamb3 + lamb4 &
- 2*lamb5R)*R31**2 + 2*lamb7R*R31*R32 + 2*lamb2*R32**2 - &
6*lamb7I*R31*R33 + 6*lamb2*R33**2) + R21**2*(6*lamb1*R31**2 + &
(lamb3 + lamb4 + 2*lamb5R)*R32**2 - 4*lamb5I*R32*R33 + &
(lamb3 + lamb4 - 2*lamb5R)*R33**2 + R31*(6*lamb6R*R32 - &
6*lamb6I*R33)) + 2*R21*(R22*(3*lamb6R*R31**2 + 3*lamb7R*R32**2 - &
2*lamb7I*R32*R33 + lamb7R*R33**2 + 2*R31*((lamb3 + lamb4 + &
2*lamb5R)*R32 - 2*lamb5I*R33)) - R23*(3*lamb6I*R31**2 + &
lamb7I*R32**2 - 2*lamb7R*R32*R33 + 3*lamb7I*R33**2 - &
2*R31*(-2*lamb5I*R32 + (lamb3 + lamb4 - 2*lamb5R)*R33))))/4.0
gh2333 = (R23*(-(lamb6I*R31**3) + (lamb3 + lamb4 - &
2*lamb5R)*R31**2*R33 - 3*lamb7I*R31*R33**2 + 2*lamb2*R33**3 &
+ R32**2*(-(lamb7I*R31) + 2*lamb2*R33) + 2*R31*R32*(-(lamb5I*R31) &
+ lamb7R*R33)) + R21*(2*lamb1*R31**3 + R31**2*(3*lamb6R*R32 &
- 3*lamb6I*R33) + (lamb7R*R32 - lamb7I*R33)*(R32**2 + R33**2) + &
R31*((lamb3 + lamb4 + 2*lamb5R)*R32**2 - 4*lamb5I*R32*R33 + &
(lamb3 + lamb4 - 2*lamb5R)*R33**2)) + R22*(3*lamb7R*R31*R32**2 &
+ 2*lamb2*R32**3 + R32*((lamb3 + lamb4 + 2*lamb5R)*R31**2 - &
2*lamb7I*R31*R33 + 2*lamb2*R33**2) + R31*(lamb6R*R31**2 - &
2*lamb5I*R31*R33 + lamb7R*R33**2)))/2.0
gh3333 = (lamb1*R31**4 + R31**3*(2*lamb6R*R32 - 2*lamb6I*R33) + &
2*R31*(lamb7R*R32 - lamb7I*R33)*(R32**2 + R33**2) + lamb2*(R32**2 + &
R33**2)**2 + R31**2*((lamb3 + lamb4 + 2*lamb5R)*R32**2 - &
4*lamb5I*R32*R33 + (lamb3 + lamb4 - 2*lamb5R)*R33**2))/4.
!!! Yukawa couplings. Up-type quark, up-type quark
gh1uu(1) = -((2*R11 + (R12 + imago*R13)*zetau + (R12 - &
imago*R13)*conjg(zetau))*mass(2))/(2.*vev)
gh2uu(1) = -((2*R21 + (R22 + imago*R23)*zetau + (R22 - &
imago*R23)*conjg(zetau))*mass(2))/(2.*vev)
gh3uu(1) = -((2*R31 + (R32 + imago*R33)*zetau + (R32 - &
imago*R33)*conjg(zetau))*mass(2))/(2.*vev)
gh1cc(1) = -((2*R11 + (R12 + imago*R13)*zetau + (R12 - &
imago*R13)*conjg(zetau))*mass(4))/(2.*vev)
gh2cc(1) = -((2*R21 + (R22 + imago*R23)*zetau + (R22 - &
imago*R23)*conjg(zetau))*mass(4))/(2.*vev)
gh3cc(1) = -((2*R31 + (R32 + imago*R33)*zetau + (R32 - &
imago*R33)*conjg(zetau))*mass(4))/(2.*vev)
gh1tt(1) = -((2*R11 + (R12 + imago*R13)*zetau + (R12 - &
imago*R13)*conjg(zetau))*mass(6))/(2.*vev)
gh2tt(1) = -((2*R21 + (R22 + imago*R23)*zetau + (R22 - &
imago*R23)*conjg(zetau))*mass(6))/(2.*vev)
gh3tt(1) = -((2*R31 + (R32 + imago*R33)*zetau + (R32 - &
imago*R33)*conjg(zetau))*mass(6))/(2.*vev)
gh1uu(2) = ((R12*zetau + imago*R13*zetau - R12*conjg(zetau) + &
imago*R13*conjg(zetau))*mass(2))/(2.*vev)
gh2uu(2) = ((R22*zetau + imago*R23*zetau - R22*conjg(zetau) + &
imago*R23*conjg(zetau))*mass(2))/(2.*vev)
gh3uu(2) = ((R32*zetau + imago*R33*zetau - R32*conjg(zetau) + &
imago*R33*conjg(zetau))*mass(2))/(2.*vev)
gh1cc(2) = ((R12*zetau + imago*R13*zetau - R12*conjg(zetau) + &
imago*R13*conjg(zetau))*mass(4))/(2.*vev)
gh2cc(2) = ((R22*zetau + imago*R23*zetau - R22*conjg(zetau) + &
imago*R23*conjg(zetau))*mass(4))/(2.*vev)
gh3cc(2) = ((R32*zetau + imago*R33*zetau - R32*conjg(zetau) + &
imago*R33*conjg(zetau))*mass(4))/(2.*vev)
gh1tt(2) = ((R12*zetau + imago*R13*zetau - R12*conjg(zetau) + &
imago*R13*conjg(zetau))*mass(6))/(2.*vev)
gh2tt(2) = ((R22*zetau + imago*R23*zetau - R22*conjg(zetau) + &
imago*R23*conjg(zetau))*mass(6))/(2.*vev)
gh3tt(2) = ((R32*zetau + imago*R33*zetau - R32*conjg(zetau) + &
imago*R33*conjg(zetau))*mass(6))/(2.*vev)
!!! Yukawa couplings. Down-type quark, down-type quark
gh1dd(1) = -((2*R11 + (R12 + imago*R13)*zetad + (R12 - &
imago*R13)*conjg(zetad))*mass(1))/(2.*vev)
gh2dd(1) = -((2*R21 + (R22 + imago*R23)*zetad + (R22 - &
imago*R23)*conjg(zetad))*mass(1))/(2.*vev)
gh3dd(1) = -((2*R31 + (R32 + imago*R33)*zetad + (R32 - &
imago*R33)*conjg(zetad))*mass(1))/(2.*vev)
gh1ss(1) = -((2*R11 + (R12 + imago*R13)*zetad + (R12 - &
imago*R13)*conjg(zetad))*mass(3))/(2.*vev)
gh2ss(1) = -((2*R21 + (R22 + imago*R23)*zetad + (R22 - &
imago*R23)*conjg(zetad))*mass(3))/(2.*vev)
gh3ss(1) = -((2*R31 + (R32 + imago*R33)*zetad + (R32 - &
imago*R33)*conjg(zetad))*mass(3))/(2.*vev)
gh1bb(1) = -((2*R11 + (R12 + imago*R13)*zetad + (R12 - &
imago*R13)*conjg(zetad))*mass(5))/(2.*vev)
gh2bb(1) = -((2*R21 + (R22 + imago*R23)*zetad + (R22 - &
imago*R23)*conjg(zetad))*mass(5))/(2.*vev)
gh3bb(1) = -((2*R31 + (R32 + imago*R33)*zetad + (R32 - &
imago*R33)*conjg(zetad))*mass(5))/(2.*vev)
gh1dd(2) = ((-((R12 + imago*R13)*zetad) + (R12 - &
imago*R13)*conjg(zetad))*mass(1))/(2.*vev)
gh2dd(2) = ((-((R22 + imago*R23)*zetad) + (R22 - &
imago*R23)*conjg(zetad))*mass(1))/(2.*vev)
gh3dd(2) = ((-((R32 + imago*R33)*zetad) + (R32 - &
imago*R33)*conjg(zetad))*mass(1))/(2.*vev)
gh1ss(2) = ((-((R12 + imago*R13)*zetad) + (R12 - &
imago*R13)*conjg(zetad))*mass(3))/(2.*vev)
gh2ss(2) = ((-((R22 + imago*R23)*zetad) + (R22 - &
imago*R23)*conjg(zetad))*mass(3))/(2.*vev)
gh3ss(2) = ((-((R32 + imago*R33)*zetad) + (R32 - &
imago*R33)*conjg(zetad))*mass(3))/(2.*vev)
gh1bb(2) = ((-((R12 + imago*R13)*zetad) + (R12 - &
imago*R13)*conjg(zetad))*mass(5))/(2.*vev)
gh2bb(2) = ((-((R22 + imago*R23)*zetad) + (R22 - &
imago*R23)*conjg(zetad))*mass(5))/(2.*vev)
gh3bb(2) = ((-((R32 + imago*R33)*zetad) + (R32 - &
imago*R33)*conjg(zetad))*mass(5))/(2.*vev)
!!! Yukawa couplings: up quark down quark Higgs charged
ghud(1) = -(sqrt2*(zetad*mass(1) - zetau*mass(2)))/(2.*vev)
ghcs(1) = -(sqrt2*(zetad*mass(3) - zetau*mass(4)))/(2.*vev)
ghtb(1) = -(sqrt2*(zetad*mass(5) - zetau*mass(6)))/(2.*vev)
ghud(2) = -(sqrt2*(zetad*mass(1) + zetau*mass(2)))/(2.*vev)
ghcs(2) = -(sqrt2*(zetad*mass(3) + zetau*mass(4)))/(2.*vev)
ghtb(2) = -(sqrt2*(zetad*mass(5) + zetau*mass(6)))/(2.*vev)
ghdu(1) = (sqrt2*(-(conjg(zetad)*mass(1)) + &
conjg(zetau)*mass(2)))/(2.*vev)
ghsc(1) = (sqrt2*(-(conjg(zetad)*mass(3)) + &
conjg(zetau)*mass(4)))/(2.*vev)
ghbt(1) = (sqrt2*(-(conjg(zetad)*mass(5)) + &
conjg(zetau)*mass(6)))/(2.*vev)
ghdu(2) = (sqrt2*(conjg(zetad)*mass(1) + &
conjg(zetau)*mass(2)))/(2.*vev)
ghsc(2) = (sqrt2*(conjg(zetad)*mass(3) + &
conjg(zetau)*mass(4)))/(2.*vev)
ghbt(2) = (sqrt2*(conjg(zetad)*mass(5) + &
conjg(zetau)*mass(6)))/(2.*vev)
!!! Yukawa couplings: Charged Higgs + Neutrino + Electron
ghe1n1 = -((conjg(zetal)*mass(11))/(sqrt2*vev))
ghe2n2 = -((conjg(zetal)*mass(13))/(sqrt2*vev))
ghe3n3 = -((conjg(zetal)*mass(15))/(sqrt2*vev))
ghn1e1 = -((zetal*mass(11))/(sqrt2*vev))
ghn2e2 = -((zetal*mass(13))/(sqrt2*vev))
ghn3e3 = -((zetal*mass(15))/(sqrt2*vev))
!!! Yukawa couplings: Neutral Higgs + Electron + Elentron
gh1e1e1(1) = -((2*R11 + (R12 + imago*R13)*zetal + (R12 - &
imago*R13)*conjg(zetal))*mass(11))/(2.*vev)
gh2e1e1(1) = -((2*R21 + (R22 + imago*R23)*zetal + (R22 - &
imago*R23)*conjg(zetal))*mass(11))/(2.*vev)
gh3e1e1(1) = -((2*R31 + (R32 + imago*R33)*zetal + (R32 - &
imago*R33)*conjg(zetal))*mass(11))/(2.*vev)
gh1e2e2(1) = -((2*R11 + (R12 + imago*R13)*zetal + (R12 - &
imago*R13)*conjg(zetal))*mass(13))/(2.*vev)
gh2e2e2(1) = -((2*R21 + (R22 + imago*R23)*zetal + (R22 - &
imago*R23)*conjg(zetal))*mass(13))/(2.*vev)
gh3e2e2(1) = -((2*R31 + (R32 + imago*R33)*zetal + (R32 - &
imago*R33)*conjg(zetal))*mass(13))/(2.*vev)
gh1e3e3(1) = -((2*R11 + (R12 + imago*R13)*zetal + (R12 - &
imago*R13)*conjg(zetal))*mass(15))/(2.*vev)
gh2e3e3(1) = -((2*R21 + (R22 + imago*R23)*zetal + (R22 - &
imago*R23)*conjg(zetal))*mass(15))/(2.*vev)
gh3e3e3(1) = -((2*R31 + (R32 + imago*R33)*zetal + (R32 - &
imago*R33)*conjg(zetal))*mass(15))/(2.*vev)
gh1e1e1(2) = ((-((R12 + imago*R13)*zetal) + (R12 - &
imago*R13)*conjg(zetal))*mass(11))/(2.*vev)
gh2e1e1(2) = ((-((R22 + imago*R23)*zetal) + (R22 - &
imago*R23)*conjg(zetal))*mass(11))/(2.*vev)
gh3e1e1(2) = ((-((R32 + imago*R33)*zetal) + (R32 - &
imago*R33)*conjg(zetal))*mass(11))/(2.*vev)
gh1e2e2(2) = ((-((R12 + imago*R13)*zetal) + (R12 - &
imago*R13)*conjg(zetal))*mass(13))/(2.*vev)
gh2e2e2(2) = ((-((R22 + imago*R23)*zetal) + (R22 - &
imago*R23)*conjg(zetal))*mass(13))/(2.*vev)
gh3e2e2(2) = ((-((R32 + imago*R33)*zetal) + (R32 - &
imago*R33)*conjg(zetal))*mass(13))/(2.*vev)
gh1e3e3(2) = ((-((R12 + imago*R13)*zetal) + (R12 - &
imago*R13)*conjg(zetal))*mass(15))/(2.*vev)
gh2e3e3(2) = ((-((R22 + imago*R23)*zetal) + (R22 - &
imago*R23)*conjg(zetal))*mass(15))/(2.*vev)
gh3e3e3(2) = ((-((R32 + imago*R33)*zetal) + (R32 - &
imago*R33)*conjg(zetal))*mass(15))/(2.*vev)
!!! Color flow basis, divide by sqrt(2)
gs = sqrt(2.0_default*PI*par%alphas)
igs = cmplx (0.0_default, 1.0_default, kind=default) * gs
end subroutine import_from_whizard
subroutine model_update_alpha_s (alpha_s)
real(default), intent(in) :: alpha_s
gs = sqrt(2.0_default*PI*alpha_s)
igs = cmplx (0.0_default, 1.0_default, kind=default) * gs
end subroutine model_update_alpha_s
end module parameters_2hdm

File Metadata

Mime Type
text/plain
Expires
Sat, May 3, 6:34 AM (1 d, 17 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4973366
Default Alt Text
parameters.2HDM.f90 (35 KB)

Event Timeline