Extensions to Evt data format. More...

Classes
struct	JAtmosphericNeutrinoFlux
	Implementation of atmospheric neutrino flux using official KM3NeT atmospheric flux function. More...

struct	JConstantFlux
	Function object for constant fluxes. More...

struct	JDiffuseFlux
	Low-level interface for diffuse fluxes. More...

struct	JDiffuseFluxHelper
	Helper class for diffuse flux factor. More...

struct	JEvtCategory
	Low-level interface for event categories. More...

struct	JEvtCategoryHelper
	Helper class for event categories. More...

class	JEvtCategoryMap
	Auxiliary class for reading a map of event categories. More...

class	JEvtCategoryMap< JEvtWeightFactorProduct >
	Template specialisation for a map between event categories and event-weight factor products. More...

class	JEvtCategoryMap< JFluxHelper >
	Template specialisation for a map between event categories and flux factors. More...

struct	JEvtCategorySet
	Container for a set of event categories. More...

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

struct	JEvtWeight
	Abstract base class for event weighing. More...

struct	JEvtWeightCorsika
	Implementation of event weighting for Corsika data. More...

struct	JEvtWeightDAQ
	Implementation of event weighing for DAQ data. More...

struct	JEvtWeightFactor
	Abstract base class for specifiable event-weight factors. More...

struct	JEvtWeightFactorConstant
	Class for constant event-weight factors. More...

struct	JEvtWeightFactorFunction
	Implementation of event-weight factor interface. More...

struct	JEvtWeightFactorFunction< JDiffuseFluxFunction_t, JDiffuseFlux >
	Implementation of event-weight factor interface for diffuse flux objects. More...

struct	JEvtWeightFactorFunction< Flux_Atmospheric, JDiffuseFlux >
	Specialisation of event-weight factor interface for atmospheric neutrino flux. More...

struct	JEvtWeightFactorFunction< pEvtWeightFactor, JEvtWeightFactor_t >
	Implementation of C-style event-weight factor. More...

struct	JEvtWeightFactorFunction< pDiffuseFlux, JDiffuseFlux >
	Implementation of C-style diffuse flux event-weight factor. More...

class	JEvtWeightFactorGSeaGen
	Implementation of reweighting factor for simulated neutrino interactions according to a specifiable ROOT TFormula. More...

struct	JEvtWeightFactorHelper
	Helper class for event-weight factor. More...

struct	JEvtWeightFactorMupage
	Implementation of reweighting factor for mupage events according to a specifiable ROOT TFormula. More...

struct	JEvtWeightFactorProduct
	Class for product of event-weight factors. More...

class	JEvtWeightFactorTFormula
	Base class implementation for reweighting factor for simulated events according to a specifiable ROOT TFormula. More...

struct	JEvtWeightGenhen
	Implementation of event weighting for Genhen data. More...

struct	JEvtWeightGSeaGen
	Implementation of event weighting for GSeaGen data. More...

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

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

struct	JEvtWeightKM3BUU
	Implementation of event weighting for KM3BUU data. More...

struct	JEvtWeightMiscellaneous
	Implementation of event weighing for miscellaneous data such as a merged offline file containing neutrinos and atmospheric muons. More...

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

struct	JEvtWeightNormalisation
	Auxiliary data structure for storing pairs of header UUIDs and event-weight normalisations. More...

struct	JEvtWeighter
	Look-up table for event weighters. More...

struct	JFlux
	Low-level interface for retrieving the flux corresponding to a given event. More...

struct	JFluxHelper
	Helper class for flux function. More...

struct	JRange_t
	Type definition of range. More...

struct	start_run
	Start of run record. More...

struct	String
	General purpose string class. More...

struct	detector
	Detector file. More...

struct	muon_desc_file
	Muon descriptor file. More...

struct	target
	Target. More...

struct	XSecFile
	Neutrino cross section file. More...

struct	drawing
	Drawing. More...

struct	calibration
	Calibration. More...

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

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

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

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

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

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

struct	physics
	Physics information. More...

struct	simul
	Generator for simulation. More...

struct	spectrum
	Neutrino energy spectrum. More...

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

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

struct	genvol
	Neutrino vertex volume. More...

struct	coord_origin
	Coordinate origin. More...

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

struct	norma
	Normlisation of CORSIKA events. More...

struct	livetime
	Normalisation of MUPAGE events. More...

struct	flux
	Neutrino flux. More...

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

struct	depth
	Depth. More...

struct	DAQ
	Livetime of DAQ data. More...

struct	K40
	Livetime of noise data. More...

struct	tgen
	Time duration of event generation. More...

struct	time_interval
	UTC time interval for event generation. More...

struct	primary
	Primary particle. More...

struct	end_event
	End of event record. More...

class	JHead
	Monte Carlo run header. More...

struct	JHeadHelper
	Auxiliary map of application to check method. More...

struct	JHeadWriter
	Implementation for Head output of JHead objects with ROOT dictionary. More...

struct	JHondaBinRange
	Auxiliary data structure for reading Honda bin ranges. More...

struct	JHondaAngularBinSpecs
	Auxiliary data structure for reading angular binning specifications of Honda flux table. More...

class	JHondaFluxInterpolator
	Template definition for Honda flux table interpolator. More...

class	JHondaFluxInterpolator< JConstantFunction1D< double, JArray< 4, double > >, JCoszFunctionalMap_t, JEnergyFunctionalMap_t >
	Template specialisation for interpolator of azimuth-averaged Honda flux table. More...

struct	JLorentzBoost
	Auxiliary class for performing Lorentz boosts. More...

struct	JMultiHead
	Auxiliary data structure to store multiple headers and bookkeep event-weight normalisations. More...

struct	JMuonBundleCategory
	Class for muon bundle categories. More...

struct	JNeutrinoInteractionCategory
	Class for neutrino interaction categories. More...

struct	JOscFlux
	Implementation of oscillated neutrino flux. More...

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

struct	JPDB
	Collection of particles. More...

struct	JPowerLawFlux
	Example function object for computing power-law flux. More...

struct	JVolume
	Auxiliary class for histogramming of effective volume. More...

struct	JEventSelector
	Event selector. More...

Typedefs
using	JEvtWeightFactorMap = JEvtCategoryMap< JEvtWeightFactorProduct >

using	JFluxMap = JEvtCategoryMap< JFluxHelper >
	Type alias for flux maps. More...

using	JFluxAtmospheric = JEvtWeightFactorFunction< Flux_Atmospheric, JDiffuseFlux >
	Type alias for atmospheric flux function interface. More...

typedef double(*	pEvtWeightFactor) (const Evt &)
	Type definition of event-weight factor pointer. More...

using	pFlux = pEvtWeightFactor
	Type definition of flux function pointer. More...

typedef double(*	pDiffuseFlux) (int, double, double)
	Type definition of pointer to diffuse flux function. More...

typedef bool(*	is_head) (const JHead &)
	Type definition of test function of header. More...

template<template< class, class, class > class JCoszFunctionalMap_t = JPolint1FunctionalMap, template< class, class, class > class JEnergyFunctionalMap_t = JPolint1FunctionalMap>
using	JHondaFluxInterpolator2D = JHondaFluxInterpolator< JConstantFunction1D< double, JArray< 4, double > >, JCoszFunctionalMap_t, JEnergyFunctionalMap_t >
	Alias template definition for 2-dimensional Honda flux interpolator. 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 class	JInteractionTypeGENIE_t { INTERACTION_TYPES_START = -1 , UNDEFINED , ELECTROMAGNETIC , WEAK_CHARGED_CURRENT , WEAK_NEUTRAL_CURRENT , WEAK_MIX , DARK_MATTER , NUCLEON_DECAY , NEUTRON_OSCILLATION , NEUTRAL_HEAVY_LEPTON , DARK_NEUTRAL_CURRENT , INTERACTION_TYPES_END }
	Enumeration of GENIE interaction types. More...

enum class	JScatteringTypeGENIE_t { UNKNOWN = -100 , NEUTRINO_SCATTERING_TYPES_START = -1 , UNDEFINED = 0 , QUASI_ELASTIC , SINGLE_KAON , DEEP_INELASTIC , RESONANT , COHERENT_PRODUCTION , DIFFRACTIVE , ELECTRON_NEUTRINO_ELASTIC , INVERSE_MUON_DECAY , ANOMALY_MEDIATED_NEUTRINO_GAMMA , MESON_EXCHANGE_CURRENT , COHERENT_ELASTIC , INVERSE_BETA_DECAY , GLASHOW_RESONANCE , INVERSE_MUON_DECAY_ANNIHILATION , NEUTRINO_SCATTERING_TYPES_END , DARK_MATTER_SCATTERING_TYPES_START = 100 , DARK_MATTER_ELASTIC , DARK_MATTER_DEEP_INELASTIC , DARK_MATTER_ELECTRON , DARK_MATTER_SCATTERING_TYPES_END }
	Enumeration of GENIE scattering types. 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_D_PLUS = 35 , GEANT4_TYPE_D_MINUS = 36 , GEANT4_TYPE_D_ZERO = 37 , GEANT4_TYPE_ANTID_ZERO = 38 , GEANT4_TYPE_D_PLUS_S = 39 , GEANT4_TYPE_D_MINUS_S = 40 , GEANT4_TYPE_LAMBDA_PLUS_C = 41 , 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 , GEANT4_TYPE_XI_PLUS_C = 0 , GEANT4_TYPE_XI_MINUS_C = 0 , GEANT4_TYPE_XI_ZERO_C = 0 , GEANT4_TYPE_ANTIXI_ZERO_C = 0 , GEANT4_TYPE_OMEGA_PLUS_B = 0 , GEANT4_TYPE_ANTIOMEGA_ZERO_C = 0 , GEANT4_TYPE_OMEGA_MINUS_B = 0 , GEANT4_TYPE_OMEGA_ZERO_C = 0 , GEANT4_TYPE_B_PLUS = 0 , GEANT4_TYPE_B_MINUS = 0 , GEANT4_TYPE_XI_MINUS_B = 0 , GEANT4_TYPE_XI_PLUS_B = 0 , GEANT4_TYPE_B_ZERO = 0 , GEANT4_TYPE_ANTIB_ZERO = 0 , GEANT4_TYPE_B_ZERO_S = 0 , GEANT4_TYPE_ANTIB_ZERO_S = 0 , GEANT4_TYPE_XI_ZERO_B = 0 , GEANT4_TYPE_ANTIXI_ZERO_B = 0 , GEANT4_TYPE_LAMBDA_B = 0 , GEANT4_TYPE_ANTILAMBDA_B = 0 , GEANT4_TYPE_B_PLUS_C = 0 , GEANT4_TYPE_B_MINUS_C = 0 }
	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_D_PLUS = 411 , TRACK_TYPE_D_ZERO = 421 , TRACK_TYPE_D_PLUS_S = 431 , TRACK_TYPE_B_ZERO = 511 , TRACK_TYPE_B_PLUS = 521 , TRACK_TYPE_B_ZERO_S = 531 , TRACK_TYPE_B_PLUS_C = 541 , 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_LAMBDA_PLUS_C = 4122 , TRACK_TYPE_XI_ZERO_C = 4132 , TRACK_TYPE_XI_PLUS_C = 4232 , TRACK_TYPE_OMEGA_ZERO_C = 4332 , TRACK_TYPE_LAMBDA_B = 5122 , TRACK_TYPE_XI_MINUS_B = 5132 , TRACK_TYPE_XI_ZERO_B = 5232 , TRACK_TYPE_OMEGA_MINUS_B = 5332 , 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_D_MINUS = -411 , TRACK_TYPE_D_MINUS_S = -431 , TRACK_TYPE_ANTID_ZERO = -421 , TRACK_TYPE_ANTIB_ZERO = -511 , TRACK_TYPE_B_MINUS = -521 , TRACK_TYPE_ANTIB_ZERO_S = -531 , TRACK_TYPE_B_MINUS_C = -541 , TRACK_TYPE_ANTIPROTON = -2212 , TRACK_TYPE_ANTINEUTRON = -2112 , TRACK_TYPE_ANTILAMBDA = -3122 , TRACK_TYPE_ANTISIGMA_PLUS = -3222 , TRACK_TYPE_NEUTRAL_ANTISIGMA = -3212 , TRACK_TYPE_ANTISIGMA_MINUS = -3112 , TRACK_TYPE_NEUTRAL_ANTIXI = -3322 , TRACK_TYPE_ANTIXI_MINUS = -3312 , TRACK_TYPE_ANTIOMEGA_MINUS = -3334 , TRACK_TYPE_ANTIXI_ZERO_C = -4132 , TRACK_TYPE_XI_MINUS_C = -4232 , TRACK_TYPE_ANTIOMEGA_ZERO_C = -4332 , TRACK_TYPE_ANTILAMBDA_B = -5122 , TRACK_TYPE_XI_PLUS_B = -5132 , TRACK_TYPE_ANTIXI_ZERO_B = -5232 , TRACK_TYPE_OMEGA_PLUS_B = -5332 }
	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. time between earliest and latest hit) of Monte Carlo event. More...

