Error when setting low particle cuts for em showers
An error that terminates the shower occurs when setting low particle cuts (10 MeV) and electromagnetic interaction. The error is not present when setting high cuts (10 GeV). The error message received was:
terminate called after throwing an instance of 'boost::wrapexcept<boost::math::evaluation_error>' what(): Error in function boost::math::tools::bisect<double>: No change of sign in boost::math::tools::bisect, either there is no root to find, or there are multiple roots in the interval (f(min) = -9.8713280286356291). Aborted (core dumped)
Previously several messages of:
No node was found where the function value was greater than zero. Either there are no continuous losses for the cut, or something went wrong.
coming from crosssection/CrossSectionDEDX/AxisBuilderDEDX.h are shown.
The proposal_new_installation_and_bugfixes branch was used (latest commit 60f00801) and PROPOSAL 7.0.2 is running.
I leave here the examples used and the ouput to reproduce it. Vertical EAS (modified to start with an electron and have electromagnetic interaction) was run with ./bin/vertical_EAS 0 1 1000 1 and em shower with ./bin/em_shower 1000. em_shower.cpp vertical_EAS.cpp output_em_shower output_vertical_EAS