cvn Namespace Reference

Defines an enumeration for prong classification. More...

Classes
class	Boundary
class	CVNAddTrainingData
class	CVNCosmicEventDump
class	CVNCosmicMapper
class	CVNEventDump
class	CVNEventProngDump
class	CVNMapper
class	EventLabeledPMaps
	The EventLabeledPMaps objects contains a PixelMap and the output class type, and any other bit that goes into the ANN. More...
class	Features
	Features, basic output of CVN neural net. More...
class	ParticlesTrainingData
	The TrainingData objects contains a PixelMap and the output class type, and any other bit that goes into the ANN. More...
class	PixelMap
	PixelMap, basic input to CVN neural net. More...
class	PixelMapProducer
	Producer algorithm for PixelMap, input to CVN neural net. More...
class	ProngSummary
class	ProngTrainingData
	The TrainingData objects contains a PixelMap and the output class type, and any other bit that goes into the ANN. More...
class	RegCVNImageUtils
	Class containing some utility functions for all things RegCVN. More...
class	RegCVNMapper
class	RegHadronResult
	Hadronic energy output from RegHadCVN. More...
class	RegNuResult
	Neutrino energy output from RegCVN. More...
class	RegPixelMap
	RegPixelMap, basic input to CVN neural net. More...
class	RegPixelMapProducer
	Producer algorithm for RegPixelMap, input to CVN neural net. More...
class	RegProngResult
	Electron energy output from RegCVN. More...
class	RegResult
	Result, basic output of Regression CVN neural net. More...
class	Result
	Result, basic output of CVN neural net. More...
class	Results
	Result, basic output of Regression CVN neural net. More...
class	TrainingData
	The TrainingData objects contains a PixelMap and the output class type, and any other bit that goes into the ANN. More...

Typedefs
typedef enum cvn::FinalStateProngs	FinalStateProngsType
typedef enum cvn::FinalState	FinalStateType
typedef enum cvn::HType	HitType
typedef enum cvn::Interaction	InteractionType
typedef enum cvn::ParentParticle	ParentParticleType
typedef enum cvn::Particles	ParticlesType
typedef enum cvn::PType	ProngType
typedef std::vector< std::vector< float > >	ViewVectorF
	Useful typedefs. More...
typedef std::vector< ViewVectorF >	ImageVectorF
typedef enum cvn::SimpleNeutrino	SimpleNeutrinoType

Enumerations
enum	FinalStateProngs {   kNumu1prong, kNumu2prong, kNumu3prong, kNumuMprong,   kNue1prong, kNue2prong, kNue3prong, kNueMprong,   kNC1prong, kNC2prong, kNC3prong, kNCMprong,   kCosmicFSp, kOtherFSp, kNFSpType }
enum	FinalState {   kNumu0tr0sh, kNumu0tr1sh, kNumu0tr2sh, kNumu0trMsh,   kNumu1tr0sh, kNumu1tr1sh, kNumu1tr2sh, kNumu1trMsh,   kNumu2tr0sh, kNumu2tr1sh, kNumu2tr2sh, kNumu2trMsh,   kNumuMtr0sh, kNumuMtr1sh, kNumuMtr2sh, kNumuMtrMsh,   kNue0tr0sh, kNue0tr1sh, kNue0tr2sh, kNue0trMsh,   kNue1tr0sh, kNue1tr1sh, kNue1tr2sh, kNue1trMsh,   kNue2tr0sh, kNue2tr1sh, kNue2tr2sh, kNue2trMsh,   kNueMtr0sh, kNueMtr1sh, kNueMtr2sh, kNueMtrMsh,   kNC0tr0sh, kNC0tr1sh, kNC0tr2sh, kNC0trMsh,   kNC1tr0sh, kNC1tr1sh, kNC1tr2sh, kNC1trMsh,   kNC2tr0sh, kNC2tr1sh, kNC2tr2sh, kNC2trMsh,   kNCMtr0sh, kNCMtr1sh, kNCMtr2sh, kNCMtrMsh,   kCosmicFS, kOtherFS, kNFStType }
enum	HType {   kEmptyHit, kElectronHit, kMuonHit, kProtonHit,   kNeutronHit, kPionHit, kPiZeroHit, kGammaHit,   kMostlyHA, kMostlyEM, kMostlyMU, kUnknownHit }
enum	Interaction {   kNumuQE, kNumuRes, kNumuDIS, kNumuOther,   kNueQE, kNueRes, kNueDIS, kNueOther,   kNutauQE, kNutauRes, kNutauDIS, kNutauOther,   kNuElectronElastic, kNC, kCosmic, kOther,   kNIntType }
enum	ParentParticle {   kPiNuCCPP, kKaNuCCPP, kMuNuCCPP, kPiNuNCPP,   kKaNuNCPP, kMuNuNCPP, kCosmicPP, kOtherPP,   kUnknownPP, kNPPtType }
enum	Particles {   nE_e_1_mu_0_pi0_0_pi_0_n_0_p_1, nM_e_0_mu_1_pi0_0_pi_0_n_0_p_1, nM_e_0_mu_1_pi0_0_pi_1_n_0_p_1, nU_e_0_mu_0_pi0_1_pi_0_n_0_p_1,   nE_e_1_mu_0_pi0_0_pi_1_n_0_p_1, nE_e_1_mu_0_pi0_0_pi_0_n_2_p_2, nU_e_0_mu_0_pi0_0_pi_1_n_0_p_1, nE_e_1_mu_0_pi0_0_pi_1_n_1_p_0,   nM_e_0_mu_1_pi0_0_pi_1_n_1_p_0, nU_e_0_mu_0_pi0_0_pi_1_n_1_p_0, nU_e_0_mu_0_pi0_1_pi_0_n_1_p_0, nM_e_0_mu_1_pi0_0_pi_0_n_2_p_2,   nU_e_0_mu_0_pi0_1_pi_0_n_1_p_1, nM_e_0_mu_1_pi0_0_pi_0_n_1_p_1, nE_e_1_mu_0_pi0_0_pi_0_n_1_p_1, nU_e_0_mu_0_pi0_1_pi_0_n_2_p_2,   nE_e_1_mu_0_pi0_0_pi_0_n_1_p_2, nE_e_1_mu_0_pi0_0_pi_1_n_2_p_2, nM_e_0_mu_1_pi0_0_pi_1_n_1_p_1, nE_e_1_mu_0_pi0_0_pi_0_n_0_p_2,   nU_e_0_mu_0_pi0_0_pi_2_n_0_p_1, nM_e_0_mu_1_pi0_0_pi_0_n_1_p_2, nU_e_0_mu_0_pi0_0_pi_0_n_2_p_2, nM_e_0_mu_1_pi0_0_pi_0_n_0_p_2,   nE_e_1_mu_0_pi0_1_pi_0_n_0_p_1, nE_e_1_mu_0_pi0_0_pi_1_n_1_p_1, nU_e_0_mu_0_pi0_1_pi_1_n_1_p_0, nU_e_0_mu_0_pi0_0_pi_1_n_1_p_1,   nM_e_0_mu_1_pi0_1_pi_0_n_0_p_1, nU_e_0_mu_0_pi0_0_pi_2_n_1_p_0, nU_e_0_mu_0_pi0_0_pi_1_n_2_p_2, nE_e_1_mu_0_pi0_1_pi_1_n_0_p_1,   nU_e_0_mu_0_pi0_1_pi_2_n_0_p_1, nM_e_0_mu_1_pi0_0_pi_1_n_2_p_2, nE_e_1_mu_0_pi0_0_pi_2_n_0_p_1, nU_e_0_mu_0_pi0_1_pi_2_n_2_p_2,   nU_e_0_mu_0_pi0_1_pi_2_n_1_p_0, nU_e_0_mu_0_pi0_1_pi_1_n_0_p_1, nU_e_0_mu_0_pi0_2_pi_2_n_0_p_1, nU_e_0_mu_0_pi0_0_pi_0_n_0_p_1,   nU_e_0_mu_0_pi0_0_pi_2_n_2_p_2, nU_e_0_mu_0_pi0_0_pi_2_n_1_p_1, nM_e_0_mu_1_pi0_0_pi_1_n_0_p_2, nU_e_0_mu_0_pi0_1_pi_0_n_0_p_0,   nE_e_1_mu_0_pi0_0_pi_1_n_0_p_2, nU_e_0_mu_0_pi0_0_pi_1_n_0_p_2, nT_e_0_mu_0_pi0_0_pi_2_n_0_p_1, nU_e_0_mu_0_pi0_1_pi_0_n_2_p_0,   nU_e_0_mu_0_pi0_0_pi_1_n_2_p_0, nM_e_0_mu_1_pi0_1_pi_1_n_0_p_1, nE_e_1_mu_0_pi0_1_pi_2_n_0_p_1, nU_e_0_mu_0_pi0_2_pi_2_n_1_p_0,   nM_e_0_mu_1_pi0_0_pi_1_n_2_p_1, nU_e_0_mu_0_pi0_1_pi_0_n_0_p_2, nU_e_0_mu_0_pi0_0_pi_1_n_2_p_1, nM_e_0_mu_1_pi0_1_pi_2_n_0_p_1,   nM_e_0_mu_1_pi0_0_pi_2_n_1_p_0, nE_e_1_mu_0_pi0_0_pi_2_n_1_p_0, nU_e_0_mu_0_pi0_0_pi_0_n_1_p_0, nM_e_0_mu_1_pi0_0_pi_2_n_0_p_1,   nM_e_0_mu_1_pi0_0_pi_1_n_1_p_2, nU_e_0_mu_0_pi0_2_pi_0_n_0_p_1, nU_e_0_mu_0_pi0_1_pi_1_n_2_p_2, nU_e_0_mu_0_pi0_1_pi_0_n_2_p_1,   nU_e_0_mu_0_pi0_0_pi_1_n_1_p_2, nE_e_1_mu_0_pi0_1_pi_2_n_1_p_0, nE_e_1_mu_0_pi0_0_pi_0_n_2_p_1, nT_e_0_mu_0_pi0_0_pi_2_n_1_p_0,   nM_e_0_mu_1_pi0_1_pi_0_n_0_p_2, nE_e_1_mu_0_pi0_0_pi_1_n_2_p_1, nU_e_0_mu_0_pi0_2_pi_2_n_1_p_1, nM_e_0_mu_1_pi0_1_pi_1_n_1_p_0,   nM_e_0_mu_1_pi0_0_pi_2_n_1_p_1, nU_e_0_mu_0_pi0_0_pi_0_n_2_p_1, nM_e_0_mu_1_pi0_1_pi_0_n_2_p_2, nE_e_1_mu_0_pi0_0_pi_2_n_1_p_1,   nU_e_0_mu_0_pi0_2_pi_2_n_2_p_2, nT_e_0_mu_0_pi0_2_pi_2_n_0_p_1, nM_e_0_mu_1_pi0_1_pi_1_n_1_p_1, nU_e_0_mu_0_pi0_1_pi_2_n_1_p_1,   nU_e_0_mu_0_pi0_0_pi_2_n_2_p_1, nU_e_0_mu_0_pi0_0_pi_2_n_0_p_2, nM_e_0_mu_1_pi0_0_pi_0_n_2_p_1, nE_e_1_mu_0_pi0_1_pi_0_n_1_p_1,   nM_e_0_mu_1_pi0_1_pi_0_n_1_p_1, nE_e_1_mu_0_pi0_1_pi_1_n_1_p_0, nE_e_1_mu_0_pi0_0_pi_2_n_2_p_2, nU_e_0_mu_0_pi0_1_pi_1_n_1_p_1,   nE_e_1_mu_0_pi0_2_pi_2_n_2_p_2, nE_e_1_mu_0_pi0_2_pi_2_n_0_p_1, nE_e_1_mu_0_pi0_1_pi_1_n_1_p_1, nU_e_0_mu_0_pi0_1_pi_1_n_2_p_0,   nU_e_0_mu_0_pi0_1_pi_1_n_1_p_2, nU_e_0_mu_0_pi0_0_pi_2_n_1_p_2, nU_e_0_mu_0_pi0_0_pi_0_n_1_p_2, nT_e_0_mu_0_pi0_2_pi_2_n_2_p_2,   nT_e_0_mu_0_pi0_2_pi_2_n_1_p_0, nT_e_0_mu_0_pi0_1_pi_2_n_0_p_1, nM_e_0_mu_1_pi0_1_pi_2_n_1_p_0, nE_e_1_mu_0_pi0_1_pi_0_n_0_p_2,   nM_e_0_mu_1_pi0_1_pi_2_n_1_p_1, nM_e_0_mu_1_pi0_0_pi_1_n_2_p_0, nU_e_0_mu_0_pi0_1_pi_2_n_2_p_0, nU_e_0_mu_0_pi0_0_pi_0_n_0_p_0,   nT_e_0_mu_0_pi0_1_pi_2_n_1_p_0, nM_e_0_mu_1_pi0_2_pi_2_n_0_p_1, nE_e_1_mu_0_pi0_0_pi_2_n_1_p_2, nE_e_1_mu_0_pi0_0_pi_1_n_0_p_0,   nU_e_0_mu_0_pi0_2_pi_1_n_2_p_2, nU_e_0_mu_0_pi0_2_pi_1_n_0_p_1, nT_e_0_mu_0_pi0_0_pi_2_n_1_p_1, nM_e_0_mu_1_pi0_0_pi_2_n_0_p_2,   nE_e_1_mu_0_pi0_2_pi_2_n_1_p_0, nE_e_1_mu_0_pi0_0_pi_1_n_1_p_2, nU_e_0_mu_0_pi0_2_pi_0_n_2_p_2, nU_e_0_mu_0_pi0_0_pi_2_n_2_p_0,   nU_e_0_mu_0_pi0_0_pi_0_n_1_p_1, nM_e_0_mu_1_pi0_1_pi_2_n_0_p_2, nM_e_0_mu_1_pi0_0_pi_1_n_0_p_0, nE_e_1_mu_0_pi0_1_pi_0_n_2_p_2,   nU_e_0_mu_0_pi0_2_pi_2_n_0_p_2, nU_e_0_mu_0_pi0_1_pi_1_n_0_p_2, nM_e_0_mu_1_pi0_1_pi_2_n_2_p_0, nM_e_0_mu_1_pi0_0_pi_2_n_2_p_2,   nE_e_1_mu_0_pi0_0_pi_1_n_2_p_0, nU_e_0_mu_0_pi0_2_pi_0_n_1_p_1, nM_e_0_mu_1_pi0_2_pi_2_n_2_p_2, nE_e_1_mu_0_pi0_1_pi_2_n_2_p_2,   nE_e_1_mu_0_pi0_1_pi_2_n_2_p_1, nE_e_1_mu_0_pi0_0_pi_0_n_1_p_0, nU_e_0_mu_0_pi0_2_pi_2_n_2_p_0, nU_e_0_mu_0_pi0_2_pi_0_n_1_p_0,   nU_e_0_mu_0_pi0_1_pi_2_n_2_p_1, nU_e_0_mu_0_pi0_1_pi_2_n_1_p_2, nU_e_0_mu_0_pi0_1_pi_1_n_2_p_1, nU_e_0_mu_0_pi0_0_pi_0_n_2_p_0,   nT_e_0_mu_0_pi0_0_pi_1_n_0_p_1, nM_e_0_mu_1_pi0_1_pi_2_n_2_p_2, nM_e_0_mu_1_pi0_0_pi_0_n_1_p_0, nE_e_1_mu_0_pi0_2_pi_1_n_0_p_1,   nE_e_1_mu_0_pi0_1_pi_1_n_2_p_2, nE_e_1_mu_0_pi0_1_pi_0_n_1_p_2, nU_e_0_mu_0_pi0_1_pi_0_n_1_p_2, nU_e_0_mu_0_pi0_0_pi_2_n_0_p_0,   nT_e_0_mu_0_pi0_1_pi_2_n_2_p_2, nT_e_0_mu_0_pi0_0_pi_2_n_2_p_2, nM_e_0_mu_1_pi0_2_pi_2_n_1_p_0, nM_e_0_mu_1_pi0_1_pi_1_n_0_p_2,   nM_e_0_mu_1_pi0_1_pi_0_n_1_p_2, nE_e_1_mu_0_pi0_1_pi_2_n_1_p_2, nE_e_1_mu_0_pi0_1_pi_1_n_0_p_2, nU_e_0_mu_0_pi0_2_pi_1_n_1_p_0,   nT_e_0_mu_1_pi0_0_pi_0_n_0_p_1, nM_e_0_mu_1_pi0_2_pi_2_n_1_p_1, nM_e_0_mu_1_pi0_2_pi_1_n_0_p_1, nM_e_0_mu_1_pi0_1_pi_1_n_2_p_2,   nM_e_0_mu_1_pi0_1_pi_1_n_1_p_2, nM_e_0_mu_1_pi0_0_pi_2_n_2_p_1, nE_e_1_mu_0_pi0_2_pi_1_n_2_p_2, nE_e_1_mu_0_pi0_0_pi_2_n_0_p_2,   nE_e_1_mu_0_pi0_0_pi_0_n_0_p_0, nU_e_0_mu_0_pi0_2_pi_2_n_2_p_1, nU_e_0_mu_0_pi0_2_pi_2_n_1_p_2, nU_e_0_mu_0_pi0_2_pi_1_n_2_p_1,   nU_e_0_mu_0_pi0_2_pi_0_n_1_p_2, nU_e_0_mu_0_pi0_1_pi_2_n_0_p_2, nU_e_0_mu_0_pi0_0_pi_0_n_0_p_2, nM_e_0_mu_1_pi0_2_pi_0_n_0_p_1,   nE_e_1_mu_0_pi0_2_pi_2_n_1_p_1, nE_e_1_mu_0_pi0_2_pi_0_n_0_p_1, nE_e_1_mu_0_pi0_1_pi_1_n_2_p_1, nT_e_0_mu_1_pi0_0_pi_1_n_0_p_1,   nM_e_0_mu_1_pi0_2_pi_0_n_0_p_2, nM_e_0_mu_1_pi0_0_pi_2_n_2_p_0, nM_e_0_mu_1_pi0_0_pi_2_n_1_p_2, nE_e_1_mu_0_pi0_2_pi_2_n_2_p_1,   nE_e_1_mu_0_pi0_2_pi_2_n_0_p_2, nE_e_1_mu_0_pi0_1_pi_2_n_0_p_2, nE_e_1_mu_0_pi0_1_pi_1_n_1_p_2, nT_e_0_mu_0_pi0_1_pi_1_n_0_p_1,   nT_e_0_mu_0_pi0_0_pi_2_n_0_p_2, nM_e_0_mu_1_pi0_1_pi_2_n_2_p_1, nE_e_1_mu_0_pi0_2_pi_0_n_0_p_2, nE_e_1_mu_0_pi0_1_pi_0_n_2_p_1,   nE_e_1_mu_0_pi0_0_pi_2_n_2_p_0, nE_e_1_mu_0_pi0_0_pi_0_n_2_p_0, nU_e_0_mu_0_pi0_2_pi_1_n_1_p_1, nT_e_0_mu_1_pi0_0_pi_1_n_1_p_0,   nT_e_0_mu_0_pi0_0_pi_2_n_2_p_0, nM_e_0_mu_1_pi0_1_pi_2_n_1_p_2, nE_e_1_mu_0_pi0_2_pi_2_n_1_p_2, nE_e_1_mu_0_pi0_2_pi_1_n_1_p_0,   nU_e_0_mu_0_pi0_2_pi_0_n_2_p_1, nU_e_0_mu_0_pi0_1_pi_2_n_0_p_0, nT_e_1_mu_0_pi0_0_pi_1_n_1_p_0, nT_e_1_mu_0_pi0_0_pi_1_n_0_p_1,   nT_e_0_mu_0_pi0_2_pi_2_n_0_p_2, nT_e_0_mu_0_pi0_2_pi_1_n_0_p_1, nT_e_0_mu_0_pi0_1_pi_2_n_1_p_1, nM_e_0_mu_1_pi0_2_pi_2_n_0_p_2,   nM_e_0_mu_1_pi0_2_pi_1_n_1_p_0, nE_e_1_mu_0_pi0_1_pi_2_n_1_p_1, nE_e_1_mu_0_pi0_0_pi_2_n_0_p_0, nU_e_0_mu_0_pi0_2_pi_1_n_2_p_0,   nU_e_0_mu_0_pi0_2_pi_1_n_0_p_2, nU_e_0_mu_0_pi0_2_pi_0_n_2_p_0, nT_e_1_mu_0_pi0_1_pi_2_n_1_p_0, nT_e_1_mu_0_pi0_1_pi_2_n_0_p_1,   nT_e_1_mu_0_pi0_0_pi_0_n_0_p_1, nT_e_0_mu_1_pi0_0_pi_2_n_2_p_2, nT_e_0_mu_1_pi0_0_pi_2_n_0_p_1, nT_e_0_mu_0_pi0_2_pi_2_n_1_p_1,   nT_e_0_mu_0_pi0_1_pi_2_n_2_p_1, nT_e_0_mu_0_pi0_1_pi_2_n_1_p_2, nT_e_0_mu_0_pi0_0_pi_2_n_1_p_2, nT_e_0_mu_0_pi0_0_pi_2_n_0_p_0,   nT_e_0_mu_0_pi0_0_pi_1_n_2_p_2, nM_e_0_mu_1_pi0_2_pi_1_n_0_p_2, nM_e_0_mu_1_pi0_1_pi_1_n_2_p_1, nM_e_0_mu_1_pi0_1_pi_1_n_2_p_0,   nM_e_0_mu_1_pi0_1_pi_0_n_1_p_0, nM_e_0_mu_1_pi0_0_pi_2_n_0_p_0, nE_e_1_mu_0_pi0_2_pi_1_n_1_p_2, nE_e_1_mu_0_pi0_2_pi_1_n_1_p_1,   nE_e_1_mu_0_pi0_2_pi_1_n_0_p_2, nE_e_1_mu_0_pi0_2_pi_0_n_1_p_1, nE_e_1_mu_0_pi0_0_pi_2_n_2_p_1, nU_e_0_mu_0_pi0_2_pi_2_n_0_p_0,   nT_e_1_mu_0_pi0_0_pi_2_n_0_p_1, nT_e_1_mu_0_pi0_0_pi_0_n_1_p_1, nT_e_0_mu_0_pi0_2_pi_2_n_0_p_0, nM_e_0_mu_1_pi0_2_pi_2_n_2_p_0,   nM_e_0_mu_1_pi0_2_pi_1_n_2_p_2, nM_e_0_mu_1_pi0_2_pi_1_n_1_p_2, nM_e_0_mu_1_pi0_2_pi_1_n_1_p_1, nM_e_0_mu_1_pi0_2_pi_0_n_2_p_2,   nM_e_0_mu_1_pi0_2_pi_0_n_1_p_2, nE_e_1_mu_0_pi0_2_pi_2_n_2_p_0, nE_e_1_mu_0_pi0_2_pi_0_n_2_p_2, nE_e_1_mu_0_pi0_2_pi_0_n_2_p_1,   nE_e_1_mu_0_pi0_2_pi_0_n_1_p_2, nE_e_1_mu_0_pi0_1_pi_2_n_2_p_0, nU_e_0_mu_0_pi0_2_pi_1_n_1_p_2, nU_e_0_mu_0_pi0_1_pi_1_n_0_p_0,   nU_e_0_mu_0_pi0_0_pi_1_n_0_p_0, nT_e_1_mu_0_pi0_1_pi_2_n_1_p_2, nT_e_1_mu_0_pi0_0_pi_2_n_2_p_0, nT_e_1_mu_0_pi0_0_pi_1_n_2_p_2,   nT_e_0_mu_1_pi0_2_pi_2_n_1_p_0, nT_e_0_mu_1_pi0_2_pi_2_n_0_p_1, nT_e_0_mu_1_pi0_1_pi_2_n_0_p_1, nT_e_0_mu_1_pi0_1_pi_0_n_0_p_1,   nT_e_0_mu_1_pi0_0_pi_2_n_1_p_0, nT_e_0_mu_1_pi0_0_pi_2_n_0_p_2, nT_e_0_mu_1_pi0_0_pi_0_n_2_p_2, nT_e_0_mu_1_pi0_0_pi_0_n_0_p_2,   nT_e_0_mu_0_pi0_2_pi_2_n_2_p_0, nT_e_0_mu_0_pi0_2_pi_2_n_1_p_2, nT_e_0_mu_0_pi0_2_pi_1_n_0_p_2, nT_e_0_mu_0_pi0_1_pi_2_n_0_p_2,   nT_e_0_mu_0_pi0_1_pi_1_n_1_p_2, nT_e_0_mu_0_pi0_1_pi_1_n_1_p_1, nT_e_0_mu_0_pi0_1_pi_1_n_0_p_2, nT_e_0_mu_0_pi0_0_pi_1_n_1_p_2,   nT_e_0_mu_0_pi0_0_pi_1_n_0_p_0, nM_e_0_mu_1_pi0_2_pi_2_n_2_p_1, nM_e_0_mu_1_pi0_2_pi_2_n_1_p_2, nM_e_0_mu_1_pi0_2_pi_0_n_1_p_1,   nM_e_0_mu_1_pi0_1_pi_0_n_2_p_1, nM_e_0_mu_1_pi0_1_pi_0_n_0_p_0, nM_e_0_mu_1_pi0_0_pi_0_n_0_p_0, nE_e_1_mu_0_pi0_1_pi_2_n_0_p_0,   nU_e_0_mu_0_pi0_2_pi_0_n_0_p_2, nT_e_1_mu_0_pi0_1_pi_0_n_0_p_1, nT_e_1_mu_0_pi0_0_pi_0_n_0_p_0, nT_e_0_mu_1_pi0_1_pi_2_n_2_p_2,   nT_e_0_mu_1_pi0_1_pi_1_n_1_p_0, nT_e_0_mu_0_pi0_2_pi_2_n_2_p_1, nT_e_0_mu_0_pi0_1_pi_1_n_2_p_2, nT_e_0_mu_0_pi0_1_pi_1_n_1_p_0,   nT_e_0_mu_0_pi0_0_pi_1_n_1_p_1, nT_e_0_mu_0_pi0_0_pi_1_n_0_p_2, nM_e_0_mu_1_pi0_2_pi_1_n_2_p_1, nM_e_0_mu_1_pi0_2_pi_1_n_2_p_0,   nM_e_0_mu_1_pi0_1_pi_2_n_0_p_0, nM_e_0_mu_1_pi0_0_pi_0_n_2_p_0, nE_e_1_mu_0_pi0_2_pi_1_n_2_p_1, nE_e_1_mu_0_pi0_1_pi_1_n_2_p_0,   nE_e_1_mu_0_pi0_1_pi_1_n_0_p_0, nE_e_1_mu_0_pi0_1_pi_0_n_0_p_0, nU_e_1_mu_0_pi0_0_pi_2_n_1_p_1, nU_e_0_mu_0_pi0_2_pi_1_n_0_p_0,   nT_e_1_mu_0_pi0_2_pi_2_n_1_p_1, nT_e_1_mu_0_pi0_2_pi_2_n_1_p_0, nT_e_1_mu_0_pi0_2_pi_2_n_0_p_1, nT_e_1_mu_0_pi0_2_pi_1_n_1_p_1,   nT_e_1_mu_0_pi0_1_pi_2_n_1_p_1, nT_e_1_mu_0_pi0_1_pi_2_n_0_p_2, nT_e_1_mu_0_pi0_1_pi_1_n_1_p_1, nT_e_1_mu_0_pi0_1_pi_1_n_0_p_1,   nT_e_1_mu_0_pi0_0_pi_0_n_2_p_2, nT_e_0_mu_1_pi0_2_pi_2_n_2_p_2, nT_e_0_mu_1_pi0_2_pi_2_n_2_p_1, nT_e_0_mu_1_pi0_2_pi_0_n_1_p_2,   nT_e_0_mu_1_pi0_1_pi_2_n_1_p_1, nT_e_0_mu_1_pi0_1_pi_2_n_1_p_0, nT_e_0_mu_1_pi0_1_pi_0_n_2_p_2, nT_e_0_mu_1_pi0_1_pi_0_n_2_p_1,   nT_e_0_mu_1_pi0_0_pi_2_n_1_p_1, nT_e_0_mu_1_pi0_0_pi_1_n_2_p_2, nT_e_0_mu_1_pi0_0_pi_1_n_2_p_0, nT_e_0_mu_1_pi0_0_pi_1_n_1_p_2,   nT_e_0_mu_1_pi0_0_pi_1_n_1_p_1, nT_e_0_mu_1_pi0_0_pi_0_n_1_p_1, nT_e_0_mu_0_pi0_2_pi_1_n_2_p_2, nT_e_0_mu_0_pi0_1_pi_2_n_2_p_0,   nT_e_0_mu_0_pi0_1_pi_2_n_0_p_0, nT_e_0_mu_0_pi0_1_pi_1_n_2_p_0, nT_e_0_mu_0_pi0_0_pi_2_n_2_p_1, nT_e_0_mu_0_pi0_0_pi_1_n_2_p_1,   nT_e_0_mu_0_pi0_0_pi_1_n_2_p_0, nT_e_0_mu_0_pi0_0_pi_1_n_1_p_0, nM_e_1_mu_1_pi0_0_pi_1_n_0_p_1, nM_e_0_mu_1_pi0_1_pi_1_n_0_p_0,   nE_e_1_mu_0_pi0_2_pi_1_n_0_p_0, nU_e_1_mu_0_pi0_0_pi_0_n_0_p_0, nT_e_1_mu_0_pi0_2_pi_2_n_2_p_2, nT_e_1_mu_0_pi0_2_pi_2_n_0_p_2,   nT_e_1_mu_0_pi0_2_pi_1_n_1_p_2, nT_e_1_mu_0_pi0_2_pi_1_n_1_p_0, nT_e_1_mu_0_pi0_2_pi_1_n_0_p_1, nT_e_1_mu_0_pi0_2_pi_1_n_0_p_0,   nT_e_1_mu_0_pi0_2_pi_0_n_1_p_1, nT_e_1_mu_0_pi0_2_pi_0_n_0_p_1, nT_e_1_mu_0_pi0_1_pi_1_n_2_p_2, nT_e_1_mu_0_pi0_1_pi_1_n_1_p_2,   nT_e_1_mu_0_pi0_1_pi_1_n_1_p_0, nT_e_1_mu_0_pi0_1_pi_0_n_2_p_2, nT_e_1_mu_0_pi0_1_pi_0_n_2_p_1, nT_e_1_mu_0_pi0_1_pi_0_n_1_p_1,   nT_e_1_mu_0_pi0_1_pi_0_n_0_p_2, nT_e_1_mu_0_pi0_0_pi_2_n_2_p_1, nT_e_1_mu_0_pi0_0_pi_2_n_1_p_2, nT_e_1_mu_0_pi0_0_pi_2_n_1_p_0,   nT_e_1_mu_0_pi0_0_pi_2_n_0_p_0, nT_e_1_mu_0_pi0_0_pi_1_n_2_p_1, nT_e_1_mu_0_pi0_0_pi_1_n_2_p_0, nT_e_1_mu_0_pi0_0_pi_1_n_1_p_1,   nT_e_1_mu_0_pi0_0_pi_0_n_2_p_1, nT_e_1_mu_0_pi0_0_pi_0_n_1_p_2, nT_e_1_mu_0_pi0_0_pi_0_n_1_p_0, nT_e_1_mu_0_pi0_0_pi_0_n_0_p_2,   nT_e_0_mu_1_pi0_2_pi_2_n_2_p_0, nT_e_0_mu_1_pi0_2_pi_2_n_1_p_2, nT_e_0_mu_1_pi0_2_pi_2_n_1_p_1, nT_e_0_mu_1_pi0_2_pi_2_n_0_p_2,   nT_e_0_mu_1_pi0_2_pi_1_n_1_p_0, nT_e_0_mu_1_pi0_2_pi_1_n_0_p_1, nT_e_0_mu_1_pi0_2_pi_0_n_2_p_2, nT_e_0_mu_1_pi0_2_pi_0_n_0_p_1,   nT_e_0_mu_1_pi0_1_pi_2_n_1_p_2, nT_e_0_mu_1_pi0_1_pi_2_n_0_p_2, nT_e_0_mu_1_pi0_1_pi_1_n_2_p_1, nT_e_0_mu_1_pi0_1_pi_1_n_0_p_1,   nT_e_0_mu_1_pi0_0_pi_2_n_2_p_0, nT_e_0_mu_1_pi0_0_pi_2_n_1_p_2, nT_e_0_mu_1_pi0_0_pi_1_n_2_p_1, nT_e_0_mu_1_pi0_0_pi_1_n_0_p_2,   nT_e_0_mu_1_pi0_0_pi_0_n_1_p_2, nT_e_0_mu_0_pi0_2_pi_1_n_1_p_2, nT_e_0_mu_0_pi0_2_pi_1_n_1_p_1, nT_e_0_mu_0_pi0_2_pi_1_n_0_p_0,   nT_e_0_mu_0_pi0_2_pi_0_n_1_p_2, nT_e_0_mu_0_pi0_1_pi_1_n_2_p_1, nT_e_0_mu_0_pi0_1_pi_0_n_0_p_1, nT_e_0_mu_0_pi0_0_pi_0_n_0_p_1,   nM_e_0_mu_1_pi0_2_pi_2_n_0_p_0, nM_e_0_mu_1_pi0_2_pi_0_n_2_p_1, nM_e_0_mu_1_pi0_2_pi_0_n_2_p_0, nM_e_0_mu_1_pi0_2_pi_0_n_0_p_0,   nM_e_0_mu_1_pi0_1_pi_0_n_2_p_0, nE_e_1_mu_1_pi0_1_pi_0_n_1_p_0, nE_e_1_mu_0_pi0_2_pi_2_n_0_p_0, nE_e_1_mu_0_pi0_2_pi_1_n_2_p_0,   nE_e_1_mu_0_pi0_2_pi_0_n_2_p_0, nE_e_1_mu_0_pi0_1_pi_0_n_2_p_0, nE_e_1_mu_0_pi0_1_pi_0_n_1_p_0, nM_Other,   nE_Other, nT_Other, kCosmic_PT, kOther_PT,   kNIntType_PT }
enum	PType {   kElectron, kMuon, kProton, kNeutron,   kPion, kPiZero, kGamma, kOtherPDG,   kUnknown, kEmpty, kEM, kHadron }
enum	SimpleNeutrino {   kNumu, kNue, kNutau, kNuNC,   kCos, kNSimpleNeutrinoType }

