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();
96 const std::string pdfFile,
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));
156 sort(data.begin(), data.end(),
compare);
166 if (shower->getE() > 0.1)
171 double chi2 = (*this)(vx, data.begin(), __end);
174 this->value.getT(), shower->getE());
188 static const struct {
201 using namespace KM3NETDAQ;
203 if (equal_to<JDAQPMTIdentifier>()(first, second))
204 return less<JHit>()(first, second);
206 return less<JDAQPMTIdentifier>()(first, second);
static int debug
debug level (default is off).
Template definition of a data regressor of given model.
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.
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.
then echo The file $DIR KM3NeT_00000001_00000000 root already please rename or remove it first
Regressor function object for JPoint4D fit using JGandalf minimiser.
static const int JSHOWERPOSITIONFIT
Basic data structure for L0 hit.
const JSummaryRouter & summary
Auxiliary class to extract a subset of optical modules from a detector.
static struct JRECONSTRUCTION::JShowerPositionFit::@58 compare
Auxiliary data structure for sorting of hits.
do set_variable OUTPUT_DIRECTORY $WORKDIR T
JPosition3D getPosition(const Vec &pos)
Get position.
Router for fast addressing of summary data in JDAQSummaryslice data structure as a function of the op...
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].
JFit getFit(const JHistory &history, const JTrack3D &track, const double Q, const int NDF, const double energy=0.0, const int status=0)
Get fit.
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.
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.
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 source JAcoustics sh $DETECTOR_ID typeset A TRIPODS get_tripods $WORKDIR tripod txt TRIPODS 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].