calculate beam profile stats

ARM-DOE · Jan 19, 2024 · aa0a041 · aa0a041
1 parent f6aabbf
commit aa0a041
Showing 1 changed file with 5 additions and 0 deletions.
diff --git a/tests/test_gfit.py b/tests/test_gfit.py
@@ -29,3 +29,8 @@ def test_psds():
         my_binary['ScatDiaBC50'].values[ScatRejectKey == 0]) < 1000.
     my_psds = pysp2.util.process_psds(my_binary, my_hk, my_ini)
     np.testing.assert_almost_equal(my_psds['NumConcIncan'].max(), 0.95805343)
+    np.testing.assert_almost_equal(my_binary['DeadtimeRelativeBias'].mean(), -0.00023515)
+    coeff, beam_profile = pysp2.util.beam_shape(my_binary,beam_position_from='peak maximum',Globals=pysp2.util.DMTGlobals())
+
+
+