Functions
std::ostream &	operator<< (std::ostream &os, const PixelMapProducer &p)
std::ostream &	operator<< (std::ostream &os, const RegPixelMapProducer &p)
std::map< int, int >	GetMultiplicityMap (const art::Ptr< simb::MCTruth > truth, std::vector< int > PDGs, bool *isnuint, bool *iscc, int *nupdg, int *trueMode)
FinalStateType	GetFinalStateType (const art::Ptr< simb::MCTruth > truth)
FinalStateProngsType	GetFinalStateProngsType (const art::Ptr< simb::MCTruth > truth)
InteractionType	GetInteractionType (const simb::MCTruth &truth)
InteractionType	GetInteractionTypeFromSlice (int pdg, bool iscc, int trueMode)
InteractionType	SliceClassify (const art::Ptr< simb::MCTruth > truth, const rb::Cluster &slice, float &nuEnergy, float &lepEnergy, double &vtxx, double &vtxy, double &vtxz)
ParticlesType	MakeParticleLabel (const art::Ptr< simb::MCTruth > truth, std::string NuType)
ParticlesType	ParticlesSliceClassify (const art::Ptr< simb::MCTruth > truth, const rb::Cluster &slice, float &nuEnergy, float &lepEnergy, double &vtxx, double &vtxy, double &vtxz)
HitType	GetHitTypeByPDG (int pdg)
HitType	HitClassify (art::Ptr< rb::CellHit > Hit, HitType *hType, double *hPurity)
int	HitObject (art::Ptr< rb::CellHit > Hit, int &hObjectCount, std::map< int, int > &hObjectMap)
int	HitNuIndex (art::Ptr< rb::CellHit > Hit, const rb::Cluster &evtCluster)
ParentParticleType	GetParentParticleType (const art::Ptr< simb::MCTruth > truth, const art::Event &evt)
ProngType	GetTypeByPDG (int pdg)
ProngType	ProngClassify (const rb::Prong &prong, ProngType *pType3D, ProngType *pTypeX, ProngType *pTypeY, bool *isprimary, double *purity3D, double *purityX, double *purityY, double *recE, unsigned int *ncellX, unsigned int *ncellY)
std::vector< std::pair< int, double > >	GetProngPurityByPDG (const rb::Cluster &prong)
std::vector< std::pair< int, double > >	GetProngEfficiencyByPDG (const rb::Cluster &prong, const rb::Cluster &slice)
std::vector< std::pair< int, double > >	GetProngEnergyByPDG (const rb::Cluster &prong)
std::vector< std::pair< int, double > >	Get2DProngPurityByPDG (const rb::Cluster &prong)
std::vector< std::pair< int, double > >	Get2DProngEfficiencyByPDG (const rb::Cluster &prong, const rb::Cluster &slice)
std::vector< std::pair< int, double > >	Get2DProngEnergyByPDG (const rb::Cluster &prong)
std::ostream &	operator<< (std::ostream &os, const Boundary &b)
std::ostream &	operator<< (std::ostream &os, const PixelMap &m)
int	GetPDGByPType (PType ptype)
std::ostream &	operator<< (std::ostream &os, const RegPixelMap &m)
SimpleNeutrinoType	GetSimpleNeutrinoType (const InteractionType &intType)

Detailed Description

Defines an enumeration for prong classification.

The ProngTrainingData objects contains a PixelMap and the output type within the chosen classification which will later go into the CNN.

The ProngTrainingData object contains a PixelMap and the output class type, and any other bit that goes into the ANN.

The TrainingData objects contains a PixelMap and the output class type, and any other bit that goes into the ANN.

Defines an enumeration for int type by parent particle type.

Defines an enumeration for interaction type.

Defines an enumeration for cellhit classification.

Defines an enumeration for int type by final state.

Defines an enumeration for int type by final state nprongs.

The EventLabeledPMaps objects contains a PixelMap and the output class type, and any other bit that goes into the ANN.

The EventLabeledPMaps object contains all pixel maps and corresponding labels for selected slices an prongs in the event. These are the input for the CNN.

Typedef Documentation

typedef enum cvn::FinalStateProngs cvn::FinalStateProngsType

typedef enum cvn::FinalState cvn::FinalStateType

typedef enum cvn::HType cvn::HitType

typedef std::vector<ViewVectorF> cvn::ImageVectorF

Definition at line 19 of file RegCVNImageUtils.h.

typedef enum cvn::Interaction cvn::InteractionType

typedef enum cvn::ParentParticle cvn::ParentParticleType

typedef enum cvn::Particles cvn::ParticlesType

typedef enum cvn::PType cvn::ProngType

typedef enum cvn::SimpleNeutrino cvn::SimpleNeutrinoType

typedef std::vector<std::vector<float> > cvn::ViewVectorF

Useful typedefs.

Definition at line 18 of file RegCVNImageUtils.h.

Enumeration Type Documentation

enum cvn::FinalState

Enumerator
kNumu0tr0sh	Numu CC - no track no shower.
kNumu0tr1sh	Numu CC - no track 1 shower.
kNumu0tr2sh	Numu CC - no track 2 shower.
kNumu0trMsh	Numu CC - no track 3+ shower.
kNumu1tr0sh	Numu CC - 1 track no shower.
kNumu1tr1sh	Numu CC - 1 track 1 shower.
kNumu1tr2sh	Numu CC - 1 track 2 shower.
kNumu1trMsh	Numu CC - 1 track 3+ shower.
kNumu2tr0sh	Numu CC - 2 track no shower.
kNumu2tr1sh	Numu CC - 2 track 1 shower.
kNumu2tr2sh	Numu CC - 2 track 2 shower.
kNumu2trMsh	Numu CC - 2 track 3+ shower.
kNumuMtr0sh	Numu CC - 3+ track no shower.
kNumuMtr1sh	Numu CC - 3+ track 1 shower.
kNumuMtr2sh	Numu CC - 3+ track 2 shower.
kNumuMtrMsh	Numu CC - 3+ track 3+ shower.
kNue0tr0sh	Nue CC - no track no shower.
kNue0tr1sh	Nue CC - no track 1 shower.
kNue0tr2sh	Nue CC - no track 2 shower.
kNue0trMsh	Nue CC - no track 3+ shower.
kNue1tr0sh	Nue CC - 1 track no shower.
kNue1tr1sh	Nue CC - 1 track 1 shower.
kNue1tr2sh	Nue CC - 1 track 2 shower.
kNue1trMsh	Nue CC - 1 track 3+ shower.
kNue2tr0sh	Nue CC - 2 track no shower.
kNue2tr1sh	Nue CC - 2 track 1 shower.
kNue2tr2sh	Nue CC - 2 track 2 shower.
kNue2trMsh	Nue CC - 2 track 3+ shower.
kNueMtr0sh	Nue CC - 3+ track no shower.
kNueMtr1sh	Nue CC - 3+ track 1 shower.
kNueMtr2sh	Nue CC - 3+ track 2 shower.
kNueMtrMsh	Nue CC - 3+ track 3+ shower.
kNC0tr0sh	NC CC - no track no shower.
kNC0tr1sh	NC CC - no track 1 shower.
kNC0tr2sh	NC CC - no track 2 shower.
kNC0trMsh	NC CC - no track 3+ shower.
kNC1tr0sh	NC CC - 1 track no shower.
kNC1tr1sh	NC CC - 1 track 1 shower.
kNC1tr2sh	NC CC - 1 track 2 shower.
kNC1trMsh	NC CC - 1 track 3+ shower.
kNC2tr0sh	NC CC - 2 track no shower.
kNC2tr1sh	NC CC - 2 track 1 shower.
kNC2tr2sh	NC CC - 2 track 2 shower.
kNC2trMsh	NC CC - 2 track 3+ shower.
kNCMtr0sh	NC CC - 3+ track no shower.
kNCMtr1sh	NC CC - 3+ track 1 shower.
kNCMtr2sh	NC CC - 3+ track 2 shower.
kNCMtrMsh	NC CC - 3+ track 3+ shower.
kCosmicFS	Cosmic ray background.
kOtherFS	Something else. Tau? Hopefully we don't use this.
kNFStType	Number of interaction types, used like a vector size.

Definition at line 15 of file FinalStateType.h.

  {
    kNumu0tr0sh,          ///< Numu CC - no track no shower
    kNumu0tr1sh,          ///< Numu CC - no track  1 shower
    kNumu0tr2sh,          ///< Numu CC - no track  2 shower
    kNumu0trMsh,          ///< Numu CC - no track 3+ shower
    kNumu1tr0sh,          ///< Numu CC -  1 track no shower
    kNumu1tr1sh,          ///< Numu CC -  1 track  1 shower
    kNumu1tr2sh,          ///< Numu CC -  1 track  2 shower
    kNumu1trMsh,          ///< Numu CC -  1 track 3+ shower
    kNumu2tr0sh,          ///< Numu CC -  2 track no shower
    kNumu2tr1sh,          ///< Numu CC -  2 track  1 shower
    kNumu2tr2sh,          ///< Numu CC -  2 track  2 shower
    kNumu2trMsh,          ///< Numu CC -  2 track 3+ shower
    kNumuMtr0sh,          ///< Numu CC - 3+ track no shower
    kNumuMtr1sh,          ///< Numu CC - 3+ track  1 shower
    kNumuMtr2sh,          ///< Numu CC - 3+ track  2 shower
    kNumuMtrMsh,          ///< Numu CC - 3+ track 3+ shower
    kNue0tr0sh,           ///< Nue CC - no track no shower
    kNue0tr1sh,           ///< Nue CC - no track  1 shower
    kNue0tr2sh,           ///< Nue CC - no track  2 shower
    kNue0trMsh,           ///< Nue CC - no track 3+ shower
    kNue1tr0sh,           ///< Nue CC -  1 track no shower
    kNue1tr1sh,           ///< Nue CC -  1 track  1 shower
    kNue1tr2sh,           ///< Nue CC -  1 track  2 shower
    kNue1trMsh,           ///< Nue CC -  1 track 3+ shower
    kNue2tr0sh,           ///< Nue CC -  2 track no shower
    kNue2tr1sh,           ///< Nue CC -  2 track  1 shower
    kNue2tr2sh,           ///< Nue CC -  2 track  2 shower
    kNue2trMsh,           ///< Nue CC -  2 track 3+ shower
    kNueMtr0sh,           ///< Nue CC - 3+ track no shower
    kNueMtr1sh,           ///< Nue CC - 3+ track  1 shower
    kNueMtr2sh,           ///< Nue CC - 3+ track  2 shower
    kNueMtrMsh,           ///< Nue CC - 3+ track 3+ shower
    kNC0tr0sh,           ///< NC CC - no track no shower
    kNC0tr1sh,           ///< NC CC - no track  1 shower
    kNC0tr2sh,           ///< NC CC - no track  2 shower
    kNC0trMsh,           ///< NC CC - no track 3+ shower
    kNC1tr0sh,           ///< NC CC -  1 track no shower
    kNC1tr1sh,           ///< NC CC -  1 track  1 shower
    kNC1tr2sh,           ///< NC CC -  1 track  2 shower
    kNC1trMsh,           ///< NC CC -  1 track 3+ shower
    kNC2tr0sh,           ///< NC CC -  2 track no shower
    kNC2tr1sh,           ///< NC CC -  2 track  1 shower
    kNC2tr2sh,           ///< NC CC -  2 track  2 shower
    kNC2trMsh,           ///< NC CC -  2 track 3+ shower
    kNCMtr0sh,           ///< NC CC - 3+ track no shower
    kNCMtr1sh,           ///< NC CC - 3+ track  1 shower
    kNCMtr2sh,           ///< NC CC - 3+ track  2 shower
    kNCMtrMsh,           ///< NC CC - 3+ track 3+ shower
    kCosmicFS,           ///< Cosmic ray background
    kOtherFS,            ///< Something else.  Tau?  Hopefully we don't use this
    kNFStType            ///< Number of interaction types, used like a vector size
  } FinalStateType;

enum cvn::FinalStateProngs

Enumerator
kNumu1prong	Numu CC - 1 prong.
kNumu2prong	Numu CC - 2 prong.
kNumu3prong	Numu CC - 3 prong.
kNumuMprong	Numu CC - 4+ prong.
kNue1prong	Nue CC - 1 prong.
kNue2prong	Nue CC - 2 prong.
kNue3prong	Nue CC - 3 prong.
kNueMprong	Nue CC - 4+ prong.
kNC1prong	NC - 1 prong.
kNC2prong	NC - 2 prong.
kNC3prong	NC - 3 prong.
kNCMprong	NC - 4+ prong.
kCosmicFSp	Cosmic ray background.
kOtherFSp	Something else. Tau? Hopefully we don't use this.
kNFSpType	Number of interaction types, used like a vector size.

Definition at line 15 of file FinalStateProngsType.h.

  {
    kNumu1prong,          ///< Numu CC -  1 prong
    kNumu2prong,          ///< Numu CC -  2 prong
    kNumu3prong,          ///< Numu CC -  3 prong
    kNumuMprong,          ///< Numu CC - 4+ prong
    kNue1prong,           ///< Nue CC -  1 prong
    kNue2prong,           ///< Nue CC -  2 prong
    kNue3prong,           ///< Nue CC -  3 prong
    kNueMprong,           ///< Nue CC - 4+ prong
    kNC1prong,            ///< NC  -  1 prong
    kNC2prong,            ///< NC  -  2 prong
    kNC3prong,            ///< NC  -  3 prong
    kNCMprong,            ///< NC  - 4+ prong
    kCosmicFSp,         ///< Cosmic ray background
    kOtherFSp,          ///< Something else.  Tau?  Hopefully we don't use this
    kNFSpType          ///< Number of interaction types, used like a vector size
  } FinalStateProngsType;

enum cvn::HType

Enumerator
kEmptyHit
kElectronHit
kMuonHit
kProtonHit
kNeutronHit
kPionHit
kPiZeroHit
kGammaHit
kMostlyHA
kMostlyEM
kMostlyMU
kUnknownHit

Definition at line 14 of file HitType.h.

  {
    kEmptyHit,
    kElectronHit,
    kMuonHit,
    kProtonHit,
    kNeutronHit,
    kPionHit,
    kPiZeroHit,
    kGammaHit,
    kMostlyHA,
    kMostlyEM,
    kMostlyMU,
    kUnknownHit
  } HitType;

enum cvn::Interaction

Enumerator
kNumuQE	Numu CC QE interaction.
kNumuRes	Numu CC Resonant interaction.
kNumuDIS	Numu CC DIS interaction.
kNumuOther	Numu CC, other than above.
kNueQE	Nue CC QE interaction.
kNueRes	Nue CC Resonant interaction.
kNueDIS	Nue CC DIS interaction.
kNueOther	Nue CC, other than above.
kNutauQE	Nutau CC QE interaction.
kNutauRes	Nutau CC Resonant interaction.
kNutauDIS	Nutau CC DIS interaction.
kNutauOther	Nutau CC, other than above.
kNuElectronElastic	NC Nu On E Scattering.
kNC	NC interaction.
kCosmic	Cosmic ray background.
kOther	Something else. Tau? Hopefully we don't use this.
kNIntType	Number of interaction types, used like a vector size.