JTimeRange	getTimeRange (const Evt &event, const JTimeRange &T_ns)
	Get time range (i.e. time between earliest and latest hit) of Monte Carlo event. More...

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

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

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

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

Vec	getDirection (const JDirection3D &dir)
	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...

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

double	getW (const Trk &track, const size_t index, const double value)
	Get track information. More...

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

bool	is_muonbundle (const Trk &track)
	Test whether given track is a (anti-)muon. 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_lepton (const Trk &track)
	Test whether given track is a lepton. More...

bool	is_charged_lepton (const Trk &track)
	Test whether given track is a charged lepton. More...

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

bool	is_meson (const Trk &track)
	Test whether given track is a meson (is hadron and third digit of PDG code is zero) More...

bool	is_baryon (const Trk &track)
	Test whether given track is a baryon (is hadron and third digit of PDG code is not zero) More...

bool	is_nucleus (const Trk &track)
	Test whether given track is a nucleus (PDG code corresponds to a baryon or has 10 digits, starting with '10') More...

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

bool	is_initialstate (const Trk &track)
	Test whether given track corresponds to an initial state particle. More...

bool	is_finalstate (const Trk &track)
	Test whether given track corresponds to a final state particle. 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...

bool	has_target_nucleus (const Evt &evt)
	Test whether given event has well-defined target nucleus. More...

const Trk &	get_target_nucleus (const Evt &evt)
	Get target nucleus. More...

const Trk &	get_primary (const Evt &evt)
	Get primary. More...

const Trk &	get_leading_lepton (const Evt &event)
	Auxiliary function to retrieve the leading lepton of a neutrino interaction. More...

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

JVertex3D	getVertex (const Evt &event)
	Get event vertex. 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	has_muonic_taudecay (const Evt &event)
	Test whether given event contains a tau-lepton decay into a muon. More...

bool	has_leptonic_taudecay (const Evt &event)
	Test whether given event contains a leptonic tau-decay. More...

bool	has_hadronic_taudecay (const Evt &event)
	Test whether given event contains a hadronic tau-decay. More...

int	count_taudecay_prongs (const Evt &event)
	Count the number of tau-lepton decay prongs in a given event. 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...

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

double	getKineticEnergy (const Trk &trk)
	Get track kinetic energy. More...

double	getE0 (const Evt &evt)
	Get initial state energy of a neutrino interaction. More...

double	getE1 (const Evt &evt)
	Get final state energy of a neutrino interaction. More...

Vec	getP0 (const Evt &evt)
	Get momentum vector of the initial state of a neutrino interaction. More...

Vec	getP1 (const Evt &evt)
	Get momentum vector of the final state of a neutrino interaction. More...

double	getInvariantMass (const Evt &event)
	Get final state invariant mass. More...

JEvtCategorySet	getCategories (const JHead &header)
	Auxiliary method to retrieve the unique event categories corresponding to a MC header. More...

template<class JFunction_t , class JEvtWeightFactor_t = JEvtWeightFactor>
JEvtWeightFactorFunction< JFunction_t, JEvtWeightFactor_t >	make_weightFactor (const JFunction_t &function)
	Auxiliary method for creating an interface to an event-weight factor. More...

template<class JEvtWeightFactor_t = JEvtWeightFactor>
JEvtWeightFactorFunction< pEvtWeightFactor, JEvtWeightFactor_t >	make_weightFactor (pEvtWeightFactor function)
	Auxiliary method for creating an interface to an event-weight factor. More...

template<class JFunction_t >
JEvtWeightFactorFunction< JFunction_t, JFlux >	make_fluxFunction (const JFunction_t &flux)
	Auxiliary method for creating an interface to a flux function. More...

JEvtWeightFactorFunction< pFlux, JFlux >	make_fluxFunction (pFlux flux)
	Auxiliary method for creating an interface to a flux function. More...

template<class JFunction_t >
JEvtWeightFactorFunction< JFunction_t, JDiffuseFlux >	make_diffuseFluxFunction (const JFunction_t &flux)
	Auxiliary method for creating an interface to a diffuse flux function. More...

JEvtWeightFactorFunction< pDiffuseFlux, JDiffuseFlux >	make_diffuseFluxFunction (pDiffuseFlux flux)
	Auxiliary method for creating an interface to a diffuse flux function. More...

double	getVolume (const JType< JEvtWeightGSeaGen > &type, const Evt &evt)
	Get volume of given event according given weighter. More...

double	getVolume (const JType< JEvtWeightKM3BUU > &type, const Evt &evt)
	Get volume of given event according given weighter. More...

double	getVolume (const JType< JEvtWeightGenhen > &type, const Evt &evt)
	Get volume of given event according given weighter. More...

double	getVolume (const JEvtWeight &weighter, const Evt &evt)
	Get volume of given event according given weighter. More...

JInteractionTypeGENIE_t	getInteractionType (const int intType)
	Auxiliary function to convert interaction types. More...

JScatteringTypeGENIE_t	getScatteringType (const int scatType)
	Auxiliary function to convert scattering types. 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...

std::string	getTag (const std::string &tag, const int counter)
	Get tag with aanet extension "_<counter>" for identical tags. 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_genhen (const JHead &header)
	Check for generator. More...

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

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

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

bool	is_km3buu (const JHead &header)
	Check for generator. 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_sirene (const JHead &header)
	Check for detector simulation. More...

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

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

Vec	getOffset (const JHead &header)
	Get offset. More...

Vec	getOrigin (const JHead &header)
	Get origin. More...

