JAANET Namespace Reference

Extensions to Evt data format. More...

Classes
struct	JEvtEvaluator
	Auxiliary class to determine value of Evt objects. More...
struct	JEvtWeight
	Low-level interface for event weighing. More...
struct	JEvtWeightDAQ
	Implementation of event weighing for DAQ data. More...
struct	JEvtWeightFactor
	Low-level interface for retrieving a specifiable multiplication factor corresponding to a given event. More...
struct	JEvtWeightFactorFunction
	Implementation of event-weight factor interface. More...
struct	JEvtWeightFactorFunction< pEvtWeightFactor >
	Implementation of C-style event-weight factor. More...
struct	JEvtWeightFactorHelper
	Helper class for event-weight factor. More...
struct	JEvtWeightFactorMultiParticle
	Implementation of event-weight factor for multiple particle types. More...
struct	JEvtWeightGSeaGen
	Implementation of event weighting for GSeaGen data. More...
struct	JEvtWeightHelper
	Helper class for event weighing. More...
struct	JEvtWeightKM3BUU
	Implementation of event weighting for KM3BUU data. More...
struct	JEvtWeightMupage
	Implementation of event weighing for MUPAGE data. More...
struct	JEvtWeighter
	Look-up table for event weighters. 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	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	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...
struct	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	JHit_t
	Auxiliary class to set-up Hit. More...
struct	JHits_t
	Auxiliary data structure for list of hits with hit merging capability. More...
struct	JParticle
	Auxiliary class to handle particle name, codes and mass. More...
struct	JPDB
	Collection of particles. More...
struct	JTrk_t
	Auxiliary class to set-up Trk. More...
struct	JVolume
	Auxiliary class for histogramming of effective volume. More...

Typedefs
typedef double(*	pEvtWeightFactor )(const Evt &)
	Type definition of pointer to event-weight factor. More...
using	JFlux = JEvtWeightFactor
	Alias for flux function interface. More...
using	JFluxHelper = JEvtWeightFactorHelper
	Alias for flux function interface helper. More...
using	pFlux = pEvtWeightFactor
	Alias for pointer to flux function. More...
template<class JFunction_t >
using	JFluxFunction = JEvtWeightFactorFunction< JFunction_t >
	Alias for implementation of flux function interface for functors. More...
using	JFluxMultiParticle = JEvtWeightFactorMultiParticle
	Alias for implementation of an interface to multiple flux functions. More...
typedef bool(*	is_head )(const JHead &)
	Type definition of test function of header. More...

Enumerations
enum	JHitType_t {   HIT_TYPE_MUON_DIRECT = +5, HIT_TYPE_MUON_SCATTERED = -5, HIT_TYPE_DELTARAYS_DIRECT = +4, HIT_TYPE_DELTARAYS_SCATTERED = -4,   HIT_TYPE_BREMSSTRAHLUNG_DIRECT = +3, HIT_TYPE_BREMSSTRAHLUNG_SCATTERED = -3, HIT_TYPE_SHOWER_DIRECT = +99, HIT_TYPE_SHOWER_SCATTERED = -99,   HIT_TYPE_NOISE = -1, HIT_TYPE_UNKNOWN = 0 }
	Enumeration of hit types based on km3 codes. More...
