JFIT Namespace Reference

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	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...
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...
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...
static const int	JKM3DELTAPOS = 10000
	This is not a fit, this gives position information only. More...

Detailed Description

Auxiliary classes and methods for linear and iterative data regression.

Author

mdejong

mdejong, gmaggi

gmaggi

adomi, gmaggi

adomi

lquinn

adomi

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 Documentation

JZRange

typedef JTOOLS::JRange<double> JFIT::JZRange

Definition at line 21 of file JModel.hh.

Enumeration Type Documentation

JFitStatus_t

enum JFIT::JFitStatus_t

Enumerator
ERROR
INSUFFICIENT_NUMBER_OF_HITS
INSUFFICIENT_NUMBER_OF_DOMS
OKAY

Definition at line 16 of file JFitStatus.hh.

                    {
 
    ERROR             = -100,
    INSUFFICIENT_NUMBER_OF_HITS,
    INSUFFICIENT_NUMBER_OF_DOMS,
    
    OKAY              = 1
  };

Function Documentation

getPosition()

JPosition3D JFIT::getPosition ( const JFit &  fit )

inline

Get position.

Parameters

fit

fit

Returns

position

Definition at line 63 of file JEvtToolkit.hh.

  {
    return JPosition3D(fit.getX(), fit.getY(), fit.getZ());
  }

getDirection()

JDirection3D JFIT::getDirection ( const JFit &  fit )

inline

Get direction.

Parameters

fit

fit

Returns

direction

Definition at line 75 of file JEvtToolkit.hh.

  {
    return JDirection3D(fit.getDX(), fit.getDY(), fit.getDZ());
  }

getAxis()

JAxis3D JFIT::getAxis ( const JFit &  fit )

inline

Get axis.

Parameters

fit

fit

Returns

axis

Definition at line 87 of file JEvtToolkit.hh.

  {
    return JAxis3D(getPosition(fit), getDirection(fit));
  }

getTrack()

JTrack3E JFIT::getTrack ( const JFit &  fit )

inline

Get track.

Parameters

fit

fit

Returns

track

Definition at line 99 of file JEvtToolkit.hh.

  {
    return JTrack3E(JTrack3D(getAxis(fit), fit.getT()), fit.getE());
  }

getFit() [1/2]

JFit JFIT::getFit ( const JHistory &  history, const JTrack3D &  track, const double  Q, const int  NDF, const double  energy = 0.0, const int  status = 0  )

inline

Get fit.

Parameters

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

Returns

fit

Definition at line 116 of file JEvtToolkit.hh.

  {
    return JFit(history,
                track.getX(),  track.getY(),  track.getZ(),
                track.getDX(), track.getDY(), track.getDZ(),
                track.getT(),
                Q, NDF,
                energy, status);
  }

getFit() [2/2]

JFit 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  )

inline

Get fit.

Parameters

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

Returns

fit

Definition at line 144 of file JEvtToolkit.hh.

  {
    return JFit(history,
                track.getX(),  track.getY(),  track.getZ(),
                angle.getDX(), angle.getDY(), angle.getDZ(),
                track.getT(),
                Q, NDF,
                energy, status);
  }

getDot() [1/2]

double JFIT::getDot ( const JFit &  first, const JFit &  second  )

inline

Get dot product.

Parameters

first	first fit
second	second fit

Returns

dot product

Definition at line 168 of file JEvtToolkit.hh.

  {
    return JMATH::getDot(getDirection(first), getDirection(second));
  }

getDot() [2/2]

double JFIT::getDot ( const JFit &  fit, const JDirection3D &  dir  )

inline

Get dot product.

Parameters

fit	fit
dir	direction

Returns

dot product

Definition at line 181 of file JEvtToolkit.hh.

  {
    return JMATH::getDot(getDirection(fit), dir);
  }

getQuality() [1/2]

double JFIT::getQuality ( const double  chi2, const int  NDF  )

inline

Get quality of fit.

The larger the quality, the better the fit.

Parameters

chi2	chi2
NDF	number of degrees of freedom

Returns

quality

Definition at line 195 of file JEvtToolkit.hh.

  {
    const double w = NDF  -  0.25 * chi2 / NDF;
      
    return w;
  }

getQuality() [2/2]

double JFIT::getQuality ( const double  chi2 )

inline

Get quality of fit.

The larger the quality, the better the fit.

Parameters

chi2

chi2

Returns

quality

Definition at line 210 of file JEvtToolkit.hh.

  {
    return -chi2;
  }

qualitySorter()

bool JFIT::qualitySorter ( const JFit &  first, const JFit &  second  )

inline

Comparison of fit results.

Parameters

first	first fit
second	second fit

Returns

true if first fit has better quality than second; else false