JCylinder3D	getCylinder (const JHead &header)
	Get cylinder corresponding to the positions of generated tracks (i.e. More...

JLorentzBoost	getBoostToCOM (const Evt &event)
	Get Lorentz boost to the Center of Mass (COM) frame for a given neutrino interaction. More...

void	boostToCOM (Evt &event)
	Boost event to the Center of Mass (COM) frame. More...

bool	is_muon_bundle_primary (const int type)
	Auxiliary function to check if given PDG code corresponds to a valid muon bundle primary type. More...

bool	is_neutrino_primary (const int type)
	Auxiliary function to check if given PDG code corresponds to a neutrino. More...

JInteractionTypeGENIE_t	getInteractionType (const JHead &header)
	Auxiliary function to retrieve the GENIE interaction type corresponding to a given header. More...

std::pair< JParameter< int >, JParameter< int > >	getTargetInfo (const JHead &header)
	Auxiliary function to retrieve the target atomic number and mass number. More...

bool	select (const Trk &trk, const Evt &evt)
	Event selection. More...

bool	hasW (const Trk &trk, const int i)
	Check availability of value. More...

double	getW (const Trk &trk, const int i)
	Get associated value. More...

double	getW (const Trk &trk, const int i, const double value)
	Get associated value. More...

void	setW (Trk &trk, const int i, const double value)
	Set associated value. More...

Variables
static const int	AASHOWER_RECONSTRUCTION_TYPE = 101
	AAnet shower fit reconstruction type. More...

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

JEvtWeighter	getEventWeighter
	Function object for mapping header to event weighter. More...

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

static JHeadHelper	get_is_head
	Function object to get check method for given application. More...

static const char *const	HONDA_ENERGY_KEYWORD = "Enu(GeV)"

static const char *const	HONDA_COSINE_ZENITH_ANGLE_KEYWORD = "cosZ"

static const char *const	HONDA_AZIMUTH_ANGLE_KEYWORD = "phi_Az"

static const char *const	HONDA_MUON_NEUTRINO_FLUX_KEYWORD = "NuMu"

static const char *const	HONDA_MUON_ANTINEUTRINO_FLUX_KEYWORD = "NuMubar"

static const char *const	HONDA_ELECTRON_NEUTRINO_FLUX_KEYWORD = "NuE"

static const char *const	HONDA_ELECTRON_ANTINEUTRINO_FLUX_KEYWORD = "NuEbar"

static const char *const	HONDA_BINSPECS_SEPARATOR = "--"

static const char	HONDA_BINSPECS_BEGIN = '['

static const char	HONDA_BINSPECS_END = ']'

static const JRange< double >	DEFAULT_BUNDLE_ENERGY_RANGE = JRange<double>(0, 1e11)
	Default bundle energy range [GeV]. More...

static const JRange< double >	DEFAULT_BUNDLE_COSINE_ZENITH_ANGLE_RANGE = JRange<double>(0.0, 1.0)
	Default bundle cosine zenith angle range [-]. More...

static const JRange< int >	DEFAULT_BUNDLE_MULTIPLICITY_RANGE = JRange<int> (0, 1000)
	Default bundle multiplicity range [-]. More...

static const double	DEFAULT_BUNDLE_RADIAL_EXTENT = 1e5
	Default bundle radial extent [m]. More...

static const char *const	NEUTRAL_CURRENT = "NC"
	Neutral current header string identifier. More...

static const char *const	CHARGED_CURRENT = "CC"
	Charged current header string identifier. More...

static const JRange< double >	DEFAULT_NEUTRINO_ENERGY_RANGE = JRange<double>(0, 1e11)
	Default neutrino energy range [GeV]. More...

static const JRange< double >	DEFAULT_NEUTRINO_COSINE_ZENITH_ANGLE_RANGE = JRange<double>(-1.0, 1.0)
	Default neutrino cosine zenith angle range [-]. More...

Detailed Description

Extensions to Evt data format.

Author: bjung; bjjung; mlincetto; mdejong, bjjung; mdejong; bjung, mdejong

Typedef Documentation

◆ JEvtWeightFactorMap

using JAANET::JEvtWeightFactorMap = typedef JEvtCategoryMap<JEvtWeightFactorProduct>

Definition at line 493 of file JEvtCategoryMap.hh.

◆ JFluxMap

using JAANET::JFluxMap = typedef JEvtCategoryMap<JFluxHelper>

Type alias for flux maps.

Definition at line 807 of file JEvtCategoryMap.hh.

◆ JFluxAtmospheric

using JAANET::JFluxAtmospheric = typedef JEvtWeightFactorFunction<Flux_Atmospheric, JDiffuseFlux>

Type alias for atmospheric flux function interface.

Definition at line 189 of file JEvtWeightFactorFunction.hh.

◆ pEvtWeightFactor

typedef double(* JAANET::pEvtWeightFactor) (const Evt &)

Type definition of event-weight factor pointer.

Definition at line 193 of file JEvtWeightFactorFunction.hh.

◆ pFlux

using JAANET::pFlux = typedef pEvtWeightFactor

Type definition of flux function pointer.

Definition at line 197 of file JEvtWeightFactorFunction.hh.

◆ pDiffuseFlux

typedef double(* JAANET::pDiffuseFlux) (int, double, double)

Type definition of pointer to diffuse flux function.

Definition at line 238 of file JEvtWeightFactorFunction.hh.

◆ is_head

typedef bool(* JAANET::is_head) (const JHead &)

Type definition of test function of header.

Definition at line 28 of file JHeadToolkit.hh.

◆ JHondaFluxInterpolator2D

template<template< class, class, class > class JCoszFunctionalMap_t = JPolint1FunctionalMap, template< class, class, class > class JEnergyFunctionalMap_t = JPolint1FunctionalMap>

using JAANET::JHondaFluxInterpolator2D = typedef JHondaFluxInterpolator<JConstantFunction1D<double, JArray<4, double> >, JCoszFunctionalMap_t, JEnergyFunctionalMap_t>

Alias template definition for 2-dimensional Honda flux interpolator.

Definition at line 911 of file JHondaFluxInterpolator.hh.

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

◆ JInteractionTypeGENIE_t

enum JAANET::JInteractionTypeGENIE_t

strong

Enumeration of GENIE interaction types.

Enumerator
INTERACTION_TYPES_START	Start of interaction types.
UNDEFINED	Unknown interaction type.
ELECTROMAGNETIC	Electromagnetic interaction.
WEAK_CHARGED_CURRENT	Weak charged current interaction.
WEAK_NEUTRAL_CURRENT	Weak neutral current interaction.
WEAK_MIX	CC + NC + interference.
DARK_MATTER	Dark matter interaction.
NUCLEON_DECAY	Nucleon decay.
NEUTRON_OSCILLATION	Neutron oscillation.
NEUTRAL_HEAVY_LEPTON	Heavy neutral lepton.
DARK_NEUTRAL_CURRENT	Dark neutral current.
INTERACTION_TYPES_END	End of interaction types.

Definition at line 20 of file JGENIETypes.hh.

                                      {
     INTERACTION_TYPES_START            =   -1, //!< Start of interaction types
     UNDEFINED,                                 //!< Unknown interaction type
     ELECTROMAGNETIC,                           //!< Electromagnetic interaction
     WEAK_CHARGED_CURRENT,                      //!< Weak charged current interaction
     WEAK_NEUTRAL_CURRENT,                      //!< Weak neutral current interaction
     WEAK_MIX,                                  //!< CC + NC + interference
     DARK_MATTER,                               //!< Dark matter interaction
     NUCLEON_DECAY,                             //!< Nucleon decay
     NEUTRON_OSCILLATION,                       //!< Neutron oscillation
     NEUTRAL_HEAVY_LEPTON,                      //!< Heavy neutral lepton
     DARK_NEUTRAL_CURRENT,                      //!< Dark neutral current
     INTERACTION_TYPES_END                      //!< End of interaction types
   };

◆ JScatteringTypeGENIE_t

enum JAANET::JScatteringTypeGENIE_t

strong

Enumeration of GENIE scattering types.

Enumerator
UNKNOWN	Unknown scattering type.
NEUTRINO_SCATTERING_TYPES_START	Start of standard scattering types.
UNDEFINED	Undefined scattering type.
QUASI_ELASTIC	Quasi-elastic scattering.
SINGLE_KAON	Single kaon scattering.
DEEP_INELASTIC	Deep Inelastic Scattering.
RESONANT	Resonant scattering.
COHERENT_PRODUCTION	Coherent production.
DIFFRACTIVE	Diffractive scattering.
ELECTRON_NEUTRINO_ELASTIC	Electron-neutrino elastic scattering.
INVERSE_MUON_DECAY	Inverse muon decay.
ANOMALY_MEDIATED_NEUTRINO_GAMMA	Anomaly mediated neutrino-gamma scattering.
MESON_EXCHANGE_CURRENT	Meson exchange current.
COHERENT_ELASTIC	Coherent elastic scattering.
INVERSE_BETA_DECAY	Inverse beta-decay.
GLASHOW_RESONANCE	Glashow-resonance.
INVERSE_MUON_DECAY_ANNIHILATION	Inverse muon decay annihilation.
NEUTRINO_SCATTERING_TYPES_END	End of standard scattering types.
DARK_MATTER_SCATTERING_TYPES_START	Start of dark matter scattering types.
DARK_MATTER_ELASTIC	Dark matter elastic scattering.
DARK_MATTER_DEEP_INELASTIC	Dark matter deep inelastic scattering.
DARK_MATTER_ELECTRON	Dark matter electron scattering.
DARK_MATTER_SCATTERING_TYPES_END	End of dark matter scattering types.

Definition at line 39 of file JGENIETypes.hh.

                                     {
     UNKNOWN                            = -100, //!< Unknown   scattering type
     NEUTRINO_SCATTERING_TYPES_START    =   -1, //!< Start of standard scattering types
     UNDEFINED                          =    0, //!< Undefined scattering type
     QUASI_ELASTIC,                             //!< Quasi-elastic scattering
     SINGLE_KAON,                               //!< Single kaon scattering
     DEEP_INELASTIC,                            //!< Deep Inelastic Scattering
     RESONANT,                                  //!< Resonant scattering
     COHERENT_PRODUCTION,                       //!< Coherent production
     DIFFRACTIVE,                               //!< Diffractive scattering
     ELECTRON_NEUTRINO_ELASTIC,                 //!< Electron-neutrino elastic scattering
     INVERSE_MUON_DECAY,                        //!< Inverse muon decay
     ANOMALY_MEDIATED_NEUTRINO_GAMMA,           //!< Anomaly mediated neutrino-gamma scattering
     MESON_EXCHANGE_CURRENT,                    //!< Meson exchange current
     COHERENT_ELASTIC,                          //!< Coherent elastic scattering
     INVERSE_BETA_DECAY,                        //!< Inverse beta-decay
     GLASHOW_RESONANCE,                         //!< Glashow-resonance
     INVERSE_MUON_DECAY_ANNIHILATION,           //!< Inverse muon decay annihilation
     NEUTRINO_SCATTERING_TYPES_END,             //!< End of standard scattering types
     DARK_MATTER_SCATTERING_TYPES_START = 100,  //!< Start of dark matter scattering types
     DARK_MATTER_ELASTIC,                       //!< Dark matter elastic scattering
     DARK_MATTER_DEEP_INELASTIC,                //!< Dark matter deep inelastic scattering
     DARK_MATTER_ELECTRON,                      //!< Dark matter electron scattering
     DARK_MATTER_SCATTERING_TYPES_END           //!< End of dark matter scattering types
   };

◆ 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_D_PLUS
GEANT4_TYPE_D_MINUS
GEANT4_TYPE_D_ZERO
GEANT4_TYPE_ANTID_ZERO
GEANT4_TYPE_D_PLUS_S
GEANT4_TYPE_D_MINUS_S
GEANT4_TYPE_LAMBDA_PLUS_C
GEANT4_TYPE_DEUTERON
GEANT4_TYPE_TRITON
GEANT4_TYPE_ALPHA
GEANT4_TYPE_GEANTINO
GEANT4_TYPE_HE3
GEANT4_TYPE_ANTITAU
GEANT4_TYPE_TAU
GEANT4_TYPE_XI_PLUS_C
GEANT4_TYPE_XI_MINUS_C
GEANT4_TYPE_XI_ZERO_C
GEANT4_TYPE_ANTIXI_ZERO_C
GEANT4_TYPE_OMEGA_PLUS_B
GEANT4_TYPE_ANTIOMEGA_ZERO_C
GEANT4_TYPE_OMEGA_MINUS_B
GEANT4_TYPE_OMEGA_ZERO_C
GEANT4_TYPE_B_PLUS
GEANT4_TYPE_B_MINUS
GEANT4_TYPE_XI_MINUS_B
GEANT4_TYPE_XI_PLUS_B
GEANT4_TYPE_B_ZERO
GEANT4_TYPE_ANTIB_ZERO
GEANT4_TYPE_B_ZERO_S
GEANT4_TYPE_ANTIB_ZERO_S
GEANT4_TYPE_XI_ZERO_B
GEANT4_TYPE_ANTIXI_ZERO_B
GEANT4_TYPE_LAMBDA_B
GEANT4_TYPE_ANTILAMBDA_B
GEANT4_TYPE_B_PLUS_C
GEANT4_TYPE_B_MINUS_C

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_D_PLUS            =  35,
                        GEANT4_TYPE_D_MINUS           =  36,
                        GEANT4_TYPE_D_ZERO            =  37,
                        GEANT4_TYPE_ANTID_ZERO        =  38,
                        GEANT4_TYPE_D_PLUS_S          =  39,
                        GEANT4_TYPE_D_MINUS_S         =  40,
                        GEANT4_TYPE_LAMBDA_PLUS_C     =  41,
                        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, 
                        //
                        // %Particles with PDG code but not GEANT
                        //
                        GEANT4_TYPE_XI_PLUS_C         =   0,
                        GEANT4_TYPE_XI_MINUS_C        =   0,
                        GEANT4_TYPE_XI_ZERO_C         =   0,
                        GEANT4_TYPE_ANTIXI_ZERO_C     =   0,
                        GEANT4_TYPE_OMEGA_PLUS_B      =   0,
                        GEANT4_TYPE_ANTIOMEGA_ZERO_C  =   0,
                        GEANT4_TYPE_OMEGA_MINUS_B     =   0,
                        GEANT4_TYPE_OMEGA_ZERO_C      =   0,
                        GEANT4_TYPE_B_PLUS            =   0,
                        GEANT4_TYPE_B_MINUS           =   0,
                        GEANT4_TYPE_XI_MINUS_B        =   0,
                        GEANT4_TYPE_XI_PLUS_B         =   0,
                        GEANT4_TYPE_B_ZERO            =   0,
                        GEANT4_TYPE_ANTIB_ZERO        =   0,
                        GEANT4_TYPE_B_ZERO_S          =   0,
                        GEANT4_TYPE_ANTIB_ZERO_S      =   0,
                        GEANT4_TYPE_XI_ZERO_B         =   0,
                        GEANT4_TYPE_ANTIXI_ZERO_B     =   0,
                        GEANT4_TYPE_LAMBDA_B          =   0,
                        GEANT4_TYPE_ANTILAMBDA_B      =   0,
                        GEANT4_TYPE_B_PLUS_C          =   0,
                        GEANT4_TYPE_B_MINUS_C         =   0 };

