[Rivet] Fwd: FourMomentum mass

Hendrik Hoeth hendrik.hoeth at cern.ch
Thu Nov 11 00:45:50 GMT 2010


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.

Cheers,

   Hendrik

-- 
One doesn't discover new lands without consenting to lose
sight of the shore for a very long time.    -- André Gide


More information about the Rivet mailing list