diff --git a/.idea/other.xml b/.idea/other.xml new file mode 100644 index 0000000000000000000000000000000000000000..a708ec78108573c549112f9a71e489890930c832 --- /dev/null +++ b/.idea/other.xml @@ -0,0 +1,6 @@ +<?xml version="1.0" encoding="UTF-8"?> +<project version="4"> + <component name="PySciProjectComponent"> + <option name="PY_SCI_VIEW_SUGGESTED" value="true" /> + </component> +</project> \ No newline at end of file diff --git a/sixsutils.py b/sixsutils.py index 48dec75743f4bf879e7d2ea2e5e97a9c63a39b93..661ee0cc05ebda4e87c535f2d3a2f7231f100c5f 100644 --- a/sixsutils.py +++ b/sixsutils.py @@ -57,7 +57,10 @@ def birks_correction(energy): # Birks coefficient and parameter estimates taken :return: Corrected energy in Mev """ kbirks = 0.68E-3 - return energy * (1 - scs.hyp2f1(1, 1 / 0.678, 1 + 1 / 0.678, -np.power(energy, 0.678) / (95 * kbirks))) + if energy > 100: + return energy + else: + return energy * (1 - scs.hyp2f1(1, 1 / 0.678, 1 + 1 / 0.678, -np.power(energy, 0.678) / (95 * kbirks))) def core_energy_to_64(energy):