Page Menu
Home
HEPForge
Search
Configure Global Search
Log In
Files
F7879174
No One
Temporary
Actions
View File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Flag For Later
Size
44 KB
Subscribers
None
View Options
Index: trunk/src/models/parameters.AltH.f90
===================================================================
--- trunk/src/models/parameters.AltH.f90 (revision 0)
+++ trunk/src/models/parameters.AltH.f90 (revision 5352)
@@ -0,0 +1,362 @@
+! $Id: parameters.VBS.f90, 2013/09/23 13:11:48 msekulla Exp $
+!
+! Copyright (C) 1999-2014 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>
+! Marco Sekulla <sekulla@physik.uni-siegen.de> (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_alth
+ use kinds
+ use constants
+ implicit none
+ private
+
+ real(default), dimension(55), 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, &
+ iqw, igzww, igwww, gw4, gzzww, gazww, gaaww, &
+ gszz, gszzt, gsww, gswwt, gssww, gsszz, &
+ gpnww, gpnzz, gpwz, gpww, &
+ gfww, gfzz, gfwwt, gfzzt, &
+ gtnww, gtnzz, gtwz, gtww
+ real(default), public :: vev
+ complex(default), dimension(2), public :: &
+ gncneu, gnclep, gncup, gncdwn
+ real(default), public :: a4, a5, a6, a7, a10
+ complex(default), public :: ig1a, ig1z, rg5a, rg5z, &
+ ig1pkpg4a, ig1pkpg4z, ig1pkmg4a, ig1pkmg4z, &
+ ig1mkpg4a, ig1mkpg4z, ig1mkmg4a, ig1mkmg4z, &
+ ila, ilz, il5a, il5z, ik5a, ik5z, &
+ ialww0, ialww2, ialzw0, ialzw1, ialzz, &
+ alww0, alww2, alzw0, alzw1, alzz
+ real(default), public :: lam_reg
+ real(default), public :: fudge_higgs, fudge_km, w_res
+ real(default), dimension(1:12), public :: gkm
+ real(default), dimension(1:5), public :: mkm, wkm
+ complex(default), public :: ghvva
+ integer, public :: part_r
+ integer, private :: i
+
+ public :: import_from_whizard, model_update_alpha_s
+
+contains
+ subroutine import_from_whizard (par_array)
+ real(default), dimension(64), 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) :: xi0
+ real(default) :: xipm
+ real(default) :: a4
+ real(default) :: a5
+ real(default) :: a6
+ real(default) :: a7
+ real(default) :: a10
+ real(default) :: mkm_s
+ real(default) :: mkm_p
+ real(default) :: mkm_r
+ real(default) :: mkm_f
+ real(default) :: mkm_t
+ real(default) :: gkm_s
+ real(default) :: gkm_p
+ real(default) :: gkm_r
+ real(default) :: gkm_f
+ real(default) :: gkm_t
+ real(default) :: wkm_s
+ real(default) :: wkm_p
+ real(default) :: wkm_r
+ real(default) :: wkm_f
+ real(default) :: wkm_t
+ real(default) :: g1a
+ real(default) :: g1z
+ real(default) :: g4a
+ real(default) :: g4z
+ real(default) :: g5a
+ real(default) :: g5z
+ real(default) :: ka
+ real(default) :: kz
+ real(default) :: la
+ real(default) :: lz
+ real(default) :: k5a
+ real(default) :: k5z
+ real(default) :: l5a
+ real(default) :: l5z
+ real(default) :: mreg
+ real(default) :: fkm
+ real(default) :: wres
+ real(default) :: gkm_st
+ real(default) :: gkm_pt
+ real(default) :: gkm_rt
+ real(default) :: gkm_ft
+ real(default) :: gkm_tt
+ real(default) :: fmixed
+ real(default) :: fkappa
+ 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%xi0 = par_array(15)
+ par%xipm = par_array(16)
+ par%a4 = par_array(17)
+ par%a5 = par_array(18)
+ par%a6 = par_array(19)
+ par%a7 = par_array(20)
+ par%a10 = par_array(21)
+ par%mkm_s = par_array(22)
+ par%mkm_p = par_array(23)
+ par%mkm_r = par_array(24)
+ par%mkm_f = par_array(25)
+ par%mkm_t = par_array(26)
+ par%gkm_s = par_array(27)
+ par%gkm_p = par_array(28)
+ par%gkm_r = par_array(29)
+ par%gkm_f = par_array(30)
+ par%gkm_t = par_array(31)
+ par%wkm_s = par_array(32)
+ par%wkm_p = par_array(33)
+ par%wkm_r = par_array(34)
+ par%wkm_f = par_array(35)
+ par%wkm_t = par_array(36)
+ par%g1a = par_array(37)
+ par%g1z = par_array(38)
+ par%g4a = par_array(39)
+ par%g4z = par_array(40)
+ par%g5a = par_array(41)
+ par%g5z = par_array(42)
+ par%ka = par_array(43)
+ par%kz = par_array(44)
+ par%la = par_array(45)
+ par%lz = par_array(46)
+ par%k5a = par_array(47)
+ par%k5z = par_array(48)
+ par%l5a = par_array(49)
+ par%l5z = par_array(50)
+ par%mreg = par_array(51)
+ par%fkm = par_array(52)
+ par%wres = par_array(53)
+ par%gkm_st = par_array(54)
+ par%gkm_pt = par_array(55)
+ par%gkm_rt = par_array(56)
+ par%gkm_ft = par_array(57)
+ par%gkm_tt = par_array(58)
+ par%fmixed = par_array(59)
+ par%fkappa = par_array(60)
+ par%v = par_array(61)
+ par%cw = par_array(62)
+ par%sw = par_array(63)
+ par%ee = par_array(64)
+ mass(1:27) = 0
+ width(1:27) = 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(26) = par%xi0 * mass(23)
+ width(26) = 0
+ mass(27) = par%xipm * mass(24)
+ width(27) = 0
+ mass(45) = par%mkm_s
+ width(45) = par%wkm_s
+ mass(46) = par%mkm_p
+ width(46) = par%wkm_p
+ mass(47) = par%mkm_p
+ width(47) = par%wkm_p
+ mass(48) = par%mkm_p
+ width(48) = par%wkm_p
+ mass(52) = par%mkm_f
+ width(52) = par%wkm_f
+ mass(53) = par%mkm_t
+ width(53) = par%wkm_t
+ mass(54) = par%mkm_t
+ width(54) = par%wkm_t
+ mass(55) = par%mkm_t
+ width(55) = par%wkm_t
+ mkm(1) = par%mkm_s
+ mkm(2) = par%mkm_p
+ mkm(3) = par%mkm_r
+ mkm(4) = par%mkm_f
+ mkm(5) = par%mkm_t
+ gkm(1) = par%gkm_s
+ gkm(2) = par%gkm_p
+ gkm(3) = par%gkm_r
+ gkm(4) = par%gkm_f
+ gkm(5) = par%gkm_t
+ gkm(6) = par%gkm_st
+ gkm(7) = par%gkm_pt
+ gkm(8) = par%gkm_rt
+ gkm(9) = par%gkm_ft
+ gkm(10) = par%gkm_tt
+ gkm(11) = par%fmixed
+ gkm(12) = par%fkappa
+ wkm(1) = par%wkm_s
+ wkm(2) = par%wkm_p
+ wkm(3) = par%wkm_r
+ wkm(4) = par%wkm_f
+ wkm(5) = par%wkm_t
+ e = par%ee
+ sinthw = par%sw
+ sin2thw = par%sw**2
+ costhw = par%cw
+ qelep = - 1
+ qeup = 2.0_default / 3.0_default
+ qedwn = - 1.0_default / 3.0_default
+ g = e / sinthw
+ w_res = par%wres
+ vev = par%v
+ 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
+ gsww = gkm(1) * mass(24) * g
+ gszz = gkm(1) * mass(23) * g / costhw
+ gswwt = gkm(6) * g**3 / mass(24) / (16.0 * PI)
+ gszzt = gkm(6) * g**3 / costhw**3 / mass(23) /(16.0 * PI)
+ gpnww = - gkm(2) * mass(24) * g / 2 / sqrt(3.0_default)
+ gpnzz = gkm(2) * mass(23) * g / costhw / sqrt(3.0_default)
+ gpwz = gkm(2) * mass(23) * g / 2
+ gpww = gkm(2) * mass(24) * g / sqrt(2.0_default)
+ gfww = gkm(4) * mass(24) * g / 2
+ gfzz = gkm(4) * mass(23) * g / costhw / 2
+ gfwwt = gkm(9) * g**3 / mass(24) / (32.0 * PI)
+ gfzzt = gkm(9) * g**3 / costhw**3 / mass(23) /(32.0 * PI)
+ gtnww = - gkm(5) * mass(24) * g / 4 / sqrt(3.0_default)
+ gtnzz = gkm(5) * mass(23) * g / costhw / 2 / sqrt(3.0_default)
+ gtwz = gkm(5) * mass(23) * g / 4
+ gtww = gkm(5) * mass(24) * g / 2 / sqrt(2.0_default)
+ gssww = 0
+ gsszz = 0
+ part_r = 1
+ !!! 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
+ a4 = par%a4
+ a5 = par%a5
+ a6 = par%a6
+ a7 = par%a7
+ a10 = par%a10
+ lam_reg = par%mreg
+ fudge_higgs = 1
+ ghvva = 1
+ fudge_km = par%fkm
+ ig1a = iqw * par%g1a
+ ig1z = igzww * par%g1z
+ ig1pkpg4a = iqw * (par%g1a + par%ka + par%g4a) / 2
+ ig1pkpg4z = igzww * (par%g1z + par%kz + par%g4z) / 2
+ ig1pkmg4a = iqw * (par%g1a + par%ka - par%g4a) / 2
+ ig1pkmg4z = igzww * (par%g1z + par%kz - par%g4z) / 2
+ ig1mkpg4a = iqw * (par%g1a - par%ka + par%g4a) / 2
+ ig1mkpg4z = igzww * (par%g1z - par%kz + par%g4z) / 2
+ ig1mkmg4a = iqw * (par%g1a - par%ka - par%g4a) / 2
+ ig1mkmg4z = igzww * (par%g1z - par%kz - par%g4z) / 2
+ ila = iqw * par%la / (mass(24)*mass(24))
+ ilz = igzww * par%lz / (mass(24)*mass(24))
+ rg5a = qw * par%g5a
+ rg5z = gzww * par%g5z
+ ik5a = iqw * par%k5a
+ ik5z = igzww * par%k5z
+ il5a = iqw * par%l5a / (mass(24)*mass(24))
+ il5z = igzww * par%l5z / (mass(24)*mass(24))
+ alww0 = g**4 * (a4 + 2 * a5)
+ alww2 = g**4 * 2 * a4
+ alzw1 = g**4 / costhw**2 * (a4 + a6)
+ alzw0 = g**4 / costhw**2 * 2 * (a5 + a7)
+ alzz = g**4 / costhw**4 * 2 * (a4 + a5 + (a6+a7+a10)*2)
+ ialww0 = g**2 * sqrt (-cmplx(a4 + 2 * a5, kind=default))
+ ialww2 = g**2 * sqrt (-cmplx(2 * a4, kind=default))
+ ialzw1 = g**2 / costhw * sqrt (-cmplx(a4 + a6, kind=default))
+ ialzw0 = g**2 / costhw &
+ & * sqrt (-cmplx(2 * (a5 + a7), kind=default))
+ ialzz = g**2 / (costhw*costhw) &
+ & * sqrt (-cmplx(2 * (a4 + a5 + (a6+a7+a10)*2), &
+ & kind=default))
+ 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
+ !!! The Hgg coupling should not get a running alpha_s
+ end subroutine model_update_alpha_s
+end module parameters_alth
Index: trunk/src/models/parameters.NoH_rx.f90
===================================================================
--- trunk/src/models/parameters.NoH_rx.f90 (revision 0)
+++ trunk/src/models/parameters.NoH_rx.f90 (revision 5352)
@@ -0,0 +1,324 @@
+! $Id: parameters.NoH_rx.f90,v 1.0 2014/02/21 13:31:48 msekulla Exp $
+!
+! Copyright (C) 1999-2014 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>
+! Marco Sekulla <sekulla@physik.uni-siegen.de> (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_noh_rx
+ use kinds
+ use constants
+ implicit none
+ private
+
+ real(default), dimension(27), 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, &
+ iqw, igzww, igwww, gw4, gzzww, gazww, gaaww
+ real(default), public :: vev
+ complex(default), dimension(2), public :: &
+ gncneu, gnclep, gncup, gncdwn
+ real(default), public :: a4, a5, a6, a7, a10
+ complex(default), public :: ig1a, ig1z, rg5a, rg5z, &
+ ig1pkpg4a, ig1pkpg4z, ig1pkmg4a, ig1pkmg4z, &
+ ig1mkpg4a, ig1mkpg4z, ig1mkmg4a, ig1mkmg4z, &
+ ila, ilz, il5a, il5z, ik5a, ik5z, &
+ ialww0, ialww2, ialzw0, ialzw1, ialzz, &
+ alww0, alww2, alzw0, alzw1, alzz
+ real(default), public :: lam_reg
+ real(default), public :: fudge_higgs, fudge_km, w_res
+ real(default), dimension(1:12), public :: gkm
+ real(default), dimension(1:5), public :: mkm, wkm
+ complex(default), public :: ghvva
+ integer, public :: part_r
+ integer, private :: i
+
+ public :: import_from_whizard, model_update_alpha_s
+
+contains
+ subroutine import_from_whizard (par_array)
+ real(default), dimension(64), 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) :: xi0
+ real(default) :: xipm
+ real(default) :: a4
+ real(default) :: a5
+ real(default) :: a6
+ real(default) :: a7
+ real(default) :: a10
+ real(default) :: mkm_s
+ real(default) :: mkm_p
+ real(default) :: mkm_r
+ real(default) :: mkm_f
+ real(default) :: mkm_t
+ real(default) :: gkm_s
+ real(default) :: gkm_p
+ real(default) :: gkm_r
+ real(default) :: gkm_f
+ real(default) :: gkm_t
+ real(default) :: wkm_s
+ real(default) :: wkm_p
+ real(default) :: wkm_r
+ real(default) :: wkm_f
+ real(default) :: wkm_t
+ real(default) :: g1a
+ real(default) :: g1z
+ real(default) :: g4a
+ real(default) :: g4z
+ real(default) :: g5a
+ real(default) :: g5z
+ real(default) :: ka
+ real(default) :: kz
+ real(default) :: la
+ real(default) :: lz
+ real(default) :: k5a
+ real(default) :: k5z
+ real(default) :: l5a
+ real(default) :: l5z
+ real(default) :: mreg
+ real(default) :: fkm
+ real(default) :: wres
+ real(default) :: gkm_st
+ real(default) :: gkm_pt
+ real(default) :: gkm_rt
+ real(default) :: gkm_ft
+ real(default) :: gkm_tt
+ real(default) :: fmixed
+ real(default) :: fkappa
+ 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%xi0 = par_array(15)
+ par%xipm = par_array(16)
+ par%a4 = par_array(17)
+ par%a5 = par_array(18)
+ par%a6 = par_array(19)
+ par%a7 = par_array(20)
+ par%a10 = par_array(21)
+ par%mkm_s = par_array(22)
+ par%mkm_p = par_array(23)
+ par%mkm_r = par_array(24)
+ par%mkm_f = par_array(25)
+ par%mkm_t = par_array(26)
+ par%gkm_s = par_array(27)
+ par%gkm_p = par_array(28)
+ par%gkm_r = par_array(29)
+ par%gkm_f = par_array(30)
+ par%gkm_t = par_array(31)
+ par%wkm_s = par_array(32)
+ par%wkm_p = par_array(33)
+ par%wkm_r = par_array(34)
+ par%wkm_f = par_array(35)
+ par%wkm_t = par_array(36)
+ par%g1a = par_array(37)
+ par%g1z = par_array(38)
+ par%g4a = par_array(39)
+ par%g4z = par_array(40)
+ par%g5a = par_array(41)
+ par%g5z = par_array(42)
+ par%ka = par_array(43)
+ par%kz = par_array(44)
+ par%la = par_array(45)
+ par%lz = par_array(46)
+ par%k5a = par_array(47)
+ par%k5z = par_array(48)
+ par%l5a = par_array(49)
+ par%l5z = par_array(50)
+ par%mreg = par_array(51)
+ par%fkm = par_array(52)
+ par%wres = par_array(53)
+ par%gkm_st = par_array(54)
+ par%gkm_pt = par_array(55)
+ par%gkm_rt = par_array(56)
+ par%gkm_ft = par_array(57)
+ par%gkm_tt = par_array(58)
+ par%fmixed = par_array(59)
+ par%fkappa = par_array(60)
+ par%v = par_array(61)
+ par%cw = par_array(62)
+ par%sw = par_array(63)
+ par%ee = par_array(64)
+ mass(1:27) = 0
+ width(1:27) = 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(26) = par%xi0 * mass(23)
+ width(26) = 0
+ mass(27) = par%xipm * mass(24)
+ width(27) = 0
+ mkm(1) = par%mkm_s
+ mkm(2) = par%mkm_p
+ mkm(3) = par%mkm_r
+ mkm(4) = par%mkm_f
+ mkm(5) = par%mkm_t
+ gkm(1) = par%gkm_s
+ gkm(2) = par%gkm_p
+ gkm(3) = par%gkm_r
+ gkm(4) = par%gkm_f
+ gkm(5) = par%gkm_t
+ wkm(1) = par%wkm_s
+ wkm(2) = par%wkm_p
+ wkm(3) = par%wkm_r
+ wkm(4) = par%wkm_f
+ wkm(5) = par%wkm_t
+ gkm(6) = par%gkm_st
+ gkm(7) = par%gkm_pt
+ gkm(8) = par%gkm_rt
+ gkm(9) = par%gkm_ft
+ gkm(10) = par%gkm_tt
+ gkm(11) = par%fmixed
+ gkm(12) = par%fkappa
+ vev = par%v
+ e = par%ee
+ sinthw = par%sw
+ sin2thw = par%sw**2
+ costhw = par%cw
+ 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
+ part_r = 0
+ !!! 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
+ a4 = par%a4
+ a5 = par%a5
+ a6 = par%a6
+ a7 = par%a7
+ a10 = par%a10
+ lam_reg = par%mreg
+ fudge_higgs = 1
+ ghvva = 1
+ fudge_km = par%fkm
+ w_res = par%wres
+ ig1a = iqw * par%g1a
+ ig1z = igzww * par%g1z
+ ig1pkpg4a = iqw * (par%g1a + par%ka + par%g4a) / 2
+ ig1pkpg4z = igzww * (par%g1z + par%kz + par%g4z) / 2
+ ig1pkmg4a = iqw * (par%g1a + par%ka - par%g4a) / 2
+ ig1pkmg4z = igzww * (par%g1z + par%kz - par%g4z) / 2
+ ig1mkpg4a = iqw * (par%g1a - par%ka + par%g4a) / 2
+ ig1mkpg4z = igzww * (par%g1z - par%kz + par%g4z) / 2
+ ig1mkmg4a = iqw * (par%g1a - par%ka - par%g4a) / 2
+ ig1mkmg4z = igzww * (par%g1z - par%kz - par%g4z) / 2
+ ila = iqw * par%la / (mass(24)*mass(24))
+ ilz = igzww * par%lz / (mass(24)*mass(24))
+ rg5a = qw * par%g5a
+ rg5z = gzww * par%g5z
+ ik5a = iqw * par%k5a
+ ik5z = igzww * par%k5z
+ il5a = iqw * par%l5a / (mass(24)*mass(24))
+ il5z = igzww * par%l5z / (mass(24)*mass(24))
+ alww0 = g**4 * (a4 + 2 * a5)
+ alww2 = g**4 * 2 * a4
+ alzw1 = g**4 / costhw**2 * (a4 + a6)
+ alzw0 = g**4 / costhw**2 * 2 * (a5 + a7)
+ alzz = g**4 / costhw**4 * 2 * (a4 + a5 + (a6+a7+a10)*2)
+ ialww0 = g**2 * sqrt (-cmplx(a4 + 2 * a5, kind=default))
+ ialww2 = g**2 * sqrt (-cmplx(2 * a4, kind=default))
+ ialzw1 = g**2 / costhw * sqrt (-cmplx(a4 + a6, kind=default))
+ ialzw0 = g**2 / costhw &
+ & * sqrt (-cmplx(2 * (a5 + a7), kind=default))
+ ialzz = g**2 / (costhw*costhw) &
+ & * sqrt (-cmplx(2 * (a4 + a5 + (a6+a7+a10)*2), &
+ & kind=default))
+ 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
+ !!! The Hgg coupling should not get a running alpha_s
+ end subroutine model_update_alpha_s
+end module parameters_noh_rx
Index: trunk/src/omega/src/omega_NoH_rx.ml
===================================================================
--- trunk/src/omega/src/omega_NoH_rx.ml (revision 0)
+++ trunk/src/omega/src/omega_NoH_rx.ml (revision 5352)
@@ -0,0 +1,36 @@
+(* $Id: omega_NoH_rx.ml 4926 2014-03-21 12:35:06 msekulla $
+
+ Copyright (C) 1999-2014 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>
+ Marco Sekulla <sekulla@physik.uni-siegen.de>
+
+ 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 O = Omega.Make(Fusion.Mixed23)(Targets.Fortran)
+ (Modellib_NoH.NoH(Modellib_NoH.NoH_k_matrix))
+let _ = O.main ()
+
+(*i
+ * Local Variables:
+ * mode:caml
+ * indent-tabs-mode:nil
+ * page-delimiter:"^(\\* .*\n"
+ * End:
+i*)
Index: trunk/src/omega/src/omega_AltH.ml
===================================================================
--- trunk/src/omega/src/omega_AltH.ml (revision 0)
+++ trunk/src/omega/src/omega_AltH.ml (revision 5352)
@@ -0,0 +1,36 @@
+(* $Id: omega_AltH.ml 3670 2014-03-21 17:07:07Z msekulla $
+
+ Copyright (C) 1999-2014 by
+
+ Wolfgang Kilian <kilian@physik.uni-siegen.de>
+ Thorsten Ohl <ohl@physik.uni-wuerzburg.de>
+ Juergen Reuter <juergen.reuter@physik.uni-freiburg.de>
+ with contributions from
+ Marco Sekulla <sekulla@physik.uni-siegen.de>
+ Christian Speckner <cnspeckn@googlemail.com>
+
+ 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 O = Omega.Make(Fusion.Mixed23)(Targets.Fortran)
+ (Modellib_NoH.AltH(Modellib_NoH.NoH_k_matrix))
+let _ = O.main ()
+
+(*i
+ * Local Variables:
+ * mode:caml
+ * indent-tabs-mode:nil
+ * page-delimiter:"^(\\* .*\n"
+ * End:
+i*)
Index: trunk/share/models/AltH.mdl
===================================================================
--- trunk/share/models/AltH.mdl (revision 0)
+++ trunk/share/models/AltH.mdl (revision 5352)
@@ -0,0 +1,348 @@
+########################################################################
+# Standard Model (unitary gauge) with anomalous couplings,
+# CKM matrix and K-Matrix unitarization
+# Default Higgs mass is very large
+model "AltH"
+
+# Independent parameters
+### DO NOT CHANGE THE ORDER OF THESE PARAMETERS
+parameter GF = 1.16639E-5 # Fermi constant
+parameter mZ = 91.1882 # Z-boson mass
+parameter mW = 80.419 # W-boson mass
+parameter alphas = 0.1178 # Strong coupling constant (Z point)
+parameter me = 0.000511 # electron mass
+parameter mmu = 0.1057 # muon mass
+parameter mtau = 1.777 # tau-lepton mass
+parameter ms = 0.12 # s-quark mass
+parameter mc = 1.25 # c-quark mass
+parameter mb = 4.2 # b-quark mass
+parameter mtop = 174 # t-quark mass
+parameter wtop = 1.523 # t-quark width
+parameter wZ = 2.443 # Z-boson width
+parameter wW = 2.049 # W-boson width
+parameter xi0 = 0 # R_xi parameter for Z-boson
+parameter xipm = 0 # R_xi parameter for W-boson
+parameter a4 = 0 # Coefficient of L4
+parameter a5 = 0 # Coefficient of L5
+parameter a6 = 0 # Coefficient of L6
+parameter a7 = 0 # Coefficient of L7
+parameter a10 = 0 # Coefficient of L10
+parameter mkm_s = 1.E10 # Mass for scalar isoscalar resonance
+parameter mkm_p = 1.E10 # Mass for scalar isotensor resonance
+parameter mkm_r = 1.E10 # Mass for vector isovector resonance
+parameter mkm_f = 1.E10 # Mass for tensor isotensor resonance
+parameter mkm_t = 1.E10 # Mass for tensor isotensor resonance
+parameter gkm_s = 0 # Coupling for scalar isoscalar resonance
+parameter gkm_p = 0 # Coupling for scalar isotensor resonance
+parameter gkm_r = 0 # Coupling for vector isovector resonance
+parameter gkm_f = 0 # Coupling for tensor isoscalar resonance
+parameter gkm_t = 0 # Coupling for tensor isotensor resonance
+parameter wkm_s = 1.E-4 # Width for scalar isoscalar resonance
+parameter wkm_p = 1.E-4 # Width for scalar isotensor resonance
+parameter wkm_r = 1.E-4 # Width for vector isovector resonance
+parameter wkm_f = 1.E-4 # Width for tensor isotensor resonance
+parameter wkm_t = 1.E-4 # Width for tensor isotensor resonance
+parameter g1a = 1 # anomalous TGC
+parameter g1z = 1 # anomalous TGC
+parameter g4a = 0 # anomalous TGC
+parameter g4z = 0 # anomalous TGC
+parameter g5a = 0 # anomalous TGC
+parameter g5z = 0 # anomalous TGC
+parameter ka = 1 # anomalous TGC
+parameter kz = 1 # anomalous TGC
+parameter la = 0 # anomalous TGC
+parameter lz = 0 # anomalous TGC
+parameter k5a = 0 # anomalous TGC
+parameter k5z = 0 # anomalous TGC
+parameter l5a = 0 # anomalous TGC
+parameter l5z = 0 # anomalous TGC
+parameter mreg = 1000 # regulator mass for one-loop correction
+parameter fkm = 1 # fudge factor to switch K matrix unitarization
+ # for resonances on or off
+parameter wres = 0 # switch for resonance widths without K matrix
+parameter gkm_st = 0 # Transversal coupling for scalar isoscalar resonance
+parameter gkm_pt = 0 # Transversal coupling for scalar isotensor resonance
+parameter gkm_rt = 0 # Transversal Coupling for vector isovector resonance
+parameter gkm_ft = 0 # Transversal coupling for tensor isoscalar resonance
+parameter gkm_tt = 0 # Coupling for tensor isotensor resonance
+parameter fmixed = 1 # Coupling for mixed terms of transversal and longitudinal couplings
+parameter fkappa = 0 # kappa Term for tensor isoscalar resonance (experimental)
+
+# Dependent parameters
+derived v = 1 / sqrt (sqrt (2.) * GF) # v (Higgs vev)
+derived cw = mW / mZ # cos(theta-W)
+derived sw = sqrt (1-cw**2) # sin(theta-W)
+derived ee = 2 * sw * mW / v # em-coupling (GF scheme)
+derived alpha_em_i = 4 * pi / ee**2 # inverse fine structure const
+
+########################################################################
+# Particle content
+
+# The quarks
+particle D_QUARK 1 parton
+ spin 1/2 charge -1/3 isospin -1/2 color 3
+ name d down
+ anti dbar D "d~"
+ tex_anti "\bar{d}"
+particle U_QUARK 2 parton
+ spin 1/2 charge 2/3 isospin 1/2 color 3
+ name u up
+ anti ubar U "u~"
+ tex_anti "\bar{u}"
+particle S_QUARK 3 like D_QUARK
+ name s strange
+ anti sbar S "s~"
+ tex_anti "\bar{s}"
+ mass ms
+particle C_QUARK 4 like U_QUARK
+ name c charm
+ anti cbar C "c~"
+ tex_anti "\bar{c}"
+ mass mc
+particle B_QUARK 5 like D_QUARK
+ name b bottom
+ anti bbar B "b~"
+ tex_anti "\bar{b}"
+ mass mb
+particle T_QUARK 6 like U_QUARK
+ name t top
+ anti tbar T "t~"
+ tex_anti "\bar{t}"
+ mass mtop width wtop
+
+# The leptons
+particle E_LEPTON 11
+ spin 1/2 charge -1 isospin -1/2
+ name "e-" e1 electron e
+ anti "e+" E1 positron
+ tex_name "e^-"
+ tex_anti "e^+"
+ mass me
+particle E_NEUTRINO 12 left
+ spin 1/2 isospin 1/2
+ name nue n1 "nu_e" ve "e-neutrino"
+ anti nuebar N1 "ve~"
+ tex_name "\nu_e"
+ tex_anti "\bar\nu_e"
+particle MU_LEPTON 13 like E_LEPTON
+ name "mu-" e2 mu muon
+ anti "mu+" E2
+ tex_name "\mu^-"
+ tex_anti "\mu^+"
+ mass mmu
+particle MU_NEUTRINO 14 like E_NEUTRINO
+ name numu "nu_mu" n2 vm "mu-neutrino"
+ anti numubar N2 "vm~"
+ tex_name "\nu_\mu"
+ tex_anti "\bar\nu_\mu"
+particle TAU_LEPTON 15 like E_LEPTON
+ name "tau-" e3 tau "ta-" tauon
+ anti "tau+" E3 "ta+"
+ tex_name "\tau^-"
+ tex_anti "\tau^+"
+ mass mtau
+particle TAU_NEUTRINO 16 like E_NEUTRINO
+ name nutau "nu_tau" n3 vt "tau_neutrino"
+ anti nutaubar N3 "vt~"
+ tex_name "\nu_\tau"
+ tex_anti "\bar\nu_\tau"
+
+# The vector bosons
+particle GLUON 21 parton gauge
+ spin 1 color 8
+ name gl g G gluon
+particle PHOTON 22 gauge
+ spin 1
+ name A gamma photon
+ tex_name "\gamma"
+particle Z_BOSON 23 gauge
+ spin 1
+ name Z
+ mass mZ width wZ
+particle W_BOSON 24 gauge
+ spin 1 charge 1
+ name "W+" Wp
+ anti "W-" Wm
+ tex_name "W^+"
+ tex_anti "W^-"
+ mass mW width wW
+
+# Resonance Scalar IsoScalar
+particle AltH_SIGMA 45
+ spin 0
+ name Rsigma
+ tex_name "\sigma"
+ mass mkm_s width wkm_s
+
+# Resonance Scalar IsoTensor
+particle AltH_PHI_N 46
+ spin 0
+ name "Rphi0" Rphin
+ tex_name "\phi^0"
+ mass mkm_p width wkm_p
+
+particle AltH_PHI_P 47
+ spin 0
+ name "Rphi+" Rphip
+ anti "Rphi-" Rphim
+ tex_name "\phi^+"
+ tex_anti "\phi^-"
+ mass mkm_p width wkm_p
+
+particle AltH_PHI_PP 48
+ spin 0
+ name "Rphi++" Rphipp
+ anti "Rphi--" Rphimm
+ tex_name "\phi^{++}"
+ tex_anti "\phi^{--}"
+ mass mkm_p width wkm_p
+
+#Resonance Tensor Isoscalar
+particle AltH_f 52
+ spin 2
+ name Rf
+ tex_name "f"
+ mass mkm_f width wkm_f
+
+# Resonance Tensor IsoTensor
+particle AltH_T_N 53
+ spin 2
+ name "Rt0" Rtn
+ tex_name "t^0"
+ mass mkm_t width wkm_t
+
+particle AltH_T_P 54
+ spin 2
+ name "Rt+" Rtp
+ anti "Rt-" Rtm
+ tex_name "t^+"
+ tex_anti "t^-"
+ mass mkm_t width wkm_t
+
+particle AltH_T_PP 55
+ spin 2
+ name "Rt++" Rtpp
+ anti "Rt--" Rtmm
+ tex_name "t^{++}"
+ tex_anti "t^{--}"
+ mass mkm_t width wkm_t
+
+# Hadrons
+particle PROTON 2212
+ spin 1/2 charge 1
+ name p "p+"
+ anti pbar "p-"
+
+# Beam remnants for proton colliders
+particle HADRON_REMNANT 90
+ name hr
+ tex_name "had_r"
+particle HADRON_REMNANT_SINGLET 91
+ name hr1
+ tex_name "had_r^{(1)}"
+particle HADRON_REMNANT_TRIPLET 92
+ color 3
+ name hr3
+ tex_name "had_r^{(3)}"
+ anti hr3bar
+ tex_name "had_r^{(\bar 3)}"
+particle HADRON_REMNANT_OCTET 93
+ color 8
+ name hr8
+ tex_name "had_r^{(8)}"
+
+########################################################################
+# Vertices of the Standard model with anomalous couplings
+# In graphs with identical structure, the first vertex is kept for phase space,
+# therefore, lighter particles come before heavier ones.
+#
+# Note: The electric charge has opposite sign compared to SM.mdl
+
+
+# QED
+vertex D d A
+vertex U u A
+vertex S s A
+vertex C c A
+vertex B b A
+vertex T t A
+
+vertex E1 e1 A
+vertex E2 e2 A
+vertex E3 e3 A
+
+# QCD
+vertex G G G
+
+vertex D d G
+vertex U u G
+vertex S s G
+vertex C c G
+vertex B b G
+vertex T t G
+
+# Neutral currents
+vertex D d Z
+vertex U u Z
+vertex S s Z
+vertex C c Z
+vertex B b Z
+vertex T t Z
+
+vertex E1 e1 Z
+vertex E2 e2 Z
+vertex E3 e3 Z
+vertex N1 n1 Z
+vertex N2 n2 Z
+vertex N3 n3 Z
+
+# Charged currents
+vertex U d Wp
+vertex C s Wp
+vertex T b Wp
+vertex D u Wm
+vertex S c Wm
+vertex B t Wm
+
+vertex N1 e1 Wp
+vertex N2 e2 Wp
+vertex N3 e3 Wp
+vertex E1 n1 Wm
+vertex E2 n2 Wm
+vertex E3 n3 Wm
+
+# Vector-boson self-interactions
+vertex A Wm Wp
+vertex Wm Wp Z
+
+vertex A A Wm Wp
+vertex A Wm Wp Z
+vertex Wm Wm Wp Wp
+vertex Wm Wp Z Z
+vertex Z Z Z Z
+
+# Sigma - vector boson
+vertex Rsigma Wp Wm
+vertex Rsigma Z Z
+vertex Rsigma Rsigma Z Z
+vertex Rsigma Rsigma Wp Wm
+
+# Phi - vector boson
+vertex Rphin Wp Wm
+vertex Rphin Z Z
+vertex Rphip Wm Z
+vertex Rphipp Wm Wm
+vertex Rphim Wp Z
+vertex Rphimm Wp Wp
+
+# Tensor f - vector boson
+vertex Rf Wp Wm
+vertex Rf Z Z
+
+# Tensor t - vector boson
+vertex Rtn Wp Wm
+vertex Rtn Z Z
+vertex Rtp Wm Z
+vertex Rtpp Wm Wm
+vertex Rtm Wp Z
+vertex Rtmm Wp Wp
Index: trunk/share/models/NoH_rx.mdl
===================================================================
--- trunk/share/models/NoH_rx.mdl (revision 0)
+++ trunk/share/models/NoH_rx.mdl (revision 5352)
@@ -0,0 +1,266 @@
+########################################################################
+# Higgsless Model (unitary gauge) with anomalous couplings,
+# CKM matrix and K-Matrix unitarization
+model "NoH_rx"
+
+# Independent parameters
+### DO NOT CHANGE THE ORDER OF THESE PARAMETERS
+parameter GF = 1.16639E-5 # Fermi constant
+parameter mZ = 91.1882 # Z-boson mass
+parameter mW = 80.419 # W-boson mass
+parameter alphas = 0.1178 # Strong coupling constant (Z point)
+parameter me = 0.000511 # electron mass
+parameter mmu = 0.1057 # muon mass
+parameter mtau = 1.777 # tau-lepton mass
+parameter ms = 0.12 # s-quark mass
+parameter mc = 1.25 # c-quark mass
+parameter mb = 4.2 # b-quark mass
+parameter mtop = 174 # t-quark mass
+parameter wtop = 1.523 # t-quark width
+parameter wZ = 2.443 # Z-boson width
+parameter wW = 2.049 # W-boson width
+parameter xi0 = 0 # R_xi parameter for Z-boson
+parameter xipm = 0 # R_xi parameter for W-boson
+parameter a4 = 0 # Coefficient of L4
+parameter a5 = 0 # Coefficient of L5
+parameter a6 = 0 # Coefficient of L6
+parameter a7 = 0 # Coefficient of L7
+parameter a10 = 0 # Coefficient of L10
+parameter mkm_s = 1.E10 # Mass for scalar isoscalar resonance
+parameter mkm_p = 1.E10 # Mass for scalar isotensor resonance
+parameter mkm_r = 1.E10 # Mass for vector isovector resonance
+parameter mkm_f = 1.E10 # Mass for tensor isotensor resonance
+parameter mkm_t = 1.E10 # Mass for tensor isotensor resonance
+parameter gkm_s = 0 # Coupling for scalar isoscalar resonance
+parameter gkm_p = 0 # Coupling for scalar isotensor resonance
+parameter gkm_r = 0 # Coupling for vector isovector resonance
+parameter gkm_f = 0 # Coupling for tensor isoscalar resonance
+parameter gkm_t = 0 # Coupling for tensor isotensor resonance
+parameter wkm_s = 0 # Width for scalar isoscalar resonance
+parameter wkm_p = 0 # Width for scalar isotensor resonance
+parameter wkm_r = 0 # Width for vector isovector resonance
+parameter wkm_f = 0 # Width for tensor isotensor resonance
+parameter wkm_t = 0 # Width for tensor isotensor resonance
+parameter g1a = 1 # anomalous TGC
+parameter g1z = 1 # anomalous TGC
+parameter g4a = 0 # anomalous TGC
+parameter g4z = 0 # anomalous TGC
+parameter g5a = 0 # anomalous TGC
+parameter g5z = 0 # anomalous TGC
+parameter ka = 1 # anomalous TGC
+parameter kz = 1 # anomalous TGC
+parameter la = 0 # anomalous TGC
+parameter lz = 0 # anomalous TGC
+parameter k5a = 0 # anomalous TGC
+parameter k5z = 0 # anomalous TGC
+parameter l5a = 0 # anomalous TGC
+parameter l5z = 0 # anomalous TGC
+parameter mreg = 1000 # regulator mass for one-loop correction
+parameter fkm = 1 # fudge factor to switch K matrix unitarization
+ # for resonances on or off
+parameter wres = 1 # switch for resonance widths without K matrix
+parameter gkm_st = 0 # Transversal coupling for scalar isoscalar resonance
+parameter gkm_pt = 0 # Transversal coupling for scalar isotensor resonance
+parameter gkm_rt = 0 # Transversal Coupling for vector isovector resonance
+parameter gkm_ft = 0 # Transversal coupling for tensor isoscalar resonance
+parameter gkm_tt = 0 # Coupling for tensor isotensor resonance
+parameter fmixed = 1 # Coupling for mixed terms of transversal and longitudinal couplings
+parameter fkappa = 0 # kappa Term for tensor isoscalar resonance (experimental)
+
+# Dependent parameters
+derived v = 1 / sqrt (sqrt (2.) * GF) # v (Higgs vev)
+derived cw = mW / mZ # cos(theta-W)
+derived sw = sqrt (1-cw**2) # sin(theta-W)
+derived ee = 2 * sw * mW / v # em-coupling (GF scheme)
+derived alpha_em_i = 4 * pi / ee**2 # inverse fine structure const
+
+########################################################################
+# Particle content
+
+# The quarks
+particle D_QUARK 1 parton
+ spin 1/2 charge -1/3 isospin -1/2 color 3
+ name d down
+ anti dbar D "d~"
+ tex_anti "\bar{d}"
+particle U_QUARK 2 parton
+ spin 1/2 charge 2/3 isospin 1/2 color 3
+ name u up
+ anti ubar U "u~"
+ tex_anti "\bar{u}"
+particle S_QUARK 3 like D_QUARK
+ name s strange
+ anti sbar S "s~"
+ tex_anti "\bar{s}"
+ mass ms
+particle C_QUARK 4 like U_QUARK
+ name c charm
+ anti cbar C "c~"
+ tex_anti "\bar{c}"
+ mass mc
+particle B_QUARK 5 like D_QUARK
+ name b bottom
+ anti bbar B "b~"
+ tex_anti "\bar{b}"
+ mass mb
+particle T_QUARK 6 like U_QUARK
+ name t top
+ anti tbar T "t~"
+ tex_anti "\bar{t}"
+ mass mtop width wtop
+
+# The leptons
+particle E_LEPTON 11
+ spin 1/2 charge -1 isospin -1/2
+ name "e-" e1 electron e
+ anti "e+" E1 positron
+ tex_name "e^-"
+ tex_anti "e^+"
+ mass me
+particle E_NEUTRINO 12 left
+ spin 1/2 isospin 1/2
+ name nue n1 "nu_e" ve "e-neutrino"
+ anti nuebar N1 "ve~"
+ tex_name "\nu_e"
+ tex_anti "\bar\nu_e"
+particle MU_LEPTON 13 like E_LEPTON
+ name "mu-" e2 mu muon
+ anti "mu+" E2
+ tex_name "\mu^-"
+ tex_anti "\mu^+"
+ mass mmu
+particle MU_NEUTRINO 14 like E_NEUTRINO
+ name numu "nu_mu" n2 vm "mu-neutrino"
+ anti numubar N2 "vm~"
+ tex_name "\nu_\mu"
+ tex_anti "\bar\nu_\mu"
+particle TAU_LEPTON 15 like E_LEPTON
+ name "tau-" e3 tau "ta-" tauon
+ anti "tau+" E3 "ta+"
+ tex_name "\tau^-"
+ tex_anti "\tau^+"
+ mass mtau
+particle TAU_NEUTRINO 16 like E_NEUTRINO
+ name nutau "nu_tau" n3 vt "tau_neutrino"
+ anti nutaubar N3 "vt~"
+ tex_name "\nu_\tau"
+ tex_anti "\bar\nu_\tau"
+
+# The vector bosons
+particle GLUON 21 parton gauge
+ spin 1 color 8
+ name gl g G gluon
+particle PHOTON 22 gauge
+ spin 1
+ name A gamma photon
+ tex_name "\gamma"
+particle Z_BOSON 23 gauge
+ spin 1
+ name Z
+ mass mZ width wZ
+particle W_BOSON 24 gauge
+ spin 1 charge 1
+ name "W+" Wp
+ anti "W-" Wm
+ tex_name "W^+"
+ tex_anti "W^-"
+ mass mW width wW
+
+# Hadrons
+particle PROTON 2212
+ spin 1/2 charge 1
+ name p "p+"
+ anti pbar "p-"
+
+# Beam remnants for proton colliders
+particle HADRON_REMNANT 90
+ name hr
+ tex_name "had_r"
+particle HADRON_REMNANT_SINGLET 91
+ name hr1
+ tex_name "had_r^{(1)}"
+particle HADRON_REMNANT_TRIPLET 92
+ color 3
+ name hr3
+ tex_name "had_r^{(3)}"
+ anti hr3bar
+ tex_name "had_r^{(\bar 3)}"
+particle HADRON_REMNANT_OCTET 93
+ color 8
+ name hr8
+ tex_name "had_r^{(8)}"
+
+########################################################################
+# Vertices of the Standard model with anomalous couplings
+# In graphs with identical structure, the first vertex is kept for phase space,
+# therefore, lighter particles come before heavier ones.
+#
+# Note: The electric charge has opposite sign compared to SM.mdl
+
+
+# QED
+vertex D d A
+vertex U u A
+vertex S s A
+vertex C c A
+vertex B b A
+vertex T t A
+
+vertex E1 e1 A
+vertex E2 e2 A
+vertex E3 e3 A
+
+# QCD
+vertex G G G
+
+vertex D d G
+vertex U u G
+vertex S s G
+vertex C c G
+vertex B b G
+vertex T t G
+
+# Neutral currents
+vertex D d Z
+vertex U u Z
+vertex S s Z
+vertex C c Z
+vertex B b Z
+vertex T t Z
+
+vertex E1 e1 Z
+vertex E2 e2 Z
+vertex E3 e3 Z
+vertex N1 n1 Z
+vertex N2 n2 Z
+vertex N3 n3 Z
+
+# Charged currents
+vertex U d Wp
+vertex C s Wp
+vertex T b Wp
+vertex D u Wm
+vertex S c Wm
+vertex B t Wm
+
+vertex N1 e1 Wp
+vertex N2 e2 Wp
+vertex N3 e3 Wp
+vertex E1 n1 Wm
+vertex E2 n2 Wm
+vertex E3 n3 Wm
+
+# Vector-boson self-interactions
+vertex A Wm Wp
+vertex Wm Wp Z
+
+vertex A A Wm Wp
+vertex A Wm Wp Z
+vertex Wm Wm Wp Wp
+vertex Wm Wp Z Z
+vertex Z Z Z Z
+
+
+
+
+
File Metadata
Details
Attached
Mime Type
text/x-diff
Expires
Tue, Nov 19, 7:41 PM (1 d, 9 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
3805884
Default Alt Text
(44 KB)
Attached To
rWHIZARDSVN whizardsvn
Event Timeline
Log In to Comment