blob: c0c1791cf6cf3d0fadab194203d2798f79ecfd18 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
|
/* Copyright (c) 2019, Anthony Latorre <tlatorre at uchicago>
*
* This program is free software: you can redistribute it and/or modify it
* under the terms of the GNU General Public License as published by the Free
* Software Foundation, either version 3 of the License, or (at your option)
* any later version.
* This program is distributed in the hope that it will be useful, but WITHOUT
* ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
* FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for
* more details.
* You should have received a copy of the GNU General Public License along with
* this program. If not, see <https://www.gnu.org/licenses/>.
*/
#ifndef DC_H
#define DC_H
#include "event.h"
#include "zebra.h"
#include <stdint.h>
/* Data cleaning bitmasks. */
#define DC_MUON 0x1
#define DC_JUNK 0x2
#define DC_CRATE_ISOTROPY 0x4
#define DC_QVNHIT 0x8
#define DC_NECK 0x10
#define DC_FLASHER 0x20
#define DC_ESUM 0x40
#define DC_OWL 0x80
#define DC_OWL_TRIGGER 0x100
#define DC_FTS 0x200
#define DC_ITC 0x400
#define DC_BREAKDOWN 0x800
/* Minimum number of hits required for a crate to be considered the source of a
* breakdown. */
#define MIN_NHIT_BREAKDOWN 256
/* Minimum number of hits for a crate to calculate the median TAC in the
* breakdown cut. */
#define MIN_NHIT_CRATE 20
/* Length of the sliding window used in the ITC cut (ns). */
#define ITC_TIME_WINDOW 93.0
/* Fraction of hits which must be in a sliding window of 93 ns to pass the ITC
* cut. */
#define ITC_TIME_FRACTION 0.6
/* Height of PMTs which are considered to be at the "bottom" of the detector
* for the neck cut. Comes from the SNOOP definition of a neck event in SNO. */
#define NECK_BOTTOM_HEIGHT -425.0
/* Fraction of normal PMT hits which must be below NECK_BOTTOM_HEIGHT. */
#define NECK_BOTTOM_FRACTION 0.5
/* Delta-t threshold for PMT pairs in the FTS cut (ns). */
#define FTS_DT_THRESH 25.0
/* Distance threshold for PMT pairs in the FTS cut (cm). */
#define FTS_DIST_THRESH 300.0
/* Minimum number of PMT pairs required for the FTS cut. */
#define FTS_COUNT_THRESH 15
/* For the FTS cut the median time difference between PMT pairs within 3 meters
* is computed and if this difference is greater than FTS_MEDIAN_CUT the event
* doesn't pass. */
#define FTS_MEDIAN_CUT 6.8
/* Minimum number of normal PMTs which must be hit to be tagged as an incoming
* muon. */
#define MUON_MIN_NHIT 150
/* Minimum number of OWL PMTs which must be hit to be tagged as an incoming
* muon. */
#define MUON_MIN_OWL_NHIT 5
/* Minimum value for the maximum OWL QHS value. */
#define MUON_OWL_QHS 20.0
/* Percentile of PSUP PMT times to define early OWL hits. */
#define MUON_OWL_TIME_PERCENTILE 0.1
/* QvNHIT ratio threshold. */
#define QRATIO_THRESHOLD 0.25
int is_itc(event *ev);
int is_fts(event *ev);
int is_owl_trigger(event *ev);
int is_owl(event *ev);
int is_esum(event *ev);
uint32_t get_dc_word(event *ev, zebraFile *f, zebraBank *bmast, zebraBank *bev);
int is_muon(event *ev);
int junk_cut(zebraFile *f, zebraBank *bmast, zebraBank *bev);
int crate_isotropy(event *ev);
int qvnhit(event *ev);
double get_neck_tube_cut_time_diff(event *ev);
int is_neck_event(event *ev);
int is_breakdown(event *ev);
int is_flasher(event *ev);
#endif
|