◆ 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_D_PLUS
TRACK_TYPE_D_ZERO
TRACK_TYPE_D_PLUS_S
TRACK_TYPE_B_ZERO
TRACK_TYPE_B_PLUS
TRACK_TYPE_B_ZERO_S
TRACK_TYPE_B_PLUS_C
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_LAMBDA_PLUS_C
TRACK_TYPE_XI_ZERO_C
TRACK_TYPE_XI_PLUS_C
TRACK_TYPE_OMEGA_ZERO_C
TRACK_TYPE_LAMBDA_B
TRACK_TYPE_XI_MINUS_B
TRACK_TYPE_XI_ZERO_B
TRACK_TYPE_OMEGA_MINUS_B
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_D_MINUS
TRACK_TYPE_D_MINUS_S
TRACK_TYPE_ANTID_ZERO
TRACK_TYPE_ANTIB_ZERO
TRACK_TYPE_B_MINUS
TRACK_TYPE_ANTIB_ZERO_S
TRACK_TYPE_B_MINUS_C
TRACK_TYPE_ANTIPROTON
TRACK_TYPE_ANTINEUTRON
TRACK_TYPE_ANTILAMBDA
TRACK_TYPE_ANTISIGMA_PLUS
TRACK_TYPE_NEUTRAL_ANTISIGMA
TRACK_TYPE_ANTISIGMA_MINUS
TRACK_TYPE_NEUTRAL_ANTIXI
TRACK_TYPE_ANTIXI_MINUS
TRACK_TYPE_ANTIOMEGA_MINUS
TRACK_TYPE_ANTIXI_ZERO_C
TRACK_TYPE_XI_MINUS_C
TRACK_TYPE_ANTIOMEGA_ZERO_C
TRACK_TYPE_ANTILAMBDA_B
TRACK_TYPE_XI_PLUS_B
TRACK_TYPE_ANTIXI_ZERO_B
TRACK_TYPE_OMEGA_PLUS_B