Definition at line 223 of file JEvtToolkit.hh.

  {
    if (first.getHistory().size() == second.getHistory().size())
      return first.getQ() > second.getQ();
    else
      return first.getHistory().size() > second.getHistory().size();
  }

has_history() [1/2]

bool JFIT::has_history ( const JFit &  fit, const int  type  )

inline

Test whether given fit has specified history.

Parameters

fit	fit
type	application type

Returns

true if type in history; else false

Definition at line 262 of file JEvtToolkit.hh.

  {
    return std::find_if(fit.getHistory().begin(),
                        fit.getHistory().end(),
                        JLANG::make_predicate(&JEvent::type, type)) != fit.getHistory().end();
  }

has_history() [2/2]

bool JFIT::has_history ( const JFit &  fit, const JRange< int > &  range  )

inline

Test whether given fit has specified history.

Parameters

fit	fit
range	application type range

Returns

true if type in history; else false

Definition at line 277 of file JEvtToolkit.hh.

  {
    return std::find_if(fit.getHistory().begin(),
                        fit.getHistory().end(),
                        JLANG::make_find_if(&JEvent::type, range)) != fit.getHistory().end();
  }

has_muon_prefit()

bool JFIT::has_muon_prefit ( const JFit &  fit )

inline

Test whether given fit has muon prefit in history.

Parameters

fit

fit

Returns

true if muon prefit in history; else false

Definition at line 291 of file JEvtToolkit.hh.

  {
    return has_history(fit, JMUONPREFIT);
  }

has_muon_simplex()

bool JFIT::has_muon_simplex ( const JFit &  fit )

inline

Test whether given fit has muon simplex in history.

Parameters

fit

fit

Returns

true if muon simplex in history; else false

Definition at line 303 of file JEvtToolkit.hh.

  {
    return has_history(fit, JMUONSIMPLEX);
  }

has_muon_gandalf()

bool JFIT::has_muon_gandalf ( const JFit &  fit )

inline

Test whether given fit has muon gandalf in history.

Parameters

fit

fit

Returns

true if muon gandalf in history; else false

Definition at line 315 of file JEvtToolkit.hh.

  {
    return has_history(fit, JMUONGANDALF);
  }

has_muon_energy()

bool JFIT::has_muon_energy ( const JFit &  fit )

inline

Test whether given fit has muon energy in history.

Parameters

fit

fit

Returns

true if muon energy in history; else false

Definition at line 327 of file JEvtToolkit.hh.

  {
    return has_history(fit, JMUONENERGY);
  }

has_muon_start()

bool JFIT::has_muon_start ( const JFit &  fit )

inline

Test whether given fit has muon start in history.

Parameters

fit

fit

Returns

true if muon start in history; else false

Definition at line 339 of file JEvtToolkit.hh.

  {
    return has_history(fit, JMUONSTART);
  }

has_muon_fit()

bool JFIT::has_muon_fit ( const JFit &  fit )

inline

Test whether given fit has muon fit in history.

Parameters

fit

fit

Returns

true if muon fit in history; else false

Definition at line 351 of file JEvtToolkit.hh.

  {
    return has_history(fit, JRange<int>(JMUONBEGIN, JMUONEND));
  }

has_shower_prefit()

bool JFIT::has_shower_prefit ( const JFit &  fit )

inline

Test whether given fit has shower prefit in history.

Parameters

fit

fit

Returns

true if shower prefit in history; else false

Definition at line 363 of file JEvtToolkit.hh.

  {
    return has_history(fit, JSHOWERPREFIT);
  }

has_shower_positionfit()

bool JFIT::has_shower_positionfit ( const JFit &  fit )

inline

Test whether given fit has shower position fit in history.

Parameters

fit

fit

Returns

true if shower position fit in history; else false

Definition at line 375 of file JEvtToolkit.hh.

  {
    return has_history(fit, JSHOWERPOSITIONFIT);
  }

has_shower_completefit()

bool JFIT::has_shower_completefit ( const JFit &  fit )

inline

Test whether given fit has shower complete fit in history.

Parameters

fit

fit

Returns

true if shower complete fit in history; else false

Definition at line 387 of file JEvtToolkit.hh.

  {
    return has_history(fit, JSHOWERCOMPLETEFIT);
  }

has_shower_fit()

bool JFIT::has_shower_fit ( const JFit &  fit )

inline

Test whether given fit has shower fit in history.

Parameters

fit

fit

Returns

true if shower fit in history; else false

Definition at line 399 of file JEvtToolkit.hh.

  {
    return has_history(fit, JRange<int>(JSHOWERBEGIN, JSHOWEREND));
  }

has_reconstructed_track()

template<class JTrackSelector_t >

bool JFIT::has_reconstructed_track ( const JEvt &  evt, JTrackSelector_t  selector  )

