Page 1 of 1

2HDM model file in SARAH

Posted: 11. Oct 2016, 09:39
by MailQuestion
Hi Florian,
We were trying to use the default 2HDM type-II model file from the SARAH model database https://sarah.hepforge.org/trac/wiki/THDM-II, but for some reason, the UFO output file (UFO.tar.gz) given there doesn't work with MadGraph and gives a compilation error (please see attached) which we cannot decode. Hope you can help us.

Thanks.

Re: 2HDM model file in SARAH

Posted: 11. Oct 2016, 09:46
by FStaub
Hi,

I think the crucial line is:

Code: Select all

MDL_TW = ASIN(SQRT(DCMPLX(1.000000D+00-MDL_MWM__EXP__2   
Error: 'x' argument of 'asin' intrinsic at (1) must be REAL
The compiler complaints because Asin expects a real argument. I haven't seen this problem before, so it might just come from a picky compiler. Althought, I'm surprised because I used gfortran already myself.
Anyway, maybe it already helps if you change the following block parameters.py

Code: Select all

TW = 	 Parameter(name='TW', 
	 nature = 'internal', 
	 type = 'real', 
	 value = 'cmath.asin(cmath.sqrt(1 - MWm**2/MZ**2))', 
	 texname = 'TW') 
and use math.sqrt instead of cmath.sqrt.

Cheers,
Florian