Explicit numerical instability handling #145

[jacanchaplais]

• replaced phi where pt is low with nan (and raised NumericalStabilityWarning)
• corrected _root_diff_two_squares() for case when both inputs are negative (in turn fixing mass_t)
• rapidity division by zero explicitly handled, converting nan values into appropriately signed infinities