JRECONSTRUCTION Namespace Reference

Support. More...

Classes
struct	event_type
	Auxiliary data structure to store data and fit in memory. More...
struct	JARCAMuonReconstruction
	ARCA muon reconstruction. More...
class	JAtmosphericMuon
	Auxiliary class to evaluate atmospheric muon hypothesis. More...
struct	JEditor
	Auxiliary class for editing time offset. More...
class	JEnergyCorrection
	Auxiliary class for correction of energy determined by JEnergy.cc. More...
struct	JEventSelector
	Event selector. More...
class	JHitW0
	Auxiliary class for a hit with background rate value. More...
struct	JModuleL0
	Auxiliary class for handling module response. More...
class	JMultiThreadedReconstruction
	Thread pool for event-by-event reconstruction. More...
class	JMuonEnergy
	Auxiliary class to to determine muon energy. More...
struct	JMuonEnergyParameters_t
	Data structure for fit parameters. More...
struct	JMuonFeatures
	Wrapper class to add features after the final fit of muon trajectory. 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	JORCAMuonReconstruction
	ORCA muon reconstruction. More...
struct	JORCAReconstruction
	ORCA reconstruction. More...
struct	JORCAShowerReconstruction
	ORCA shower reconstruction. More...
struct	JPerth
	Thread pool for fits to data. More...
struct	JPerth_t
	Auxiliary data structure for chi2 function object. 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. More...
struct	JQualitySorter
	General purpose sorter of fit results. More...
struct	JRange_t
	Type definition of range. More...
class	JReconstructionDictionary
	Simple wrapper class around JROOT::JRootDictionary so that other classes are included by construction,. 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	JShowerEnergy
	Auxiliary class to compare fit results with respect to a reference energy. 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...
struct	JWeight
	Auxiliary data structure to get weight of given fit. 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
typedef std::map< int, buffer_type >	map_type
	identifier -> hits
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
void	copy (const JFIT::JEvt::const_iterator __begin, const JFIT::JEvt::const_iterator __end, Evt &out)
	Copy tracks.
JPosition3D	getPosition (const JFit &fit)
	Get position.
JDirection3D	getDirection (const JFit &fit)
	Get direction.
JVertex3D	getVertex (const JFit &fit)
	Get vertex.
JAxis3D	getAxis (const JFit &fit)
	Get axis.
JTrack3E	getTrack (const JFit &fit)
	Get track.
JShower3E	getShower (const JFit &fit)
	Get shower.
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.
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.
JFit	getFit (const JHistory &history, const JShower3D &shower, const double Q, const int NDF, const double energy=0.0, const int status=SINGLE_STAGE)
	Get fit.
double	getDot (const JFit &first, const JFit &second)
	Get dot product.
double	getDot (const JFit &fit, const JDirection3D &dir)
	Get dot product.
double	getAngle (const JFit &first, const JFit &second)
	Get space angle.
double	getAngle (const JFit &fit, const JDirection3D &dir)
	Get space angle.
double	getQuality (const double chi2, const int N, const int NDF)
	Get quality of fit.
double	getQuality (const double chi2, const int NDF)
	Get quality of fit.
double	getQuality (const double chi2)
	Get quality of fit.
bool	qualitySorter (const JFit &first, const JFit &second)
	Comparison of fit results.
bool	has_history (const JFit &fit, const int type)
	Test whether given fit has specified history.
bool	has_history (const JFit &fit, const JRange< int > &range)
	Test whether given fit has specified history.
bool	has_muon_prefit (const JFit &fit)
	Test whether given fit has muon prefit in history.
bool	has_muon_simplex (const JFit &fit)
	Test whether given fit has muon simplex in history.
bool	has_muon_gandalf (const JFit &fit)
	Test whether given fit has muon gandalf in history.
bool	has_muon_energy (const JFit &fit)
	Test whether given fit has muon energy in history.
bool	has_muon_start (const JFit &fit)
	Test whether given fit has muon start in history.
