diff options
| -rw-r--r-- | src/dc.c | 4 | ||||
| -rwxr-xr-x | utils/plot-energy | 4 |
2 files changed, 4 insertions, 4 deletions
@@ -106,7 +106,7 @@ int is_fts(event *ev) /* Get the absolute value of the time difference. */ dt = fabs(ev->pmt_hits[i].t - ev->pmt_hits[j].t); - if (dt <= 25.0 && distance < 300.0) { + if (dt <= FTS_DT_THRESH && distance < FTS_DIST_THRESH) { deltat_array[count++] = dt; if (count >= MAX_PAIRS) goto end; } @@ -118,7 +118,7 @@ end: if (count > FTS_COUNT_THRESH) { gsl_sort(deltat_array,1,count); dtmedian = gsl_stats_median_from_sorted_data(deltat_array,1,count); - return dtmedian < FTS_MEDIAN_CUT; + return dtmedian > FTS_MEDIAN_CUT; } return 0; diff --git a/utils/plot-energy b/utils/plot-energy index 1e8cbcb..3c645d5 100755 --- a/utils/plot-energy +++ b/utils/plot-energy @@ -205,7 +205,7 @@ if __name__ == '__main__': if prompt_plus_muons.size and follower.size: # require Michel events to pass more of the SNO data cleaning cuts - michel = follower[follower.dc & (DC_JUNK | DC_CRATE_ISOTROPY | DC_QVNHIT | DC_FLASHER | DC_NECK | DC_ESUM | DC_OWL | DC_OWL_TRIGGER | DC_FTS) == DC_FTS] + michel = follower[follower.dc & (DC_JUNK | DC_CRATE_ISOTROPY | DC_QVNHIT | DC_FLASHER | DC_NECK | DC_ESUM | DC_OWL | DC_OWL_TRIGGER | DC_FTS) == 0] michel = michel[michel.nhit >= 100] @@ -242,7 +242,7 @@ if __name__ == '__main__': if prompt.size and follower.size: # neutron followers have to obey stricter set of data cleaning cuts - neutron = follower[follower.dc & (DC_JUNK | DC_CRATE_ISOTROPY | DC_QVNHIT | DC_FLASHER | DC_NECK | DC_ESUM | DC_OWL | DC_OWL_TRIGGER | DC_FTS) == DC_FTS] + neutron = follower[follower.dc & (DC_JUNK | DC_CRATE_ISOTROPY | DC_QVNHIT | DC_FLASHER | DC_NECK | DC_ESUM | DC_OWL | DC_OWL_TRIGGER | DC_FTS) == 0] neutron = neutron[~np.isnan(neutron.ftp_x) & ~np.isnan(neutron.rsp_energy)] r = np.sqrt(neutron.ftp_x**2 + neutron.ftp_y**2 + neutron.ftp_z**2) neutron = neutron[r < AV_RADIUS] |