blob: 89a7670021342c5e79841569dbd37aa0f0f88013 (
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
|
/* 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"
/* 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_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
|