Definition at line 15 of file InteractionType.h.

  {
    kNumuQE,           ///< Numu CC QE interaction
    kNumuRes,          ///< Numu CC Resonant interaction
    kNumuDIS,          ///< Numu CC DIS interaction
    kNumuOther,        ///< Numu CC, other than above
    kNueQE,            ///< Nue CC QE interaction
    kNueRes,           ///< Nue CC Resonant interaction
    kNueDIS,           ///< Nue CC DIS interaction
    kNueOther,         ///< Nue CC, other than above
    kNutauQE,          ///< Nutau CC QE interaction
    kNutauRes,         ///< Nutau CC Resonant interaction
    kNutauDIS,         ///< Nutau CC DIS interaction
    kNutauOther,       ///< Nutau CC, other than above
    kNuElectronElastic,///< NC Nu On E Scattering
    kNC,               ///< NC interaction
    kCosmic,           ///< Cosmic ray background
    kOther,            ///< Something else.  Tau?  Hopefully we don't use this
    kNIntType          ///< Number of interaction types, used like a vector size
  } InteractionType;

enum cvn::ParentParticle

Enumerator
kPiNuCCPP	Pi Parent Particle.
kKaNuCCPP	K Parent Particle.
kMuNuCCPP	mu minus Parent Particle
kPiNuNCPP	Pi Parent Particle.
kKaNuNCPP	K Parent Particle.
kMuNuNCPP	mu minus Parent Particle
kCosmicPP	Cosmic ray background.
kOtherPP	Something else.
kUnknownPP	Catch all term, shouldn't be used.
kNPPtType	Number of Parent Particle types, used like a vector size.

Definition at line 15 of file ParentParticleType.h.

                             {
    kPiNuCCPP,              ///< Pi Parent Particle
    kKaNuCCPP,              ///< K Parent Particle
    kMuNuCCPP,              ///< mu minus Parent Particle
    kPiNuNCPP,              ///< Pi Parent Particle
    kKaNuNCPP,              ///< K Parent Particle
    kMuNuNCPP,              ///< mu minus Parent Particle
    kCosmicPP,              ///< Cosmic ray background
    kOtherPP,               ///< Something else
    kUnknownPP,             ///< Catch all term, shouldn't be used
    kNPPtType              ///< Number of Parent Particle types, used like a vector size
  } ParentParticleType;

enum cvn::Particles

Enumerator
nE_e_1_mu_0_pi0_0_pi_0_n_0_p_1
nM_e_0_mu_1_pi0_0_pi_0_n_0_p_1
nM_e_0_mu_1_pi0_0_pi_1_n_0_p_1
nU_e_0_mu_0_pi0_1_pi_0_n_0_p_1
nE_e_1_mu_0_pi0_0_pi_1_n_0_p_1
nE_e_1_mu_0_pi0_0_pi_0_n_2_p_2
nU_e_0_mu_0_pi0_0_pi_1_n_0_p_1
nE_e_1_mu_0_pi0_0_pi_1_n_1_p_0
nM_e_0_mu_1_pi0_0_pi_1_n_1_p_0
nU_e_0_mu_0_pi0_0_pi_1_n_1_p_0
nU_e_0_mu_0_pi0_1_pi_0_n_1_p_0
nM_e_0_mu_1_pi0_0_pi_0_n_2_p_2
nU_e_0_mu_0_pi0_1_pi_0_n_1_p_1
nM_e_0_mu_1_pi0_0_pi_0_n_1_p_1
nE_e_1_mu_0_pi0_0_pi_0_n_1_p_1
nU_e_0_mu_0_pi0_1_pi_0_n_2_p_2
nE_e_1_mu_0_pi0_0_pi_0_n_1_p_2
nE_e_1_mu_0_pi0_0_pi_1_n_2_p_2
nM_e_0_mu_1_pi0_0_pi_1_n_1_p_1
nE_e_1_mu_0_pi0_0_pi_0_n_0_p_2
nU_e_0_mu_0_pi0_0_pi_2_n_0_p_1
nM_e_0_mu_1_pi0_0_pi_0_n_1_p_2
nU_e_0_mu_0_pi0_0_pi_0_n_2_p_2
nM_e_0_mu_1_pi0_0_pi_0_n_0_p_2
nE_e_1_mu_0_pi0_1_pi_0_n_0_p_1
nE_e_1_mu_0_pi0_0_pi_1_n_1_p_1
nU_e_0_mu_0_pi0_1_pi_1_n_1_p_0
nU_e_0_mu_0_pi0_0_pi_1_n_1_p_1
nM_e_0_mu_1_pi0_1_pi_0_n_0_p_1
nU_e_0_mu_0_pi0_0_pi_2_n_1_p_0
nU_e_0_mu_0_pi0_0_pi_1_n_2_p_2
nE_e_1_mu_0_pi0_1_pi_1_n_0_p_1
nU_e_0_mu_0_pi0_1_pi_2_n_0_p_1
nM_e_0_mu_1_pi0_0_pi_1_n_2_p_2
nE_e_1_mu_0_pi0_0_pi_2_n_0_p_1
nU_e_0_mu_0_pi0_1_pi_2_n_2_p_2
nU_e_0_mu_0_pi0_1_pi_2_n_1_p_0
nU_e_0_mu_0_pi0_1_pi_1_n_0_p_1
nU_e_0_mu_0_pi0_2_pi_2_n_0_p_1
nU_e_0_mu_0_pi0_0_pi_0_n_0_p_1
nU_e_0_mu_0_pi0_0_pi_2_n_2_p_2
nU_e_0_mu_0_pi0_0_pi_2_n_1_p_1
nM_e_0_mu_1_pi0_0_pi_1_n_0_p_2
nU_e_0_mu_0_pi0_1_pi_0_n_0_p_0
nE_e_1_mu_0_pi0_0_pi_1_n_0_p_2
nU_e_0_mu_0_pi0_0_pi_1_n_0_p_2
nT_e_0_mu_0_pi0_0_pi_2_n_0_p_1
nU_e_0_mu_0_pi0_1_pi_0_n_2_p_0
nU_e_0_mu_0_pi0_0_pi_1_n_2_p_0
nM_e_0_mu_1_pi0_1_pi_1_n_0_p_1
nE_e_1_mu_0_pi0_1_pi_2_n_0_p_1
nU_e_0_mu_0_pi0_2_pi_2_n_1_p_0
nM_e_0_mu_1_pi0_0_pi_1_n_2_p_1
nU_e_0_mu_0_pi0_1_pi_0_n_0_p_2
nU_e_0_mu_0_pi0_0_pi_1_n_2_p_1
nM_e_0_mu_1_pi0_1_pi_2_n_0_p_1
nM_e_0_mu_1_pi0_0_pi_2_n_1_p_0
nE_e_1_mu_0_pi0_0_pi_2_n_1_p_0
nU_e_0_mu_0_pi0_0_pi_0_n_1_p_0
nM_e_0_mu_1_pi0_0_pi_2_n_0_p_1
nM_e_0_mu_1_pi0_0_pi_1_n_1_p_2
nU_e_0_mu_0_pi0_2_pi_0_n_0_p_1
nU_e_0_mu_0_pi0_1_pi_1_n_2_p_2
nU_e_0_mu_0_pi0_1_pi_0_n_2_p_1
nU_e_0_mu_0_pi0_0_pi_1_n_1_p_2
nE_e_1_mu_0_pi0_1_pi_2_n_1_p_0
nE_e_1_mu_0_pi0_0_pi_0_n_2_p_1
nT_e_0_mu_0_pi0_0_pi_2_n_1_p_0
nM_e_0_mu_1_pi0_1_pi_0_n_0_p_2
nE_e_1_mu_0_pi0_0_pi_1_n_2_p_1
nU_e_0_mu_0_pi0_2_pi_2_n_1_p_1
nM_e_0_mu_1_pi0_1_pi_1_n_1_p_0
nM_e_0_mu_1_pi0_0_pi_2_n_1_p_1
nU_e_0_mu_0_pi0_0_pi_0_n_2_p_1
nM_e_0_mu_1_pi0_1_pi_0_n_2_p_2
nE_e_1_mu_0_pi0_0_pi_2_n_1_p_1
nU_e_0_mu_0_pi0_2_pi_2_n_2_p_2
nT_e_0_mu_0_pi0_2_pi_2_n_0_p_1
nM_e_0_mu_1_pi0_1_pi_1_n_1_p_1
nU_e_0_mu_0_pi0_1_pi_2_n_1_p_1
nU_e_0_mu_0_pi0_0_pi_2_n_2_p_1
nU_e_0_mu_0_pi0_0_pi_2_n_0_p_2
nM_e_0_mu_1_pi0_0_pi_0_n_2_p_1
nE_e_1_mu_0_pi0_1_pi_0_n_1_p_1
nM_e_0_mu_1_pi0_1_pi_0_n_1_p_1
nE_e_1_mu_0_pi0_1_pi_1_n_1_p_0
nE_e_1_mu_0_pi0_0_pi_2_n_2_p_2
nU_e_0_mu_0_pi0_1_pi_1_n_1_p_1
nE_e_1_mu_0_pi0_2_pi_2_n_2_p_2
nE_e_1_mu_0_pi0_2_pi_2_n_0_p_1
nE_e_1_mu_0_pi0_1_pi_1_n_1_p_1
nU_e_0_mu_0_pi0_1_pi_1_n_2_p_0
nU_e_0_mu_0_pi0_1_pi_1_n_1_p_2
nU_e_0_mu_0_pi0_0_pi_2_n_1_p_2
nU_e_0_mu_0_pi0_0_pi_0_n_1_p_2
nT_e_0_mu_0_pi0_2_pi_2_n_2_p_2
nT_e_0_mu_0_pi0_2_pi_2_n_1_p_0
nT_e_0_mu_0_pi0_1_pi_2_n_0_p_1
nM_e_0_mu_1_pi0_1_pi_2_n_1_p_0
nE_e_1_mu_0_pi0_1_pi_0_n_0_p_2
nM_e_0_mu_1_pi0_1_pi_2_n_1_p_1
nM_e_0_mu_1_pi0_0_pi_1_n_2_p_0
nU_e_0_mu_0_pi0_1_pi_2_n_2_p_0
nU_e_0_mu_0_pi0_0_pi_0_n_0_p_0
nT_e_0_mu_0_pi0_1_pi_2_n_1_p_0
nM_e_0_mu_1_pi0_2_pi_2_n_0_p_1
nE_e_1_mu_0_pi0_0_pi_2_n_1_p_2
nE_e_1_mu_0_pi0_0_pi_1_n_0_p_0
nU_e_0_mu_0_pi0_2_pi_1_n_2_p_2
nU_e_0_mu_0_pi0_2_pi_1_n_0_p_1
nT_e_0_mu_0_pi0_0_pi_2_n_1_p_1
nM_e_0_mu_1_pi0_0_pi_2_n_0_p_2
nE_e_1_mu_0_pi0_2_pi_2_n_1_p_0
nE_e_1_mu_0_pi0_0_pi_1_n_1_p_2
nU_e_0_mu_0_pi0_2_pi_0_n_2_p_2
nU_e_0_mu_0_pi0_0_pi_2_n_2_p_0
nU_e_0_mu_0_pi0_0_pi_0_n_1_p_1
nM_e_0_mu_1_pi0_1_pi_2_n_0_p_2
nM_e_0_mu_1_pi0_0_pi_1_n_0_p_0
nE_e_1_mu_0_pi0_1_pi_0_n_2_p_2
nU_e_0_mu_0_pi0_2_pi_2_n_0_p_2
nU_e_0_mu_0_pi0_1_pi_1_n_0_p_2
nM_e_0_mu_1_pi0_1_pi_2_n_2_p_0
nM_e_0_mu_1_pi0_0_pi_2_n_2_p_2
nE_e_1_mu_0_pi0_0_pi_1_n_2_p_0
nU_e_0_mu_0_pi0_2_pi_0_n_1_p_1
nM_e_0_mu_1_pi0_2_pi_2_n_2_p_2
nE_e_1_mu_0_pi0_1_pi_2_n_2_p_2
nE_e_1_mu_0_pi0_1_pi_2_n_2_p_1
nE_e_1_mu_0_pi0_0_pi_0_n_1_p_0
nU_e_0_mu_0_pi0_2_pi_2_n_2_p_0
nU_e_0_mu_0_pi0_2_pi_0_n_1_p_0
nU_e_0_mu_0_pi0_1_pi_2_n_2_p_1
nU_e_0_mu_0_pi0_1_pi_2_n_1_p_2
nU_e_0_mu_0_pi0_1_pi_1_n_2_p_1
nU_e_0_mu_0_pi0_0_pi_0_n_2_p_0
nT_e_0_mu_0_pi0_0_pi_1_n_0_p_1
nM_e_0_mu_1_pi0_1_pi_2_n_2_p_2
nM_e_0_mu_1_pi0_0_pi_0_n_1_p_0
nE_e_1_mu_0_pi0_2_pi_1_n_0_p_1
nE_e_1_mu_0_pi0_1_pi_1_n_2_p_2
nE_e_1_mu_0_pi0_1_pi_0_n_1_p_2
nU_e_0_mu_0_pi0_1_pi_0_n_1_p_2
nU_e_0_mu_0_pi0_0_pi_2_n_0_p_0
nT_e_0_mu_0_pi0_1_pi_2_n_2_p_2
nT_e_0_mu_0_pi0_0_pi_2_n_2_p_2
nM_e_0_mu_1_pi0_2_pi_2_n_1_p_0
nM_e_0_mu_1_pi0_1_pi_1_n_0_p_2
nM_e_0_mu_1_pi0_1_pi_0_n_1_p_2
nE_e_1_mu_0_pi0_1_pi_2_n_1_p_2
nE_e_1_mu_0_pi0_1_pi_1_n_0_p_2
nU_e_0_mu_0_pi0_2_pi_1_n_1_p_0
nT_e_0_mu_1_pi0_0_pi_0_n_0_p_1
nM_e_0_mu_1_pi0_2_pi_2_n_1_p_1
nM_e_0_mu_1_pi0_2_pi_1_n_0_p_1
nM_e_0_mu_1_pi0_1_pi_1_n_2_p_2
nM_e_0_mu_1_pi0_1_pi_1_n_1_p_2
nM_e_0_mu_1_pi0_0_pi_2_n_2_p_1
nE_e_1_mu_0_pi0_2_pi_1_n_2_p_2
nE_e_1_mu_0_pi0_0_pi_2_n_0_p_2
nE_e_1_mu_0_pi0_0_pi_0_n_0_p_0
nU_e_0_mu_0_pi0_2_pi_2_n_2_p_1
nU_e_0_mu_0_pi0_2_pi_2_n_1_p_2
nU_e_0_mu_0_pi0_2_pi_1_n_2_p_1
nU_e_0_mu_0_pi0_2_pi_0_n_1_p_2
nU_e_0_mu_0_pi0_1_pi_2_n_0_p_2
nU_e_0_mu_0_pi0_0_pi_0_n_0_p_2
nM_e_0_mu_1_pi0_2_pi_0_n_0_p_1
nE_e_1_mu_0_pi0_2_pi_2_n_1_p_1
nE_e_1_mu_0_pi0_2_pi_0_n_0_p_1
nE_e_1_mu_0_pi0_1_pi_1_n_2_p_1
nT_e_0_mu_1_pi0_0_pi_1_n_0_p_1
nM_e_0_mu_1_pi0_2_pi_0_n_0_p_2
nM_e_0_mu_1_pi0_0_pi_2_n_2_p_0
nM_e_0_mu_1_pi0_0_pi_2_n_1_p_2
nE_e_1_mu_0_pi0_2_pi_2_n_2_p_1
nE_e_1_mu_0_pi0_2_pi_2_n_0_p_2
nE_e_1_mu_0_pi0_1_pi_2_n_0_p_2
nE_e_1_mu_0_pi0_1_pi_1_n_1_p_2
nT_e_0_mu_0_pi0_1_pi_1_n_0_p_1
nT_e_0_mu_0_pi0_0_pi_2_n_0_p_2
nM_e_0_mu_1_pi0_1_pi_2_n_2_p_1
nE_e_1_mu_0_pi0_2_pi_0_n_0_p_2
nE_e_1_mu_0_pi0_1_pi_0_n_2_p_1
nE_e_1_mu_0_pi0_0_pi_2_n_2_p_0
nE_e_1_mu_0_pi0_0_pi_0_n_2_p_0
nU_e_0_mu_0_pi0_2_pi_1_n_1_p_1
nT_e_0_mu_1_pi0_0_pi_1_n_1_p_0
nT_e_0_mu_0_pi0_0_pi_2_n_2_p_0
nM_e_0_mu_1_pi0_1_pi_2_n_1_p_2
nE_e_1_mu_0_pi0_2_pi_2_n_1_p_2
nE_e_1_mu_0_pi0_2_pi_1_n_1_p_0
nU_e_0_mu_0_pi0_2_pi_0_n_2_p_1
nU_e_0_mu_0_pi0_1_pi_2_n_0_p_0
nT_e_1_mu_0_pi0_0_pi_1_n_1_p_0
nT_e_1_mu_0_pi0_0_pi_1_n_0_p_1
nT_e_0_mu_0_pi0_2_pi_2_n_0_p_2
nT_e_0_mu_0_pi0_2_pi_1_n_0_p_1
nT_e_0_mu_0_pi0_1_pi_2_n_1_p_1
nM_e_0_mu_1_pi0_2_pi_2_n_0_p_2
nM_e_0_mu_1_pi0_2_pi_1_n_1_p_0
nE_e_1_mu_0_pi0_1_pi_2_n_1_p_1
nE_e_1_mu_0_pi0_0_pi_2_n_0_p_0
nU_e_0_mu_0_pi0_2_pi_1_n_2_p_0
nU_e_0_mu_0_pi0_2_pi_1_n_0_p_2
nU_e_0_mu_0_pi0_2_pi_0_n_2_p_0
nT_e_1_mu_0_pi0_1_pi_2_n_1_p_0
nT_e_1_mu_0_pi0_1_pi_2_n_0_p_1
nT_e_1_mu_0_pi0_0_pi_0_n_0_p_1
nT_e_0_mu_1_pi0_0_pi_2_n_2_p_2
nT_e_0_mu_1_pi0_0_pi_2_n_0_p_1
nT_e_0_mu_0_pi0_2_pi_2_n_1_p_1
nT_e_0_mu_0_pi0_1_pi_2_n_2_p_1
nT_e_0_mu_0_pi0_1_pi_2_n_1_p_2
nT_e_0_mu_0_pi0_0_pi_2_n_1_p_2
nT_e_0_mu_0_pi0_0_pi_2_n_0_p_0
nT_e_0_mu_0_pi0_0_pi_1_n_2_p_2
nM_e_0_mu_1_pi0_2_pi_1_n_0_p_2
nM_e_0_mu_1_pi0_1_pi_1_n_2_p_1
nM_e_0_mu_1_pi0_1_pi_1_n_2_p_0
nM_e_0_mu_1_pi0_1_pi_0_n_1_p_0
nM_e_0_mu_1_pi0_0_pi_2_n_0_p_0
nE_e_1_mu_0_pi0_2_pi_1_n_1_p_2
nE_e_1_mu_0_pi0_2_pi_1_n_1_p_1
nE_e_1_mu_0_pi0_2_pi_1_n_0_p_2
nE_e_1_mu_0_pi0_2_pi_0_n_1_p_1
nE_e_1_mu_0_pi0_0_pi_2_n_2_p_1
nU_e_0_mu_0_pi0_2_pi_2_n_0_p_0
nT_e_1_mu_0_pi0_0_pi_2_n_0_p_1
nT_e_1_mu_0_pi0_0_pi_0_n_1_p_1
nT_e_0_mu_0_pi0_2_pi_2_n_0_p_0
nM_e_0_mu_1_pi0_2_pi_2_n_2_p_0
nM_e_0_mu_1_pi0_2_pi_1_n_2_p_2
nM_e_0_mu_1_pi0_2_pi_1_n_1_p_2
nM_e_0_mu_1_pi0_2_pi_1_n_1_p_1
nM_e_0_mu_1_pi0_2_pi_0_n_2_p_2
nM_e_0_mu_1_pi0_2_pi_0_n_1_p_2
nE_e_1_mu_0_pi0_2_pi_2_n_2_p_0
nE_e_1_mu_0_pi0_2_pi_0_n_2_p_2
nE_e_1_mu_0_pi0_2_pi_0_n_2_p_1
nE_e_1_mu_0_pi0_2_pi_0_n_1_p_2
nE_e_1_mu_0_pi0_1_pi_2_n_2_p_0
nU_e_0_mu_0_pi0_2_pi_1_n_1_p_2
nU_e_0_mu_0_pi0_1_pi_1_n_0_p_0
nU_e_0_mu_0_pi0_0_pi_1_n_0_p_0
nT_e_1_mu_0_pi0_1_pi_2_n_1_p_2
nT_e_1_mu_0_pi0_0_pi_2_n_2_p_0
nT_e_1_mu_0_pi0_0_pi_1_n_2_p_2
nT_e_0_mu_1_pi0_2_pi_2_n_1_p_0
nT_e_0_mu_1_pi0_2_pi_2_n_0_p_1
nT_e_0_mu_1_pi0_1_pi_2_n_0_p_1
nT_e_0_mu_1_pi0_1_pi_0_n_0_p_1
nT_e_0_mu_1_pi0_0_pi_2_n_1_p_0
nT_e_0_mu_1_pi0_0_pi_2_n_0_p_2
nT_e_0_mu_1_pi0_0_pi_0_n_2_p_2
nT_e_0_mu_1_pi0_0_pi_0_n_0_p_2
nT_e_0_mu_0_pi0_2_pi_2_n_2_p_0
nT_e_0_mu_0_pi0_2_pi_2_n_1_p_2
nT_e_0_mu_0_pi0_2_pi_1_n_0_p_2
nT_e_0_mu_0_pi0_1_pi_2_n_0_p_2
nT_e_0_mu_0_pi0_1_pi_1_n_1_p_2
nT_e_0_mu_0_pi0_1_pi_1_n_1_p_1
nT_e_0_mu_0_pi0_1_pi_1_n_0_p_2
nT_e_0_mu_0_pi0_0_pi_1_n_1_p_2
nT_e_0_mu_0_pi0_0_pi_1_n_0_p_0
nM_e_0_mu_1_pi0_2_pi_2_n_2_p_1
nM_e_0_mu_1_pi0_2_pi_2_n_1_p_2
nM_e_0_mu_1_pi0_2_pi_0_n_1_p_1
nM_e_0_mu_1_pi0_1_pi_0_n_2_p_1
nM_e_0_mu_1_pi0_1_pi_0_n_0_p_0
nM_e_0_mu_1_pi0_0_pi_0_n_0_p_0
nE_e_1_mu_0_pi0_1_pi_2_n_0_p_0
nU_e_0_mu_0_pi0_2_pi_0_n_0_p_2
nT_e_1_mu_0_pi0_1_pi_0_n_0_p_1
nT_e_1_mu_0_pi0_0_pi_0_n_0_p_0
nT_e_0_mu_1_pi0_1_pi_2_n_2_p_2
nT_e_0_mu_1_pi0_1_pi_1_n_1_p_0
nT_e_0_mu_0_pi0_2_pi_2_n_2_p_1
nT_e_0_mu_0_pi0_1_pi_1_n_2_p_2
nT_e_0_mu_0_pi0_1_pi_1_n_1_p_0
nT_e_0_mu_0_pi0_0_pi_1_n_1_p_1
nT_e_0_mu_0_pi0_0_pi_1_n_0_p_2
nM_e_0_mu_1_pi0_2_pi_1_n_2_p_1
nM_e_0_mu_1_pi0_2_pi_1_n_2_p_0
nM_e_0_mu_1_pi0_1_pi_2_n_0_p_0
nM_e_0_mu_1_pi0_0_pi_0_n_2_p_0
nE_e_1_mu_0_pi0_2_pi_1_n_2_p_1
nE_e_1_mu_0_pi0_1_pi_1_n_2_p_0
nE_e_1_mu_0_pi0_1_pi_1_n_0_p_0
nE_e_1_mu_0_pi0_1_pi_0_n_0_p_0
nU_e_1_mu_0_pi0_0_pi_2_n_1_p_1
nU_e_0_mu_0_pi0_2_pi_1_n_0_p_0
nT_e_1_mu_0_pi0_2_pi_2_n_1_p_1
nT_e_1_mu_0_pi0_2_pi_2_n_1_p_0
nT_e_1_mu_0_pi0_2_pi_2_n_0_p_1
nT_e_1_mu_0_pi0_2_pi_1_n_1_p_1
nT_e_1_mu_0_pi0_1_pi_2_n_1_p_1
nT_e_1_mu_0_pi0_1_pi_2_n_0_p_2
nT_e_1_mu_0_pi0_1_pi_1_n_1_p_1
nT_e_1_mu_0_pi0_1_pi_1_n_0_p_1
nT_e_1_mu_0_pi0_0_pi_0_n_2_p_2
nT_e_0_mu_1_pi0_2_pi_2_n_2_p_2
nT_e_0_mu_1_pi0_2_pi_2_n_2_p_1
nT_e_0_mu_1_pi0_2_pi_0_n_1_p_2
nT_e_0_mu_1_pi0_1_pi_2_n_1_p_1
nT_e_0_mu_1_pi0_1_pi_2_n_1_p_0
nT_e_0_mu_1_pi0_1_pi_0_n_2_p_2
nT_e_0_mu_1_pi0_1_pi_0_n_2_p_1
nT_e_0_mu_1_pi0_0_pi_2_n_1_p_1
nT_e_0_mu_1_pi0_0_pi_1_n_2_p_2
nT_e_0_mu_1_pi0_0_pi_1_n_2_p_0
nT_e_0_mu_1_pi0_0_pi_1_n_1_p_2
nT_e_0_mu_1_pi0_0_pi_1_n_1_p_1
nT_e_0_mu_1_pi0_0_pi_0_n_1_p_1
nT_e_0_mu_0_pi0_2_pi_1_n_2_p_2
nT_e_0_mu_0_pi0_1_pi_2_n_2_p_0
nT_e_0_mu_0_pi0_1_pi_2_n_0_p_0
nT_e_0_mu_0_pi0_1_pi_1_n_2_p_0
nT_e_0_mu_0_pi0_0_pi_2_n_2_p_1
nT_e_0_mu_0_pi0_0_pi_1_n_2_p_1
nT_e_0_mu_0_pi0_0_pi_1_n_2_p_0
nT_e_0_mu_0_pi0_0_pi_1_n_1_p_0
nM_e_1_mu_1_pi0_0_pi_1_n_0_p_1
nM_e_0_mu_1_pi0_1_pi_1_n_0_p_0
nE_e_1_mu_0_pi0_2_pi_1_n_0_p_0
nU_e_1_mu_0_pi0_0_pi_0_n_0_p_0
nT_e_1_mu_0_pi0_2_pi_2_n_2_p_2
nT_e_1_mu_0_pi0_2_pi_2_n_0_p_2
nT_e_1_mu_0_pi0_2_pi_1_n_1_p_2
nT_e_1_mu_0_pi0_2_pi_1_n_1_p_0
nT_e_1_mu_0_pi0_2_pi_1_n_0_p_1
nT_e_1_mu_0_pi0_2_pi_1_n_0_p_0
nT_e_1_mu_0_pi0_2_pi_0_n_1_p_1
nT_e_1_mu_0_pi0_2_pi_0_n_0_p_1
nT_e_1_mu_0_pi0_1_pi_1_n_2_p_2
nT_e_1_mu_0_pi0_1_pi_1_n_1_p_2
nT_e_1_mu_0_pi0_1_pi_1_n_1_p_0
nT_e_1_mu_0_pi0_1_pi_0_n_2_p_2
nT_e_1_mu_0_pi0_1_pi_0_n_2_p_1
nT_e_1_mu_0_pi0_1_pi_0_n_1_p_1
nT_e_1_mu_0_pi0_1_pi_0_n_0_p_2
nT_e_1_mu_0_pi0_0_pi_2_n_2_p_1
nT_e_1_mu_0_pi0_0_pi_2_n_1_p_2
nT_e_1_mu_0_pi0_0_pi_2_n_1_p_0
nT_e_1_mu_0_pi0_0_pi_2_n_0_p_0
nT_e_1_mu_0_pi0_0_pi_1_n_2_p_1
nT_e_1_mu_0_pi0_0_pi_1_n_2_p_0
nT_e_1_mu_0_pi0_0_pi_1_n_1_p_1
nT_e_1_mu_0_pi0_0_pi_0_n_2_p_1
nT_e_1_mu_0_pi0_0_pi_0_n_1_p_2
nT_e_1_mu_0_pi0_0_pi_0_n_1_p_0
nT_e_1_mu_0_pi0_0_pi_0_n_0_p_2
nT_e_0_mu_1_pi0_2_pi_2_n_2_p_0
nT_e_0_mu_1_pi0_2_pi_2_n_1_p_2
nT_e_0_mu_1_pi0_2_pi_2_n_1_p_1
nT_e_0_mu_1_pi0_2_pi_2_n_0_p_2
nT_e_0_mu_1_pi0_2_pi_1_n_1_p_0
nT_e_0_mu_1_pi0_2_pi_1_n_0_p_1
nT_e_0_mu_1_pi0_2_pi_0_n_2_p_2
nT_e_0_mu_1_pi0_2_pi_0_n_0_p_1
nT_e_0_mu_1_pi0_1_pi_2_n_1_p_2
nT_e_0_mu_1_pi0_1_pi_2_n_0_p_2
nT_e_0_mu_1_pi0_1_pi_1_n_2_p_1
nT_e_0_mu_1_pi0_1_pi_1_n_0_p_1
nT_e_0_mu_1_pi0_0_pi_2_n_2_p_0
nT_e_0_mu_1_pi0_0_pi_2_n_1_p_2
nT_e_0_mu_1_pi0_0_pi_1_n_2_p_1
nT_e_0_mu_1_pi0_0_pi_1_n_0_p_2
nT_e_0_mu_1_pi0_0_pi_0_n_1_p_2
nT_e_0_mu_0_pi0_2_pi_1_n_1_p_2
nT_e_0_mu_0_pi0_2_pi_1_n_1_p_1
nT_e_0_mu_0_pi0_2_pi_1_n_0_p_0
nT_e_0_mu_0_pi0_2_pi_0_n_1_p_2
nT_e_0_mu_0_pi0_1_pi_1_n_2_p_1
nT_e_0_mu_0_pi0_1_pi_0_n_0_p_1
nT_e_0_mu_0_pi0_0_pi_0_n_0_p_1
nM_e_0_mu_1_pi0_2_pi_2_n_0_p_0
nM_e_0_mu_1_pi0_2_pi_0_n_2_p_1
nM_e_0_mu_1_pi0_2_pi_0_n_2_p_0
nM_e_0_mu_1_pi0_2_pi_0_n_0_p_0
nM_e_0_mu_1_pi0_1_pi_0_n_2_p_0
nE_e_1_mu_1_pi0_1_pi_0_n_1_p_0
nE_e_1_mu_0_pi0_2_pi_2_n_0_p_0
nE_e_1_mu_0_pi0_2_pi_1_n_2_p_0
nE_e_1_mu_0_pi0_2_pi_0_n_2_p_0
nE_e_1_mu_0_pi0_1_pi_0_n_2_p_0
nE_e_1_mu_0_pi0_1_pi_0_n_1_p_0
nM_Other
nE_Other
nT_Other
kCosmic_PT
kOther_PT	Whatever is left.
kNIntType_PT	Number of interaction types, used like a vector size.

