now cutting out low energy photons

fusion-energy · Nov 20, 2023 · f69bf1b · f69bf1b
1 parent 2beafd0
commit f69bf1b
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