Model fits to data. More...

Classes
struct	JQualitySorter
	General purpose sorter of fit results. More...

class	JPointing
	Auxiliary class to compare fit results with respect to a reference direction (e.g. true muon). More...

class	JPosition
	Auxiliary class to compare fit results with respect to a reference position (e.g. true muon). More...

class	JShowerEnergy

class	JAtmosphericMuon
	Auxiliary class to evaluate atmospheric muon hypothesis. More...

class	JHitW0
	Auxiliary class for a hit with background rate value. More...

struct	JMuonEnergyParameters_t
	Data structure for fit parameters. More...

struct	JMuonGandalf
	Wrapper class to make final fit of muon trajectory. More...

struct	JMuonGandalfParameters_t
	Data structure for fit parameters. More...

struct	JMuonParameters
	Muon reconstruction parameters. More...

struct	JMuonPrefit
	Wrapper class to make pre-fit of muon trajectory. More...

struct	JMuonPrefitParameters_t
	Data structure for fit parameters. More...

struct	JMuonSimplex
	Wrapper class to make intermediate fit of muon trajectory. More...

struct	JMuonSimplexParameters_t
	Data structure for fit parameters. More...

struct	JMuonStartParameters_t
	Data structure for fit parameters. More...

class	JShowerBjorkenY
	class to handle the direction fit of the shower reconstruction, mainly dedicated for ORCA More...

struct	JShowerBjorkenYParameters_t
	Data structure for fit parameters. More...

class	JShowerDirectionPrefit
	class to handle the direction fit of the shower reconstruction, mainly dedicated for ORCA More...

struct	JShowerDirectionPrefitParameters_t
	Data structure for fit parameters. More...

class	JShowerEnergyPrefit
	class to handle the third step of the shower reconstruction, mainly dedicated for ORCA More...

struct	JShowerEnergyPrefitParameters_t
	Data structure for fit parameters. More...

class	JShowerFit
	class to handle the direction fit of the shower reconstruction, mainly dedicated for ORCA More...

struct	JShowerFitParameters_t
	Data structure for fit parameters. More...

struct	JShowerParameters
	Shower reconstruction parameters. More...

class	JShowerPointSimplex
	class to handle the second position fit of the shower reconstruction, mainly dedicated for ORCA More...

struct	JShowerPointSimplexParameters_t

class	JShowerPositionFit
	class to handle the second position fit of the shower reconstruction, mainly dedicated for ORCA More...

struct	JShowerPositionFitParameters_t

class	JShowerPrefit
	class to handle first step of the shower reconstruction in ORCA: it reconstructs the shower vertex, intended as the shower brightest point, as the barycenter of the hits More...

struct	JShowerPrefitParameters_t
	Data structure for fit parameters. More...

struct	JStart
	Auxiliary class for start or end point evaluation. More...

Typedefs
typedef JFIT::JFit	JFit

typedef JFIT::JEvt	JEvt

typedef JFIT::JEvent	JEvent

typedef JFIT::JHistory	JHistory

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 N, const int NDF)
	Get quality of fit. 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...

Detailed Description

Model fits to data.

Author: mdejong; mdejong, gmaggi; gmaggi; adomi; adomi, gmaggi

Typedef Documentation

typedef JFIT::JFit JRECONSTRUCTION::JFit

Definition at line 453 of file JReconstruction/JEvt.hh.

typedef JFIT::JEvt JRECONSTRUCTION::JEvt

Definition at line 454 of file JReconstruction/JEvt.hh.

typedef JFIT::JEvent JRECONSTRUCTION::JEvent

Definition at line 282 of file JHistory.hh.

typedef JFIT::JHistory JRECONSTRUCTION::JHistory

Definition at line 283 of file JHistory.hh.

Function Documentation

JPosition3D JRECONSTRUCTION::getPosition ( const JFit & fit )

inline

Get position.

Parameters

fit fit

Returns: position

Definition at line 70 of file JReconstruction/JEvtToolkit.hh.

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

JDirection3D JRECONSTRUCTION::getDirection ( const JFit & fit )

inline

Get direction.

Parameters

fit fit

Returns: direction

Definition at line 82 of file JReconstruction/JEvtToolkit.hh.

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

JAxis3D JRECONSTRUCTION::getAxis ( const JFit & fit )

inline

Get axis.

Parameters

fit fit

Returns: axis

Definition at line 94 of file JReconstruction/JEvtToolkit.hh.

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