enum	JGeant4Type_t {   GEANT4_TYPE_PHOTON = 1, GEANT4_TYPE_ANTIELECTRON = 2, GEANT4_TYPE_ELECTRON = 3, GEANT4_TYPE_NEUTRINO = 4,   GEANT4_TYPE_ANTIMUON = 5, GEANT4_TYPE_MUON = 6, GEANT4_TYPE_NEUTRAL_PION = 7, GEANT4_TYPE_PION_PLUS = 8,   GEANT4_TYPE_PION_MINUS = 9, GEANT4_TYPE_KAON_LONG = 10, GEANT4_TYPE_KAON_PLUS = 11, GEANT4_TYPE_KAON_MINUS = 12,   GEANT4_TYPE_NEUTRON = 13, GEANT4_TYPE_PROTON = 14, GEANT4_TYPE_ANTIPROTON = 15, GEANT4_TYPE_KAON_SHORT = 16,   GEANT4_TYPE_ETA = 17, GEANT4_TYPE_LAMBDA = 18, GEANT4_TYPE_SIGMA_PLUS = 19, GEANT4_TYPE_NEUTRAL_SIGMA = 20,   GEANT4_TYPE_SIGMA_MINUS = 21, GEANT4_TYPE_NEUTRAL_XI = 22, GEANT4_TYPE_XI_MINUS = 23, GEANT4_TYPE_OMEGA_MINUS = 24,   GEANT4_TYPE_ANTINEUTRON = 25, GEANT4_TYPE_ANTILAMBDA = 26, GEANT4_TYPE_ANTISIGMA_MINUS = 27, GEANT4_TYPE_NEUTRAL_ANTISIGMA = 28,   GEANT4_TYPE_ANTISIGMA_PLUS = 29, GEANT4_TYPE_NEUTRAL_ANTIXI = 30, GEANT4_TYPE_ANTIXI_PLUS = 31, GEANT4_TYPE_ANTIOMEGA_PLUS = 32,   GEANT4_TYPE_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_ANTINEUTRAL_SIGMA = -3212, TRACK_TYPE_ANTISIGMA_MINUS = -3112, TRACK_TYPE_ANTINEUTRAL_XI = -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...
JVertex3D	getVertex (const Trk &track)
	Get vertex. 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_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	has_particleID (const Trk &track, const int type)
	Test whether given track corresponds to given particle type. More...
bool	has_neutrino (const Evt &evt)
	Test whether given event has an incoming neutrino. More...
const Trk &	get_neutrino (const Evt &evt)
	Get incoming neutrino. More...
int	count_electrons (const Evt &evt)
	Count the number of electrons in a given event. More...
bool	has_electron (const Evt &evt)
	Test whether given event has an electron. More...
const Trk &	get_electron (const Evt &evt)
	Get first electron from the event tracklist. More...
bool	from_electron (const Hit &hit)
	Test whether given hit was produced by an electron. More...
int	count_muons (const Evt &evt)
	Count the number of muons in a given event. More...
bool	has_muon (const Evt &evt)
	Test whether given event has a muon. More...
const Trk &	get_muon (const Evt &evt)
	Get first muon from the event tracklist. More...
bool	from_muon (const Hit &hit)
	Test whether given hit was produced by a muon. More...
int	count_taus (const Evt &evt)
	Count the number of taus in a given event. More...
bool	has_tau (const Evt &evt)
	Test whether given event has a tau. More...
const Trk &	get_tau (const Evt &evt)
	Get first tau from the event tracklist. More...
bool	from_tau (const Hit &hit)
	Test whether given hit was produced by a tau. More...
int	count_hadrons (const Evt &evt)
	Count the number of hadrons in a given event (not including neutral pions) More...
bool	from_hadron (const Hit &hit)
	Test whether given hit was produced by a hadronic shower. More...
Trk	get_hadronic_cascade (const Evt &evt)
	Auxiliary function to get average direction and total energy of a hadronic cascade. More...
void	add_time_offset (Evt &evt, const double tOff)
	Add time offset to mc event; according to current definition, the absolute time of the event is defined by the track "t" attribute; this could change in the future if the global attribute mc_t is assigned to this purpose. More...
template<class JFunction_t >
JEvtWeightFactorFunction < JFunction_t >	make_weightFactor (const JFunction_t &function)
	Auxiliary method for creating an interface to an event-weight factor. More...
JEvtWeightFactorFunction < pEvtWeightFactor >	make_weightFactor (pEvtWeightFactor function)
	Auxiliary method for creating an interface to an event-weight factor. More...
template<class JFunction_t >
JFluxFunction< JFunction_t >	make_fluxFunction (const JFunction_t &function)
	Auxiliary method for creating an interface to a flux function. More...
JFluxFunction< pFlux >	make_fluxFunction (pFlux function)
	Auxiliary method for creating an interface to a flux function. 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_daq (const JHead &header)
	Check for real data. More...
template<class T >
T	get (const JHead &header)
	Get object from header. More...
template<>
Vec	get (const JHead &header)
	Get position offset of detector due to generator. More...
template<>
JPosition3D	get (const JHead &header)
	Get position offset of detector due to generator. More...
template<>
JCylinder3D	get (const JHead &header)
	Get cylinder corresponding to can. More...

Variables
static const int	RESCALED_WEIGHT_INDEX = 3
	Index of rescaled weight. More...
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...

Detailed Description

Extensions to Evt data format.

Author

mlincetto

mdejong

bjung

Typedef Documentation

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

Type definition of pointer to event-weight factor.

Definition at line 57 of file JEvtWeightFactorFunction.hh.

using JAANET::JFlux = typedef JEvtWeightFactor

Alias for flux function interface.

Definition at line 18 of file JFlux.hh.

using JAANET::JFluxHelper = typedef JEvtWeightFactorHelper

Alias for flux function interface helper.

Definition at line 23 of file JFlux.hh.

using JAANET::pFlux = typedef pEvtWeightFactor

Alias for pointer to flux function.

Definition at line 28 of file JFlux.hh.

template<class JFunction_t >

using JAANET::JFluxFunction = typedef JEvtWeightFactorFunction<JFunction_t>

Alias for implementation of flux function interface for functors.

Definition at line 34 of file JFlux.hh.

using JAANET::JFluxMultiParticle = typedef JEvtWeightFactorMultiParticle

Alias for implementation of an interface to multiple flux functions.

Definition at line 39 of file JFlux.hh.

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

Type definition of test function of header.

Definition at line 27 of file JHeadToolkit.hh.

Enumeration Type Documentation

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

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

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_ANTINEUTRAL_SIGMA
TRACK_TYPE_ANTISIGMA_MINUS
TRACK_TYPE_ANTINEUTRAL_XI
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_ANTINEUTRAL_SIGMA  = -3212,
                      TRACK_TYPE_ANTISIGMA_MINUS    = -3112,
                      TRACK_TYPE_ANTINEUTRAL_XI     = -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

double JAANET::getTime ( const Hit &  hit )

inline

Get true time of hit.

Parameters

hit

hit

Returns

true time of hit

Definition at line 84 of file JAAnetToolkit.hh.