Definition at line 15 of file ParticlesType.h.

  {
    //kSignal_ND,           /// At least one muon and one pion
    //kOther_ND,            ///< Background
    //kNIntType_ND          ///< Number of interaction types, used like a vector size

nE_e_1_mu_0_pi0_0_pi_0_n_0_p_1,
nM_e_0_mu_1_pi0_0_pi_0_n_0_p_1,
nM_e_0_mu_1_pi0_0_pi_1_n_0_p_1,
nU_e_0_mu_0_pi0_1_pi_0_n_0_p_1,
nE_e_1_mu_0_pi0_0_pi_1_n_0_p_1,
nE_e_1_mu_0_pi0_0_pi_0_n_2_p_2,
nU_e_0_mu_0_pi0_0_pi_1_n_0_p_1,
nE_e_1_mu_0_pi0_0_pi_1_n_1_p_0,
nM_e_0_mu_1_pi0_0_pi_1_n_1_p_0,
nU_e_0_mu_0_pi0_0_pi_1_n_1_p_0,
nU_e_0_mu_0_pi0_1_pi_0_n_1_p_0,
nM_e_0_mu_1_pi0_0_pi_0_n_2_p_2,
nU_e_0_mu_0_pi0_1_pi_0_n_1_p_1,
nM_e_0_mu_1_pi0_0_pi_0_n_1_p_1,
nE_e_1_mu_0_pi0_0_pi_0_n_1_p_1,
nU_e_0_mu_0_pi0_1_pi_0_n_2_p_2,
nE_e_1_mu_0_pi0_0_pi_0_n_1_p_2,
nE_e_1_mu_0_pi0_0_pi_1_n_2_p_2,
nM_e_0_mu_1_pi0_0_pi_1_n_1_p_1,
nE_e_1_mu_0_pi0_0_pi_0_n_0_p_2,
nU_e_0_mu_0_pi0_0_pi_2_n_0_p_1,
nM_e_0_mu_1_pi0_0_pi_0_n_1_p_2,
nU_e_0_mu_0_pi0_0_pi_0_n_2_p_2,
nM_e_0_mu_1_pi0_0_pi_0_n_0_p_2,
nE_e_1_mu_0_pi0_1_pi_0_n_0_p_1,
nE_e_1_mu_0_pi0_0_pi_1_n_1_p_1,
nU_e_0_mu_0_pi0_1_pi_1_n_1_p_0,
nU_e_0_mu_0_pi0_0_pi_1_n_1_p_1,
nM_e_0_mu_1_pi0_1_pi_0_n_0_p_1,
nU_e_0_mu_0_pi0_0_pi_2_n_1_p_0,
nU_e_0_mu_0_pi0_0_pi_1_n_2_p_2,
nE_e_1_mu_0_pi0_1_pi_1_n_0_p_1,
nU_e_0_mu_0_pi0_1_pi_2_n_0_p_1,
nM_e_0_mu_1_pi0_0_pi_1_n_2_p_2,
nE_e_1_mu_0_pi0_0_pi_2_n_0_p_1,
nU_e_0_mu_0_pi0_1_pi_2_n_2_p_2,
nU_e_0_mu_0_pi0_1_pi_2_n_1_p_0,
nU_e_0_mu_0_pi0_1_pi_1_n_0_p_1,
nU_e_0_mu_0_pi0_2_pi_2_n_0_p_1,
nU_e_0_mu_0_pi0_0_pi_0_n_0_p_1,
nU_e_0_mu_0_pi0_0_pi_2_n_2_p_2,
nU_e_0_mu_0_pi0_0_pi_2_n_1_p_1,
nM_e_0_mu_1_pi0_0_pi_1_n_0_p_2,
nU_e_0_mu_0_pi0_1_pi_0_n_0_p_0,
nE_e_1_mu_0_pi0_0_pi_1_n_0_p_2,
nU_e_0_mu_0_pi0_0_pi_1_n_0_p_2,
nT_e_0_mu_0_pi0_0_pi_2_n_0_p_1,
nU_e_0_mu_0_pi0_1_pi_0_n_2_p_0,
nU_e_0_mu_0_pi0_0_pi_1_n_2_p_0,
nM_e_0_mu_1_pi0_1_pi_1_n_0_p_1,
nE_e_1_mu_0_pi0_1_pi_2_n_0_p_1,
nU_e_0_mu_0_pi0_2_pi_2_n_1_p_0,
nM_e_0_mu_1_pi0_0_pi_1_n_2_p_1,
nU_e_0_mu_0_pi0_1_pi_0_n_0_p_2,
nU_e_0_mu_0_pi0_0_pi_1_n_2_p_1,
nM_e_0_mu_1_pi0_1_pi_2_n_0_p_1,
nM_e_0_mu_1_pi0_0_pi_2_n_1_p_0,
nE_e_1_mu_0_pi0_0_pi_2_n_1_p_0,
nU_e_0_mu_0_pi0_0_pi_0_n_1_p_0,
nM_e_0_mu_1_pi0_0_pi_2_n_0_p_1,
nM_e_0_mu_1_pi0_0_pi_1_n_1_p_2,
nU_e_0_mu_0_pi0_2_pi_0_n_0_p_1,
nU_e_0_mu_0_pi0_1_pi_1_n_2_p_2,
nU_e_0_mu_0_pi0_1_pi_0_n_2_p_1,
nU_e_0_mu_0_pi0_0_pi_1_n_1_p_2,
nE_e_1_mu_0_pi0_1_pi_2_n_1_p_0,
nE_e_1_mu_0_pi0_0_pi_0_n_2_p_1,
nT_e_0_mu_0_pi0_0_pi_2_n_1_p_0,
nM_e_0_mu_1_pi0_1_pi_0_n_0_p_2,
nE_e_1_mu_0_pi0_0_pi_1_n_2_p_1,
nU_e_0_mu_0_pi0_2_pi_2_n_1_p_1,
nM_e_0_mu_1_pi0_1_pi_1_n_1_p_0,
nM_e_0_mu_1_pi0_0_pi_2_n_1_p_1,
nU_e_0_mu_0_pi0_0_pi_0_n_2_p_1,
nM_e_0_mu_1_pi0_1_pi_0_n_2_p_2,
nE_e_1_mu_0_pi0_0_pi_2_n_1_p_1,
nU_e_0_mu_0_pi0_2_pi_2_n_2_p_2,
nT_e_0_mu_0_pi0_2_pi_2_n_0_p_1,
nM_e_0_mu_1_pi0_1_pi_1_n_1_p_1,
nU_e_0_mu_0_pi0_1_pi_2_n_1_p_1,
nU_e_0_mu_0_pi0_0_pi_2_n_2_p_1,
nU_e_0_mu_0_pi0_0_pi_2_n_0_p_2,
nM_e_0_mu_1_pi0_0_pi_0_n_2_p_1,
nE_e_1_mu_0_pi0_1_pi_0_n_1_p_1,
nM_e_0_mu_1_pi0_1_pi_0_n_1_p_1,
nE_e_1_mu_0_pi0_1_pi_1_n_1_p_0,
nE_e_1_mu_0_pi0_0_pi_2_n_2_p_2,
nU_e_0_mu_0_pi0_1_pi_1_n_1_p_1,
nE_e_1_mu_0_pi0_2_pi_2_n_2_p_2,
nE_e_1_mu_0_pi0_2_pi_2_n_0_p_1,
nE_e_1_mu_0_pi0_1_pi_1_n_1_p_1,
nU_e_0_mu_0_pi0_1_pi_1_n_2_p_0,
nU_e_0_mu_0_pi0_1_pi_1_n_1_p_2,
nU_e_0_mu_0_pi0_0_pi_2_n_1_p_2,
nU_e_0_mu_0_pi0_0_pi_0_n_1_p_2,
nT_e_0_mu_0_pi0_2_pi_2_n_2_p_2,
nT_e_0_mu_0_pi0_2_pi_2_n_1_p_0,
nT_e_0_mu_0_pi0_1_pi_2_n_0_p_1,
nM_e_0_mu_1_pi0_1_pi_2_n_1_p_0,
nE_e_1_mu_0_pi0_1_pi_0_n_0_p_2,
nM_e_0_mu_1_pi0_1_pi_2_n_1_p_1,
nM_e_0_mu_1_pi0_0_pi_1_n_2_p_0,
nU_e_0_mu_0_pi0_1_pi_2_n_2_p_0,
nU_e_0_mu_0_pi0_0_pi_0_n_0_p_0,
nT_e_0_mu_0_pi0_1_pi_2_n_1_p_0,
nM_e_0_mu_1_pi0_2_pi_2_n_0_p_1,
nE_e_1_mu_0_pi0_0_pi_2_n_1_p_2,
nE_e_1_mu_0_pi0_0_pi_1_n_0_p_0,
nU_e_0_mu_0_pi0_2_pi_1_n_2_p_2,
nU_e_0_mu_0_pi0_2_pi_1_n_0_p_1,
nT_e_0_mu_0_pi0_0_pi_2_n_1_p_1,
nM_e_0_mu_1_pi0_0_pi_2_n_0_p_2,
nE_e_1_mu_0_pi0_2_pi_2_n_1_p_0,
nE_e_1_mu_0_pi0_0_pi_1_n_1_p_2,
nU_e_0_mu_0_pi0_2_pi_0_n_2_p_2,
nU_e_0_mu_0_pi0_0_pi_2_n_2_p_0,
nU_e_0_mu_0_pi0_0_pi_0_n_1_p_1,
nM_e_0_mu_1_pi0_1_pi_2_n_0_p_2,
nM_e_0_mu_1_pi0_0_pi_1_n_0_p_0,
nE_e_1_mu_0_pi0_1_pi_0_n_2_p_2,
nU_e_0_mu_0_pi0_2_pi_2_n_0_p_2,
nU_e_0_mu_0_pi0_1_pi_1_n_0_p_2,
nM_e_0_mu_1_pi0_1_pi_2_n_2_p_0,
nM_e_0_mu_1_pi0_0_pi_2_n_2_p_2,
nE_e_1_mu_0_pi0_0_pi_1_n_2_p_0,
nU_e_0_mu_0_pi0_2_pi_0_n_1_p_1,
nM_e_0_mu_1_pi0_2_pi_2_n_2_p_2,
nE_e_1_mu_0_pi0_1_pi_2_n_2_p_2,
nE_e_1_mu_0_pi0_1_pi_2_n_2_p_1,
nE_e_1_mu_0_pi0_0_pi_0_n_1_p_0,
nU_e_0_mu_0_pi0_2_pi_2_n_2_p_0,
nU_e_0_mu_0_pi0_2_pi_0_n_1_p_0,
nU_e_0_mu_0_pi0_1_pi_2_n_2_p_1,
nU_e_0_mu_0_pi0_1_pi_2_n_1_p_2,
nU_e_0_mu_0_pi0_1_pi_1_n_2_p_1,
nU_e_0_mu_0_pi0_0_pi_0_n_2_p_0,
nT_e_0_mu_0_pi0_0_pi_1_n_0_p_1,
nM_e_0_mu_1_pi0_1_pi_2_n_2_p_2,
nM_e_0_mu_1_pi0_0_pi_0_n_1_p_0,
nE_e_1_mu_0_pi0_2_pi_1_n_0_p_1,
nE_e_1_mu_0_pi0_1_pi_1_n_2_p_2,
nE_e_1_mu_0_pi0_1_pi_0_n_1_p_2,
nU_e_0_mu_0_pi0_1_pi_0_n_1_p_2,
nU_e_0_mu_0_pi0_0_pi_2_n_0_p_0,
nT_e_0_mu_0_pi0_1_pi_2_n_2_p_2,
nT_e_0_mu_0_pi0_0_pi_2_n_2_p_2,
nM_e_0_mu_1_pi0_2_pi_2_n_1_p_0,
nM_e_0_mu_1_pi0_1_pi_1_n_0_p_2,
nM_e_0_mu_1_pi0_1_pi_0_n_1_p_2,
nE_e_1_mu_0_pi0_1_pi_2_n_1_p_2,
nE_e_1_mu_0_pi0_1_pi_1_n_0_p_2,
nU_e_0_mu_0_pi0_2_pi_1_n_1_p_0,
nT_e_0_mu_1_pi0_0_pi_0_n_0_p_1,
nM_e_0_mu_1_pi0_2_pi_2_n_1_p_1,
nM_e_0_mu_1_pi0_2_pi_1_n_0_p_1,
nM_e_0_mu_1_pi0_1_pi_1_n_2_p_2,
nM_e_0_mu_1_pi0_1_pi_1_n_1_p_2,
nM_e_0_mu_1_pi0_0_pi_2_n_2_p_1,
nE_e_1_mu_0_pi0_2_pi_1_n_2_p_2,
nE_e_1_mu_0_pi0_0_pi_2_n_0_p_2,
nE_e_1_mu_0_pi0_0_pi_0_n_0_p_0,
nU_e_0_mu_0_pi0_2_pi_2_n_2_p_1,
nU_e_0_mu_0_pi0_2_pi_2_n_1_p_2,
nU_e_0_mu_0_pi0_2_pi_1_n_2_p_1,
nU_e_0_mu_0_pi0_2_pi_0_n_1_p_2,
nU_e_0_mu_0_pi0_1_pi_2_n_0_p_2,
nU_e_0_mu_0_pi0_0_pi_0_n_0_p_2,
nM_e_0_mu_1_pi0_2_pi_0_n_0_p_1,
nE_e_1_mu_0_pi0_2_pi_2_n_1_p_1,
nE_e_1_mu_0_pi0_2_pi_0_n_0_p_1,
nE_e_1_mu_0_pi0_1_pi_1_n_2_p_1,
nT_e_0_mu_1_pi0_0_pi_1_n_0_p_1,
nM_e_0_mu_1_pi0_2_pi_0_n_0_p_2,
nM_e_0_mu_1_pi0_0_pi_2_n_2_p_0,
nM_e_0_mu_1_pi0_0_pi_2_n_1_p_2,
nE_e_1_mu_0_pi0_2_pi_2_n_2_p_1,
nE_e_1_mu_0_pi0_2_pi_2_n_0_p_2,
nE_e_1_mu_0_pi0_1_pi_2_n_0_p_2,
nE_e_1_mu_0_pi0_1_pi_1_n_1_p_2,
nT_e_0_mu_0_pi0_1_pi_1_n_0_p_1,
nT_e_0_mu_0_pi0_0_pi_2_n_0_p_2,
nM_e_0_mu_1_pi0_1_pi_2_n_2_p_1,
nE_e_1_mu_0_pi0_2_pi_0_n_0_p_2,
nE_e_1_mu_0_pi0_1_pi_0_n_2_p_1,
nE_e_1_mu_0_pi0_0_pi_2_n_2_p_0,
nE_e_1_mu_0_pi0_0_pi_0_n_2_p_0,
nU_e_0_mu_0_pi0_2_pi_1_n_1_p_1,
nT_e_0_mu_1_pi0_0_pi_1_n_1_p_0,
nT_e_0_mu_0_pi0_0_pi_2_n_2_p_0,
nM_e_0_mu_1_pi0_1_pi_2_n_1_p_2,
nE_e_1_mu_0_pi0_2_pi_2_n_1_p_2,
nE_e_1_mu_0_pi0_2_pi_1_n_1_p_0,
nU_e_0_mu_0_pi0_2_pi_0_n_2_p_1,
nU_e_0_mu_0_pi0_1_pi_2_n_0_p_0,
nT_e_1_mu_0_pi0_0_pi_1_n_1_p_0,
nT_e_1_mu_0_pi0_0_pi_1_n_0_p_1,
nT_e_0_mu_0_pi0_2_pi_2_n_0_p_2,
nT_e_0_mu_0_pi0_2_pi_1_n_0_p_1,
nT_e_0_mu_0_pi0_1_pi_2_n_1_p_1,
nM_e_0_mu_1_pi0_2_pi_2_n_0_p_2,
nM_e_0_mu_1_pi0_2_pi_1_n_1_p_0,
nE_e_1_mu_0_pi0_1_pi_2_n_1_p_1,
nE_e_1_mu_0_pi0_0_pi_2_n_0_p_0,
nU_e_0_mu_0_pi0_2_pi_1_n_2_p_0,
nU_e_0_mu_0_pi0_2_pi_1_n_0_p_2,
nU_e_0_mu_0_pi0_2_pi_0_n_2_p_0,
nT_e_1_mu_0_pi0_1_pi_2_n_1_p_0,
nT_e_1_mu_0_pi0_1_pi_2_n_0_p_1,
nT_e_1_mu_0_pi0_0_pi_0_n_0_p_1,
nT_e_0_mu_1_pi0_0_pi_2_n_2_p_2,
nT_e_0_mu_1_pi0_0_pi_2_n_0_p_1,
nT_e_0_mu_0_pi0_2_pi_2_n_1_p_1,
nT_e_0_mu_0_pi0_1_pi_2_n_2_p_1,
nT_e_0_mu_0_pi0_1_pi_2_n_1_p_2,
nT_e_0_mu_0_pi0_0_pi_2_n_1_p_2,
nT_e_0_mu_0_pi0_0_pi_2_n_0_p_0,
nT_e_0_mu_0_pi0_0_pi_1_n_2_p_2,
nM_e_0_mu_1_pi0_2_pi_1_n_0_p_2,
nM_e_0_mu_1_pi0_1_pi_1_n_2_p_1,
nM_e_0_mu_1_pi0_1_pi_1_n_2_p_0,
nM_e_0_mu_1_pi0_1_pi_0_n_1_p_0,
nM_e_0_mu_1_pi0_0_pi_2_n_0_p_0,
nE_e_1_mu_0_pi0_2_pi_1_n_1_p_2,
nE_e_1_mu_0_pi0_2_pi_1_n_1_p_1,
nE_e_1_mu_0_pi0_2_pi_1_n_0_p_2,
nE_e_1_mu_0_pi0_2_pi_0_n_1_p_1,
nE_e_1_mu_0_pi0_0_pi_2_n_2_p_1,
nU_e_0_mu_0_pi0_2_pi_2_n_0_p_0,
nT_e_1_mu_0_pi0_0_pi_2_n_0_p_1,
nT_e_1_mu_0_pi0_0_pi_0_n_1_p_1,
nT_e_0_mu_0_pi0_2_pi_2_n_0_p_0,
nM_e_0_mu_1_pi0_2_pi_2_n_2_p_0,
nM_e_0_mu_1_pi0_2_pi_1_n_2_p_2,
nM_e_0_mu_1_pi0_2_pi_1_n_1_p_2,
nM_e_0_mu_1_pi0_2_pi_1_n_1_p_1,
nM_e_0_mu_1_pi0_2_pi_0_n_2_p_2,
nM_e_0_mu_1_pi0_2_pi_0_n_1_p_2,
nE_e_1_mu_0_pi0_2_pi_2_n_2_p_0,
nE_e_1_mu_0_pi0_2_pi_0_n_2_p_2,
nE_e_1_mu_0_pi0_2_pi_0_n_2_p_1,
nE_e_1_mu_0_pi0_2_pi_0_n_1_p_2,
nE_e_1_mu_0_pi0_1_pi_2_n_2_p_0,
nU_e_0_mu_0_pi0_2_pi_1_n_1_p_2,
nU_e_0_mu_0_pi0_1_pi_1_n_0_p_0,
nU_e_0_mu_0_pi0_0_pi_1_n_0_p_0,
nT_e_1_mu_0_pi0_1_pi_2_n_1_p_2,
nT_e_1_mu_0_pi0_0_pi_2_n_2_p_0,
nT_e_1_mu_0_pi0_0_pi_1_n_2_p_2,
nT_e_0_mu_1_pi0_2_pi_2_n_1_p_0,
nT_e_0_mu_1_pi0_2_pi_2_n_0_p_1,
nT_e_0_mu_1_pi0_1_pi_2_n_0_p_1,
nT_e_0_mu_1_pi0_1_pi_0_n_0_p_1,
nT_e_0_mu_1_pi0_0_pi_2_n_1_p_0,
nT_e_0_mu_1_pi0_0_pi_2_n_0_p_2,
nT_e_0_mu_1_pi0_0_pi_0_n_2_p_2,
nT_e_0_mu_1_pi0_0_pi_0_n_0_p_2,
nT_e_0_mu_0_pi0_2_pi_2_n_2_p_0,
nT_e_0_mu_0_pi0_2_pi_2_n_1_p_2,
nT_e_0_mu_0_pi0_2_pi_1_n_0_p_2,
nT_e_0_mu_0_pi0_1_pi_2_n_0_p_2,
nT_e_0_mu_0_pi0_1_pi_1_n_1_p_2,
nT_e_0_mu_0_pi0_1_pi_1_n_1_p_1,
nT_e_0_mu_0_pi0_1_pi_1_n_0_p_2,
nT_e_0_mu_0_pi0_0_pi_1_n_1_p_2,
nT_e_0_mu_0_pi0_0_pi_1_n_0_p_0,
nM_e_0_mu_1_pi0_2_pi_2_n_2_p_1,
nM_e_0_mu_1_pi0_2_pi_2_n_1_p_2,
nM_e_0_mu_1_pi0_2_pi_0_n_1_p_1,
nM_e_0_mu_1_pi0_1_pi_0_n_2_p_1,
nM_e_0_mu_1_pi0_1_pi_0_n_0_p_0,
nM_e_0_mu_1_pi0_0_pi_0_n_0_p_0,
nE_e_1_mu_0_pi0_1_pi_2_n_0_p_0,
nU_e_0_mu_0_pi0_2_pi_0_n_0_p_2,
nT_e_1_mu_0_pi0_1_pi_0_n_0_p_1,
nT_e_1_mu_0_pi0_0_pi_0_n_0_p_0,
nT_e_0_mu_1_pi0_1_pi_2_n_2_p_2,
nT_e_0_mu_1_pi0_1_pi_1_n_1_p_0,
nT_e_0_mu_0_pi0_2_pi_2_n_2_p_1,
nT_e_0_mu_0_pi0_1_pi_1_n_2_p_2,
nT_e_0_mu_0_pi0_1_pi_1_n_1_p_0,
nT_e_0_mu_0_pi0_0_pi_1_n_1_p_1,
nT_e_0_mu_0_pi0_0_pi_1_n_0_p_2,
nM_e_0_mu_1_pi0_2_pi_1_n_2_p_1,
nM_e_0_mu_1_pi0_2_pi_1_n_2_p_0,
nM_e_0_mu_1_pi0_1_pi_2_n_0_p_0,
nM_e_0_mu_1_pi0_0_pi_0_n_2_p_0,
nE_e_1_mu_0_pi0_2_pi_1_n_2_p_1,
nE_e_1_mu_0_pi0_1_pi_1_n_2_p_0,
nE_e_1_mu_0_pi0_1_pi_1_n_0_p_0,
nE_e_1_mu_0_pi0_1_pi_0_n_0_p_0,
nU_e_1_mu_0_pi0_0_pi_2_n_1_p_1,
nU_e_0_mu_0_pi0_2_pi_1_n_0_p_0,
nT_e_1_mu_0_pi0_2_pi_2_n_1_p_1,
nT_e_1_mu_0_pi0_2_pi_2_n_1_p_0,
nT_e_1_mu_0_pi0_2_pi_2_n_0_p_1,
nT_e_1_mu_0_pi0_2_pi_1_n_1_p_1,
nT_e_1_mu_0_pi0_1_pi_2_n_1_p_1,
nT_e_1_mu_0_pi0_1_pi_2_n_0_p_2,
nT_e_1_mu_0_pi0_1_pi_1_n_1_p_1,
nT_e_1_mu_0_pi0_1_pi_1_n_0_p_1,
nT_e_1_mu_0_pi0_0_pi_0_n_2_p_2,
nT_e_0_mu_1_pi0_2_pi_2_n_2_p_2,
nT_e_0_mu_1_pi0_2_pi_2_n_2_p_1,
nT_e_0_mu_1_pi0_2_pi_0_n_1_p_2,
nT_e_0_mu_1_pi0_1_pi_2_n_1_p_1,
nT_e_0_mu_1_pi0_1_pi_2_n_1_p_0,
nT_e_0_mu_1_pi0_1_pi_0_n_2_p_2,
nT_e_0_mu_1_pi0_1_pi_0_n_2_p_1,
nT_e_0_mu_1_pi0_0_pi_2_n_1_p_1,
nT_e_0_mu_1_pi0_0_pi_1_n_2_p_2,
nT_e_0_mu_1_pi0_0_pi_1_n_2_p_0,
nT_e_0_mu_1_pi0_0_pi_1_n_1_p_2,
nT_e_0_mu_1_pi0_0_pi_1_n_1_p_1,
nT_e_0_mu_1_pi0_0_pi_0_n_1_p_1,
nT_e_0_mu_0_pi0_2_pi_1_n_2_p_2,
nT_e_0_mu_0_pi0_1_pi_2_n_2_p_0,
nT_e_0_mu_0_pi0_1_pi_2_n_0_p_0,
nT_e_0_mu_0_pi0_1_pi_1_n_2_p_0,
nT_e_0_mu_0_pi0_0_pi_2_n_2_p_1,
nT_e_0_mu_0_pi0_0_pi_1_n_2_p_1,
nT_e_0_mu_0_pi0_0_pi_1_n_2_p_0,
nT_e_0_mu_0_pi0_0_pi_1_n_1_p_0,
nM_e_1_mu_1_pi0_0_pi_1_n_0_p_1,
nM_e_0_mu_1_pi0_1_pi_1_n_0_p_0,
nE_e_1_mu_0_pi0_2_pi_1_n_0_p_0,
nU_e_1_mu_0_pi0_0_pi_0_n_0_p_0,
nT_e_1_mu_0_pi0_2_pi_2_n_2_p_2,
nT_e_1_mu_0_pi0_2_pi_2_n_0_p_2,
nT_e_1_mu_0_pi0_2_pi_1_n_1_p_2,
nT_e_1_mu_0_pi0_2_pi_1_n_1_p_0,
nT_e_1_mu_0_pi0_2_pi_1_n_0_p_1,
nT_e_1_mu_0_pi0_2_pi_1_n_0_p_0,
nT_e_1_mu_0_pi0_2_pi_0_n_1_p_1,
nT_e_1_mu_0_pi0_2_pi_0_n_0_p_1,
nT_e_1_mu_0_pi0_1_pi_1_n_2_p_2,
nT_e_1_mu_0_pi0_1_pi_1_n_1_p_2,
nT_e_1_mu_0_pi0_1_pi_1_n_1_p_0,
nT_e_1_mu_0_pi0_1_pi_0_n_2_p_2,
nT_e_1_mu_0_pi0_1_pi_0_n_2_p_1,
nT_e_1_mu_0_pi0_1_pi_0_n_1_p_1,
nT_e_1_mu_0_pi0_1_pi_0_n_0_p_2,
nT_e_1_mu_0_pi0_0_pi_2_n_2_p_1,
nT_e_1_mu_0_pi0_0_pi_2_n_1_p_2,
nT_e_1_mu_0_pi0_0_pi_2_n_1_p_0,
nT_e_1_mu_0_pi0_0_pi_2_n_0_p_0,
nT_e_1_mu_0_pi0_0_pi_1_n_2_p_1,
nT_e_1_mu_0_pi0_0_pi_1_n_2_p_0,
nT_e_1_mu_0_pi0_0_pi_1_n_1_p_1,
nT_e_1_mu_0_pi0_0_pi_0_n_2_p_1,
nT_e_1_mu_0_pi0_0_pi_0_n_1_p_2,
nT_e_1_mu_0_pi0_0_pi_0_n_1_p_0,
nT_e_1_mu_0_pi0_0_pi_0_n_0_p_2,
nT_e_0_mu_1_pi0_2_pi_2_n_2_p_0,
nT_e_0_mu_1_pi0_2_pi_2_n_1_p_2,
nT_e_0_mu_1_pi0_2_pi_2_n_1_p_1,
nT_e_0_mu_1_pi0_2_pi_2_n_0_p_2,
nT_e_0_mu_1_pi0_2_pi_1_n_1_p_0,
nT_e_0_mu_1_pi0_2_pi_1_n_0_p_1,
nT_e_0_mu_1_pi0_2_pi_0_n_2_p_2,
nT_e_0_mu_1_pi0_2_pi_0_n_0_p_1,
nT_e_0_mu_1_pi0_1_pi_2_n_1_p_2,
nT_e_0_mu_1_pi0_1_pi_2_n_0_p_2,
nT_e_0_mu_1_pi0_1_pi_1_n_2_p_1,
nT_e_0_mu_1_pi0_1_pi_1_n_0_p_1,
nT_e_0_mu_1_pi0_0_pi_2_n_2_p_0,
nT_e_0_mu_1_pi0_0_pi_2_n_1_p_2,
nT_e_0_mu_1_pi0_0_pi_1_n_2_p_1,
nT_e_0_mu_1_pi0_0_pi_1_n_0_p_2,
nT_e_0_mu_1_pi0_0_pi_0_n_1_p_2,
nT_e_0_mu_0_pi0_2_pi_1_n_1_p_2,
nT_e_0_mu_0_pi0_2_pi_1_n_1_p_1,
nT_e_0_mu_0_pi0_2_pi_1_n_0_p_0,
nT_e_0_mu_0_pi0_2_pi_0_n_1_p_2,
nT_e_0_mu_0_pi0_1_pi_1_n_2_p_1,
nT_e_0_mu_0_pi0_1_pi_0_n_0_p_1,
nT_e_0_mu_0_pi0_0_pi_0_n_0_p_1,
nM_e_0_mu_1_pi0_2_pi_2_n_0_p_0,
nM_e_0_mu_1_pi0_2_pi_0_n_2_p_1,
nM_e_0_mu_1_pi0_2_pi_0_n_2_p_0,
nM_e_0_mu_1_pi0_2_pi_0_n_0_p_0,
nM_e_0_mu_1_pi0_1_pi_0_n_2_p_0,
nE_e_1_mu_1_pi0_1_pi_0_n_1_p_0,
nE_e_1_mu_0_pi0_2_pi_2_n_0_p_0,
nE_e_1_mu_0_pi0_2_pi_1_n_2_p_0,
nE_e_1_mu_0_pi0_2_pi_0_n_2_p_0,
nE_e_1_mu_0_pi0_1_pi_0_n_2_p_0,
nE_e_1_mu_0_pi0_1_pi_0_n_1_p_0,
nM_Other,
nE_Other,
nT_Other,
kCosmic_PT,
kOther_PT,            ///< Whatever is left
kNIntType_PT          ///< Number of interaction types, used like a vector size
  } ParticlesType;

