That would indeed be a very nice possibility, but I am not sure you can always find an analytical solution. Martin's solution only works for the MSSM. As soon as you complicate the model a little you start finding difficulties to solve the tadpole equations (unless you do it for squared mass terms)....

Solving the tadpole equations for parameters other than squared mass terms is usually problematic. The reason is that Mathematica may not be able to find an analytical solution to the equations. I am not sure there is a way around this...

Great!

That's strange. Have you previously created a SM SPheno module with SARAH and later copied it (with cp -r) into SPheno?

Dear Robin, I believe that the files that you downloaded from https://sarah.hepforge.org/trac/wiki/NMSSM is pretty much outdated. My recommendation is to regenerate it using the latest SARAH version, that will produce an output for SPheno that works perfectly with SPheno 4.0.3.

Hi and sorry for the delay in my response. Now I see your approach. Since you already identified the loops leading to the discrepancy and know how to find their analytical expressions in the Fortran code (as you write, they are in the FlavorKit_LFV_modelname.f90 file), I would compare these to your ...

For my understanding: how do you compare your analytical results to those obtained by FlavorKit?

Hi Saiyad, sometimes there are discrepancies for individual loop diagrams that cancel out in the sum. However, I understand that is not your case. You also get a different result in the sum. Is this correct? If so, I can think of two possible reasons: 1) Please check whether the involved loop functi...

As a far as I know, introducing "fictitious" mixings is usually a bad idea in SARAH. I tried to do it several times but finally realized that it leads to many problems. My guess is that that's precisely the origin of your issues. However, you mention that not including the mixing leads to problems a...

I just tried this implementation of the Scotogenic model with SARAH-4.13.0 and SPheno-4.0.3. The SPheno code is generated correctly, and I only get some warnings... Non-zero result for {etI,etp,conj(Hp)} vertex which violates charge. (This might just be a problem with simplifying the vertex, but cou...