inline

Test whether given event has a track according selection.

The track selector corresponds to the function operator bool selector(const Trk&);.

Parameters

evt	event
selector	track selector

Returns

true if at least one corresponding track; else false

Definition at line 414 of file JEvtToolkit.hh.

  {
    return std::find_if(evt.begin(), evt.end(), selector) != evt.end();
  }

has_reconstructed_muon()

bool JFIT::has_reconstructed_muon ( const JEvt &  evt )

inline

Test whether given event has a track with muon reconstruction.

Parameters

evt

event

Returns

true if at least one reconstructed muon; else false

Definition at line 426 of file JEvtToolkit.hh.

  {
    return has_reconstructed_track(evt, has_muon_fit);
  }

has_reconstructed_shower()

bool JFIT::has_reconstructed_shower ( const JEvt &  evt )

inline

Test whether given event has a track with shower reconstruction.

Parameters

evt

event

Returns

true if at least one reconstructed shower; else false

Definition at line 438 of file JEvtToolkit.hh.

  {
    return has_reconstructed_track(evt, has_shower_fit);
  }

get_best_reconstructed_track()

template<class JTrackSelector_t , class JQualitySorter_t >

const JFit& JFIT::get_best_reconstructed_track ( const JEvt &  evt, JTrackSelector_t  selector, JQualitySorter_t  comparator  )

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&);.

Parameters

evt	event
selector	track selector
comparator	track comparator

Returns

track

Definition at line 455 of file JEvtToolkit.hh.

  {
    JEvt::const_iterator p = std::find_if(evt.begin(), evt.begin(), selector);
 
    for (JEvt::const_iterator i = p; i != evt.end(); ++i) {
      if (selector(*i) && comparator(*i, *p)) {
        p = i;
      }
    }
 
    if (p != evt.end())
      return *p;
    else
      THROW(JIndexOutOfRange, "This event has no fit with given selector.");
  }

get_best_reconstructed_muon()

const JFit& JFIT::get_best_reconstructed_muon ( const JEvt &  evt )

inline

Get best reconstructed muon.

Parameters

evt

event

Returns

track

Definition at line 480 of file JEvtToolkit.hh.

  {
    return get_best_reconstructed_track(evt, has_muon_fit, qualitySorter);
  }

get_best_reconstructed_shower()

const JFit& JFIT::get_best_reconstructed_shower ( const JEvt &  evt )

inline

Get best reconstructed shower.

Parameters

evt

event

Returns

track

Definition at line 492 of file JEvtToolkit.hh.

  {
    return get_best_reconstructed_track(evt, has_shower_fit, qualitySorter);
  }

getCount() [1/4]

int JFIT::getCount ( const JHitL0 &  hit )

inline

Get hit count.

Parameters

hit

L0 hit

Returns

1

Definition at line 728 of file JEvtToolkit.hh.

  {
    return 1;
  }

getCount() [2/4]

int JFIT::getCount ( const JHitL1 &  hit )

inline

Get hit count.

Parameters

hit

hit

Returns

count

Definition at line 740 of file JEvtToolkit.hh.

  {
    return hit.size();
  }

getCount() [3/4]

int JFIT::getCount ( const JHitR1 &  hit )

inline

Get hit count.

Parameters

hit

hit

Returns

count

Definition at line 752 of file JEvtToolkit.hh.

  {
    return hit.getN();
  }

getCount() [4/4]

template<class T >

int JFIT::getCount ( T  __begin, T  __end  )

inline

Get hit count.

Parameters

__begin	begin of L0/L1 data
__end	end of L0/L1 data

Returns

hit count

Definition at line 766 of file JEvtToolkit.hh.

  {
    int count = 0;
    
    for (T hit = __begin; hit !=__end; ++hit) {
      count += getCount(*hit);
    }
    
    return count;
  }

getP() [1/2]

double JFIT::getP ( const double  expval, bool  hit  )

inline

Get Poisson probability to observe a hit or not for given expectation value for the number of hits.

Parameters

expval	expectation value
hit	hit

Returns

probability

Definition at line 37 of file JFitToolkit.hh.

  {
    if (hit)
      return 1.0 - exp(-expval);    // 1 - P(0)
    else
      return exp(-expval);          // P(0)
  }

getP() [2/2]

double JFIT::getP ( const double  expval, const unsigned int  numberOfHits, const double  precision = 1.0e-5  )

inline

Get Poisson probability to observe given number of hits or more for the given expectation value for the number of hits.

Parameters

expval	expectation value
numberOfHits	number of hits
precision	precision

Returns

propability

Definition at line 55 of file JFitToolkit.hh.

  {
    double P = 0.0;
 
    for (int i = numberOfHits; ; ++i) {
      
      const double p = TMath::PoissonI(i, expval);
      
      P += p;
      
      if ((double) i > expval && p < precision) {
        break;
      }
    }
 
    return P;
  }