JTrack3E JRECONSTRUCTION::getTrack ( const JFit & fit )

inline

Get track.

Parameters

fit fit

Returns: track

Definition at line 106 of file JReconstruction/JEvtToolkit.hh.

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

JFit JRECONSTRUCTION::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 123 of file JReconstruction/JEvtToolkit.hh.

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

JFit JRECONSTRUCTION::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 151 of file JReconstruction/JEvtToolkit.hh.

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

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

inline

Get dot product.

Parameters

first	first fit
second	second fit

Returns: dot product

Definition at line 175 of file JReconstruction/JEvtToolkit.hh.

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

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

inline

Get dot product.

Parameters

fit	fit
dir	direction

Returns: dot product

Definition at line 188 of file JReconstruction/JEvtToolkit.hh.

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

double JRECONSTRUCTION::getQuality	(	const double	chi2,
		const int	N,
		const int	NDF
	)

inline

Get quality of fit.

The larger the quality, the better the fit.

Parameters

chi2	chi2
N	number of hits
NDF	number of degrees of freedom

Returns: quality

Definition at line 203 of file JReconstruction/JEvtToolkit.hh.

   {
     return N  -  0.25 * chi2 / NDF;
   }

double JRECONSTRUCTION::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 217 of file JReconstruction/JEvtToolkit.hh.

   {
     return NDF  -  0.25 * chi2 / NDF;
   }

double JRECONSTRUCTION::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 230 of file JReconstruction/JEvtToolkit.hh.

   {
     return -chi2;
   }

bool JRECONSTRUCTION::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 243 of file JReconstruction/JEvtToolkit.hh.

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

bool JRECONSTRUCTION::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 282 of file JReconstruction/JEvtToolkit.hh.

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

bool JRECONSTRUCTION::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 297 of file JReconstruction/JEvtToolkit.hh.

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

bool JRECONSTRUCTION::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 311 of file JReconstruction/JEvtToolkit.hh.

   {
     return has_history(fit, JMUONPREFIT);
   }

bool JRECONSTRUCTION::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 323 of file JReconstruction/JEvtToolkit.hh.

   {
     return has_history(fit, JMUONSIMPLEX);
   }

bool JRECONSTRUCTION::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 335 of file JReconstruction/JEvtToolkit.hh.

   {
     return has_history(fit, JMUONGANDALF);
   }

bool JRECONSTRUCTION::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 347 of file JReconstruction/JEvtToolkit.hh.

   {
     return has_history(fit, JMUONENERGY);
   }

bool JRECONSTRUCTION::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 359 of file JReconstruction/JEvtToolkit.hh.

   {
     return has_history(fit, JMUONSTART);
   }

bool JRECONSTRUCTION::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 371 of file JReconstruction/JEvtToolkit.hh.

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

bool JRECONSTRUCTION::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 383 of file JReconstruction/JEvtToolkit.hh.

   {
     return has_history(fit, JSHOWERPREFIT);
   }

bool JRECONSTRUCTION::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 395 of file JReconstruction/JEvtToolkit.hh.

   {
     return has_history(fit, JSHOWERPOSITIONFIT);
   }

bool JRECONSTRUCTION::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 407 of file JReconstruction/JEvtToolkit.hh.

   {
     return has_history(fit, JSHOWERCOMPLETEFIT);
   }

bool JRECONSTRUCTION::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 419 of file JReconstruction/JEvtToolkit.hh.

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

template<class JTrackSelector_t >

bool JRECONSTRUCTION::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 434 of file JReconstruction/JEvtToolkit.hh.

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

bool JRECONSTRUCTION::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 446 of file JReconstruction/JEvtToolkit.hh.

   {
     return has_reconstructed_track(evt, has_muon_fit);
   }

bool JRECONSTRUCTION::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 458 of file JReconstruction/JEvtToolkit.hh.

   {
     return has_reconstructed_track(evt, has_shower_fit);
   }

template<class JTrackSelector_t , class JQualitySorter_t >

const JFit& JRECONSTRUCTION::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 475 of file JReconstruction/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.");
   }

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

inline

Get best reconstructed muon.

Parameters

evt event

Returns: track

Definition at line 500 of file JReconstruction/JEvtToolkit.hh.

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

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

inline

Get best reconstructed shower.

Parameters

evt event

Returns: track

Definition at line 512 of file JReconstruction/JEvtToolkit.hh.

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

Classes

Typedefs

Functions

Detailed Description

Typedef Documentation

Function Documentation