Extensions to AAnet data format. More...

Classes
struct	can
	The cylinder used for photon tracking. More...

struct	coord_origin
	Coordinate origin. More...

struct	cut
	General purpose class of phase space generation. More...

struct	cut_in
	Phase space of incoming particle. More...

struct	cut_nu
	Phase space of incident neutrino. More...

struct	cut_primary
	Phase space of primary particle. More...

struct	cut_seamuon
	Phase space of atmospheric muon generation. More...

struct	DAQ
	Livetime of DAQ data. More...

struct	detector
	Detector file. More...

struct	end_event
	End of event record. More...

struct	fixedcan
	The fixe cylinder used for photon tracking. More...

struct	flux
	Neutrino flux. More...

struct	generator
	Description of Monte Carlo event generation applications. More...

struct	genhencut
	Phase space for incident neutrino. More...

struct	genvol
	Neutrino vertex volume. More...

struct	getCorsikaHeader
	Match header for Corsika. More...

struct	getDAQHeader
	Match header for DAQ. More...

struct	getGenhenHeader
	Match header for genhen. More...

struct	getGenieHeader
	Match header for genie. More...

struct	getGSeaGenHeader
	Match header for gseagen. More...

struct	getKM3Header
	Match header for km3. More...

struct	getKM3SimHeader
	Match header for km3. More...

struct	getMUPAGEHeader
	Match header for MUPAGE. More...

struct	getSireneHeader
	Match header for JSirene.cc. More...

struct	has_history
	Auxiliary class to test whether given track has specified history. More...

struct	JEvtEvaluator
	Auxiliary class to determine value of Evt objects. More...

struct	JHead
	Monte Carlo run header. More...

struct	JHit_t
	Auxiliary class to set-up Hit. More...

struct	JHits_t
	Auxiliary data structure for list of hits with hit merging capability. More...

struct	JParticle
	Auxiliary class to handle particle name, codes and mass. More...

struct	JPDB
	Collection of particles. More...

struct	JTrk_t
	Auxiliary class to set-up Trk. More...

struct	JWeightEvent
	Low-level interface for event weighing. More...

struct	JWeightEventHelper
	Helper class for event weighing. More...

struct	JWeightFileScanner
	Template file scanner with event weight. More...

struct	JWeightMupage
	Implementation of event weighing for MUPAGE data. More...

struct	livetime
	Normalisation of MUPAGE events. More...

struct	muon_desc_file
	Muon descriptor file. More...

struct	norma
	Normlisation of CORSIKA events. More...

struct	physics
	Generator for neutrino interaction. More...

struct	primary
	Primary particle. More...

struct	quality_sorter
	Reconstruction type dependent comparison of track quality. More...

struct	seabottom
	The bottom of the sea. More...

struct	simul
	Generator for simulation. More...

struct	spectrum
	Neutrino energy spectrum. More...

struct	start_run
	Start of run record. More...

struct	String
	General purpose string class. More...

struct	target
	Target. More...

struct	tgen
	Scaling time of neutrino interaction generators. More...

struct	XSecFile
	Neutrino cross section file. More...

Enumerations
enum	JHitType_t { HIT_TYPE_MUON_DIRECT = +5, HIT_TYPE_MUON_SCATTERED = -5, HIT_TYPE_DELTARAYS_DIRECT = +4, HIT_TYPE_DELTARAYS_SCATTERED = -4, HIT_TYPE_BREMSSTRAHLUNG_DIRECT = +3, HIT_TYPE_BREMSSTRAHLUNG_SCATTERED = -3, HIT_TYPE_SHOWER_DIRECT = +99, HIT_TYPE_SHOWER_SCATTERED = -99, HIT_TYPE_NOISE = -1, HIT_TYPE_UNKNOWN = 0 }
	Enumeration of hit types based on km3 codes. More...

enum	JGeant4Type_t { GEANT4_TYPE_PHOTON = 1, GEANT4_TYPE_ANTIELECTRON = 2, GEANT4_TYPE_ELECTRON = 3, GEANT4_TYPE_NEUTRINO = 4, GEANT4_TYPE_ANTIMUON = 5, GEANT4_TYPE_MUON = 6, GEANT4_TYPE_NEUTRAL_PION = 7, GEANT4_TYPE_PION_PLUS = 8, GEANT4_TYPE_PION_MINUS = 9, GEANT4_TYPE_KAON_LONG = 10, GEANT4_TYPE_KAON_PLUS = 11, GEANT4_TYPE_KAON_MINUS = 12, GEANT4_TYPE_NEUTRON = 13, GEANT4_TYPE_PROTON = 14, GEANT4_TYPE_ANTIPROTON = 15, GEANT4_TYPE_KAON_SHORT = 16, GEANT4_TYPE_ETA = 17, GEANT4_TYPE_LAMBDA = 18, GEANT4_TYPE_SIGMA_PLUS = 19, GEANT4_TYPE_NEUTRAL_SIGMA = 20, GEANT4_TYPE_SIGMA_MINUS = 21, GEANT4_TYPE_NEUTRAL_XI = 22, GEANT4_TYPE_XI_MINUS = 23, GEANT4_TYPE_OMEGA_MINUS = 24, GEANT4_TYPE_ANTINEUTRON = 25, GEANT4_TYPE_ANTILAMBDA = 26, GEANT4_TYPE_ANTISIGMA_MINUS = 27, GEANT4_TYPE_NEUTRAL_ANTISIGMA = 28, GEANT4_TYPE_ANTISIGMA_PLUS = 29, GEANT4_TYPE_NEUTRAL_ANTIXI = 30, GEANT4_TYPE_ANTIXI_PLUS = 31, GEANT4_TYPE_ANTIOMEGA_PLUS = 32, GEANT4_TYPE_DEUTERON = 45, GEANT4_TYPE_TRITON = 46, GEANT4_TYPE_ALPHA = 47, GEANT4_TYPE_GEANTINO = 48, GEANT4_TYPE_HE3 = 49, GEANT4_TYPE_ANTITAU = 33, GEANT4_TYPE_TAU = 34 }
	Enumeration of hit types based on Geant4 codes, for compatbility with KM3Sim. More...

enum	JTrackType_t { TRACK_TYPE_ELECTRON = 11, TRACK_TYPE_NUE = 12, TRACK_TYPE_MUON = 13, TRACK_TYPE_NUMU = 14, TRACK_TYPE_TAU = 15, TRACK_TYPE_NUTAU = 16, TRACK_TYPE_PHOTON = 22, TRACK_TYPE_NEUTRAL_PION = 111, TRACK_TYPE_CHARGED_PION_PLUS = 211, TRACK_TYPE_PION_PLUS = 211, TRACK_TYPE_K_LONG = 130, TRACK_TYPE_K_SHORT = 310, TRACK_TYPE_K_PLUS = 321, TRACK_TYPE_PROTON = 2212, TRACK_TYPE_NEUTRON = 2112, TRACK_TYPE_LAMBDA = 3122, TRACK_TYPE_SIGMA_PLUS = 3222, TRACK_TYPE_NEUTRAL_SIGMA = 3212, TRACK_TYPE_SIGMA_MINUS = 3112, TRACK_TYPE_NEUTRAL_XI = 3322, TRACK_TYPE_XI_MINUS = 3312, TRACK_TYPE_OMEGA_MINUS = 3334, TRACK_TYPE_ANTIELECTRON = -11, TRACK_TYPE_ANTINUE = -12, TRACK_TYPE_ANTIMUON = -13, TRACK_TYPE_ANTINUMU = -14, TRACK_TYPE_ANTITAU = -15, TRACK_TYPE_ANTINUTAU = -16, TRACK_TYPE_NEUTRAL_ANTIPION = -111, TRACK_TYPE_CHARGED_PION_MINUS = -211, TRACK_TYPE_PION_MINUS = -211, TRACK_TYPE_ANTIK_LONG = -130, TRACK_TYPE_ANTIK_SHORT = -310, TRACK_TYPE_K_MINUS = -321, TRACK_TYPE_ANTIPROTON = -2212, TRACK_TYPE_ANTINEUTRON = -2112, TRACK_TYPE_ANTILAMBDA = -3122, TRACK_TYPE_ANTISIGMA_PLUS = -3222, TRACK_TYPE_ANTINEUTRAL_SIGMA = -3212, TRACK_TYPE_ANTISIGMA_MINUS = -3112, TRACK_TYPE_ANTINEUTRAL_XI = -3322, TRACK_TYPE_ANTIXI_MINUS = -3312, TRACK_TYPE_ANTIOMEGA_MINUS = -3334 }
	Enumeration of track types based on PDG codes. More...

