From 5d98fdc1bd1b0516e49a05dd8f55febec587b931 Mon Sep 17 00:00:00 2001 From: tlatorre Date: Mon, 16 Nov 2020 09:02:20 -0600 Subject: calculate livetime based on pulse gt and 10 MHz clock This commit updates get_events() to calculate the livetime based on both the number of pulse gt events and the 10 MHz clock and to return it in a dictionary stored with the dataframe. I also update dm-search so that the results are now reported as a function of events/cm^3/s. Also updated radius cut to be the AV radius. --- utils/chi2 | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) (limited to 'utils/chi2') diff --git a/utils/chi2 b/utils/chi2 index 0326e1a..39db081 100755 --- a/utils/chi2 +++ b/utils/chi2 @@ -480,10 +480,10 @@ if __name__ == '__main__': ev_mc = ev_mc[~np.isnan(ev_mc.fmin)] muon_mc = muon_mc[~np.isnan(muon_mc.fmin)] - # require (r/r_psup)^3 < 0.9 - ev = ev[ev.r_psup < 0.9] - ev_mc = ev_mc[ev_mc.r_psup < 0.9] - muon_mc = muon_mc[muon_mc.r_psup < 0.9] + # require (r < av radius) + ev = ev[ev.r < AV_RADIUS] + ev_mc = ev_mc[ev_mc.r < AV_RADIUS] + muon_mc = muon_mc[muon_mc.r < AV_RADIUS] # require psi < 6 ev = ev[ev.psi < 6] -- cgit