|
[Rivet] Fwd: FourMomentum massAndy Buckley andy.buckley at ed.ac.ukThu 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 |