  {
    return hit.t;
  }

double JAANET::getNPE ( const Hit &  hit )

inline

Get true charge of hit.

Parameters

hit

hit

Returns

true charge of hit

Definition at line 96 of file JAAnetToolkit.hh.

  {
    return hit.a;
  }

bool JAANET::is_noise ( const Hit &  hit )

inline

Verify hit origin.

Parameters

hit

hit

Returns

true if noise; else false

Definition at line 108 of file JAAnetToolkit.hh.

  {
    return hit.type == HIT_TYPE_NOISE;
  }

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

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

  {
    JTimeRange time_range = getTimeRange(event);

    if (time_range.getLength() < T_ns.getLength()) {
      time_range.setLength(T_ns.getLength());
    }

    return time_range;
  }

JPosition3D JAANET::getPosition ( const Vec &  pos )

inline

Get position.

Parameters

pos

position

Returns

position

Definition at line 162 of file JAAnetToolkit.hh.

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

Vec JAANET::getPosition ( const JPosition3D &  pos )

inline

Get position.

Parameters

pos

position

Returns

position

Definition at line 174 of file JAAnetToolkit.hh.

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

JPosition3D JAANET::getPosition ( const Trk &  track )

inline

Get position.

Parameters

track

track

Returns

position

Definition at line 186 of file JAAnetToolkit.hh.

  {
    return getPosition(track.pos);
  }

JDirection3D JAANET::getDirection ( const Vec &  dir )

inline

Get direction.

Parameters

dir

direction

Returns

direction

Definition at line 198 of file JAAnetToolkit.hh.

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

Vec JAANET::getDirection ( const JDirection3D &  dir )

inline

Get direction.

Parameters

dir

direction

Returns

direction

Definition at line 210 of file JAAnetToolkit.hh.

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

JDirection3D JAANET::getDirection ( const Trk &  track )

inline

Get direction.

Parameters

track

track

Returns

direction

Definition at line 222 of file JAAnetToolkit.hh.