bool	has_muon_fit (const JFit &fit)
	Test whether given fit has muon fit in history.
bool	has_shower_prefit (const JFit &fit)
	Test whether given fit has shower prefit in history.
bool	has_shower_positionfit (const JFit &fit)
	Test whether given fit has shower position fit in history.
bool	has_shower_completefit (const JFit &fit)
	Test whether given fit has shower complete fit in history.
bool	has_shower_fit (const JFit &fit)
	Test whether given fit has shower fit in history.
template<class JTrackSelector_t >
bool	has_reconstructed_track (const JEvt &evt, JTrackSelector_t selector)
	Test whether given event has a track according selection.
bool	has_reconstructed_muon (const JEvt &evt)
	Test whether given event has a track with muon reconstruction.
bool	has_reconstructed_shower (const JEvt &evt)
	Test whether given event has a track with shower reconstruction.
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.
const JFit &	get_best_reconstructed_muon (const JEvt &evt)
	Get best reconstructed muon.
const JFit &	get_best_reconstructed_shower (const JEvt &evt)
	Get best reconstructed shower.
JEvt::iterator	gridify (JEvt::iterator __begin, JEvt::iterator __end, const int N)
	Gridify set of fits.
bool	filter (const JDAQEvent &tev, const JEvt &evt, const Evt *const pE)
	Event selection.

Variables
JRECONSTRUCTION::JWeight	getWeight
const char *const	h2_t = "h2"
	Name of histogram with results from JMuonCompass.cc.

Detailed Description

Support.

Model fits to data.

Author

mdejong

mdejong, azegarelli, scelli

eleguirriec, mdejong

mdejong, gmaggi

gmaggi

mdejong, gmaggi, azegarelli

adomi

adomi, vcarretero

adomi, gmaggi, vcarretero

adomi, gmaggi

Typedef Documentation

JFit

typedef JFIT::JFit JRECONSTRUCTION::JFit

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

JEvt

typedef JFIT::JEvt JRECONSTRUCTION::JEvt

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

JEvent

typedef JFIT::JEvent JRECONSTRUCTION::JEvent

Definition at line 454 of file JHistory.hh.

JHistory

typedef JFIT::JHistory JRECONSTRUCTION::JHistory

Definition at line 455 of file JHistory.hh.

buffer_type

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

hits

Definition at line 70 of file JPerth.cc.

map_type

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

identifier -> hits

Definition at line 71 of file JPerth.cc.

data_type

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

Definition at line 81 of file JPerth.cc.

input_type

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

Definition at line 82 of file JPerth.cc.

JRegressorStorage_t

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

Definition at line 131 of file JPerth.cc.

JRegressor_t

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

Definition at line 132 of file JPerth.cc.

Enumeration Type Documentation

JFitStatus_t

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

copy()

void JRECONSTRUCTION::copy ( const JFIT::JEvt::const_iterator __begin, const JFIT::JEvt::const_iterator __end, Evt & out )

inline

Copy tracks.

Parameters

__begin	begin of input tracks
__end	end of input tracks
out	output tracks

