1 #ifndef JSHOWERPOINTSIMPLEX_INCLUDE
2 #define JSHOWERPOINTSIMPLEX_INCLUDE
53 namespace JRECONSTRUCTION {}
54 namespace JPP {
using namespace JRECONSTRUCTION; }
56 namespace JRECONSTRUCTION {
74 using JRegressor_t::operator();
124 for (JEvt::const_iterator shower = in.begin(); shower != in.end(); ++shower) {
130 data.reserve(dataL0.size() + dataL1.size());
150 for (buffer_type::const_iterator i = dataL1.begin(); i != dataL1.end(); ++i) {
157 buffer_type::iterator __end = data.end();
159 for (buffer_type::const_iterator i = dataL0.begin(); i != dataL0.end(); ++i) {
161 if (find_if(data.begin(), __end, bind2nd(equal_to<JDAQModuleIdentifier>(), i->getModuleID())) == __end) {
169 const int NDF =
getCount(data.begin(), data.end()) - this->step.size();
178 JPoint4D vxs(
JVector3D(vx.getX() + x, vx.getY() + y, vx.getZ() + z), vx.getT() + t);
180 this->step.resize(4);
186 double chi2 = (*this)(vxs, data.begin(), data.end());
189 this->value.getT(), shower->getE());
static int debug
debug level (default is off).
Template definition of a data regressor of given model.
class to handle the second position fit of the shower reconstruction, mainly dedicated for ORCA ...
const JModuleRouter & router
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.
Data structure for vertex fit.
JTOOLS::JRange< double > JTimeRange
Type definition for time range (unit [ns]).
Router for direct addressing of module data in detector data structure.
static const int JSHOWERPOINTSIMPLEX
*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
Basic data structure for time and time over threshold information of hit.
Data structure for detector geometry and calibration.
JTOOLS::JRange< double > time_grid_ns
JRegressor< JPoint4D, JSimplex > JRegressor_t
Auxiliary class to extract a subset of optical modules from a detector.
JEvt operator()(const KM3NETDAQ::JDAQEvent &event, const JFIT::JEvt &in)
Declaration of the member function that actually performs the reconstruction.
Regressor function object for JPoint4D fit using JSimplex minimiser.
std::vector< hit_type > buffer_type
Data structure for vector in three dimensions.
JTOOLS::JRange< double > pos_grid_m
JPosition3D getPosition(const Vec &pos)
Get position.
JTRIGGER::JHitR1 hit_type
JShowerPointSimplex(const JShowerPointSimplexParameters_t ¶meters, const JModuleRouter &router, const int debug=0)
Parameterized constructor.
Direct access to module in detector data structure.
Data structure for L2 parameters.
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.
Data structure for set of track fit results.
Auxiliary class to define a range between two values.
Simple fit method based on Powell's algorithm, see reference: Numerical Recipes in C++...
int getCount(const T &hit)
Get hit count.
Reduced data structure for L1 hit.
static int MAXIMUM_ITERATIONS
maximal number of iterations
JTOOLS::JRange< double > TWindow_ns
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
Tukey's biweight M-estimator.
Match operator for Cherenkov light from shower in any direction.