Declaration of the member function that actually performs the reconstruction.
136 for (JEvt::const_iterator shower = in.begin(); shower != in.end(); ++shower) {
146 if(sh.getE() >
E_GeV) {
156 for (JDataL0_t::const_iterator i = dataL0.begin(); i != dataL0.end(); ++i) {
160 top.insert(i->getPMTIdentifier());
162 start_dir.add(
JDirection3D(i->getPosition() - sh.getPosition()));
171 JOmega3D scan_directions(sh.getDirection(),
173 scan_angle *
PI/180.0);
175 for (JOmega3D_t::const_iterator dir = scan_directions.begin(); dir != scan_directions.end(); ++dir) {
182 for (JDetectorSubset_t::iterator module = subdetector.begin();
183 module != subdetector.end(); ++module) {
188 for (
unsigned int i = 0; i != dom.size(); ++i) {
194 buffer.push_back(
JPMTW0(dom.getPMT(i), rate_Hz, top.count(
id)));
199 JEnergy(log10(sh.getE()))), buffer.begin(), buffer.end());
203 JShower3E sh_fit(this->value.getPosition(), this->value.getDirection(),
204 this->value.getT(), this->value.getE());
208 sh_fit.
add(sh.getPosition());
211 NDF, sh_fit.getE()));
Data structure for direction in three dimensions.
Data structure for a composite optical module.
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.
Auxiliary class for handling PMT geometry, rate and response.
Data structure for vertex fit.
const JSummaryRouter & summary
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
Direction set covering (part of) solid angle.
double MaxAngle_deg_lowE
maximum scanning angle in [deg] for low energy events
std::pair< double, double > range_type
double roadWidth_m
road width [m]
double getQuality(const double chi2, const int N, const int NDF)
Get quality of fit.
double scanAngle_deg_highE
scanning angle step in [deg] for high energy events
double scanAngle_deg_lowE
scanning angle step in [deg] for low energy events
double DMin_m
radius of the hit selection for high energy events (used to speedup the code)
Data structure for fit of straight line in positive z-direction with energy.
Data structure for vector in three dimensions.
double MaxAngle_deg_highE
maximum scanning angle in [deg] for high energy events
static const int JSHOWERDIRECTIONPREFIT
JDirection3D getDirection(const JFit &fit)
Get direction.
JAxis3D & rotate_back(const JRotation3D &R)
Rotate back axis.
static const double PI
Mathematical constants.
then usage $script[distance] fi case set_variable R
Detector subset without binary search functionality.
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.
const JClass_t & getReference() const
Get reference to object.
int getCount(const T &hit)
Get hit count.
const JModuleRouter & router
Data structure for fit of energy.
JPosition3D getPosition(const JFit &fit)
Get position.
static JTimeRange T_ns
Time window with respect to Cherenkov hypothesis [ns].
Template specialisation of class JModel to match hit with bright point.
Data structure for normalised vector in positive z-direction.
JVector3D & add(const JVector3D &vector)
Add vector.