1 #ifndef JSHOWERPOSITIONFIT_INCLUDE
2 #define JSHOWERPOSITIONFIT_INCLUDE
58 namespace JRECONSTRUCTION {}
59 namespace JPP {
using namespace JRECONSTRUCTION; }
61 namespace JRECONSTRUCTION {
80 using JRegressor_t::operator();
110 this->parameters.resize(4);
137 for (JEvt::const_iterator shower = in.begin(); shower != in.end(); ++shower) {
145 for (JDataL0_t::const_iterator
i = dataL0.begin();
i != dataL0.end(); ++
i) {
150 data.push_back(
JHitW0(*
i, rate_Hz));
166 if (shower->getE() > 0.1)
171 double chi2 = (*this)(vx, data.begin(), __end);
174 this->value.getT(), shower->getE());
static int debug
debug level (default is off).
Template definition of a data regressor of given model.
int NMax
maximum number of iterations
Data structure for direction in three dimensions.
double getQuality(const double chi2, const int NDF)
Get quality of fit.
Algorithms for hit clustering and sorting.
std::vector< T >::difference_type distance(typename std::vector< T >::const_iterator first, typename PhysicsEvent::const_iterator< T > second)
Specialisation of STL distance.
Template specialisation of L0 builder for JHitL0 data type.
Data structure for vertex fit.
Router for direct addressing of module data in detector data structure.
double roadWidth_m
road width [m]
static parameter_type pT()
double getRate() const
Get default rate.
JRegressor< JPoint4D, JGandalf > JRegressor_t
*fatal Wrong number of arguments esac JCookie sh typeset Z DETECTOR typeset Z SOURCE_RUN typeset Z TARGET_RUN set_variable PARAMETERS_FILE $WORKDIR parameters
class to handle the second position fit of the shower reconstruction, mainly dedicated for ORCA ...
double E_GeV
Energy of the shower [GeV].
Basic data structure for time and time over threshold information of hit.
const JModuleRouter & router
static parameter_type pZ()
JEvt operator()(const KM3NETDAQ::JDAQEvent &event, const JFIT::JEvt &in)
Declaration of the member function that actually performs the reconstruction.
Data structure for detector geometry and calibration.
double TMin_ns
minimum time for local coincidences [ns]
static struct JTRIGGER::JHitL0::compare compare
Regressor function object for JPoint4D fit using JGandalf minimiser.
static const int JSHOWERPOSITIONFIT
double VMax_npe
maximum number of of photo-electrons
Basic data structure for L0 hit.
const JSummaryRouter & summary
Auxiliary class to extract a subset of optical modules from a detector.
double TMax_ns
maximum time for local coincidences [ns]
double TTS_ns
transition-time spread [ns]
JPosition3D getPosition(const Vec &pos)
Get position.
Router for fast addressing of summary data in KM3NETDAQ::JDAQSummaryslice data structure as a functio...
static parameter_type pX()
static parameter_type pY()
Direct access to module in detector data structure.
Fit method based on the Levenberg-Marquardt method.
static double Vmax_npe
Maximal integral of PDF [npe].
Reduced data structure for L1 hit.
Auxiliary class for a hit with background rate value.
Data regression method for JFIT::JPoint4D from a bright point isoptropic emission PDF...
Data structure for set of track fit results.
JFit getFit(const int id, const JMODEL::JString &string)
Get fit parameters of string.
Auxiliary class to define a range between two values.
JShowerPositionFit(const JShowerPositionFitParameters_t ¶meters, const JModuleRouter &router, const JSummaryRouter &summary, const std::string pdfFile, const int debug=0)
Parameterized 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
static int MAXIMUM_ITERATIONS
maximal number of iterations
Template specialisation of class JModel to match hit with bright point.
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2] 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
Match operator for Cherenkov light from shower in any direction.
Basic data structure for L1 hit.
static JTimeRange T_ns
Time window with respect to Cherenkov hypothesis [ns].