Definition at line 35 of file JEvtSupportkit.hh.

  {
    using namespace std;
    using namespace JPP;
 
    struct : public set<JUUID> {
 
      inline int get_index(const JUUID& element) const
      {
        const_iterator i = this->find(element);
 
        if (i != this->end())
          return std::distance(this->begin(), i);
        else
          return -1;
      }
    } uuid;
 
    for (JFIT::JEvt::const_iterator fit = __begin; fit != __end; ++fit) {
      uuid.insert(fit->getUUID());
    }
 
    for (JFIT::JEvt::const_iterator fit = __begin; fit != __end; ++fit) {     // import reconstruction data
 
      Trk trk;
 
      trk.id        = uuid.get_index(fit->getUUID());
      trk.pos       = Vec(fit->getX(),  fit->getY(),  fit->getZ());
      trk.dir       = Vec(fit->getDX(), fit->getDY(), fit->getDZ());
      trk.t         = fit->getT();
      trk.E         = fit->getE();
      trk.lik       = fit->getQ();
      trk.rec_type  = JPP_RECONSTRUCTION_TYPE;
      trk.status    = (fit->getStatus() == COMPLETE_CHAIN ? TRK_ST_FINALSTATE : TRK_ST_UNDEFINED);
 
      if (fit->hasParentUUID()) {
        trk.mother_id = uuid.get_index(fit->getParentUUID());
      }
 
      for (JHistory::const_iterator i = fit->getHistory().begin(); i != fit->getHistory().end(); ++i) {
        trk.rec_stages.push_back(i->type);
      }
 
      for (int i = 0; i != fit->getN(); ++i) {
        trk.fitinf.push_back(fit->getW(i));
      }
 
      trk.error_matrix = fit->getV();
 
      out.trks.push_back(trk);
    }
  }

getPosition()

JPosition3D JRECONSTRUCTION::getPosition ( const JFit & fit )

inline

Get position.

Parameters

fit

fit

Returns

position

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

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

getDirection()

JDirection3D JRECONSTRUCTION::getDirection ( const JFit & fit )

inline

Get direction.

Parameters

fit

fit

Returns

direction

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

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

getVertex()

JVertex3D JRECONSTRUCTION::getVertex ( const JFit & fit )

inline

Get vertex.

Parameters

fit

fit

Returns

vertex

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

  {
    return JVertex3D(getPosition(fit), fit.getT());
  }

getAxis()

JAxis3D JRECONSTRUCTION::getAxis ( const JFit & fit )

inline

Get axis.

Parameters

fit

fit

Returns

axis

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

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

getTrack()

JTrack3E JRECONSTRUCTION::getTrack ( const JFit & fit )

inline

Get track.

Parameters

fit

fit

Returns

track

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

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

getShower()

JShower3E JRECONSTRUCTION::getShower ( const JFit & fit )

inline

Get shower.

Parameters

fit

fit

Returns

shower

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

  {
    return JShower3E(JShower3D(getVertex(fit), getDirection(fit)), fit.getE());
  }

getFit() [1/3]

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 231 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);
  }

getFit() [2/3]

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 259 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);
  }

getFit() [3/3]

JFit JRECONSTRUCTION::getFit ( const JHistory & history, const JShower3D & shower, const double Q, const int NDF, const double energy = 0.0, const int status = SINGLE_STAGE )

inline

Get fit.

Parameters

history	history
shower	shower
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 287 of file JReconstruction/JEvtToolkit.hh.

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

getDot() [1/2]

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 310 of file JReconstruction/JEvtToolkit.hh.

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

getDot() [2/2]

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

inline

Get dot product.

Parameters

fit	fit
dir	direction

Returns

dot product

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

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

getAngle() [1/2]

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

inline

Get space angle.

Parameters

first	first fit
second	second fit

Returns

angle [deg]

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

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

getAngle() [2/2]

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

inline

Get space angle.

Parameters

fit	fit
dir	direction

Returns

angle [deg]

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

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

getQuality() [1/3]

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 364 of file JReconstruction/JEvtToolkit.hh.

  {
    return N  -  0.25 * chi2 / NDF;
  }

getQuality() [2/3]

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 378 of file JReconstruction/JEvtToolkit.hh.

  {
    return NDF  -  0.25 * chi2 / NDF;
  }

getQuality() [3/3]

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 391 of file JReconstruction/JEvtToolkit.hh.

  {
    return -chi2;
  }

qualitySorter()

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 404 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();
  }

has_history() [1/2]

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 443 of file JReconstruction/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 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 458 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();
  }

has_muon_prefit()

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 472 of file JReconstruction/JEvtToolkit.hh.

  {
    return has_history(fit, JMUONPREFIT);
  }

