Model fits to data. More...

Classes
struct	JEventSelector
	Event selector. More...

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...

class	JMuonEnergy
	Auxiliary class to to determine muon energy. 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...

class	JMuonStart
	Auxiliary class to determine start and end position of muon trajectory. More...

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

struct	event_type
	Auxiliary data structure to store data and fit in memory. More...

struct	JStringEditor
	Auxiliary class to edit time offset of data per string. More...

struct	JPerth
	Thread pool for fits to data. More...

struct	JPerth_t
	Auxiliary data structure for chi2 function object. More...

struct	JRange_t
	Type definition of range. 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

typedef std::vector< JHitW0 >	buffer_type
	hits More...

typedef std::map< int, buffer_type >	map_type
	string -> hits More...

typedef std::vector< event_type >	data_type

typedef JLANG::JSTDObjectReader< const event_type >	input_type

typedef JFIT::JRegressorStorage < JFIT::JLine3Z, JFIT::JGandalf >	JRegressorStorage_t

typedef JFIT::JRegressor < JFIT::JLine3Z, JFIT::JGandalf >	JRegressor_t

Enumerations
enum	JFitStatus_t { ERROR = -100, INSUFFICIENT_NUMBER_OF_HITS, INSUFFICIENT_NUMBER_OF_DOMS, UNDEFINED = -1, SINGLE_STAGE = 0, INCOMPLETE_CHAIN = +1, COMPLETE_CHAIN = +100 }

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=SINGLE_STAGE)
	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=SINGLE_STAGE)
	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...

double	getProbability (const JModule &module, const JDAQSummaryFrame &frame, const JRateL1_t &R_Hz, const double T_ns, const std::multiset< int > &top)
	Get probability of given response in optical module due to random background. More...

bool	filter (const JDAQEvent &tev, const JEvt &evt, const Evt *const pE)
	Event selection. More...

Detailed Description

Model fits to data.

Author: mdejong; mdejong, azegarelli, scelli; mdejong, gmaggi; gmaggi; mdejong, gmaggi, azegarelli; adomi; adomi, gmaggi

Typedef Documentation

typedef JFIT::JFit JRECONSTRUCTION::JFit

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

typedef JFIT::JEvt JRECONSTRUCTION::JEvt

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

typedef JFIT::JEvent JRECONSTRUCTION::JEvent

Definition at line 353 of file JHistory.hh.

typedef JFIT::JHistory JRECONSTRUCTION::JHistory

Definition at line 354 of file JHistory.hh.

typedef std::vector<JHitW0> JRECONSTRUCTION::buffer_type

hits

Definition at line 67 of file JPerth.cc.

typedef std::map<int, buffer_type> JRECONSTRUCTION::map_type

string -> hits

Definition at line 68 of file JPerth.cc.

typedef std::vector<event_type> JRECONSTRUCTION::data_type

Definition at line 78 of file JPerth.cc.

typedef JLANG::JSTDObjectReader<const event_type> JRECONSTRUCTION::input_type

Definition at line 79 of file JPerth.cc.

typedef JFIT::JRegressorStorage<JFIT::JLine3Z, JFIT::JGandalf> JRECONSTRUCTION::JRegressorStorage_t

Definition at line 128 of file JPerth.cc.

typedef JFIT::JRegressor<JFIT::JLine3Z, JFIT::JGandalf> JRECONSTRUCTION::JRegressor_t

Definition at line 129 of file JPerth.cc.

Enumeration Type Documentation

enum JRECONSTRUCTION::JFitStatus_t

Enumerator
ERROR
INSUFFICIENT_NUMBER_OF_HITS
INSUFFICIENT_NUMBER_OF_DOMS
UNDEFINED
SINGLE_STAGE
INCOMPLETE_CHAIN
COMPLETE_CHAIN

Definition at line 16 of file JFitStatus.hh.

                     {
 
     ERROR             = -100,
     INSUFFICIENT_NUMBER_OF_HITS,
     INSUFFICIENT_NUMBER_OF_DOMS,
     UNDEFINED         = -1,
     SINGLE_STAGE      =  0,
     INCOMPLETE_CHAIN  = +1,
     COMPLETE_CHAIN    = +100
   };

Function Documentation

JPosition3D JRECONSTRUCTION::getPosition ( const JFit & fit )

inline

Get position.

Parameters

fit fit

Returns: position

Definition at line 71 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 83 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 95 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 107 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 = `SINGLE_STAGE`
	)

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 124 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 = `SINGLE_STAGE`
	)

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 152 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 176 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 189 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 204 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 218 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 231 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 244 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 283 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 298 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 312 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 324 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 336 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 348 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 360 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 372 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 384 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 396 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 408 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 420 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 435 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 447 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 459 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 476 of file JReconstruction/JEvtToolkit.hh.

   {
     JEvt::const_iterator p = std::find_if(evt.begin(), evt.end(), 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 501 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 513 of file JReconstruction/JEvtToolkit.hh.

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

double JRECONSTRUCTION::getProbability	(	const JModule &	module,
		const JDAQSummaryFrame &	frame,
		const JRateL1_t &	R_Hz,
		const double	T_ns,
		const std::multiset< int > &	top
	)

inline

Get probability of given response in optical module due to random background.

Parameters

module	module
frame	summary frame
R_Hz	multiples rates [Hz]
T_ns	time window [ns]
top	list of indentifiers of PMTs with hit in given module

Returns: probability

Definition at line 75 of file JMuonStart.hh.

   {
     using namespace std;
     using namespace JPP;
     using namespace KM3NETDAQ;
 
     size_t N = 0;         // number of active PMTs
     size_t M = 0;         // multiplicity
     double R = 0.0;       // total rate
 
     for (size_t i = 0; i != module.size(); ++i) {
 
       if (getDAQStatus(frame, module, i)            &&
           getPMTStatus(frame, module, i)            &&
           !module.getPMT(i).has(PMT_DISABLE)) {
 
         N += 1;
         M += top.count(i);
         R += frame.getRate(i);
       }
     }
 
     if (N > 0)
       return JFIT::getProbability(N, M, JK40Rates(R/N, R_Hz), T_ns);
     else
       return (M == 0 ? 1.0 : 0.0);
   }

bool JRECONSTRUCTION::filter	(	const JDAQEvent &	tev,
		const JEvt &	evt,
		const Evt *const	pE
	)

Event selection.

Parameters

tev	DAQ event
evt	input event
pE	pointer to Monte Carlo event

Returns: true if selected; else false

Definition at line 21 of file JReconstruction/event_selector.cc.

   {
     using namespace std;
     using namespace JPP;
     /*
     if (!evt.empty() && evt[0].getQ() >= 100) {
 
       if (pE != NULL) {
 
         Trk muon;
 
         for (const auto& t1 : pE->mc_trks) {
           if (is_muon(t1)) {
             if (t1.E > muon.E) {
               muon = t1;
             }
           }
         }
 
         return (is_muon(muon) && getAngle(getDirection(event[0]), getDirection(muon)) > 5.0);
       }
     }
     */
     return true;
   }

Classes

Typedefs

Enumerations

Functions

Detailed Description

Typedef Documentation

Enumeration Type Documentation

Function Documentation