1 #ifndef __JACOUSTICS__JEVT__
2 #define __JACOUSTICS__JEVT__
23 namespace JACOUSTICS {}
24 namespace JPP {
using namespace JACOUSTICS; }
26 namespace JACOUSTICS {
94 out << setw(4) << fit.
id <<
' '
193 UNIXTimeStart(range.getLowerLimit()),
194 UNIXTimeStop (range.getUpperLimit()),
220 in >>
object.UNIXTimeStart;
221 in >>
object.UNIXTimeStop;
242 out <<
object.UNIXTimeStart;
243 out <<
object.UNIXTimeStop;
278 if (first.
oid == second.
oid) {
291 return first.
oid < second.
oid;
325 return event.UNIXTimeStart;
366 out <<
event.oid << endl
367 <<
FIXED(20,5) <<
event.UNIXTimeStart << endl
368 <<
FIXED(20,5) <<
event.UNIXTimeStop << endl
369 << setw(5) <<
event.nhit <<
' '
370 << setw(5) <<
event.nfit <<
' '
371 << setw(4) <<
event.npar << endl
372 <<
FIXED(12,3) <<
event.chi2 <<
'/'
373 <<
FIXED(7,1) <<
event.ndf << endl;
375 for (JEvt::const_iterator fit = event.begin(); fit !=
event.end(); ++fit) {
391 in >>
static_cast<JHead&
> (*this);
406 out << static_cast<const JHead&> (*this);
407 out << static_cast<const std::vector<JFit>&>(*this);
friend JReader & operator>>(JReader &in, JHead &object)
Read head from input.
int npar
number of fit parameters
Interface for binary output.
ClassDefOverride(JFit, 2)
int nfit
number of hits used in fit (after outlier removal)
JEvaluator()
Default constructor.
JHead(const std::string &oid, const JTimeRange &range, const int nhit, const int nfit, const int npar, const double ndf, const double chi2)
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
friend JWriter & operator<<(JWriter &out, const JHead &object)
Write head to output.
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.
friend JReader & operator>>(JReader &in, JFit &object)
Read fit from input.
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.
Forward declaration of binary 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.
Interface for binary input.
double ty2
2nd order correction of slope dy/dz
double vs
stretching factor
#define ClassDef(name, version)
z range($ZMAX-$ZMIN)< $MINIMAL_DZ." fi fi typeset -Z 4 STRING typeset -Z 2 FLOOR JPlot1D -f $
double UNIXTimeStart
start time
then fatal The output file must have the wildcard in the e g root fi eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY JAcoustics sh $DETECTOR_ID source JAcousticsToolkit sh CHECK_EXIT_CODE typeset A EMITTERS get_tripods $WORKDIR tripod txt EMITTERS get_transmitters $WORKDIR transmitter txt EMITTERS for EMITTER in
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
friend JWriter & operator<<(JWriter &out, const JFit &object)
Write fit to output.
virtual ~JHead()
Virtual destructor.
virtual JWriter & write(JWriter &out) const override
Write to output.
virtual JReader & read(JReader &in) override
Read from input.
Auxiliary data structure for floating point format specification.
JEvt()
Default constructor.
STD extensions for binary I/O.
JFit()
Default constructor.