Excellent! I had defined them as "Dirac Neutrinos", but once changing them to "Neutrinos", and fixing some issues with the ElectronYukawa variable appearing every here and there, the code compiled.

Many thanks, let's hope that I find no problems when doing the tests. :-O

Ok, it didn't work. :-D Now SARAH does recognize that the nu_tau is massless. However, I now get MNeutrino and MNeutrino2 variables within LoopCouplings, FlavorKit_QFV and SPheno subroutines. These variables are generally undefined, except in FlavorKit_QFV, where they appear as: Complex(dp),Intent(i...

Right, but the problem there is that the reconstructed mass matrix ends up being Majorana, and I need it to be Dirac.

What about not using the Flavour option, and separating the lmu states into two families from the start?

Sure! I did have to split generations to make only the nu_e interact. This is how I did it (I only show the parts involving leptons): FermionFields[[2]] = {le, 1, {veL0, eL0}, -1/2, 2, 1}; FermionFields[[6]] = {v1c, 1, vR1, 0, 1, 1}; FermionFields[[7]] = {v2, 1, conj[vR2], 0, 1, 1}; FermionFields[[8...

Hi Florian, I am building a model with two heavy neutrinos N_i (which do not decouple). Here, the SM neutrinos are massless, the heavy neutrinos are Dirac, but there do exist Yukawa couplings between the SM nus and the N_i (such that they mix). As I want the final state fermions to be Dirac after EW...

Hi Florian, I am working with a non-decoupled version of Seesaw1, and noticed that the mass matrices were not being correctly reproduced in SPheno. I later connected this with the definition of nuR, I had to define it with an additional conj: SuperFields[[8]] = {v, 3, conj[vR], 0, 1, 1, RpM}; Out of...

Hi all, I have spotted a potential problem in RunSM_MSSM. Let me first summarize what I understand the procedure does. The procedure runs SUSY and SM parameters from M_SUSY to Q_in. This can be controlled with flags 62 and 63. With flag 62 = 1, both SUSY and SM parameters run from M_SUSY to Q_in. Wi...