Page 1 of 1
Confusing output of FermionQ
Posted: 14. Oct 2019, 14:45
by AlexanderVoigt
Dear SARAH developers,
in FlexibleSUSY we stumbled upon a confusing output of the FermionQ function:
Code: Select all
Needs["SARAH`"];
Start["SM"];
FermionQ[Fe] (* -> False *)
FermionQ[Fv] (* -> True *)
Is this a bug or are we misinterpreting the output of FermionQ?
Best regards from Flensburg
Alex
Re: Confusing output of FermionQ
Posted: 15. Oct 2019, 19:26
by FStaub
Hi Alex,
I guess, somewhere in between. Without checking (because I don't have Mathematica anymore):
make sure that SA`CurrentStates is set correct at this stage (i.e. EWSB). It is usually set when calling other functions to perform some calculations, but just after running Start is might be GaugeES (and the components of Fe are not known there).
Cheers,
Florian
Re: Confusing output of FermionQ
Posted: 16. Oct 2019, 15:43
by AlexanderVoigt
Hi Florian,
many thanks for the hint, this was in deed the problem!
Code: Select all
Needs["SARAH`"];
Start["SM"];
?SA`CurrentStates (* -> GaugeES *)
FermionQ[Fe] (* -> False *)
FermionQ[Fv] (* -> True *)
SA`CurrentStates = SARAH`EWSB
FermionQ[Fe] (* -> True *)
FermionQ[Fv] (* -> True *)
Best regards,
Alex