Page MenuHomeHEPForge
Files F7879182

No OneTemporary
Actions

View Options

Index: trunk/src/models/parameters.AltH.f90
===================================================================
--- trunk/src/models/parameters.AltH.f90 (revision 5354)
+++ trunk/src/models/parameters.AltH.f90 (revision 5355)
@@ -1,362 +1,362 @@
-! $Id: parameters.VBS.f90, 2013/09/23 13:11:48 msekulla Exp $
+! $Id: parameters.AltH.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/Makefile.am
===================================================================
--- trunk/src/models/Makefile.am (revision 5354)
+++ trunk/src/models/Makefile.am (revision 5355)
@@ -1,316 +1,320 @@
## Makefile.am -- Makefile for model-parameter modules in WHIZARD
##
## Process this file with automake to produce Makefile.in
#
# 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>
# 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.
#
########################################################################
## Build the threeshl bundle required for the 3SHL models.
SUBDIRS = \
threeshl_bundle
## Build WHIZARD/O'Mega model interface library
## These modules transfer WHIZARD's parameters to the matrix-element code
## as generated by O'Mega
noinst_LTLIBRARIES = libmodels.la
libmodels_la_SOURCES = \
parameters.2HDM.f90 \
parameters.2HDM_CKM.f90 \
parameters.AltH.f90 \
parameters.GravTest.f90 \
parameters.Littlest.f90 \
parameters.Littlest_Eta.f90 \
parameters.Littlest_Tpar.f90 \
parameters.MSSM.f90 \
parameters.MSSM_CKM.f90 \
parameters.MSSM_Grav.f90 \
parameters.MSSM_Hgg.f90 \
parameters.NMSSM.f90 \
parameters.NMSSM_CKM.f90 \
parameters.NMSSM_Hgg.f90 \
parameters.NoH_rx.f90 \
parameters.QED.f90 \
parameters.QCD.f90 \
parameters.PSSSM.f90 \
parameters.Simplest.f90 \
parameters.Simplest_univ.f90 \
parameters.SM.f90 \
parameters.SM_ac.f90 \
parameters.SM_ac_CKM.f90 \
parameters.SM_CKM.f90 \
parameters.SM_top.f90 \
parameters.SM_top_anom.f90 \
parameters.SM_Higgs.f90 \
parameters.SM_rx.f90 \
parameters.Template.f90 \
parameters.UED.f90 \
parameters.SSC.f90 \
parameters.Xdim.f90 \
parameters.Zprime.f90 \
parameters.Threeshl.f90 \
parameters.Threeshl_nohf.f90 \
parameters.Test.f90
## Install the modules used by generated matrix element code
execmoddir = $(pkglibdir)/mod/models
nodist_execmod_HEADERS = \
parameters_2hdm.$(FC_MODULE_EXT) \
parameters_2hdm_ckm.$(FC_MODULE_EXT) \
+ parameters_alth.$(FC_MODULE_EXT) \
parameters_gravtest.$(FC_MODULE_EXT) \
parameters_littlest.$(FC_MODULE_EXT) \
parameters_littlest_eta.$(FC_MODULE_EXT) \
parameters_littlest_tpar.$(FC_MODULE_EXT) \
parameters_mssm.$(FC_MODULE_EXT) \
parameters_mssm_ckm.$(FC_MODULE_EXT) \
parameters_mssm_hgg.$(FC_MODULE_EXT) \
parameters_mssm_grav.$(FC_MODULE_EXT) \
parameters_nmssm.$(FC_MODULE_EXT) \
parameters_nmssm_ckm.$(FC_MODULE_EXT) \
parameters_nmssm_hgg.$(FC_MODULE_EXT) \
parameters_noh_rx.$(FC_MODULE_EXT) \
parameters_qed.$(FC_MODULE_EXT) \
parameters_qcd.$(FC_MODULE_EXT) \
parameters_psssm.$(FC_MODULE_EXT) \
parameters_simplest.$(FC_MODULE_EXT) \
parameters_simplest_univ.$(FC_MODULE_EXT) \
parameters_sm.$(FC_MODULE_EXT) \
parameters_sm_ac.$(FC_MODULE_EXT) \
parameters_sm_ac_ckm.$(FC_MODULE_EXT) \
parameters_sm_ckm.$(FC_MODULE_EXT) \
parameters_sm_top.$(FC_MODULE_EXT) \
parameters_sm_top_anom.$(FC_MODULE_EXT) \
parameters_sm_higgs.$(FC_MODULE_EXT) \
parameters_sm_rx.$(FC_MODULE_EXT) \
parameters_template.$(FC_MODULE_EXT) \
parameters_ued.$(FC_MODULE_EXT) \
parameters_ssc.$(FC_MODULE_EXT) \
parameters_xdim.$(FC_MODULE_EXT) \
parameters_zprime.$(FC_MODULE_EXT) \
parameters_threeshl.$(FC_MODULE_EXT) \
parameters_threeshl_nohf.$(FC_MODULE_EXT) \
parameters_test.$(FC_MODULE_EXT)
## Build WHIZARD auxiliary dynamic libraries for model access
## These modules handle initialization operations that are not accounted
## for by definitions in the model files (external parameters).
modellibdir = $(pkglibdir)/models
modellib_LTLIBRARIES = \
external.Threeshl.la \
external.Threeshl_nohf.la \
external.Test.la
external_Threeshl_la_SOURCES = external.Threeshl.f90
external_Threeshl_la_LDFLAGS = -module
external_Threeshl_nohf_la_SOURCES = external.Threeshl_nohf.f90
external_Threeshl_nohf_la_LDFLAGS = -module
external_Test_la_SOURCES = external.Test.f90
external_Test_la_LDFLAGS = -module
## Collect the 3SHL code.
external_Threeshl_la_LIBADD = threeshl_bundle/libthreeshl_bundle.la
libmodels_la_LIBADD = threeshl_bundle/libthreeshl_bundle.la
external_Threeshl_nohf_la_LIBADD = threeshl_bundle/libthreeshl_bundle.la
## Switch off optimization if requested
## Optimization is useless here, it just consumes compile time.
if !OPTIMIZATION_FOR_PARFILES
NO_OPT=-O0
endif
.f90.o:
$(FCCOMPILE) $(NO_OPT) -c -o $@ $<
.f90.obj:
$(FCCOMPILE) $(NO_OPT) -c -o $@ `$(CYGPATH_W) '$<'`
.f90.lo:
$(LTFCCOMPILE) $(NO_OPT) -c -o $@ $<
## The mismatch of filenames and module names requires a lot of repetition:
parameters_2hdm.$(FC_MODULE_EXT): parameters.2HDM.lo
@:
parameters_2hdm_ckm.$(FC_MODULE_EXT): parameters.2HDM_CKM.lo
@:
+parameters_alth.$(FC_MODULE_EXT): parameters.2HDM_CKM.lo
+ @:
+
parameters_gravtest.$(FC_MODULE_EXT): parameters.GravTest.lo
@:
parameters_littlest.$(FC_MODULE_EXT): parameters.Littlest.lo
@:
parameters_littlest_eta.$(FC_MODULE_EXT): parameters.Littlest_Eta.lo
@:
parameters_littlest_tpar.$(FC_MODULE_EXT): parameters.Littlest_Tpar.lo
@:
parameters_mssm.$(FC_MODULE_EXT): parameters.MSSM.lo
@:
parameters_mssm_4.$(FC_MODULE_EXT): parameters.MSSM_4.lo
@:
parameters_mssm_ckm.$(FC_MODULE_EXT): parameters.MSSM_CKM.lo
@:
parameters_mssm_grav.$(FC_MODULE_EXT): parameters.MSSM_Grav.lo
@:
parameters_mssm_hgg.$(FC_MODULE_EXT): parameters.MSSM_Hgg.lo
@:
parameters_nmssm.$(FC_MODULE_EXT): parameters.NMSSM.lo
@:
parameters_nmssm_ckm.$(FC_MODULE_EXT): parameters.NMSSM_CKM.lo
@:
parameters_nmssm_hgg.$(FC_MODULE_EXT): parameters.NMSSM_Hgg.lo
@:
parameters_noh_rx.$(FC_MODULE_EXT): parameters.NoH_rx.lo
@:
parameters_psssm.$(FC_MODULE_EXT): parameters.PSSSM.lo
@:
parameters_qcd.$(FC_MODULE_EXT): parameters.QCD.lo
@:
parameters_qed.$(FC_MODULE_EXT): parameters.QED.lo
@:
parameters_sm.$(FC_MODULE_EXT): parameters.SM.lo
@:
parameters_sm_ckm.$(FC_MODULE_EXT): parameters.SM_CKM.lo
@:
parameters_sm_ac.$(FC_MODULE_EXT): parameters.SM_ac.lo
@:
parameters_sm_ac_ckm.$(FC_MODULE_EXT): parameters.SM_ac_CKM.lo
@:
parameters_sm_top.$(FC_MODULE_EXT): parameters.SM_top.lo
@:
parameters_sm_top_anom.$(FC_MODULE_EXT): parameters.SM_top_anom.lo
@:
parameters_sm_higgs.$(FC_MODULE_EXT): parameters.SM_Higgs.lo
@:
parameters_sm_rx.$(FC_MODULE_EXT): parameters.SM_rx.lo
@:
parameters_simplest.$(FC_MODULE_EXT): parameters.Simplest.lo
@:
parameters_simplest_univ.$(FC_MODULE_EXT): parameters.Simplest_univ.lo
@:
parameters_template.$(FC_MODULE_EXT): parameters.Template.lo
@:
parameters_test.$(FC_MODULE_EXT): parameters.Test.lo
@:
parameters_threeshl.$(FC_MODULE_EXT): parameters.Threeshl.lo
@:
parameters_threeshl_nohf.$(FC_MODULE_EXT): parameters.Threeshl_nohf.lo
@:
parameters_ued.$(FC_MODULE_EXT): parameters.UED.lo
@:
parameters_ssc.$(FC_MODULE_EXT): parameters.SSC.lo
@:
parameters_xdim.$(FC_MODULE_EXT): parameters.Xdim.lo
@:
parameters_zprime.$(FC_MODULE_EXT): parameters.Zprime.lo
@:
########################################################################
## Dependency on O'Mega modules
auxdir = $(top_builddir)/src/misc
KINDS_MOD = kinds.$(FC_MODULE_EXT)
CONSTANTS_MOD = constants.$(FC_MODULE_EXT)
DIAGNOSTICS_MOD = diagnostics.$(FC_MODULE_EXT)
SMPHYSICS_MOD = sm_physics.$(FC_MODULE_EXT)
FILEUTILS_MOD = file_utils.$(FC_MODULE_EXT)
### To be changed when moved in an external physics library
$(libmodels_la_OBJECTS) $(external_Threeshl_la_OBJECTS): \
$(auxdir)/$(KINDS_MOD) \
$(auxdir)/$(CONSTANTS_MOD) \
$(auxdir)/$(SMPHYSICS_MOD) \
$(auxdir)/$(DIAGNOSTICS_MOD) \
$(auxdir)/$(FILEUTILS_MOD) \
threeshl_bundle/threeshl.$(FC_MODULE_EXT) \
threeshl_bundle/tglue.$(FC_MODULE_EXT)
AM_FCFLAGS = -I$(auxdir) -Ithreeshl_bundle
########################################################################
## Default Fortran compiler options
## Profiling
if FC_USE_PROFILING
AM_FCFLAGS += $(FCFLAGS_PROFILING)
endif
## OpenMP
if FC_USE_OPENMP
AM_FCFLAGS += $(FCFLAGS_OPENMP)
endif
########################################################################
## Threeshl and Threeshl_nohf are identical
parameters.Threeshl_nohf.f90: parameters.Threeshl.f90
sed -e 's/\(module parameters_threeshl\)/\1_nohf/' < $< > $@
external.Threeshl_nohf.f90: external.Threeshl.f90
cp $< $@
########################################################################
## Non-standard cleanup tasks
CLEANFILES = \
parameters.Threeshl_nohf.f90 \
external.Threeshl_nohf.f90
clean-local:
-rm -f *.$(FC_MODULE_EXT)
## Remove backup files
maintainer-clean-local:
-rm -f *~
Index: trunk/share/models/SSC.mdl
===================================================================
--- trunk/share/models/SSC.mdl (revision 5354)
+++ trunk/share/models/SSC.mdl (revision 5355)
@@ -1,366 +1,366 @@
########################################################################
# Standard Model (unitary gauge) with anomalous couplings,
# CKM matrix and K-Matrix unitarization
# Default Higgs mass is very large
model "SSC"
# 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 mH = 1.E10 # Higgs mass
+parameter mH = 125 # Higgs 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 wH = 0 # Higgs width
+parameter wH = 0.004143 # Higgs 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
# The Higgs
particle HIGGS 25
spin 0
name H h Higgs
mass mH width wH
# Resonance Scalar IsoScalar
particle SSC_SIGMA 45
spin 0
name Rsigma
tex_name "\sigma"
mass mkm_s width wkm_s
# Resonance Scalar IsoTensor
particle SSC_PHI_N 46
spin 0
name "Rphi0" Rphin
tex_name "\phi^0"
mass mkm_p width wkm_p
particle SSC_PHI_P 47
spin 0
name "Rphi+" Rphip
anti "Rphi-" Rphim
tex_name "\phi^+"
tex_anti "\phi^-"
mass mkm_p width wkm_p
particle SSC_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 SSC_f 52
spin 2
name Rf
tex_name "f"
mass mkm_f width wkm_f
# Resonance Tensor IsoTensor
particle SSC_T_N 53
spin 2
name "Rt0" Rtn
tex_name "t^0"
mass mkm_t width wkm_t
particle SSC_T_P 54
spin 2
name "Rt+" Rtp
anti "Rt-" Rtm
tex_name "t^+"
tex_anti "t^-"
mass mkm_t width wkm_t
particle SSC_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
# Higgs - vector boson
vertex H Wp Wm
vertex H Z Z
vertex H H Z Z
vertex H H Wp Wm
# 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
# Higgs self-interactions
vertex H H H
vertex H H H H

File Metadata

Mime Type
text/x-diff
Expires
Tue, Nov 19, 7:42 PM (1 d, 7 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
3805887
Default Alt Text
(32 KB)

Event Timeline