1 #ifndef __JACOUSTICS__JEVT__
2 #define __JACOUSTICS__JEVT__
22 namespace JACOUSTICS {}
23 namespace JPP {
using namespace JACOUSTICS; }
25 namespace JACOUSTICS {
89 out << setw(4) << fit.
id <<
' '
146 UNIXTimeStart(range.getLowerLimit()),
147 UNIXTimeStop (range.getUpperLimit()),
186 if (first.
oid == second.
oid) {
199 return first.
oid < second.
oid;
232 return event.UNIXTimeStart;
273 out <<
event.oid << endl
274 <<
FIXED(20,5) <<
event.UNIXTimeStart << endl
275 <<
FIXED(20,5) <<
event.UNIXTimeStop << endl
276 << setw(5) <<
event.nhit <<
' '
277 << setw(4) <<
event.npar << endl
278 <<
FIXED(12,3) <<
event.chi2 <<
' '
279 <<
FIXED(7,1) <<
event.ndf << endl;
281 for (JEvt::const_iterator fit = event.begin(); fit !=
event.end(); ++fit) {
int npar
number of fit parameters
JEvaluator()
Default constructor.
JFit(const int id, const double tx, const double ty, const double tx2, const double ty2, const double vs)
Constructor.
std::string oid
detector identifier
bool operator<(const Head &first, const Head &second)
Less than operator.
friend std::ostream & operator<<(std::ostream &out, const JFit &fit)
Write fit to output.
Auxiliary data structure for floating point format specification.
then echo The file $DIR KM3NeT_00000001_00000000 root already please rename or remove it first
virtual ~JFit()
Virtual destructor.
Auxiliary class to determine value of acoustic events.
friend std::ostream & operator<<(std::ostream &out, const JEvt &event)
Write event to output.
double UNIXTimeStop
stop time
double operator()(const JEvt &event) const
Get value of object.
double ndf
weighed number of degrees of freedom
JEvt(const JHead &header)
Constructor.
virtual ~JEvt()
Virtual destructor.
JHead()
Default constructor.
double ty2
2nd order correction of slope dy/dz
double vs
stretching factor
z range($ZMAX-$ZMIN)< $MINIMAL_DZ." fi fi typeset -Z 4 STRING typeset -Z 2 FLOOR JPlot1D -f $
double UNIXTimeStart
start time
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
virtual ~JHead()
Virtual destructor.
Auxiliary data structure for floating point format specification.
JHead(const std::string &oid, const JTimeRange &range, const int nhit, const int npar, const double ndf, const double chi2)
Constructor.
JEvt()
Default constructor.
JFit()
Default constructor.