has_muon_simplex()

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 484 of file JReconstruction/JEvtToolkit.hh.

  {
    return has_history(fit, JMUONSIMPLEX);
  }

has_muon_gandalf()

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 496 of file JReconstruction/JEvtToolkit.hh.

  {
    return has_history(fit, JMUONGANDALF);
  }

has_muon_energy()

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 508 of file JReconstruction/JEvtToolkit.hh.

  {
    return has_history(fit, JMUONENERGY);
  }

has_muon_start()

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 520 of file JReconstruction/JEvtToolkit.hh.

  {
    return has_history(fit, JMUONSTART);
  }

has_muon_fit()

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 532 of file JReconstruction/JEvtToolkit.hh.

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

has_shower_prefit()

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 544 of file JReconstruction/JEvtToolkit.hh.

  {
    return has_history(fit, JSHOWERPREFIT);
  }

has_shower_positionfit()

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 556 of file JReconstruction/JEvtToolkit.hh.

  {
    return has_history(fit, JSHOWERPOSITIONFIT);
  }

has_shower_completefit()

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 568 of file JReconstruction/JEvtToolkit.hh.

  {
    return has_history(fit, JSHOWERCOMPLETEFIT);
  }

has_shower_fit()

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 580 of file JReconstruction/JEvtToolkit.hh.

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

has_reconstructed_track()

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

Parameters

evt	event
selector	track selector

Returns

true if at least one corresponding track; else false

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

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

has_reconstructed_muon()

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 607 of file JReconstruction/JEvtToolkit.hh.

  {
    return has_reconstructed_track(evt, has_muon_fit);
  }

has_reconstructed_shower()

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 619 of file JReconstruction/JEvtToolkit.hh.

  {
    return has_reconstructed_track(evt, has_shower_fit);
  }

get_best_reconstructed_track()

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 comparator(const Trk&, const Trk&);.

Parameters

evt	event
selector	track selector
comparator	track comparator

Returns

track

Definition at line 636 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.");
  }

get_best_reconstructed_muon()

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

inline

Get best reconstructed muon.

Parameters

evt

event

Returns

track

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

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

get_best_reconstructed_shower()

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

inline

Get best reconstructed shower.

Parameters

evt

event

Returns

track

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

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

gridify()

JEvt::iterator JRECONSTRUCTION::gridify ( JEvt::iterator __begin, JEvt::iterator __end, const int N )

inline

Gridify set of fits.

Parameters

__begin	begin of fits
__end	end of fits
N	number of directions in grid

Returns

end of fits sorted by quality

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

  {
    using namespace std;
    using namespace JPP;
 
    static const double MINIMAL_ANGLE_DEG  =  60.0;
 
    // estimate grid angle
 
    double angle_deg = 10.0 * sqrt((double) 240 / (double) N);
 
    if (angle_deg > MINIMAL_ANGLE_DEG) {
      angle_deg = MINIMAL_ANGLE_DEG;
    }
 
    const JOmega3D ps(angle_deg * PI/180.0);
 
    JEvt::iterator out = __begin;
 
    for (JOmega3D_t::const_iterator dir = ps.begin(); dir != ps.end() && out != __end; ++dir) {
 
      const JPointing pointing(*dir);
 
      sort(out, __end, pointing);
 
      JEvt::iterator p = out;
 
      for (JEvt::iterator i = out; i != __end && pointing.getAngle(*i) <= angle_deg; ++i) {
        if (qualitySorter(*i, *p)) {
          p = i;
        }
      }
 
      swap(*(out++), *p);
    }
 
    sort(__begin, out, qualitySorter);
 
    return out;
  }

filter()

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;
  }

Variable Documentation

getWeight

JRECONSTRUCTION::JWeight JRECONSTRUCTION::getWeight

static

h2_t

const char* const JRECONSTRUCTION::h2_t = "h2"

Name of histogram with results from JMuonCompass.cc.

Definition at line 12 of file JMuonCompass.hh.