1 #ifndef __JACOUSTICS__JEVTTOOLKIT__
2 #define __JACOUSTICS__JEVTTOOLKIT__
18 namespace JACOUSTICS {}
19 namespace JPP {
using namespace JACOUSTICS; }
21 namespace JACOUSTICS {
33 return NDF - 0.25 * chi2 / NDF;
51 for (
T i = __begin;
i != __end; ++
i) {
53 buffer.insert(
i->getID());
75 for (
T i = __begin;
i != __end; ++
i) {
76 buffer[
i->getString()].insert(
i->getID());
79 if (!buffer.empty()) {
81 size_t n = numeric_limits<size_t>::max();
83 for (
map<
int,
set<int> >::const_iterator
i = buffer.begin();
i != buffer.end(); ++
i) {
84 if (
i->second.size() <
n) {
117 return JFit(
id,
string.tx,
string.ty,
string.tx2,
string.ty2,
string.
vs);
144 for (JEvt::const_iterator
i = evt.begin();
i != evt.end(); ++
i) {
167 this->push_back(
getFit(
i->first,
i->second));
size_t getMinimumNumberOfEmitters(T __begin, T __end)
Get minimum number of emitters for any string in data.
double getQuality(const double chi2, const int NDF)
Get quality of fit.
getModel(const T &object)
Constructor.
JMODEL::JString getString(const JFit &fit)
Get model parameters of string.
size_t getNumberOfEmitters(T __begin, T __end)
Get number of emitters.
Model for fit to acoustics data.
Auxiliary data structure to convert model to event.
do set_variable OUTPUT_DIRECTORY $WORKDIR T
double ty2
2nd order correction of slope dy/dz
double vs
stretching factor
JACOUSTICS::JModel::string_type string
JFit getFit(const int id, const JMODEL::JString &string)
Get fit parameters of string.
getEvt(const JHead &header, const JModel &model)
Constructor.
then if[[!-f $DETECTOR]] then JDetector sh $DETECTOR fi cat $WORKDIR trigger_parameters txt<< EOFtrigger3DMuon.enabled=1;trigger3DMuon.numberOfHits=5;trigger3DMuon.gridAngle_deg=1;ctMin=0.0;TMaxLocal_ns=15.0;EOF set_variable TRIGGEREFFICIENCY_TRIGGERED_EVENTS_ONLY INPUT_FILES=() for((i=1;$i<=$NUMBER_OF_RUNS;++i));do JSirene.sh $DETECTOR $JPP_DATA/genhen.km3net_wpd_V2_0.evt.gz $WORKDIR/sirene_ ${i}.root JTriggerEfficiency.sh $DETECTOR $DETECTOR $WORKDIR/sirene_ ${i}.root $WORKDIR/trigger_efficiency_ ${i}.root $WORKDIR/trigger_parameters.txt $JPP_DATA/PMT_parameters.txt INPUT_FILES+=($WORKDIR/trigger_efficiency_ ${i}.root) done for ANGLE_DEG in $ANGLES_DEG[*];do set_variable SIGMA_NS 3.0 set_variable OUTLIERS 3 set_variable OUTPUT_FILE $WORKDIR/matrix\[${ANGLE_DEG}\deg\].root $JPP_DIR/examples/JReconstruction-f"$INPUT_FILES[*]"-o $OUTPUT_FILE-S ${SIGMA_NS}-A ${ANGLE_DEG}-O ${OUTLIERS}-d ${DEBUG}--!fiif[[$OPTION=="plot"]];then if((0));then for H1 in h0 h1;do JPlot1D-f"$WORKDIR/matrix["${^ANGLES_DEG}" deg].root:${H1}"-y"1 2e3"-Y-L TR-T""-\^"number of events [a.u.]"-> o chi2
double tx2
2nd order correction of slope dx/dz
Auxiliary data structure to convert event to model.
Model for fit to acoutsics data.