Jpp
|
Auxiliary classes and methods for linear and iterative data regression. More...
Classes | |
struct | JAbstractMinimiser |
Abstract minimiser. More... | |
struct | JAbstractRegressor |
Abstract class for global fit method. More... | |
class | JAtmosphericMuon |
Auxiliary class to evaluate atmospheric muon hypothesis. More... | |
class | JEnergy |
Data structure for fit of energy. More... | |
class | JEnergyCorrection |
Auxiliary class for correction of energy determined by JEnergy.cc. More... | |
struct | JEstimator |
Template definition of linear fit. More... | |
class | JEstimator< JLine1Z > |
Linear fit of straight line parallel to z-axis to set of hits (objects with position and time). More... | |
class | JEstimator< JPoint3D > |
Linear fit of crossing point (position) between axes (objects with position and direction). More... | |
class | JEstimator< JPoint4D > |
Linear fit of bright point (position and time) between hits (objects with position and time). More... | |
struct | JEvent |
Auxiliary class for historical event. More... | |
class | JEvt |
Data structure for set of track fit results. More... | |
class | JFit |
Data structure for track fit results. More... | |
class | JGandalf |
Fit method based on the Levenberg-Marquardt method. More... | |
struct | JHistory |
Container for historical events. More... | |
struct | JHitL1Comparator |
Auxiliary class for permutations of L1 hits. More... | |
class | JHitW0 |
Auxiliary class for a hit with background rate value. More... | |
struct | JK40 |
Auxiliary class for converting various rates to expectation values of the number of hits within a given time interval. More... | |
struct | JK40Hit |
Auxiliary class for simultaneously handling light yields and response of module. More... | |
class | JLine1Z |
Data structure for fit of straight line paralel to z-axis. More... | |
class | JLine3EZ |
Data structure for fit of straight line in positive z-direction with energy. More... | |
class | JLine3Z |
Data structure for fit of straight line in positive z-direction. More... | |
class | JMatrixNZ |
Determination of the co-variance matrix of hits for a track along z-axis (JFIT::JLine1Z). More... | |
struct | JMEstimator |
Interface for maximum likelihood estimator (M-estimator). More... | |
struct | JMEstimatorLinear |
Linear M-estimator. More... | |
struct | JMEstimatorLorentzian |
Lorentzian M-estimator. More... | |
struct | JMEstimatorNormal |
Normal M-estimator. More... | |
struct | JMEstimatorNull |
Null M-estimator. More... | |
struct | JMEstimatorTukey |
Tukey's biweight M-estimator. More... | |
struct | JModel |
Auxiliary class to match data points with given model. More... | |
struct | JModel< JEnergy > |
Template specialisation of class JModel to match hit with muon energy. More... | |
struct | JModel< JLine1Z > |
Template specialisation of class JModel to match hit with muon trajectory along z-axis. More... | |
struct | JModel< JPoint4D > |
Template specialisation of class JModel to match hit with bright point. More... | |
class | JMuonGandalf |
class to handle Muon Gandalf reconstruction first, two other angular reconstructionS should be ran : JMuonPrefit and JMuonSimplex More... | |
struct | JMuonGandalfParameters_t |
struct | JMuonParameters |
class | JMuonPrefit |
class to handle Prefit angular reconstruction. More... | |
struct | JMuonPrefitParameters_t |
Data struct that holds the parameters for JMuonPrefit angular reconstruction This is part of the constructor parameters of JMuonPrefit More... | |
class | JMuonSimplex |
class to handle the Simplex angular reconstruction. More... | |
struct | JMuonSimplexParameters_t |
struct | JNPE |
Auxiliary class for handling various light yields. More... | |
struct | JNPEHit |
Auxiliary class for simultaneously handling light yields and response of PMT. More... | |
class | JORCAShowerFit |
class | JORCAShowerPositionFit |
class to handle the second step of the shower reconstruction, mainly dedicated for ORCA More... | |
class | JORCAShowerPrefit |
class to handle first step of the shower reconstruction, mainly dedicated for ORCA More... | |
struct | JPMTW0 |
Auxiliary class for handling PMT geometry, rate and response. More... | |
class | JPoint3D |
Data structure for position fit. More... | |
class | JPoint4D |
Data structure for vertex fit. More... | |
class | JPointing |
Auxiliary class to compare fit results with respect to a reference direction (e.g. More... | |
struct | JQualitySorter |
General purpose sorter of fit results. More... | |
struct | JRegressor |
Template definition of a data regressor of given model. More... | |
struct | JRegressor< JEnergy > |
Regressor function object for JEnergy fit. More... | |
struct | JRegressor< JLine3Z, JGandalf > |
Regressor function object for JLine3Z fit using JGandalf minimiser. More... | |
struct | JRegressor< JLine3Z, JSimplex > |
Regressor function object for JLine3Z fit using JSimplex minimiser. More... | |
struct | JRegressor< JPoint4D, JSimplex > |
Regressor function object for JPoint4D fit using JSimplex minimiser. More... | |
struct | JRegressor< JShower3EZ, JSimplex > |
Regressor function object for JShower3EZ fit using JSimplex minimiser. More... | |
class | JShower3EZ |
Data structure for fit of straight line in positive z-direction with energy. More... | |
class | JShower3Z |
Data structure for cascade in positive z-direction. More... | |
struct | JShowerFitParameters_t |
struct that holds the Parameters used for JORCAShowerFit More... | |
struct | JShowerParameters |
struct | JShowerPositionFitParameters_t |
struct | JShowerPrefitParameters_t |
class | JSimplex |
Simple fit method based on Powell's algorithm, see reference: Numerical Recipes in C++, W.H. More... | |
struct | JStart |
Auxiliary class for start or end point evaluation. More... | |
class | JVectorNZ |
Determination of the time residual vector of hits for a track along z-axis (JFIT::JLine1Z). More... | |
Typedefs | |
typedef JTOOLS::JRange< double > | JZRange |
Enumerations | |
enum | JFitApplication_t { JMUONBEGIN = 0, JMUONPREFIT, JMUONSIMPLEX, JMUONGANDALF, JMUONENERGY, JMUONSTART, JLINEFIT, JMUONEND = 99, JSHOWERBEGIN = 100, JSHOWERPREFIT, JSHOWERPOSITIONFIT, JSHOWERCOMPLETEFIT, JSHOWER_BJORKEN_Y, JSHOWEREND = 199, JUSERBEGIN = 1000, JMUONVETO, JPRESIM, JMUONPATH, JMCEVT, JUSEREND = 1099, JKM3DELTAPOS = 10000 } |
enum | JFitParameter_t { JGANDALF_BETA0_RAD = 0, JGANDALF_BETA1_RAD = 1, JGANDALF_CHI2 = 2, JGANDALF_NUMBER_OF_HITS = 3, JENERGY_ENERGY = 4, JENERGY_CHI2 = 5, JGANDALF_LAMBDA = 6, JGANDALF_NUMBER_OF_ITERATIONS = 7, JSTART_NPE_MIP = 8, JSTART_NPE_MIP_TOTAL = 9, JSTART_LENGTH_METRES = 10, JVETO_NPE = 11, JVETO_NUMBER_OF_HITS = 12, JENERGY_MUON_RANGE_METRES = 13, JENERGY_NOISE_LIKELIHOOD = 14, JENERGY_NDF = 15, JENERGY_NUMBER_OF_HITS = 16, JCOPY_Z_M = 17 } |
enum | JFitStatus_t { ERROR = -100, INSUFFICIENT_NUMBER_OF_HITS, INSUFFICIENT_NUMBER_OF_DOMS, OKAY = 1 } |
Functions | |
JPosition3D | getPosition (const JFit &fit) |
Get position. More... | |
JDirection3D | getDirection (const JFit &fit) |
Get direction. More... | |
JAxis3D | getAxis (const JFit &fit) |
Get axis. More... | |
JTrack3E | getTrack (const JFit &fit) |
Get track. More... | |
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. More... | |
JFit | getFit (const JHistory &history, const JLine1Z &track, const JAngle3D &angle, const double Q, const int NDF, const double energy=0.0, const int status=0) |
Get fit. More... | |
double | getDot (const JFit &first, const JFit &second) |
Get dot product. More... | |
double | getDot (const JFit &fit, const JDirection3D &dir) |
Get dot product. More... | |
double | getQuality (const double chi2, const int NDF) |
Get quality of fit. More... | |
double | getQuality (const double chi2) |
Get quality of fit. More... | |
bool | qualitySorter (const JFit &first, const JFit &second) |
Comparison of fit results. More... | |
bool | has_history (const JFit &fit, const int type) |
Test whether given fit has specified history. More... | |
bool | has_history (const JFit &fit, const JRange< int > &range) |
Test whether given fit has specified history. More... | |
bool | has_muon_prefit (const JFit &fit) |
Test whether given fit has muon prefit in history. More... | |
bool | has_muon_simplex (const JFit &fit) |
Test whether given fit has muon simplex in history. More... | |
bool | has_muon_gandalf (const JFit &fit) |
Test whether given fit has muon gandalf in history. More... | |
bool | has_muon_energy (const JFit &fit) |
Test whether given fit has muon energy in history. More... | |
bool | has_muon_start (const JFit &fit) |
Test whether given fit has muon start in history. More... | |
bool | has_muon_fit (const JFit &fit) |
Test whether given fit has muon fit in history. More... | |
bool | has_shower_prefit (const JFit &fit) |
Test whether given fit has shower prefit in history. More... | |
bool | has_shower_positionfit (const JFit &fit) |
Test whether given fit has shower position fit in history. More... | |
bool | has_shower_completefit (const JFit &fit) |
Test whether given fit has shower complete fit in history. More... | |
bool | has_shower_fit (const JFit &fit) |
Test whether given fit has shower fit in history. More... | |
template<class JTrackSelector_t > | |
bool | has_reconstructed_track (const JEvt &evt, JTrackSelector_t selector) |
Test whether given event has a track according selection. More... | |
bool | has_reconstructed_muon (const JEvt &evt) |
Test whether given event has a track with muon reconstruction. More... | |
bool | has_reconstructed_shower (const JEvt &evt) |
Test whether given event has a track with shower reconstruction. More... | |
template<class JTrackSelector_t , class JQualitySorter_t > | |
const JFit & | get_best_reconstructed_track (const JEvt &evt, JTrackSelector_t selector, JQualitySorter_t comparator) |
Get best reconstructed track. More... | |
const JFit & | get_best_reconstructed_muon (const JEvt &evt) |
Get best reconstructed muon. More... | |
const JFit & | get_best_reconstructed_shower (const JEvt &evt) |
Get best reconstructed shower. More... | |
int | getCount (const JHitL0 &hit) |
Get hit count. More... | |
int | getCount (const JHitL1 &hit) |
Get hit count. More... | |
int | getCount (const JHitR1 &hit) |
Get hit count. More... | |
template<class T > | |
int | getCount (T __begin, T __end) |
Get hit count. More... | |
double | getP (const double expval, bool hit) |
Get Poisson probability to observe a hit or not for given expectation value for the number of hits. More... | |
double | getP (const double expval, const unsigned int numberOfHits, const double precision=1.0e-5) |
Get Poisson probability to observe given number of hits or more for the given expectation value for the number of hits. More... | |
double | getChi2 (const double P) |
Get chi2 corresponding to given probability. More... | |
double | getChi2 (const double expval, bool hit) |
Get chi2 to observe a hit or not for given expectation value for the number of hits. More... | |
template<class JModel_t , class JHit_t > | |
double | getChi2 (const JModel_t &model, const JHit_t &hit, const double sigma) |
Determine chi2 of a hit for a given model and time resolution. More... | |
template<class JModel_t , class T > | |
double | getChi2 (const JModel_t &model, T __begin, T __end, const double sigma) |
Determine chi2 of data for given model and time resolution. More... | |
template<class T > | |
double | getChi2 (T __begin, T __end, const JMatrixNZ &V) |
Determine chi2 of data using full covariance matrix. More... | |
double | getChi2 (const JVectorNZ &Y, const JMatrixNZ &V) |
Determine chi2 using full covariance matrix. More... | |
template<class T > | |
double | getChi2 (const JLine1Z &track, T __begin, T __end, const JMatrixNZ &V) |
Determine chi2 of data for given track using full covariance matrix. More... | |
template<class T > | |
double | getChi2 (const JLine1Z &track, T __begin, T __end, const double alpha, const double sigma) |
Determine chi2 of data for given track and angular and time resolution. More... | |
template<class T > | |
double | getChi2 (T __begin, T __end, const JMatrixNZ &V, const int i) |
Determine difference between chi2 with and without hit using full covariance matrix. More... | |
double | getChi2 (const JVectorNZ &Y, const JMatrixNZ &V, const int i) |
Determine difference between chi2 with and without hit using full covariance matrix. More... | |
template<class JModel_t , class JFit_t , class T > | |
double | getChi2 (const JModel_t &model, const JFit_t &fit, T __begin, T __end) |
Get chi2 of data for given model and fit function. More... | |
Variables | |
static const int | JPP_RECONSTRUCTION_TYPE = 4000 |
Jpp reconstruction type for AAnet. More... | |
Auxiliary classes and methods for linear and iterative data regression.
Class used for EM shower reconstruction in ORCA. It sets the optimal values for the shower position reconstruction. These values can be changed from command line if needed.
typedef JTOOLS::JRange<double> JFIT::JZRange |
Definition at line 20 of file JFitApplications.hh.
Enumerator | |
---|---|
JGANDALF_BETA0_RAD | angular resolution [rad] from JGandalf.cc |
JGANDALF_BETA1_RAD | angular resolution [rad] from JGandalf.cc |
JGANDALF_CHI2 | chi2 from JGandalf.cc |
JGANDALF_NUMBER_OF_HITS | number of hits from JGandalf.cc |
JENERGY_ENERGY | uncorrected energy [GeV] from JEnergy.cc |
JENERGY_CHI2 | chi2 from JEnergy.cc |
JGANDALF_LAMBDA | control parameter from JGandalf.cc |
JGANDALF_NUMBER_OF_ITERATIONS | number of iterations from JGandalf.cc |
JSTART_NPE_MIP | number of photo-electrons up to the barycentre from JStart.cc |
JSTART_NPE_MIP_TOTAL | number of photo-electrons along the whole track from JStart.cc |
JSTART_LENGTH_METRES | distance between first and last hits in metres from JStart.cc |
JVETO_NPE | number of photo-electrons from JVeto.cc |
JVETO_NUMBER_OF_HITS | number of hits from JVeto.cc |
JENERGY_MUON_RANGE_METRES | range of a muon with the reconstructed energy [m] from JEnergy.cc |
JENERGY_NOISE_LIKELIHOOD | log likelihood of every hit being K40 from JEnergy.cc |
JENERGY_NDF | number of degrees of freedom from JEnergy.cc |
JENERGY_NUMBER_OF_HITS | number of hits from JEnergy.cc |
JCOPY_Z_M | true vertex position along track [m] from JCopy.cc |
Definition at line 15 of file JFitParameters.hh.
enum JFIT::JFitStatus_t |
Enumerator | |
---|---|
ERROR | |
INSUFFICIENT_NUMBER_OF_HITS | |
INSUFFICIENT_NUMBER_OF_DOMS | |
OKAY |
Definition at line 16 of file JFitStatus.hh.
|
inline |
|
inline |
|
inline |
Get fit.
history | history |
track | track |
Q | quality |
NDF | number of degrees of freedom |
energy | Energy, which would be set as 0, if the fit does not reconstruct energy |
status | status of the fit as defined in enum JFitStatus.hh |
Definition at line 116 of file JEvtToolkit.hh.
|
inline |
Get fit.
history | history |
track | track |
angle | angle |
Q | quality |
NDF | number of degrees of freedom |
energy | Energy, which would be set as 0, if the fit does not reconstruct energy |
status | status of the fit as defined in JFitStatus.hh |
Definition at line 144 of file JEvtToolkit.hh.
Get dot product.
first | first fit |
second | second fit |
Definition at line 168 of file JEvtToolkit.hh.
|
inline |
Get dot product.
fit | fit |
dir | direction |
Definition at line 181 of file JEvtToolkit.hh.
|
inline |
Get quality of fit.
The larger the quality, the better the fit.
chi2 | chi2 |
NDF | number of degrees of freedom |
Definition at line 195 of file JEvtToolkit.hh.
|
inline |
Get quality of fit.
The larger the quality, the better the fit.
chi2 | chi2 |
Definition at line 210 of file JEvtToolkit.hh.
Comparison of fit results.
first | first fit |
second | second fit |
Definition at line 223 of file JEvtToolkit.hh.
|
inline |
Test whether given fit has specified history.
fit | fit |
type | application type |
Definition at line 262 of file JEvtToolkit.hh.
|
inline |
Test whether given fit has specified history.
fit | fit |
range | application type range |
Definition at line 277 of file JEvtToolkit.hh.
|
inline |
Test whether given fit has muon prefit in history.
fit | fit |
Definition at line 291 of file JEvtToolkit.hh.
|
inline |
Test whether given fit has muon simplex in history.
fit | fit |
Definition at line 303 of file JEvtToolkit.hh.
|
inline |
Test whether given fit has muon gandalf in history.
fit | fit |
Definition at line 315 of file JEvtToolkit.hh.
|
inline |
Test whether given fit has muon energy in history.
fit | fit |
Definition at line 327 of file JEvtToolkit.hh.
|
inline |
Test whether given fit has muon start in history.
fit | fit |
Definition at line 339 of file JEvtToolkit.hh.
|
inline |
Test whether given fit has muon fit in history.
fit | fit |
Definition at line 351 of file JEvtToolkit.hh.
|
inline |
Test whether given fit has shower prefit in history.
fit | fit |
Definition at line 363 of file JEvtToolkit.hh.
|
inline |
Test whether given fit has shower position fit in history.
fit | fit |
Definition at line 375 of file JEvtToolkit.hh.
|
inline |
Test whether given fit has shower complete fit in history.
fit | fit |
Definition at line 387 of file JEvtToolkit.hh.
|
inline |
Test whether given fit has shower fit in history.
fit | fit |
Definition at line 399 of file JEvtToolkit.hh.
|
inline |
Test whether given event has a track according selection.
The track selector corresponds to the function operator bool selector(const Trk&);
.
evt | event |
selector | track selector |
Definition at line 414 of file JEvtToolkit.hh.
|
inline |
Test whether given event has a track with muon reconstruction.
evt | event |
Definition at line 426 of file JEvtToolkit.hh.
|
inline |
Test whether given event has a track with shower reconstruction.
evt | event |
Definition at line 438 of file JEvtToolkit.hh.
|
inline |
Get best reconstructed track.
The track selector corresponds to the function operator bool selector(const Trk&);
and the track comparator to bool comprator(const Trk&, const Trk&);
.
evt | event |
selector | track selector |
comparator | track comparator |
Definition at line 455 of file JEvtToolkit.hh.
Get best reconstructed muon.
evt | event |
Definition at line 480 of file JEvtToolkit.hh.
Get best reconstructed shower.
evt | event |
Definition at line 492 of file JEvtToolkit.hh.
|
inline |
|
inline |
|
inline |
|
inline |
Get hit count.
__begin | begin of L0/L1 data |
__end | end of L0/L1 data |
Definition at line 766 of file JEvtToolkit.hh.
|
inline |
Get Poisson probability to observe a hit or not for given expectation value for the number of hits.
expval | expectation value |
hit | hit |
Definition at line 37 of file JFitToolkit.hh.
|
inline |
Get Poisson probability to observe given number of hits or more for the given expectation value for the number of hits.
expval | expectation value |
numberOfHits | number of hits |
precision | precision |
Definition at line 55 of file JFitToolkit.hh.
|
inline |
Get chi2 corresponding to given probability.
P | probability |
Definition at line 80 of file JFitToolkit.hh.
|
inline |
Get chi2 to observe a hit or not for given expectation value for the number of hits.
expval | expectation value |
hit | hit |
Definition at line 93 of file JFitToolkit.hh.
|
inline |
Determine chi2 of a hit for a given model and time resolution.
The template argument JModel_t
refers to a data structure which should provide for the following method:
double getT(const JHit_t& hit) const; // get expected time of hit
model | model |
hit | hit |
sigma | sigma |
Definition at line 117 of file JFitToolkit.hh.
|
inline |
Determine chi2 of data for given model and time resolution.
The template argument JModel_t
refers to a data structure which should provide for the following method:
double getT(const value_type& hit) const; // expected time of hit
where value_type
corresponds to the value type if the input data.
model | model |
__begin | begin of data |
__end | end of data |
sigma | time resolution [ns] |
Definition at line 142 of file JFitToolkit.hh.
|
inline |
Determine chi2 of data using full covariance matrix.
__begin | begin of residuals |
__end | end of residuals |
V | covariance matrix |
Definition at line 163 of file JFitToolkit.hh.
Determine chi2 using full covariance matrix.
Y | residuals |
V | covariance matrix |
Definition at line 176 of file JFitToolkit.hh.
|
inline |
Determine chi2 of data for given track using full covariance matrix.
track | track |
__begin | begin of data |
__end | end of data |
V | covariance matrix |
Definition at line 193 of file JFitToolkit.hh.
|
inline |
Determine chi2 of data for given track and angular and time resolution.
track | track |
__begin | begin of data |
__end | end of data |
alpha | angular resolution [deg] |
sigma | time resolution [ns] |
Definition at line 215 of file JFitToolkit.hh.
|
inline |
Determine difference between chi2 with and without hit using full covariance matrix.
__begin | begin of residuals |
__end | end of residuals |
V | covariance matrix |
i | index of excluded hit |
Definition at line 235 of file JFitToolkit.hh.
Determine difference between chi2 with and without hit using full covariance matrix.
Y | residuals |
V | covariance matrix |
i | index of excluded hit |
Definition at line 257 of file JFitToolkit.hh.
|
inline |
Get chi2 of data for given model and fit function.
The template argument JFit_t
refers to a data structure which should provide for the function object operator:
double operator()(const JModel_t& model, const value_type&) const; // chi2
where JModel_t
refers to the given model and value_type
to the value type if the input data. The return value should correspond to the chi2 of the hit.
model | model |
fit | fit function |
__begin | begin of data |
__end | end of data |
Definition at line 284 of file JFitToolkit.hh.
|
static |
Jpp reconstruction type for AAnet.
Definition at line 18 of file JFitApplications.hh.