Hi Florian,
I have implemented a model extending the MSSM with three generations of right-neutrino superfields.
Because of the presence of a Majorana term (M N^c N^c) in the superpotential, and the corresponding
soft-term, there are six CP-even and six CP-odd sneutrino mass eigenstates in the sneutrino sector.
With an older version of SARAH, I managed to generate the SPheno files successfully. However,
since we needed a loop-induced decay, I tried to run the model file with the latest version of SARAH
(4.12.2). While I get no (serious) warning/error in CheckModel[], running MakeSPheno[] gives an error
message related to the (two-loop) Higgs mass computation; please find the error message in the
attached.
Although I am able to run the SPheno model file, thus generated, with SPheno (v4), after a little edit, but
the Higgs mass is lower (for a few benchmarks) compared to the SPheno files generated by an older version
of SARAH. Kindly let me know if you have any comments.
Further, is it possible to implement chargino_1--> neutralino_1 + pion(s) in the SPheno files generated by SARAH?
I find that such decays are included in the MSSM SPheno version.
Many thanks for your help in advance.
right-sneutrino
-
MDGoodsell
- Posts: 17
- Joined: 25. Aug 2017, 08:41
Re: right-sneutrino
Dear Arin,
This error looks identical to one that I have seen elsewhere, where the problem was the Mathematica version. Are you using a version prior to 10? In which case the latest routines will not work, I am sorry to say. If you use a more recent Mathematica version, however, then it should work without problems. If this is impossible for you, then you can always turn off the two loop routines by defining
Include2LoopCorrections=False
before typing MakeSPheno.
Possibly in future I'll make this backwards compatible, but in general this is difficult to do!
This error looks identical to one that I have seen elsewhere, where the problem was the Mathematica version. Are you using a version prior to 10? In which case the latest routines will not work, I am sorry to say. If you use a more recent Mathematica version, however, then it should work without problems. If this is impossible for you, then you can always turn off the two loop routines by defining
Include2LoopCorrections=False
before typing MakeSPheno.
Possibly in future I'll make this backwards compatible, but in general this is difficult to do!
Re: right-sneutrino
Thank you for your prompt response. Yes, indeed I am using Mathematica 9.
I will run the file with Mathematica 10.
Could you also kindly address my second query in the previous post, if possible :
" Further, is it possible to implement chargino_1--> neutralino_1 + pion(s) in the SPheno files generated by SARAH?
I find that such decays are included in the MSSM SPheno version. "
Thanks once again.
I will run the file with Mathematica 10.
Could you also kindly address my second query in the previous post, if possible :
" Further, is it possible to implement chargino_1--> neutralino_1 + pion(s) in the SPheno files generated by SARAH?
I find that such decays are included in the MSSM SPheno version. "
Thanks once again.
Re: right-sneutrino
Hi,
1) If you used before a SARAH version before 4.10.0, then the changes in the HIggs mass are most likely due to the changes explained here: 1703.03267
2) You would need to hack the code to include the pion decay. That's nothing what is planned to be added on SARAH side. You find the routine (FermionToFermionPi) in SPheno/src/DecayFunctions.f90. Of course, it would need some modifications in the generated code to get this implemented consistently. An easier approach might to calculate this decay within Mathematica for each point using the information provided in the spc file.
Cheers,
Florian
1) If you used before a SARAH version before 4.10.0, then the changes in the HIggs mass are most likely due to the changes explained here: 1703.03267
2) You would need to hack the code to include the pion decay. That's nothing what is planned to be added on SARAH side. You find the routine (FermionToFermionPi) in SPheno/src/DecayFunctions.f90. Of course, it would need some modifications in the generated code to get this implemented consistently. An easier approach might to calculate this decay within Mathematica for each point using the information provided in the spc file.
Cheers,
Florian