Page 1 of 1

New SPheno Bug in SARAH 4.10.2

Posted: 19. Mar 2017, 05:16
by JPEllis
Hi,

I updated to SARAH 4.10.2 and SPheno 4.0.2 today, and there appears to be a regression. The BranchingRatios file no longer compiles do to the following lines:

Code: Select all

gPVP = 0._dp 
gTVP = 0._dp 
BRVP = 0._dp 
Call VPTwoBodyDecay(-1,DeltaM,BMC,BMN,FMLd,FMLe,FMLU,FLxd,FLxe,FMRd,FMRe,             & 
& FMRu,FRxd,FRxe,FMv,MFd,MFd2,MFe,MFe2,MFu,MFu2,MFv,MFv2,MFxd,MFxd2,MFxe,MFxe2,          & 
& MSA,MSA2,MSC,MSC2,MSS,MSS2,MVWm,MVWm2,MVZ,MVZ2,MVZp,MVZp2,MVZpp,MVZpp2,SMC,            & 
& SMP,SMS,ThetaPZ,g1,g2,g3,gU4,gU5,Ye,Yd,Yxd,Yv,Yu,Yxu,Y1d,Y1l,Y2,LAMdd,LAMdduu,         & 
& LAMd1,LAMd2,LAMduud,LAMuu,LAMu1,LAMu2,LAM11,LAM12,LAM22,kap,MUd,MUu,MU1,               & 
& MU2,gPVP,gTVP,BRVP)

! FIXME
! Do i1=1,1
! 0. =Sum(0.(i1,:))
! If (0..Gt.0._dp) 0.(i1,: ) =0.(i1,:)/0.
! End Do
 
Here, VP is the photon and I suspect this bug might be due do to the photon having its width set to 0 in the particles.m file.

Re: New SPheno Bug in SARAH 4.10.2

Posted: 19. Mar 2017, 12:26
by FStaub
Hi,

I'm confused why there is code to generate the photon decays at all. Is the particle marked by

Code: Select all

Description -> "Photon"
in particles.m? Anyway, I'll change it that decays for massless states are never included.

Cheers,
Florian

Re: New SPheno Bug in SARAH 4.10.2

Posted: 19. Mar 2017, 12:45
by JPEllis
I indeed have

Code: Select all

    Description -> "Photon",
in the definition of the symbol for the photon. Probably safer to not have massless states decay in general.