/* 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