enum cvn::PType

Enumerator
kElectron
kMuon
kProton
kNeutron
kPion
kPiZero
kGamma
kOtherPDG
kUnknown
kEmpty
kEM
kHadron

Definition at line 18 of file ProngType.h.

  {
    kElectron,
    kMuon,
    kProton,
    kNeutron,
    kPion,
    kPiZero,
    kGamma,
    kOtherPDG,
    kUnknown,
    kEmpty,
    kEM,
    kHadron
  } ProngType;

enum cvn::SimpleNeutrino

Enumerator
kNumu	Numu CC interaction.
kNue	Nue CC QE interaction.
kNutau	Nutau CC QE interaction.
kNuNC	NC interaction.
kCos	Cosmics.
kNSimpleNeutrinoType	Number of types, used like a vector size.

Definition at line 16 of file SimpleNeutrinoType.h.

  {
    kNumu,                ///< Numu CC interaction
    kNue,                 ///< Nue CC QE interaction
    kNutau,               ///< Nutau CC QE interaction
    kNuNC,                ///< NC interaction
    kCos,                 ///< Cosmics
    kNSimpleNeutrinoType  ///< Number of types, used like a vector size
  } SimpleNeutrinoType;

Function Documentation

std::vector< std::pair< int, double > > cvn::Get2DProngEfficiencyByPDG	(	const rb::Cluster &	prong,
		const rb::Cluster &	slice
	)

Definition at line 1401 of file AssignLabels.cxx.

References rb::Cluster::AllCells(), dE, cheat::BackTracker::HitCollectionEfficiency(), cheat::BackTracker::HitsToTrackIDE(), kEmpty, kUnknown, geo::kXorY, make_pair(), rb::Cluster::NXCell(), rb::Cluster::NYCell(), simb::MCParticle::PdgCode(), and cheat::BackTracker::TrackIDToParticle().

Referenced by bpf::BPFCVNAna::analyze().

  {
    
    art::ServiceHandle<cheat::BackTracker> bt;
    
    std::vector< std::pair<int, double> > efficiencybypdg;
    
    // Check nhits > 0 and number of associated trackIDs > 0 for the prong
    ProngType pType3D = kUnknown;
    if( prong.NXCell() == 0 && 
        prong.NYCell() == 0) pType3D = kEmpty;
    
    std::vector<cheat::TrackIDE> ParticleTrackIDs;
    if( pType3D != kEmpty )
      ParticleTrackIDs  = bt->HitsToTrackIDE(prong.AllCells());
    
    if( ParticleTrackIDs.size()  < 1 ) pType3D = kEmpty;
    
    if( pType3D != kEmpty ){
      
      for (unsigned int iID=0; iID < ParticleTrackIDs.size(); ++iID){
        std::set<int> TrackID;
        std::map<int, double> purMap;
        double dE, tE;
        std::map<int, int> parents;
        
        int iPartID = ParticleTrackIDs[iID].trackID;
        TrackID.insert(iPartID);
        int iPartPDG = bt->TrackIDToParticle(iPartID)->PdgCode();
        
        double iPartEff = bt->HitCollectionEfficiency(TrackID, prong.AllCells(),slice.AllCells(),
                                                      geo::kXorY, &purMap,true, &dE, &tE);

        efficiencybypdg.push_back(std::make_pair(iPartPDG, iPartEff));
        
      }//iTrackID      
    }
    
    return efficiencybypdg;
  }// 2DProngEfficiencyByPDG

std::vector< std::pair< int, double > > cvn::Get2DProngEnergyByPDG

(

const rb::Cluster &

prong

)

Definition at line 1444 of file AssignLabels.cxx.

References rb::Cluster::AllCells(), simb::MCParticle::E(), cheat::BackTracker::HitsToTrackIDE(), kEmpty, kUnknown, make_pair(), rb::Cluster::NXCell(), rb::Cluster::NYCell(), simb::MCParticle::PdgCode(), and cheat::BackTracker::TrackIDToParticle().

Referenced by bpf::BPFCVNAna::analyze().

  {
    
    art::ServiceHandle<cheat::BackTracker> bt;
    
    std::vector< std::pair<int, double> > energybypdg;
    
    // Check nhits > 0 and number of associated trackIDs > 0 for the prong
    ProngType pType3D = kUnknown;
    if( prong.NXCell() == 0 && 
        prong.NYCell() == 0) pType3D = kEmpty;
    
    std::vector<cheat::TrackIDE> ParticleTrackIDs;
    if( pType3D != kEmpty )
      ParticleTrackIDs  = bt->HitsToTrackIDE(prong.AllCells());
    
    if( ParticleTrackIDs.size()  < 1 ) pType3D = kEmpty;
    
    if( pType3D != kEmpty ){
      
      for (unsigned int iID=0; iID < ParticleTrackIDs.size(); ++iID){
        std::set<int> TrackID;
        std::map<int, double> purMap;
        std::map<int, int> parents;
        
        int iPartID = ParticleTrackIDs[iID].trackID;
        TrackID.insert(iPartID);
        int  iPartPDG = bt->TrackIDToParticle(iPartID)->PdgCode();
        double iPartE = bt->TrackIDToParticle(iPartID)->E();
        
        energybypdg.push_back(std::make_pair(iPartPDG, iPartE));
        
      }//iTrackID      
    }
    
    return energybypdg;
  }// 2DProngEnergyByPDG  

std::vector< std::pair< int, double > > cvn::Get2DProngPurityByPDG

(

const rb::Cluster &

prong

)

Definition at line 1361 of file AssignLabels.cxx.

References rb::Cluster::AllCells(), cheat::BackTracker::HitCollectionPurity(), cheat::BackTracker::HitsToTrackIDE(), kEmpty, kUnknown, make_pair(), rb::Cluster::NXCell(), rb::Cluster::NYCell(), simb::MCParticle::PdgCode(), and cheat::BackTracker::TrackIDToParticle().

Referenced by bpf::BPFCVNAna::analyze().

  {
    
    art::ServiceHandle<cheat::BackTracker> bt;

    std::vector< std::pair<int, double> > puritybypdg;

    // Check nhits > 0 and number of associated trackIDs > 0 for the prong
    ProngType   pType3D = kUnknown;    
    if( prong.NXCell() == 0 && 
        prong.NYCell() == 0) pType3D = kEmpty;
    
    
    std::vector<cheat::TrackIDE> ParticleTrackIDs;    
    if( pType3D != kEmpty )
      ParticleTrackIDs  = bt->HitsToTrackIDE(prong.AllCells());
    
    if( ParticleTrackIDs.size()  < 1 ) pType3D = kEmpty;
    
    if( pType3D != kEmpty ){
      
      for (unsigned int iID=0; iID < ParticleTrackIDs.size(); ++iID){
        std::set<int> TrackID;
        std::map<int, double> purMap;
        std::map<int, int> parents;
        
        int iPartID = ParticleTrackIDs[iID].trackID;
        TrackID.insert(iPartID);
        int iPartPDG = bt->TrackIDToParticle(iPartID)->PdgCode();
        double iPartPur = bt->HitCollectionPurity(TrackID, prong.AllCells(),
                                                  &purMap, &parents,true);
        puritybypdg.push_back(std::make_pair(iPartPDG, iPartPur));
        
      }//iTrackID
    }
    
    return puritybypdg;
  }// 2DProngPurityByPDG

FinalStateProngsType cvn::GetFinalStateProngsType

(

const art::Ptr< simb::MCTruth >

truth

)

Definition at line 137 of file AssignLabels.cxx.

References GetMultiplicityMap(), MECModelEnuComparisons::i, iscc, kCosmicFSp, kNC1prong, kNC2prong, kNC3prong, kNCMprong, kNue1prong, kNue2prong, kNueMprong, kNumu1prong, kNumu2prong, kNumuMprong, and kOtherFSp.

Referenced by cvn::CVNEventDump::analyze(), cvn::CVNCosmicEventDump::analyze(), cvn::CVNEventProngDump::analyze(), and cvn::CVNAddTrainingData::produce().

  {

    int nupdg, truemode;
    bool isnuint, iscc;
    std::vector<int>  PDGs={11,13,2212,211,111,2112,0}; // e,mu,p,pi,pi0,n,other
    std::map<int,int> NPrimariesPDG = GetMultiplicityMap(truth,PDGs,&isnuint
                                                         ,&iscc,&nupdg,&truemode);


    // Count number of prongs = charged primaries + 2*(pi0 primaries)
    int Nprongs=0;
    for ( unsigned int i = 0; i < PDGs.size(); ++i){

      if ( PDGs[i] == 2112 ) continue;
      if ( PDGs[i] == 111 )      
        Nprongs = Nprongs+(2*NPrimariesPDG[PDGs[i]]);
      else 
        Nprongs = Nprongs+( NPrimariesPDG[PDGs[i]] );

    }


    if ( isnuint ){
      if( iscc ){
        if ( nupdg == 14 && Nprongs == 1 ) return kNumu1prong;
        if ( nupdg == 14 && Nprongs == 2 ) return kNumu2prong;
        if ( nupdg == 14 && Nprongs == 3 ) return kNumu2prong;
        if ( nupdg == 14 && Nprongs > 3  ) return kNumuMprong;
        if ( nupdg == 12 && Nprongs == 1 ) return kNue1prong;
        if ( nupdg == 12 && Nprongs == 2 ) return kNue2prong;
        if ( nupdg == 12 && Nprongs == 3 ) return kNue2prong;
        if ( nupdg == 12 && Nprongs > 3  ) return kNueMprong;
        if ( nupdg != 12 && nupdg != 14 ) return kOtherFSp;     
      }//cc
      else{
        if ( Nprongs == 1 ) return kNC1prong;
        if ( Nprongs == 2 ) return kNC2prong;
        if ( Nprongs == 3 ) return kNC3prong;
        if ( Nprongs > 3  ) return kNCMprong;
      }//nc
    }
    else return kCosmicFSp;
    
    return kOtherFSp;
  }

FinalStateType cvn::GetFinalStateType

(

const art::Ptr< simb::MCTruth >

truth

)

Definition at line 64 of file AssignLabels.cxx.

References GetMultiplicityMap(), iscc, kCosmicFS, kNC0tr0sh, kNC0tr1sh, kNC0tr2sh, kNC0trMsh, kNC1tr0sh, kNC1tr1sh, kNC1tr2sh, kNC1trMsh, kNC2tr0sh, kNC2tr1sh, kNC2tr2sh, kNC2trMsh, kNCMtr0sh, kNCMtr1sh, kNCMtr2sh, kNCMtrMsh, kNue0tr0sh, kNue0tr1sh, kNue0tr2sh, kNue0trMsh, kNue1tr0sh, kNue1tr1sh, kNue1tr2sh, kNue1trMsh, kNue2tr0sh, kNue2tr1sh, kNue2tr2sh, kNue2trMsh, kNueMtr0sh, kNueMtr1sh, kNueMtr2sh, kNueMtrMsh, kNumu0tr0sh, kNumu0tr1sh, kNumu0tr2sh, kNumu0trMsh, kNumu1tr0sh, kNumu1tr1sh, kNumu1tr2sh, kNumu1trMsh, kNumu2tr0sh, kNumu2tr1sh, kNumu2tr2sh, kNumu2trMsh, kNumuMtr0sh, kNumuMtr1sh, kNumuMtr2sh, kNumuMtrMsh, and kOtherFS.

