Merge pull request #255 from fusion-energy/add_photon_cut

cutting photons with negligible contribution to dose to speed up the r2s simulation
fusion-energy · Nov 20, 2023 · a31bf92 · a31bf92
2 parents 2beafd0 + f69bf1b
commit a31bf92
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 61 deletions.
diff --git a/tasks/task_07_CSG_cell_tally_spectra/5_neutron_energy_loss_plot.py b/tasks/task_07_CSG_cell_tally_spectra/5_neutron_energy_loss_plot.py
diff --git a/tasks/task_11_CSG_shut_down_dose_tallies/1_cell_based_shut_down_dose_rate_example.py b/tasks/task_11_CSG_shut_down_dose_tallies/1_cell_based_shut_down_dose_rate_example.py
@@ -179,7 +179,10 @@
             # gets the activated material using the material id
             activated_mat = results[i_cool].get_material(str(material_id))
             # gets the energy and probabilities for the 
-            energy = activated_mat.get_decay_photon_energy()
+            energy = activated_mat.get_decay_photon_energy(
+                clip_tolerance = 1e-6,  # cuts out a small fraction of the very low energy (and hence negligible dose contribution) photons
+                units = 'Bq',
+            )
             strength = energy.integral()
 
             if strength > 0.:  # only makes sources for