1 #ifndef __JRECONSTRUCTION__JEVT__
2 #define __JRECONSTRUCTION__JEVT__
26 namespace JRECONSTRUCTION {}
27 namespace JPP {
using namespace JRECONSTRUCTION; }
46 __x(std::numeric_limits<double>::lowest()),
47 __y(std::numeric_limits<double>::lowest()),
48 __z(std::numeric_limits<double>::lowest()),
52 __t(std::numeric_limits<double>::lowest()),
53 __Q(std::numeric_limits<double>::lowest()),
180 void move(
const double step,
const double velocity)
185 __t += step / velocity;
207 return ((
size_t) sqrt(1.0 + 8.0*
V.size()) - 1) / 2;
231 double getV(
const size_t row,
const size_t col)
const
235 const size_t __row = max(row, col) + 1;
236 const size_t __col = min(row, col);
238 return V.at((__row * (__row + 1)) / 2 - __row + __col);
257 for (
size_t row = 0; row != size; ++row) {
258 for (
size_t col = 0; col <= row; ++col) {
259 V.push_back(data[row][col]);
306 return (i >= 0 && i < (
int)
W.size());
329 double getW(
const int i,
const double value)
const
344 void setW(
const int i,
const double value)
346 if (i >= (
int)
W.size()) {
347 W.resize(i + 1, 0.0);
393 template<
class JSelector_t>
400 iterator __end = partition(this->begin(), this->end(), selector);
402 this->erase(__end, this->end());
413 template<
class JComparator_t>
415 const JComparator_t& comparator)
423 if (number_of_fits > 0 && number_of_fits < this->size()) {
425 advance(__end = this->begin(), number_of_fits);
427 partial_sort(this->begin(), __end, this->end(), comparator);
431 sort(this->begin(), __end, comparator);
434 this->erase(__end, this->end());
450 out <<
"Event: " << endl;
452 for (JEvt::const_iterator fit = event.begin(); fit !=
event.end(); ++fit) {
465 namespace JRECONSTRUCTION {
Q(UTCMax_s-UTCMin_s)-livetime_s
JHistory & add(const int type)
Add event to history.
friend std::ostream & operator<<(std::ostream &out, const JEvt &event)
Write event to output.
then usage $script< input file >[option[primary[working directory]]] nWhere option can be E
Definition for fit results.
void setW(const std::vector< double > &W)
Set associated values.
void setE(const double E)
Set energy.
void select(const size_t number_of_fits, const JComparator_t &comparator)
Select fits.
int getNDF() const
Get number of degrees of freedom.
void setW(const int i, const double value)
Set associated value.
int getN() const
Get number of associated values.
friend std::ostream & operator<<(std::ostream &out, const JHistory &history)
Write history to output stream.
Container for historical events.
This include file serves the purpose of hiding ROOT dependencies and circumphere namespace problems w...
Auxiliary class for a type holder.
void setV(const size_t size, const T &data)
Set error matrix.
JFit & add(const int type)
Add event to history.
size_t getDimensionOfErrorMatrix() const
Get dimension of error matrix.
Data structure for track fit results with history and optional associated values. ...
void setStatus(const int status)
Set status of the fit.
double getV(const size_t row, const size_t col) const
Get element of error matrix.
bool hasW(const int i) const
Check availability of value.
double getT() const
Get time.
double getX() const
Get X-position.
do set_variable OUTPUT_DIRECTORY $WORKDIR T
JFit(const JHistory &history, const double x, const double y, const double z, const int status=JRECONSTRUCTION::SINGLE_STAGE)
Constructor for storing position only.
counter_type advance(counter_type &counter, const counter_type value, const counter_type limit=std::numeric_limits< counter_type >::max())
Advance counter.
double getY() const
Get Y-position.
int getStatus() const
Get status of the fit; negative values should refer to a bad fit.
Data structure for TTree parameters.
JEvt()
Default constructor.
double getDZ() const
Get Z-slope.
JROOT::JTreeParameters getTreeParameters(JLANG::JType< JRECONSTRUCTION::JEvt >)
Get TTree parameters for given data type.
double getQ() const
Get quality.
Data structure for set of track fit results.
double getZ() const
Get Z-position.
const std::vector< double > & getV() const
Get error matrix.
double getW(const int i) const
Get associated value.
double getW(const int i, const double value) const
Get associated value.
double getDY() const
Get Y-slope.
void move(const double step, const double velocity)
Move vertex along this track with given velocity.
const std::vector< double > & getW() const
Get associated values.
const JHistory & getHistory() const
Get history.
JFit(const JHistory &history, const double x, const double y, const double z, const double dx, const double dy, const double dz, const double t, const double Q, const int NDF, const double E=0.0, const int status=JRECONSTRUCTION::SINGLE_STAGE)
Constructor.
void select(const JSelector_t &selector)
Select fits.
double getE() const
Get energy.
double getDX() const
Get X-slope.
JFit()
Default constructor.