Referenced by cvn::CVNEventDump::analyze(), cvn::CVNCosmicEventDump::analyze(), cvn::CVNEventProngDump::analyze(), and cvn::CVNAddTrainingData::produce().

  {

    int nupdg, truemode;
    bool isnuint, iscc;
    std::vector<int>  PDGs={11,13,2212,211,111,0}; // e,mu,p,pi,pi0,other

    std::map<int,int> NPrimariesPDG = GetMultiplicityMap(truth,PDGs,&isnuint
                                                         ,&iscc,&nupdg,&truemode);

    int Ntrack = NPrimariesPDG[2212] + NPrimariesPDG[211];
    int Nshowr = NPrimariesPDG[111];

    if ( isnuint ){
      if( iscc ){
        if ( nupdg == 14 && Ntrack == 0 && Nshowr == 0 ) return kNumu0tr0sh;
        if ( nupdg == 14 && Ntrack == 0 && Nshowr == 1 ) return kNumu0tr1sh;
        if ( nupdg == 14 && Ntrack == 0 && Nshowr == 2 ) return kNumu0tr2sh;
        if ( nupdg == 14 && Ntrack == 0 && Nshowr > 2 )  return kNumu0trMsh;
        if ( nupdg == 14 && Ntrack == 1 && Nshowr == 0 ) return kNumu1tr0sh;
        if ( nupdg == 14 && Ntrack == 1 && Nshowr == 1 ) return kNumu1tr1sh;
        if ( nupdg == 14 && Ntrack == 1 && Nshowr == 2 ) return kNumu1tr2sh;
        if ( nupdg == 14 && Ntrack == 1 && Nshowr > 2 )  return kNumu1trMsh;
        if ( nupdg == 14 && Ntrack == 2 && Nshowr == 0 ) return kNumu2tr0sh;
        if ( nupdg == 14 && Ntrack == 2 && Nshowr == 1 ) return kNumu2tr1sh;
        if ( nupdg == 14 && Ntrack == 2 && Nshowr == 2 ) return kNumu2tr2sh;
        if ( nupdg == 14 && Ntrack == 2 && Nshowr > 2 )  return kNumu2trMsh;
        if ( nupdg == 14 && Ntrack > 2  && Nshowr == 0 ) return kNumuMtr0sh;
        if ( nupdg == 14 && Ntrack > 2  && Nshowr == 1 ) return kNumuMtr1sh;
        if ( nupdg == 14 && Ntrack > 2  && Nshowr == 2 ) return kNumuMtr2sh;
        if ( nupdg == 14 && Ntrack > 2  && Nshowr > 2 )  return kNumuMtrMsh;
        if ( nupdg == 12 && Ntrack == 0 && Nshowr == 0 ) return kNue0tr0sh;
        if ( nupdg == 12 && Ntrack == 0 && Nshowr == 1 ) return kNue0tr1sh;
        if ( nupdg == 12 && Ntrack == 0 && Nshowr == 2 ) return kNue0tr2sh;
        if ( nupdg == 12 && Ntrack == 0 && Nshowr > 2 )  return kNue0trMsh;
        if ( nupdg == 12 && Ntrack == 1 && Nshowr == 0 ) return kNue1tr0sh;
        if ( nupdg == 12 && Ntrack == 1 && Nshowr == 1 ) return kNue1tr1sh;
        if ( nupdg == 12 && Ntrack == 1 && Nshowr == 2 ) return kNue1tr2sh;
        if ( nupdg == 12 && Ntrack == 1 && Nshowr > 2 )  return kNue1trMsh;
        if ( nupdg == 12 && Ntrack == 2 && Nshowr == 0 ) return kNue2tr0sh;
        if ( nupdg == 12 && Ntrack == 2 && Nshowr == 1 ) return kNue2tr1sh;
        if ( nupdg == 12 && Ntrack == 2 && Nshowr == 2 ) return kNue2tr2sh;
        if ( nupdg == 12 && Ntrack == 2 && Nshowr > 2 )  return kNue2trMsh;
        if ( nupdg == 12 && Ntrack > 2  && Nshowr == 0 ) return kNueMtr0sh;
        if ( nupdg == 12 && Ntrack > 2  && Nshowr == 1 ) return kNueMtr1sh;
        if ( nupdg == 12 && Ntrack > 2  && Nshowr == 2 ) return kNueMtr2sh;
        if ( nupdg == 12 && Ntrack > 2  && Nshowr > 2 )  return kNueMtrMsh;
        if ( nupdg != 12 && nupdg != 14 ) return kOtherFS;
      }//cc
      else{
        if ( Ntrack == 0 && Nshowr == 0 ) return kNC0tr0sh;
        if ( Ntrack == 0 && Nshowr == 1 ) return kNC0tr1sh;
        if ( Ntrack == 0 && Nshowr == 2 ) return kNC0tr2sh;
        if ( Ntrack == 0 && Nshowr > 2 )  return kNC0trMsh;
        if ( Ntrack == 1 && Nshowr == 0 ) return kNC1tr0sh;
        if ( Ntrack == 1 && Nshowr == 1 ) return kNC1tr1sh;
        if ( Ntrack == 1 && Nshowr == 2 ) return kNC1tr2sh;
        if ( Ntrack == 1 && Nshowr > 2 )  return kNC1trMsh;
        if ( Ntrack == 2 && Nshowr == 0 ) return kNC2tr0sh;
        if ( Ntrack == 2 && Nshowr == 1 ) return kNC2tr1sh;
        if ( Ntrack == 2 && Nshowr == 2 ) return kNC2tr2sh;
        if ( Ntrack == 2 && Nshowr > 2 )  return kNC2trMsh;
        if ( Ntrack > 2  && Nshowr == 0 ) return kNCMtr0sh;
        if ( Ntrack > 2  && Nshowr == 1 ) return kNCMtr1sh;
        if ( Ntrack > 2  && Nshowr == 2 ) return kNCMtr2sh;
        if ( Ntrack > 2  && Nshowr > 2 )  return kNCMtrMsh;
      }//nc
    }
    else return kCosmicFS;
    
 return kOtherFS;
}

HitType cvn::GetHitTypeByPDG

(

int

pdg

)

Definition at line 793 of file AssignLabels.cxx.

References kElectronHit, kGammaHit, kMostlyHA, kMuonHit, kNeutronHit, kPionHit, kPiZeroHit, and kProtonHit.

Referenced by HitClassify().

  {
    switch (pdg) {
    case  11:
    case -11:
      return kElectronHit;
    case  13:
    case -13:
      return kMuonHit;
    case  2212:
      return kProtonHit;
    case  2112:
      return kNeutronHit;
    case  111:
      return kPiZeroHit;
    case  211:
    case -211:
      return kPionHit;
    case  22:
      return kGammaHit;
    default:
      return kMostlyHA;
    }
    
  }//GetHitTypeByPDG

InteractionType cvn::GetInteractionType

(

const simb::MCTruth &

truth

)

Get Interaction_t from pdg, mode and iscc. Setting pdg and mode to zero triggers cosmic ray

Definition at line 187 of file AssignLabels.cxx.

References abs(), simb::MCNeutrino::CCNC(), simb::MCTruth::GetNeutrino(), iscc, simb::kCC, kCosmic, simb::kDIS, kNC, kNueDIS, kNuElectronElastic, simb::kNuElectronElastic, kNueOther, kNueQE, kNueRes, kNumuDIS, kNumuOther, kNumuQE, kNumuRes, kNutauDIS, kNutauOther, kNutauQE, kNutauRes, kOther, simb::kQE, simb::kRes, simb::MCNeutrino::Mode(), simb::MCTruth::NeutrinoSet(), simb::MCNeutrino::Nu(), make_root_from_grid_output::pdg, and simb::MCParticle::PdgCode().

  {
    
    if(truth.NeutrinoSet())
      {
        int pdg = truth.GetNeutrino().Nu().PdgCode();
        bool iscc = truth.GetNeutrino().CCNC() == simb::kCC;
        int trueMode = truth.GetNeutrino().Mode();
        
        if(iscc)
          {
            if(abs(pdg) == 14)
              {
                switch(trueMode){
                case simb::kQE: return kNumuQE; break;
                case simb::kRes: return kNumuRes; break;
                case simb::kDIS: return kNumuDIS; break;
                default: return kNumuOther;
                }
              }
            else if(abs(pdg) == 12)
              {
                switch(trueMode){
                case simb::kQE: return kNueQE; break;
                case simb::kRes: return kNueRes; break;
                case simb::kDIS: return kNueDIS; break;
                default: return kNueOther;
                }
              }
            else if(abs(pdg) == 16)
              {
                switch(trueMode){
                case simb::kQE: return kNutauQE; break;
                case simb::kRes: return kNutauRes; break;
                case simb::kDIS: return kNutauDIS; break;
                default: return kNutauOther;
                }
              }
          }
        else if(trueMode==simb::kNuElectronElastic){
          return kNuElectronElastic;
        }
        else return kNC;
      }
    else return kCosmic;
    
    return kOther;
  }

InteractionType cvn::GetInteractionTypeFromSlice	(	int	pdg,
		bool	iscc,
		int	trueMode
	)

Definition at line 236 of file AssignLabels.cxx.

References abs(), simb::kDIS, kNC, kNueDIS, kNuElectronElastic, simb::kNuElectronElastic, kNueOther, kNueQE, kNueRes, kNumuDIS, kNumuOther, kNumuQE, kNumuRes, kNutauDIS, kNutauOther, kNutauQE, kNutauRes, kOther, simb::kQE, and simb::kRes.

Referenced by SliceClassify().

  {
    
    if(iscc){
      
      if(abs(pdg) == 14){
        switch(trueMode){
        case simb::kQE: return kNumuQE; break;
        case simb::kRes: return kNumuRes; break;
        case simb::kDIS: return kNumuDIS; break;
        default: return kNumuOther;
        }
      }
      else if(abs(pdg) == 12){
        switch(trueMode){
        case simb::kQE: return kNueQE; break;
        case simb::kRes: return kNueRes; break;
        case simb::kDIS: return kNueDIS; break;
        default: return kNueOther;
        }
      }
      else if(abs(pdg) == 16){
        switch(trueMode){
        case simb::kQE: return kNutauQE; break;
        case simb::kRes: return kNutauRes; break;
        case simb::kDIS: return kNutauDIS; break;
        default: return kNutauOther;
        }
      }
      
    }
    else if(trueMode==simb::kNuElectronElastic){
      return kNuElectronElastic;
    }
    else return kNC;
    
    return kOther;
  }

std::map< int, int > cvn::GetMultiplicityMap	(	const art::Ptr< simb::MCTruth >	truth,
		std::vector< int >	PDGs,
		bool *	isnuint,
		bool *	iscc,
		int *	nupdg,
		int *	trueMode
	)

Definition at line 24 of file AssignLabels.cxx.

References abs(), simb::MCNeutrino::CCNC(), simb::MCParticle::E(), simb::MCTruth::GetNeutrino(), simb::MCTruth::GetParticle(), MECModelEnuComparisons::i, simb::kCC, simb::MCNeutrino::Mode(), simb::MCTruth::NeutrinoSet(), simb::MCTruth::NParticles(), simb::MCNeutrino::Nu(), make_root_from_grid_output::pdg, simb::MCParticle::PdgCode(), and simb::MCParticle::StatusCode().

Referenced by GetFinalStateProngsType(), and GetFinalStateType().

  {
    //
    art::ServiceHandle<cheat::BackTracker> bt;

    std::map<int,int> NumberOfPrimariesPDG;
    
    for ( unsigned int i = 0; i < PDGs.size(); ++i)
      NumberOfPrimariesPDG.insert ( std::pair<int,int>(PDGs[i],0) );
    
    if(truth->NeutrinoSet()){

      *isnuint  = true;
      *iscc     = truth->GetNeutrino().CCNC() == simb::kCC;
      *nupdg    = truth->GetNeutrino().Nu().PdgCode();
      *trueMode = truth->GetNeutrino().Mode();
      
      for( int iPart = 0; iPart < truth->NParticles(); iPart++){
        
        const simb::MCParticle daughter = truth->GetParticle(iPart);
        int pdg = abs(daughter.PdgCode());        
        
        if ( daughter.E() < 0.03 )        continue;
        if ( daughter.StatusCode() != 1 ) continue;

        if(NumberOfPrimariesPDG.find(pdg) != NumberOfPrimariesPDG.end())
          ++NumberOfPrimariesPDG[pdg];          
        else 
          ++NumberOfPrimariesPDG[0];
        
      }//Primaries      

    }
    else 
      *isnuint  = false;;
    
    return NumberOfPrimariesPDG;
  }

ParentParticleType cvn::GetParentParticleType	(	const art::Ptr< simb::MCTruth >	truth,
		const art::Event &	evt
	)

Definition at line 1011 of file AssignLabels.cxx.

References simb::MCNeutrino::CCNC(), evt, flux, simb::MCFlux::fptype, simb::MCTruth::GetNeutrino(), iscc, simb::kCC, kCosmicPP, kKaNuCCPP, kKaNuNCPP, kMuNuCCPP, kMuNuNCPP, kOtherPP, kPiNuCCPP, kPiNuNCPP, kUnknownPP, and simb::MCTruth::NeutrinoSet().

Referenced by cvn::CVNEventDump::analyze(), cvn::CVNCosmicEventDump::analyze(), and cvn::CVNAddTrainingData::produce().

  {
      int parentpdg = -1;
      if(truth->NeutrinoSet()){
        bool iscc  = truth->GetNeutrino().CCNC() == simb::kCC;
        art::FindManyP<simb::MCFlux> fmFlux({truth}, evt, art::InputTag("generator"));
        if(fmFlux.isValid()){
          std::vector< art::Ptr<simb::MCFlux> > fluxes = fmFlux.at(0);
          if(!fluxes.empty()){
            const simb::MCFlux& flux = *fluxes[0];
            parentpdg = flux.fptype;
          }
        }
        if(iscc){
          switch(parentpdg) {
          case 211:
            return kPiNuCCPP;
          case -211: 
            return kPiNuCCPP;
          case 321:
            return kKaNuCCPP;
          case -321: 
            return kKaNuCCPP;
          case 130:
          case 310:
          case 311:
            return kKaNuCCPP;
          case 13:
            return kMuNuCCPP;
          case -13:
            return kMuNuCCPP;   
          } 
        }
        else if(!iscc){
          switch(parentpdg) {
          case 211:
            return kPiNuNCPP;
          case -211:
            return kPiNuNCPP;
          case 321:
            return kKaNuNCPP;
          case -321:
            return kKaNuNCPP;
          case 130:
          case 310:
          case 311:
            return kKaNuNCPP;
          case 13:
            return kMuNuNCPP;
          case -13:
            return kMuNuNCPP;
          } 
        }
        return kOtherPP;
      }
      else return kCosmicPP;
      return kUnknownPP;
  
  }// GetParentParticleType

int cvn::GetPDGByPType

(

PType

ptype

)

Definition at line 5 of file ProngType.cxx.

References kElectron, kEM, kGamma, kHadron, kMuon, kNeutron, kPion, kPiZero, and kProton.

Referenced by cvnprongtf::CVNProngTF::produce(), and cvntf::CVNProngEvaluatorTF::produce().

  {

    switch (ptype) {
    case  kElectron:
      return 11;
    case  kMuon:
      return  13;
    case  kProton:
      return  2212;
    case  kNeutron:
      return  2112;
    case  kPion:
      return  211;
    case  kPiZero:
      return  111;
    case  kGamma:
      return  22;
      
      // In the case of composite CVN scores like EMID and HadronID,
      // represent the PDG code of those types as the product of the two 
      // pdg codes. 
      // Eg. pdgEM = 22 * 11 = 242

      // This way we can make comparisons with %
      //     pdgEM % 11 == 0
    case kEM:
      // pdg 22*11 = 242 describes a bound state of an up and charm quark with spin 1/2
      // This breaks quantum mechanics, so we're safe
      return 22 * 11;
    case kHadron:
      // pdg 2212*211 = 466732 describes a tb's baryon in the 4th radial eigen state
      // and an unphysical orbital momentum state.
      // Pretty sure we won't have to deal with those
      return 2212 * 211;
    default:
      return 0;
    }

  }// GetPDGByPType

std::vector< std::pair< int, double > > cvn::GetProngEfficiencyByPDG	(	const rb::Cluster &	prong,
		const rb::Cluster &	slice
	)

Definition at line 1278 of file AssignLabels.cxx.

References rb::Cluster::AllCells(), dE, cheat::BackTracker::HitCollectionEfficiency(), cheat::BackTracker::HitsToTrackIDE(), kEmpty, kUnknown, geo::kXorY, make_pair(), rb::Cluster::NXCell(), rb::Cluster::NYCell(), simb::MCParticle::PdgCode(), and cheat::BackTracker::TrackIDToParticle().

Referenced by bpf::BPFCVNAna::analyze().

  {
    
    art::ServiceHandle<cheat::BackTracker> bt;
    
    std::vector< std::pair<int, double> > efficiencybypdg;
    
    // Check nhits > 0 and number of associated trackIDs > 0 for the prong
    ProngType pType3D = kUnknown;
    if( prong.NXCell() == 0 || 
        prong.NYCell() == 0) pType3D = kEmpty;
    
    std::vector<cheat::TrackIDE> ParticleTrackIDs;
    if( pType3D != kEmpty )
      ParticleTrackIDs  = bt->HitsToTrackIDE(prong.AllCells());
    
    if( ParticleTrackIDs.size()  < 1 ) pType3D = kEmpty;
    
    if( pType3D != kEmpty ){
      
      for (unsigned int iID=0; iID < ParticleTrackIDs.size(); ++iID){
        std::set<int> TrackID;
        std::map<int, double> purMap;
        double dE, tE;
        std::map<int, int> parents;
        
        int iPartID = ParticleTrackIDs[iID].trackID;
        TrackID.insert(iPartID);
        int iPartPDG = bt->TrackIDToParticle(iPartID)->PdgCode();
        
        double iPartEff = bt->HitCollectionEfficiency(TrackID, prong.AllCells(),slice.AllCells(),
                                                      geo::kXorY, &purMap,true, &dE, &tE);

        efficiencybypdg.push_back(std::make_pair(iPartPDG, iPartEff));
        
      }//iTrackID      
    }
    
    return efficiencybypdg;
  }// ProngEfficiencyByPDG

std::vector< std::pair< int, double > > cvn::GetProngEnergyByPDG

(

const rb::Cluster &

prong

)

Definition at line 1321 of file AssignLabels.cxx.

References rb::Cluster::AllCells(), simb::MCParticle::E(), cheat::BackTracker::HitsToTrackIDE(), kEmpty, kUnknown, make_pair(), rb::Cluster::NXCell(), rb::Cluster::NYCell(), simb::MCParticle::PdgCode(), and cheat::BackTracker::TrackIDToParticle().

Referenced by bpf::BPFCVNAna::analyze().

  {
    
    art::ServiceHandle<cheat::BackTracker> bt;
    
    std::vector< std::pair<int, double> > energybypdg;
    
    // Check nhits > 0 and number of associated trackIDs > 0 for the prong
    ProngType pType3D = kUnknown;
    if( prong.NXCell() == 0 || 
        prong.NYCell() == 0) pType3D = kEmpty;
    
    std::vector<cheat::TrackIDE> ParticleTrackIDs;
    if( pType3D != kEmpty )
      ParticleTrackIDs  = bt->HitsToTrackIDE(prong.AllCells());
    
    if( ParticleTrackIDs.size()  < 1 ) pType3D = kEmpty;
    
    if( pType3D != kEmpty ){
      
      for (unsigned int iID=0; iID < ParticleTrackIDs.size(); ++iID){
        std::set<int> TrackID;
        std::map<int, double> purMap;
        std::map<int, int> parents;
        
        int iPartID = ParticleTrackIDs[iID].trackID;
        TrackID.insert(iPartID);
        int  iPartPDG = bt->TrackIDToParticle(iPartID)->PdgCode();
        double iPartE = bt->TrackIDToParticle(iPartID)->E();
        
        energybypdg.push_back(std::make_pair(iPartPDG, iPartE));
        
      }//iTrackID      
    }
    
    return energybypdg;
  }// ProngEnergyByPDG  

std::vector< std::pair< int, double > > cvn::GetProngPurityByPDG

(

const rb::Cluster &

prong

)

Definition at line 1238 of file AssignLabels.cxx.

References rb::Cluster::AllCells(), cheat::BackTracker::HitCollectionPurity(), cheat::BackTracker::HitsToTrackIDE(), kEmpty, kUnknown, make_pair(), rb::Cluster::NXCell(), rb::Cluster::NYCell(), simb::MCParticle::PdgCode(), and cheat::BackTracker::TrackIDToParticle().

Referenced by bpf::BPFCVNAna::analyze().

  {
    
    art::ServiceHandle<cheat::BackTracker> bt;

    std::vector< std::pair<int, double> > puritybypdg;

    // Check nhits > 0 and number of associated trackIDs > 0 for the prong
    ProngType   pType3D = kUnknown;    
    if( prong.NXCell() == 0 || 
        prong.NYCell() == 0) pType3D = kEmpty;
    
    
    std::vector<cheat::TrackIDE> ParticleTrackIDs;    
    if( pType3D != kEmpty )
      ParticleTrackIDs  = bt->HitsToTrackIDE(prong.AllCells());
    
    if( ParticleTrackIDs.size()  < 1 ) pType3D = kEmpty;
    
    if( pType3D != kEmpty ){
      
      for (unsigned int iID=0; iID < ParticleTrackIDs.size(); ++iID){
        std::set<int> TrackID;
        std::map<int, double> purMap;
        std::map<int, int> parents;
        
        int iPartID = ParticleTrackIDs[iID].trackID;
        TrackID.insert(iPartID);
        int iPartPDG = bt->TrackIDToParticle(iPartID)->PdgCode();
        double iPartPur = bt->HitCollectionPurity(TrackID, prong.AllCells(),
                                                  &purMap, &parents,true);
        puritybypdg.push_back(std::make_pair(iPartPDG, iPartPur));
        
      }//iTrackID
    }
    
    return puritybypdg;
  }// ProngPurityByPDG

SimpleNeutrinoType cvn::GetSimpleNeutrinoType

(

const InteractionType &

intType

)

Get Interaction_t from pdg, mode and iscc. Setting pdg and mode to zero triggers cosmic ray

Definition at line 9 of file SimpleNeutrinoType.cxx.

References kCos, kCosmic, kNC, kNSimpleNeutrinoType, kNue, kNueOther, kNumu, kNumuOther, kNuNC, kNutau, and kNutauOther.

Referenced by cvn::TrainingData::FillSimpleOutputVector().

  {

      if (intType <= kNumuOther)  return kNumu;
      if (intType <= kNueOther)   return kNue;
      if (intType <= kNutauOther) return kNutau;
      if (intType == kNC)         return kNuNC;
      if (intType == kCosmic)     return kCos;

      return kNSimpleNeutrinoType;

  }

ProngType cvn::GetTypeByPDG

(

int

pdg

)

Definition at line 1072 of file AssignLabels.cxx.

References kElectron, kGamma, kMuon, kNeutron, kOtherPDG, kPion, kPiZero, and kProton.

Referenced by ProngClassify().

  {
    
    switch (pdg) {
    case  11:
    case -11:
      return kElectron;
    case  13:
    case -13:
      return kMuon;
    case  2212:
      return kProton;
    case  2112:
      return kNeutron;
    case  111:
      return kPiZero;
    case  211:
    case -211:
      return kPion;
    case  22:
      return kGamma;
    default:
      return kOtherPDG;
    }
    
  }// GetTypeByInt

HitType cvn::HitClassify	(	art::Ptr< rb::CellHit >	Hit,
		HitType *	hType,
		double *	hPurity
	)

Definition at line 820 of file AssignLabels.cxx.

References GetHitTypeByPDG(), cheat::BackTracker::HaveTruthInfo(), cheat::BackTracker::HitCollectionPurity(), cheat::BackTracker::HitToTrackIDE(), kEmptyHit, kMostlyEM, kMostlyHA, kMostlyMU, kUnknownHit, simb::MCParticle::PdgCode(), art::PtrVector< T >::push_back(), and cheat::BackTracker::TrackIDToParticle().

Referenced by cvn::RegPixelMapProducer::CreateMapGivenBoundary(), cvn::PixelMapProducer::CreateMapGivenBoundary(), cvn::RegPixelMapProducer::CreateMapGivenShowerVertex(), and cvn::RegPixelMapProducer::CreateMapGivenVertex().

  {
    
    art::ServiceHandle<cheat::BackTracker> bt;
    std::vector<cheat::TrackIDE>           ParticleTrackIDs;
    
    *hType = kUnknownHit;
    *hPurity = 0.0;

    HitType typeEmpty = *hType;
    if ( !(bt->HaveTruthInfo()) ) return typeEmpty;
    
    // Check the hit is associated to any true particles
    art::PtrVector<rb::CellHit> hit;
    hit.push_back( Hit );
    ParticleTrackIDs  = bt->HitToTrackIDE( Hit );
    if( ParticleTrackIDs.size()  < 1 ) *hType = kEmptyHit;
    
    if( *hType != kEmptyHit ){
      double totalPur   = 0.0;
      double highestPur = 0.0;
      int highestPurPDG;
      double fracEm = 0.0;
      double fracHa = 0.0;
      double fracMu = 0.0;
      
      for (unsigned int iID=0; iID < ParticleTrackIDs.size(); ++iID){
        std::set<int> TrackID;
        std::map<int, double> purMap;
        std::map<int, int> parents;
        
        // Ensure we label it as the highest purity true hit
        int iPartID = ParticleTrackIDs[iID].trackID;
        TrackID.insert(iPartID);
        double iPartPur = bt->HitCollectionPurity(TrackID, hit,
                                                  &purMap, &parents,true);
        int    iPartPDG = bt->TrackIDToParticle(iPartID)->PdgCode(); 

        // Keep track of what fraction is EM, Hadronic or Muon
        switch (iPartPDG) {
        case  11:
        case -11:
        case  22:
        case  111:
          fracEm += iPartPur;
          break;
        case  13:
        case -13:
          fracMu += iPartPur;
          break;
        case  2212:
        case  2112:
        case  211:
        case -211:
          fracHa += iPartPur;
          break;
        default:
          fracHa += iPartPur;
        }

        // Label hits by highest Pur contributor....*
        if( iPartPur > highestPur ){
          highestPur    = iPartPur;
          highestPurPDG = iPartPDG;
          *hType  = GetHitTypeByPDG(highestPurPDG);
          *hPurity = highestPur;
        }
        
        totalPur+=iPartPur; //don't need this
        
      }// for trackIDs 

      // .... *unless that "highest" is not that high
      if ( highestPur < 0.5 ){
        if ( fracMu > fracEm && fracMu > fracHa )
          *hType = kMostlyMU;
        else{
          if ( fracEm > fracHa ) 
            *hType = kMostlyEM;
          else
            *hType = kMostlyHA;       
        }
      } // ....
      
    }// type != empty
    
    HitType type = *hType;
    
    return type;
    
  }//HitClassify

