Problems with one loop decays

[PhilipDiessner]

Hey Florian,

I am having some rather strange issues with the one loop decay widths where some the channels seem not to appear
for me. I attach my SPheno output for the MSSM default using freshly downloaded SARAH 4.12.3 and SPheno 4.0.3.

Most notable are the decay of the top quark missing the main final state and the Higgs missing bbar.
Have you seen this before and/or any idea what I am doing wrong?

Cheers Philip

[FStaub]

Hi Philipp,

that sounds strange. Did you try to look at the debug output? Putting flag 1103 to one, you see the values of the different pieces. Maybe, that helps to locate the problem.

Cheers,
Florian

[PhilipDiessner]

It seems like something is going wrong with the loop contributions as they are larger than the tree-level ones,
but are not included because they are quite often negative.
The top decay looks e.g. like this

Code: Select all

 ------------------ Fu->Fd conj[VWm] -----------------------
           3           1
 TREE x TREE:   0.28218357849087300     
 TREE x WAVE:   0.34934785430167420     
 TREE x VERTEX:   -4.3353411475606718     
 TREE x (TREE+WAVE+VERTEX):   -3.7038097147681239     
           3           2
 TREE x TREE:    13.418818578938646     
 TREE x WAVE:    18.671175059019522     
 TREE x VERTEX:   -155.62547295809065     
 TREE x (TREE+WAVE+VERTEX):   -123.53547932013248     
           3           3
 TREE x TREE:    8188.0507180939558     
 TREE x WAVE:    13030.747761126197     
 TREE x VERTEX:   -48911.484168844581     
 TREE x (TREE+WAVE+VERTEX):   -27692.685689624424  
 

[FStaub]

Hi,

the big negative numbers are actually not a problem, because they stll contain the IR divergent pieces. However, in one of the last updates an "eq" got exchanged by "le" what caused the problem. The patch is attached.

Cheers,
Florian