blob: 4635806e91c5cd16891a809eae7c695817b7d721 (
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
|
/* 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
/* 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
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);
int is_neck_event(event *ev);
int is_flasher(event *ev);
#endif
|
