1 #ifndef JSHOWERPOSITIONFIT_INCLUDE
2 #define JSHOWERPOSITIONFIT_INCLUDE
58 namespace JRECONSTRUCTION {}
59 namespace JPP {
using namespace JRECONSTRUCTION; }
61 namespace JRECONSTRUCTION {
81 using JRegressor_t::operator();
97 const std::string pdfFile,
118 for (
int i = 0;
i !=
En; ++
i) {
122 this->parameters.resize(4);
149 for (JEvt::const_iterator shower = in.begin(); shower != in.end(); ++shower) {
157 for (JDataL0_t::const_iterator
i = dataL0.begin();
i != dataL0.end(); ++
i) {
162 data.push_back(
JHitW0(*
i, rate_Hz));
179 double chi2 = (*this)(sh, data.begin(), __end);
182 this->value.getT(), this->value.getE());
static int debug
debug level (default is off).
Template definition of a data regressor of given model.
int NMax
maximum number of iterations
static double EPSILON
maximal distance to minimum
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.
JRegressor< JPoint4E, JGandalf > JRegressor_t
Router for direct addressing of module data in detector data structure.
double Emin_GeV
minimum energy to scan
#define THROW(JException_t, A)
Marco for throwing exception with std::ostream compatible message.
double getRate() const
Get default rate.
*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 ...
static double Vmax_npe
Maximal integral of PDF [npe].
Basic data structure for time and time over threshold information of hit.
const JModuleRouter & router
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
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.
static JTimeRange T_ns
Time window with respect to Cherenkov hypothesis [ns].
double TMax_ns
maximum time for local coincidences [ns]
double TTS_ns
transition-time spread [ns]
int En
number of points to scan in energy range
Data structure for vertex fit.
JPosition3D getPosition(const Vec &pos)
Get position.
double DMax_m
maximal distance to optical module [m]
T pow(const T &x, const double y)
Power .
Router for fast addressing of summary data in KM3NETDAQ::JDAQSummaryslice data structure as a functio...
Direct access to module in detector data structure.
Fit method based on the Levenberg-Marquardt method.
Reduced data structure for L1 hit.
Auxiliary class for a hit with background rate value.
Data regression method for JFIT::JPoint4E from a bright point isoptropic emission PDF...
static parameter_type pT()
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.
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
static parameter_type pX()
static parameter_type pY()
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
double Emax_GeV
maximum energy to scan
static int MAXIMUM_ITERATIONS
maximal number of iterations
Template specialisation of class JModel to match hit with bright point.
Match operator for Cherenkov light from shower in any direction.
Basic data structure for L1 hit.
static parameter_type pZ()
Regressor function object for JPoint4E fit using JGandalf minimiser.