[12gevmoller_sim] [JeffersonLab/remoll] b51040: Bugfix to #233 where generated pion energy fails a...
Wouter Deconinck
noreply at github.com
Tue Mar 19 15:57:44 EDT 2019
Branch: refs/heads/bugfix-test-pion-single-threaded
Home: https://gcc01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2FJeffersonLab%2Fremoll&data=02%7C01%7C12gevmoller_sim%40jlab.org%7Cc2086b937c9b4780a9c008d6aca52d1f%7Cb4d7ee1f4fb34f0690372b5b522042ab%7C1%7C1%7C636886222773087191&sdata=%2B2IT4kl6L%2F%2FB%2FcQvrVZ57XK3KQE3GGf6xEqgLVPJVTw%3D&reserved=0
Commit: b5104065c368fc9fb1d97eb07f8f500d9cf15f1e
https://gcc01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2FJeffersonLab%2Fremoll%2Fcommit%2Fb5104065c368fc9fb1d97eb07f8f500d9cf15f1e&data=02%7C01%7C12gevmoller_sim%40jlab.org%7Cc2086b937c9b4780a9c008d6aca52d1f%7Cb4d7ee1f4fb34f0690372b5b522042ab%7C1%7C1%7C636886222773087191&sdata=6aUX76TCbAuTvl0EHnMzgPAGp75dG2FgW4oH%2B70R1vA%3D&reserved=0
Author: Wouter Deconinck <wdconinc at gmail.com>
Date: 2019-03-19 (Tue, 19 Mar 2019)
Changed paths:
M src/remollGenPion.cc
Log Message:
-----------
Bugfix to #233 where generated pion energy fails assertion
The problem was that the beam can radiate so much that the energy is
lower than the minimum energy with which we wish to generate. In that
case G4RandFlat::shoot(emin, emax) receives a larger emin than emax,
and interprets this as a backwards interval. This subsequently fails
the e > emin test.
In this case, when emax is the radiated energy, and emax < emin, we
now set emin = emax and generate in that interval. This also changes
the phase space allocated to this event to zero (so properly weighted
outputs should not see this event anymore).
More information about the 12gevmoller_sim
mailing list