int cvn::HitNuIndex	(	art::Ptr< rb::CellHit >	Hit,
		const rb::Cluster &	evtCluster
	)

Definition at line 968 of file AssignLabels.cxx.

References rb::Cluster::Cell(), cheat::BackTracker::HaveTruthInfo(), cheat::BackTracker::HitToTrackIDE(), rb::Cluster::NCell(), art::PtrVector< T >::push_back(), and cheat::BackTracker::SliceToNeutrinoIndex().

Referenced by cvn::PixelMapProducer::CreateMapGivenBoundary().

  {
    
    art::ServiceHandle<cheat::BackTracker> bt;
    std::vector<cheat::TrackIDE>           ParticleTrackIDs;

    if ( !(bt->HaveTruthInfo()) ) return -5;

    // Check the hit is associated to any true particles
    art::PtrVector<rb::CellHit> hit;
    hit.push_back( Hit );
    ParticleTrackIDs  = bt->HitToTrackIDE( Hit );
    if( ParticleTrackIDs.size() < 1 ) return 0;

    // Make vectors of the one hit and all the hits
    const rb::CellHit* oneHit = &*Hit;
    std::vector<const rb::CellHit*> cell_vector;

    cell_vector.push_back(oneHit);

    std::vector<const rb::CellHit*> cluster_vector;

    for ( size_t iHit = 0; iHit < evtCluster.NCell(); ++iHit )
      {
        art::Ptr< rb::CellHit >evtHit = evtCluster.Cell(iHit);

        const rb::CellHit* oneevtHit = &*evtHit;

        cluster_vector.push_back(oneevtHit);

      }

    // Get indices of nus contributing to cell_vector (ordered by pur.)
    std::vector<int> nuIndices = bt->SliceToNeutrinoIndex(cell_vector, cluster_vector, true);
    if(nuIndices.empty()) return 0;

    return nuIndices[0] + 1; // reserve 0 for noise
    
  }//HitNuIndex

int cvn::HitObject	(	art::Ptr< rb::CellHit >	Hit,
		int &	hObjectCount,
		std::map< int, int > &	hObjectMap
	)

Definition at line 914 of file AssignLabels.cxx.

References cheat::BackTracker::HaveTruthInfo(), cheat::BackTracker::HitCollectionPurity(), cheat::BackTracker::HitToTrackIDE(), maxID(), and art::PtrVector< T >::push_back().

Referenced by cvn::PixelMapProducer::CreateMapGivenBoundary().

  {
    int object = 0;
    
    art::ServiceHandle<cheat::BackTracker> bt;
    std::vector<cheat::TrackIDE>           ParticleTrackIDs;

    if ( !(bt->HaveTruthInfo()) ) return 0;

    // Check the hit is associated to any true particles
    art::PtrVector<rb::CellHit> hit;
    hit.push_back( Hit );
    ParticleTrackIDs  = bt->HitToTrackIDE( Hit );
    if( ParticleTrackIDs.size() < 1 ) return 0;
    
    double maxPur = 0.0;
    int    maxID  = 0;
    
    for (unsigned int iID=0; iID < ParticleTrackIDs.size(); ++iID){
      std::set<int> TrackID;
      std::map<int, double> purMap;
      std::map<int, int> parents;
      
      int iPartID = ParticleTrackIDs[iID].trackID;
      TrackID.insert(iPartID);
      double iPartPur = bt->HitCollectionPurity(TrackID, hit,
                                                &purMap, &parents,true);

      if(iPartPur > maxPur){
        maxPur = iPartPur;
        maxID = iPartID;
      }
      
    }// for trackIDs 

    // Label hits by highest Pur contributor....
    if( hObjectMap.count(maxID) > 0){
      object = hObjectMap[maxID];
    }
    else {
      hObjectCount++;
      hObjectMap[maxID] = hObjectCount;
      object = hObjectMap[maxID];
    }

    return object;
    
  }//HitObject

ParticlesType cvn::MakeParticleLabel	(	const art::Ptr< simb::MCTruth >	truth,
		std::string	NuType
	)

Definition at line 303 of file AssignLabels.cxx.

References abs(), simb::MCParticle::E(), energy, simb::MCTruth::GetParticle(), kOther_PT, simb::MCParticle::Mass(), simb::MCTruth::NParticles(), simb::MCParticle::PdgCode(), elec2geo::pos, fillBadChanDBTables::result, simb::MCParticle::StatusCode(), string, and registry_explorer::v.

Referenced by ParticlesSliceClassify().

                                                   {

    ParticlesType sType = kOther_PT;
    
    //int count_e,count_mu,count_g,count_pi0,count_pi,count_n,count_p;
    //count_e=count_mu=count_g=count_pi0=count_pi=count_n=count_p=0;
    int count_e,count_mu,count_pi0,count_pi,count_n,count_p;
    count_e=count_mu=count_pi0=count_pi=count_n=count_p=0;
    std::string result;
    std::stringstream sstm;
    
    for (int iPart = 0; iPart < truth->NParticles(); ++iPart){
      if (truth->GetParticle(iPart).StatusCode() != 1) continue;
      int pdgcode = abs(truth->GetParticle(iPart).PdgCode());
      double mass = truth->GetParticle(iPart).Mass();
      double energy = truth->GetParticle(iPart).E(0);

      if(energy-mass<0.04) continue;

      if (pdgcode==11) count_e++;
      if (pdgcode==13) count_mu++;
      if (pdgcode==111) count_pi0++;
      if (pdgcode==211) count_pi++;
      if (pdgcode==2112) count_n++;
      if (pdgcode==2212) count_p++;
    }
    
    if (count_e>1) count_e=1; // Do not care to distinguish between events with 1 and 2 electrons
    if (count_mu>1) count_mu=1; // Do not care to distinguish between events with 1 and 2 muons
    if (count_pi0>1) count_pi0=2; // Above 1 pi0 it is one category
    if (count_pi>1) count_pi=2; // Above 1 pi it is one category
    if (count_p>1) count_p=2; // Above 1 p it is one category
    if (count_n>1) count_n=2; // Above 1 n it is one category
    
    sstm << NuType << "_e_" << count_e << "_mu_" << count_mu
         << "_pi0_" << count_pi0 << "_pi_" << count_pi << "_n_"
         << count_n << "_p_" << count_p;
    result = sstm.str();
    
  std::vector<std::string> v = {
"nE_e_1_mu_0_pi0_0_pi_0_n_0_p_1",
"nM_e_0_mu_1_pi0_0_pi_0_n_0_p_1",
"nM_e_0_mu_1_pi0_0_pi_1_n_0_p_1",
"nU_e_0_mu_0_pi0_1_pi_0_n_0_p_1",
"nE_e_1_mu_0_pi0_0_pi_1_n_0_p_1",
"nE_e_1_mu_0_pi0_0_pi_0_n_2_p_2",
"nU_e_0_mu_0_pi0_0_pi_1_n_0_p_1",
"nE_e_1_mu_0_pi0_0_pi_1_n_1_p_0",
"nM_e_0_mu_1_pi0_0_pi_1_n_1_p_0",
"nU_e_0_mu_0_pi0_0_pi_1_n_1_p_0",
"nU_e_0_mu_0_pi0_1_pi_0_n_1_p_0",
"nM_e_0_mu_1_pi0_0_pi_0_n_2_p_2",
"nU_e_0_mu_0_pi0_1_pi_0_n_1_p_1",
"nM_e_0_mu_1_pi0_0_pi_0_n_1_p_1",
"nE_e_1_mu_0_pi0_0_pi_0_n_1_p_1",
"nU_e_0_mu_0_pi0_1_pi_0_n_2_p_2",
"nE_e_1_mu_0_pi0_0_pi_0_n_1_p_2",
"nE_e_1_mu_0_pi0_0_pi_1_n_2_p_2",
"nM_e_0_mu_1_pi0_0_pi_1_n_1_p_1",
"nE_e_1_mu_0_pi0_0_pi_0_n_0_p_2",
"nU_e_0_mu_0_pi0_0_pi_2_n_0_p_1",
"nM_e_0_mu_1_pi0_0_pi_0_n_1_p_2",
"nU_e_0_mu_0_pi0_0_pi_0_n_2_p_2",
"nM_e_0_mu_1_pi0_0_pi_0_n_0_p_2",
"nE_e_1_mu_0_pi0_1_pi_0_n_0_p_1",
"nE_e_1_mu_0_pi0_0_pi_1_n_1_p_1",
"nU_e_0_mu_0_pi0_1_pi_1_n_1_p_0",
"nU_e_0_mu_0_pi0_0_pi_1_n_1_p_1",
"nM_e_0_mu_1_pi0_1_pi_0_n_0_p_1",
"nU_e_0_mu_0_pi0_0_pi_2_n_1_p_0",
"nU_e_0_mu_0_pi0_0_pi_1_n_2_p_2",
"nE_e_1_mu_0_pi0_1_pi_1_n_0_p_1",
"nU_e_0_mu_0_pi0_1_pi_2_n_0_p_1",
"nM_e_0_mu_1_pi0_0_pi_1_n_2_p_2",
"nE_e_1_mu_0_pi0_0_pi_2_n_0_p_1",
"nU_e_0_mu_0_pi0_1_pi_2_n_2_p_2",
"nU_e_0_mu_0_pi0_1_pi_2_n_1_p_0",
"nU_e_0_mu_0_pi0_1_pi_1_n_0_p_1",
"nU_e_0_mu_0_pi0_2_pi_2_n_0_p_1",
"nU_e_0_mu_0_pi0_0_pi_0_n_0_p_1",
"nU_e_0_mu_0_pi0_0_pi_2_n_2_p_2",
"nU_e_0_mu_0_pi0_0_pi_2_n_1_p_1",
"nM_e_0_mu_1_pi0_0_pi_1_n_0_p_2",
"nU_e_0_mu_0_pi0_1_pi_0_n_0_p_0",
"nE_e_1_mu_0_pi0_0_pi_1_n_0_p_2",
"nU_e_0_mu_0_pi0_0_pi_1_n_0_p_2",
"nT_e_0_mu_0_pi0_0_pi_2_n_0_p_1",
"nU_e_0_mu_0_pi0_1_pi_0_n_2_p_0",
"nU_e_0_mu_0_pi0_0_pi_1_n_2_p_0",
"nM_e_0_mu_1_pi0_1_pi_1_n_0_p_1",
"nE_e_1_mu_0_pi0_1_pi_2_n_0_p_1",
"nU_e_0_mu_0_pi0_2_pi_2_n_1_p_0",
"nM_e_0_mu_1_pi0_0_pi_1_n_2_p_1",
"nU_e_0_mu_0_pi0_1_pi_0_n_0_p_2",
"nU_e_0_mu_0_pi0_0_pi_1_n_2_p_1",
"nM_e_0_mu_1_pi0_1_pi_2_n_0_p_1",
"nM_e_0_mu_1_pi0_0_pi_2_n_1_p_0",
"nE_e_1_mu_0_pi0_0_pi_2_n_1_p_0",
"nU_e_0_mu_0_pi0_0_pi_0_n_1_p_0",
"nM_e_0_mu_1_pi0_0_pi_2_n_0_p_1",
"nM_e_0_mu_1_pi0_0_pi_1_n_1_p_2",
"nU_e_0_mu_0_pi0_2_pi_0_n_0_p_1",
"nU_e_0_mu_0_pi0_1_pi_1_n_2_p_2",
"nU_e_0_mu_0_pi0_1_pi_0_n_2_p_1",
"nU_e_0_mu_0_pi0_0_pi_1_n_1_p_2",
"nE_e_1_mu_0_pi0_1_pi_2_n_1_p_0",
"nE_e_1_mu_0_pi0_0_pi_0_n_2_p_1",
"nT_e_0_mu_0_pi0_0_pi_2_n_1_p_0",
"nM_e_0_mu_1_pi0_1_pi_0_n_0_p_2",
"nE_e_1_mu_0_pi0_0_pi_1_n_2_p_1",
"nU_e_0_mu_0_pi0_2_pi_2_n_1_p_1",
"nM_e_0_mu_1_pi0_1_pi_1_n_1_p_0",
"nM_e_0_mu_1_pi0_0_pi_2_n_1_p_1",
"nU_e_0_mu_0_pi0_0_pi_0_n_2_p_1",
"nM_e_0_mu_1_pi0_1_pi_0_n_2_p_2",
"nE_e_1_mu_0_pi0_0_pi_2_n_1_p_1",
"nU_e_0_mu_0_pi0_2_pi_2_n_2_p_2",
"nT_e_0_mu_0_pi0_2_pi_2_n_0_p_1",
"nM_e_0_mu_1_pi0_1_pi_1_n_1_p_1",
"nU_e_0_mu_0_pi0_1_pi_2_n_1_p_1",
"nU_e_0_mu_0_pi0_0_pi_2_n_2_p_1",
"nU_e_0_mu_0_pi0_0_pi_2_n_0_p_2",
"nM_e_0_mu_1_pi0_0_pi_0_n_2_p_1",
"nE_e_1_mu_0_pi0_1_pi_0_n_1_p_1",
"nM_e_0_mu_1_pi0_1_pi_0_n_1_p_1",
"nE_e_1_mu_0_pi0_1_pi_1_n_1_p_0",
"nE_e_1_mu_0_pi0_0_pi_2_n_2_p_2",
"nU_e_0_mu_0_pi0_1_pi_1_n_1_p_1",
"nE_e_1_mu_0_pi0_2_pi_2_n_2_p_2",
"nE_e_1_mu_0_pi0_2_pi_2_n_0_p_1",
"nE_e_1_mu_0_pi0_1_pi_1_n_1_p_1",
"nU_e_0_mu_0_pi0_1_pi_1_n_2_p_0",
"nU_e_0_mu_0_pi0_1_pi_1_n_1_p_2",
"nU_e_0_mu_0_pi0_0_pi_2_n_1_p_2",
"nU_e_0_mu_0_pi0_0_pi_0_n_1_p_2",
"nT_e_0_mu_0_pi0_2_pi_2_n_2_p_2",
"nT_e_0_mu_0_pi0_2_pi_2_n_1_p_0",
"nT_e_0_mu_0_pi0_1_pi_2_n_0_p_1",
"nM_e_0_mu_1_pi0_1_pi_2_n_1_p_0",
"nE_e_1_mu_0_pi0_1_pi_0_n_0_p_2",
"nM_e_0_mu_1_pi0_1_pi_2_n_1_p_1",
"nM_e_0_mu_1_pi0_0_pi_1_n_2_p_0",
"nU_e_0_mu_0_pi0_1_pi_2_n_2_p_0",
"nU_e_0_mu_0_pi0_0_pi_0_n_0_p_0",
"nT_e_0_mu_0_pi0_1_pi_2_n_1_p_0",
"nM_e_0_mu_1_pi0_2_pi_2_n_0_p_1",
"nE_e_1_mu_0_pi0_0_pi_2_n_1_p_2",
"nE_e_1_mu_0_pi0_0_pi_1_n_0_p_0",
"nU_e_0_mu_0_pi0_2_pi_1_n_2_p_2",
"nU_e_0_mu_0_pi0_2_pi_1_n_0_p_1",
"nT_e_0_mu_0_pi0_0_pi_2_n_1_p_1",
"nM_e_0_mu_1_pi0_0_pi_2_n_0_p_2",
"nE_e_1_mu_0_pi0_2_pi_2_n_1_p_0",
"nE_e_1_mu_0_pi0_0_pi_1_n_1_p_2",
"nU_e_0_mu_0_pi0_2_pi_0_n_2_p_2",
"nU_e_0_mu_0_pi0_0_pi_2_n_2_p_0",
"nU_e_0_mu_0_pi0_0_pi_0_n_1_p_1",
"nM_e_0_mu_1_pi0_1_pi_2_n_0_p_2",
"nM_e_0_mu_1_pi0_0_pi_1_n_0_p_0",
"nE_e_1_mu_0_pi0_1_pi_0_n_2_p_2",
"nU_e_0_mu_0_pi0_2_pi_2_n_0_p_2",
"nU_e_0_mu_0_pi0_1_pi_1_n_0_p_2",
"nM_e_0_mu_1_pi0_1_pi_2_n_2_p_0",
"nM_e_0_mu_1_pi0_0_pi_2_n_2_p_2",
"nE_e_1_mu_0_pi0_0_pi_1_n_2_p_0",
"nU_e_0_mu_0_pi0_2_pi_0_n_1_p_1",
"nM_e_0_mu_1_pi0_2_pi_2_n_2_p_2",
"nE_e_1_mu_0_pi0_1_pi_2_n_2_p_2",
"nE_e_1_mu_0_pi0_1_pi_2_n_2_p_1",
"nE_e_1_mu_0_pi0_0_pi_0_n_1_p_0",
"nU_e_0_mu_0_pi0_2_pi_2_n_2_p_0",
"nU_e_0_mu_0_pi0_2_pi_0_n_1_p_0",
"nU_e_0_mu_0_pi0_1_pi_2_n_2_p_1",
"nU_e_0_mu_0_pi0_1_pi_2_n_1_p_2",
"nU_e_0_mu_0_pi0_1_pi_1_n_2_p_1",
"nU_e_0_mu_0_pi0_0_pi_0_n_2_p_0",
"nT_e_0_mu_0_pi0_0_pi_1_n_0_p_1",
"nM_e_0_mu_1_pi0_1_pi_2_n_2_p_2",
"nM_e_0_mu_1_pi0_0_pi_0_n_1_p_0",
"nE_e_1_mu_0_pi0_2_pi_1_n_0_p_1",
"nE_e_1_mu_0_pi0_1_pi_1_n_2_p_2",
"nE_e_1_mu_0_pi0_1_pi_0_n_1_p_2",
"nU_e_0_mu_0_pi0_1_pi_0_n_1_p_2",
"nU_e_0_mu_0_pi0_0_pi_2_n_0_p_0",
"nT_e_0_mu_0_pi0_1_pi_2_n_2_p_2",
"nT_e_0_mu_0_pi0_0_pi_2_n_2_p_2",
"nM_e_0_mu_1_pi0_2_pi_2_n_1_p_0",
"nM_e_0_mu_1_pi0_1_pi_1_n_0_p_2",
"nM_e_0_mu_1_pi0_1_pi_0_n_1_p_2",
"nE_e_1_mu_0_pi0_1_pi_2_n_1_p_2",
"nE_e_1_mu_0_pi0_1_pi_1_n_0_p_2",
"nU_e_0_mu_0_pi0_2_pi_1_n_1_p_0",
"nT_e_0_mu_1_pi0_0_pi_0_n_0_p_1",
"nM_e_0_mu_1_pi0_2_pi_2_n_1_p_1",
"nM_e_0_mu_1_pi0_2_pi_1_n_0_p_1",
"nM_e_0_mu_1_pi0_1_pi_1_n_2_p_2",
"nM_e_0_mu_1_pi0_1_pi_1_n_1_p_2",
"nM_e_0_mu_1_pi0_0_pi_2_n_2_p_1",
"nE_e_1_mu_0_pi0_2_pi_1_n_2_p_2",
"nE_e_1_mu_0_pi0_0_pi_2_n_0_p_2",
"nE_e_1_mu_0_pi0_0_pi_0_n_0_p_0",
"nU_e_0_mu_0_pi0_2_pi_2_n_2_p_1",
"nU_e_0_mu_0_pi0_2_pi_2_n_1_p_2",
"nU_e_0_mu_0_pi0_2_pi_1_n_2_p_1",
"nU_e_0_mu_0_pi0_2_pi_0_n_1_p_2",
"nU_e_0_mu_0_pi0_1_pi_2_n_0_p_2",
"nU_e_0_mu_0_pi0_0_pi_0_n_0_p_2",
"nM_e_0_mu_1_pi0_2_pi_0_n_0_p_1",
"nE_e_1_mu_0_pi0_2_pi_2_n_1_p_1",
"nE_e_1_mu_0_pi0_2_pi_0_n_0_p_1",
"nE_e_1_mu_0_pi0_1_pi_1_n_2_p_1",
"nT_e_0_mu_1_pi0_0_pi_1_n_0_p_1",
"nM_e_0_mu_1_pi0_2_pi_0_n_0_p_2",
"nM_e_0_mu_1_pi0_0_pi_2_n_2_p_0",
"nM_e_0_mu_1_pi0_0_pi_2_n_1_p_2",
"nE_e_1_mu_0_pi0_2_pi_2_n_2_p_1",
"nE_e_1_mu_0_pi0_2_pi_2_n_0_p_2",
"nE_e_1_mu_0_pi0_1_pi_2_n_0_p_2",
"nE_e_1_mu_0_pi0_1_pi_1_n_1_p_2",
"nT_e_0_mu_0_pi0_1_pi_1_n_0_p_1",
"nT_e_0_mu_0_pi0_0_pi_2_n_0_p_2",
"nM_e_0_mu_1_pi0_1_pi_2_n_2_p_1",
"nE_e_1_mu_0_pi0_2_pi_0_n_0_p_2",
"nE_e_1_mu_0_pi0_1_pi_0_n_2_p_1",
"nE_e_1_mu_0_pi0_0_pi_2_n_2_p_0",
"nE_e_1_mu_0_pi0_0_pi_0_n_2_p_0",
"nU_e_0_mu_0_pi0_2_pi_1_n_1_p_1",
"nT_e_0_mu_1_pi0_0_pi_1_n_1_p_0",
"nT_e_0_mu_0_pi0_0_pi_2_n_2_p_0",
"nM_e_0_mu_1_pi0_1_pi_2_n_1_p_2",
"nE_e_1_mu_0_pi0_2_pi_2_n_1_p_2",
"nE_e_1_mu_0_pi0_2_pi_1_n_1_p_0",
"nU_e_0_mu_0_pi0_2_pi_0_n_2_p_1",
"nU_e_0_mu_0_pi0_1_pi_2_n_0_p_0",
"nT_e_1_mu_0_pi0_0_pi_1_n_1_p_0",
"nT_e_1_mu_0_pi0_0_pi_1_n_0_p_1",
"nT_e_0_mu_0_pi0_2_pi_2_n_0_p_2",
"nT_e_0_mu_0_pi0_2_pi_1_n_0_p_1",
"nT_e_0_mu_0_pi0_1_pi_2_n_1_p_1",
"nM_e_0_mu_1_pi0_2_pi_2_n_0_p_2",
"nM_e_0_mu_1_pi0_2_pi_1_n_1_p_0",
"nE_e_1_mu_0_pi0_1_pi_2_n_1_p_1",
"nE_e_1_mu_0_pi0_0_pi_2_n_0_p_0",
"nU_e_0_mu_0_pi0_2_pi_1_n_2_p_0",
"nU_e_0_mu_0_pi0_2_pi_1_n_0_p_2",
"nU_e_0_mu_0_pi0_2_pi_0_n_2_p_0",
"nT_e_1_mu_0_pi0_1_pi_2_n_1_p_0",
"nT_e_1_mu_0_pi0_1_pi_2_n_0_p_1",
"nT_e_1_mu_0_pi0_0_pi_0_n_0_p_1",
"nT_e_0_mu_1_pi0_0_pi_2_n_2_p_2",
"nT_e_0_mu_1_pi0_0_pi_2_n_0_p_1",
"nT_e_0_mu_0_pi0_2_pi_2_n_1_p_1",
"nT_e_0_mu_0_pi0_1_pi_2_n_2_p_1",
"nT_e_0_mu_0_pi0_1_pi_2_n_1_p_2",
"nT_e_0_mu_0_pi0_0_pi_2_n_1_p_2",
"nT_e_0_mu_0_pi0_0_pi_2_n_0_p_0",
"nT_e_0_mu_0_pi0_0_pi_1_n_2_p_2",
"nM_e_0_mu_1_pi0_2_pi_1_n_0_p_2",
"nM_e_0_mu_1_pi0_1_pi_1_n_2_p_1",
"nM_e_0_mu_1_pi0_1_pi_1_n_2_p_0",
"nM_e_0_mu_1_pi0_1_pi_0_n_1_p_0",
"nM_e_0_mu_1_pi0_0_pi_2_n_0_p_0",
"nE_e_1_mu_0_pi0_2_pi_1_n_1_p_2",
"nE_e_1_mu_0_pi0_2_pi_1_n_1_p_1",
"nE_e_1_mu_0_pi0_2_pi_1_n_0_p_2",
"nE_e_1_mu_0_pi0_2_pi_0_n_1_p_1",
"nE_e_1_mu_0_pi0_0_pi_2_n_2_p_1",
"nU_e_0_mu_0_pi0_2_pi_2_n_0_p_0",
"nT_e_1_mu_0_pi0_0_pi_2_n_0_p_1",
"nT_e_1_mu_0_pi0_0_pi_0_n_1_p_1",
"nT_e_0_mu_0_pi0_2_pi_2_n_0_p_0",
"nM_e_0_mu_1_pi0_2_pi_2_n_2_p_0",
"nM_e_0_mu_1_pi0_2_pi_1_n_2_p_2",
"nM_e_0_mu_1_pi0_2_pi_1_n_1_p_2",
"nM_e_0_mu_1_pi0_2_pi_1_n_1_p_1",
"nM_e_0_mu_1_pi0_2_pi_0_n_2_p_2",
"nM_e_0_mu_1_pi0_2_pi_0_n_1_p_2",
"nE_e_1_mu_0_pi0_2_pi_2_n_2_p_0",
"nE_e_1_mu_0_pi0_2_pi_0_n_2_p_2",
"nE_e_1_mu_0_pi0_2_pi_0_n_2_p_1",
"nE_e_1_mu_0_pi0_2_pi_0_n_1_p_2",
"nE_e_1_mu_0_pi0_1_pi_2_n_2_p_0",
"nU_e_0_mu_0_pi0_2_pi_1_n_1_p_2",
"nU_e_0_mu_0_pi0_1_pi_1_n_0_p_0",
"nU_e_0_mu_0_pi0_0_pi_1_n_0_p_0",
"nT_e_1_mu_0_pi0_1_pi_2_n_1_p_2",
"nT_e_1_mu_0_pi0_0_pi_2_n_2_p_0",
"nT_e_1_mu_0_pi0_0_pi_1_n_2_p_2",
"nT_e_0_mu_1_pi0_2_pi_2_n_1_p_0",
"nT_e_0_mu_1_pi0_2_pi_2_n_0_p_1",
"nT_e_0_mu_1_pi0_1_pi_2_n_0_p_1",
"nT_e_0_mu_1_pi0_1_pi_0_n_0_p_1",
"nT_e_0_mu_1_pi0_0_pi_2_n_1_p_0",
"nT_e_0_mu_1_pi0_0_pi_2_n_0_p_2",
"nT_e_0_mu_1_pi0_0_pi_0_n_2_p_2",
"nT_e_0_mu_1_pi0_0_pi_0_n_0_p_2",
"nT_e_0_mu_0_pi0_2_pi_2_n_2_p_0",
"nT_e_0_mu_0_pi0_2_pi_2_n_1_p_2",
"nT_e_0_mu_0_pi0_2_pi_1_n_0_p_2",
"nT_e_0_mu_0_pi0_1_pi_2_n_0_p_2",
"nT_e_0_mu_0_pi0_1_pi_1_n_1_p_2",
"nT_e_0_mu_0_pi0_1_pi_1_n_1_p_1",
"nT_e_0_mu_0_pi0_1_pi_1_n_0_p_2",
"nT_e_0_mu_0_pi0_0_pi_1_n_1_p_2",
"nT_e_0_mu_0_pi0_0_pi_1_n_0_p_0",
"nM_e_0_mu_1_pi0_2_pi_2_n_2_p_1",
"nM_e_0_mu_1_pi0_2_pi_2_n_1_p_2",
"nM_e_0_mu_1_pi0_2_pi_0_n_1_p_1",
"nM_e_0_mu_1_pi0_1_pi_0_n_2_p_1",
"nM_e_0_mu_1_pi0_1_pi_0_n_0_p_0",
"nM_e_0_mu_1_pi0_0_pi_0_n_0_p_0",
"nE_e_1_mu_0_pi0_1_pi_2_n_0_p_0",
"nU_e_0_mu_0_pi0_2_pi_0_n_0_p_2",
"nT_e_1_mu_0_pi0_1_pi_0_n_0_p_1",
"nT_e_1_mu_0_pi0_0_pi_0_n_0_p_0",
"nT_e_0_mu_1_pi0_1_pi_2_n_2_p_2",
"nT_e_0_mu_1_pi0_1_pi_1_n_1_p_0",
"nT_e_0_mu_0_pi0_2_pi_2_n_2_p_1",
"nT_e_0_mu_0_pi0_1_pi_1_n_2_p_2",
"nT_e_0_mu_0_pi0_1_pi_1_n_1_p_0",
"nT_e_0_mu_0_pi0_0_pi_1_n_1_p_1",
"nT_e_0_mu_0_pi0_0_pi_1_n_0_p_2",
"nM_e_0_mu_1_pi0_2_pi_1_n_2_p_1",
"nM_e_0_mu_1_pi0_2_pi_1_n_2_p_0",
"nM_e_0_mu_1_pi0_1_pi_2_n_0_p_0",
"nM_e_0_mu_1_pi0_0_pi_0_n_2_p_0",
"nE_e_1_mu_0_pi0_2_pi_1_n_2_p_1",
"nE_e_1_mu_0_pi0_1_pi_1_n_2_p_0",
"nE_e_1_mu_0_pi0_1_pi_1_n_0_p_0",
"nE_e_1_mu_0_pi0_1_pi_0_n_0_p_0",
"nU_e_1_mu_0_pi0_0_pi_2_n_1_p_1",
"nU_e_0_mu_0_pi0_2_pi_1_n_0_p_0",
"nT_e_1_mu_0_pi0_2_pi_2_n_1_p_1",
"nT_e_1_mu_0_pi0_2_pi_2_n_1_p_0",
"nT_e_1_mu_0_pi0_2_pi_2_n_0_p_1",
"nT_e_1_mu_0_pi0_2_pi_1_n_1_p_1",
"nT_e_1_mu_0_pi0_1_pi_2_n_1_p_1",
"nT_e_1_mu_0_pi0_1_pi_2_n_0_p_2",
"nT_e_1_mu_0_pi0_1_pi_1_n_1_p_1",
"nT_e_1_mu_0_pi0_1_pi_1_n_0_p_1",
"nT_e_1_mu_0_pi0_0_pi_0_n_2_p_2",
"nT_e_0_mu_1_pi0_2_pi_2_n_2_p_2",
"nT_e_0_mu_1_pi0_2_pi_2_n_2_p_1",
"nT_e_0_mu_1_pi0_2_pi_0_n_1_p_2",
"nT_e_0_mu_1_pi0_1_pi_2_n_1_p_1",
"nT_e_0_mu_1_pi0_1_pi_2_n_1_p_0",
"nT_e_0_mu_1_pi0_1_pi_0_n_2_p_2",
"nT_e_0_mu_1_pi0_1_pi_0_n_2_p_1",
"nT_e_0_mu_1_pi0_0_pi_2_n_1_p_1",
"nT_e_0_mu_1_pi0_0_pi_1_n_2_p_2",
"nT_e_0_mu_1_pi0_0_pi_1_n_2_p_0",
"nT_e_0_mu_1_pi0_0_pi_1_n_1_p_2",
"nT_e_0_mu_1_pi0_0_pi_1_n_1_p_1",
"nT_e_0_mu_1_pi0_0_pi_0_n_1_p_1",
"nT_e_0_mu_0_pi0_2_pi_1_n_2_p_2",
"nT_e_0_mu_0_pi0_1_pi_2_n_2_p_0",
"nT_e_0_mu_0_pi0_1_pi_2_n_0_p_0",
"nT_e_0_mu_0_pi0_1_pi_1_n_2_p_0",
"nT_e_0_mu_0_pi0_0_pi_2_n_2_p_1",
"nT_e_0_mu_0_pi0_0_pi_1_n_2_p_1",
"nT_e_0_mu_0_pi0_0_pi_1_n_2_p_0",
"nT_e_0_mu_0_pi0_0_pi_1_n_1_p_0",
"nM_e_1_mu_1_pi0_0_pi_1_n_0_p_1",
"nM_e_0_mu_1_pi0_1_pi_1_n_0_p_0",
"nE_e_1_mu_0_pi0_2_pi_1_n_0_p_0",
"nU_e_1_mu_0_pi0_0_pi_0_n_0_p_0",
"nT_e_1_mu_0_pi0_2_pi_2_n_2_p_2",
"nT_e_1_mu_0_pi0_2_pi_2_n_0_p_2",
"nT_e_1_mu_0_pi0_2_pi_1_n_1_p_2",
"nT_e_1_mu_0_pi0_2_pi_1_n_1_p_0",
"nT_e_1_mu_0_pi0_2_pi_1_n_0_p_1",
"nT_e_1_mu_0_pi0_2_pi_1_n_0_p_0",
"nT_e_1_mu_0_pi0_2_pi_0_n_1_p_1",
"nT_e_1_mu_0_pi0_2_pi_0_n_0_p_1",
"nT_e_1_mu_0_pi0_1_pi_1_n_2_p_2",
"nT_e_1_mu_0_pi0_1_pi_1_n_1_p_2",
"nT_e_1_mu_0_pi0_1_pi_1_n_1_p_0",
"nT_e_1_mu_0_pi0_1_pi_0_n_2_p_2",
"nT_e_1_mu_0_pi0_1_pi_0_n_2_p_1",
"nT_e_1_mu_0_pi0_1_pi_0_n_1_p_1",
"nT_e_1_mu_0_pi0_1_pi_0_n_0_p_2",
"nT_e_1_mu_0_pi0_0_pi_2_n_2_p_1",
"nT_e_1_mu_0_pi0_0_pi_2_n_1_p_2",
"nT_e_1_mu_0_pi0_0_pi_2_n_1_p_0",
"nT_e_1_mu_0_pi0_0_pi_2_n_0_p_0",
"nT_e_1_mu_0_pi0_0_pi_1_n_2_p_1",
"nT_e_1_mu_0_pi0_0_pi_1_n_2_p_0",
"nT_e_1_mu_0_pi0_0_pi_1_n_1_p_1",
"nT_e_1_mu_0_pi0_0_pi_0_n_2_p_1",
"nT_e_1_mu_0_pi0_0_pi_0_n_1_p_2",
"nT_e_1_mu_0_pi0_0_pi_0_n_1_p_0",
"nT_e_1_mu_0_pi0_0_pi_0_n_0_p_2",
"nT_e_0_mu_1_pi0_2_pi_2_n_2_p_0",
"nT_e_0_mu_1_pi0_2_pi_2_n_1_p_2",
"nT_e_0_mu_1_pi0_2_pi_2_n_1_p_1",
"nT_e_0_mu_1_pi0_2_pi_2_n_0_p_2",
"nT_e_0_mu_1_pi0_2_pi_1_n_1_p_0",
"nT_e_0_mu_1_pi0_2_pi_1_n_0_p_1",
"nT_e_0_mu_1_pi0_2_pi_0_n_2_p_2",
"nT_e_0_mu_1_pi0_2_pi_0_n_0_p_1",
"nT_e_0_mu_1_pi0_1_pi_2_n_1_p_2",
"nT_e_0_mu_1_pi0_1_pi_2_n_0_p_2",
"nT_e_0_mu_1_pi0_1_pi_1_n_2_p_1",
"nT_e_0_mu_1_pi0_1_pi_1_n_0_p_1",
"nT_e_0_mu_1_pi0_0_pi_2_n_2_p_0",
"nT_e_0_mu_1_pi0_0_pi_2_n_1_p_2",
"nT_e_0_mu_1_pi0_0_pi_1_n_2_p_1",
"nT_e_0_mu_1_pi0_0_pi_1_n_0_p_2",
"nT_e_0_mu_1_pi0_0_pi_0_n_1_p_2",
"nT_e_0_mu_0_pi0_2_pi_1_n_1_p_2",
"nT_e_0_mu_0_pi0_2_pi_1_n_1_p_1",
"nT_e_0_mu_0_pi0_2_pi_1_n_0_p_0",
"nT_e_0_mu_0_pi0_2_pi_0_n_1_p_2",
"nT_e_0_mu_0_pi0_1_pi_1_n_2_p_1",
"nT_e_0_mu_0_pi0_1_pi_0_n_0_p_1",
"nT_e_0_mu_0_pi0_0_pi_0_n_0_p_1",
"nM_e_0_mu_1_pi0_2_pi_2_n_0_p_0",
"nM_e_0_mu_1_pi0_2_pi_0_n_2_p_1",
"nM_e_0_mu_1_pi0_2_pi_0_n_2_p_0",
"nM_e_0_mu_1_pi0_2_pi_0_n_0_p_0",
"nM_e_0_mu_1_pi0_1_pi_0_n_2_p_0",
"nE_e_1_mu_1_pi0_1_pi_0_n_1_p_0",
"nE_e_1_mu_0_pi0_2_pi_2_n_0_p_0",
"nE_e_1_mu_0_pi0_2_pi_1_n_2_p_0",
"nE_e_1_mu_0_pi0_2_pi_0_n_2_p_0",
"nE_e_1_mu_0_pi0_1_pi_0_n_2_p_0",
"nE_e_1_mu_0_pi0_1_pi_0_n_1_p_0",
"nM_Other",
"nE_Other",
"nT_Other",
"kCosmic_PT",
"kOther_PT"};

  int pos;
  if (std::find(v.begin(), v.end(), result) != v.end()){
    pos = std::find(v.begin(), v.end(), result) - v.begin();
    sType = ParticlesType(pos);
  }
  else{
    result = NuType + "_Other";
    if (std::find(v.begin(), v.end(), result) != v.end()){
      pos = std::find(v.begin(), v.end(), result) - v.begin();
      sType = ParticlesType(pos);
    }
  }

  return sType;
  }

