[Rivet] Fwd: FourMomentum mass

Andy Buckley andy.buckley at ed.ac.uk
Thu Nov 11 11:42:49 GMT 2010


On 11/11/10 00:45, Hendrik Hoeth wrote:
> Thus spake Hendrik Hoeth (hendrik.hoeth at cern.ch):
> 
>> Guess we should catch that. In fact I thought we had already fixed
>> this several months ago.
> 
> Looking at the code, we didn't. We inserted an assert(), but left it to
> the caller to check if he gets a reasonable result.
> 
> Can someone remind me why we decided against doing at least something
> like
> 
> if (Rivet::isZero(mass2())
>   return 0.0;
> else
>   return sqrt(mass2());
> 
> This should catch the photons but will still blow up in case of a real
> problem.

No idea why we didn't do that. Because we were dealing with a different
failure mode, I guess. I'm happy for this to go in: can you make the
change if no-one else objects?

Andy

-- 
Dr Andy Buckley
SUPA Advanced Research Fellow
Particle Physics Experiment Group, University of Edinburgh


More information about the Rivet mailing list