Confusing output of FermionQ

[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

[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

[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