FeynHiggs Generated SLHA file cannot be used by MG5

Is this not a question for FeynHiggs author?

Not sure how I can help you on MG5 side.

Cheers,

Olivier

On 06 Jan 2015, at 08:31, Lingfeng Li <email address hidden> wrote:

> New question #260233 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/260233
>
> What I am doing is using FeynHiggs to generate spectrum and branching ratios in MSSM. So, when I try to use the "OutputSLHA" subroutine in the package, it generated a small SLHA file, which doesn't include information about decay and Branching Ratios (coming from subroutine FHCouplings in FeynHIggs). MadGraph calculator accepted this file and generated a param_card.dat.
>
> Needless to say, the file cannot pass the MG5, because the decay blocks are missing. So I use FeynHiggs again, this time using:
>
> call OutputSLHA(error,slhadata,key)
>
> where key= 1+2+4+8+16+32 (which means including decay information with the 'fast' flag)
>
> Now, there seems to be a small bug here, and the program gets in a infinite loop. I will explian below:
>
> in subroutine OutputSLHA, everything is fine before calling the subroutine SLHANewDecay. but in SLHANewDecay, it said:
>
>
> par_id = abs(parent_id)
> decay = 1
>
> do while( Decay(decay) .ne. invalid )
> next = Decay_Next(decay)
>
> if( abs(Decay_Id(decay)) .eq. par_id ) then
> do while( Decay(decay + next) .ne. invalid )
> Decay(decay) = Decay(decay + next)
> decay = decay + 1
> enddo
> goto 1
> endif
>
> decay = decay + next
> enddo
> ........
>
> I started with par_id =25. And function Decay(decay)=Decay(1)(the initial value)=0!
> Also, next= Decau_NextDecay(decay)=0 too!
> Then the if block has never triggered, and we go to decay=decay+next=decay+0=decay. This gives a infinite loop.
>
> I am really confused if I get wronged somehow or there is a loophole.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.