Update gammapkt.cc

artis-mcrt · Sep 16, 2024 · 7d0a27e · 7d0a27e
1 parent c353346
commit 7d0a27e
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/gammapkt.cc b/gammapkt.cc
@@ -696,9 +696,12 @@ void pair_prod(Packet &pkt) {
   }
 
   if (rng_uniform() > prob_gamma) {
+    // Convert it to an e-minus packet - actually it could be positron EK too, but this works
+    // for consistency with compton_scatter.
+    pkt.type = TYPE_NONTHERMAL_PREDEPOSIT_BETAMINUS;
     // Convert it to an e-minus or positron kinetic energy packet
+    // pkt.type = (rng_uniform() > 0.5) ? TYPE_NONTHERMAL_PREDEPOSIT_BETAMINUS : TYPE_NONTHERMAL_PREDEPOSIT_BETAPLUS;
     // nu_cmf stays the same as the gamma energy becomes the kinetic energy of the electron
-    pkt.type = (rng_uniform() > 0.5) ? TYPE_NONTHERMAL_PREDEPOSIT_BETAMINUS : TYPE_NONTHERMAL_PREDEPOSIT_BETAPLUS;
     pkt.absorptiontype = -5;
     stats::increment(stats::COUNTER_NT_STAT_FROM_GAMMA);
   } else {