1 #ifndef __JACOUSTICS__JEVTTOOLKIT__
2 #define __JACOUSTICS__JEVTTOOLKIT__
16 namespace JACOUSTICS {}
17 namespace JPP {
using namespace JACOUSTICS; }
19 namespace JACOUSTICS {
31 return NDF - 0.25 * chi2 / NDF;
47 for (
T i = __begin; i != __end; ++i) {
49 buffer.insert(i->getID());
78 return JFit(
id,
string.tx,
string.ty,
string.tx2,
string.ty2,
string.
vs);
92 for (JEvt::const_iterator i = evt.begin(); i != evt.end(); ++i) {
113 evt.push_back(
getFit(i->first, i->second));
double getQuality(const double chi2, const int NDF)
Get quality of fit.
JEvt getEvt(const JHead &header, const JModel &model)
Get event.
size_t getNumberOfEmitters(T __begin, T __end)
Get number of emitters.
Model for fit to acoustics data.
do set_variable OUTPUT_DIRECTORY $WORKDIR T
JMODEL::JString getModel(const JFit &fit)
Get model parameters of string.
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.
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
Model for fit to acoutsics data.