Declaration of the member function that actually performs the reconstruction.
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();
173 for(
double x = posGrid_m.getLowerLimit(); x <= posGrid_m.getUpperLimit(); x +=
pos_step_m){
174 for(
double y = posGrid_m.getLowerLimit(); y <= posGrid_m.getUpperLimit(); y +=
pos_step_m){
175 for(
double z = posGrid_m.getLowerLimit(); z <= posGrid_m.getUpperLimit(); z +=
pos_step_m){
176 for(
double t = timeGrid_ns.getLowerLimit(); t <= timeGrid_ns.getUpperLimit(); t +=
time_step_ns){
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());
const JModuleRouter & router
Data structure for direction in three dimensions.
JTOOLS::JTimeRange time_grid_ns
Data structure for vertex fit.
static const int JSHOWERPOINTSIMPLEX
double getQuality(const double chi2, const int N, const int NDF)
Get quality of fit.
std::vector< hit_type > buffer_type
Data structure for vector in three dimensions.
JTOOLS::JTimeRange TWindow_ns
JTOOLS::JRange< double > pos_grid_m
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.
Data structure for set of track fit results.
int getCount(const T &hit)
Get hit count.
JPosition3D getPosition(const JFit &fit)
Get position.
Template specialisation of class JModel to match hit with bright point.