Update update_packets.cc

artis-mcrt · Sep 1, 2024 · 7a66ef1 · 7a66ef1
1 parent 6aaee0e
commit 7a66ef1
Showing 1 changed file with 1 addition and 5 deletions.
diff --git a/update_packets.cc b/update_packets.cc
@@ -83,11 +83,7 @@ void do_nonthermal_predeposit(Packet &pkt, const int nts, const double t2) {
     // for endot independent of energy, the next line is trival (for E dependent endot, an integral would be needed)
 
     const double t_enzero = ts + (particle_en / endot);  // time at which zero energy is reached
-    if (t_enzero > t2) {
-      en_deposited = pkt.e_cmf * (t2 - ts) / (particle_en / endot);
-    } else {
-      en_deposited = pkt.e_cmf * (t_enzero - ts) / (particle_en / endot);
-    }
+    en_deposited = pkt.e_cmf * (std::min(t2, t_enzero) - ts) / (particle_en / endot);
 
     // A discrete absorption event should occur somewhere along the
     // continuous track from initial kinetic energy to zero KE.