std::ostream & cvn::operator<<	(	std::ostream &	os,
		const Boundary &	b
	)

Definition at line 49 of file Boundary.cxx.

References cvn::Boundary::FirstCell(), cvn::Boundary::FirstPlane(), cvn::Boundary::LastCell(), and cvn::Boundary::LastPlane().

  {
    os<<"Boundary with "
      <<"(first,last) plane: (" << b.FirstPlane()<<", "<< b.LastPlane()
      <<"), (first,last) cell even: ("<<b.FirstCell(0)<<", "<<b.LastCell(0)<<")"
      <<"), (first,last) cell odd: ("<<b.FirstCell(1)<<", "<<b.LastCell(1)<<")";

    return os;
  }

std::ostream& cvn::operator<<	(	std::ostream &	os,
		const PixelMapProducer &	p
	)

Definition at line 115 of file PixelMapProducer.cxx.

References cvn::PixelMapProducer::NCell(), and cvn::PixelMapProducer::NPlane().

  {
    os << "PixelMapProducer: "
       << p.NCell()  <<" cells X  " <<  p.NPlane() << " planes";
    return os;
  }

std::ostream& cvn::operator<<	(	std::ostream &	os,
		const RegPixelMapProducer &	p
	)

Definition at line 226 of file RegPixelMapProducer.cxx.

References cvn::RegPixelMapProducer::NCell(), and cvn::RegPixelMapProducer::NPlane().

  {
    os << "RegPixelMapProducer: "
       << p.NCell()  <<" cells X  " <<  p.NPlane() << " planes";
    return os;
  }

std::ostream & cvn::operator<<	(	std::ostream &	os,
		const RegPixelMap &	m
	)

Definition at line 298 of file RegPixelMap.cxx.

References cvn::RegPixelMap::NCell(), cvn::RegPixelMap::NPixel(), and cvn::RegPixelMap::NPlane().

  {
    os << "RegPixelMap with " << m.NPixel() << " pixels, "
                           << m.NPlane() << " planes"
                << " by "  << m.NCell()  << " cells" ;
    return os;
  }

std::ostream & cvn::operator<<	(	std::ostream &	os,
		const PixelMap &	m
	)

Definition at line 387 of file PixelMap.cxx.

References cvn::PixelMap::NCell(), cvn::PixelMap::NPixel(), and cvn::PixelMap::NPlane().

  {
    os << "PixelMap with " << m.NPixel() << " pixels, "
                           << m.NPlane() << " planes"
                << " by "  << m.NCell()  << " cells" ;
    return os;
  }

ParticlesType cvn::ParticlesSliceClassify	(	const art::Ptr< simb::MCTruth >	truth,
		const rb::Cluster &	slice,
		float &	nuEnergy,
		float &	lepEnergy,
		double &	vtxx,
		double &	vtxy,
		double &	vtxz
	)

Definition at line 753 of file AssignLabels.cxx.

References abs(), simb::MCNeutrino::CCNC(), simb::MCParticle::E(), simb::MCTruth::GetNeutrino(), simb::kCC, kCosmic_PT, simb::MCNeutrino::Lepton(), MakeParticleLabel(), simb::MCTruth::NeutrinoSet(), simb::MCNeutrino::Nu(), simb::MCParticle::PdgCode(), string, simb::MCParticle::Vx(), simb::MCParticle::Vy(), and simb::MCParticle::Vz().

Referenced by cvn::CVNEventDump::analyze(), cvn::CVNCosmicEventDump::analyze(), and cvn::CVNAddTrainingData::produce().

                                                                                {

    //std::cout << "I am here" << std::endl;
  if(!truth->NeutrinoSet()) return kCosmic_PT;

  // get nu energy
  const simb::MCNeutrino& nu = truth->GetNeutrino();

  std::string NuType;
  bool nuIsCC = nu.CCNC() == simb::kCC;
  if (nuIsCC){

    if (abs(nu.Nu().PdgCode()) == 14){
      NuType = "nM";      
    }
    else if (abs(nu.Nu().PdgCode()) == 12){
      NuType = "nE";
    }
    else if (abs(nu.Nu().PdgCode()) == 16){ 
      NuType = "nT";
    }

  } // close CC if check

  else{
    NuType = "nU";
  } // close NC if check

  nuEnergy = nu.Nu().E();
  lepEnergy = nu.Lepton().E();
  vtxx = nu.Nu().Vx();
  vtxy = nu.Nu().Vy();
  vtxz = nu.Nu().Vz();
  return MakeParticleLabel(truth, NuType);
}

ProngType cvn::ProngClassify	(	const rb::Prong &	prong,
		ProngType *	pType3D,
		ProngType *	pTypeX,
		ProngType *	pTypeY,
		bool *	isprimary,
		double *	purity3D,
		double *	purityX,
		double *	purityY,
		double *	recE,
		unsigned int *	ncellX,
		unsigned int *	ncellY
	)

Definition at line 1101 of file AssignLabels.cxx.

References rb::Cluster::AllCells(), GetTypeByPDG(), cheat::BackTracker::HitCollectionPurity(), cheat::BackTracker::HitsToTrackIDE(), sim::ParticleNavigator::IsPrimary(), kEmpty, kUnknown, rb::Cluster::NXCell(), rb::Cluster::NYCell(), cheat::BackTracker::ParticleNavigator(), simb::MCParticle::PdgCode(), rb::Cluster::TotalGeV(), cheat::BackTracker::TrackIDToParticle(), rb::Cluster::XCells(), and rb::Cluster::YCells().

Referenced by cvn::CVNEventProngDump::analyze(), bpf::BPFCVNAna::analyze(), and cvn::CVNAddTrainingData::produce().

  {
    
    *isprimary = false;
    
    *pType3D = kUnknown;
    *pTypeX  = kUnknown;
    *pTypeY  = kUnknown;
    
    *purity3D = 0.0;
    *purityX  = 0.0;
    *purityY  = 0.0;
    *purityY  = 0.0;
    
    *ncellX = prong.NXCell();
    *ncellY = prong.NYCell();
    
    *recE = prong.TotalGeV();
    
    if( prong.NXCell() == 0 || 
        prong.NYCell() == 0) *pType3D = kEmpty;
    if( prong.NXCell() == 0 ) *pTypeX  = kEmpty;
    if( prong.NYCell() == 0 ) *pTypeY  = kEmpty;
    
    art::ServiceHandle<cheat::BackTracker> bt;
    
    const sim::ParticleNavigator& pnav = bt->ParticleNavigator();
    
    std::vector<cheat::TrackIDE> ParticleTrackIDs, ParticleTrackIDsX,
      ParticleTrackIDsY;
    
    if( *pType3D != kEmpty )
      ParticleTrackIDs  = bt->HitsToTrackIDE(prong.AllCells());
    if( *pTypeX  != kEmpty )
      ParticleTrackIDsX = bt->HitsToTrackIDE(prong.XCells());
    if( *pTypeY != kEmpty )
      ParticleTrackIDsY = bt->HitsToTrackIDE(prong.YCells());
    
    if( ParticleTrackIDs.size()  < 1 ) *pType3D = kEmpty;
    if( ParticleTrackIDsX.size() < 1 ) *pTypeX  = kEmpty;
    if( ParticleTrackIDsY.size() < 1 ) *pTypeY  = kEmpty;
    
    if( *pType3D != kEmpty ){
      double totalPur   = 0.0;
      double highestPur = 0.0;
      //int highestPurID;
      int highestPurPDG;
      for (unsigned int iID=0; iID < ParticleTrackIDs.size(); ++iID){
        std::set<int> TrackID;
        std::map<int, double> purMap;
        std::map<int, int> parents;
        
        int iPartID = ParticleTrackIDs[iID].trackID;
        TrackID.insert(iPartID);
        double iPartPur = bt->HitCollectionPurity(TrackID, prong.AllCells(),
                                                  &purMap, &parents,true);
        
        if( iPartPur > highestPur ){
          highestPur    = iPartPur;
          //highestPurID  = iPartID;
          highestPurPDG = bt->TrackIDToParticle(iPartID)->PdgCode();
          *pType3D   = GetTypeByPDG(highestPurPDG);
          *purity3D  = highestPur;
          *isprimary = pnav.IsPrimary(iPartID);
        }
        
        totalPur+=iPartPur;
      }
      
    }
    
    if( *pTypeX  != kEmpty ){
      double totalPurX   = 0.0;
      double highestPurX = 0.0;
      //      int highestPurIDX;
      int  highestPurPDGx;
      for (unsigned int ixID=0; ixID < ParticleTrackIDsX.size(); ++ixID){
        std::set<int> TrackID;
        std::map<int, double> purMap;
        std::map<int, int> parents;
        
        int ixPartID = ParticleTrackIDsX[ixID].trackID;
        TrackID.insert(ixPartID);
        double ixPartPur = bt->HitCollectionPurity(TrackID, prong.XCells(),
                                                   &purMap, &parents,true);
        
        if( ixPartPur > highestPurX ){
          highestPurX   = ixPartPur;
          //highestPurIDX = ixPartID;
          highestPurPDGx = bt->TrackIDToParticle(ixPartID)->PdgCode();
          *pTypeX  = GetTypeByPDG(highestPurPDGx);
          *purityX = highestPurX;
        }
        
        totalPurX+=ixPartPur;
      }
      
    }
    if( *pTypeY != kEmpty ){
      double totalPurY   = 0.0;
      double highestPurY = 0.0;
      //      int highestPurIDY;
      int highestPurPDGy;
      for (unsigned int iyID=0; iyID < ParticleTrackIDsY.size(); ++iyID){
        std::set<int> TrackID;
        std::map<int, double> purMap;
        std::map<int, int> parents;
        
        int iyPartID = ParticleTrackIDsY[iyID].trackID;
        TrackID.insert(iyPartID);
        double iyPartPur = bt->HitCollectionPurity(TrackID, prong.YCells(), &purMap,
                                                   &parents,true);
        
        if( iyPartPur > highestPurY ){
          highestPurY   = iyPartPur;
          //highestPurIDY = iyPartID;
          highestPurPDGy = bt->TrackIDToParticle(iyPartID)->PdgCode();
          *pTypeY = GetTypeByPDG(highestPurPDGy);
          *purityY = highestPurY;
        }
        
        totalPurY+=iyPartPur;
      }
      
    }
    
    ProngType pType;
    
    pType = *pType3D;
    
    return pType;
  }// ProngClassify.

InteractionType cvn::SliceClassify	(	const art::Ptr< simb::MCTruth >	truth,
		const rb::Cluster &	slice,
		float &	nuEnergy,
		float &	lepEnergy,
		double &	vtxx,
		double &	vtxy,
		double &	vtxz
	)

Definition at line 275 of file AssignLabels.cxx.

References simb::MCNeutrino::CCNC(), simb::MCParticle::E(), GetInteractionTypeFromSlice(), simb::MCTruth::GetNeutrino(), simb::kCC, kOther, simb::MCNeutrino::Lepton(), simb::MCNeutrino::Mode(), simb::MCTruth::NeutrinoSet(), simb::MCNeutrino::Nu(), simb::MCParticle::PdgCode(), simb::MCParticle::Vx(), simb::MCParticle::Vy(), and simb::MCParticle::Vz().

Referenced by cvn::CVNEventDump::analyze(), cvn::CVNCosmicEventDump::analyze(), cvn::CVNEventProngDump::analyze(), and cvn::CVNAddTrainingData::produce().

                                                                         {
    
    //std::cout << "Slice NCells = " << slice.NCell() << std::endl;
    
    if(!truth->NeutrinoSet()) return kOther;
    
    // get nu energy
    const simb::MCNeutrino& nu = truth->GetNeutrino();
    int  nuPDG       = nu.Nu().PdgCode();
    bool nuIsCC      = nu.CCNC() == simb::kCC;
    int  trueMode    = nu.Mode();
    
    InteractionType sType = kOther;
    
    sType = GetInteractionTypeFromSlice(nuPDG, nuIsCC, trueMode);
    
    nuEnergy = nu.Nu().E();
    lepEnergy = nu.Lepton().E();
    vtxx = nu.Nu().Vx();
    vtxy = nu.Nu().Vy();
    vtxz = nu.Nu().Vz();