getChi2() [1/9]

double JFIT::getChi2 ( const double  P )

inline

Get chi2 corresponding to given probability.

Parameters

P	probability

Returns

chi2

Definition at line 80 of file JFitToolkit.hh.

  {
    return -log(P);
  }

getChi2() [2/9]

double JFIT::getChi2 ( const double  expval, bool  hit  )

inline

Get chi2 to observe a hit or not for given expectation value for the number of hits.

Parameters

expval	expectation value
hit	hit

Returns

probability

Definition at line 93 of file JFitToolkit.hh.

  {
    if (hit)
      return -log(1.0 - exp(-expval));
    else
      return expval;
  }

getChi2() [3/9]

template<class JModel_t , class JHit_t >

double JFIT::getChi2 ( const JModel_t &  model, const JHit_t &  hit, const double  sigma  )

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

Parameters

model	model
hit	hit
sigma	sigma

Returns

chi2

Definition at line 117 of file JFitToolkit.hh.

  {
    const double ds = (hit.getT() - model.getT(hit)) / sigma;
 
    return ds * ds;
  }

getChi2() [4/9]

template<class JModel_t , class T >

double JFIT::getChi2 ( const JModel_t &  model, T  __begin, T  __end, const double  sigma  )

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.

Parameters

model	model
__begin	begin of data
__end	end of data
sigma	time resolution [ns]

Returns

chi2

Definition at line 142 of file JFitToolkit.hh.

  {
    double chi2 = 0.0;
 
    for (T hit = __begin; hit != __end; ++hit) {
      chi2 += getChi2(model, *hit, sigma);
    }
 
    return chi2;
  }

getChi2() [5/9]

double JFIT::getChi2 ( const JVectorNZ &  Y, const JMatrixNZ &  V  )

inline

Determine chi2 using full covariance matrix.

Parameters

Y	residuals
V	covariance matrix

Returns

chi2

Definition at line 161 of file JFitToolkit.hh.

  {
    return V.getDot(Y);
  }

getChi2() [6/9]

template<class T >

double JFIT::getChi2 ( const JLine1Z &  track, T  __begin, T  __end, const JMatrixNZ &  V  )

inline

Determine chi2 of data for given track using full covariance matrix.

Parameters

track	track
__begin	begin of data
__end	end of data
V	covariance matrix

Returns

chi2

Definition at line 178 of file JFitToolkit.hh.

  {
    const JVectorNZ Y(track, __begin, __end);
 
    return getChi2(Y, V);
  }

getChi2() [7/9]

template<class T >

double JFIT::getChi2 ( const JLine1Z &  track, T  __begin, T  __end, const double  alpha, const double  sigma  )

inline

Determine chi2 of data for given track and angular and time resolution.

Parameters

track	track
__begin	begin of data
__end	end of data
alpha	angular resolution [deg]
sigma	time resolution [ns]

Returns

chi2

Definition at line 200 of file JFitToolkit.hh.

  {
    JMatrixNZ V(track, __begin, __end, alpha, sigma);
 
    V.invert();
 
    return getChi2(track, __begin, __end, V);
  }

getChi2() [8/9]

double JFIT::getChi2 ( const JVectorNZ &  Y, const JMatrixNZ &  V, const int  i  )

inline

Determine difference between chi2 with and without hit using full covariance matrix.

Parameters

Y	residuals
V	covariance matrix
i	index of excluded hit

Returns

chi2

Definition at line 218 of file JFitToolkit.hh.

  {
    double chi2 = 0.0;
 
    for (size_t j = 0; j != V.size(); ++j) {
      chi2 += V(i,j) * Y[j];
    }
 
    return chi2*chi2 / V(i,i);
  }

getChi2() [9/9]

template<class JModel_t , class JFit_t , class T >

double JFIT::getChi2 ( const JModel_t &  model, const JFit_t &  fit, T  __begin, T  __end  )

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.

Parameters

model	model
fit	fit function
__begin	begin of data
__end	end of data

Definition at line 251 of file JFitToolkit.hh.

  {
    double chi2 = 0.0;
 
    for (T hit = __begin; hit != __end; ++hit) {
      chi2 += fit(model, *hit);
    }
 
    return chi2;
  }

Variable Documentation

JPP_RECONSTRUCTION_TYPE

const int JFIT::JPP_RECONSTRUCTION_TYPE = 4000

static

Jpp reconstruction type for AAnet.

Definition at line 19 of file JFitApplications.hh.

JKM3DELTAPOS

const int JFIT::JKM3DELTAPOS = 10000

static

This is not a fit, this gives position information only.

Definition at line 21 of file JFitApplications.hh.