  {
    return getDirection(track.dir);
  }

JAxis3D JAANET::getAxis ( const Trk &  track )

inline

Get axis.

Parameters

track

track

Returns

axis

Definition at line 234 of file JAAnetToolkit.hh.

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

JTrack3E JAANET::getTrack ( const Trk &  track )

inline

Get track.

Parameters

track

track

Returns

track

Definition at line 246 of file JAAnetToolkit.hh.

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

JVertex3D JAANET::getVertex ( const Trk &  track )

inline

Get vertex.

Parameters

track

track

Returns

vertex

Definition at line 258 of file JAAnetToolkit.hh.

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

JTransformation3D JAANET::getTransformation ( const Trk &  track )

inline

Get transformation.

Parameters

track

track

Returns

transformation

Definition at line 270 of file JAAnetToolkit.hh.

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

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

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

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

                                            { return (    track.type  ==  TRACK_TYPE_PHOTON || 
                                                      abs(track.type) ==  TRACK_TYPE_NEUTRAL_PION); }

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                            { return (abs(track.type) != TRACK_TYPE_PHOTON &&
                                                      !is_lepton(track)); }

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

                                            { return (is_hadron(track) &&
                                                      ((int)(track.type / 1000)) % 10 == 0); }

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

                                            { return (is_hadron(track) &&
                                                      ((int)(track.type / 1000)) % 10 != 0); }

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

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

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

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

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

inline

Get incoming neutrino.

Parameters

evt

event

Returns

neutrino

Definition at line 434 of file JAAnetToolkit.hh.

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

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

  {
    return std::count_if(evt.mc_trks.begin(), evt.mc_trks.end(), is_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 459 of file JAAnetToolkit.hh.

  {
    return count_electrons(evt) != 0;
  }

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

inline

Get first electron from the event tracklist.

Parameters

evt

event

Returns

electron

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

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

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

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

  {
    return std::count_if(evt.mc_trks.begin(), evt.mc_trks.end(), is_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 509 of file JAAnetToolkit.hh.

  {
    return count_muons(evt) != 0;
  }

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

inline

Get first muon from the event tracklist.

Parameters

evt

event

Returns

muon

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

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

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

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

  {
    return std::count_if(evt.mc_trks.begin(), evt.mc_trks.end(), is_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 559 of file JAAnetToolkit.hh.

  {
    return count_taus(evt) != 0;
  }

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

inline

Get first tau from the event tracklist.

Parameters

evt

event

Returns

tau

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

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

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

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

  {
    return std::count_if(evt.mc_trks.begin(), evt.mc_trks.end(), is_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 612 of file JAAnetToolkit.hh.

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

Trk JAANET::get_hadronic_cascade ( const Evt &  evt )

inline

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

Parameters

evt

event

Returns

cascade (default if there are no hadrons)

Definition at line 623 of file JAAnetToolkit.hh.

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

    cascade.dir.normalize();

    return cascade;
  }

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

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

template<class JFunction_t >

JEvtWeightFactorFunction<JFunction_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 101 of file JEvtWeightFactorFunction.hh.

                                                                                              {
    return JEvtWeightFactorFunction<JFunction_t>(function);
  }  

JEvtWeightFactorFunction<pEvtWeightFactor> 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 112 of file JEvtWeightFactorFunction.hh.

                                                                                                 {
    return JEvtWeightFactorFunction<pEvtWeightFactor>(function);
  }

template<class JFunction_t >

JFluxFunction<JFunction_t> JAANET::make_fluxFunction ( const JFunction_t &  function )

inline

Auxiliary method for creating an interface to a flux function.

Parameters

function

function object

Returns

event-weight factor interface

Definition at line 48 of file JFlux.hh.

                                                                                   {
    return JFluxFunction<JFunction_t>(function);
  }  

JFluxFunction<pFlux> JAANET::make_fluxFunction ( pFlux  function )

inline

Auxiliary method for creating an interface to a flux function.

Parameters

function

function pointer

Returns

event-weight factor interface

Definition at line 58 of file JFlux.hh.

                                                                {
    return JFluxFunction<pFlux>(function);
  }

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

Copy header from from to to.

Parameters

from	header
to	header

Definition at line 139 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
      }
    }
  }

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

Copy header from from to to.

Parameters

from	header
to	header

Definition at line 178 of file JHead.cc.

  {
    using namespace std;
    using namespace JPP;

    to = Head();

    JHeadWriter writer(to, JHead::getEquationParameters(), JAAnetDictionary::getInstance());

    JRootWritableClass cls(from);

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

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

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

  {
    std::ostringstream os;

    os << tag << AANET_TAG_SEPARATOR << counter;

    return os.str();
  }

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

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

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

inline

Less than operator.

Parameters

first	first header
second	second header

Returns

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

Definition at line 1603 of file JHead.hh.

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

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

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

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

    return false;
  }

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

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

    return false;
  }

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

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

    return false;
  }

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

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

    return false;
  }

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

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

    return false;
  }

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

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

