JRECONSTRUCTION Namespace Reference

Model fits to data. More...

Classes
struct	JQualitySorter
	General purpose sorter of fit results. More...
struct	JHitL0Comparator
	Auxiliary class for permutations of L0 hits. More...
struct	JHitL1Comparator
	Auxiliary class for permutations of L1 hits. 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 76 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 88 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 100 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 112 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 129 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 157 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 181 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 194 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 209 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 223 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 236 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 249 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 288 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 303 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 317 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 329 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 341 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 353 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 365 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 377 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 389 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 401 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 413 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 425 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 440 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 452 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 464 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 481 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 506 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 518 of file JReconstruction/JEvtToolkit.hh.

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