Declaration of the member function that actually performs the reconstruction.
136 for (JEvt::const_iterator shower = in.begin(); shower != in.end(); ++shower) {
144 for (JDataL0_t::const_iterator
i = dataL0.begin();
i != dataL0.end(); ++
i) {
147 top.insert(
i->getPMTIdentifier());
158 for (JDetectorSubset_t::iterator module = subdetector.begin();
159 module != subdetector.end(); ++module) {
167 for (
size_t i = 0;
i != dom.size(); ++
i) {
171 frame[
i].is_valid() &&
177 const size_t count = top.count(
id);
179 buffer.push_back(
JPMTW0(dom.getPMT(
i), rate_Hz, count));
184 this->step.resize(2);
189 JEnergy(
log10(sh.getE()))), buffer.begin(), buffer.end());
191 double NDF =
getCount(buffer.begin(), buffer.end()) - this->step.size();
193 JShower3E sh_fit(this->value.getPosition(), this->value.getDirection(),
194 this->value.getT(),
correct(this->value.getE()));
196 sh_fit.rotate_back(
R);
198 sh_fit.add(sh.getPosition());
201 NDF, sh_fit.getE()));
double roadWidth_m
road width [m]
Data structure for direction in three dimensions.
Data structure for a composite optical module.
Template specialisation of L0 builder for JHitL0 data type.
Auxiliary class for handling PMT geometry, rate and response.
Data structure for vertex fit.
double getRate() const
Get default rate.
JFit getFit(const JHistory &history, const JTrack3D &track, const double Q, const int NDF, const double energy=0.0, const int status=SINGLE_STAGE)
Get fit.
const JModuleRouter & router
const JDAQSummaryFrame & getSummaryFrame() const
Get default summary frame.
double getQuality(const double chi2, const int N, const int NDF)
Get quality of fit.
const JShowerEnergyCorrection & correct
const JSummaryRouter & summary
static const int JSHOWERFIT_ENERGY
uncorrected energy [GeV] from JShowerFit.cc
Data structure for fit of straight line in positive z-direction with energy.
Data structure for vector in three dimensions.
set_variable E_E log10(E_{fit}/E_{#mu})"
Data storage class for rate measurements of all PMTs in one module.
static const int PMT_DISABLE
KM3NeT Data Definitions v3.3.0-2-g5cc95cf https://git.km3net.de/common/km3net-dataformat.
JDirection3D getDirection(const JFit &fit)
Get direction.
double fit_step
simplex initial step
Detector subset without binary search functionality.
static const int JSHOWERCOMPLETEFIT
then JCookie sh JDataQuality D $DETECTOR_ID R
bool getPMTStatus(const JStatus &status)
Test status of PMT.
const JClass_t & getReference() const
Get reference to object.
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
int getCount(const T &hit)
Get hit count.
Data structure for fit of energy.
JPosition3D getPosition(const JFit &fit)
Get position.
bool getDAQStatus(const JDAQFrameStatus &frame, const JStatus &status)
Test status of DAQ.
Template specialisation of class JModel to match hit with bright point.
Data structure for normalised vector in positive z-direction.
static JTimeRange T_ns
Time window with respect to Cherenkov hypothesis [ns].
1.8.5