    return false;
  }

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

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

    return false;
  }

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

  {
    return header.DAQ.livetime_s > 0.0;
  }

template<class T >

T JAANET::get ( const JHead &  header )

inline

Get object from header.

Parameters

header

header

Returns

object

Get object from header.

Parameters

header

header

Returns

position

Get object from header.

Parameters

header

header

Returns

cylinder

Definition at line 254 of file JHeadToolkit.hh.

  {
    if (is_sirene(header) || is_km3(header)) {

      return header.coord_origin;

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

template<>

Vec JAANET::get ( const JHead &  header )

inline

Get position offset of detector due to generator.

Get object from header.

Parameters

header

header

Returns

position

Definition at line 254 of file JHeadToolkit.hh.

  {
    if (is_sirene(header) || is_km3(header)) {

      return header.coord_origin;

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

template<>

JPosition3D JAANET::get ( const JHead &  header )

inline

Get position offset of detector due to generator.

Get object from header.

Parameters

header

header

Returns

position

Definition at line 279 of file JHeadToolkit.hh.

  {
    const Vec pos = get<Vec>(header);

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

template<>

JCylinder3D JAANET::get ( const JHead &  header )

inline

Get cylinder corresponding to can.

Get object from header.

Parameters

header

header

Returns

cylinder

Definition at line 294 of file JHeadToolkit.hh.

  {
    using namespace JGEOMETRY2D;

    return JCylinder3D(JCircle2D(JVector2D(), header.can.r),
                       header.can.zmin,
                       header.can.zmax);
  }

Variable Documentation

const int JAANET::RESCALED_WEIGHT_INDEX = 3

static

Index of rescaled weight.

Definition at line 54 of file JAAnetToolkit.hh.

const int JAANET::AASHOWER_RECONSTRUCTION_TYPE = 101

static

AAnet shower fit reconstruction type.

Definition at line 59 of file JAAnetToolkit.hh.

const JEvtEvaluator JAANET::getEvtValue

static

Function object for evaluation of DAQ objects.

Definition at line 48 of file JEvtEvaluator.hh.

JEvtWeighter JAANET::getEventWeighter

Function object for mapping header to event weighter.

Definition at line 5 of file JEvtWeightToolkit.cc.

const char JAANET::AANET_TAG_SEPARATOR = '_'

static

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

Definition at line 35 of file JHead.hh.

JHeadHelper JAANET::get_is_head

static

Initial value:

= {
    std::make_pair(APPLICATION_GENHEN,   is_genhen),
    std::make_pair(APPLICATION_GSEAGEN,  is_gseagen),
    std::make_pair(APPLICATION_MUPAGE,   is_mupage),
    std::make_pair(APPLICATION_CORSIKA,  is_corsika),
    std::make_pair(APPLICATION_KM3BUU,   is_km3buu),
    std::make_pair(APPLICATION_KM3,      is_km3),
    std::make_pair(APPLICATION_KM3SIM,   is_km3sim),
    std::make_pair(APPLICATION_JSIRENE,  is_sirene),
    std::make_pair("DAQ",                is_daq)
  }

Function object to get check method for given application.

Definition at line 224 of file JHeadToolkit.hh.