Functions
double	getTime (const Hit &hit)
	Get true time of hit. More...

double	getNPE (const Hit &hit)
	Get true charge of hit. More...

bool	is_noise (const Hit &hit)
	Verify hit origin. More...

JTimeRange	getTimeRange (const Evt &event)
	Get time range (i.e. More...

JTimeRange	getTimeRange (const Evt &event, const JTimeRange &T_ns)
	Get time range (i.e. More...

JPosition3D	getPosition (const Vec &v)
	Get position. More...

JPosition3D	getPosition (const Trk &track)
	Get position. More...

JDirection3D	getDirection (const Vec &v)
	Get direction. More...

JDirection3D	getDirection (const Trk &track)
	Get direction. More...

JAxis3D	getAxis (const Trk &track)
	Get axis. More...

JTrack3E	getTrack (const Trk &track)
	Get track. More...

JVertex3D	getVertex (const Trk &track)
	Get vertex. More...

JTransformation3D	getTransformation (const Trk &track)
	Get transformation. More...

JHitL0	getHit (const Hit &hit)
	Get transformation. More...

JHitL0	getHit (const Hit &hit, const JModuleRouter &router)
	Get transformation. More...

JHitL0	getHit (const Hit &hit, const JPMTRouter &router)
	Get transformation. More...

bool	is_photon (const Trk &track)
	Test whether given track is a photon or neutral pion. More...

bool	is_neutrino (const Trk &track)
	Test whether given track is a neutrino. More...

bool	is_electron (const Trk &track)
	Test whether given track is a (anti-)electron. More...

bool	is_muon (const Trk &track)
	Test whether given track is a (anti-)muon. More...

bool	is_tau (const Trk &track)
	Test whether given track is a (anti-)tau. More...

bool	is_pion (const Trk &track)
	Test whether given track is a charged pion. More...

bool	is_proton (const Trk &track)
	Test whether given track is a (anti-)proton. More...

bool	is_neutron (const Trk &track)
	Test whether given track is a (anti-)neutron. More...

bool	is_hadron (const Trk &track)
	Test whether given track is a hadron. More...

bool	has_particleID (const Trk &track, const int type)
	Test whether given track corresponds to given particle type. More...

bool	has_neutrino (const Evt &evt)
	Test whether given event has an incoming neutrino. More...

const Trk &	get_neutrino (const Evt &evt)
	Get incoming neutrino. More...

int	count_electrons (const Evt &evt)
	Count the number of electrons in a given event. More...

bool	has_electron (const Evt &evt)
	Test whether given event has an electron. More...

const Trk &	get_electron (const Evt &evt)
	Get first electron from the event tracklist. More...

bool	from_electron (const Hit &hit)
	Test whether given hit was produced by an electron. More...

int	count_muons (const Evt &evt)
	Count the number of muons in a given event. More...

bool	has_muon (const Evt &evt)
	Test whether given event has a muon. More...

const Trk &	get_muon (const Evt &evt)
	Get first muon from the event tracklist. More...

bool	from_muon (const Hit &hit)
	Test whether given hit was produced by a muon. More...

int	count_taus (const Evt &evt)
	Count the number of taus in a given event. More...

bool	has_tau (const Evt &evt)
	Test whether given event has a tau. More...

const Trk &	get_tau (const Evt &evt)
	Get first tau from the event tracklist. More...

bool	from_tau (const Hit &hit)
	Test whether given hit was produced by a tau. More...

int	count_hadrons (const Evt &evt)
	Count the number of hadrons in a given event (not including neutral pions) More...

bool	from_hadron (const Hit &hit)
	Test whether given hit was produced by a hadronic shower. More...

Trk	get_hadronic_cascade (const Evt &evt)
	Auxiliary function to get average direction and total energy of a hadronic cascade. More...

void	add_time_offset (Evt &evt, const double tOff)
	Add time offset to mc event; according to current definition, the absolute time of the event is defined by the track "t" attribute; this could change in the future if the global attribute mc_t is assigned to this purpose. More...

bool	has_muon_prefit (const Trk &track)
	Test whether given track has muon prefit in history. More...

bool	has_muon_simplex (const Trk &track)
	Test whether given track has muon simplex fit in history. More...

bool	has_muon_gandalf (const Trk &track)
	Test whether given track has muon gandalf fit in history. More...

bool	has_muon_energy (const Trk &track)
	Test whether given track has muon energy fit in history. More...

bool	has_muon_start (const Trk &track)
	Test whether given track has muon start fit in history. More...

bool	has_muon_fit (const Trk &track)
	Test whether given track has default muon fit in history. More...

bool	has_shower_prefit (const Trk &track)
	Test whether given track has shower prefit in history. More...

bool	has_shower_positionfit (const Trk &track)
	Test whether given track has shower position fit in history. More...

bool	has_shower_completefit (const Trk &track)
	Test whether given track has shower complete fit in history. More...

bool	has_shower_fit (const Trk &track)
	Test whether given track has default shower fit in history. More...

bool	has_aashower_fit (const Trk &track)
	Test whether given track has default shower fit in history. More...

template<class JTrackSelector_t >
bool	has_reconstructed_track (const Evt &evt, JTrackSelector_t selector)
	Test whether given event has a track according selection. More...

template<int reconstruction_type>
bool	has_reconstructed_track (const Evt &evt, const JRange< int > range=JRange< int >())
	Test whether given event has a track according selection. More...

bool	has_reconstructed_muon (const Evt &evt)
	Test whether given event has a track with muon reconstruction. More...

bool	has_reconstructed_shower (const Evt &evt)
	Test whether given event has a track with shower reconstruction. More...

bool	has_reconstructed_aashower (const Evt &evt)
	Test whether given event has a track with aashower reconstruction. More...

template<class JTrackSelector_t , class JQualitySorter_t >
const Trk &	get_best_reconstructed_track (const Evt &evt, JTrackSelector_t selector, JQualitySorter_t comparator)
	Get best reconstructed track. More...

template<int reconstruction_type>
const Trk &	get_best_reconstructed_track (const Evt &evt, const JRange< int > range=JRange< int >())
	Get best reconstructed track. More...

const Trk &	get_best_reconstructed_muon (const Evt &evt)
	Get best reconstructed muon. More...

const Trk &	get_best_reconstructed_shower (const Evt &evt)
	Get best reconstructed shower. More...

const Trk &	get_best_reconstructed_aashower (const Evt &evt)
	Get best reconstructed aashower. More...

void	copy (const Head &from, JHead &to)
	Copy header from `from` to `to`. More...

void	copy (const JHead &from, Head &to)
	Copy header from `from` to `to`. More...

std::string	getTag (const std::string &tag)
	Get tag without aanet extension "_<counter>" for identical tags. More...

template<class T >
bool	is_valid (const T &value)
	Check validity of given value. More...

bool	operator== (const Head &first, const Head &second)
	Equal operator. More...

bool	operator< (const Head &first, const Head &second)
	Less than operator. More...

bool	is_mupage (const JHead &header)
	Check for generator. More...

bool	is_genhen (const JHead &header)
	Check for generator. More...

bool	is_genie (const JHead &header)
	Check for generator. More...

bool	is_sirene (const JHead &header)
	Check for detector simulation. More...

bool	is_km3 (const JHead &header)
	Check for detector simulation. More...

bool	is_km3sim (const JHead &header)
	Check for detector simulation. More...

bool	is_daq (const JHead &header)
	Check for real data. More...

template<class T >
T	get (const JHead &header)
	Get object from header. More...

template<>
Vec	get (const JHead &header)
	Get position offset of detector due to generator. More...

template<>
JPosition3D	get (const JHead &header)
	Get position offset of detector due to generator. More...

template<>
JCylinder3D	get (const JHead &header)
	Get cylinder corresponding to can. More...

Variables
static const int	AASHOWER_RECONSTRUCTION_TYPE = 101
	AAshower reconstruction type for AAnet. More...

static const JEvtEvaluator	getEvtValue
	Function object for evaluation of DAQ objects. More...

static const char	AANET_TAG_SEPARATOR = '_'
	Separator for AAnet tag extension for multiple tags ("_<counter>"). More...

Detailed Description

Extensions to AAnet data format.

Author: mlincetto; mdejong

Enumeration Type Documentation

◆ JHitType_t

enum JAANET::JHitType_t

Enumeration of hit types based on km3 codes.

Enumerator
HIT_TYPE_MUON_DIRECT	Direct light from muon.
HIT_TYPE_MUON_SCATTERED	Scattered light from muon.
HIT_TYPE_DELTARAYS_DIRECT	Direct light from delta-rays.
HIT_TYPE_DELTARAYS_SCATTERED	Scattered light from delta-rays.
HIT_TYPE_BREMSSTRAHLUNG_DIRECT	Direct light from Bremsstrahlung.
HIT_TYPE_BREMSSTRAHLUNG_SCATTERED	Scattered light from Bremsstrahlung.
HIT_TYPE_SHOWER_DIRECT	Direct light from primary shower.
HIT_TYPE_SHOWER_SCATTERED	Scattered light from primary shower.
HIT_TYPE_NOISE	Random noise.
HIT_TYPE_UNKNOWN	Unknown source.

Definition at line 67 of file JAAnetToolkit.hh.

                   {
     HIT_TYPE_MUON_DIRECT              =   +5,    //!< Direct    light from muon
     HIT_TYPE_MUON_SCATTERED           =   -5,    //!< Scattered light from muon
     HIT_TYPE_DELTARAYS_DIRECT         =   +4,    //!< Direct    light from delta-rays
     HIT_TYPE_DELTARAYS_SCATTERED      =   -4,    //!< Scattered light from delta-rays
     HIT_TYPE_BREMSSTRAHLUNG_DIRECT    =   +3,    //!< Direct    light from Bremsstrahlung
     HIT_TYPE_BREMSSTRAHLUNG_SCATTERED =   -3,    //!< Scattered light from Bremsstrahlung
     HIT_TYPE_SHOWER_DIRECT            =  +99,    //!< Direct    light from primary shower
     HIT_TYPE_SHOWER_SCATTERED         =  -99,    //!< Scattered light from primary shower
     HIT_TYPE_NOISE                    =   -1,    //!< Random noise
     HIT_TYPE_UNKNOWN                  =    0     //!< Unknown source
   };

◆ JGeant4Type_t

enum JAANET::JGeant4Type_t

Enumeration of hit types based on Geant4 codes, for compatbility with KM3Sim.

Enumerator
GEANT4_TYPE_PHOTON
GEANT4_TYPE_ANTIELECTRON
GEANT4_TYPE_ELECTRON
GEANT4_TYPE_NEUTRINO
GEANT4_TYPE_ANTIMUON
GEANT4_TYPE_MUON
GEANT4_TYPE_NEUTRAL_PION
GEANT4_TYPE_PION_PLUS
GEANT4_TYPE_PION_MINUS
GEANT4_TYPE_KAON_LONG
GEANT4_TYPE_KAON_PLUS
GEANT4_TYPE_KAON_MINUS
GEANT4_TYPE_NEUTRON
GEANT4_TYPE_PROTON
GEANT4_TYPE_ANTIPROTON
GEANT4_TYPE_KAON_SHORT
GEANT4_TYPE_ETA
GEANT4_TYPE_LAMBDA
GEANT4_TYPE_SIGMA_PLUS
GEANT4_TYPE_NEUTRAL_SIGMA
GEANT4_TYPE_SIGMA_MINUS
GEANT4_TYPE_NEUTRAL_XI
GEANT4_TYPE_XI_MINUS
GEANT4_TYPE_OMEGA_MINUS
GEANT4_TYPE_ANTINEUTRON
GEANT4_TYPE_ANTILAMBDA
GEANT4_TYPE_ANTISIGMA_MINUS
GEANT4_TYPE_NEUTRAL_ANTISIGMA
GEANT4_TYPE_ANTISIGMA_PLUS
GEANT4_TYPE_NEUTRAL_ANTIXI
GEANT4_TYPE_ANTIXI_PLUS
GEANT4_TYPE_ANTIOMEGA_PLUS
GEANT4_TYPE_DEUTERON
GEANT4_TYPE_TRITON
GEANT4_TYPE_ALPHA
GEANT4_TYPE_GEANTINO
GEANT4_TYPE_HE3
GEANT4_TYPE_ANTITAU
GEANT4_TYPE_TAU

Definition at line 18 of file JParticleTypes.hh.

                      { GEANT4_TYPE_PHOTON            =   1,
                        GEANT4_TYPE_ANTIELECTRON      =   2,
                        GEANT4_TYPE_ELECTRON          =   3,
                        GEANT4_TYPE_NEUTRINO          =   4,
                        GEANT4_TYPE_ANTIMUON          =   5,
                        GEANT4_TYPE_MUON              =   6,
                        GEANT4_TYPE_NEUTRAL_PION      =   7,
                        GEANT4_TYPE_PION_PLUS         =   8,
                        GEANT4_TYPE_PION_MINUS        =   9,
                        GEANT4_TYPE_KAON_LONG         =  10,
                        GEANT4_TYPE_KAON_PLUS         =  11,
                        GEANT4_TYPE_KAON_MINUS        =  12,
                        GEANT4_TYPE_NEUTRON           =  13,
                        GEANT4_TYPE_PROTON            =  14,
                        GEANT4_TYPE_ANTIPROTON        =  15,
                        GEANT4_TYPE_KAON_SHORT        =  16,
                        GEANT4_TYPE_ETA               =  17,
                        GEANT4_TYPE_LAMBDA            =  18,
                        GEANT4_TYPE_SIGMA_PLUS        =  19,
                        GEANT4_TYPE_NEUTRAL_SIGMA     =  20,
                        GEANT4_TYPE_SIGMA_MINUS       =  21,
                        GEANT4_TYPE_NEUTRAL_XI        =  22,
                        GEANT4_TYPE_XI_MINUS          =  23,
                        GEANT4_TYPE_OMEGA_MINUS       =  24,
                        GEANT4_TYPE_ANTINEUTRON       =  25,
                        GEANT4_TYPE_ANTILAMBDA        =  26,
                        GEANT4_TYPE_ANTISIGMA_MINUS   =  27,
                        GEANT4_TYPE_NEUTRAL_ANTISIGMA =  28,
                        GEANT4_TYPE_ANTISIGMA_PLUS    =  29,
                        GEANT4_TYPE_NEUTRAL_ANTIXI    =  30,
                        GEANT4_TYPE_ANTIXI_PLUS       =  31,
                        GEANT4_TYPE_ANTIOMEGA_PLUS    =  32,
                        GEANT4_TYPE_DEUTERON          =  45,
                        GEANT4_TYPE_TRITON            =  46,
                        GEANT4_TYPE_ALPHA             =  47,
                        GEANT4_TYPE_GEANTINO          =  48,
                        GEANT4_TYPE_HE3               =  49,                   
                        //
                        // KM3NeT specific codes
                        //
                        GEANT4_TYPE_ANTITAU           =  33,
                        GEANT4_TYPE_TAU               =  34 };

◆ JTrackType_t

enum JAANET::JTrackType_t

Enumeration of track types based on PDG codes.

Enumerator
TRACK_TYPE_ELECTRON
TRACK_TYPE_NUE
TRACK_TYPE_MUON
TRACK_TYPE_NUMU
TRACK_TYPE_TAU
TRACK_TYPE_NUTAU
TRACK_TYPE_PHOTON
TRACK_TYPE_NEUTRAL_PION
TRACK_TYPE_CHARGED_PION_PLUS
TRACK_TYPE_PION_PLUS
TRACK_TYPE_K_LONG
TRACK_TYPE_K_SHORT
TRACK_TYPE_K_PLUS
TRACK_TYPE_PROTON
TRACK_TYPE_NEUTRON
TRACK_TYPE_LAMBDA
TRACK_TYPE_SIGMA_PLUS
TRACK_TYPE_NEUTRAL_SIGMA
TRACK_TYPE_SIGMA_MINUS
TRACK_TYPE_NEUTRAL_XI
TRACK_TYPE_XI_MINUS
TRACK_TYPE_OMEGA_MINUS
TRACK_TYPE_ANTIELECTRON
TRACK_TYPE_ANTINUE
TRACK_TYPE_ANTIMUON
TRACK_TYPE_ANTINUMU
TRACK_TYPE_ANTITAU
TRACK_TYPE_ANTINUTAU
TRACK_TYPE_NEUTRAL_ANTIPION
TRACK_TYPE_CHARGED_PION_MINUS
TRACK_TYPE_PION_MINUS
TRACK_TYPE_ANTIK_LONG
TRACK_TYPE_ANTIK_SHORT
TRACK_TYPE_K_MINUS
TRACK_TYPE_ANTIPROTON
TRACK_TYPE_ANTINEUTRON
TRACK_TYPE_ANTILAMBDA
TRACK_TYPE_ANTISIGMA_PLUS
TRACK_TYPE_ANTINEUTRAL_SIGMA
TRACK_TYPE_ANTISIGMA_MINUS
TRACK_TYPE_ANTINEUTRAL_XI
TRACK_TYPE_ANTIXI_MINUS
TRACK_TYPE_ANTIOMEGA_MINUS

Definition at line 64 of file JParticleTypes.hh.

                     { TRACK_TYPE_ELECTRON           =  11,
                       TRACK_TYPE_NUE                =  12,
                       TRACK_TYPE_MUON               =  13,
                       TRACK_TYPE_NUMU               =  14,
                       TRACK_TYPE_TAU                =  15,
                       TRACK_TYPE_NUTAU              =  16,
                       TRACK_TYPE_PHOTON             =  22,
                       TRACK_TYPE_NEUTRAL_PION       =  111,
                       TRACK_TYPE_CHARGED_PION_PLUS  =  211,
                       TRACK_TYPE_PION_PLUS          =  211,
                       TRACK_TYPE_K_LONG             =  130,
                       TRACK_TYPE_K_SHORT            =  310,
                       TRACK_TYPE_K_PLUS             =  321,
                       TRACK_TYPE_PROTON             =  2212,
                       TRACK_TYPE_NEUTRON            =  2112,
                       TRACK_TYPE_LAMBDA             =  3122,
                       TRACK_TYPE_SIGMA_PLUS         =  3222,
                       TRACK_TYPE_NEUTRAL_SIGMA      =  3212,
                       TRACK_TYPE_SIGMA_MINUS        =  3112,
                       TRACK_TYPE_NEUTRAL_XI         =  3322,
                       TRACK_TYPE_XI_MINUS           =  3312,
                       TRACK_TYPE_OMEGA_MINUS        =  3334,
  
                       TRACK_TYPE_ANTIELECTRON       = -11,
                       TRACK_TYPE_ANTINUE            = -12,
                       TRACK_TYPE_ANTIMUON           = -13,
                       TRACK_TYPE_ANTINUMU           = -14,
                       TRACK_TYPE_ANTITAU            = -15,
                       TRACK_TYPE_ANTINUTAU          = -16,
                       TRACK_TYPE_NEUTRAL_ANTIPION   = -111,
                       TRACK_TYPE_CHARGED_PION_MINUS = -211,
                       TRACK_TYPE_PION_MINUS         = -211,
                       TRACK_TYPE_ANTIK_LONG         = -130,
                       TRACK_TYPE_ANTIK_SHORT        = -310,
                       TRACK_TYPE_K_MINUS            = -321,
                       TRACK_TYPE_ANTIPROTON         = -2212,
                       TRACK_TYPE_ANTINEUTRON        = -2112,
                       TRACK_TYPE_ANTILAMBDA         = -3122,
                       TRACK_TYPE_ANTISIGMA_PLUS     = -3222,
                       TRACK_TYPE_ANTINEUTRAL_SIGMA  = -3212,
                       TRACK_TYPE_ANTISIGMA_MINUS    = -3112,
                       TRACK_TYPE_ANTINEUTRAL_XI     = -3322,
                       TRACK_TYPE_ANTIXI_MINUS       = -3312,
                       TRACK_TYPE_ANTIOMEGA_MINUS    = -3334 };

Function Documentation

◆ getTime()

double JAANET::getTime ( const Hit & hit )

inline

Get true time of hit.

Parameters

hit hit

Returns: true time of hit

Definition at line 87 of file JAAnetToolkit.hh.

   {
     // if IGAIN = 0 in km3 program, then pure_xx values are not set.
     
     if (hit.pure_a == 0.0)
       return hit.t;
     else
       return hit.pure_t;
   }

◆ getNPE()

double JAANET::getNPE ( const Hit & hit )

inline

Get true charge of hit.

Parameters

hit hit

Returns: true charge of hit

Definition at line 104 of file JAAnetToolkit.hh.

   {
     // if IGAIN = 0 in km3 program, then pure_xx values are not set.
     
     if (hit.pure_a == 0.0)
       return hit.a;
     else
       return hit.pure_a;
   }

◆ is_noise()

bool JAANET::is_noise ( const Hit & hit )

inline

Verify hit origin.

Parameters

hit hit

Returns: true if noise; else false

Definition at line 121 of file JAAnetToolkit.hh.

   {
     return hit.type == HIT_TYPE_NOISE;
   }

◆ getTimeRange() [1/2]

JTimeRange JAANET::getTimeRange ( const Evt & event )

inline

Get time range (i.e.

earlist and latest hit time) of Monte Carlo event. Note that the global event time in not taken into account.

Parameters

event event

Returns: time range

Definition at line 134 of file JAAnetToolkit.hh.

   {
     JTimeRange timeRange(JTimeRange::DEFAULT_RANGE);
  
     for (std::vector<Hit>::const_iterator hit = event.mc_hits.begin(); hit != event.mc_hits.end(); ++hit) {
       if (!is_noise(*hit)) {
         timeRange.include(getTime(*hit));
       }
     }
  
     return timeRange;
   }

◆ getTimeRange() [2/2]

JTimeRange JAANET::getTimeRange	(	const Evt &	event,
		const JTimeRange &	T_ns
	)

inline

Get time range (i.e.

earlist and latest hit time) of Monte Carlo event. The given time window is offset so that the lower limit is equal to the earliest hit time. The time window is then applied to all hits. Note that the global event time in not taken into account.

Parameters

event	event
T_ns	time window [ns]

Returns: time range

Definition at line 158 of file JAAnetToolkit.hh.

   {
     using namespace std;
  
     JTimeRange timeRange(JTimeRange::DEFAULT_RANGE);
  
     double t0 = numeric_limits<double>::max();
  
     for (std::vector<Hit>::const_iterator hit = event.mc_hits.begin(); hit != event.mc_hits.end(); ++hit) {
       if (!is_noise(*hit)) {
         if (getTime(*hit) < t0) {
           t0 = getTime(*hit);
         }
       }
     }
  
     if (t0 != numeric_limits<double>::max()) {
  
       JTimeRange window(T_ns);
     
       window.fixLowerLimit(t0);
     
       for (std::vector<Hit>::const_iterator hit = event.mc_hits.begin(); hit != event.mc_hits.end(); ++hit) {
         if (!is_noise(*hit) && window(getTime(*hit))) {
           timeRange.include(getTime(*hit));
         }
       }
     }
  
     return timeRange;
   }

◆ getPosition() [1/2]

JPosition3D JAANET::getPosition ( const Vec & v )

inline

Get position.

Parameters

v vector

Returns: position

Definition at line 197 of file JAAnetToolkit.hh.

   {
     return JPosition3D(v.x, v.y, v.z);
   }

◆ getPosition() [2/2]

JPosition3D JAANET::getPosition ( const Trk & track )

inline

Get position.

Parameters

track track

Returns: position

Definition at line 209 of file JAAnetToolkit.hh.

   {
     return getPosition(track.pos);
   }

◆ getDirection() [1/2]

JDirection3D JAANET::getDirection ( const Vec & v )

inline

Get direction.

Parameters

v vector

Returns: direction

Definition at line 221 of file JAAnetToolkit.hh.

   {
     return JDirection3D(v.x, v.y, v.z);
   }

◆ getDirection() [2/2]

JDirection3D JAANET::getDirection ( const Trk & track )

inline

Get direction.

Parameters

track track

Returns: direction

Definition at line 232 of file JAAnetToolkit.hh.

   {
     return getDirection(track.dir);
   }

◆ getAxis()

JAxis3D JAANET::getAxis ( const Trk & track )

inline

Get axis.

Parameters

track track

Returns: axis

Definition at line 244 of file JAAnetToolkit.hh.

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

◆ getTrack()

JTrack3E JAANET::getTrack ( const Trk & track )

inline

Get track.

Parameters

track track

Returns: track

Definition at line 256 of file JAAnetToolkit.hh.

   {
     return JTrack3E(JTrack3D(getAxis(track), track.t), track.E);
   }

◆ getVertex()

JVertex3D JAANET::getVertex ( const Trk & track )

inline

Get vertex.

Parameters

track track

Returns: vertex

Definition at line 268 of file JAAnetToolkit.hh.

   {
     return JVertex3D(getPosition(track), track.t);
   }

◆ getTransformation()

JTransformation3D JAANET::getTransformation ( const Trk & track )

inline

Get transformation.

Parameters

track track

Returns: transformation

Definition at line 280 of file JAAnetToolkit.hh.

   {
     return JTransformation3D(getPosition(track), getDirection(track));
   }

◆ getHit() [1/3]

JHitL0 JAANET::getHit ( const Hit & hit )

inline

Get transformation.

Parameters

hit hit

Returns: hit

Definition at line 292 of file JAAnetToolkit.hh.

   {
     using namespace JPP;
  
     return JHitL0(JDAQPMTIdentifier(hit.dom_id, hit.channel_id), JAxis3D(getPosition(hit.pos), getDirection(hit.dir)), JHit(hit.t, hit.tot) );
   }

◆ getHit() [2/3]

JHitL0 JAANET::getHit	(	const Hit &	hit,
		const JModuleRouter &	router
	)

inline

Get transformation.

Parameters

hit	hit
router	module router

Returns: hit

Definition at line 307 of file JAAnetToolkit.hh.

   {
     using namespace JPP;
  
     const JPMT& pmt = router.getModule(hit.dom_id).getPMT(hit.channel_id);
  
     return JHitL0(JDAQPMTIdentifier(hit.dom_id, hit.channel_id), pmt.getAxis(), JHit(getTime(hit.tdc, pmt.getCalibration()), getToT(hit.tot, pmt.getCalibration())));
   }

◆ getHit() [3/3]

JHitL0 JAANET::getHit	(	const Hit &	hit,
		const JPMTRouter &	router
	)

inline

Get transformation.

Parameters

hit	hit
router	pmt router

Returns: hit

Definition at line 324 of file JAAnetToolkit.hh.

   {
     using namespace JPP;
  
     const JPMT& pmt = router.getPMT(hit.pmt_id);
  
     return JHitL0(JDAQPMTIdentifier(hit.dom_id, hit.channel_id), pmt.getAxis(), JHit(getTime(hit.tdc, pmt.getCalibration()), getToT(hit.tot, pmt.getCalibration())));
   }

◆ is_photon()

bool JAANET::is_photon ( const Trk & track )

inline

Test whether given track is a photon or neutral pion.

Parameters

track track

Returns: true if photon or neutral pion; else false

Definition at line 340 of file JAAnetToolkit.hh.

340 { return ( track.type == TRACK_TYPE_PHOTON ||

341 abs(track.type) == TRACK_TYPE_NEUTRAL_PION); }

◆ is_neutrino()

bool JAANET::is_neutrino ( const Trk & track )

inline

Test whether given track is a neutrino.

Parameters

track track

Returns: true if neutrino; else false

Definition at line 349 of file JAAnetToolkit.hh.

                                             { return (abs(track.type) == TRACK_TYPE_NUE  ||
                                                       abs(track.type) == TRACK_TYPE_NUMU ||
                                                       abs(track.type) == TRACK_TYPE_NUTAU); }

◆ is_electron()

bool JAANET::is_electron ( const Trk & track )

inline

Test whether given track is a (anti-)electron.

Parameters

track track

Returns: true if (anti-)electron; else false

Definition at line 359 of file JAAnetToolkit.hh.

359 { return (abs(track.type) == TRACK_TYPE_ELECTRON); }

◆ is_muon()

bool JAANET::is_muon ( const Trk & track )

inline

Test whether given track is a (anti-)muon.

Parameters

track track

Returns: true if (anti-)muon; else false

Definition at line 367 of file JAAnetToolkit.hh.

367 { return (abs(track.type) == TRACK_TYPE_MUON); }

◆ is_tau()

bool JAANET::is_tau ( const Trk & track )

inline

Test whether given track is a (anti-)tau.

Parameters

track track

Returns: true if (anti-)tau; else false

Definition at line 375 of file JAAnetToolkit.hh.

375 { return (abs(track.type) == TRACK_TYPE_TAU); }

◆ is_pion()

bool JAANET::is_pion ( const Trk & track )

inline

Test whether given track is a charged pion.

Parameters

track track

Returns: true if charged pion; else false

Definition at line 383 of file JAAnetToolkit.hh.

383 { return (abs(track.type) == TRACK_TYPE_CHARGED_PION_PLUS); }

◆ is_proton()

bool JAANET::is_proton ( const Trk & track )

inline

Test whether given track is a (anti-)proton.

Parameters

track track

Returns: true if (anti-)proton; else false

Definition at line 391 of file JAAnetToolkit.hh.

391 { return (abs(track.type) == TRACK_TYPE_PROTON); }

◆ is_neutron()

bool JAANET::is_neutron ( const Trk & track )

inline

Test whether given track is a (anti-)neutron.

Parameters

track track

Returns: true if (anti-)neutron; else false

Definition at line 399 of file JAAnetToolkit.hh.

399 { return (abs(track.type) == TRACK_TYPE_NEUTRON); }

◆ is_hadron()

bool JAANET::is_hadron ( const Trk & track )

inline

Test whether given track is a hadron.

Parameters

track track

Returns: true if hadron; else fails

Definition at line 407 of file JAAnetToolkit.hh.

                                             { return (!is_neutrino(track) &&
                                                       !is_electron(track) &&
                                                       !is_muon    (track) &&
                                                       !is_tau     (track)); }

◆ has_particleID()

bool JAANET::has_particleID	(	const Trk &	track,
		const int	type
	)

inline

Test whether given track corresponds to given particle type.

Parameters

track	track
type	particle type

Returns: true if matches the corresponding particle; else false

Definition at line 419 of file JAAnetToolkit.hh.

419 { return (track.type == type); }

◆ has_neutrino()

bool JAANET::has_neutrino ( const Evt & evt )

inline

Test whether given event has an incoming neutrino.

Parameters

evt event

Returns: true if neutrino; else false

Definition at line 427 of file JAAnetToolkit.hh.

   { 
     return !evt.mc_trks.empty() && is_neutrino(evt.mc_trks[0]);
   }

◆ get_neutrino()

const Trk& JAANET::get_neutrino ( const Evt & evt )

inline

Get incoming neutrino.

Parameters

evt event

Returns: neutrino

Definition at line 438 of file JAAnetToolkit.hh.

   {
     if (has_neutrino(evt))
       return evt.mc_trks[0];
     else
       THROW(JIndexOutOfRange, "This event has no neutrino.");
   }

◆ count_electrons()

int JAANET::count_electrons ( const Evt & evt )

inline

Count the number of electrons in a given event.

Parameters

evt event

Returns: number of electrons

Definition at line 452 of file JAAnetToolkit.hh.

   {
     return std::count_if(evt.mc_trks.begin(), evt.mc_trks.end(), is_electron);
   }

◆ has_electron()

bool JAANET::has_electron ( const Evt & evt )

inline

Test whether given event has an electron.

Parameters

evt event

Returns: true if event has electron; else false

Definition at line 463 of file JAAnetToolkit.hh.

   {
     return count_electrons(evt) != 0;
   }

◆ get_electron()

const Trk& JAANET::get_electron ( const Evt & evt )

inline

Get first electron from the event tracklist.

Parameters

evt event

Returns: electron

Definition at line 474 of file JAAnetToolkit.hh.

   {
     if (count_electrons(evt) > 0)
       return *std::find_if(evt.mc_trks.begin(), evt.mc_trks.end(), is_electron);
     else
       THROW(JIndexOutOfRange, "This event has no electron.");
   }

◆ from_electron()

bool JAANET::from_electron ( const Hit & hit )

inline

Test whether given hit was produced by an electron.

Warning: This method will only work with the output of KM3Sim. For JSirene or KM3, you should check if track.id == hit.origin instead.

Parameters

hit hit

Returns: true if hit was produced by an electron; else false

Definition at line 491 of file JAAnetToolkit.hh.

   {
     return hit.type == GEANT4_TYPE_ELECTRON || hit.type == GEANT4_TYPE_ANTIELECTRON;
   }

◆ count_muons()

int JAANET::count_muons ( const Evt & evt )

inline

Count the number of muons in a given event.

Parameters

evt event

Returns: number of muons

Definition at line 502 of file JAAnetToolkit.hh.

   {
     return std::count_if(evt.mc_trks.begin(), evt.mc_trks.end(), is_muon);
   }

◆ has_muon()

bool JAANET::has_muon ( const Evt & evt )

inline

Test whether given event has a muon.

Parameters

evt event

Returns: true if event has muon; else false

Definition at line 513 of file JAAnetToolkit.hh.

   {
     return count_muons(evt) != 0;
   }

◆ get_muon()

const Trk& JAANET::get_muon ( const Evt & evt )

inline

Get first muon from the event tracklist.

Parameters

evt event

Returns: muon

Definition at line 524 of file JAAnetToolkit.hh.

   {
     if (count_muons(evt) > 0)
       return *std::find_if(evt.mc_trks.begin(), evt.mc_trks.end(), is_muon);
     else
       THROW(JIndexOutOfRange, "This event has no muon.");
   }

◆ from_muon()

bool JAANET::from_muon ( const Hit & hit )

inline

Test whether given hit was produced by a muon.

Warning: This method will only work with the output of KM3Sim. For JSirene or KM3, you should check if track.id == hit.origin instead.

Parameters

hit hit

Returns: true if hit was produced by a muon; else false

Definition at line 541 of file JAAnetToolkit.hh.

   {
     return hit.type == GEANT4_TYPE_MUON || hit.type == GEANT4_TYPE_ANTIMUON;
   }

◆ count_taus()

int JAANET::count_taus ( const Evt & evt )

inline

Count the number of taus in a given event.

Parameters

evt event

Returns: number of taus

Definition at line 552 of file JAAnetToolkit.hh.

   {
     return std::count_if(evt.mc_trks.begin(), evt.mc_trks.end(), is_tau);
   }

◆ has_tau()

bool JAANET::has_tau ( const Evt & evt )

inline

Test whether given event has a tau.

Parameters

evt event

Returns: true if event has tau; else false

Definition at line 563 of file JAAnetToolkit.hh.

   {
     return count_taus(evt) != 0;
   }

◆ get_tau()

const Trk& JAANET::get_tau ( const Evt & evt )

inline

Get first tau from the event tracklist.

Parameters

evt event

Returns: tau

Definition at line 574 of file JAAnetToolkit.hh.

   {
     if (count_taus(evt) > 0)
       return *std::find_if(evt.mc_trks.begin(), evt.mc_trks.end(), is_tau);
     else
       THROW(JIndexOutOfRange, "This event has no tau.");
   }

◆ from_tau()

bool JAANET::from_tau ( const Hit & hit )

inline

Test whether given hit was produced by a tau.

Warning: This method will only work with the output of KM3Sim. For JSirene or KM3, you should check if track.id == hit.origin instead.

Parameters

hit hit

Returns: true if hit was produced by a tau; else false

Definition at line 591 of file JAAnetToolkit.hh.

   {
     return hit.type == GEANT4_TYPE_TAU || hit.type == GEANT4_TYPE_ANTITAU;
   }

◆ count_hadrons()

int JAANET::count_hadrons ( const Evt & evt )

inline

Count the number of hadrons in a given event (not including neutral pions)

Parameters

evt event

Returns: number of hadrons

Definition at line 602 of file JAAnetToolkit.hh.

   {
     return std::count_if(evt.mc_trks.begin(), evt.mc_trks.end(), is_hadron);
   }

◆ from_hadron()

bool JAANET::from_hadron ( const Hit & hit )

inline

Test whether given hit was produced by a hadronic shower.

Warning: This method will only work with the output of KM3Sim. For JSirene or KM3, you should check if track.id == hit.origin instead.

Parameters

hit hit

Returns: true if hit was produced by a hadron; else false

Definition at line 616 of file JAAnetToolkit.hh.

   {
     return (!from_electron(hit) && !from_muon(hit) && !from_tau(hit));
   }

◆ get_hadronic_cascade()

Trk JAANET::get_hadronic_cascade ( const Evt & evt )

inline

Auxiliary function to get average direction and total energy of a hadronic cascade.

Parameters

evt event

Returns: cascade (default if there are no hadrons)

Definition at line 627 of file JAAnetToolkit.hh.

   {
     Trk cascade;
     
     for (vector<Trk>::const_iterator track = evt.mc_trks.begin(); track != evt.mc_trks.end(); ++track) {
       if (is_hadron(*track)) {
         cascade.E   += track->E;
         cascade.dir += track->dir * track->E;
         cascade.t    = track->t;
         cascade.pos  = track->pos;
       }
     }
  
     cascade.dir.normalize();
  
     return cascade;
   }

◆ add_time_offset()

void JAANET::add_time_offset	(	Evt &	evt,
		const double	tOff
	)

inline

Add time offset to mc event; according to current definition, the absolute time of the event is defined by the track "t" attribute; this could change in the future if the global attribute mc_t is assigned to this purpose.

Parameters

evt	event
tOff	time offset [ns]

Definition at line 653 of file JAAnetToolkit.hh.

   {
     for (vector<Trk>::iterator p = evt.mc_trks.begin(); p != evt.mc_trks.end(); p++) {
       p->t += tOff;
     }
   }

◆ has_muon_prefit()

bool JAANET::has_muon_prefit ( const Trk & track )

inline

Test whether given track has muon prefit in history.

Parameters

track track

Returns: true if muon prefit in history; else false

Definition at line 727 of file JAAnetToolkit.hh.

   {
     return has_history(JPP_RECONSTRUCTION_TYPE, JFIT::JMUONPREFIT)(track);
   }

◆ has_muon_simplex()

bool JAANET::has_muon_simplex ( const Trk & track )

inline

Test whether given track has muon simplex fit in history.

Parameters

track track

Returns: true if muon simplex fit in history; else false

Definition at line 738 of file JAAnetToolkit.hh.

   {
     return has_history(JPP_RECONSTRUCTION_TYPE, JFIT::JMUONSIMPLEX)(track);
   }

◆ has_muon_gandalf()

bool JAANET::has_muon_gandalf ( const Trk & track )

inline

Test whether given track has muon gandalf fit in history.

Parameters

track track

Returns: true if muon gandalf fit in history; else false

Definition at line 749 of file JAAnetToolkit.hh.

   {
     return has_history(JPP_RECONSTRUCTION_TYPE, JFIT::JMUONGANDALF)(track);
   }

◆ has_muon_energy()

bool JAANET::has_muon_energy ( const Trk & track )

inline

Test whether given track has muon energy fit in history.

Parameters

track track

Returns: true if muon energy fit in history; else false

Definition at line 760 of file JAAnetToolkit.hh.

   {
     return has_history(JPP_RECONSTRUCTION_TYPE, JFIT::JMUONENERGY)(track);
   }

◆ has_muon_start()

bool JAANET::has_muon_start ( const Trk & track )

inline

Test whether given track has muon start fit in history.

Parameters

track track

Returns: true if muon start fit in history; else false

Definition at line 771 of file JAAnetToolkit.hh.

   {
     return has_history(JPP_RECONSTRUCTION_TYPE, JFIT::JMUONSTART)(track);
   }

◆ has_muon_fit()

bool JAANET::has_muon_fit ( const Trk & track )

inline

Test whether given track has default muon fit in history.

Parameters

track track

Returns: true if muon fit in history; else false

Definition at line 782 of file JAAnetToolkit.hh.

   {
     return has_history(JPP_RECONSTRUCTION_TYPE, JRange<int>(JFIT::JMUONBEGIN, JFIT::JMUONEND))(track);
   }  

◆ has_shower_prefit()

bool JAANET::has_shower_prefit ( const Trk & track )

inline

Test whether given track has shower prefit in history.

Parameters

track track

Returns: true if shower prefit in history; else false

Definition at line 793 of file JAAnetToolkit.hh.

   {
     return has_history(JPP_RECONSTRUCTION_TYPE, JFIT::JSHOWERPREFIT)(track);
   }

◆ has_shower_positionfit()

bool JAANET::has_shower_positionfit ( const Trk & track )

inline

Test whether given track has shower position fit in history.

Parameters

track track

Returns: true if shower position fit in history; else false

Definition at line 804 of file JAAnetToolkit.hh.

   {
     return has_history(JPP_RECONSTRUCTION_TYPE, JFIT::JSHOWERPOSITIONFIT)(track);
   }  

◆ has_shower_completefit()

bool JAANET::has_shower_completefit ( const Trk & track )

inline

Test whether given track has shower complete fit in history.

Parameters

track track

Returns: true if shower complete fit in history; else false

Definition at line 815 of file JAAnetToolkit.hh.

   {
     return has_history(JPP_RECONSTRUCTION_TYPE, JFIT::JSHOWERCOMPLETEFIT)(track);
   }

◆ has_shower_fit()

bool JAANET::has_shower_fit ( const Trk & track )

inline

Test whether given track has default shower fit in history.

Parameters

track track

Returns: true if shower fit in history; else false

Definition at line 826 of file JAAnetToolkit.hh.

   {
     return has_history(JPP_RECONSTRUCTION_TYPE, JRange<int>(JFIT::JSHOWERBEGIN, JFIT::JSHOWEREND))(track);
   }

◆ has_aashower_fit()

bool JAANET::has_aashower_fit ( const Trk & track )

inline

Test whether given track has default shower fit in history.

Parameters

track track

Returns: true if shower fit in history; else false

Definition at line 837 of file JAAnetToolkit.hh.

   {
     return has_history(AASHOWER_RECONSTRUCTION_TYPE, JRange<int>())(track);
   }

◆ has_reconstructed_track() [1/2]

template<class JTrackSelector_t >

bool JAANET::has_reconstructed_track	(	const Evt &	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 851 of file JAAnetToolkit.hh.

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

◆ has_reconstructed_track() [2/2]

template<int reconstruction_type>

bool JAANET::has_reconstructed_track	(	const Evt &	evt,
		const JRange< int >	range = `JRange<int>()`
	)

inline

Test whether given event has a track according selection.

Parameters

evt	event
range	range of application types

Returns: true if at least one corresponding track; else false

Definition at line 864 of file JAAnetToolkit.hh.

   {
     return has_reconstructed_track(evt, has_history(reconstruction_type, range));
   }

◆ has_reconstructed_muon()

bool JAANET::has_reconstructed_muon ( const Evt & 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 875 of file JAAnetToolkit.hh.

   {
     return has_reconstructed_track(evt, has_muon_fit);
   }

◆ has_reconstructed_shower()

bool JAANET::has_reconstructed_shower ( const Evt & 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 886 of file JAAnetToolkit.hh.

   {
     return has_reconstructed_track(evt, has_shower_fit);
   }

◆ has_reconstructed_aashower()

bool JAANET::has_reconstructed_aashower ( const Evt & evt )

inline

Test whether given event has a track with aashower reconstruction.

Parameters

evt event

Returns: true if at least one reconstructed shower; else false

Definition at line 897 of file JAAnetToolkit.hh.

   {
     return has_reconstructed_track(evt, has_aashower_fit);
   }

◆ get_best_reconstructed_track() [1/2]

template<class JTrackSelector_t , class JQualitySorter_t >

const Trk& JAANET::get_best_reconstructed_track	(	const Evt &	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 913 of file JAAnetToolkit.hh.

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

◆ get_best_reconstructed_track() [2/2]

template<int reconstruction_type>

const Trk& JAANET::get_best_reconstructed_track	(	const Evt &	evt,
		const JRange< int >	range = `JRange<int>()`
	)

inline

Get best reconstructed track.

Parameters

evt	event
range	range of application types

Returns: track

Definition at line 939 of file JAAnetToolkit.hh.

   {
     return get_best_reconstructed_track(evt, has_history(reconstruction_type, range), quality_sorter<reconstruction_type>());
   }

◆ get_best_reconstructed_muon()

const Trk& JAANET::get_best_reconstructed_muon ( const Evt & evt )

inline

Get best reconstructed muon.

Parameters

evt event

Returns: track

Definition at line 951 of file JAAnetToolkit.hh.

   {
     return get_best_reconstructed_track(evt, has_muon_fit, quality_sorter<JPP_RECONSTRUCTION_TYPE>());
   }

◆ get_best_reconstructed_shower()

const Trk& JAANET::get_best_reconstructed_shower ( const Evt & evt )

inline

Get best reconstructed shower.

Parameters

evt event

Returns: track

Definition at line 962 of file JAAnetToolkit.hh.

   {
     return get_best_reconstructed_track(evt, has_shower_fit, quality_sorter<JPP_RECONSTRUCTION_TYPE>());
   }

◆ get_best_reconstructed_aashower()

const Trk& JAANET::get_best_reconstructed_aashower ( const Evt & evt )

inline

Get best reconstructed aashower.

Parameters

evt event

Returns: track

Definition at line 973 of file JAAnetToolkit.hh.

   {
     return get_best_reconstructed_track(evt, has_aashower_fit, quality_sorter<AASHOWER_RECONSTRUCTION_TYPE>());
   }

◆ copy() [1/2]

void JAANET::copy	(	const Head &	from,
		JHead &	to
	)

Copy header from from to to.

Parameters

from	header
to	header

Definition at line 152 of file JHead.cc.

   {
     using namespace std;
     using namespace JPP;
  
     JRootReader        reader(null, JHead::getEquationParameters(), JAAnetDictionary::getInstance());
     JRootReadableClass cls(to);
  
     for (Head::const_iterator i = from.begin(); i != from.end(); ++i) {
  
       const JRootReadableClass& abc = cls.find(getTag(i->first).c_str());
  
       if (abc.is_valid()) {
  
         JRedirectString redirect(reader, i->second);
  
         reader.getObject(abc);
       }
  
       if (i->first == getTag(i->first)) {
         to.insert(*i);                                     // keep track of parsed tags
       }
     }
   }

◆ copy() [2/2]

void JAANET::copy	(	const JHead &	from,
		Head &	to
	)

Copy header from from to to.

Parameters

from	header
to	header

Definition at line 184 of file JHead.cc.

   {
     using namespace std;
  
     to = Head();
  
     stringstream io;
  
     io << from;
  
     read(to, io);
     
     // copy pending data
  
     for (JHead::const_iterator i = from.begin(); i != from.end(); ++i) {
       if (to.find(i->first) == to.end()) {
         to.insert(*i);
       }
     }
   }

◆ getTag()

std::string JAANET::getTag ( const std::string & tag )

inline

Get tag without aanet extension "_<counter>" for identical tags.

Parameters

tag tag

Returns: tag

Definition at line 61 of file JHead.hh.

   {
     using namespace std;
     
     const string::size_type pos = tag.find(AANET_TAG_SEPARATOR);
     
     if (pos != string::npos) {
       
       for (string::size_type i = pos + 1; i != tag.size(); ++i) {
         if (!isdigit(tag[i])) {
           return tag;
         }
       }
       
       return tag.substr(0, pos);
     }
     
     return tag;
   }

◆ is_valid()

template<class T >

bool JAANET::is_valid ( const T & value )

inline

Check validity of given value.

Parameters

value value

Returns: true if given value does not match the default value; else false

Definition at line 823 of file JHead.hh.

   {
     return !value.match(T());
   }

◆ operator==()

bool JAANET::operator==	(	const Head &	first,
		const Head &	second
	)

inline

Equal operator.

Note that this operator uses the JHead::match method.

Parameters

first	first header
second	second header

Returns: true if two headers are equal; else false

Definition at line 1231 of file JHead.hh.

   {
     return JHead(first).match(JHead(second));
   }

◆ operator<()

bool JAANET::operator<	(	const Head &	first,
		const Head &	second
	)

inline

Less than operator.

Parameters

first	first header
second	second header

Returns: true if first header is less than second header; else false

Definition at line 1245 of file JHead.hh.

   {
     return JHead(first).less(JHead(second));
   }

◆ is_mupage()

bool JAANET::is_mupage ( const JHead & header )

inline

Check for generator.

Parameters

header header

Returns: true if this Monte Carlo is produced by mupage; else false

Definition at line 199 of file JHeadToolkit.hh.

   {
     return header.match(getMUPAGEHeader(), false);
   }

◆ is_genhen()

bool JAANET::is_genhen ( const JHead & header )

inline

Check for generator.

Parameters

header header

Returns: true if this Monte Carlo is produced by genhen; else false

Definition at line 211 of file JHeadToolkit.hh.

   {
     return header.match(getGenhenHeader(), false);
   }

◆ is_genie()

bool JAANET::is_genie ( const JHead & header )

inline

Check for generator.

Parameters

header header

Returns: true if this Monte Carlo is produced by genie; else false

Definition at line 223 of file JHeadToolkit.hh.

   {
     return (header.match(getGenieHeader(),   false) ||
             header.match(getGSeaGenHeader(), false));
   }

◆ is_sirene()

bool JAANET::is_sirene ( const JHead & header )

inline

Check for detector simulation.

Parameters

header header

Returns: true if this Monte Carlo is processed with JSirene; else false

Definition at line 236 of file JHeadToolkit.hh.

   {
     return header.match(getSireneHeader(), false);
   }

◆ is_km3()

bool JAANET::is_km3 ( const JHead & header )

inline

Check for detector simulation.

Parameters

header header

Returns: true if this Monte Carlo is processed with km3; else false

Definition at line 248 of file JHeadToolkit.hh.

   {
     return header.match(getKM3Header(), false);
   }

◆ is_km3sim()

bool JAANET::is_km3sim ( const JHead & header )

inline

Check for detector simulation.

Parameters

header header

Returns: true if this Monte Carlo is processed with KM3Sim; else false

Definition at line 260 of file JHeadToolkit.hh.

   {
     return header.match(getKM3SimHeader(), false);
   }

◆ is_daq()

bool JAANET::is_daq ( const JHead & header )

inline

Check for real data.

Parameters

header header

Returns: true if this header corresponds to real data; else false

Definition at line 272 of file JHeadToolkit.hh.

   {
     return header.match(getDAQHeader(), false);
   }

◆ get() [1/4]

template<class T >

T JAANET::get ( const JHead & header )

inline

Get object from header.

Parameters

header header

Returns: object

Get object from header.

Parameters

header header

Returns: position

Get object from header.

Parameters

header header

Returns: cylinder

Definition at line 295 of file JHeadToolkit.hh.

   {
     if (is_sirene(header) || is_km3(header)) {
  
       return header.coord_origin;
  
     } else {
  
       if      (is_valid(header.can))
         return Vec(0.0, 0.0, -header.can.zmin);
       else if (is_valid(header.coord_origin))
         return header.coord_origin;
       else
         return Vec(0.0, 0.0, 0.0);
     }
   }

◆ get() [2/4]

template<>

Vec JAANET::get ( const JHead & header )

inline

Get position offset of detector due to generator.

Get object from header.

Parameters

header header

Returns: position

Definition at line 295 of file JHeadToolkit.hh.

   {
     if (is_sirene(header) || is_km3(header)) {
  
       return header.coord_origin;
  
     } else {
  
       if      (is_valid(header.can))
         return Vec(0.0, 0.0, -header.can.zmin);
       else if (is_valid(header.coord_origin))
         return header.coord_origin;
       else
         return Vec(0.0, 0.0, 0.0);
     }
   }

◆ get() [3/4]

template<>

JPosition3D JAANET::get ( const JHead & header )

inline

Get position offset of detector due to generator.

Get object from header.

Parameters

header header

Returns: position

Definition at line 320 of file JHeadToolkit.hh.

   {
     const Vec pos = get<Vec>(header);
  
     return JPosition3D(pos.x, pos.y, pos.z);
   }

◆ get() [4/4]

template<>

JCylinder3D JAANET::get ( const JHead & header )

inline

Get cylinder corresponding to can.

Get object from header.

Parameters

header header

Returns: cylinder

Definition at line 335 of file JHeadToolkit.hh.

   {
     using namespace JGEOMETRY2D;
  
     return JCylinder3D(JCircle2D(JVector2D(), header.can.r),
                        header.can.zmin,
                        header.can.zmax);
   }

Variable Documentation

◆ AASHOWER_RECONSTRUCTION_TYPE

const int JAANET::AASHOWER_RECONSTRUCTION_TYPE = 101

static

AAshower reconstruction type for AAnet.

Definition at line 62 of file JAAnetToolkit.hh.

◆ getEvtValue

const JEvtEvaluator JAANET::getEvtValue

static

Function object for evaluation of DAQ objects.

Definition at line 48 of file JEvtEvaluator.hh.

◆ AANET_TAG_SEPARATOR

const char JAANET::AANET_TAG_SEPARATOR = '_'

static

Separator for AAnet tag extension for multiple tags ("_<counter>").

Definition at line 32 of file JHead.hh.

Classes

Enumerations

Functions

Variables

Detailed Description

Enumeration Type Documentation

◆ JHitType_t

◆ JGeant4Type_t

◆ JTrackType_t

Function Documentation

◆ getTime()

◆ getNPE()

◆ is_noise()

◆ getTimeRange() [1/2]

◆ getTimeRange() [2/2]

◆ getPosition() [1/2]

◆ getPosition() [2/2]

◆ getDirection() [1/2]

◆ getDirection() [2/2]

◆ getAxis()

◆ getTrack()

◆ getVertex()

◆ getTransformation()

◆ getHit() [1/3]

◆ getHit() [2/3]

◆ getHit() [3/3]

◆ is_photon()

◆ is_neutrino()

◆ is_electron()

◆ is_muon()

◆ is_tau()

◆ is_pion()

◆ is_proton()

◆ is_neutron()

◆ is_hadron()

◆ has_particleID()

◆ has_neutrino()

◆ get_neutrino()

◆ count_electrons()

◆ has_electron()

◆ get_electron()

◆ from_electron()

◆ count_muons()

◆ has_muon()

◆ get_muon()

◆ from_muon()

◆ count_taus()

◆ has_tau()

◆ get_tau()

◆ from_tau()

◆ count_hadrons()

◆ from_hadron()

◆ get_hadronic_cascade()

◆ add_time_offset()

◆ has_muon_prefit()

◆ has_muon_simplex()

◆ has_muon_gandalf()

◆ has_muon_energy()

◆ has_muon_start()

◆ has_muon_fit()

◆ has_shower_prefit()

◆ has_shower_positionfit()

◆ has_shower_completefit()

◆ has_shower_fit()

◆ has_aashower_fit()

◆ has_reconstructed_track() [1/2]

◆ has_reconstructed_track() [2/2]

◆ has_reconstructed_muon()

◆ has_reconstructed_shower()

◆ has_reconstructed_aashower()

◆ get_best_reconstructed_track() [1/2]

◆ get_best_reconstructed_track() [2/2]

◆ get_best_reconstructed_muon()

◆ get_best_reconstructed_shower()

◆ get_best_reconstructed_aashower()

◆ copy() [1/2]

◆ copy() [2/2]

◆ getTag()

◆ is_valid()

◆ operator==()