Definition at line 96 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_D_PLUS             =  411,
                       TRACK_TYPE_D_ZERO             =  421,
                       TRACK_TYPE_D_PLUS_S           =  431,
                       TRACK_TYPE_B_ZERO             =  511,
                       TRACK_TYPE_B_PLUS             =  521,
                       TRACK_TYPE_B_ZERO_S           =  531,
                       TRACK_TYPE_B_PLUS_C           =  541,
                       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_LAMBDA_PLUS_C      =  4122,
                       TRACK_TYPE_XI_ZERO_C          =  4132,
                       TRACK_TYPE_XI_PLUS_C          =  4232,
                       TRACK_TYPE_OMEGA_ZERO_C       =  4332,
                       TRACK_TYPE_LAMBDA_B           =  5122,
                       TRACK_TYPE_XI_MINUS_B         =  5132,
                       TRACK_TYPE_XI_ZERO_B          =  5232,
                       TRACK_TYPE_OMEGA_MINUS_B      =  5332,
  
                       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_D_MINUS            = -411,
                       TRACK_TYPE_D_MINUS_S          = -431,
                       TRACK_TYPE_ANTID_ZERO         = -421,
                       TRACK_TYPE_ANTIB_ZERO         = -511,
                       TRACK_TYPE_B_MINUS            = -521,
                       TRACK_TYPE_ANTIB_ZERO_S       = -531,
                       TRACK_TYPE_B_MINUS_C          = -541,
                       TRACK_TYPE_ANTIPROTON         = -2212,
                       TRACK_TYPE_ANTINEUTRON        = -2112,
                       TRACK_TYPE_ANTILAMBDA         = -3122,
                       TRACK_TYPE_ANTISIGMA_PLUS     = -3222,
                       TRACK_TYPE_NEUTRAL_ANTISIGMA  = -3212,
                       TRACK_TYPE_ANTISIGMA_MINUS    = -3112,
                       TRACK_TYPE_NEUTRAL_ANTIXI     = -3322,
                       TRACK_TYPE_ANTIXI_MINUS       = -3312,
                       TRACK_TYPE_ANTIOMEGA_MINUS    = -3334,
                       TRACK_TYPE_ANTIXI_ZERO_C      = -4132,
                       TRACK_TYPE_XI_MINUS_C         = -4232,
                       TRACK_TYPE_ANTIOMEGA_ZERO_C   = -4332,
                       TRACK_TYPE_ANTILAMBDA_B       = -5122,
                       TRACK_TYPE_XI_PLUS_B          = -5132,
                       TRACK_TYPE_ANTIXI_ZERO_B      = -5232,
                       TRACK_TYPE_OMEGA_PLUS_B       = -5332 };

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 96 of file JAAnetToolkit.hh.

   {
     return hit.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 108 of file JAAnetToolkit.hh.

   {
     return hit.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 120 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. time between earliest and latest hit) of Monte Carlo event.

Note that the global event time in not taken into account.

Parameters

event event

Returns: time range

Definition at line 133 of file JAAnetToolkit.hh.

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

◆ getTimeRange() [2/2]

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

inline

Get time range (i.e. time between earliest and latest hit) of Monte Carlo event.

The resulting time range is larger than or equal to the given time window.
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 156 of file JAAnetToolkit.hh.

   {
     JTimeRange time_range = getTimeRange(event);
  
     if (!time_range.is_valid()) {
       time_range = T_ns;
     } 
  
     if (time_range.getLength() < T_ns.getLength()) {
       time_range.setLength(T_ns.getLength());
     }
  
     return time_range;
   }

◆ getPosition() [1/3]

JPosition3D JAANET::getPosition ( const Vec & pos )

inline

Get position.

Parameters

pos position

Returns: position

Definition at line 178 of file JAAnetToolkit.hh.

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

◆ getPosition() [2/3]

Vec JAANET::getPosition ( const JPosition3D & pos )

inline

Get position.

Parameters

pos position

Returns: position

Definition at line 190 of file JAAnetToolkit.hh.

   {
     return Vec(pos.getX(), pos.getY(), pos.getZ());
   }

◆ getPosition() [3/3]

JPosition3D JAANET::getPosition ( const Trk & track )

inline

Get position.

Parameters

track track

Returns: position

Definition at line 202 of file JAAnetToolkit.hh.

   {
     return getPosition(track.pos);
   }

◆ getDirection() [1/3]

JDirection3D JAANET::getDirection ( const Vec & dir )

inline

Get direction.

Parameters

dir direction

Returns: direction

Definition at line 214 of file JAAnetToolkit.hh.

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

◆ getDirection() [2/3]

Vec JAANET::getDirection ( const JDirection3D & dir )

inline

Get direction.

Parameters

dir direction

Returns: direction

Definition at line 226 of file JAAnetToolkit.hh.

   {
     return Vec(dir.getDX(), dir.getDY(), dir.getDZ());
   }

◆ getDirection() [3/3]

JDirection3D JAANET::getDirection ( const Trk & track )

inline

Get direction.

Parameters

track track

Returns: direction

Definition at line 238 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 250 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 262 of file JAAnetToolkit.hh.

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

◆ getTransformation()

JTransformation3D JAANET::getTransformation ( const Trk & track )

inline

Get transformation.

Parameters

track track

Returns: transformation

Definition at line 274 of file JAAnetToolkit.hh.

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

◆ getW() [1/3]

double JAANET::getW	(	const Trk &	track,
		const size_t	index,
		const double	value
	)

inline

Get track information.

Parameters

track	track
index	index
value	default value

Returns: actual value

Definition at line 288 of file JAAnetToolkit.hh.

   {
     if (index < track.fitinf.size())
       return track.fitinf[index];
     else
       return value;
   }

◆ 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 303 of file JAAnetToolkit.hh.

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

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

Trk::type

int type

MC: particle type in PDG encoding.

Definition: Trk.hh:24

◆ is_muonbundle()

bool JAANET::is_muonbundle ( 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 312 of file JAAnetToolkit.hh.

312 { return (abs(track.type) == PDG_MUONBUNDLE); }

PDG_MUONBUNDLE

static const int PDG_MUONBUNDLE

muon bundle reached the can level (mupage)

Definition: trkmembers.hh:38

◆ 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 320 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 330 of file JAAnetToolkit.hh.

330 { 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 338 of file JAAnetToolkit.hh.

338 { 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 346 of file JAAnetToolkit.hh.

346 { 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 354 of file JAAnetToolkit.hh.

354 { 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 362 of file JAAnetToolkit.hh.

362 { 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 370 of file JAAnetToolkit.hh.

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

◆ is_lepton()

bool JAANET::is_lepton ( const Trk & track )

inline

Test whether given track is a lepton.

Parameters

track track

Returns: true if lepton; else fails

Definition at line 378 of file JAAnetToolkit.hh.

                                             { return (is_neutrino(track) ||
                                                       is_electron(track) ||
                                                       is_muon    (track) ||
                                                       is_tau     (track)); }

◆ is_charged_lepton()

bool JAANET::is_charged_lepton ( const Trk & track )

inline

Test whether given track is a charged lepton.

Parameters

track track

Returns: true if lepton; else fails

Definition at line 389 of file JAAnetToolkit.hh.

                                                     { return (is_electron(track) ||
                                                               is_muon    (track) ||
                                                               is_tau     (track)); }

◆ 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 399 of file JAAnetToolkit.hh.

399 { return (abs(track.type) != TRACK_TYPE_PHOTON &&

400 !is_lepton(track)); }

JAANET::is_lepton

bool is_lepton(const Trk &track)

Test whether given track is a lepton.

Definition: JAAnetToolkit.hh:378

◆ is_meson()

bool JAANET::is_meson ( const Trk & track )

inline

Test whether given track is a meson (is hadron and third digit of PDG code is zero)

Parameters

track track

Returns: true if hadron; else fails

Definition at line 408 of file JAAnetToolkit.hh.

408 { return (is_hadron(track) &&

409 ((int)(track.type / 1000)) % 10 == 0); }

JAANET::is_hadron

bool is_hadron(const Trk &track)

Test whether given track is a hadron.

Definition: JAAnetToolkit.hh:399

◆ is_baryon()

bool JAANET::is_baryon ( const Trk & track )

inline

Test whether given track is a baryon (is hadron and third digit of PDG code is not zero)

Parameters

track track

Returns: true if hadron; else fails

Definition at line 417 of file JAAnetToolkit.hh.

417 { return (is_hadron(track) &&

418 ((int)(track.type / 1000)) % 10 != 0); }

◆ is_nucleus()

bool JAANET::is_nucleus ( const Trk & track )

inline

Test whether given track is a nucleus (PDG code corresponds to a baryon or has 10 digits, starting with '10')

Parameters

track track

Returns: true if hadron; else fails

Definition at line 426 of file JAAnetToolkit.hh.

   {
     using namespace std;
     
     const string typestr = to_string(abs(track.type));
     
     return (is_baryon(track) || (typestr.size() == 10 && typestr.substr(0,2) == "10"));
   }

◆ 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 442 of file JAAnetToolkit.hh.

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

◆ is_initialstate()

bool JAANET::is_initialstate ( const Trk & track )

inline

Test whether given track corresponds to an initial state particle.

Parameters

track track

Returns: true if track corresponds to an initial state particle; else false

Definition at line 450 of file JAAnetToolkit.hh.

   {
     if        (Evt::ROOT_IO_VERSION >= 14) {
       
       return (track.status == TRK_ST_PRIMARYNEUTRINO ||
               track.status == TRK_ST_PRIMARYCOSMIC   ||
               track.status == TRK_ST_MUONBUNDLE      ||
               track.status == TRK_ST_ININUCLEI);
       
     } else if (Evt::ROOT_IO_VERSION >= 11) {
  
       return (track.mother_id == TRK_MOTHER_UNDEFINED ||
               track.mother_id == TRK_MOTHER_NONE);
  
     } else {
       
       return is_neutrino(track) && track.id == 0;
     }
   }

◆ is_finalstate()

bool JAANET::is_finalstate ( const Trk & track )

inline

Test whether given track corresponds to a final state particle.

Parameters

track track

Returns: true if track corresponds to a final state particle; else false

Definition at line 476 of file JAAnetToolkit.hh.

   {
     if (Evt::ROOT_IO_VERSION >= 15) {
       return  track.status == TRK_ST_FINALSTATE;
     } else {
       return !is_initialstate(track);
     }
   }

◆ 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 491 of file JAAnetToolkit.hh.

   {
     using namespace std;
     
     if (Evt::ROOT_IO_VERSION >= 14) {
  
       vector<Trk>::const_iterator i = find_if(evt.mc_trks.cbegin(), evt.mc_trks.cend(),
                                                    make_predicate(&Trk::status, TRK_ST_PRIMARYNEUTRINO));
       return i != evt.mc_trks.cend();
       
     } else {
       
       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 513 of file JAAnetToolkit.hh.

   {
     using namespace std;
     
     if (Evt::ROOT_IO_VERSION >= 14) {
       
       vector<Trk>::const_iterator i = find_if(evt.mc_trks.cbegin(), evt.mc_trks.cend(),
                                                    make_predicate(&Trk::status, TRK_ST_PRIMARYNEUTRINO));
  
       if (i != evt.mc_trks.cend()) { return *i; }
       
     } else if (!evt.mc_trks.empty() && is_neutrino(evt.mc_trks[0])) {
       
       return evt.mc_trks[0];
     }
       
     THROW(JIndexOutOfRange, "JAANET::get_neutrino(): Event " << evt.id << " has no neutrino."); 
   }

◆ has_target_nucleus()

bool JAANET::has_target_nucleus ( const Evt & evt )

inline

Test whether given event has well-defined target nucleus.

Parameters

evt event

Returns: true if target nucleus is defined; else false

Definition at line 538 of file JAAnetToolkit.hh.

   {
     using namespace std;
     
     if (Evt::ROOT_IO_VERSION >= 14) {
  
       vector<Trk>::const_iterator i = find_if(evt.mc_trks.cbegin(), evt.mc_trks.cend(),
                                                    make_predicate(&Trk::status, TRK_ST_ININUCLEI));
       
       return i != evt.mc_trks.cend();
       
     } else {
       
       return false;
     }
   }

◆ get_target_nucleus()

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

inline

Get target nucleus.

Parameters

evt event

Returns: target nucleus

Definition at line 561 of file JAAnetToolkit.hh.

   {
     using namespace std;
     using namespace JPP;
  
     vector<Trk>::const_iterator i = find_if(evt.mc_trks.cbegin(), evt.mc_trks.cend(),
                                             make_predicate(&Trk::status, TRK_ST_ININUCLEI));
     
     if (i != evt.mc_trks.cend()) {
       return *i;
     } else {
       THROW(JIndexOutOfRange, "JAANET::get_target_nucleus(): Cannot retrieve target nucleus for event " << evt.id << ".");
     }
   }

◆ get_primary()

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

inline

Get primary.

Parameters

evt event

Returns: primary track

Definition at line 583 of file JAAnetToolkit.hh.

   {
     using namespace std;
     using namespace JPP;
  
     for (vector<Trk>::const_iterator i = evt.mc_trks.cbegin(); i != evt.mc_trks.cend(); ++i) {
       if (is_initialstate(*i) && i->status != TRK_ST_ININUCLEI) {
         return *i;
       }
     }
  
     THROW(JIndexOutOfRange, "JAANET::get_primary(): Cannot retrieve primary track for event " << evt.id << ".");
   }

◆ get_leading_lepton()

const Trk& JAANET::get_leading_lepton ( const Evt & event )

inline

Auxiliary function to retrieve the leading lepton of a neutrino interaction.

Parameters

event event

Returns: leading lepton

Definition at line 604 of file JAAnetToolkit.hh.

   {
     using namespace std;
     using namespace JPP;
  
      const Trk& primary = get_primary(event);
     
      if (is_neutrino(primary)) {
       
       for (vector<Trk>::const_iterator i = event.mc_trks.cbegin(); i != event.mc_trks.cend(); ++i) {
         if (is_lepton(*i) && i->mother_id == primary.id) {
           return *i;
         }
       }
  
       THROW(JValueOutOfRange, "get_leading_lepton(): Could not find leading lepton for given event");
       
     } else {
       THROW(JValueOutOfRange, "get_leading_lepton(): Given event does not correspond to a neutrino interaction.");
     }
   }

◆ getVertex() [1/2]

JVertex3D JAANET::getVertex ( const Trk & track )

inline

Get vertex.

Parameters

track track

Returns: vertex

Definition at line 633 of file JAAnetToolkit.hh.

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

◆ getVertex() [2/2]

JVertex3D JAANET::getVertex ( const Evt & event )

inline

Get event vertex.

Parameters

event event

Returns: event vertex

Definition at line 644 of file JAAnetToolkit.hh.

   {
     using namespace std;
     using namespace JPP;
     
     if (has_neutrino(event)) {
  
       const Trk& neutrino = get_neutrino(event);
  
       return getVertex(neutrino);
       
     } else if (!event.mc_trks.empty()) {
  
       vector<Trk>::const_iterator i = find_if(event.mc_trks.begin(), event.mc_trks.end(), is_initialstate);
  
       if (i != event.mc_trks.end()) {
         return getVertex(*i);
       } else {
         THROW(JValueOutOfRange, "getVertex(): No initial state track found.");
       }
       
     } else if (!event.trks.empty()) { // For reconstructed DAQ events
       
       const Trk& track = get_best_reconstructed_track<JPP_RECONSTRUCTION_TYPE>(event);
  
       return getVertex(track);
       
     } else {
  
       THROW(JValueOutOfRange, "getVertex(): Could not find a valid event vertex.");
     }
   }

◆ 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 683 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 694 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 705 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 722 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 733 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 744 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 755 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 772 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 783 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 794 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 805 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.");
   }

◆ has_muonic_taudecay()

bool JAANET::has_muonic_taudecay ( const Evt & event )

inline

Test whether given event contains a tau-lepton decay into a muon.

Parameters

event event

Returns: true if event contains a muonic tau-lepton decay; else false

Definition at line 820 of file JAAnetToolkit.hh.

   {
     using namespace std;
     
     if (has_tau(event)) {
       
       const Trk& tau = get_tau(event);
  
       for (vector<Trk>::const_iterator i = event.mc_trks.cbegin(); i != event.mc_trks.cend(); ++i) {
         if (is_muon(*i) && i->mother_id == tau.id) {
           return true;
         }
       }
     }
  
     return false;
   }

◆ has_leptonic_taudecay()

bool JAANET::has_leptonic_taudecay ( const Evt & event )

inline

Test whether given event contains a leptonic tau-decay.

Parameters

event event

Returns: true if event contains a leptonic tau-decay; else false

Definition at line 844 of file JAAnetToolkit.hh.

   {
     using namespace std;
     
     if (has_tau(event)) {
       
       const Trk& tau = get_tau(event);
  
       for (vector<Trk>::const_iterator i = event.mc_trks.cbegin(); i != event.mc_trks.cend(); ++i) {
         if ((is_electron(*i) || is_muon(*i)) && i->mother_id == tau.id) {
           return true;
         }
       }
     }
  
     return false;
   }

◆ has_hadronic_taudecay()

bool JAANET::has_hadronic_taudecay ( const Evt & event )

inline

Test whether given event contains a hadronic tau-decay.

Parameters

event event

Returns: true if event contains a hadronic tau-decay; else false

Definition at line 868 of file JAAnetToolkit.hh.

   {
     return has_tau(event) && !has_leptonic_taudecay(event);
   }

◆ count_taudecay_prongs()

int JAANET::count_taudecay_prongs ( const Evt & event )

inline

Count the number of tau-lepton decay prongs in a given event.

The number of prongs is defined as the number of charged tau-lepton decay products.

Parameters

event event

Returns: number of tau-lepton decay prongs

Definition at line 880 of file JAAnetToolkit.hh.

   {
     using namespace std;
     
     int n = 0;
  
     if (has_tau(event)) {
  
       const Trk& tau = get_tau(event);
  
       for (vector<Trk>::const_iterator i = event.mc_trks.cbegin(); i != event.mc_trks.cend(); ++i) {
         if (i->mother_id == tau.id && JPDB::getInstance().getPDG(i->type).charge != 0) {
           ++n;
         }
       }
     }
  
     return n;
   }

◆ 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 909 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 920 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 934 of file JAAnetToolkit.hh.

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

◆ 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 948 of file JAAnetToolkit.hh.

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

◆ getKineticEnergy()

double JAANET::getKineticEnergy ( const Trk & trk )

inline

Get track kinetic energy.

Parameters

trk track

Returns: kinetic energy [GeV]

Definition at line 962 of file JAAnetToolkit.hh.

   {
     const double energy = trk.E;
     const double mass   = JPDB::getInstance().getPDG(trk.type).mass;
  
     return getKineticEnergy(energy, mass);
   }

◆ getE0()

double JAANET::getE0 ( const Evt & evt )

inline

Get initial state energy of a neutrino interaction.

This method includes baryon number conservation.

Parameters

evt event

Returns: energy [GeV]

Definition at line 978 of file JAAnetToolkit.hh.

   {
     using namespace std;
     
     const Trk& neutrino = get_neutrino(evt);
     
     double E0 = neutrino.E;
     
     for (vector<Trk>::const_iterator track = evt.mc_trks.cbegin(); track != evt.mc_trks.end(); ++track) {
       
       if (!is_finalstate(*track)) { continue; }
       
       if        (track->type == TRACK_TYPE_NEUTRON     ||
                  track->type == TRACK_TYPE_SIGMA_MINUS ||
                  track->type == TRACK_TYPE_NEUTRAL_SIGMA) {
           
         E0 += MASS_NEUTRON;
           
       } else if (track->type == TRACK_TYPE_PROTON     ||
                  track->type == TRACK_TYPE_LAMBDA     ||
                  track->type == TRACK_TYPE_SIGMA_PLUS) {
         
         E0 += MASS_PROTON;
         
       } else if (track->type == TRACK_TYPE_ANTINEUTRON) {
         
         E0 -= MASS_NEUTRON;
         
       } else if (track->type == TRACK_TYPE_ANTIPROTON ||
                  track->type == TRACK_TYPE_ANTILAMBDA) {
           
         E0 -= MASS_PROTON;
       }
     }
  
     return E0;
   }

◆ getE1()

double JAANET::getE1 ( const Evt & evt )

inline

Get final state energy of a neutrino interaction.

This method includes muon energy loss.

Parameters

evt event

Returns: energy [GeV]

Definition at line 1024 of file JAAnetToolkit.hh.

   {
     using namespace std;
     using namespace JPP;
     
     double E1 = 0.0;
     
     const Trk& neutrino = get_neutrino(evt);
  
     for (vector<Trk>::const_iterator track = evt.mc_trks.cbegin(); track != evt.mc_trks.cend(); ++track) {
  
       if (!is_finalstate(*track)) { continue; }      
         
       if (is_muon(*track)) {
  
         const Trk& mother = ( track->mother_id > TRK_MOTHER_UNDEFINED ?
                               evt.mc_trks[track->mother_id] : neutrino );
  
         const double distance = (track->pos - mother.pos).len();
  
         E1 += gWater(track->E, -distance);
           
       } else {
  
         E1 += track->E;
       }
     }
  
     return E1;
   }

◆ getP0()

Vec JAANET::getP0 ( const Evt & evt )

inline

Get momentum vector of the initial state of a neutrino interaction.

This method assumes neutrino DIS on a stationary nucleus

Parameters

evt event

Returns: energy [GeV]

Definition at line 1063 of file JAAnetToolkit.hh.

   {
     const Trk& neutrino = get_neutrino(evt);
  
     return neutrino.E * neutrino.dir;
   }

◆ getP1()

Vec JAANET::getP1 ( const Evt & evt )

inline

Get momentum vector of the final state of a neutrino interaction.

This method includes muon energy losses.

Parameters

evt event

Returns: final state momentum vector [GeV]

Definition at line 1078 of file JAAnetToolkit.hh.

   {
     using namespace std;
     using namespace JPP;
     
     Vec P1(0,0,0);
  
     const Trk& neutrino = get_neutrino(evt);
  
     for (vector<Trk>::const_iterator track = evt.mc_trks.cbegin(); track != evt.mc_trks.cend(); ++track) {
  
       if (!is_finalstate(*track)) { continue; }
       
       double kineticEnergy = 0.0;
  
       if (is_muon(*track)) {
  
         const Trk& mother = ( track->mother_id > TRK_MOTHER_UNDEFINED ?
                               evt.mc_trks[track->mother_id] : neutrino );
  
         const double distance = (track->pos - mother.pos).len();
         const double energy   = gWater(track->E, -distance);
  
         kineticEnergy = getKineticEnergy(energy, MASS_MUON);
           
       } else {
  
         kineticEnergy = getKineticEnergy(*track);
       }
  
       P1 += kineticEnergy * track->dir;
     }
     
     return P1;
   }

◆ getInvariantMass()

double JAANET::getInvariantMass ( const Evt & event )

inline

Get final state invariant mass.

Parameters

event event

Returns: invariant mass [GeV]

Definition at line 1121 of file JAAnetToolkit.hh.

   {
     using namespace std;
     using namespace JPP;
  
     double Minv = 0.0;
  
     for (vector<Trk>::const_iterator track = event.mc_trks.cbegin(); track != event.mc_trks.cend(); ++track) {
  
       if (!is_finalstate(*track)) { continue; }
  
       const double mass = JPDB::getInstance().getPDG(track->type).mass;
  
       Minv += mass;
     }
  
     return Minv;
   }

◆ getCategories()

JEvtCategorySet JAANET::getCategories ( const JHead & header )

inline

Auxiliary method to retrieve the unique event categories corresponding to a MC header.

Parameters

header MC header

Returns: vector of unique event categories

Definition at line 110 of file JEvtCategoryToolkit.hh.

   {
     using namespace std;
     using namespace JPP;
     
     JEvtCategorySet categories;
     
     if        (is_gseagen(header)) {
  
       const JNeutrinoInteractionCategory cat(header);
       const JEvtCategoryHelper           helper(cat);
  
       categories.insert(helper);
  
       for (vector<JAANET::flux>::const_iterator i = header.flux.cbegin(); i != header.flux.cend(); ++i) {
         
         const JNeutrinoInteractionCategory cat2(header, i->type);
         const JEvtCategoryHelper           helper2(cat2);
  
         categories.insert(helper2);
       }
       
     } else if (is_mupage(header) || is_corsika(header)) {
  
       const JMuonBundleCategory cat(header);
       const JEvtCategoryHelper  helper(cat);
       
       categories.insert(helper);
     }
  
     return categories;
   }

◆ make_weightFactor() [1/2]

template<class JFunction_t , class JEvtWeightFactor_t = JEvtWeightFactor>

JEvtWeightFactorFunction<JFunction_t, JEvtWeightFactor_t> JAANET::make_weightFactor ( const JFunction_t & function )

inline

Auxiliary method for creating an interface to an event-weight factor.

Parameters

function function object

Returns: event-weight factor interface

Definition at line 288 of file JEvtWeightFactorFunction.hh.

                                                                                                                   {
     return JEvtWeightFactorFunction<JFunction_t, JEvtWeightFactor_t>(function);
   }  

◆ make_weightFactor() [2/2]

template<class JEvtWeightFactor_t = JEvtWeightFactor>

JEvtWeightFactorFunction<pEvtWeightFactor, JEvtWeightFactor_t> JAANET::make_weightFactor ( pEvtWeightFactor function )

inline

Auxiliary method for creating an interface to an event-weight factor.

Parameters

function function pointer

Returns: event-weight factor interface

Definition at line 300 of file JEvtWeightFactorFunction.hh.

                                                                                                                      {
     return JEvtWeightFactorFunction<pEvtWeightFactor, JEvtWeightFactor_t>(function);
   }

◆ make_fluxFunction() [1/2]

template<class JFunction_t >

JEvtWeightFactorFunction<JFunction_t, JFlux> JAANET::make_fluxFunction ( const JFunction_t & flux )

inline

Auxiliary method for creating an interface to a flux function.

Parameters

flux	flux function object

Returns: flux function interface

Definition at line 312 of file JEvtWeightFactorFunction.hh.

                                                                                                  {
     return JEvtWeightFactorFunction<JFunction_t, JFlux>(flux);
   }

◆ make_fluxFunction() [2/2]

JEvtWeightFactorFunction<pFlux, JFlux> JAANET::make_fluxFunction ( pFlux flux )

inline

Auxiliary method for creating an interface to a flux function.

Parameters

flux	flux function pointer

Returns: flux function interface

Definition at line 323 of file JEvtWeightFactorFunction.hh.

                                                                               {
     return JEvtWeightFactorFunction<pFlux, JFlux>(flux);
   }

◆ make_diffuseFluxFunction() [1/2]

template<class JFunction_t >

JEvtWeightFactorFunction<JFunction_t, JDiffuseFlux> JAANET::make_diffuseFluxFunction ( const JFunction_t & flux )

inline

Auxiliary method for creating an interface to a diffuse flux function.

Parameters

flux	diffuse flux function object

Returns: diffuse flux function interface

Definition at line 335 of file JEvtWeightFactorFunction.hh.

                                                                                                                {
     return JEvtWeightFactorFunction<JFunction_t, JDiffuseFlux>(flux);
   }

◆ make_diffuseFluxFunction() [2/2]

JEvtWeightFactorFunction<pDiffuseFlux, JDiffuseFlux> JAANET::make_diffuseFluxFunction ( pDiffuseFlux flux )

inline

Auxiliary method for creating an interface to a diffuse flux function.

Parameters

flux	diffuse flux function pointer

Returns: diffuse flux function interface

Definition at line 346 of file JEvtWeightFactorFunction.hh.

                                                                                                           {
     return JEvtWeightFactorFunction<pDiffuseFlux, JDiffuseFlux>(flux);
   }

◆ getVolume() [1/4]

double JAANET::getVolume	(	const JType< JEvtWeightGSeaGen > &	type,
		const Evt &	evt
	)

inline

Get volume of given event according given weighter.

Parameters

type	type
evt	event

Returns: volume [m^3 Gev sr]

Definition at line 99 of file JEvtWeightToolkit.hh.

   {
     return (evt.w[WEIGHTLIST_DIFFERENTIAL_EVENT_RATE] *
             evt.w2list[W2LIST_GSEAGEN_WATER_INT_LEN]  /
             evt.w2list[W2LIST_GSEAGEN_P_EARTH]);
   }

◆ getVolume() [2/4]

double JAANET::getVolume	(	const JType< JEvtWeightKM3BUU > &	type,
		const Evt &	evt
	)

inline

Get volume of given event according given weighter.

Parameters

type	type
evt	event

Returns: volume [m^3 Gev sr]

Definition at line 114 of file JEvtWeightToolkit.hh.

   {
     return (evt.w[WEIGHTLIST_DIFFERENTIAL_EVENT_RATE] *
             evt.w2list[W2LIST_KM3BUU_WATER_INT_LEN]   /
             evt.w2list[W2LIST_KM3BUU_P_EARTH]);
   }

◆ getVolume() [3/4]

double JAANET::getVolume	(	const JType< JEvtWeightGenhen > &	type,
		const Evt &	evt
	)

inline

Get volume of given event according given weighter.

Parameters

type	type
evt	event

Returns: volume [m^3 Gev sr]

Definition at line 129 of file JEvtWeightToolkit.hh.

   {
     using namespace JPHYSICS;
  
     const double l_int = 1.0e-6 * NUCLEON_MOLAR_MASS / (evt.w2list[W2LIST_GENHEN_SIG] * DENSITY_SEA_WATER * AVOGADRO);
     const double year  = 60*60*24*365;       // [s]
  
     return (evt.w[WEIGHTLIST_DIFFERENTIAL_EVENT_RATE] *
             l_int                                     /
             year                                      /
             evt.w2list[W2LIST_GENHEN_P_EARTH]);
   }

◆ getVolume() [4/4]

double JAANET::getVolume	(	const JEvtWeight &	weighter,
		const Evt &	evt
	)

inline

Get volume of given event according given weighter.

The return value should be normalised to

number of generated events;
solid angle covered by the generation; and
bin width of the histogram.

Parameters

weighter	weighter
evt	event

Returns: volume [m^3 GeV sr]

Definition at line 155 of file JEvtWeightToolkit.hh.

   {
     if (dynamic_cast<const JEvtWeightGSeaGen*>(&weighter) != NULL) { return getVolume(JType<JEvtWeightGSeaGen>(), evt); }
     if (dynamic_cast<const JEvtWeightGenhen*> (&weighter) != NULL) { return getVolume(JType<JEvtWeightGenhen> (), evt); }
     if (dynamic_cast<const JEvtWeightKM3BUU*> (&weighter) != NULL) { return getVolume(JType<JEvtWeightKM3BUU> (), evt); }
  
     THROW(JCastException, "No valid weighter.");
   }

◆ getInteractionType() [1/2]

JInteractionTypeGENIE_t JAANET::getInteractionType ( const int intType )

inline

Auxiliary function to convert interaction types.

Parameters

intType GENIE interaction type

Returns: GENIE interaction type

Definition at line 72 of file JGENIETypes.hh.

   {
     if (intType > (int) JInteractionTypeGENIE_t::INTERACTION_TYPES_START &&
         intType < (int) JInteractionTypeGENIE_t::INTERACTION_TYPES_END) {
       return static_cast<JInteractionTypeGENIE_t>(intType);
     } else {
       return JInteractionTypeGENIE_t::UNDEFINED;
     }
   }

◆ getScatteringType()

JScatteringTypeGENIE_t JAANET::getScatteringType ( const int scatType )

inline

Auxiliary function to convert scattering types.

Parameters

scatType GENIE scattering type

Returns: GENIE scattering type

Definition at line 89 of file JGENIETypes.hh.

   {
     if ((scatType > (int) JScatteringTypeGENIE_t::NEUTRINO_SCATTERING_TYPES_START    &&
          scatType < (int) JScatteringTypeGENIE_t::NEUTRINO_SCATTERING_TYPES_END)     ||
         (scatType > (int) JScatteringTypeGENIE_t::DARK_MATTER_SCATTERING_TYPES_START &&
          scatType > (int) JScatteringTypeGENIE_t::DARK_MATTER_SCATTERING_TYPES_END)  ||
         (scatType == (int) JScatteringTypeGENIE_t::UNKNOWN)) {
       return static_cast<JScatteringTypeGENIE_t>(scatType);
     } else {
       return JScatteringTypeGENIE_t::UNDEFINED;
     }
   }

◆ 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 162 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());
  
       const string buffer = trim(i->second);
  
       if (abc.is_valid() && buffer != "") {
  
         JRedirectString redirect(reader, buffer);
  
         reader.getObject(abc);
  
         if (i->first == getTag(i->first)) {
           to.insert(*i);                                     // keep track of parsed tags
         }
  
       } else {
  
         to.insert(*i);                                       // store data of unknown 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 201 of file JHead.cc.

   {
     using namespace std;
     using namespace JPP;
  
     to = Head();
  
     JHeadWriter writer(to, JHead::getEquationParameters(), JAAnetDictionary::getInstance());
  
     JRootWritableClass cls(from);
  
     unique_ptr<TIterator> i(cls.getClass()->GetListOfDataMembers()->MakeIterator());
  
     for (const TDataMember* p; (p = (const TDataMember*) i->Next()) != NULL; ) {
       if (!JRootClass::is_static(*p)) {
         if (from.find(p->GetName())  != from.end()                     ||
             cls.get(*p)              == JRootClass(&JHead::start_run)  ||
             cls.get(*p)              == JRootClass(&JHead::end_event)) {
           writer.put(p->GetName(), cls.get(*p), true);
         }
       }
     }
  
     // 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() [1/2]

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

◆ getTag() [2/2]

std::string JAANET::getTag	(	const std::string &	tag,
		const int	counter
	)

inline

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

Parameters

tag	tag
counter	counter

Returns: tag

Definition at line 122 of file JHead.hh.

   {
     std::ostringstream os;
  
     os << tag << AANET_TAG_SEPARATOR << counter;
  
     return os.str();
   }

◆ 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 1801 of file JHead.hh.

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

◆ operator<()

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

inline

Less than operator.

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

Parameters

first	first header
second	second header

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

Definition at line 1817 of file JHead.hh.

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

◆ is_genhen()

bool JAANET::is_genhen ( const JHead & header )

inline

Check for generator.

Parameters

header header

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

Definition at line 37 of file JHeadToolkit.hh.

   {
     for (const auto& i : header.simul) {
       if (i.program == APPLICATION_GENHEN) {
         return true;
       }
     }
  
     for (const auto& i : header.physics) { // legacy
       if (i.buffer.find(APPLICATION_GENHEN) != std::string::npos) {
         return true;
       }
     }
  
     return false;
   }

◆ is_gseagen()

bool JAANET::is_gseagen ( const JHead & header )

inline

Check for generator.

Parameters

header header

Returns: true if this header is produced by gSeaGen; else false

Definition at line 61 of file JHeadToolkit.hh.

   {
     for (const auto& i : header.simul) {
       if (i.program == APPLICATION_GSEAGEN) {
         return true;
       }
     }
  
     for (const auto& i : header.physics) { // legacy
       if (i.buffer.find(APPLICATION_GSEAGEN) != std::string::npos) {
         return true;
       }
     }
  
     return false;
   }

◆ is_mupage()

bool JAANET::is_mupage ( const JHead & header )

inline

Check for generator.

Parameters

header header

Returns: true if this header is produced by MUPAGE; else false

Definition at line 85 of file JHeadToolkit.hh.

   {
     for (const auto& i : header.simul) {
       if (i.program == APPLICATION_MUPAGE) {
         return true;
       }
     }
  
     return false;
   }

◆ is_corsika()

bool JAANET::is_corsika ( const JHead & header )

inline

Check for generator.

Parameters

header header

Returns: true if this header is produced by Corsika; else false

Definition at line 103 of file JHeadToolkit.hh.

   {
     for (const auto& i : header.simul) {
       if (i.program == APPLICATION_CORSIKA) {
         return true;
       }
     }
  
     return false;
   }

◆ is_km3buu()

bool JAANET::is_km3buu ( const JHead & header )

inline

Check for generator.

Parameters

header header

Returns: true if this header is produced by km3buu; else false

Definition at line 121 of file JHeadToolkit.hh.

   {
     for (const auto& i : header.simul) {
       if (i.program == APPLICATION_KM3BUU) {
         return true;
       }
     }
  
     return false;
   }

◆ is_km3()

bool JAANET::is_km3 ( const JHead & header )

inline

Check for detector simulation.

Parameters

header header

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

Definition at line 139 of file JHeadToolkit.hh.

   {
     for (const auto& i : header.simul) {
       if (i.program == APPLICATION_KM3) {
         return true;
       }
     }
  
     return false;
   }

◆ is_km3sim()

bool JAANET::is_km3sim ( const JHead & header )

inline

Check for detector simulation.

Parameters

header header

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

Definition at line 157 of file JHeadToolkit.hh.

   {
     for (const auto& i : header.simul) {
       if (i.program == APPLICATION_KM3SIM) {
         return true;
       }
     }
  
     return false;
   }

◆ is_sirene()

bool JAANET::is_sirene ( const JHead & header )

inline

Check for detector simulation.

Parameters

header header

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

Definition at line 175 of file JHeadToolkit.hh.

   {
     for (const auto& i : header.simul) {
       if (i.program == APPLICATION_JSIRENE) {
         return true;
       }
     }
  
     return false;
   }

◆ is_pure_noise()

bool JAANET::is_pure_noise ( const JHead & header )

inline

Check for generator.

Parameters

header header

Returns: true if this header is produced by pure noise; else false

Definition at line 193 of file JHeadToolkit.hh.

   {
     return header.K40.livetime_s > 0.0;
   }

◆ 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 205 of file JHeadToolkit.hh.

   {
     return header.DAQ.livetime_s > 0.0;
   }

◆ getOffset()

Vec JAANET::getOffset ( const JHead & header )

inline

Get offset.

For consistency with the position of the detector, a pending offset should be added to the positions of generated tracks (i.e. Evt::mc_trks).

The main logic can be summarised as follows.

if JHead::fixedcan is defined, the values correspond to an offset that has already been applied, (0,0,0) is returned;
else if JHead::can is defined, the values correspond to a z-range that has already been applied, (0,0,-zmin) is returned;
else if JHead::coord_origin is defined, the values correspond to an offset that has not yet been applied, (x,y,z) is returned;
else (0,0,0) is returned;

Parameters

header header

Returns: position

Definition at line 275 of file JHeadToolkit.hh.

   {
     if      (header.is_valid(&JHead::fixedcan))
       return Vec(0.0, 0.0, 0.0);
     else if (header.is_valid(&JHead::can))
       return Vec(0.0, 0.0, -header.can.zmin);
     else if (header.is_valid(&JHead::coord_origin))
       return header.coord_origin;
     else
       return Vec(0.0, 0.0, 0.0);
   }

◆ getOrigin()

Vec JAANET::getOrigin ( const JHead & header )

inline

Get origin.

The origin corresponds to the offset that has already been applied to the positions of generated tracks (i.e. Evt::mc_trks).

The main logic can be summarised as follows.

if JHead::fixedcan is defined, the values correspond to the offset that has already been applied, (x,y,0) is returned;
else if JHead::can is defined, the values correspond to a z-range that has already been applied, (0,0,zmin) is returned;
else if JHead::coord_origin is defined, the values correspond to the offset that has not yet been applied, (0,0,0) is returned;
else (0,0,0) is returned;

Parameters

header header

Returns: position

Definition at line 309 of file JHeadToolkit.hh.

   {
     if      (header.is_valid(&JHead::fixedcan))
       return Vec(header.fixedcan.xcenter, header.fixedcan.ycenter, 0.0);
     else if (header.is_valid(&JHead::can))
       return Vec(0.0, 0.0, header.can.zmin);
     else if (header.is_valid(&JHead::coord_origin))
       return Vec(0.0, 0.0, 0.0);
     else
       return Vec(0.0, 0.0, 0.0);
   }

◆ getCylinder()

JCylinder3D JAANET::getCylinder ( const JHead & header )

inline

Get cylinder corresponding to the positions of generated tracks (i.e.

Evt::mc_trks).

For consistency with the position of the detector, the pending offset (see method JAANET::getOffset) should be added to the position of cylinder.

Parameters

header header

Returns: cylinder

Definition at line 331 of file JHeadToolkit.hh.

   {
     using namespace JPP;
  
     if        (header.is_valid(&JHead::fixedcan)) {
  
       return JCylinder3D(JCircle2D(JVector2D(header.fixedcan.xcenter,
                                              header.fixedcan.ycenter), 
                                    header.fixedcan.radius),
                          header.fixedcan.zmin,
                          header.fixedcan.zmax);
  
     } else if (header.is_valid(&JHead::can)) {
  
       return JCylinder3D(JCircle2D(JVector2D(), header.can.r),
                          header.can.zmin,
                          header.can.zmax);
     } else {
  
       return JCylinder3D();
     }
   }

◆ getBoostToCOM()

JLorentzBoost JAANET::getBoostToCOM ( const Evt & event )

Get Lorentz boost to the Center of Mass (COM) frame for a given neutrino interaction.

Parameters

event event

Definition at line 248 of file JLorentzBoost.hh.

   {
     using namespace JPP;
     
     if (has_neutrino(event)) {
  
       const Trk& nu = get_neutrino(event);
  
       const Vec beta = -nu.E / getE0(event) * nu.dir;
  
       return Boost(beta.x, beta.y, beta.z);
       
     } else {
       THROW(JValueOutOfRange, "getBoostToCOM(): Given event does not correspond to a neutrino interaction.");
     }
   }

◆ boostToCOM()

void JAANET::boostToCOM ( Evt & event )

Boost event to the Center of Mass (COM) frame.

Parameters

event event

Definition at line 271 of file JLorentzBoost.hh.

   {
     const JLorentzBoost boost = getBoostToCOM(event);
     boost(event);
   }

◆ is_muon_bundle_primary()

bool JAANET::is_muon_bundle_primary ( const int type )

inline

Auxiliary function to check if given PDG code corresponds to a valid muon bundle primary type.

Parameters

type PDG type

Returns: true if given PDG type corresponds to a valid muon bundle primary type; else false

Definition at line 56 of file JMuonBundleCategory.hh.

   {
     using namespace std;
     
     const int    absType = abs(type);
     const string typeStr = to_string(absType);
       
     return ( (absType == PDG_MUONBUNDLE         ||                    // MUPAGE
               absType == TRACK_TYPE_MUON)                          ||
              (typeStr.size() == 10 && typeStr.substr(0,2) == "10") || // CORSIKA nuclei
              (absType != TRACK_TYPE_PHOTON      &&                    // CORSIKA single baryon primaries (e.g. proton or neutron)
               absType != TRACK_TYPE_ELECTRON    &&
               absType != TRACK_TYPE_MUON        &&
               absType != TRACK_TYPE_TAU         &&
               absType != TRACK_TYPE_NUE         &&
               absType != TRACK_TYPE_NUMU        &&
               absType != TRACK_TYPE_NUTAU       &&
               ((int)(absType / 1000)) % 10 != 0) );
   }

◆ is_neutrino_primary()

bool JAANET::is_neutrino_primary ( const int type )

inline

Auxiliary function to check if given PDG code corresponds to a neutrino.

Parameters

type PDG type

Returns: true if given PDG type corresponds to a neutrino; else false

Definition at line 57 of file JNeutrinoInteractionCategory.hh.

   {
     const int absType = abs(type);
  
     return (absType == TRACK_TYPE_NUE  ||
             absType == TRACK_TYPE_NUMU ||
             absType == TRACK_TYPE_NUTAU);
   }

◆ getInteractionType() [2/2]

JInteractionTypeGENIE_t JAANET::getInteractionType ( const JHead & header )

inline

Auxiliary function to retrieve the GENIE interaction type corresponding to a given header.

Parameters

header MC header

Returns: GENIE interaction type

Definition at line 73 of file JNeutrinoInteractionCategory.hh.

   {
     using namespace std;
     using namespace JPP;
     
     if (is_gseagen(header)) {
       
       for (vector<JAANET::physics>::const_iterator i = header.physics.cbegin(); i != header.physics.cend(); ++i) {
         
         istringstream iss(i->buffer);
         
         for (string buffer; iss >> buffer ;) {
           
           if (buffer        == string(NEUTRAL_CURRENT)) {
             return JInteractionTypeGENIE_t::WEAK_NEUTRAL_CURRENT;
           } else if (buffer == string(CHARGED_CURRENT)) {
             return JInteractionTypeGENIE_t::WEAK_CHARGED_CURRENT;
           }
         }
       }
     }
  
     return JInteractionTypeGENIE_t::UNDEFINED;
   }

◆ getTargetInfo()

std::pair<JParameter<int>, JParameter<int> > JAANET::getTargetInfo ( const JHead & header )

inline

Auxiliary function to retrieve the target atomic number and mass number.

Parameters

header MC header

Returns: pair of atomic number (Z) and mass number (A)

Definition at line 106 of file JNeutrinoInteractionCategory.hh.

   {
     using namespace std;
     using namespace JPP;
  
     JParameter<int> Z;
     JParameter<int> A;
     
     if (!header.target.buffer.empty()) {
       
       istringstream iss(header.target.buffer);
  
       if (! (iss.peek() == (int) 'A'  &&
              iss.ignore() && iss >> A &&
              iss.peek() == (int) 'Z'  &&
              iss.ignore() && iss >> Z &&
              Z > 0 && A >= Z) ) {
         THROW(JValueOutOfRange, "getTargetInfo(): Invalid target header-field: " << header.target.buffer << endl);
       }
     }
     
     return make_pair(Z,A);
   }

◆ select()

bool JAANET::select	(	const Trk &	trk,
		const Evt &	evt
	)

Event selection.

Parameters

trk	track
evt	event

Returns: true if selected; else false

Definition at line 16 of file JAAnet/event_selector.cc.

   {
     using namespace std;
     using namespace JPP;
     /*
     if (has_muon(evt)) {
  
       Trk muon;
  
       for (const auto& t1 : evt.mc_trks) {
         if (is_muon(t1)) {
           if (t1.E > muon.E) {
             muon = t1;
           }
         }
       }
  
       return (is_muon(muon) && getAngle(getDirection(trk), getDirection(muon)) > 5.0);
     }
     */
     return true;
   }

◆ hasW()

bool JAANET::hasW	(	const Trk &	trk,
		const int	i
	)

inline

Check availability of value.

Parameters

trk	track
i	index

Returns: true if available; else false

Definition at line 102 of file JAAnet/JEvD.cc.

   {
     return (i >= 0 && i < (int) trk.fitinf.size());
   }

◆ getW() [2/3]

double JAANET::getW	(	const Trk &	trk,
		const int	i
	)

inline

Get associated value.

Parameters

trk	track
i	index

Returns: value

Definition at line 115 of file JAAnet/JEvD.cc.

   {
     return trk.fitinf.at(i);
   }

◆ getW() [3/3]

double JAANET::getW	(	const Trk &	trk,
		const int	i,
		const double	value
	)

inline

Get associated value.

Parameters

trk	track
i	index
value	default value

Returns: value

Definition at line 129 of file JAAnet/JEvD.cc.

   {
     if (hasW(trk,i))
       return trk.fitinf.at(i);
     else
       return value;
   }

◆ setW()

void JAANET::setW	(	Trk &	trk,
		const int	i,
		const double	value
	)

Set associated value.

Parameters

trk	track
i	index
value	value

Definition at line 145 of file JAAnet/JEvD.cc.

   {
     if (i >= (int) trk.fitinf.size()) {
       trk.fitinf.resize(i + 1, 0.0);
     }
  
     trk.fitinf[i] = value;
   }

Variable Documentation

◆ AASHOWER_RECONSTRUCTION_TYPE

const int JAANET::AASHOWER_RECONSTRUCTION_TYPE = 101

static

AAnet shower fit reconstruction type.

Definition at line 70 of file JAAnetToolkit.hh.

◆ getEvtValue

const JEvtEvaluator JAANET::getEvtValue

static

Function object for evaluation of DAQ objects.

Definition at line 54 of file JEvtEvaluator.hh.

◆ getEventWeighter

JEvtWeighter JAANET::getEventWeighter

Function object for mapping header to event weighter.

Definition at line 12 of file JEvtWeightToolkit.cc.

◆ AANET_TAG_SEPARATOR

const char JAANET::AANET_TAG_SEPARATOR = '_'

static

Separator for tag extension of multiple tags in Head ("_<counter>").

Definition at line 65 of file JHead.hh.

◆ get_is_head

JHeadHelper JAANET::get_is_head

static

Function object to get check method for given application.

Definition at line 250 of file JHeadToolkit.hh.

◆ HONDA_ENERGY_KEYWORD

const char* const JAANET::HONDA_ENERGY_KEYWORD = "Enu(GeV)"

static

Definition at line 49 of file JHondaFluxInterpolator.hh.

◆ HONDA_COSINE_ZENITH_ANGLE_KEYWORD

const char* const JAANET::HONDA_COSINE_ZENITH_ANGLE_KEYWORD = "cosZ"

static

Definition at line 50 of file JHondaFluxInterpolator.hh.

◆ HONDA_AZIMUTH_ANGLE_KEYWORD

const char* const JAANET::HONDA_AZIMUTH_ANGLE_KEYWORD = "phi_Az"

static

Definition at line 51 of file JHondaFluxInterpolator.hh.

◆ HONDA_MUON_NEUTRINO_FLUX_KEYWORD

const char* const JAANET::HONDA_MUON_NEUTRINO_FLUX_KEYWORD = "NuMu"

static

Definition at line 53 of file JHondaFluxInterpolator.hh.

◆ HONDA_MUON_ANTINEUTRINO_FLUX_KEYWORD

const char* const JAANET::HONDA_MUON_ANTINEUTRINO_FLUX_KEYWORD = "NuMubar"

static

Definition at line 54 of file JHondaFluxInterpolator.hh.

◆ HONDA_ELECTRON_NEUTRINO_FLUX_KEYWORD

const char* const JAANET::HONDA_ELECTRON_NEUTRINO_FLUX_KEYWORD = "NuE"

static

Definition at line 55 of file JHondaFluxInterpolator.hh.

◆ HONDA_ELECTRON_ANTINEUTRINO_FLUX_KEYWORD

const char* const JAANET::HONDA_ELECTRON_ANTINEUTRINO_FLUX_KEYWORD = "NuEbar"

static

Definition at line 56 of file JHondaFluxInterpolator.hh.

◆ HONDA_BINSPECS_SEPARATOR

const char* const JAANET::HONDA_BINSPECS_SEPARATOR = "--"

static

Definition at line 58 of file JHondaFluxInterpolator.hh.

◆ HONDA_BINSPECS_BEGIN

const char JAANET::HONDA_BINSPECS_BEGIN = '['

static

Definition at line 59 of file JHondaFluxInterpolator.hh.

◆ HONDA_BINSPECS_END

const char JAANET::HONDA_BINSPECS_END = ']'

static

Definition at line 60 of file JHondaFluxInterpolator.hh.

◆ DEFAULT_BUNDLE_ENERGY_RANGE

const JRange<double> JAANET::DEFAULT_BUNDLE_ENERGY_RANGE = JRange<double>(0, 1e11)

static

Default bundle energy range [GeV].

Definition at line 44 of file JMuonBundleCategory.hh.

◆ DEFAULT_BUNDLE_COSINE_ZENITH_ANGLE_RANGE

const JRange<double> JAANET::DEFAULT_BUNDLE_COSINE_ZENITH_ANGLE_RANGE = JRange<double>(0.0, 1.0)

static

Default bundle cosine zenith angle range [-].

Definition at line 45 of file JMuonBundleCategory.hh.

◆ DEFAULT_BUNDLE_MULTIPLICITY_RANGE

const JRange<int> JAANET::DEFAULT_BUNDLE_MULTIPLICITY_RANGE = JRange<int> (0, 1000)

static

Default bundle multiplicity range [-].

Definition at line 46 of file JMuonBundleCategory.hh.

◆ DEFAULT_BUNDLE_RADIAL_EXTENT

const double JAANET::DEFAULT_BUNDLE_RADIAL_EXTENT = 1e5

static

Default bundle radial extent [m].

Definition at line 47 of file JMuonBundleCategory.hh.

◆ NEUTRAL_CURRENT

const char* const JAANET::NEUTRAL_CURRENT = "NC"

static

Neutral current header string identifier.

Definition at line 43 of file JNeutrinoInteractionCategory.hh.

◆ CHARGED_CURRENT

const char* const JAANET::CHARGED_CURRENT = "CC"

static

Charged current header string identifier.

Definition at line 44 of file JNeutrinoInteractionCategory.hh.

◆ DEFAULT_NEUTRINO_ENERGY_RANGE

const JRange<double> JAANET::DEFAULT_NEUTRINO_ENERGY_RANGE = JRange<double>(0, 1e11)

static

Default neutrino energy range [GeV].

Definition at line 47 of file JNeutrinoInteractionCategory.hh.

◆ DEFAULT_NEUTRINO_COSINE_ZENITH_ANGLE_RANGE

const JRange<double> JAANET::DEFAULT_NEUTRINO_COSINE_ZENITH_ANGLE_RANGE = JRange<double>(-1.0, 1.0)

static

Default neutrino cosine zenith angle range [-].

Definition at line 48 of file JNeutrinoInteractionCategory.hh.

Classes

Typedefs

Enumerations

Functions

Variables

Detailed Description

Typedef Documentation

◆ JEvtWeightFactorMap

◆ JFluxMap

◆ JFluxAtmospheric

◆ pEvtWeightFactor

◆ pFlux

◆ pDiffuseFlux

◆ is_head

◆ JHondaFluxInterpolator2D

Enumeration Type Documentation

◆ JHitType_t

◆ JInteractionTypeGENIE_t

◆ JScatteringTypeGENIE_t

◆ JGeant4Type_t

◆ JTrackType_t

Function Documentation

◆ getTime()

◆ getNPE()

◆ is_noise()

◆ getTimeRange() [1/2]

◆ getTimeRange() [2/2]

◆ getPosition() [1/3]

◆ getPosition() [2/3]

◆ getPosition() [3/3]

◆ getDirection() [1/3]

◆ getDirection() [2/3]

◆ getDirection() [3/3]

◆ getAxis()

◆ getTrack()

◆ getTransformation()

◆ getW() [1/3]

◆ is_photon()

◆ is_muonbundle()

◆ is_neutrino()

◆ is_electron()

◆ is_muon()

◆ is_tau()

◆ is_pion()

◆ is_proton()

◆ is_neutron()

◆ is_lepton()

◆ is_charged_lepton()

◆ is_hadron()

◆ is_meson()

◆ is_baryon()

◆ is_nucleus()

◆ has_particleID()

◆ is_initialstate()

◆ is_finalstate()

◆ has_neutrino()

◆ get_neutrino()

◆ has_target_nucleus()

◆ get_target_nucleus()

◆ get_primary()

◆ get_leading_lepton()

◆ getVertex() [1/2]

◆ getVertex() [2/2]

◆ count_electrons()

◆ has_electron()

◆ get_electron()

◆ from_electron()

◆ count_muons()

◆ has_muon()

◆ get_muon()

◆ from_muon()

◆ count_taus()

◆ has_tau()

◆ get_tau()

◆ has_muonic_taudecay()

◆ has_leptonic_taudecay()

◆ has_hadronic_taudecay()

◆ count_taudecay_prongs()

◆ from_tau()

◆ count_hadrons()