Classes | Typedefs | Enumerations | Functions | Variables
caf Namespace Reference

This module creates Common Analysis Files. More...

Classes

class  ArrayVectorProxyBase
 
class  CAFMaker
 Module to create Common Analysis Files from ART files. More...
 
struct  CAFMakerParams
 
struct  CVNFinalState
 
class  Proxy
 
class  Proxy< caf::SRBeam >
 Proxy for caf::SRBeam. More...
 
class  Proxy< caf::SRBpf >
 Proxy for caf::SRBpf. More...
 
class  Proxy< caf::SRBPFEnergy >
 Proxy for caf::SRBPFEnergy. More...
 
class  Proxy< caf::SRBpfId >
 Proxy for caf::SRBpfId. More...
 
class  Proxy< caf::SRBpfTrack >
 Proxy for caf::SRBpfTrack. More...
 
class  Proxy< caf::SRContain >
 Proxy for caf::SRContain. More...
 
class  Proxy< caf::SRCosmic >
 Proxy for caf::SRCosmic. More...
 
class  Proxy< caf::SRCosmicCVN >
 Proxy for caf::SRCosmicCVN. More...
 
class  Proxy< caf::SRCosRej >
 Proxy for caf::SRCosRej. More...
 
class  Proxy< caf::SRCVNFeatures >
 Proxy for caf::SRCVNFeatures. More...
 
class  Proxy< caf::SRCVNNeutronDaughterResult >
 Proxy for caf::SRCVNNeutronDaughterResult. More...
 
class  Proxy< caf::SRCVNParticleResult >
 Proxy for caf::SRCVNParticleResult. More...
 
class  Proxy< caf::SRCVNResult >
 Proxy for caf::SRCVNResult. More...
 
class  Proxy< caf::SRElastic >
 Proxy for caf::SRElastic. More...
 
class  Proxy< caf::SRELid >
 Proxy for caf::SRELid. More...
 
class  Proxy< caf::SREnergyBranch >
 Proxy for caf::SREnergyBranch. More...
 
class  Proxy< caf::SRFluxWeights >
 Proxy for caf::SRFluxWeights. More...
 
class  Proxy< caf::SRFuzzyK >
 Proxy for caf::SRFuzzyK. More...
 
class  Proxy< caf::SRFuzzyKProng >
 Proxy for caf::SRFuzzyKProng. More...
 
class  Proxy< caf::SRGeant4Weights >
 Proxy for caf::SRGeant4Weights. More...
 
class  Proxy< caf::SRGenieWeights >
 Proxy for caf::SRGenieWeights. More...
 
class  Proxy< caf::SRGlobalTruth >
 Proxy for caf::SRGlobalTruth. More...
 
class  Proxy< caf::SRHadClust >
 Proxy for caf::SRHadClust. More...
 
class  Proxy< caf::SRHeader >
 Proxy for caf::SRHeader. More...
 
class  Proxy< caf::SRHoughVertex >
 Proxy for caf::SRHoughVertex. More...
 
class  Proxy< caf::SRIDBranch >
 Proxy for caf::SRIDBranch. More...
 
class  Proxy< caf::SRKalman >
 Proxy for caf::SRKalman. More...
 
class  Proxy< caf::SRKalmanTrack >
 Proxy for caf::SRKalmanTrack. More...
 
class  Proxy< caf::SRLem >
 Proxy for caf::SRLem. More...
 
class  Proxy< caf::SRLorentzVector >
 Proxy for caf::SRLorentzVector. More...
 
class  Proxy< caf::SRMCReweight >
 Proxy for caf::SRMCReweight. More...
 
class  Proxy< caf::SRMichelE >
 Proxy for caf::SRMichelE. More...
 
class  Proxy< caf::SRMRCCParent >
 Proxy for caf::SRMRCCParent. More...
 
class  Proxy< caf::SRMRProperties >
 Proxy for caf::SRMRProperties. More...
 
class  Proxy< caf::SRMuonID >
 Proxy for caf::SRMuonID. More...
 
class  Proxy< caf::SRNCCosRej >
 Proxy for caf::SRNCCosRej. More...
 
class  Proxy< caf::SRNCPi0BkgRej >
 Proxy for caf::SRNCPi0BkgRej. More...
 
class  Proxy< caf::SRNeutrino >
 Proxy for caf::SRNeutrino. More...
 
class  Proxy< caf::SRNueCosRej >
 Proxy for caf::SRNueCosRej. More...
 
class  Proxy< caf::SRNueEnergy >
 Proxy for caf::SRNueEnergy. More...
 
class  Proxy< caf::SRNuGenParticle >
 Proxy for caf::SRNuGenParticle. More...
 
class  Proxy< caf::SRNumuEnergy >
 Proxy for caf::SRNumuEnergy. More...
 
class  Proxy< caf::SRNuonEResult >
 Proxy for caf::SRNuonEResult. More...
 
class  Proxy< caf::SRParentBranch >
 Proxy for caf::SRParentBranch. More...
 
class  Proxy< caf::SRParticleTruth >
 Proxy for caf::SRParticleTruth. More...
 
class  Proxy< caf::SRPixelMap >
 Proxy for caf::SRPixelMap. More...
 
class  Proxy< caf::SRPixelObjMap >
 Proxy for caf::SRPixelObjMap. More...
 
class  Proxy< caf::SRPresel >
 Proxy for caf::SRPresel. More...
 
class  Proxy< caf::SRProng >
 Proxy for caf::SRProng. More...
 
class  Proxy< caf::SRProngTrainingData >
 Proxy for caf::SRProngTrainingData. More...
 
class  Proxy< caf::SRRegCVNResult >
 Proxy for caf::SRRegCVNResult. More...
 
class  Proxy< caf::SRRemid >
 Proxy for caf::SRRemid. More...
 
class  Proxy< caf::SRRvp >
 Proxy for caf::SRRvp. More...
 
class  Proxy< caf::SRShower >
 Proxy for caf::SRShower. More...
 
class  Proxy< caf::SRShowerLID >
 Proxy for caf::SRShowerLID. More...
 
class  Proxy< caf::SRSlcME >
 Proxy for caf::SRSlcME. More...
 
class  Proxy< caf::SRSlice >
 Proxy for caf::SRSlice. More...
 
class  Proxy< caf::SRSliceLID >
 Proxy for caf::SRSliceLID. More...
 
class  Proxy< caf::SRSliceMap >
 Proxy for caf::SRSliceMap. More...
 
class  Proxy< caf::SRSLid >
 Proxy for caf::SRSLid. More...
 
class  Proxy< caf::SRSLidEnergy >
 Proxy for caf::SRSLidEnergy. More...
 
class  Proxy< caf::SRSpill >
 Proxy for caf::SRSpill. More...
 
class  Proxy< caf::SRTrack >
 Proxy for caf::SRTrack. More...
 
class  Proxy< caf::SRTrackBase >
 Proxy for caf::SRTrackBase. More...
 
class  Proxy< caf::SRTrackBranch >
 Proxy for caf::SRTrackBranch. More...
 
class  Proxy< caf::SRTrainingBranch >
 Proxy for caf::SRTrainingBranch. More...
 
class  Proxy< caf::SRTrainingData >
 Proxy for caf::SRTrainingData. More...
 
class  Proxy< caf::SRTrkME >
 Proxy for caf::SRTrkME. More...
 
class  Proxy< caf::SRTrueMichelE >
 Proxy for caf::SRTrueMichelE. More...
 
class  Proxy< caf::SRTrueNumuEnergy >
 Proxy for caf::SRTrueNumuEnergy. More...
 
class  Proxy< caf::SRTrueParticle >
 Proxy for caf::SRTrueParticle. More...
 
class  Proxy< caf::SRTruthBranch >
 Proxy for caf::SRTruthBranch. More...
 
class  Proxy< caf::SRVector3D >
 Proxy for caf::SRVector3D. More...
 
class  Proxy< caf::SRVertexBranch >
 Proxy for caf::SRVertexBranch. More...
 
class  Proxy< caf::SRVertexDT >
 Proxy for caf::SRVertexDT. More...
 
class  Proxy< caf::SRVeto >
 Proxy for caf::SRVeto. More...
 
class  Proxy< caf::SRXnue >
 Proxy for caf::SRXnue. More...
 
class  Proxy< caf::StandardRecord >
 Proxy for caf::StandardRecord. More...
 
class  Proxy< std::vector< T > >
 
class  Proxy< T[N]>
 
class  Restorer
 
class  RestorerT
 
class  SRBeam
 Information about the neutrino production. Docs from http://www.hep.utexas.edu/~zarko/wwwgnumi/v19/v19/output_gnumi.html. More...
 
class  SRBpf
 
class  SRBPFEnergy
 BPF energy estimator output. More...
 
class  SRBpfId
 Breakpoint ID (BpfId) output. More...
 
class  SRBpfTrack
 
class  SRBranchRegistry
 
class  SRContain
 Containment variables. More...
 
class  SRCosmic
 Truth information for cosmic rays. More...
 
class  SRCosmicCVN
 
class  SRCosRej
 Output from Cosmic Rejection (CosRej) module. More...
 
class  SRCVNFeatures
 CVN features. More...
 
class  SRCVNNeutronDaughterResult
 CVN PID output for a single particle. More...
 
class  SRCVNParticleResult
 CVN PID output for a single particle. More...
 
class  SRCVNResult
 CVN PID output. More...
 
class  SRElastic
 A potential interaction point from the ElasticArms algorithm. More...
 
class  SRELid
 Output of the LIDBuilder module (slid::lid objects). More...
 
class  SREnergy
 
class  SREnergyBranch
 An SREnergyBranch contains vectors of energy objects. More...
 
class  SRFluxWeights
 Reweight information for flux systematic. More...
 
class  SRFuzzyK
 
class  SRFuzzyKProng
 
class  SRGeant4Weights
 Reweight information for geant4 systematic. More...
 
class  SRGenieWeights
 Reweight information for a single GENIE systematic. More...
 
class  SRGlobalTruth
 
class  SRHadClust
 Overarching information for a numu hadronic cluster. More...
 
class  SRHeader
 Header representing overview information for the current event/slice. More...
 
class  SRHoughVertex
 A potential interaction point found by the HoughVertex algorithm. More...
 
class  SRIDBranch
 Event ID and selection variables. More...
 
class  SRJMEid
 Output of the jmshower::NueSel module. More...
 
class  SRJMEnergy
 DO NOT USE. ONLY HERE FOR BACKWARDS COMPATIBILITY. More...
 
class  SRJMShower
 A reconstructed shower from the JMShower module. More...
 
class  SRKalman
 
class  SRKalmanTrack
 
class  SRLem
 This class contains the LEM PID output. More...
 
class  SRLid
 DO NOT USE. ONLY HERE FOR BACKWARDS COMPATIBILITY. More...
 
class  SRLorentzVector
 4-vector with more efficient storage than TLorentzVector More...
 
class  SRMCReweight
 Various weights for systematic reweights of MC. More...
 
class  SRMichelE
 
class  SRMRCCParent
 An SRMRCCParent holds information about the slice that was parent to the current slice. It is currently being used by Muon Removed Charged Current Analysis. More...
 
class  SRMRProperties
 A reconstructed shower from the MRProperties module. More...
 
class  SRMuId
 This class contains the LLH muon PID output. More...
 
class  SRMuonID
 Contains the reco muon PID (ReMId) output. More...
 
class  SRNCCosRej
 Output from Cosmic Rejection (Nuecosrej) module. More...
 
class  SRNCPi0BkgRej
 
class  SRNDSandbox
 
class  SRNeutrino
 The SRNeutrino is a representation of neutrino interaction information. More...
 
class  SRNueCosRej
 Output from Cosmic Rejection (Nuecosrej) module. More...
 
class  SRNueEnergy
 Nue energy estimator output in the standard record. More...
 
class  SRNuePresel
 Nue preselection information. More...
 
class  SRNueSandbox
 
class  SRNuGenParticle
 SRNuGenParticle represents a particle within the simulation of a neutrino interaction. More...
 
class  SRNumuEnergy
 Numu energy estimator output. More...
 
class  SRNumuSandbox
 
class  SRNuonEResult
 NuonE output. More...
 
class  SRNusSandbox
 
class  SRParentBranch
 
class  SRParticleTruth
 The truth information of reco objects within a slice. More...
 
class  SRPixelMap
 Variables describing Michel E's found around the end of a track. More...
 
class  SRPixelObjMap
 Variables describing Michel E's found around the end of a track. More...
 
class  SRPresel
 preselection information More...
 
class  SRProng
 
class  SRProngTrainingData
 
class  SRProngXSec
 
class  SRProxySystController
 
class  SRQepid
 Contains the quasielastic muon PID (QePId) output. More...
 
class  SRRegCVNResult
 Regression CVN output. More...
 
class  SRRemid
 Contains the reco muon PID (ReMId) output. More...
 
class  SRRvp
 Contains the RVP PID output. More...
 
class  SRSandbox
 
class  SRShower
 
class  SRShowerBranch
 
class  SRShowerLID
 
class  SRShowerPID
 DO NOT USE. ONLY HERE FOR BACKWARDS COMPATIBILITY. More...
 
class  SRSlcME
 Represents output from SlcMEFilter. More...
 
class  SRSlice
 An SRSlice contains overarching information for a slice. More...
 
class  SRSliceLID
 
class  SRSliceMap
 Variables describing Michel E's found around the end of a track. More...
 
class  SRSLid
 This class contains the LID pid information for a shower (slid::ShowerLID objects). More...
 
class  SRSLidEnergy
 This is a class for the NueSel energy estimate. More...
 
class  SRSPid
 Contains the SPID pid information for a shower (slid::ShowerLID objects). More...
 
class  SRSpill
 
class  SRSpillTruthBranch
 Truth info for all neutrinos in the spill. More...
 
class  SRTrack
 
class  SRTrackBase
 
class  SRTrackBranch
 Reconstructed tracks found by various algorithms. More...
 
class  SRTrainingBranch
 Event ID training variables. More...
 
class  SRTrainingData
 
class  SRTrkME
 Variables describing Michel E's found around the end of a track. More...
 
class  SRTrueMichelE
 Truth information for a Michel electron. More...
 
class  SRTrueNumuEnergy
 Truth information for numu energy fitting. More...
 
class  SRTrueParticle
 
class  SRTruth
 Represents a true neutrino. More...
 
class  SRTruthBranch
 Contains truth information for the slice for the parent neutrino/cosmic. More...
 
class  SRVector3D
 A 3-vector with more efficient storage than TVector3. More...
 
class  SRVertex
 Time and position of a reconstructed vertex. More...
 
class  SRVertexBranch
 Vectors of reconstructed vertices found by various algorithms. More...
 
class  SRVertexDT
 A vertex found by the VertexDT algorithm. More...
 
class  SRVeto
 Details of processing cuts made by the veto modules. More...
 
class  SRXnue
 Store BDT variables for the short-baseline oscillation study. More...
 
class  StandardRecord
 The StandardRecord is the primary top-level object in the Common Analysis File trees. More...
 
class  VectorProxyBase
 

Typedefs

using FlatSRSLidEnergy = flat::Flat< caf::SRSLidEnergy >
 
using FlatSRNueEnergy = flat::Flat< caf::SRNueEnergy >
 
using FlatSRBPFEnergy = flat::Flat< caf::SRBPFEnergy >
 
using FlatSRVector3D = flat::Flat< caf::SRVector3D >
 
using FlatSRHadClust = flat::Flat< caf::SRHadClust >
 
using FlatSRTrueNumuEnergy = flat::Flat< caf::SRTrueNumuEnergy >
 
using FlatSRNumuEnergy = flat::Flat< caf::SRNumuEnergy >
 
using FlatSREnergyBranch = flat::Flat< caf::SREnergyBranch >
 
using FlatSRHeader = flat::Flat< caf::SRHeader >
 
using FlatSRLorentzVector = flat::Flat< caf::SRLorentzVector >
 
using FlatSRTrueMichelE = flat::Flat< caf::SRTrueMichelE >
 
using FlatSRCosmic = flat::Flat< caf::SRCosmic >
 
using FlatSRBeam = flat::Flat< caf::SRBeam >
 
using FlatSRNuGenParticle = flat::Flat< caf::SRNuGenParticle >
 
using FlatSRTrueParticle = flat::Flat< caf::SRTrueParticle >
 
using FlatSRGeant4Weights = flat::Flat< caf::SRGeant4Weights >
 
using FlatSRGenieWeights = flat::Flat< caf::SRGenieWeights >
 
using FlatSRFluxWeights = flat::Flat< caf::SRFluxWeights >
 
using FlatSRMCReweight = flat::Flat< caf::SRMCReweight >
 
using FlatSRNeutrino = flat::Flat< caf::SRNeutrino >
 
using FlatSRGlobalTruth = flat::Flat< caf::SRGlobalTruth >
 
using FlatSRTruthBranch = flat::Flat< caf::SRTruthBranch >
 
using FlatSRParticleTruth = flat::Flat< caf::SRParticleTruth >
 
using FlatSRSlcME = flat::Flat< caf::SRSlcME >
 
using FlatSRTrkME = flat::Flat< caf::SRTrkME >
 
using FlatSRMichelE = flat::Flat< caf::SRMichelE >
 
using FlatSRMRCCParent = flat::Flat< caf::SRMRCCParent >
 
using FlatSRParentBranch = flat::Flat< caf::SRParentBranch >
 
using FlatSRBpfId = flat::Flat< caf::SRBpfId >
 
using FlatSRContain = flat::Flat< caf::SRContain >
 
using FlatSRCosRej = flat::Flat< caf::SRCosRej >
 
using FlatSRCVNResult = flat::Flat< caf::SRCVNResult >
 
using FlatSRLem = flat::Flat< caf::SRLem >
 
using FlatSRELid = flat::Flat< caf::SRELid >
 
using FlatSRMuonID = flat::Flat< caf::SRMuonID >
 
using FlatSRNCCosRej = flat::Flat< caf::SRNCCosRej >
 
using FlatSRNCPi0BkgRej = flat::Flat< caf::SRNCPi0BkgRej >
 
using FlatSRNueCosRej = flat::Flat< caf::SRNueCosRej >
 
using FlatSRPresel = flat::Flat< caf::SRPresel >
 
using FlatSRNuonEResult = flat::Flat< caf::SRNuonEResult >
 
using FlatSRRemid = flat::Flat< caf::SRRemid >
 
using FlatSRRvp = flat::Flat< caf::SRRvp >
 
using FlatSRSliceLID = flat::Flat< caf::SRSliceLID >
 
using FlatSRVeto = flat::Flat< caf::SRVeto >
 
using FlatSRXnue = flat::Flat< caf::SRXnue >
 
using FlatSRIDBranch = flat::Flat< caf::SRIDBranch >
 
using FlatSRSlice = flat::Flat< caf::SRSlice >
 
using FlatSRCosmicCVN = flat::Flat< caf::SRCosmicCVN >
 
using FlatSRSpill = flat::Flat< caf::SRSpill >
 
using FlatSRCVNFeatures = flat::Flat< caf::SRCVNFeatures >
 
using FlatSRPixelObjMap = flat::Flat< caf::SRPixelObjMap >
 
using FlatSRSliceMap = flat::Flat< caf::SRSliceMap >
 
using FlatSRTrainingData = flat::Flat< caf::SRTrainingData >
 
using FlatSRTrainingBranch = flat::Flat< caf::SRTrainingBranch >
 
using FlatSRShower = flat::Flat< caf::SRShower >
 
using FlatSRSLid = flat::Flat< caf::SRSLid >
 
using FlatSRShowerLID = flat::Flat< caf::SRShowerLID >
 
using FlatSRMRProperties = flat::Flat< caf::SRMRProperties >
 
using FlatSRTrack = flat::Flat< caf::SRTrack >
 
using FlatSRTrackBase = flat::Flat< caf::SRTrackBase >
 
using FlatSRKalmanTrack = flat::Flat< caf::SRKalmanTrack >
 
using FlatSRKalman = flat::Flat< caf::SRKalman >
 
using FlatSRTrackBranch = flat::Flat< caf::SRTrackBranch >
 
using FlatSRPixelMap = flat::Flat< caf::SRPixelMap >
 
using FlatSRCVNNeutronDaughterResult = flat::Flat< caf::SRCVNNeutronDaughterResult >
 
using FlatSRCVNParticleResult = flat::Flat< caf::SRCVNParticleResult >
 
using FlatSRProngTrainingData = flat::Flat< caf::SRProngTrainingData >
 
using FlatSRProng = flat::Flat< caf::SRProng >
 
using FlatSRBpfTrack = flat::Flat< caf::SRBpfTrack >
 
using FlatSRBpf = flat::Flat< caf::SRBpf >
 
using FlatSRRegCVNResult = flat::Flat< caf::SRRegCVNResult >
 
using FlatSRFuzzyKProng = flat::Flat< caf::SRFuzzyKProng >
 
using FlatSRFuzzyK = flat::Flat< caf::SRFuzzyK >
 
using FlatSRElastic = flat::Flat< caf::SRElastic >
 
using FlatSRHoughVertex = flat::Flat< caf::SRHoughVertex >
 
using FlatSRVertexDT = flat::Flat< caf::SRVertexDT >
 
using FlatSRVertexBranch = flat::Flat< caf::SRVertexBranch >
 
using FlatStandardRecord = flat::Flat< caf::StandardRecord >
 
using SRProxy = caf::Proxy< caf::StandardRecord >
 
using SRSLidEnergyProxy = caf::Proxy< caf::SRSLidEnergy >
 
using SRNueEnergyProxy = caf::Proxy< caf::SRNueEnergy >
 
using SRBPFEnergyProxy = caf::Proxy< caf::SRBPFEnergy >
 
using SRVector3DProxy = caf::Proxy< caf::SRVector3D >
 
using SRHadClustProxy = caf::Proxy< caf::SRHadClust >
 
using SRTrueNumuEnergyProxy = caf::Proxy< caf::SRTrueNumuEnergy >
 
using SRNumuEnergyProxy = caf::Proxy< caf::SRNumuEnergy >
 
using SREnergyBranchProxy = caf::Proxy< caf::SREnergyBranch >
 
using SRHeaderProxy = caf::Proxy< caf::SRHeader >
 
using SRLorentzVectorProxy = caf::Proxy< caf::SRLorentzVector >
 
using SRTrueMichelEProxy = caf::Proxy< caf::SRTrueMichelE >
 
using SRCosmicProxy = caf::Proxy< caf::SRCosmic >
 
using SRBeamProxy = caf::Proxy< caf::SRBeam >
 
using SRNuGenParticleProxy = caf::Proxy< caf::SRNuGenParticle >
 
using SRTrueParticleProxy = caf::Proxy< caf::SRTrueParticle >
 
using SRGeant4WeightsProxy = caf::Proxy< caf::SRGeant4Weights >
 
using SRGenieWeightsProxy = caf::Proxy< caf::SRGenieWeights >
 
using SRFluxWeightsProxy = caf::Proxy< caf::SRFluxWeights >
 
using SRMCReweightProxy = caf::Proxy< caf::SRMCReweight >
 
using SRNeutrinoProxy = caf::Proxy< caf::SRNeutrino >
 
using SRGlobalTruthProxy = caf::Proxy< caf::SRGlobalTruth >
 
using SRTruthBranchProxy = caf::Proxy< caf::SRTruthBranch >
 
using SRParticleTruthProxy = caf::Proxy< caf::SRParticleTruth >
 
using SRSlcMEProxy = caf::Proxy< caf::SRSlcME >
 
using SRTrkMEProxy = caf::Proxy< caf::SRTrkME >
 
using SRMichelEProxy = caf::Proxy< caf::SRMichelE >
 
using SRMRCCParentProxy = caf::Proxy< caf::SRMRCCParent >
 
using SRParentBranchProxy = caf::Proxy< caf::SRParentBranch >
 
using SRBpfIdProxy = caf::Proxy< caf::SRBpfId >
 
using SRContainProxy = caf::Proxy< caf::SRContain >
 
using SRCosRejProxy = caf::Proxy< caf::SRCosRej >
 
using SRCVNResultProxy = caf::Proxy< caf::SRCVNResult >
 
using SRLemProxy = caf::Proxy< caf::SRLem >
 
using SRELidProxy = caf::Proxy< caf::SRELid >
 
using SRMuonIDProxy = caf::Proxy< caf::SRMuonID >
 
using SRNCCosRejProxy = caf::Proxy< caf::SRNCCosRej >
 
using SRNCPi0BkgRejProxy = caf::Proxy< caf::SRNCPi0BkgRej >
 
using SRNueCosRejProxy = caf::Proxy< caf::SRNueCosRej >
 
using SRPreselProxy = caf::Proxy< caf::SRPresel >
 
using SRNuonEResultProxy = caf::Proxy< caf::SRNuonEResult >
 
using SRRemidProxy = caf::Proxy< caf::SRRemid >
 
using SRRvpProxy = caf::Proxy< caf::SRRvp >
 
using SRSliceLIDProxy = caf::Proxy< caf::SRSliceLID >
 
using SRVetoProxy = caf::Proxy< caf::SRVeto >
 
using SRXnueProxy = caf::Proxy< caf::SRXnue >
 
using SRIDBranchProxy = caf::Proxy< caf::SRIDBranch >
 
using SRSliceProxy = caf::Proxy< caf::SRSlice >
 
using SRCosmicCVNProxy = caf::Proxy< caf::SRCosmicCVN >
 
using SRSpillProxy = caf::Proxy< caf::SRSpill >
 
using SRCVNFeaturesProxy = caf::Proxy< caf::SRCVNFeatures >
 
using SRPixelObjMapProxy = caf::Proxy< caf::SRPixelObjMap >
 
using SRSliceMapProxy = caf::Proxy< caf::SRSliceMap >
 
using SRTrainingDataProxy = caf::Proxy< caf::SRTrainingData >
 
using SRTrainingBranchProxy = caf::Proxy< caf::SRTrainingBranch >
 
using SRShowerProxy = caf::Proxy< caf::SRShower >
 
using SRSLidProxy = caf::Proxy< caf::SRSLid >
 
using SRShowerLIDProxy = caf::Proxy< caf::SRShowerLID >
 
using SRMRPropertiesProxy = caf::Proxy< caf::SRMRProperties >
 
using SRTrackProxy = caf::Proxy< caf::SRTrack >
 
using SRTrackBaseProxy = caf::Proxy< caf::SRTrackBase >
 
using SRKalmanTrackProxy = caf::Proxy< caf::SRKalmanTrack >
 
using SRKalmanProxy = caf::Proxy< caf::SRKalman >
 
using SRTrackBranchProxy = caf::Proxy< caf::SRTrackBranch >
 
using SRPixelMapProxy = caf::Proxy< caf::SRPixelMap >
 
using SRCVNNeutronDaughterResultProxy = caf::Proxy< caf::SRCVNNeutronDaughterResult >
 
using SRCVNParticleResultProxy = caf::Proxy< caf::SRCVNParticleResult >
 
using SRProngTrainingDataProxy = caf::Proxy< caf::SRProngTrainingData >
 
using SRProngProxy = caf::Proxy< caf::SRProng >
 
using SRBpfTrackProxy = caf::Proxy< caf::SRBpfTrack >
 
using SRBpfProxy = caf::Proxy< caf::SRBpf >
 
using SRRegCVNResultProxy = caf::Proxy< caf::SRRegCVNResult >
 
using SRFuzzyKProngProxy = caf::Proxy< caf::SRFuzzyKProng >
 
using SRFuzzyKProxy = caf::Proxy< caf::SRFuzzyK >
 
using SRElasticProxy = caf::Proxy< caf::SRElastic >
 
using SRHoughVertexProxy = caf::Proxy< caf::SRHoughVertex >
 
using SRVertexDTProxy = caf::Proxy< caf::SRVertexDT >
 
using SRVertexBranchProxy = caf::Proxy< caf::SRVertexBranch >
 
using StandardRecordProxy = caf::Proxy< caf::StandardRecord >
 
template<class T >
using VectorProxy = Proxy< std::vector< T >>
 
template<class T , unsigned int N>
using ArrayProxy = Proxy< T[N]>
 

Enumerations

enum  Det_t {
  kUNKNOWN, kNEARDET, kFARDET, kNDOS,
  kNDSBTEST, kTESTBEAM, kNDetector, kFCCDAQ
}
 Which NOvA detector? More...
 
enum  View_t { kX, kY, kXorY }
 Detector view, following GeometryObjects/PlaneGeo.h. More...
 
enum  generator_ { kUnknownGenerator = 0, kGENIE = 1, kGIBUU = 2 }
 Known generators of neutrino interactions (extend as other generators are used) More...
 
enum  gen_process_t {
  kModeUnknown = -1, kPrimary = 0, kHadElastic = 1, kDecay = 2,
  kPionInelastic = 3, kProtonInelastic = 4, kNeutronInelastic = 5, kOther = 6
}
 Interaction type responsible for particle production. More...
 
enum  mode_type_ {
  kUnknownMode = -1, kQE = 0, kRes = 1, kDIS = 2,
  kCoh = 3, kCohElastic = 4, kElectronScattering = 5, kIMDAnnihilation = 6,
  kInverseBetaDecay = 7, kGlashowResonance = 8, kAMNuGamma = 9, kMEC = 10,
  kDiffractive = 11, kEM = 12, kWeakMix = 13
}
 Neutrino interaction categories. More...
 
enum  int_type_ {
  kUnknownInteraction = -1, kNuanceOffset = 1000, kCCQE = kNuanceOffset + 1, kNCQE = kNuanceOffset + 2,
  kResCCNuProtonPiPlus = kNuanceOffset + 3, kResCCNuNeutronPi0 = kNuanceOffset + 4, kResCCNuNeutronPiPlus = kNuanceOffset + 5, kResNCNuProtonPi0 = kNuanceOffset + 6,
  kResNCNuProtonPiPlus = kNuanceOffset + 7, kResNCNuNeutronPi0 = kNuanceOffset + 8, kResNCNuNeutronPiMinus = kNuanceOffset + 9, kResCCNuBarNeutronPiMinus = kNuanceOffset + 10,
  kResCCNuBarProtonPi0 = kNuanceOffset + 11, kResCCNuBarProtonPiMinus = kNuanceOffset + 12, kResNCNuBarProtonPi0 = kNuanceOffset + 13, kResNCNuBarProtonPiPlus = kNuanceOffset + 14,
  kResNCNuBarNeutronPi0 = kNuanceOffset + 15, kResNCNuBarNeutronPiMinus = kNuanceOffset + 16, kResCCNuDeltaPlusPiPlus = kNuanceOffset + 17, kResCCNuDelta2PlusPiMinus = kNuanceOffset + 21,
  kResCCNuBarDelta0PiMinus = kNuanceOffset + 28, kResCCNuBarDeltaMinusPiPlus = kNuanceOffset + 32, kResCCNuProtonRhoPlus = kNuanceOffset + 39, kResCCNuNeutronRhoPlus = kNuanceOffset + 41,
  kResCCNuBarNeutronRhoMinus = kNuanceOffset + 46, kResCCNuBarNeutronRho0 = kNuanceOffset + 48, kResCCNuSigmaPlusKaonPlus = kNuanceOffset + 53, kResCCNuSigmaPlusKaon0 = kNuanceOffset + 55,
  kResCCNuBarSigmaMinusKaon0 = kNuanceOffset + 60, kResCCNuBarSigma0Kaon0 = kNuanceOffset + 62, kResCCNuProtonEta = kNuanceOffset + 67, kResCCNuBarNeutronEta = kNuanceOffset + 70,
  kResCCNuKaonPlusLambda0 = kNuanceOffset + 73, kResCCNuBarKaon0Lambda0 = kNuanceOffset + 76, kResCCNuProtonPiPlusPiMinus = kNuanceOffset + 79, kResCCNuProtonPi0Pi0 = kNuanceOffset + 80,
  kResCCNuBarNeutronPiPlusPiMinus = kNuanceOffset + 85, kResCCNuBarNeutronPi0Pi0 = kNuanceOffset + 86, kResCCNuBarProtonPi0Pi0 = kNuanceOffset + 90, kCCDIS = kNuanceOffset + 91,
  kNCDIS = kNuanceOffset + 92, kUnUsed1 = kNuanceOffset + 93, kUnUsed2 = kNuanceOffset + 94, kCCQEHyperon = kNuanceOffset + 95,
  kNCCOH = kNuanceOffset + 96, kCCCOH = kNuanceOffset + 97, kNuElectronElastic = kNuanceOffset + 98, kInverseMuDecay = kNuanceOffset + 99
}
 Neutrino interaction type. More...
 
enum  CAFType { kNested, kFlatMultiTree, kFlatSingleTree, kCopiedRecord }
 

Functions

void BlindThisSlice (SRSlice &slc)
 
void BlindThisTrack (SRTrack &trk)
 
void BlindThisJMShower (SRJMShower &jm)
 
void BlindThisElasticVtx (SRElastic &vtx)
 
void BlindThisCosRej (SRCosRej &cosrej)
 
void BlindThisRecord (StandardRecord *rec)
 
void FillNumuEnergyVars (const numue::NumuE &e, caf::SRNumuEnergy &sre)
 
void FillNumuLSTMEnergyVars (const LSTME::LSTMEnergy &artE, float &lstmmuon, float &lstmnu)
 
void FillNueEnergyVars (const rb::Cluster &slice, const SRVertexBranch &srVertexBranch, caf::SRNueEnergy &sre)
 
void FillHadClustVars (const rb::Cluster &hadclust, caf::SRHadClust &srhadclust)
 
void FillMRCCParentInfo (const murem::MRCCParent &parent, caf::SRMRCCParent &srparent, const art::Event &evt)
 
void FillNumuInfo (const art::Ptr< cvn::Result > &cvn, const std::vector< art::Ptr< rb::Track >> &trks, const art::Ptr< numue::NumuE > &numuE, const art::Ptr< cosrej::CosRejObj > &cosrejobj, const art::Event &evt, SRMRCCParent &parent, const art::Ptr< rb::Cluster > &slice)
 
void FillSliceInfo (const art::Ptr< rb::Cluster > &slice, SRMRCCParent &parent)
 
void FillMuIdVars (const remid::ReMId &remid, SRKalmanTrack &srTrk)
 
void FillMuonIDVars (const muonid::MuonID &muid, SRKalmanTrack &srTrk)
 
void FillJMEIDVars (const jmshower::EID &eid, caf::SRJMEid &sreid)
 
void FillCVNNeutronDaughterResultVars (const std::vector< art::Ptr< rb::PID > > &cvnneutrons, caf::SRCVNNeutronDaughterResult &cvnneutron)
 
void FillCVNParticleResultVars (const std::vector< art::Ptr< rb::PID > > &cvnparts, caf::SRCVNParticleResult &cvnpart)
 
void FillShowerVars (const rb::Shower &shw, caf::SRFuzzyKProng &srshw, const std::vector< rb::Cluster > &sliceList, std::vector< cheat::TrackIDE > &allTracks, std::vector< cheat::TrackIDE > &sliceTracks, std::vector< cheat::TrackIDE > &allTracksBirks, std::vector< cheat::TrackIDE > &sliceTracksBirks, const int &sliceIdx)
 
void FillSlidVars (const slid::ShowerLID &slid, SRShowerLID &shwlid)
 
void FillSpidVars (const slid::ShowerPID &spid, caf::SRSPid &srspid)
 
void FillLIDEnergyVars (const slid::ShowerLID &slid, caf::SRFuzzyKProng &png)
 
void FillLEMVars (const lem::PIDDetails &lem, caf::SRLem &srlem)
 
void FillRVPVars (const rvp::RVP &rvp, caf::SRRvp &srrvp)
 
void FillXnueVars (const xnue::Xnue &xnue, caf::SRXnue &srxnue)
 
void FillCosRejVars (const cosrej::CosRejObj &cosrej, caf::SRCosRej &srcosrej, caf::SRContain &srcontain)
 
void FillNueCosRejVars (const cosrej::NueCosRej &nuecosrej, caf::SRNueCosRej &srnuecosrej)
 
void FillNCCosRejVars (const ncid::NCCosRej &nccosrej, caf::SRNCCosRej &srnccosrej)
 
void FillNCPi0BkgRejVars (const ncpi0::NCPi0BkgRej &ncpi0bkgrej, caf::SRNCPi0BkgRej &srncpi0bkgrej)
 
void FillNuePreselVars (const presel::PreselObj &nuepre, caf::SRPresel &srnuepre)
 
void FillRockPreselVars (const presel::PreselObj &rockpre, caf::SRPresel &srrockpre)
 
void FillVetoVars (const presel::Veto &veto, const presel::Veto &nueveto, caf::SRVeto &srveto)
 
void FillNuonEResultVars (const cvn::Result &result, caf::SRNuonEResult &srnuone)
 
void FillCVNResultVars (const cvn::Result &result, caf::SRCVNResult &srcvn, bool NoCosmics)
 
void FillCVNPixelMaps (const cvn::PixelMap &map, caf::SRPixelObjMap &srcvnmap, bool useGeV)
 
void FillCVNPixelMaps (const cvn::PixelMap &map, caf::SRPixelMap &srcvnmap, bool useGeV)
 
void FillCVNTrainingData (const cvn::TrainingData &cvntdata, caf::SRTrainingData &srtdata)
 
void FillCVNProngTrainingData (const cvn::ProngTrainingData &cvnpdata, caf::SRProngTrainingData &srpdata)
 
void FillCVNFeaturesVars (const cvn::Features &features, caf::SRCVNFeatures &srcvnfeatures, int maxcomponents)
 
void CopyRemidVars (const SRKalmanTrack &srTrk, caf::SRRemid &remid)
 
void CopyMuonIDVars (const SRKalmanTrack &srTrk, caf::SRMuonID &muid)
 
void FillSliceLID (const SliceLID::Prediction &artSliceLID, caf::SRSliceLID &cafSliceLID)
 
void AddTrkMEToVec (const me::TrkME &michel, std::vector< SRTrkME > *vec, SRTrack &srTrk, const std::vector< rb::Cluster > &sliceList, std::vector< cheat::TrackIDE > &allTracks, std::vector< cheat::TrackIDE > &sliceTracks, std::vector< cheat::TrackIDE > &allTracksBirks, std::vector< cheat::TrackIDE > &sliceTracksBirks, int sliceIdx)
 
bool IsTrueDiF (const rb::Cluster &slice)
 Helper for FillDiFVars and FillDiFShowerVars. More...
 
void FillDiFVars (const rb::Cluster &slice, const std::vector< rb::Cluster > &sliceList, std::vector< cheat::TrackIDE > &allTracks, std::vector< cheat::TrackIDE > &sliceTracks, std::vector< cheat::TrackIDE > &allTracksBirks, std::vector< cheat::TrackIDE > &sliceTracksBirks, int sliceIdx, SRMRProperties &srMR)
 
void AddSlcMEToVec (const me::SlcME &michel, std::vector< SRSlcME > *vec, const std::vector< rb::Cluster > &sliceList, std::vector< cheat::TrackIDE > &allTracks, std::vector< cheat::TrackIDE > &sliceTracks, std::vector< cheat::TrackIDE > &allTracksBirks, std::vector< cheat::TrackIDE > &sliceTracksBirks, int sliceIdx)
 
template<class T >
void FillProngVars (const rb::Prong &prng, T &srPrng, const std::vector< rb::Cluster > &sliceList, std::vector< cheat::TrackIDE > &allTracks, std::vector< cheat::TrackIDE > &sliceTracks, std::vector< cheat::TrackIDE > &allTracksBirks, std::vector< cheat::TrackIDE > &sliceTracksBirks, int sliceIdx)
 
template void FillProngVars< SRProng > (const rb::Prong &, SRProng &, const std::vector< rb::Cluster > &, std::vector< cheat::TrackIDE > &allTracks, std::vector< cheat::TrackIDE > &sliceTracks, std::vector< cheat::TrackIDE > &allTracksBirks, std::vector< cheat::TrackIDE > &sliceTracksBirks, int)
 
template void FillProngVars< SRFuzzyKProng > (const rb::Prong &, SRFuzzyKProng &, const std::vector< rb::Cluster > &, std::vector< cheat::TrackIDE > &allTracks, std::vector< cheat::TrackIDE > &sliceTracks, std::vector< cheat::TrackIDE > &allTracksBirks, std::vector< cheat::TrackIDE > &sliceTracksBirks, int)
 
template void FillProngVars< SRTrack > (const rb::Prong &, SRTrack &, const std::vector< rb::Cluster > &, std::vector< cheat::TrackIDE > &allTracks, std::vector< cheat::TrackIDE > &sliceTracks, std::vector< cheat::TrackIDE > &allTracksBirks, std::vector< cheat::TrackIDE > &sliceTracksBirks, int)
 
void FillTrackVars (const rb::Track &trk, SRTrack &srTrk, const std::vector< rb::Cluster > &sliceList, std::vector< cheat::TrackIDE > &allTracks, std::vector< cheat::TrackIDE > &sliceTracks, std::vector< cheat::TrackIDE > &allTracksBirks, std::vector< cheat::TrackIDE > &sliceTracksBirks, int sliceIdx)
 
void FillTrackVarsWithOverlapE (const rb::Track &trk, const rb::Energy &energy, SRTrack &srTrk, const std::vector< rb::Cluster > &sliceList, std::vector< cheat::TrackIDE > &allTracks, std::vector< cheat::TrackIDE > &sliceTracks, std::vector< cheat::TrackIDE > &allTracksBirks, std::vector< cheat::TrackIDE > &sliceTracksBirks, int sliceIdx)
 Function to fill SRTrack branch if overlapping energy object is present. More...
 
void FillTrackVarsBpfFitSum (const rb::FitSum &fitsum, SRBpfTrack &srBpfTrk)
 
void FillTrackVarsBpfPid (const bpfit::BPFPId &bpfpid, SRBpfTrack &srBpfTrk)
 
void FillSliceVars (const rb::Cluster &slice, caf::SRSlice &srslice, bool allowEmpty)
 
void FillShowerVars (const rb::Shower &shw, caf::SRFuzzyKProng &srshw, const std::vector< rb::Cluster > &sliceList, std::vector< cheat::TrackIDE > &allTracks, std::vector< cheat::TrackIDE > &sliceTracks, std::vector< cheat::TrackIDE > &allTracksBirks, std::vector< cheat::TrackIDE > &sliceTracksBirks, int sliceIdx)
 
void FillDiFShowerVars (const rb::Shower &shw, const std::vector< rb::Cluster > &sliceList, std::vector< cheat::TrackIDE > &allTracks, std::vector< cheat::TrackIDE > &sliceTracks, std::vector< cheat::TrackIDE > &allTracksBirks, std::vector< cheat::TrackIDE > &sliceTracksBirks, int sliceIdx, SRMRProperties &srMR)
 
void FillTrackContainmentVars (const cosrej::TrkCntObj &, SRTrack &)
 Function to fill SRTrack containment information. More...
 
void FillNDRecoTrkVars (const ndreco::NDRecoTrkObj &, SRKalmanTrack &)
 Functions to fill SRKalmanTrack pion energy reconstruction information. More...
 
void FillNDRecoBPFVars (const ndreco::NDRecoBpfObj &pionreco, SRBpfTrack &srTrk)
 
void FillNDRecoPngVars (const ndreco::NDRecoPngObj &pionreco, SRFuzzyKProng &srPng)
 
void FillTrackInfoVars (const trackinfo::TrackInfoObj &, SRTrack &)
 Function to fill additional SRTrack information. More...
 
void FillNueSandVars (const nuesand::NueSandObj &sand, caf::SRNueSandbox &srsand)
 
void FillNumuSandVars (const numusand::NumuSandObj &sand, caf::SRNumuSandbox &srsand, caf::SRContain &srcontain)
 
void FillNusSandVars (const nussand::NusSandObj &sand, caf::SRNusSandbox &srsand)
 
SRParticleTruth FillParticleTruth (const std::vector< rb::Cluster > &sliceList, const art::PtrVector< rb::CellHit > &hits, std::vector< cheat::TrackIDE > &allTracks, std::vector< cheat::TrackIDE > &sliceTracks, std::vector< cheat::TrackIDE > &allTracksBirks, std::vector< cheat::TrackIDE > &sliceTracksBirks, int sliceIdx)
 Function to fill the particle truth for a set of hits. More...
 
void AddParticleToVec (const sim::Particle &part, std::vector< cheat::TrackIDE > &allTracks, std::vector< cheat::TrackIDE > &sliceTracks, std::vector< cheat::TrackIDE > &allTracksBirks, std::vector< cheat::TrackIDE > &sliceTracksBirks, std::vector< SRTrueParticle > *vec, const std::vector< sim::TrueEnergy > &TrueEnergies)
 
float FindDaughterVisE (const sim::Particle &part, std::vector< cheat::TrackIDE > &tracks)
 
float FindDaughterVisENonRecur (const sim::Particle &part, std::vector< cheat::TrackIDE > &tracks)
 
void AddCosmicTruthToVec (const cheat::NeutrinoEffPur *effPur, std::vector< cheat::TrackIDE > &allTracks, std::vector< cheat::TrackIDE > &sliceTracks, std::vector< cheat::TrackIDE > &allTracksBirks, std::vector< cheat::TrackIDE > &sliceTracksBirks, std::vector< SRCosmic > *vec)
 
void FindAndAddMichels (std::vector< const sim::Particle * > particles, std::vector< cheat::TrackIDE > &allTracks, std::vector< SRTrueMichelE > *michelVec)
 
void AddTrueMichelEToVec (const sim::Particle &michel, const sim::Particle &motherMuon, std::vector< cheat::TrackIDE > &allTracks, std::vector< SRTrueMichelE > *vec)
 
void AddPreFSI (const art::Ptr< simb::MCTruth > &truth, SRNeutrino &nu)
 
double TrueNeutrinoDistance (novadaq::cnv::DetId det, const SRNeutrino &nu)
 
std::vector< SRGenieWeights > GenieReweightTable (const rwgt::GENIEReweightTable &table)
 
SRFluxWeights FluxReweights (const fxwgt::FluxWeights &flxwgts)
 
SRGeant4Weights Geant4Reweights (const g4rwgt::G4WeightTable &g4wgts)
 
caf::generator_ CAFGeneratorEnum (simb::Generator_t simbGeneratorEnum)
 
std::vector< unsigned int > DecodeGeneratorVersion (const std::string &versionString)
 
template<class T >
std::vector< T > PtrVecToVec (const art::PtrVector< T > &xs)
 
CAFType GetCAFType (const TDirectory *dir, TTree *tr)
 
std::string StripSubscripts (const std::string &s)
 
int NSubscripts (const std::string &name)
 Count the subscripts in the name. More...
 
const long & AdjustBase (const long &base, CAFType type, const std::string &name)
 
template<class T >
void GetTypedValueWrapper (TLeaf *leaf, T &x, int subidx)
 
void GetTypedValueWrapper (TLeaf *leaf, std::string &x, int subidx)
 
template<class T >
void EvalInstanceWrapper (TTreeFormula *ttf, T &x)
 
void EvalInstanceWrapper (TTreeFormula *ttf, std::string &x)
 
template<class T >
bool operator< (const Proxy< std::vector< T >> &a, const std::vector< T > &b)
 Used in comparison of GENIE version numbers. More...
 

Variables

static const int blind_val = -5
 
static const unsigned int blind_val_unsigned = 0
 
const std::string whattracker = "bpf"
 
const std::string slicelabel = "slicer"
 
const std::string trklabelkalman = "kalmantrackmerge"
 
const std::string trklabelbpf = "breakpoint"
 
const std::string numuelabel = "numue"
 
const std::string cosrejlabel = "cosrej"
 
const std::string remidlabel = "remid"
 
const std::string cvnlabel = "cvnevaloldpresel"
 
const std::string vtxlabel = "elasticarmshs"
 
const std::string pronglabel = "fuzzykvertex:Prongs3D"
 
const std::string cvnpronglabel = "cvnprongeval"
 
const int kNumCVNFinalStates = 387
 
const CVNFinalState cvnStates [kNumCVNFinalStates]
 
const int kCVN_nM_Other = kNumCVNFinalStates
 
const int kCVN_nE_Other = kNumCVNFinalStates+1
 
const int kCVN_nT_Other = kNumCVNFinalStates+2
 
const int kCVN_Cosmic_PT = kNumCVNFinalStates+3
 
const int kCVN_Other_PT = kNumCVNFinalStates+4
 
const long kZero = 0
 
const long kDummyBase = -1
 

Detailed Description

This module creates Common Analysis Files.

Common Analysis Files.

Class to make all of the CAFCuts for the RecoValidation modules.

Author
rocco.nosp@m.@phy.nosp@m.sics..nosp@m.umn..nosp@m.edu
Date
July 2012
Authors
mbair.nosp@m.d42@.nosp@m.india.nosp@m.na.e.nosp@m.du
Date
2017/02/14

Typedef Documentation

template<class T , unsigned int N>
using caf::ArrayProxy = typedef Proxy<T[N]>

Definition at line 301 of file BasicTypesProxy.h.

using caf::FlatSRBeam = typedef flat::Flat<caf::SRBeam>

Definition at line 86 of file FwdDeclare.h.

using caf::FlatSRBpf = typedef flat::Flat<caf::SRBpf>

Definition at line 428 of file FwdDeclare.h.

using caf::FlatSRBPFEnergy = typedef flat::Flat<caf::SRBPFEnergy>

Definition at line 26 of file FwdDeclare.h.

using caf::FlatSRBpfId = typedef flat::Flat<caf::SRBpfId>

Definition at line 182 of file FwdDeclare.h.

using caf::FlatSRBpfTrack = typedef flat::Flat<caf::SRBpfTrack>

Definition at line 422 of file FwdDeclare.h.

using caf::FlatSRContain = typedef flat::Flat<caf::SRContain>

Definition at line 188 of file FwdDeclare.h.

using caf::FlatSRCosmic = typedef flat::Flat<caf::SRCosmic>

Definition at line 80 of file FwdDeclare.h.

using caf::FlatSRCosmicCVN = typedef flat::Flat<caf::SRCosmicCVN>

Definition at line 296 of file FwdDeclare.h.

using caf::FlatSRCosRej = typedef flat::Flat<caf::SRCosRej>

Definition at line 194 of file FwdDeclare.h.

using caf::FlatSRCVNFeatures = typedef flat::Flat<caf::SRCVNFeatures>

Definition at line 308 of file FwdDeclare.h.

using caf::FlatSRCVNNeutronDaughterResult = typedef flat::Flat<caf::SRCVNNeutronDaughterResult>

Definition at line 398 of file FwdDeclare.h.

using caf::FlatSRCVNParticleResult = typedef flat::Flat<caf::SRCVNParticleResult>

Definition at line 404 of file FwdDeclare.h.

using caf::FlatSRCVNResult = typedef flat::Flat<caf::SRCVNResult>

Definition at line 200 of file FwdDeclare.h.

using caf::FlatSRElastic = typedef flat::Flat<caf::SRElastic>

Definition at line 452 of file FwdDeclare.h.

using caf::FlatSRELid = typedef flat::Flat<caf::SRELid>

Definition at line 212 of file FwdDeclare.h.

using caf::FlatSREnergyBranch = typedef flat::Flat<caf::SREnergyBranch>

Definition at line 56 of file FwdDeclare.h.

using caf::FlatSRFluxWeights = typedef flat::Flat<caf::SRFluxWeights>

Definition at line 116 of file FwdDeclare.h.

using caf::FlatSRFuzzyK = typedef flat::Flat<caf::SRFuzzyK>

Definition at line 446 of file FwdDeclare.h.

using caf::FlatSRFuzzyKProng = typedef flat::Flat<caf::SRFuzzyKProng>

Definition at line 440 of file FwdDeclare.h.

using caf::FlatSRGeant4Weights = typedef flat::Flat<caf::SRGeant4Weights>

Definition at line 104 of file FwdDeclare.h.

using caf::FlatSRGenieWeights = typedef flat::Flat<caf::SRGenieWeights>

Definition at line 110 of file FwdDeclare.h.

using caf::FlatSRGlobalTruth = typedef flat::Flat<caf::SRGlobalTruth>

Definition at line 134 of file FwdDeclare.h.

using caf::FlatSRHadClust = typedef flat::Flat<caf::SRHadClust>

Definition at line 38 of file FwdDeclare.h.

using caf::FlatSRHeader = typedef flat::Flat<caf::SRHeader>

Definition at line 62 of file FwdDeclare.h.

using caf::FlatSRHoughVertex = typedef flat::Flat<caf::SRHoughVertex>

Definition at line 458 of file FwdDeclare.h.

using caf::FlatSRIDBranch = typedef flat::Flat<caf::SRIDBranch>

Definition at line 284 of file FwdDeclare.h.

using caf::FlatSRKalman = typedef flat::Flat<caf::SRKalman>

Definition at line 380 of file FwdDeclare.h.

using caf::FlatSRKalmanTrack = typedef flat::Flat<caf::SRKalmanTrack>

Definition at line 374 of file FwdDeclare.h.

using caf::FlatSRLem = typedef flat::Flat<caf::SRLem>

Definition at line 206 of file FwdDeclare.h.

using caf::FlatSRLorentzVector = typedef flat::Flat<caf::SRLorentzVector>

Definition at line 68 of file FwdDeclare.h.

using caf::FlatSRMCReweight = typedef flat::Flat<caf::SRMCReweight>

Definition at line 122 of file FwdDeclare.h.

using caf::FlatSRMichelE = typedef flat::Flat<caf::SRMichelE>

Definition at line 164 of file FwdDeclare.h.

using caf::FlatSRMRCCParent = typedef flat::Flat<caf::SRMRCCParent>

Definition at line 170 of file FwdDeclare.h.

using caf::FlatSRMRProperties = typedef flat::Flat<caf::SRMRProperties>

Definition at line 356 of file FwdDeclare.h.

using caf::FlatSRMuonID = typedef flat::Flat<caf::SRMuonID>

Definition at line 218 of file FwdDeclare.h.

using caf::FlatSRNCCosRej = typedef flat::Flat<caf::SRNCCosRej>

Definition at line 224 of file FwdDeclare.h.

using caf::FlatSRNCPi0BkgRej = typedef flat::Flat<caf::SRNCPi0BkgRej>

Definition at line 230 of file FwdDeclare.h.

using caf::FlatSRNeutrino = typedef flat::Flat<caf::SRNeutrino>

Definition at line 128 of file FwdDeclare.h.

using caf::FlatSRNueCosRej = typedef flat::Flat<caf::SRNueCosRej>

Definition at line 236 of file FwdDeclare.h.

using caf::FlatSRNueEnergy = typedef flat::Flat<caf::SRNueEnergy>

Definition at line 20 of file FwdDeclare.h.

using caf::FlatSRNuGenParticle = typedef flat::Flat<caf::SRNuGenParticle>

Definition at line 92 of file FwdDeclare.h.

using caf::FlatSRNumuEnergy = typedef flat::Flat<caf::SRNumuEnergy>

Definition at line 50 of file FwdDeclare.h.

using caf::FlatSRNuonEResult = typedef flat::Flat<caf::SRNuonEResult>

Definition at line 248 of file FwdDeclare.h.

using caf::FlatSRParentBranch = typedef flat::Flat<caf::SRParentBranch>

Definition at line 176 of file FwdDeclare.h.

using caf::FlatSRParticleTruth = typedef flat::Flat<caf::SRParticleTruth>

Definition at line 146 of file FwdDeclare.h.

using caf::FlatSRPixelMap = typedef flat::Flat<caf::SRPixelMap>

Definition at line 392 of file FwdDeclare.h.

using caf::FlatSRPixelObjMap = typedef flat::Flat<caf::SRPixelObjMap>

Definition at line 314 of file FwdDeclare.h.

using caf::FlatSRPresel = typedef flat::Flat<caf::SRPresel>

Definition at line 242 of file FwdDeclare.h.

using caf::FlatSRProng = typedef flat::Flat<caf::SRProng>

Definition at line 416 of file FwdDeclare.h.

using caf::FlatSRProngTrainingData = typedef flat::Flat<caf::SRProngTrainingData>

Definition at line 410 of file FwdDeclare.h.

using caf::FlatSRRegCVNResult = typedef flat::Flat<caf::SRRegCVNResult>

Definition at line 434 of file FwdDeclare.h.

using caf::FlatSRRemid = typedef flat::Flat<caf::SRRemid>

Definition at line 254 of file FwdDeclare.h.

using caf::FlatSRRvp = typedef flat::Flat<caf::SRRvp>

Definition at line 260 of file FwdDeclare.h.

using caf::FlatSRShower = typedef flat::Flat<caf::SRShower>

Definition at line 338 of file FwdDeclare.h.

using caf::FlatSRShowerLID = typedef flat::Flat<caf::SRShowerLID>

Definition at line 350 of file FwdDeclare.h.

using caf::FlatSRSlcME = typedef flat::Flat<caf::SRSlcME>

Definition at line 152 of file FwdDeclare.h.

using caf::FlatSRSlice = typedef flat::Flat<caf::SRSlice>

Definition at line 290 of file FwdDeclare.h.

using caf::FlatSRSliceLID = typedef flat::Flat<caf::SRSliceLID>

Definition at line 266 of file FwdDeclare.h.

using caf::FlatSRSliceMap = typedef flat::Flat<caf::SRSliceMap>

Definition at line 320 of file FwdDeclare.h.

using caf::FlatSRSLid = typedef flat::Flat<caf::SRSLid>

Definition at line 344 of file FwdDeclare.h.

using caf::FlatSRSLidEnergy = typedef flat::Flat<caf::SRSLidEnergy>

Definition at line 14 of file FwdDeclare.h.

using caf::FlatSRSpill = typedef flat::Flat<caf::SRSpill>

Definition at line 302 of file FwdDeclare.h.

using caf::FlatSRTrack = typedef flat::Flat<caf::SRTrack>

Definition at line 362 of file FwdDeclare.h.

using caf::FlatSRTrackBase = typedef flat::Flat<caf::SRTrackBase>

Definition at line 368 of file FwdDeclare.h.

using caf::FlatSRTrackBranch = typedef flat::Flat<caf::SRTrackBranch>

Definition at line 386 of file FwdDeclare.h.

using caf::FlatSRTrainingBranch = typedef flat::Flat<caf::SRTrainingBranch>

Definition at line 332 of file FwdDeclare.h.

using caf::FlatSRTrainingData = typedef flat::Flat<caf::SRTrainingData>

Definition at line 326 of file FwdDeclare.h.

using caf::FlatSRTrkME = typedef flat::Flat<caf::SRTrkME>

Definition at line 158 of file FwdDeclare.h.

using caf::FlatSRTrueMichelE = typedef flat::Flat<caf::SRTrueMichelE>

Definition at line 74 of file FwdDeclare.h.

using caf::FlatSRTrueNumuEnergy = typedef flat::Flat<caf::SRTrueNumuEnergy>

Definition at line 44 of file FwdDeclare.h.

using caf::FlatSRTrueParticle = typedef flat::Flat<caf::SRTrueParticle>

Definition at line 98 of file FwdDeclare.h.

using caf::FlatSRTruthBranch = typedef flat::Flat<caf::SRTruthBranch>

Definition at line 140 of file FwdDeclare.h.

using caf::FlatSRVector3D = typedef flat::Flat<caf::SRVector3D>

Definition at line 32 of file FwdDeclare.h.

using caf::FlatSRVertexBranch = typedef flat::Flat<caf::SRVertexBranch>

Definition at line 470 of file FwdDeclare.h.

using caf::FlatSRVertexDT = typedef flat::Flat<caf::SRVertexDT>

Definition at line 464 of file FwdDeclare.h.

using caf::FlatSRVeto = typedef flat::Flat<caf::SRVeto>

Definition at line 272 of file FwdDeclare.h.

using caf::FlatSRXnue = typedef flat::Flat<caf::SRXnue>

Definition at line 278 of file FwdDeclare.h.

using caf::FlatStandardRecord = typedef flat::Flat<caf::StandardRecord>

Definition at line 476 of file FwdDeclare.h.

using caf::SRBeamProxy = typedef caf::Proxy<caf::SRBeam>

Definition at line 86 of file FwdDeclare.h.

using caf::SRBPFEnergyProxy = typedef caf::Proxy<caf::SRBPFEnergy>

Definition at line 26 of file FwdDeclare.h.

using caf::SRBpfIdProxy = typedef caf::Proxy<caf::SRBpfId>

Definition at line 182 of file FwdDeclare.h.

using caf::SRBpfProxy = typedef caf::Proxy<caf::SRBpf>

Definition at line 428 of file FwdDeclare.h.

using caf::SRBpfTrackProxy = typedef caf::Proxy<caf::SRBpfTrack>

Definition at line 422 of file FwdDeclare.h.

using caf::SRContainProxy = typedef caf::Proxy<caf::SRContain>

Definition at line 188 of file FwdDeclare.h.

using caf::SRCosmicCVNProxy = typedef caf::Proxy<caf::SRCosmicCVN>

Definition at line 296 of file FwdDeclare.h.

using caf::SRCosmicProxy = typedef caf::Proxy<caf::SRCosmic>

Definition at line 80 of file FwdDeclare.h.

using caf::SRCosRejProxy = typedef caf::Proxy<caf::SRCosRej>

Definition at line 194 of file FwdDeclare.h.

using caf::SRCVNFeaturesProxy = typedef caf::Proxy<caf::SRCVNFeatures>

Definition at line 308 of file FwdDeclare.h.

using caf::SRCVNNeutronDaughterResultProxy = typedef caf::Proxy<caf::SRCVNNeutronDaughterResult>

Definition at line 398 of file FwdDeclare.h.

using caf::SRCVNParticleResultProxy = typedef caf::Proxy<caf::SRCVNParticleResult>

Definition at line 404 of file FwdDeclare.h.

using caf::SRCVNResultProxy = typedef caf::Proxy<caf::SRCVNResult>

Definition at line 200 of file FwdDeclare.h.

using caf::SRElasticProxy = typedef caf::Proxy<caf::SRElastic>

Definition at line 452 of file FwdDeclare.h.

using caf::SRELidProxy = typedef caf::Proxy<caf::SRELid>

Definition at line 212 of file FwdDeclare.h.

using caf::SREnergyBranchProxy = typedef caf::Proxy<caf::SREnergyBranch>

Definition at line 56 of file FwdDeclare.h.

using caf::SRFluxWeightsProxy = typedef caf::Proxy<caf::SRFluxWeights>

Definition at line 116 of file FwdDeclare.h.

using caf::SRFuzzyKProngProxy = typedef caf::Proxy<caf::SRFuzzyKProng>

Definition at line 440 of file FwdDeclare.h.

using caf::SRFuzzyKProxy = typedef caf::Proxy<caf::SRFuzzyK>

Definition at line 446 of file FwdDeclare.h.

using caf::SRGeant4WeightsProxy = typedef caf::Proxy<caf::SRGeant4Weights>

Definition at line 104 of file FwdDeclare.h.

using caf::SRGenieWeightsProxy = typedef caf::Proxy<caf::SRGenieWeights>

Definition at line 110 of file FwdDeclare.h.

using caf::SRGlobalTruthProxy = typedef caf::Proxy<caf::SRGlobalTruth>

Definition at line 134 of file FwdDeclare.h.

using caf::SRHadClustProxy = typedef caf::Proxy<caf::SRHadClust>

Definition at line 38 of file FwdDeclare.h.

using caf::SRHeaderProxy = typedef caf::Proxy<caf::SRHeader>

Definition at line 62 of file FwdDeclare.h.

using caf::SRHoughVertexProxy = typedef caf::Proxy<caf::SRHoughVertex>

Definition at line 458 of file FwdDeclare.h.

using caf::SRIDBranchProxy = typedef caf::Proxy<caf::SRIDBranch>

Definition at line 284 of file FwdDeclare.h.

using caf::SRKalmanProxy = typedef caf::Proxy<caf::SRKalman>

Definition at line 380 of file FwdDeclare.h.

using caf::SRKalmanTrackProxy = typedef caf::Proxy<caf::SRKalmanTrack>

Definition at line 374 of file FwdDeclare.h.

using caf::SRLemProxy = typedef caf::Proxy<caf::SRLem>

Definition at line 206 of file FwdDeclare.h.

using caf::SRLorentzVectorProxy = typedef caf::Proxy<caf::SRLorentzVector>

Definition at line 68 of file FwdDeclare.h.

using caf::SRMCReweightProxy = typedef caf::Proxy<caf::SRMCReweight>

Definition at line 122 of file FwdDeclare.h.

using caf::SRMichelEProxy = typedef caf::Proxy<caf::SRMichelE>

Definition at line 164 of file FwdDeclare.h.

using caf::SRMRCCParentProxy = typedef caf::Proxy<caf::SRMRCCParent>

Definition at line 170 of file FwdDeclare.h.

using caf::SRMRPropertiesProxy = typedef caf::Proxy<caf::SRMRProperties>

Definition at line 356 of file FwdDeclare.h.

using caf::SRMuonIDProxy = typedef caf::Proxy<caf::SRMuonID>

Definition at line 218 of file FwdDeclare.h.

using caf::SRNCCosRejProxy = typedef caf::Proxy<caf::SRNCCosRej>

Definition at line 224 of file FwdDeclare.h.

using caf::SRNCPi0BkgRejProxy = typedef caf::Proxy<caf::SRNCPi0BkgRej>

Definition at line 230 of file FwdDeclare.h.

typedef Proxy< SRNeutrino > caf::SRNeutrinoProxy

Definition at line 128 of file FwdDeclare.h.

using caf::SRNueCosRejProxy = typedef caf::Proxy<caf::SRNueCosRej>

Definition at line 236 of file FwdDeclare.h.

using caf::SRNueEnergyProxy = typedef caf::Proxy<caf::SRNueEnergy>

Definition at line 20 of file FwdDeclare.h.

using caf::SRNuGenParticleProxy = typedef caf::Proxy<caf::SRNuGenParticle>

Definition at line 92 of file FwdDeclare.h.

using caf::SRNumuEnergyProxy = typedef caf::Proxy<caf::SRNumuEnergy>

Definition at line 50 of file FwdDeclare.h.

using caf::SRNuonEResultProxy = typedef caf::Proxy<caf::SRNuonEResult>

Definition at line 248 of file FwdDeclare.h.

using caf::SRParentBranchProxy = typedef caf::Proxy<caf::SRParentBranch>

Definition at line 176 of file FwdDeclare.h.

using caf::SRParticleTruthProxy = typedef caf::Proxy<caf::SRParticleTruth>

Definition at line 146 of file FwdDeclare.h.

using caf::SRPixelMapProxy = typedef caf::Proxy<caf::SRPixelMap>

Definition at line 392 of file FwdDeclare.h.

using caf::SRPixelObjMapProxy = typedef caf::Proxy<caf::SRPixelObjMap>

Definition at line 314 of file FwdDeclare.h.

using caf::SRPreselProxy = typedef caf::Proxy<caf::SRPresel>

Definition at line 242 of file FwdDeclare.h.

using caf::SRProngProxy = typedef caf::Proxy<caf::SRProng>

Definition at line 416 of file FwdDeclare.h.

using caf::SRProngTrainingDataProxy = typedef caf::Proxy<caf::SRProngTrainingData>

Definition at line 410 of file FwdDeclare.h.

typedef Proxy< StandardRecord > caf::SRProxy

Definition at line 2 of file EpilogFwd.h.

using caf::SRRegCVNResultProxy = typedef caf::Proxy<caf::SRRegCVNResult>

Definition at line 434 of file FwdDeclare.h.

using caf::SRRemidProxy = typedef caf::Proxy<caf::SRRemid>

Definition at line 254 of file FwdDeclare.h.

using caf::SRRvpProxy = typedef caf::Proxy<caf::SRRvp>

Definition at line 260 of file FwdDeclare.h.

using caf::SRShowerLIDProxy = typedef caf::Proxy<caf::SRShowerLID>

Definition at line 350 of file FwdDeclare.h.

using caf::SRShowerProxy = typedef caf::Proxy<caf::SRShower>

Definition at line 338 of file FwdDeclare.h.

using caf::SRSlcMEProxy = typedef caf::Proxy<caf::SRSlcME>

Definition at line 152 of file FwdDeclare.h.

using caf::SRSliceLIDProxy = typedef caf::Proxy<caf::SRSliceLID>

Definition at line 266 of file FwdDeclare.h.

using caf::SRSliceMapProxy = typedef caf::Proxy<caf::SRSliceMap>

Definition at line 320 of file FwdDeclare.h.

using caf::SRSliceProxy = typedef caf::Proxy<caf::SRSlice>

Definition at line 290 of file FwdDeclare.h.

using caf::SRSLidEnergyProxy = typedef caf::Proxy<caf::SRSLidEnergy>

Definition at line 14 of file FwdDeclare.h.

using caf::SRSLidProxy = typedef caf::Proxy<caf::SRSLid>

Definition at line 344 of file FwdDeclare.h.

typedef Proxy< SRSpill > caf::SRSpillProxy

Definition at line 302 of file FwdDeclare.h.

using caf::SRTrackBaseProxy = typedef caf::Proxy<caf::SRTrackBase>

Definition at line 368 of file FwdDeclare.h.

using caf::SRTrackBranchProxy = typedef caf::Proxy<caf::SRTrackBranch>

Definition at line 386 of file FwdDeclare.h.

using caf::SRTrackProxy = typedef caf::Proxy<caf::SRTrack>

Definition at line 362 of file FwdDeclare.h.

using caf::SRTrainingBranchProxy = typedef caf::Proxy<caf::SRTrainingBranch>

Definition at line 332 of file FwdDeclare.h.

using caf::SRTrainingDataProxy = typedef caf::Proxy<caf::SRTrainingData>

Definition at line 326 of file FwdDeclare.h.

using caf::SRTrkMEProxy = typedef caf::Proxy<caf::SRTrkME>

Definition at line 158 of file FwdDeclare.h.

using caf::SRTrueMichelEProxy = typedef caf::Proxy<caf::SRTrueMichelE>

Definition at line 74 of file FwdDeclare.h.

using caf::SRTrueNumuEnergyProxy = typedef caf::Proxy<caf::SRTrueNumuEnergy>

Definition at line 44 of file FwdDeclare.h.

using caf::SRTrueParticleProxy = typedef caf::Proxy<caf::SRTrueParticle>

Definition at line 98 of file FwdDeclare.h.

using caf::SRTruthBranchProxy = typedef caf::Proxy<caf::SRTruthBranch>

Definition at line 140 of file FwdDeclare.h.

using caf::SRVector3DProxy = typedef caf::Proxy<caf::SRVector3D>

Definition at line 32 of file FwdDeclare.h.

using caf::SRVertexBranchProxy = typedef caf::Proxy<caf::SRVertexBranch>

Definition at line 470 of file FwdDeclare.h.

using caf::SRVertexDTProxy = typedef caf::Proxy<caf::SRVertexDT>

Definition at line 464 of file FwdDeclare.h.

using caf::SRVetoProxy = typedef caf::Proxy<caf::SRVeto>

Definition at line 272 of file FwdDeclare.h.

using caf::SRXnueProxy = typedef caf::Proxy<caf::SRXnue>

Definition at line 278 of file FwdDeclare.h.

using caf::StandardRecordProxy = typedef caf::Proxy<caf::StandardRecord>

Definition at line 476 of file FwdDeclare.h.

template<class T >
using caf::VectorProxy = typedef Proxy<std::vector<T>>

Definition at line 235 of file BasicTypesProxy.h.

Enumeration Type Documentation

enum caf::CAFType
Enumerator
kNested 
kFlatMultiTree 
kFlatSingleTree 
kCopiedRecord 

Definition at line 31 of file BasicTypesProxy.h.

32  {
33  kNested,
34  kFlatMultiTree,
35  kFlatSingleTree,
36  kCopiedRecord // Assigned into, not associated with a file
37  };
enum caf::Det_t

Which NOvA detector?

Enumerator
kUNKNOWN 

Unknown detector.

kNEARDET 

Near Detector underground.

kFARDET 

Far Detector at Ash River.

kNDOS 

Prototype Near Detector on the Surface.

kNDSBTEST 
kTESTBEAM 
kNDetector 
kFCCDAQ 

Definition at line 7 of file SREnums.h.

8  {
9  kUNKNOWN, ///< Unknown detector
10  kNEARDET, ///< Near Detector underground
11  kFARDET, ///< Far Detector at Ash River
12  kNDOS, ///< Prototype Near Detector on the Surface
13  kNDSBTEST,
14  kTESTBEAM,
15  kNDetector,
16  kFCCDAQ
17  };
caf::kNEARDET
Near Detector underground.
Definition: SREnums.h:10
caf::kUNKNOWN
Unknown detector.
Definition: SREnums.h:9
caf::kFARDET
Far Detector at Ash River.
Definition: SREnums.h:11
caf::kNDOS
Prototype Near Detector on the Surface.
Definition: SREnums.h:12
enum caf::gen_process_t

Interaction type responsible for particle production.

Enumerator
kModeUnknown 
kPrimary 
kHadElastic 
kDecay 
kPionInelastic 
kProtonInelastic 
kNeutronInelastic 
kOther 

Definition at line 35 of file SREnums.h.

35  {
36  kModeUnknown = -1,
37  kPrimary = 0,
38  kHadElastic = 1,
39  kDecay = 2,
40  kPionInelastic = 3,
41  kProtonInelastic = 4,
42  kNeutronInelastic = 5,
43  kOther = 6
44  };
enum caf::generator_

Known generators of neutrino interactions (extend as other generators are used)

Enumerator
kUnknownGenerator 
kGENIE 
kGIBUU 

Definition at line 28 of file SREnums.h.

28  {
29  kUnknownGenerator = 0,
30  kGENIE = 1,
31  kGIBUU = 2
32  };
enum caf::int_type_

Neutrino interaction type.

Enumerator
kUnknownInteraction 
kNuanceOffset 

offset to account for adding in Nuance codes to this enum

kCCQE 

charged current quasi-elastic

kNCQE 

neutral current quasi-elastic

kResCCNuProtonPiPlus 

resonant charged current, $\nu p \to l^- p \pi^+$

kResCCNuNeutronPi0 

resonant charged current, $\nu n \to l^- p \pi^0$

kResCCNuNeutronPiPlus 

resonant charged current, $\nu n \to l^- n \pi^+$

kResNCNuProtonPi0 

resonant neutral current, $\nu p \to \nu p \pi^0$

kResNCNuProtonPiPlus 

resonant neutral current, $\nu p \to \nu p \pi^+$

kResNCNuNeutronPi0 

resonant neutral current, $\nu n \to \nu n \pi^0$

kResNCNuNeutronPiMinus 

resonant neutral current, $\nu n \to \nu p \pi^-$

kResCCNuBarNeutronPiMinus 

resonant charged current, $\bar\nu n \to l^+ n \pi^-$

kResCCNuBarProtonPi0 

resonant charged current, $\bar\nu p \to l^+ n \pi^0$

kResCCNuBarProtonPiMinus 

resonant charged current, $\bar\nu p \to l^+ p \pi^-$

kResNCNuBarProtonPi0 

resonant charged current, $\bar\nu p \to \bar\nu p \pi^0$

kResNCNuBarProtonPiPlus 

resonant charged current, $\bar\nu p \to \bar\nu n \pi^+$

kResNCNuBarNeutronPi0 

resonant charged current, $\bar\nu n \to \bar\nu n \pi^0$

kResNCNuBarNeutronPiMinus 

resonant charged current, $\bar\nu n \to \bar\nu p \pi^-$

kResCCNuDeltaPlusPiPlus 
kResCCNuDelta2PlusPiMinus 
kResCCNuBarDelta0PiMinus 
kResCCNuBarDeltaMinusPiPlus 
kResCCNuProtonRhoPlus 
kResCCNuNeutronRhoPlus 
kResCCNuBarNeutronRhoMinus 
kResCCNuBarNeutronRho0 
kResCCNuSigmaPlusKaonPlus 
kResCCNuSigmaPlusKaon0 
kResCCNuBarSigmaMinusKaon0 
kResCCNuBarSigma0Kaon0 
kResCCNuProtonEta 
kResCCNuBarNeutronEta 
kResCCNuKaonPlusLambda0 
kResCCNuBarKaon0Lambda0 
kResCCNuProtonPiPlusPiMinus 
kResCCNuProtonPi0Pi0 
kResCCNuBarNeutronPiPlusPiMinus 
kResCCNuBarNeutronPi0Pi0 
kResCCNuBarProtonPi0Pi0 
kCCDIS 

charged current deep inelastic scatter

kNCDIS 

charged current deep inelastic scatter

kUnUsed1 
kUnUsed2 
kCCQEHyperon 
kNCCOH 
kCCCOH 

charged current coherent pion

kNuElectronElastic 

neutrino electron elastic scatter

kInverseMuDecay 

inverse muon decay

Definition at line 66 of file SREnums.h.

66  {
67  kUnknownInteraction = -1,
68  kNuanceOffset = 1000, ///< offset to account for adding in Nuance codes to this enum
69  kCCQE = kNuanceOffset + 1, ///< charged current quasi-elastic
70  kNCQE = kNuanceOffset + 2, ///< neutral current quasi-elastic
71  kResCCNuProtonPiPlus = kNuanceOffset + 3, ///< resonant charged current, \f$\nu p \to l^- p \pi^+\f$
72  kResCCNuNeutronPi0 = kNuanceOffset + 4, ///< resonant charged current, \f$\nu n \to l^- p \pi^0\f$
73  kResCCNuNeutronPiPlus = kNuanceOffset + 5, ///< resonant charged current, \f$\nu n \to l^- n \pi^+\f$
74  kResNCNuProtonPi0 = kNuanceOffset + 6, ///< resonant neutral current, \f$\nu p \to \nu p \pi^0\f$
75  kResNCNuProtonPiPlus = kNuanceOffset + 7, ///< resonant neutral current, \f$\nu p \to \nu p \pi^+\f$
76  kResNCNuNeutronPi0 = kNuanceOffset + 8, ///< resonant neutral current, \f$\nu n \to \nu n \pi^0\f$
77  kResNCNuNeutronPiMinus = kNuanceOffset + 9, ///< resonant neutral current, \f$\nu n \to \nu p \pi^-\f$
78  kResCCNuBarNeutronPiMinus = kNuanceOffset + 10, ///< resonant charged current, \f$\bar\nu n \to l^+ n \pi^-\f$
79  kResCCNuBarProtonPi0 = kNuanceOffset + 11, ///< resonant charged current, \f$\bar\nu p \to l^+ n \pi^0\f$
80  kResCCNuBarProtonPiMinus = kNuanceOffset + 12, ///< resonant charged current, \f$\bar\nu p \to l^+ p \pi^-\f$
81  kResNCNuBarProtonPi0 = kNuanceOffset + 13, ///< resonant charged current, \f$\bar\nu p \to \bar\nu p \pi^0\f$
82  kResNCNuBarProtonPiPlus = kNuanceOffset + 14, ///< resonant charged current, \f$\bar\nu p \to \bar\nu n \pi^+\f$
83  kResNCNuBarNeutronPi0 = kNuanceOffset + 15, ///< resonant charged current, \f$\bar\nu n \to \bar\nu n \pi^0\f$
84  kResNCNuBarNeutronPiMinus = kNuanceOffset + 16, ///< resonant charged current, \f$\bar\nu n \to \bar\nu p \pi^-\f$
85  kResCCNuDeltaPlusPiPlus = kNuanceOffset + 17,
86  kResCCNuDelta2PlusPiMinus = kNuanceOffset + 21,
87  kResCCNuBarDelta0PiMinus = kNuanceOffset + 28,
88  kResCCNuBarDeltaMinusPiPlus= kNuanceOffset + 32,
89  kResCCNuProtonRhoPlus = kNuanceOffset + 39,
90  kResCCNuNeutronRhoPlus = kNuanceOffset + 41,
91  kResCCNuBarNeutronRhoMinus = kNuanceOffset + 46,
92  kResCCNuBarNeutronRho0 = kNuanceOffset + 48,
93  kResCCNuSigmaPlusKaonPlus = kNuanceOffset + 53,
94  kResCCNuSigmaPlusKaon0 = kNuanceOffset + 55,
95  kResCCNuBarSigmaMinusKaon0 = kNuanceOffset + 60,
96  kResCCNuBarSigma0Kaon0 = kNuanceOffset + 62,
97  kResCCNuProtonEta = kNuanceOffset + 67,
98  kResCCNuBarNeutronEta = kNuanceOffset + 70,
99  kResCCNuKaonPlusLambda0 = kNuanceOffset + 73,
100  kResCCNuBarKaon0Lambda0 = kNuanceOffset + 76,
101  kResCCNuProtonPiPlusPiMinus= kNuanceOffset + 79,
102  kResCCNuProtonPi0Pi0 = kNuanceOffset + 80,
103  kResCCNuBarNeutronPiPlusPiMinus = kNuanceOffset + 85,
104  kResCCNuBarNeutronPi0Pi0 = kNuanceOffset + 86,
105  kResCCNuBarProtonPi0Pi0 = kNuanceOffset + 90,
106  kCCDIS = kNuanceOffset + 91, ///< charged current deep inelastic scatter
107  kNCDIS = kNuanceOffset + 92, ///< charged current deep inelastic scatter
108  kUnUsed1 = kNuanceOffset + 93,
109  kUnUsed2 = kNuanceOffset + 94,
110  kCCQEHyperon = kNuanceOffset + 95,
111  kNCCOH = kNuanceOffset + 96,
112  kCCCOH = kNuanceOffset + 97, ///< charged current coherent pion
113  kNuElectronElastic = kNuanceOffset + 98, ///< neutrino electron elastic scatter
114  kInverseMuDecay = kNuanceOffset + 99 ///< inverse muon decay
115  };
caf::kResNCNuProtonPiPlus
resonant neutral current,
Definition: SREnums.h:75
caf::kResCCNuNeutronPi0
resonant charged current,
Definition: SREnums.h:72
caf::kResNCNuBarProtonPiPlus
resonant charged current,
Definition: SREnums.h:82
caf::kResCCNuBarProtonPiMinus
resonant charged current,
Definition: SREnums.h:80
caf::kCCCOH
charged current coherent pion
Definition: SREnums.h:112
caf::kResNCNuBarProtonPi0
resonant charged current,
Definition: SREnums.h:81
caf::kCCDIS
charged current deep inelastic scatter
Definition: SREnums.h:106
caf::kNCDIS
charged current deep inelastic scatter
Definition: SREnums.h:107
caf::kNCQE
neutral current quasi-elastic
Definition: SREnums.h:70
caf::kResNCNuProtonPi0
resonant neutral current,
Definition: SREnums.h:74
caf::kNuanceOffset
offset to account for adding in Nuance codes to this enum
Definition: SREnums.h:68
caf::kResCCNuBarProtonPi0
resonant charged current,
Definition: SREnums.h:79
caf::kResCCNuProtonPiPlus
resonant charged current,
Definition: SREnums.h:71
caf::kResCCNuBarNeutronPiMinus
resonant charged current,
Definition: SREnums.h:78
caf::kResNCNuBarNeutronPiMinus
resonant charged current,
Definition: SREnums.h:84
caf::kCCQE
charged current quasi-elastic
Definition: SREnums.h:69
caf::kNuElectronElastic
neutrino electron elastic scatter
Definition: SREnums.h:113
caf::kResCCNuNeutronPiPlus
resonant charged current,
Definition: SREnums.h:73
caf::kInverseMuDecay
inverse muon decay
Definition: SREnums.h:114
caf::kResNCNuNeutronPi0
resonant neutral current,
Definition: SREnums.h:76
caf::kResNCNuNeutronPiMinus
resonant neutral current,
Definition: SREnums.h:77
caf::kResNCNuBarNeutronPi0
resonant charged current,
Definition: SREnums.h:83
enum caf::mode_type_

Neutrino interaction categories.

Enumerator
kUnknownMode 
kQE 
kRes 
kDIS 
kCoh 
kCohElastic 
kElectronScattering 
kIMDAnnihilation 
kInverseBetaDecay 
kGlashowResonance 
kAMNuGamma 
kMEC 
kDiffractive 
kEM 
kWeakMix 

Definition at line 47 of file SREnums.h.

47  {
48  kUnknownMode = -1,
49  kQE = 0,
50  kRes = 1,
51  kDIS = 2,
52  kCoh = 3,
53  kCohElastic = 4,
54  kElectronScattering = 5,
55  kIMDAnnihilation = 6,
56  kInverseBetaDecay = 7,
57  kGlashowResonance = 8,
58  kAMNuGamma = 9,
59  kMEC = 10,
60  kDiffractive = 11,
61  kEM = 12,
62  kWeakMix = 13
63  };
enum caf::View_t

Detector view, following GeometryObjects/PlaneGeo.h.

Enumerator
kX 

Vertical planes which measure X.

kY 

Horizontal planes which measure Y.

kXorY 

X or Y views.

Definition at line 20 of file SREnums.h.

21  {
22  kX, ///< Vertical planes which measure X
23  kY, ///< Horizontal planes which measure Y
24  kXorY ///< X or Y views
25  };
kXorY
Definition: geo.h:1
caf::kY
Horizontal planes which measure Y.
Definition: SREnums.h:23
kX
Definition: geo.h:1

Function Documentation

void caf::AddCosmicTruthToVec ( const cheat::NeutrinoEffPur *  effPur,
std::vector< cheat::TrackIDE > &  allTracks,
std::vector< cheat::TrackIDE > &  sliceTracks,
std::vector< cheat::TrackIDE > &  allTracksBirks,
std::vector< cheat::TrackIDE > &  sliceTracksBirks,
std::vector< SRCosmic > *  vec 
)

Definition at line 535 of file FillTruth.cxx.

References std::atan2(), caf::SRCosmic::azimuth, geo::GeometryBase::DetHalfHeight(), geo::GeometryBase::DetHalfWidth(), geo::GeometryBase::DetLength(), dir, caf::SRCosmic::E, caf::SRCosmic::eff, cheat::NeutrinoEffPur::efficiency, cheat::NeutrinoEffPur::energySlice, cheat::NeutrinoEffPur::energyTotal, caf::SRCosmic::enter, caf::SRCosmic::exit, FindAndAddMichels(), geom(), simb::MCTruth::GetParticle(), cheat::BackTracker::MCTruthToParticles(), caf::SRCosmic::michel, simb::MCTrajectory::Momentum(), simb::MCParticle::Momentum(), cheat::NeutrinoEffPur::neutrinoInt, caf::SRCosmic::nhitslc, caf::SRCosmic::nhittot, cheat::NeutrinoEffPur::nSliceHits, cheat::NeutrinoEffPur::nTotalHits, caf::SRCosmic::p, caf::SRCosmic::pdg, simb::MCParticle::PdgCode(), caf::SRCosmic::penter, simb::MCTrajectory::Position(), simb::MCParticle::Position(), geo::ProjectToBoxEdge(), gen_hdf5record::pt, caf::SRCosmic::pur, cheat::NeutrinoEffPur::purity, caf::SRVector3D::SetXYZ(), simb::MCTrajectory::size(), caf::SRCosmic::stop, caf::SRCosmic::time, Unit(), Vect(), caf::SRCosmic::visE, caf::SRCosmic::visEinslc, caf::SRCosmic::vtx, and caf::SRCosmic::zenith.

Referenced by caf::CAFMaker::produce().

541  {
542  static art::ServiceHandle<cheat::BackTracker> bt;
543  static art::ServiceHandle<geo::Geometry> geom;
544 
545  // Get information about true cosmic
546  vec->push_back(SRCosmic());
547  SRCosmic& srTruth = vec->back();
548 
549  const art::Ptr<simb::MCTruth> &tru = effPur->neutrinoInt;
550 
551  // BUT GetParticle(int i) requires looping over list of particles in
552  // single MCTruth object. So we have made an MCTruth object for each
553  // "cosmic" and GetParticle(0) corresponds to the only particle then
554  srTruth.pdg = tru->GetParticle(0).PdgCode();
555 
556  srTruth.E = tru->GetParticle(0).Momentum().E();
557  srTruth.visE = effPur->energyTotal;
558  srTruth.visEinslc = effPur->energySlice;
559  srTruth.eff = effPur->efficiency;
560  srTruth.pur = effPur->purity;
561  srTruth.nhitslc = effPur->nSliceHits;
562  srTruth.nhittot = effPur->nTotalHits;
563 
564  srTruth.time = tru->GetParticle(0).Position().T();
565  srTruth.vtx = tru->GetParticle(0).Position().Vect();
566  srTruth.p = tru->GetParticle(0).Momentum();
567 
568 
569  srTruth.azimuth = atan2(tru->GetParticle(0).Momentum().Pz(),
570  tru->GetParticle(0).Momentum().Px());
571  srTruth.zenith = -(tru->GetParticle(0).Momentum().Py())/tru->GetParticle(0).Momentum().P();
572 
573  std::vector<const sim::Particle*> particles = bt->MCTruthToParticles(tru);
574 
575 
576  const simb::MCTrajectory& traj = particles[0]->Trajectory();
577 
578 
579  bool entered = 0;
580  bool exited = 0;
581 
582  for(size_t pt = 0; pt < traj.size(); ++pt){
583  TVector3 point = traj.Position(pt).Vect();
584 
585  // Check to see if trajectory point is in the detector
586  bool inDet = fabs(point.X()) <= geom->DetHalfWidth() // Check X inside detector
587  && fabs(point.Y()) <= geom->DetHalfHeight() // then check Y
588  && point.Z() >= 0 && point.Z() <= geom->DetLength(); // then check Z
589 
590  if(inDet && !entered){
591  // Then we just entered the detector
592  entered = true;
593 
594  if(pt > 0){
595  // Put xyz from last point into an array for geo function call
596  TVector3 lastPoint = traj.Position(pt - 1).Vect();
597  double pointArray[3] = {point.X(), point.Y(), point.Z()};
598 
599  // Put xyz from direction vector into an array.
600  TVector3 dir = (lastPoint - point).Unit();
601  double dirArray[3] = {dir.X(), dir.Y(), dir.Z()};
602 
603  double enterArray[3] = {0,0,0};
604 
605  geo::ProjectToBoxEdge(pointArray, dirArray,
606  - geom->DetHalfWidth() , geom->DetHalfWidth(),
607  - geom->DetHalfHeight() , geom->DetHalfHeight(),
608  0 , geom->DetLength(),
609  enterArray);
610 
611  srTruth.enter.SetXYZ(enterArray[0], enterArray[1], enterArray[2]);
612  // No good way to interpolate between momentum at pt-1 and pt.
613  srTruth.penter = traj.Momentum(pt-1);
614  }
615  else{ // If this is the case, we started inside the detector.
616  srTruth.enter = point;
617  srTruth.penter = traj.Momentum(pt);
618  }
619  }
620 
621  if(!inDet && entered && !exited){
622  // Then we just exited the detector
623  exited= true;
624 
625  // Put xyz from last point into an array for geo function call
626  TVector3 lastPoint = traj.Position(pt - 1).Vect();
627  double lastPointArray[3] = {lastPoint.X(), lastPoint.Y(), lastPoint.Z()};
628 
629  // Put xyz from direction vector into an array.
630  TVector3 dir = (point - lastPoint).Unit();
631  double dirArray[3] = {dir.X(), dir.Y(), dir.Z()};
632 
633  double exitArray[3] = {0,0,0};
634 
635  geo::ProjectToBoxEdge(lastPointArray, dirArray,
636  - geom->DetHalfWidth() , geom->DetHalfWidth(),
637  - geom->DetHalfHeight() , geom->DetHalfHeight(),
638  0 , geom->DetLength(),
639  exitArray);
640 
641  srTruth.exit.SetXYZ(exitArray[0], exitArray[1], exitArray[2]);
642  }
643  }
644  srTruth.stop = traj.Position(traj.size() - 1).Vect();
645 
646  // Make sure the particle actually entered the detector
647  if(!entered){
648  // If the primary never entered the detector, set the
649  // coordinates to -5*10^9
650 
651  srTruth.enter.SetXYZ(-5e9, -5e9, -5e9);
652  srTruth.exit.SetXYZ(-5e9, -5e9, -5e9);
653  }
654 
655  // If the particle entered and did not exit, set the exit position to stop
656  if(entered && !exited){
657  srTruth.exit = srTruth.stop;
658  }
659  // Find the Michel electrons
660  FindAndAddMichels(particles, allTracks, &srTruth.michel);
661  }
simb::MCParticle::Position
const TLorentzVector & Position(const int i=0) const
Definition: MCParticle.h:218
simb::MCParticle::PdgCode
int PdgCode() const
Definition: MCParticle.h:211
geo::GeometryBase::DetLength
double DetLength() const
Definition: GeometryBase.cxx:481
cheat::NeutrinoEffPur::nTotalHits
int nTotalHits
Total number of hits the neutrino left in the detector.
Definition: BackTracker.h:56
cheat::NeutrinoEffPur::energyTotal
double energyTotal
Sum of FLS hits from the neutrino contributing to any hit in the event.
Definition: BackTracker.h:54
cheat::NeutrinoEffPur::efficiency
double efficiency
Efficiency (based on FLS energy) of neutrino interaction relative to slice.
Definition: BackTracker.h:50
caf::FindAndAddMichels
void FindAndAddMichels(std::vector< const sim::Particle * > particles, std::vector< cheat::TrackIDE > &allTracks, std::vector< SRTrueMichelE > *michelVec)
Definition: FillTruth.cxx:664
geo::GeometryBase::DetHalfHeight
double DetHalfHeight() const
Definition: GeometryBase.cxx:474
Unit
TVector3 Unit() const
Definition: SRVector3DAddons.h:14
simb::MCTrajectory::Position
const TLorentzVector & Position(const size_type) const
The accessor methods described above.
Definition: MCTrajectory.cxx:36
gen_hdf5record.pt
pt
Definition: gen_hdf5record.py:385
simb::MCTruth::GetParticle
const simb::MCParticle & GetParticle(int i) const
Definition: MCTruth.h:76
cheat::NeutrinoEffPur::nSliceHits
int nSliceHits
Number of hits from this neutrino in this slice.
Definition: BackTracker.h:55
geo::GeometryBase::DetHalfWidth
double DetHalfWidth() const
Definition: GeometryBase.cxx:467
simb::MCTrajectory::size
size_type size() const
Definition: MCTrajectory.h:165
dir
TDirectory * dir
Definition: macro.C:5
cheat::NeutrinoEffPur::energySlice
double energySlice
Sum of FLS hits from the neutrino contributing to hits included in the slice.
Definition: BackTracker.h:53
simb::MCParticle::Momentum
const TLorentzVector & Momentum(const int i=0) const
Definition: MCParticle.h:219
geom
void geom(int which=0)
Definition: geom.C:163
geo::ProjectToBoxEdge
void ProjectToBoxEdge(const double xyz[], const double dxyz[], double xlo, double xhi, double ylo, double yhi, double zlo, double zhi, double xyzout[])
Project along a direction from a particular starting point to the edge of a box.
Definition: Geo.cxx:38
simb::MCTrajectory::Momentum
const TLorentzVector & Momentum(const size_type) const
Definition: MCTrajectory.cxx:44
Vect
TVector3 Vect() const
Definition: SRLorentzVectorAddons.h:14
cheat::NeutrinoEffPur::neutrinoInt
art::Ptr< simb::MCTruth > neutrinoInt
Truth information about neutrino interaction.
Definition: BackTracker.h:49
std::atan2
T atan2(T number)
Definition: d0nt_math.hpp:72
cheat::BackTracker::MCTruthToParticles
std::vector< const sim::Particle * > MCTruthToParticles(art::Ptr< simb::MCTruth > const &mct) const
Definition: BackTracker_service.cc:534
cheat::NeutrinoEffPur::purity
double purity
Purity (based on FLS energy) of neutrino interaction relative to slice.
Definition: BackTracker.h:51
void caf::AddParticleToVec ( const sim::Particle &  part,
std::vector< cheat::TrackIDE > &  allTracks,
std::vector< cheat::TrackIDE > &  sliceTracks,
std::vector< cheat::TrackIDE > &  allTracksBirks,
std::vector< cheat::TrackIDE > &  sliceTracksBirks,
std::vector< SRTrueParticle > *  vec,
const std::vector< sim::TrueEnergy > &  TrueEnergies 
)

Definition at line 258 of file FillTruth.cxx.

References simb::MCParticle::Daughter(), caf::SRTrueParticle::daughterEnergies, caf::SRTrueParticle::daughterlist, caf::SRTrueParticle::daughterVisE, caf::SRTrueParticle::daughterVisEBirks, caf::SRTrueParticle::daughterVisEinslc, caf::SRTrueParticle::daughterVisEinslcBirks, simb::MCParticle::E(), caf::SRTrueParticle::elasticProtonMaxVisE, caf::SRTrueParticle::elasticProtonMaxVisEinslc, caf::SRTrueParticle::elasticProtonSumVisE, caf::SRTrueParticle::elasticProtonSumVisEinslc, caf::SRTrueParticle::enteringE, FindDaughterVisE(), FindDaughterVisENonRecur(), MECModelEnuComparisons::i, caf::SRTrueParticle::inelasticPhotonMaxVisE, caf::SRTrueParticle::inelasticPhotonMaxVisEinslc, caf::SRTrueParticle::inelasticPhotonSumVisE, caf::SRTrueParticle::inelasticPhotonSumVisEinslc, caf::SRTrueParticle::inelasticProtonMaxVisE, caf::SRTrueParticle::inelasticProtonMaxVisEinslc, caf::SRTrueParticle::inelasticProtonSumVisE, caf::SRTrueParticle::inelasticProtonSumVisEinslc, it, caf::SRTrueParticle::maxElasticProtonTrueE, caf::SRTrueParticle::maxInelasticPhotonTrueE, caf::SRTrueParticle::maxInelasticProtonTrueE, simb::MCParticle::Momentum(), simb::MCParticle::NumberDaughters(), caf::SRTrueParticle::p, caf::SRTrueParticle::pdg, genie::utils::res::PdgCode(), simb::MCParticle::PdgCode(), simb::MCParticle::Process(), caf::SRTrueParticle::totEscE, simb::MCParticle::TrackId(), cheat::BackTracker::TrackIDToParticle(), caf::SRTrueParticle::trkID, caf::SRTrueParticle::visE, caf::SRTrueParticle::visEBirks, caf::SRTrueParticle::visEinslc, and caf::SRTrueParticle::visEinslcBirks.

Referenced by caf::CAFMaker::AddMCTruthToVec().

265  {
266  static art::ServiceHandle<cheat::BackTracker> bt;
267  int trackId = part.TrackId();
268  // Find visible energy
269  double visE = 0;
270  for (std::vector<cheat::TrackIDE>::iterator it = allTracks.begin();
271  it != allTracks.end(); ++it)
272  {
273  if(trackId == it->trackID)
274  {
275  visE += it->energy;
276  break;
277  }
278 
279  }
280  // Find visible energy in slice
281  double inSliceVisE = 0;
282  for (std::vector<cheat::TrackIDE>::iterator it = sliceTracks.begin();
283  it != sliceTracks.end(); ++it)
284  {
285  if(trackId == it->trackID)
286  {
287  inSliceVisE += it->energy;
288  break;
289  }
290  }
291  //now calculate visible energies with birks suppression
292  // Find visible energy
293  double visEBirks = 0;
294  for (std::vector<cheat::TrackIDE>::iterator it = allTracksBirks.begin();
295  it != allTracksBirks.end(); ++it)
296  {
297  if(trackId == it->trackID)
298  {
299  visEBirks += it->energy;
300  break;
301  }
302 
303  }
304  // Find visible energy in slice
305  double inSliceVisEBirks = 0;
306  for (std::vector<cheat::TrackIDE>::iterator it = sliceTracksBirks.begin();
307  it != sliceTracksBirks.end(); ++it)
308  {
309  if(trackId == it->trackID)
310  {
311  inSliceVisEBirks += it->energy;
312  break;
313  }
314  }
315 
316 
317  double daughterVisE = 0;
318  double daughterInSliceVisE = 0;
319  double daughterVisEBirks = 0;
320  double daughterInSliceVisEBirks = 0;
321 
322 
323 
324  for (int i=0; i<part.NumberDaughters();++i){
325  const sim::Particle * p = bt->TrackIDToParticle(part.Daughter(i));
326 
327  daughterVisE += FindDaughterVisE(*p,allTracks);
328  daughterInSliceVisE += FindDaughterVisE(*p,sliceTracks);
329  daughterVisEBirks += FindDaughterVisE(*p,allTracksBirks);
330  daughterInSliceVisEBirks += FindDaughterVisE(*p,sliceTracksBirks);
331  }
332 
333 
334  // Find the total escaping energy in slice.
335  float TotEscE = -5;
336  float EnterEn = -5;
337  for(const sim::TrueEnergy& trueE: TrueEnergies){
338  if ( trackId == trueE.TrackId() ) {
339  EnterEn = trueE.EnteringEnergy();
340  TotEscE = trueE.TotalEscEnergy();
341  break;
342  }
343  }
344 
345  // Make a new SRTrueParticle and get an accessor
346  vec->push_back(SRTrueParticle());
347  SRTrueParticle& srPart = vec->back();
348  // Set fields
349  srPart.pdg = part.PdgCode();
350  srPart.visE = visE;
351  srPart.visEinslc = inSliceVisE;
352  srPart.daughterVisE = daughterVisE;
353  srPart.daughterVisEinslc = daughterInSliceVisE;
354  srPart.visEBirks = visEBirks;
355  srPart.visEinslcBirks = inSliceVisEBirks;
356  srPart.daughterVisEBirks = daughterVisEBirks;
357  srPart.daughterVisEinslcBirks = daughterInSliceVisEBirks;
358  srPart.totEscE = TotEscE;
359  srPart.enteringE = EnterEn;
360  srPart.p = part.Momentum();
361  srPart.trkID = trackId;
362 
363  float maxProInel = 0;
364  float maxProInelS =0;
365  float maxProElS =0;
366  float maxProEl =0;
367  float maxPhoInel = 0;
368  float maxPhoInelS = 0;
369 
370  float maxProInelTrueE = 0;
371  float maxProElTrueE = 0;
372  float maxPhoInelTrueE =0;
373 
374  double inelasticProtonSum = 0;
375  double inelasticPhotonSum = 0 ;
376  double elasticProtonSum = 0;
377  double inelasticProtonInslcSum = 0;
378  double inelasticPhotonInslcSum = 0 ;
379  double elasticProtonInslcSum = 0;
380 
381 
382  double inelasticProtonMax = 0;
383  double inelasticPhotonMax = 0 ;
384  double elasticProtonMax = 0;
385  double inelasticProtonInslcMax = 0;
386  double inelasticPhotonInslcMax = 0 ;
387  double elasticProtonInslcMax = 0;
388 
389 
390  for (int i=0; i<part.NumberDaughters();++i){
391  const sim::Particle * p = bt->TrackIDToParticle(part.Daughter(i));
392  srPart.daughterlist.push_back(bt->TrackIDToParticle(part.Daughter(i))->PdgCode());
393  srPart.daughterEnergies.push_back(p->E());
394 
395 
396  if(p->Process() == "neutronInelastic"){
397  if(p->PdgCode() ==2212){
398 
399  inelasticProtonSum += FindDaughterVisENonRecur(*p,allTracks);
400  inelasticProtonInslcSum += FindDaughterVisENonRecur(*p,sliceTracks);
401 
402  if( FindDaughterVisENonRecur(*p,allTracks) > maxProInel){
403  inelasticProtonMax = FindDaughterVisENonRecur(*p,allTracks);
404  maxProInel = inelasticProtonMax;
405 
406  }
407  if( FindDaughterVisENonRecur(*p,sliceTracks) > maxProInelS){
408  inelasticProtonInslcMax = FindDaughterVisENonRecur(*p,sliceTracks);
409  maxProInelS = inelasticProtonInslcMax;
410  maxProInelTrueE = p->E();
411  }
412 
413 
414  }
415  if(p->PdgCode() ==22){
416 
417  inelasticPhotonSum += FindDaughterVisENonRecur(*p,allTracks);
418  inelasticPhotonInslcSum += FindDaughterVisENonRecur(*p,sliceTracks);
419 
420  if( FindDaughterVisENonRecur(*p,allTracks) > maxPhoInel){
421  inelasticPhotonMax = FindDaughterVisENonRecur(*p,allTracks);
422  maxPhoInel = inelasticPhotonMax;
423  }
424  if( FindDaughterVisENonRecur(*p,sliceTracks) > maxPhoInelS){
425  inelasticPhotonInslcMax = FindDaughterVisENonRecur(*p,sliceTracks);
426  maxPhoInelS = inelasticPhotonInslcMax;
427  maxPhoInelTrueE = p->E();
428  }
429 
430 
431  }
432 
433  }
434  if(p->Process() == "hadElastic"){
435 
436  if(p->PdgCode() == 2212){
437 
438  elasticProtonSum += FindDaughterVisENonRecur(*p,allTracks);
439  elasticProtonInslcSum += FindDaughterVisENonRecur(*p,sliceTracks);
440 
441  if( FindDaughterVisENonRecur(*p,sliceTracks) > maxProElS){
442  elasticProtonInslcMax = FindDaughterVisENonRecur(*p,sliceTracks);
443  maxProElS = elasticProtonInslcMax;
444  maxProElTrueE = p->E();
445  }
446  if( FindDaughterVisENonRecur(*p,allTracks) > maxProEl){
447  elasticProtonMax = FindDaughterVisENonRecur(*p,allTracks);
448  maxProEl = elasticProtonMax;
449  }
450 
451  }
452 
453  }
454 
455  }
456 
457  srPart.inelasticProtonSumVisE = inelasticProtonSum;
458  srPart.inelasticPhotonSumVisE = inelasticPhotonSum;
459  srPart.elasticProtonSumVisE = elasticProtonSum;
460 
461 
462  srPart.inelasticProtonSumVisEinslc = inelasticProtonInslcSum;
463  srPart.inelasticPhotonSumVisEinslc = inelasticPhotonInslcSum;
464  srPart.elasticProtonSumVisEinslc = elasticProtonInslcSum;
465 
466 
467  srPart.inelasticProtonMaxVisE = inelasticProtonMax;
468  srPart.inelasticPhotonMaxVisE = inelasticPhotonMax;
469  srPart.elasticProtonMaxVisE = elasticProtonMax;
470 
471 
472  srPart.inelasticProtonMaxVisEinslc = inelasticProtonInslcMax;
473  srPart.inelasticPhotonMaxVisEinslc = inelasticPhotonInslcMax;
474  srPart.elasticProtonMaxVisEinslc = elasticProtonInslcMax;
475 
476 
477  srPart.maxElasticProtonTrueE = maxProElTrueE;
478  srPart.maxInelasticProtonTrueE = maxProInelTrueE;
479  srPart.maxInelasticPhotonTrueE = maxPhoInelTrueE;
480 
481 
482 
483 
484  }
simb::MCParticle::E
double E(const int i=0) const
Definition: MCParticle.h:232
simb::MCParticle::PdgCode
int PdgCode() const
Definition: MCParticle.h:211
caf::FindDaughterVisE
float FindDaughterVisE(const sim::Particle &part, std::vector< cheat::TrackIDE > &tracks)
Definition: FillTruth.cxx:487
it
set< int >::iterator it
Definition: ProduceBadRunsByFailureMode.C:101
p
const char * p
Definition: xmltok.h:285
simb::MCParticle::Process
std::string Process() const
Definition: MCParticle.h:214
simb::MCParticle::NumberDaughters
int NumberDaughters() const
Definition: MCParticle.h:216
caf::FindDaughterVisENonRecur
float FindDaughterVisENonRecur(const sim::Particle &part, std::vector< cheat::TrackIDE > &tracks)
Definition: FillTruth.cxx:512
simb::MCParticle::TrackId
int TrackId() const
Definition: MCParticle.h:209
simb::MCParticle::Daughter
int Daughter(const int i) const
Definition: MCParticle.cxx:112
genie::utils::res::PdgCode
int PdgCode(Resonance_t res, int Q)
(resonance id, charge) -> PDG code
Definition: BaryonResUtils.cxx:187
cheat::BackTracker::TrackIDToParticle
const sim::Particle * TrackIDToParticle(int const &id) const
Returns a pointer to the sim::Particle object corresponding to the given TrackID. ...
Definition: BackTracker_service.cc:488
simb::MCParticle::Momentum
const TLorentzVector & Momentum(const int i=0) const
Definition: MCParticle.h:219
void caf::AddPreFSI ( const art::Ptr< simb::MCTruth > &  truth,
SRNeutrino &  nu 
)

Definition at line 726 of file FillTruth.cxx.

References om::cout, e, allTimeWatchdog::endl, simb::MCTruth::GetParticle(), makeTrainCVSamples::int, simb::MCParticle::Momentum(), simb::MCParticle::Mother(), genie::utils::fragmrec::NParticles(), simb::MCTruth::NParticles(), caf::SRNuGenParticle::p, caf::SRNuGenParticle::pdg, period_run_lists::pfd, caf::SRNeutrino::prefsi, caf::SRNeutrino::prim, caf::SRNuGenParticle::primList, caf::SRNuGenParticle::status, and simb::MCParticle::StatusCode().

Referenced by caf::CAFMaker::AddMCTruthToVec().

727  {
728  // Add information about "pre-FSI hadrons", which are recoil hadrons inside the nucleus after
729  // resonance decays and hadronization and before FSI. This excludes hadrons from hadronization
730  // and resonance decays occuring outside the nucleus, as well as hadrons from coherent scattering
731  // and scattering on hydrogen.
732 
733  // Also included is the mapping between pre-FSI hadrons and their final state descendants.
734  // Final state particle information in the CAFs comes from sim::Particle objects from the
735  // BackTracker service and is stored in SRTrueParticle objects in SRNeutrino. Information of
736  // final state descendants of pre-FSI hadrons is stored in simb::MCParticle objects from
737  // simb::MCTruth. This situation requires matching the pre-FSI hadron final state descendant
738  // simb::MCParticle objects to the corresponding sim::Particle or SRTrueParticle objects by
739  // pdg code and 4-momentum.
740 
741  // First, retrieve the pre-FSI hadrons and their final state descendants
742 
743  std::map<const simb::MCParticle*,std::vector<const simb::MCParticle*>> prefsi_final_descendants;
744 
745  for ( int iprefsi = 0; iprefsi < truth->NParticles(); ++iprefsi )
746  {
747  const simb::MCParticle* prefsi_part = &truth->GetParticle( iprefsi );
748  // Require hadron undergoing intranuclear transport (status code = 14)
749  if ( prefsi_part->StatusCode() != 14 ) continue;
750  // Require hadron not a product of FSI
751  if ( truth->GetParticle( prefsi_part->Mother() ).StatusCode() == 14 ) continue;
752  // Find final state descendants
753  std::vector<const simb::MCParticle*> final_descendants;
754  for ( int ifinal = 0; ifinal < truth->NParticles(); ++ifinal )
755  {
756  const simb::MCParticle* final_part = &truth->GetParticle( ifinal );
757  // Require final state particle (particle (status code = 1)
758  if ( final_part->StatusCode() != 1 ) continue;
759  // Exclude at-rest final state particles produced by GENIE presumably for baryon accounting.
760  // Iddentical paricles with zero momentum are not unique and cannot be matched as described above.
761  if ( final_part->Momentum().Vect().Mag() == 0 ) continue;
762  int iancestor = final_part->Mother();
763  int search_count = 0;
764  while ( iancestor != iprefsi && iancestor > 0 && search_count < truth->NParticles() )
765  {
766  iancestor = truth->GetParticle( iancestor ).Mother();
767  search_count++;
768  }
769  if ( iancestor == iprefsi ) final_descendants.push_back( final_part );
770  }
771  prefsi_final_descendants[ prefsi_part ] = final_descendants;
772  }
773 
774  // Determine mapping between pre-FSI hadrons and their final state descendants in the nu.prim
775  // vector and add the pre-FSI hadron information to the record.
776 
777  for ( auto pfd : prefsi_final_descendants )
778  {
779  SRNuGenParticle nugen_part;
780  nugen_part.pdg = pfd.first->PdgCode();
781  nugen_part.status = pfd.first->StatusCode();
782  nugen_part.p = pfd.first->Momentum();
783 
784  const int prefsiID = (int)nu.prefsi.size();
785 
786  for ( unsigned int iprim = 0; iprim < nu.prim.size(); ++iprim )
787  {
788  for ( auto final_part : pfd.second )
789  {
790  double max_diff = 1e-4;
791  if ( final_part->PdgCode() == nu.prim[ iprim ].pdg &&
792  fabs( final_part->Px() - nu.prim[ iprim ].p.px ) < max_diff &&
793  fabs( final_part->Py() - nu.prim[ iprim ].p.py ) < max_diff &&
794  fabs( final_part->Pz() - nu.prim[ iprim ].p.pz ) < max_diff &&
795  fabs( final_part->E() - nu.prim[ iprim ].p.E ) < max_diff )
796  {
797  nugen_part.primList.push_back( iprim );
798  nu.prim[ iprim ].prefsiID = prefsiID;
799  }
800  }
801  }
802 
803  nu.prefsi.push_back( nugen_part );
804 
805  int nfinal = pfd.second.size();
806  int nprim = nu.prefsi[ prefsiID ].primList.size();
807  if ( nfinal != nprim )
808  {
809  std::cout << Form( "Mismatch: pre-FSI N final state descendants = %d, N primary matches = %d", nfinal, nprim ) << std::endl;
810  //if ( nfinal != nprim ) LOG_WARNING( "FillTruth" ) << Form( "Mismatch: pre-FSI N final state descendants = %d, N primary matches = %d", nfinal, nprim );
811  }
812  }
813 
814  }
simb::MCParticle::E
double E(const int i=0) const
Definition: MCParticle.h:232
simb::MCParticle::PdgCode
int PdgCode() const
Definition: MCParticle.h:211
simb::MCParticle::Py
double Py(const int i=0) const
Definition: MCParticle.h:230
simb::MCParticle::Mother
int Mother() const
Definition: MCParticle.h:212
simb::MCParticle::Px
double Px(const int i=0) const
Definition: MCParticle.h:229
simb::MCParticle::StatusCode
int StatusCode() const
Definition: MCParticle.h:210
simb::MCTruth::NParticles
int NParticles() const
Definition: MCTruth.h:75
om::cout
OStream cout
Definition: OStream.cxx:6
period_run_lists.pfd
pfd
Definition: period_run_lists.py:53
simb::MCTruth::GetParticle
const simb::MCParticle & GetParticle(int i) const
Definition: MCTruth.h:76
simb::MCParticle::Momentum
const TLorentzVector & Momentum(const int i=0) const
Definition: MCParticle.h:219
simb::MCParticle::Pz
double Pz(const int i=0) const
Definition: MCParticle.h:231
e
Float_t e
Definition: plot.C:35
genie::utils::fragmrec::NParticles
int NParticles(int pdg_code, const TClonesArray *const particle_list)
Definition: FragmRecUtils.cxx:32
void caf::AddSlcMEToVec ( const me::SlcME &  michel,
std::vector< SRSlcME > *  vec,
const std::vector< rb::Cluster > &  sliceList,
std::vector< cheat::TrackIDE > &  allTracks,
std::vector< cheat::TrackIDE > &  sliceTracks,
std::vector< cheat::TrackIDE > &  allTracksBirks,
std::vector< cheat::TrackIDE > &  sliceTracksBirks,
int  sliceIdx 
)

Definition at line 148 of file FillReco.cxx.

References caf::SRSlcME::adc, rb::Cluster::AllCells(), caf::SRSlcME::calE, rb::Cluster::CalorimetricEnergy(), me::SlcME::DeltaT(), caf::SRSlcME::deltat, me::SlcME::DistToSlc(), caf::SRSlcME::disttoslc, FillParticleTruth(), caf::SRSlcME::meanpos, rb::Cluster::MeanX(), rb::Cluster::MeanY(), rb::Cluster::MeanZ(), me::SlcME::MID(), caf::SRSlcME::mid, caf::SRSlcME::nhitx, caf::SRSlcME::nhity, rb::Cluster::NXCell(), rb::Cluster::NYCell(), rb::Cluster::TotalADC(), and caf::SRSlcME::truth.

Referenced by caf::CAFMaker::produce().

156  {
157  vec->push_back(SRSlcME());
158  SRSlcME& srME = vec->back();
159 
160  srME.mid = michel.MID();
161  srME.nhitx = michel.NXCell();
162  srME.nhity = michel.NYCell();
163  srME.adc = michel.TotalADC();
164  srME.calE = michel.CalorimetricEnergy();
165  srME.deltat = michel.DeltaT();
166  srME.disttoslc = michel.DistToSlc();
167  srME.meanpos = TVector3(michel.MeanX(),michel.MeanY(),michel.MeanZ());
168 
169  srME.truth = FillParticleTruth(sliceList,
170  michel.AllCells(),
171  allTracks,
172  sliceTracks,
173  allTracksBirks,
174  sliceTracksBirks,
175  sliceIdx);
176  }
me::SlcME::DistToSlc
double DistToSlc() const
Definition: MEClusters.cxx:12
rb::Cluster::MeanZ
double MeanZ(rb::AveragingScheme scheme=kDefaultScheme) const
Definition: Cluster.h:233
rb::Cluster::AllCells
art::PtrVector< rb::CellHit > AllCells() const
Get all cells from both views.
Definition: Cluster.cxx:180
rb::Cluster::MeanX
double MeanX(rb::AveragingScheme scheme=kDefaultScheme) const
Definition: Cluster.h:231
caf::FillParticleTruth
SRParticleTruth FillParticleTruth(const std::vector< rb::Cluster > &sliceList, const art::PtrVector< rb::CellHit > &hits, std::vector< cheat::TrackIDE > &allTracks, std::vector< cheat::TrackIDE > &sliceTracks, std::vector< cheat::TrackIDE > &allTracksBirks, std::vector< cheat::TrackIDE > &sliceTracksBirks, int sliceIdx)
Function to fill the particle truth for a set of hits.
Definition: FillTruth.cxx:28
rb::Cluster::TotalADC
double TotalADC() const
Sum of the ADC of all the contained hits.
Definition: Cluster.cxx:360
me::SlcME::DeltaT
double DeltaT() const
Definition: MEClusters.cxx:10
rb::Cluster::MeanY
double MeanY(rb::AveragingScheme scheme=kDefaultScheme) const
Definition: Cluster.h:232
rb::Cluster::CalorimetricEnergy
double CalorimetricEnergy(EEnergyCalcScheme escheme=kRecomputeEnergy) const
Simple estimate of neutrino energy.
Definition: Cluster.cxx:439
me::SlcME::MID
double MID() const
Definition: MEClusters.cxx:8
rb::Cluster::NYCell
unsigned int NYCell() const
Number of cells in the y-view.
Definition: Cluster.h:108
rb::Cluster::NXCell
unsigned int NXCell() const
Number of cells in the x-view.
Definition: Cluster.h:106
void caf::AddTrkMEToVec ( const me::TrkME &  michel,
std::vector< SRTrkME > *  vec,
SRTrack &  srTrk,
const std::vector< rb::Cluster > &  sliceList,
std::vector< cheat::TrackIDE > &  allTracks,
std::vector< cheat::TrackIDE > &  sliceTracks,
std::vector< cheat::TrackIDE > &  allTracksBirks,
std::vector< cheat::TrackIDE > &  sliceTracksBirks,
int  sliceIdx 
)

Definition at line 34 of file FillReco.cxx.

References caf::SRTrkME::adc, rb::Cluster::AllCells(), caf::SRTrkME::calE, rb::Cluster::CalorimetricEnergy(), me::SlcME::DeltaT(), caf::SRTrkME::deltat, caf::SRTrkME::disttotrack, me::TrkME::DistToTrk(), FillParticleTruth(), caf::SRTrack::me, caf::SRTrkME::meanpos, rb::Cluster::MeanX(), rb::Cluster::MeanY(), rb::Cluster::MeanZ(), me::SlcME::MID(), caf::SRTrkME::mid, caf::SRTrkME::nhitx, caf::SRTrkME::nhity, rb::Cluster::NXCell(), rb::Cluster::NYCell(), rb::Cluster::TotalADC(), and caf::SRTrkME::truth.

Referenced by caf::CAFMaker::produce().

43  {
44  vec->push_back(SRTrkME());
45  SRTrkME& srME = vec->back();
46 
47  srME.nhitx = michel.NXCell();
48  srME.nhity = michel.NYCell();
49  srME.mid = michel.MID();
50  srME.calE = michel.CalorimetricEnergy();
51  srME.adc = michel.TotalADC();
52  srME.deltat = michel.DeltaT();
53  srME.disttotrack = michel.DistToTrk();
54  srME.meanpos = TVector3(michel.MeanX(),michel.MeanY(),michel.MeanZ());
55  srME.truth = FillParticleTruth(sliceList,
56  michel.AllCells(),
57  allTracks,
58  sliceTracks,
59  allTracksBirks,
60  sliceTracksBirks,
61  sliceIdx);
62 
63  srTrk.me.push_back(srME); // Also copy this TrkME to parent track
64  }
rb::Cluster::MeanZ
double MeanZ(rb::AveragingScheme scheme=kDefaultScheme) const
Definition: Cluster.h:233
rb::Cluster::AllCells
art::PtrVector< rb::CellHit > AllCells() const
Get all cells from both views.
Definition: Cluster.cxx:180
rb::Cluster::MeanX
double MeanX(rb::AveragingScheme scheme=kDefaultScheme) const
Definition: Cluster.h:231
caf::FillParticleTruth
SRParticleTruth FillParticleTruth(const std::vector< rb::Cluster > &sliceList, const art::PtrVector< rb::CellHit > &hits, std::vector< cheat::TrackIDE > &allTracks, std::vector< cheat::TrackIDE > &sliceTracks, std::vector< cheat::TrackIDE > &allTracksBirks, std::vector< cheat::TrackIDE > &sliceTracksBirks, int sliceIdx)
Function to fill the particle truth for a set of hits.
Definition: FillTruth.cxx:28
rb::Cluster::TotalADC
double TotalADC() const
Sum of the ADC of all the contained hits.
Definition: Cluster.cxx:360
me::SlcME::DeltaT
double DeltaT() const
Definition: MEClusters.cxx:10
rb::Cluster::MeanY
double MeanY(rb::AveragingScheme scheme=kDefaultScheme) const
Definition: Cluster.h:232
rb::Cluster::CalorimetricEnergy
double CalorimetricEnergy(EEnergyCalcScheme escheme=kRecomputeEnergy) const
Simple estimate of neutrino energy.
Definition: Cluster.cxx:439
me::SlcME::MID
double MID() const
Definition: MEClusters.cxx:8
rb::Cluster::NYCell
unsigned int NYCell() const
Number of cells in the y-view.
Definition: Cluster.h:108
rb::Cluster::NXCell
unsigned int NXCell() const
Number of cells in the x-view.
Definition: Cluster.h:106
me::TrkME::DistToTrk
double DistToTrk() const
Definition: MEClusters.cxx:34
void caf::AddTrueMichelEToVec ( const sim::Particle &  michel,
const sim::Particle &  motherMuon,
std::vector< cheat::TrackIDE > &  allTracks,
std::vector< SRTrueMichelE > *  vec 
)

Definition at line 698 of file FillTruth.cxx.

References caf::SRTrueMichelE::E, simb::MCParticle::E(), simb::MCParticle::EndPosition(), it, simb::MCParticle::Momentum(), caf::SRTrueMichelE::mustop, caf::SRTrueMichelE::p, simb::MCParticle::T(), caf::SRTrueMichelE::time, simb::MCParticle::TrackId(), and caf::SRTrueMichelE::visE.

Referenced by FindAndAddMichels().

702  {
703  vec->push_back(SRTrueMichelE());
704  SRTrueMichelE& srMichel = vec->back();
705 
706  // Find the visible energy from the FLS hits
707  float visE = 0;
708  for (std::vector<cheat::TrackIDE>::iterator it = allTracks.begin();
709  it != allTracks.end(); ++it)
710  {
711  if(michel.TrackId() == it->trackID)
712  {
713  visE = it->energy;
714  break;
715  }
716  }
717 
718  srMichel.E = michel.E();
719  srMichel.visE = visE;
720  srMichel.time = michel.T();
721  srMichel.mustop = motherMuon.EndPosition();
722  srMichel.p = michel.Momentum();
723  }
simb::MCParticle::E
double E(const int i=0) const
Definition: MCParticle.h:232
it
set< int >::iterator it
Definition: ProduceBadRunsByFailureMode.C:101
simb::MCParticle::EndPosition
const TLorentzVector & EndPosition() const
Definition: MCParticle.h:224
simb::MCParticle::TrackId
int TrackId() const
Definition: MCParticle.h:209
simb::MCParticle::T
double T(const int i=0) const
Definition: MCParticle.h:223
simb::MCParticle::Momentum
const TLorentzVector & Momentum(const int i=0) const
Definition: MCParticle.h:219
const long& caf::AdjustBase ( const long &  base,
CAFType  type,
const std::string &  name 
)

Definition at line 94 of file BasicTypesProxy.cxx.

References kFlatSingleTree, kZero, and NSubscripts().

95  {
96  if(type == kFlatSingleTree && NSubscripts(name) == 0) return kZero;
97  return base;
98  }
name
const XML_Char * name
Definition: expat.h:151
base
const XML_Char int const XML_Char int const XML_Char * base
Definition: expat.h:331
caf::kZero
const long kZero
Definition: BasicTypesProxy.cxx:87
caf::NSubscripts
int NSubscripts(const std::string &name)
Count the subscripts in the name.
Definition: BasicTypesProxy.cxx:81
void caf::BlindThisCosRej ( SRCosRej &  cosrej)

Definition at line 156 of file Blinding.cxx.

References blind_val, caf::SRCosRej::eratio, and caf::SRCosRej::hadE.

Referenced by BlindThisRecord().

156  {
157  cosrej.eratio = blind_val;
158  cosrej.hadE = blind_val;
159  }
caf::blind_val
static const int blind_val
Definition: Blinding.cxx:16
cosrej
Cosmic Rejection PIDs for Numu analysis.
Definition: FillParentInfo.h:9
void caf::BlindThisElasticVtx ( SRElastic &  vtx)

Definition at line 58 of file Blinding.cxx.

References blind_val, blind_val_unsigned, caf::SRElastic::fuzzyk, MECModelEnuComparisons::i, caf::SRFuzzyK::png, and caf::SRFuzzyK::png2d.

Referenced by BlindThisRecord().

58  {
59  for (unsigned int i=0; i<vtx.fuzzyk.png.size(); i++) {
60  vtx.fuzzyk.png[i].nhit = blind_val_unsigned;
61  vtx.fuzzyk.png[i].nhitx = blind_val_unsigned;
62  vtx.fuzzyk.png[i].nhity = blind_val_unsigned;
63  vtx.fuzzyk.png[i].nplane = blind_val_unsigned;
64  vtx.fuzzyk.png[i].maxplanecont = blind_val_unsigned;
65  vtx.fuzzyk.png[i].calE = blind_val;
66  vtx.fuzzyk.png[i].weightedCalE = blind_val;
67  vtx.fuzzyk.png[i].len = blind_val;
68 
69  vtx.fuzzyk.png[i].shwlid.nhit = blind_val_unsigned;
70  vtx.fuzzyk.png[i].shwlid.nhitx = blind_val_unsigned;
71  vtx.fuzzyk.png[i].shwlid.nhity = blind_val_unsigned;
72  vtx.fuzzyk.png[i].shwlid.nplane = blind_val_unsigned;
73  vtx.fuzzyk.png[i].shwlid.maxplanecont = blind_val_unsigned;
74  vtx.fuzzyk.png[i].shwlid.calE = blind_val;
75  vtx.fuzzyk.png[i].shwlid.len = blind_val;
76  vtx.fuzzyk.png[i].shwlid.stop.SetZ (blind_val);
77  vtx.fuzzyk.png[i].shwlid.lid.setDefault();
78  vtx.fuzzyk.png[i].shwlid.lidE.setDefault();
79  vtx.fuzzyk.png[i].cvnpart.setDefault();
80  vtx.fuzzyk.png[i].regcvn.setDefault();
81 
82  vtx.fuzzyk.png[i].bpf.muon.nhit = blind_val_unsigned;
83  vtx.fuzzyk.png[i].bpf.muon.nhitx = blind_val_unsigned;
84  vtx.fuzzyk.png[i].bpf.muon.nhity = blind_val_unsigned;
85  vtx.fuzzyk.png[i].bpf.muon.nplane = blind_val_unsigned;
86  vtx.fuzzyk.png[i].bpf.muon.maxplanecont = blind_val_unsigned;
87  vtx.fuzzyk.png[i].bpf.muon.calE = blind_val;
88  vtx.fuzzyk.png[i].bpf.muon.len = blind_val;
89  vtx.fuzzyk.png[i].bpf.muon.lenE = blind_val;
90  vtx.fuzzyk.png[i].bpf.muon.stop. SetZ( blind_val );
91  vtx.fuzzyk.png[i].bpf.muon.pid = blind_val;
92  vtx.fuzzyk.png[i].bpf.muon.chi2T = blind_val;
93  vtx.fuzzyk.png[i].bpf.muon.dEdXLL = blind_val;
94  vtx.fuzzyk.png[i].bpf.muon.hitRatio = blind_val;
95  vtx.fuzzyk.png[i].bpf.muon.energy = blind_val;
96  vtx.fuzzyk.png[i].bpf.muon.lenE = blind_val;
97  vtx.fuzzyk.png[i].bpf.muon.overlapE = blind_val;
98  vtx.fuzzyk.png[i].bpf.muon.momentum. SetX( blind_val );
99  vtx.fuzzyk.png[i].bpf.muon.momentum. SetY( blind_val );
100  vtx.fuzzyk.png[i].bpf.muon.momentum. SetZ( blind_val );
101 
102  vtx.fuzzyk.png[i].bpf.pion.nhit = blind_val_unsigned;
103  vtx.fuzzyk.png[i].bpf.pion.nhitx = blind_val_unsigned;
104  vtx.fuzzyk.png[i].bpf.pion.nhity = blind_val_unsigned;
105  vtx.fuzzyk.png[i].bpf.pion.nplane = blind_val_unsigned;
106  vtx.fuzzyk.png[i].bpf.pion.maxplanecont = blind_val_unsigned;
107  vtx.fuzzyk.png[i].bpf.pion.calE = blind_val;
108  vtx.fuzzyk.png[i].bpf.pion.len = blind_val;
109  vtx.fuzzyk.png[i].bpf.pion.lenE = blind_val;
110  vtx.fuzzyk.png[i].bpf.pion.stop. SetZ( blind_val );
111  vtx.fuzzyk.png[i].bpf.pion.pid = blind_val;
112  vtx.fuzzyk.png[i].bpf.pion.chi2T = blind_val;
113  vtx.fuzzyk.png[i].bpf.pion.dEdXLL = blind_val;
114  vtx.fuzzyk.png[i].bpf.pion.hitRatio = blind_val;
115  vtx.fuzzyk.png[i].bpf.pion.energy = blind_val;
116  vtx.fuzzyk.png[i].bpf.pion.lenE = blind_val;
117  vtx.fuzzyk.png[i].bpf.pion.overlapE = blind_val;
118  vtx.fuzzyk.png[i].bpf.pion.momentum. SetX( blind_val );
119  vtx.fuzzyk.png[i].bpf.pion.momentum. SetY( blind_val );
120  vtx.fuzzyk.png[i].bpf.pion.momentum. SetZ( blind_val );
121 
122  vtx.fuzzyk.png[i].bpf.proton.nhit = blind_val_unsigned;
123  vtx.fuzzyk.png[i].bpf.proton.nhitx = blind_val_unsigned;
124  vtx.fuzzyk.png[i].bpf.proton.nhity = blind_val_unsigned;
125  vtx.fuzzyk.png[i].bpf.proton.nplane = blind_val_unsigned;
126  vtx.fuzzyk.png[i].bpf.proton.maxplanecont = blind_val_unsigned;
127  vtx.fuzzyk.png[i].bpf.proton.calE = blind_val;
128  vtx.fuzzyk.png[i].bpf.proton.len = blind_val;
129  vtx.fuzzyk.png[i].bpf.proton.lenE = blind_val;
130  vtx.fuzzyk.png[i].bpf.proton.stop. SetZ( blind_val );
131  vtx.fuzzyk.png[i].bpf.proton.pid = blind_val;
132  vtx.fuzzyk.png[i].bpf.proton.chi2T = blind_val;
133  vtx.fuzzyk.png[i].bpf.proton.dEdXLL = blind_val;
134  vtx.fuzzyk.png[i].bpf.proton.hitRatio = blind_val;
135  vtx.fuzzyk.png[i].bpf.proton.energy = blind_val;
136  vtx.fuzzyk.png[i].bpf.proton.lenE = blind_val;
137  vtx.fuzzyk.png[i].bpf.proton.overlapE = blind_val;
138  vtx.fuzzyk.png[i].bpf.proton.momentum. SetX( blind_val );
139  vtx.fuzzyk.png[i].bpf.proton.momentum. SetY( blind_val );
140  vtx.fuzzyk.png[i].bpf.proton.momentum. SetZ( blind_val );
141 
142  }
143  for (unsigned int i=0; i<vtx.fuzzyk.png2d.size(); i++) {
144  vtx.fuzzyk.png2d[i].nhit = blind_val_unsigned;
145  vtx.fuzzyk.png2d[i].nhitx = blind_val_unsigned;
146  vtx.fuzzyk.png2d[i].nhity = blind_val_unsigned;
147  vtx.fuzzyk.png2d[i].nplane = blind_val_unsigned;
148  vtx.fuzzyk.png2d[i].maxplanecont = blind_val_unsigned;
149  vtx.fuzzyk.png2d[i].calE = blind_val;
150  vtx.fuzzyk.png2d[i].weightedCalE = blind_val;
151  vtx.fuzzyk.png2d[i].len = blind_val;
152  }
153  }
caf::blind_val
static const int blind_val
Definition: Blinding.cxx:16
caf::blind_val_unsigned
static const unsigned int blind_val_unsigned
Definition: Blinding.cxx:17
void caf::BlindThisJMShower ( SRJMShower &  jm)

Definition at line 44 of file Blinding.cxx.

References blind_val, blind_val_unsigned, caf::SRJMShower::calE, caf::SRJMShower::len, caf::SRJMShower::maxplanecont, caf::SRJMShower::nhit, caf::SRJMShower::nhitx, caf::SRJMShower::nhity, caf::SRJMShower::nplane, caf::SRJMShower::shwE, caf::SRJMShower::stop, and caf::SRJMShower::vtxE.

44  {
45  jm.nhit = blind_val_unsigned;
46  jm.nhitx = blind_val_unsigned;
47  jm.nhity = blind_val_unsigned;
48  jm.nplane = blind_val_unsigned;
49  jm.maxplanecont = blind_val_unsigned;
50  jm.calE = blind_val;
51  jm.len = blind_val;
52  jm.shwE = blind_val;
53  jm.vtxE = blind_val;
54  jm.stop. SetZ( blind_val );
55  }
caf::blind_val
static const int blind_val
Definition: Blinding.cxx:16
caf::blind_val_unsigned
static const unsigned int blind_val_unsigned
Definition: Blinding.cxx:17
void caf::BlindThisRecord ( StandardRecord *  rec)

If required, scrub the record of any quantities that shouldn't be viewed

Definition at line 162 of file Blinding.cxx.

References caf::SRHeader::blind, blind_val, BlindThisCosRej(), BlindThisElasticVtx(), BlindThisSlice(), BlindThisTrack(), caf::SRIDBranch::bpfid, caf::SRSlice::calE, om::cerr, caf::SRTrackBranch::cosmic, caf::SRIDBranch::cosrej, caf::SRIDBranch::cvnloosepreselptp, caf::SRIDBranch::cvnnocosmics, caf::SRIDBranch::cvnoldpresel, caf::SRHeader::det, caf::SRTrackBranch::discrete, caf::SRVertexBranch::elastic, allTimeWatchdog::endl, caf::StandardRecord::energy, caf::StandardRecord::hdr, MECModelEnuComparisons::i, caf::SRHeader::ismc, caf::SRElastic::IsValid, caf::SRTrackBranch::kalman, kFARDET, caf::SRIDBranch::lem, caf::SRNueEnergy::lid, caf::SRIDBranch::lid, caf::SRSlice::meantime, caf::SREnergyBranch::nue, caf::SREnergyBranch::numu, caf::SRIDBranch::remid, caf::SRIDBranch::rvp, caf::StandardRecord::sel, caf::SRSLidEnergy::setDefault(), caf::SRNueEnergy::setDefault(), caf::SRCVNResult::setDefault(), caf::SRELid::setDefault(), caf::SRRemid::setDefault(), caf::SRRvp::setDefault(), caf::SRBpfId::setDefault(), caf::SRNumuEnergy::setDefault(), caf::SRLem::setDefault(), caf::StandardRecord::slc, caf::SRTrackBase::tracks, caf::SRKalman::tracks, caf::StandardRecord::trk, caf::StandardRecord::vtx, and caf::SRTrackBranch::window.

Referenced by caf::CAFMaker::produce().

162  {
163  /// If required, scrub the record of any quantities that shouldn't be viewed
164 
165  if(rec->hdr.det != kFARDET){
166  std::cerr << "Should only apply blinding to the Far Detector, aborting."
167  << std::endl;
168  abort();
169  }
170 
171  if(rec->hdr.ismc){
172  std::cerr << "Shouldn't apply blinding to Monte Carlo, aborting."
173  << std::endl;
174  abort();
175  }
176 
177  // Is this record in the box?
178 
179  // Blinding 0.8 to 3 GeV events within 30 microsecond window
180  // around center of beam spill (t=222000 microseconds).
181  if (rec->slc.calE > 0.8 &&
182  rec->slc.calE < 3.0 &&
183  rec->slc.meantime > 207000 &&
184  rec->slc.meantime < 237000 ) {
185 
186  // Walk through each rec branch.
187 
188  // -- start of branches --
189 
190  // -- hdr
191  // Tag the record as blind
192  rec->hdr.blind = 1;
193 
194  // -- spill
195 
196  // -- slc
197  BlindThisSlice(rec->slc);
198 
199  // -- trk
200  for (unsigned int i=0; i<rec->trk.kalman.tracks.size(); i++) {
201  BlindThisTrack(rec->trk.kalman.tracks[i]);
202  rec->trk.kalman.tracks[i].remavededx = blind_val;
203  rec->trk.kalman.tracks[i].rempid = blind_val;
204  rec->trk.kalman.tracks[i].scatllh = blind_val;
205  rec->trk.kalman.tracks[i].dedxllh = blind_val;
206  rec->trk.kalman.tracks[i].measfrac = blind_val;
207  rec->trk.kalman.tracks[i].remcont = blind_val;
208  rec->trk.kalman.tracks[i].lenE = blind_val;
209  rec->trk.kalman.tracks[i].overlapE = blind_val;
210  }
211  for (unsigned int i=0; i<rec->trk.discrete.tracks.size(); i++) {
212  BlindThisTrack(rec->trk.discrete.tracks[i]);
213  }
214  for (unsigned int i=0; i<rec->trk.cosmic.tracks.size(); i++) {
215  BlindThisTrack(rec->trk.cosmic.tracks[i]);
216  }
217  for (unsigned int i=0; i<rec->trk.window.tracks.size(); i++) {
218  BlindThisTrack(rec->trk.window.tracks[i]);
219  }
220 
221  // -- vtx
222  if (rec->vtx.elastic.IsValid == 1) {
223  BlindThisElasticVtx(rec->vtx.elastic);
224  }
225 
226  // -- energy
227  // for now, everything goes, so use setDefault() methods
228  rec->energy.numu.setDefault();
229  rec->energy.nue .setDefault();
230  rec->energy.nue.lid.setDefault();
231 
232  // for (unsigned int i=0; i<rec->energy.nue.lid.size(); i++) {
233  // rec->energy.nue.lid[i].setDefault();
234  // }
235 
236  // -- sel
237  // for now, everything goes, so use setDefault() methods
238  rec->sel.remid.setDefault();
239  rec->sel.bpfid.setDefault();
240  rec->sel.lid .setDefault();
241  // for (unsigned int i=0; i<rec->sel.lid.shwlid.size(); i++){
242  // rec->sel.lid.shwlid[i].setDefault();
243  //}
244  rec->sel.lem .setDefault();
245  rec->sel.rvp .setDefault();
246  rec->sel.cvnloosepreselptp.setDefault();
247  rec->sel.cvnoldpresel .setDefault();
248  rec->sel.cvnnocosmics .setDefault();
249 
250  // some stuff survives here:
251  BlindThisCosRej (rec->sel.cosrej);
252 
253  // -- mc
254 
255  // -- parent
256 
257  // -- end of branches --
258 
259  } // end if (rec needs blinding)
260 
261  }
caf::SRHeader::det
Det_t det
Detector, ND = 1, FD = 2, NDOS = 3.
Definition: SRHeader.h:28
caf::SRNueEnergy::setDefault
virtual void setDefault()
Definition: SRNueEnergy.cxx:23
caf::SRBpfId::setDefault
virtual void setDefault()
Definition: SRBpfId.cxx:35
caf::SRIDBranch::cvnnocosmics
SRCVNResult cvnnocosmics
Output from CVN - No cosmics ued in training (many-class PID)
Definition: SRIDBranch.h:56
caf::SRHeader::blind
bool blind
if true, record has been corrupted for blindness
Definition: SRHeader.h:29
caf::StandardRecord::hdr
SRHeader hdr
Header branch: run, subrun, etc.
Definition: StandardRecord.h:40
caf::SRIDBranch::cosrej
SRCosRej cosrej
Output from CosRej (Cosmic Rejection)
Definition: SRIDBranch.h:47
caf::blind_val
static const int blind_val
Definition: Blinding.cxx:16
caf::SRCVNResult::setDefault
virtual void setDefault()
Definition: SRCVNResult.cxx:28
om::cerr
OStream cerr
Definition: OStream.cxx:7
caf::SRTrackBranch::discrete
SRTrackBase discrete
3D tracks produced by DiscreteTrack
Definition: SRTrackBranch.h:25
rec
caf::StandardRecord * rec
Definition: tutCAFMacro.C:20
caf::SREnergyBranch::nue
SRNueEnergy nue
Nue energy variables.
Definition: SREnergyBranch.h:26
caf::SRIDBranch::cvnoldpresel
SRCVNResult cvnoldpresel
Output from CVN - Preselection used in Prod3/4 (many-class PID)
Definition: SRIDBranch.h:55
caf::SRIDBranch::cvnloosepreselptp
SRCVNResult cvnloosepreselptp
Output from CVN - Loose Presel plus PtP cut (many-class PID)
Definition: SRIDBranch.h:54
caf::SRHeader::ismc
bool ismc
data or MC? True if MC
Definition: SRHeader.h:27
caf::SRELid::setDefault
void setDefault()
Definition: SRELid.cxx:27
caf::SRIDBranch::rvp
SRRvp rvp
Output from RecoVariablePID (RVP)
Definition: SRIDBranch.h:46
caf::SRRemid::setDefault
virtual void setDefault()
Definition: SRRemid.cxx:29
caf::BlindThisCosRej
void BlindThisCosRej(SRCosRej &cosrej)
Definition: Blinding.cxx:156
caf::BlindThisTrack
void BlindThisTrack(SRTrack &trk)
Definition: Blinding.cxx:31
caf::SRNumuEnergy::setDefault
virtual void setDefault()
Definition: SRNumuEnergy.cxx:61
caf::SRTrackBranch::kalman
SRKalman kalman
Tracks produced by KalmanTrack.
Definition: SRTrackBranch.h:24
caf::SRIDBranch::remid
SRRemid remid
Output from RecoMuonID (ReMId) package.
Definition: SRIDBranch.h:39
caf::SRSlice::calE
float calE
Calorimetric energy of the cluster [GeV].
Definition: SRSlice.h:38
caf::BlindThisSlice
void BlindThisSlice(SRSlice &slc)
Definition: Blinding.cxx:20
caf::SRSLidEnergy::setDefault
virtual void setDefault()
Definition: SRSLidEnergy.cxx:27
caf::SRIDBranch::bpfid
SRBpfId bpfid
Output from the BreakPointFitter PID (BPFPIdMaker) package.
Definition: SRIDBranch.h:41
caf::BlindThisElasticVtx
void BlindThisElasticVtx(SRElastic &vtx)
Definition: Blinding.cxx:58
caf::SREnergyBranch::numu
SRNumuEnergy numu
Numu energy estimator.
Definition: SREnergyBranch.h:27
caf::SRTrackBranch::window
SRTrackBase window
Tracks produced by WindowTrack.
Definition: SRTrackBranch.h:27
caf::SRIDBranch::lid
SRELid lid
Output from LIDBuilder (LID) package.
Definition: SRIDBranch.h:42
caf::StandardRecord::sel
SRIDBranch sel
Selector (PID) branch.
Definition: StandardRecord.h:47
caf::SRVertexBranch::elastic
SRElastic elastic
Single vertex found by Elastic Arms.
Definition: SRVertexBranch.h:26
caf::StandardRecord::slc
SRSlice slc
Slice branch: nhit, extents, time, etc.
Definition: StandardRecord.h:42
caf::SRLem::setDefault
void setDefault()
Definition: SRLem.cxx:58
caf::SRIDBranch::lem
SRLem lem
Output from Library Event Matching (LEM)
Definition: SRIDBranch.h:43
caf::SRNueEnergy::lid
SRSLidEnergy lid
Energy estimate from lid for slice pid (most energetic shower for now)
Definition: SRNueEnergy.h:23
caf::StandardRecord::trk
SRTrackBranch trk
Track branch: nhit, len, etc.
Definition: StandardRecord.h:43
caf::StandardRecord::energy
SREnergyBranch energy
Energy estimator branch.
Definition: StandardRecord.h:46
caf::SRTrackBranch::cosmic
SRTrackBase cosmic
Tracks produced by CosmicTrack.
Definition: SRTrackBranch.h:26
caf::SRTrackBase::tracks
std::vector< SRTrack > tracks
Definition: SRTrackBase.h:15
caf::SRKalman::tracks
std::vector< SRKalmanTrack > tracks
3D Tracks produced by KalmanTrack
Definition: SRKalman.h:16
caf::SRSlice::meantime
float meantime
mean time, weighted by charge [ns]
Definition: SRSlice.h:41
caf::SRRvp::setDefault
void setDefault()
Definition: SRRvp.cxx:40
caf::StandardRecord::vtx
SRVertexBranch vtx
Vertex branch: location, time, etc.
Definition: StandardRecord.h:44
void caf::BlindThisSlice ( SRSlice &  slc)

Definition at line 20 of file Blinding.cxx.

References blind_val, blind_val_unsigned, caf::SRSlice::boxmax, caf::SRSlice::calE, caf::SRSlice::lastplane, std::max(), caf::SRSlice::meanpos, caf::SRSlice::ncontplanes, caf::SRVector3D::SetZ(), and submit_hadd::u.

Referenced by BlindThisRecord().

20  {
21  // Small numbers of ncontplanes are useful for cosmic rejection, large
22  // numbers can leak the event length (~ energy).
23  slc.ncontplanes = std::max(10u, slc.ncontplanes);
24  slc.lastplane = blind_val_unsigned;
25  slc.calE = blind_val;
26  slc.boxmax. SetZ( blind_val );
27  slc.meanpos.SetZ( blind_val );
28  }
std::max
T max(const caf::Proxy< T > &a, T b)
Definition: BasicTypesProxy.h:415
caf::blind_val
static const int blind_val
Definition: Blinding.cxx:16
submit_hadd.u
u
Definition: submit_hadd.py:84
caf::blind_val_unsigned
static const unsigned int blind_val_unsigned
Definition: Blinding.cxx:17
PandAna.Demos.tute_pid_validation.slc
slc
Definition: tute_pid_validation.py:182
void caf::BlindThisTrack ( SRTrack &  trk)

Definition at line 31 of file Blinding.cxx.

References blind_val, blind_val_unsigned, caf::SRTrack::calE, caf::SRTrack::len, caf::SRTrack::lenE, caf::SRTrack::maxplanecont, caf::SRTrack::nhit, caf::SRTrack::nhitx, caf::SRTrack::nhity, caf::SRTrack::nplane, and caf::SRTrack::stop.

Referenced by BlindThisRecord().

31  {
32  trk.nhit = blind_val_unsigned;
33  trk.nhitx = blind_val_unsigned;
34  trk.nhity = blind_val_unsigned;
35  trk.nplane = blind_val_unsigned;
36  trk.maxplanecont = blind_val_unsigned;
37  trk.calE = blind_val;
38  trk.len = blind_val;
39  trk.lenE = blind_val;
40  trk.stop. SetZ( blind_val );
41  }
caf::blind_val
static const int blind_val
Definition: Blinding.cxx:16
trk
Track finder for cosmic rays.
Definition: CosmicTrackAna_module.cc:182
caf::blind_val_unsigned
static const unsigned int blind_val_unsigned
Definition: Blinding.cxx:17
caf::generator_ caf::CAFGeneratorEnum ( simb::Generator_t  simbGeneratorEnum)

Definition at line 905 of file FillTruth.cxx.

References simb::kGENIE, kGENIE, simb::kGIBUU, kGIBUU, and kUnknownGenerator.

Referenced by caf::CAFMaker::AddMCTruthToVec().

906  {
907  switch (simbGeneratorEnum)
908  {
909  case simb::Generator_t::kGENIE:
910  return kGENIE;
911 
912  case simb::Generator_t::kGIBUU:
913  return kGIBUU;
914 
915  default:
916  return kUnknownGenerator;
917 
918  }
919  }
void caf::CopyMuonIDVars ( const SRKalmanTrack &  srTrk,
caf::SRMuonID &  muid 
)

Definition at line 723 of file FillPIDs.cxx.

References caf::SRKalmanTrack::muonid, and caf::SRMuonID::pid.

Referenced by caf::CAFMaker::produce().

724  {
725  muid.pid = srTrk.muonid;
726  }
caf::SRMuonID::pid
float pid
MuonID value output by ND-trained BDT.
Definition: SRMuonID.h:18
void caf::CopyRemidVars ( const SRKalmanTrack &  srTrk,
caf::SRRemid &  remid 
)

Definition at line 712 of file FillPIDs.cxx.

References caf::SRRemid::cont, caf::SRKalmanTrack::dedxllh, caf::SRRemid::dedxllh, caf::SRRemid::len, caf::SRTrack::len, caf::SRKalmanTrack::measfrac, caf::SRRemid::measfrac, caf::SRRemid::pid, caf::SRKalmanTrack::remcont, caf::SRKalmanTrack::rempid, caf::SRKalmanTrack::scatllh, and caf::SRRemid::scatllh.

Referenced by caf::CAFMaker::produce().

713  {
714  remid.pid = srTrk.rempid;
715  remid.scatllh = srTrk.scatllh;
716  remid.dedxllh = srTrk.dedxllh;
717  remid.len = srTrk.len;
718  remid.measfrac = srTrk.measfrac;
719  remid.cont = srTrk.remcont;
720  }
caf::SRRemid::scatllh
float scatllh
Log-likelihood value from scattering angle.
Definition: SRRemid.h:26
caf::SRRemid::len
float len
Track length of identified track.
Definition: SRRemid.h:28
caf::SRRemid::cont
int cont
True if track is contained as defined by ReMId.
Definition: SRRemid.h:30
caf::SRRemid::pid
float pid
PID value output by kNN.
Definition: SRRemid.h:25
caf::SRRemid::measfrac
float measfrac
Fraction of planes used to measure dE/dx.
Definition: SRRemid.h:29
caf::SRRemid::dedxllh
float dedxllh
Log-likelihood value from dE/dx.
Definition: SRRemid.h:27
std::vector< unsigned int > caf::DecodeGeneratorVersion ( const std::string &  versionString)

Definition at line 922 of file FillTruth.cxx.

References it, submit_syst::pattern, submit_cafana::regex, and runNovaSAM::ret.

Referenced by caf::CAFMaker::AddMCTruthToVec().

923  {
924  std::vector<unsigned int> ret;
925  std::regex pattern("(\\d+)");
926  auto vals_begin = std::sregex_iterator(versionString.begin(), versionString.end(), pattern);
927  auto vals_end = std::sregex_iterator();
928 
929  for (auto it = vals_begin; it != vals_end; it++)
930  ret.push_back(std::stoi(it->str()));
931 
932  return ret;
933  }
it
set< int >::iterator it
Definition: ProduceBadRunsByFailureMode.C:101
submit_syst.pattern
pattern
Definition: submit_syst.py:246
submit_cafana.regex
regex
Definition: submit_cafana.py:299
template<class T >
void caf::EvalInstanceWrapper ( TTreeFormula *  ttf,
T &  x 
)

Definition at line 236 of file BasicTypesProxy.cxx.

References T.

Referenced by caf::Proxy< T >::GetValueNested().

237  {
238  // TODO is this the safest way to cast?
239  x = (T)ttf->EvalInstance(0);
240  }
submit_syst.x
x
Definition: submit_syst.py:249
T
double T
Definition: Xdiff_gwt.C:5
void caf::EvalInstanceWrapper ( TTreeFormula *  ttf,
std::string &  x 
)

Definition at line 242 of file BasicTypesProxy.cxx.

243  {
244  x = ttf->EvalStringInstance(0);
245  }
submit_syst.x
x
Definition: submit_syst.py:249
void caf::FillCosRejVars ( const cosrej::CosRejObj &  cosrej,
caf::SRCosRej &  srcosrej,
caf::SRContain &  srcontain 
)

Definition at line 394 of file FillPIDs.cxx.

References caf::SRCosRej::anglecos, cosrej::CosRejObj::AngleCos(), caf::SRCosRej::anglekal, cosrej::CosRejObj::AngleKal(), caf::SRCosRej::cdirscore, cosrej::CosRejObj::CDirScore(), caf::SRCosRej::cfitspeed, cosrej::CosRejObj::CFitSpeed(), cosrej::CosRejObj::ConCosPID(), caf::SRContain::cosbakair, cosrej::CosRejObj::CosBakAir(), caf::SRContain::cosbakcell, cosrej::CosRejObj::CosBakCell(), caf::SRContain::cosbakcellnd, cosrej::CosRejObj::CosBakCellND(), caf::SRContain::cosbakdist, cosrej::CosRejObj::CosBakDist(), caf::SRContain::cosbaksteel, cosrej::CosRejObj::CosBakSteel(), caf::SRCosRej::coschidiff, cosrej::CosRejObj::CosChiDiff(), caf::SRCosRej::coschisq, cosrej::CosRejObj::CosChisq(), caf::SRContain::cosfwdair, cosrej::CosRejObj::CosFwdAir(), caf::SRContain::cosfwdcell, cosrej::CosRejObj::CosFwdCell(), caf::SRContain::cosfwdcellnd, cosrej::CosRejObj::CosFwdCellND(), caf::SRContain::cosfwddist, cosrej::CosRejObj::CosFwdDist(), caf::SRContain::cosfwdsteel, cosrej::CosRejObj::CosFwdSteel(), caf::SRCosRej::cosslope, cosrej::CosRejObj::CosSlope(), caf::SRCosRej::costhetatrue, cosrej::CosRejObj::CosThetaTrue(), caf::SRContain::cosyposattrans, cosrej::CosRejObj::CosYPosAtTrans(), caf::SRCosRej::cscorediff, cosrej::CosRejObj::CScoreDiff(), caf::SRContain::enddist, cosrej::CosRejObj::EndDist(), caf::SRCosRej::eratio, cosrej::CosRejObj::ERatio(), caf::SRCosRej::fscattext, cosrej::CosRejObj::FScattExt(), caf::SRCosRej::fscattmax, cosrej::CosRejObj::FScattMax(), caf::SRCosRej::fscattsig, cosrej::CosRejObj::FScattSig(), caf::SRCosRej::fscattsum, cosrej::CosRejObj::FScattSum(), caf::SRCosRej::hadE, caf::SRContain::kalbakair, cosrej::CosRejObj::KalBakAir(), caf::SRContain::kalbakcell, cosrej::CosRejObj::KalBakCell(), caf::SRContain::kalbakcellnd, cosrej::CosRejObj::KalBakCellND(), caf::SRContain::kalbakdist, cosrej::CosRejObj::KalBakDist(), caf::SRContain::kalbaksteel, cosrej::CosRejObj::KalBakSteel(), caf::SRCosRej::kalchidiff, cosrej::CosRejObj::KalChiDiff(), caf::SRCosRej::kalchisq, cosrej::CosRejObj::KalChisq(), caf::SRContain::kalfwdair, cosrej::CosRejObj::KalFwdAir(), caf::SRContain::kalfwdcell, cosrej::CosRejObj::KalFwdCell(), caf::SRContain::kalfwdcellnd, cosrej::CosRejObj::KalFwdCellND(), caf::SRContain::kalfwddist, cosrej::CosRejObj::KalFwdDist(), caf::SRContain::kalfwdsteel, cosrej::CosRejObj::KalFwdSteel(), caf::SRCosRej::kalslope, cosrej::CosRejObj::KalSlope(), caf::SRCosRej::kalthetatrue, cosrej::CosRejObj::KalThetaTrue(), caf::SRContain::kalyposattrans, cosrej::CosRejObj::KalYPosAtTrans(), caf::SRCosRej::kdirscore, cosrej::CosRejObj::KDirScore(), caf::SRCosRej::kfitspeed, cosrej::CosRejObj::KFitSpeed(), caf::SRCosRej::kscorediff, cosrej::CosRejObj::KScoreDiff(), caf::SRCosRej::mincell, cosrej::CosRejObj::MinCell(), caf::SRCosRej::mindist, cosrej::CosRejObj::MinDist(), caf::SRCosRej::nkal3d, cosrej::CosRejObj::NTracks3D(), caf::SRCosRej::numucontpid2019, caf::SRCosRej::numucontpid2020, caf::SRCosRej::numuunconttunedpid, cosrej::CosRejObj::OldCosPID(), caf::SRCosRej::pdgbest, cosrej::CosRejObj::PDGBest(), caf::SRCosRej::scatt, cosrej::CosRejObj::Scatt(), cosrej::CosRejObj::StartAct(), cosrej::CosRejObj::UnconTunedCosPID(), caf::SRContain::vtxdist, and cosrej::CosRejObj::VtxDist().

Referenced by caf::CAFMaker::produce().

397  {
398  srcosrej.anglekal = cosrej.AngleKal();
399  srcosrej.anglecos = cosrej.AngleCos();
400 
401  srcosrej.numucontpid2020 = cosrej.ConCosPID();
402  srcosrej.numucontpid2019 = cosrej.OldCosPID();
403  srcosrej.numuunconttunedpid = cosrej.UnconTunedCosPID();
404 
405  srcosrej.nkal3d = cosrej.NTracks3D();
406  srcosrej.mindist = cosrej.MinDist();
407  srcosrej.mincell = cosrej.MinCell();
408  srcosrej.eratio = cosrej.ERatio();
409  srcosrej.scatt = cosrej.Scatt();
410  srcosrej.hadE = cosrej.StartAct();
411  srcosrej.costhetatrue = cosrej.CosThetaTrue();
412  srcosrej.kalthetatrue = cosrej.KalThetaTrue();
413  srcosrej.pdgbest = cosrej.PDGBest();
414  srcosrej.fscattmax = cosrej.FScattMax();
415  srcosrej.fscattsum = cosrej.FScattSum();
416  srcosrej.fscattext = cosrej.FScattExt();
417  srcosrej.fscattsig = cosrej.FScattSig();
418  srcosrej.kfitspeed = cosrej.KFitSpeed();
419  srcosrej.cfitspeed = cosrej.CFitSpeed();
420  srcosrej.kdirscore = cosrej.KDirScore();
421  srcosrej.cdirscore = cosrej.CDirScore();
422  srcosrej.kscorediff = cosrej.KScoreDiff();
423  srcosrej.cscorediff = cosrej.CScoreDiff();
424  srcosrej.kalchisq = cosrej.KalChisq();
425  srcosrej.kalslope = cosrej.KalSlope();
426  srcosrej.kalchidiff = cosrej.KalChiDiff();
427  srcosrej.coschisq = cosrej.CosChisq();
428  srcosrej.cosslope = cosrej.CosSlope();
429  srcosrej.coschidiff = cosrej.CosChiDiff();
430 
431  srcontain.vtxdist = cosrej.VtxDist();
432  srcontain.enddist = cosrej.EndDist();
433  srcontain.cosfwddist = cosrej.CosFwdDist();
434  srcontain.cosbakdist = cosrej.CosBakDist();
435  srcontain.cosfwdcell = cosrej.CosFwdCell();
436  srcontain.cosfwdcellnd = cosrej.CosFwdCellND();
437  srcontain.cosbakcell = cosrej.CosBakCell();
438  srcontain.cosbakcellnd = cosrej.CosBakCellND();
439  srcontain.cosfwdair = cosrej.CosFwdAir();
440  srcontain.cosfwdsteel = cosrej.CosFwdSteel();
441  srcontain.cosbakair = cosrej.CosBakAir();
442  srcontain.cosbaksteel = cosrej.CosBakSteel();
443  srcontain.kalfwdair = cosrej.KalFwdAir();
444  srcontain.kalfwdsteel = cosrej.KalFwdSteel();
445  srcontain.kalbakair = cosrej.KalBakAir();
446  srcontain.kalbaksteel = cosrej.KalBakSteel();
447  srcontain.kalfwddist = cosrej.KalFwdDist();
448  srcontain.kalbakdist = cosrej.KalBakDist();
449  srcontain.kalfwdcell = cosrej.KalFwdCell();
450  srcontain.kalfwdcellnd = cosrej.KalFwdCellND();
451  srcontain.kalbakcell = cosrej.KalBakCell();
452  srcontain.kalbakcellnd = cosrej.KalBakCellND();
453  srcontain.cosyposattrans = cosrej.CosYPosAtTrans();
454  srcontain.kalyposattrans = cosrej.KalYPosAtTrans();
455  }
caf::SRContain::cosfwdcell
int cosfwdcell
cosmictrack projected # cells from end point forwards to det edge
Definition: SRContain.h:30
cosrej::CosRejObj::PDGBest
int PDGBest() const
Definition: CosRejObj.cxx:736
cosrej::CosRejObj::KalThetaTrue
float KalThetaTrue() const
Definition: CosRejObj.cxx:730
cosrej::CosRejObj::KalChisq
float KalChisq() const
Definition: CosRejObj.cxx:538
caf::SRCosRej::fscattsum
float fscattsum
sum of scattering angles (Fernanda)
Definition: SRCosRej.h:43
cosrej::CosRejObj::CosSlope
float CosSlope() const
Definition: CosRejObj.cxx:544
cosrej::CosRejObj::KalSlope
float KalSlope() const
Definition: CosRejObj.cxx:532
caf::SRCosRej::kalchisq
float kalchisq
chisq value of slope fit (kalmantrack)
Definition: SRCosRej.h:62
cosrej::CosRejObj::CosChisq
float CosChisq() const
Definition: CosRejObj.cxx:550
cosrej::CosRejObj::CosFwdDist
float CosFwdDist() const
Definition: CosRejObj.cxx:556
cosrej::CosRejObj::KalYPosAtTrans
float KalYPosAtTrans() const
Definition: CosRejObj.cxx:682
cosrej::CosRejObj::KalBakDist
float KalBakDist() const
Definition: CosRejObj.cxx:610
caf::SRCosRej::costhetatrue
float costhetatrue
cosine of angle between cosmic trk dir and true dir of most contributing particle ...
Definition: SRCosRej.h:26
caf::SRCosRej::cfitspeed
float cfitspeed
fit inverse speed (ns/cm) of Hough timing fit (TimingFit module) for cosmic track ...
Definition: SRCosRej.h:56
cosrej::CosRejObj::KalBakSteel
float KalBakSteel() const
Definition: CosRejObj.cxx:622
cosrej::CosRejObj::UnconTunedCosPID
float UnconTunedCosPID() const
Definition: CosRejObj.cxx:478
cosrej::CosRejObj::KalFwdSteel
float KalFwdSteel() const
Definition: CosRejObj.cxx:604
caf::SRContain::cosyposattrans
float cosyposattrans
Y position of cosmic track at transition (ND only, use to check if went through air gap) ...
Definition: SRContain.h:49
cosrej::CosRejObj::StartAct
float StartAct() const
Definition: CosRejObj.cxx:712
caf::SRCosRej::cosslope
float cosslope
slope of track timing fit (cosmictrack)
Definition: SRCosRej.h:64
caf::SRContain::kalyposattrans
float kalyposattrans
Y position of Kalman track and transition (ND only, use to check if went through air gap) ...
Definition: SRContain.h:50
cosrej::CosRejObj::CosBakCellND
int CosBakCellND() const
Definition: CosRejObj.cxx:646
caf::SRCosRej::pdgbest
int pdgbest
pdg code of most contributing true particle to best kalman trk
Definition: SRCosRej.h:28
cosrej::CosRejObj::Scatt
float Scatt() const
Definition: CosRejObj.cxx:706
cosrej::CosRejObj::FScattMax
float FScattMax() const
Definition: CosRejObj.cxx:767
caf::SRCosRej::fscattsig
float fscattsig
sigma of scattering variable distribution (Fernanda)
Definition: SRCosRej.h:45
cosrej::CosRejObj::KalBakCellND
int KalBakCellND() const
Definition: CosRejObj.cxx:670
cosrej::CosRejObj::KalFwdDist
float KalFwdDist() const
Definition: CosRejObj.cxx:592
caf::SRContain::vtxdist
float vtxdist
shortest distance to wall from primary Kalman track start position (proxy for vertex) ...
Definition: SRContain.h:26
cosrej::CosRejObj::AngleCos
float AngleCos() const
Definition: CosRejObj.cxx:459
caf::SRContain::cosbakcell
int cosbakcell
cosmictrack projected # cells from start point backwards to det edge
Definition: SRContain.h:35
caf::SRContain::kalfwdair
float kalfwdair
for Kalmantrack projected distance forwards how much is through air (ND only, NYI) ...
Definition: SRContain.h:42
cosrej::CosRejObj::CosBakSteel
float CosBakSteel() const
Definition: CosRejObj.cxx:586
caf::SRCosRej::scatt
float scatt
sum of all all angular deviation (looping over kaltrk traj pts) / trk len
Definition: SRCosRej.h:41
caf::SRCosRej::cscorediff
float cscorediff
chisq score difference between best +c/-c assumption and free fit in Hough timing fit (TimingFit modu...
Definition: SRCosRej.h:60
cosrej::CosRejObj::KalFwdCell
int KalFwdCell() const
Definition: CosRejObj.cxx:652
caf::SRContain::kalbakcellnd
int kalbakcellnd
Kalmantrack projected # cells from start point backwards to det edge, including muon catcher...
Definition: SRContain.h:46
cosrej::CosRejObj::KalBakAir
float KalBakAir() const
Definition: CosRejObj.cxx:616
caf::SRCosRej::kalslope
float kalslope
slope of track timing fit (kalmantrack)
Definition: SRCosRej.h:61
caf::SRContain::kalfwdsteel
float kalfwdsteel
for Kalmantrack projected distance forwards, how much is through steel (ND only, currently is just di...
Definition: SRContain.h:43
caf::SRContain::kalfwddist
float kalfwddist
Kalmantrack projected distance (cm) from end point forwards to det edge.
Definition: SRContain.h:39
caf::SRCosRej::hadE
float hadE
GeV sum of non (Kalman best ReMId) track E.
Definition: SRCosRej.h:51
caf::SRCosRej::mindist
float mindist
minimum projected distance to edge of any Kalman track with > 15 hits
Definition: SRCosRej.h:17
caf::SRContain::cosbakdist
float cosbakdist
cosmictrack projected distance (cm) from start point backwards to det edge
Definition: SRContain.h:34
caf::SRContain::kalfwdcellnd
int kalfwdcellnd
Kalmantrack projected # cells from end point forwards to det edge, including muon catcher...
Definition: SRContain.h:41
cosrej::CosRejObj::CScoreDiff
float CScoreDiff() const
Definition: CosRejObj.cxx:520
cosrej::CosRejObj::CosBakAir
float CosBakAir() const
Definition: CosRejObj.cxx:580
caf::SRCosRej::anglekal
float anglekal
cos of angle of best ReMId Kalman track
Definition: SRCosRej.h:23
cosrej::CosRejObj::CosChiDiff
float CosChiDiff() const
Definition: CosRejObj.cxx:749
caf::SRCosRej::mincell
int mincell
minimum projected cells to edge of any Kalman track with > 15 hits
Definition: SRCosRej.h:18
cosrej::CosRejObj::CosBakDist
float CosBakDist() const
Definition: CosRejObj.cxx:574
cosrej::CosRejObj::FScattSig
float FScattSig() const
Definition: CosRejObj.cxx:785
caf::SRCosRej::kdirscore
float kdirscore
chisq score difference between +c and -c assumptions in Hough timing fit (TimingFit module) for best ...
Definition: SRCosRej.h:57
cosrej::CosRejObj::VtxDist
float VtxDist() const
Definition: CosRejObj.cxx:755
cosrej::CosRejObj::ConCosPID
float ConCosPID() const
Definition: CosRejObj.cxx:465
cosrej::CosRejObj::FScattSum
float FScattSum() const
Definition: CosRejObj.cxx:773
cosrej::CosRejObj::KDirScore
float KDirScore() const
Definition: CosRejObj.cxx:502
cosrej::CosRejObj::CosFwdAir
float CosFwdAir() const
Definition: CosRejObj.cxx:562
cosrej::CosRejObj::KalBakCell
int KalBakCell() const
Definition: CosRejObj.cxx:664
caf::SRContain::cosfwddist
float cosfwddist
cosmictrack projected distance (cm) from end point forwards to det edge
Definition: SRContain.h:29
caf::SRContain::cosfwdcellnd
int cosfwdcellnd
cosmictrack projected # cells from end point forwards to det edge, including muon catcher...
Definition: SRContain.h:31
cosrej::CosRejObj::CosThetaTrue
float CosThetaTrue() const
Definition: CosRejObj.cxx:724
cosrej::CosRejObj::CosBakCell
int CosBakCell() const
Definition: CosRejObj.cxx:640
cosrej::CosRejObj::CosFwdCell
int CosFwdCell() const
Definition: CosRejObj.cxx:628
caf::SRContain::cosfwdair
float cosfwdair
for cosmictrack projected distance forwards how much is through air (ND only, NYI) ...
Definition: SRContain.h:32
cosrej::CosRejObj::EndDist
float EndDist() const
Definition: CosRejObj.cxx:761
caf::SRContain::cosfwdsteel
float cosfwdsteel
for cosmictrack projected distance forwards, how much is through steel (ND only, currently is just di...
Definition: SRContain.h:33
cosrej::CosRejObj::CDirScore
float CDirScore() const
Definition: CosRejObj.cxx:508
caf::SRContain::kalbaksteel
float kalbaksteel
for Kalmantrack projected distance backwards, how much is through steel (ND only, currently is just d...
Definition: SRContain.h:48
cosrej::CosRejObj::MinCell
int MinCell() const
Definition: CosRejObj.cxx:688
cosrej::CosRejObj::CosFwdSteel
float CosFwdSteel() const
Definition: CosRejObj.cxx:568
caf::SRCosRej::kalchidiff
float kalchidiff
chisq difference of track timing fits (kalmantrack)
Definition: SRCosRej.h:63
cosrej::CosRejObj::KScoreDiff
float KScoreDiff() const
Definition: CosRejObj.cxx:514
cosrej::CosRejObj::FScattExt
float FScattExt() const
Definition: CosRejObj.cxx:779
caf::SRCosRej::coschidiff
float coschidiff
chisq difference of track timing fits (cosmictrack)
Definition: SRCosRej.h:66
caf::SRCosRej::coschisq
float coschisq
chisq value of slope fit (cosmictrack)
Definition: SRCosRej.h:65
caf::SRContain::cosbaksteel
float cosbaksteel
for cosmictrack projected distance backwards, how much is through steel (ND only, currently is just d...
Definition: SRContain.h:38
caf::SRContain::kalfwdcell
int kalfwdcell
Kalmantrack projected # cells from end point forwards to det edge.
Definition: SRContain.h:40
caf::SRContain::enddist
float enddist
shortest distance to wall from primary Kalman track end position
Definition: SRContain.h:27
cosrej::CosRejObj::NTracks3D
int NTracks3D() const
Definition: CosRejObj.cxx:484
cosrej::CosRejObj::CosYPosAtTrans
float CosYPosAtTrans() const
Definition: CosRejObj.cxx:676
caf::SRCosRej::nkal3d
int nkal3d
number of 3D tracks as determined by Kalman tracker
Definition: SRCosRej.h:25
cosrej::CosRejObj::CFitSpeed
float CFitSpeed() const
Definition: CosRejObj.cxx:496
caf::SRCosRej::eratio
float eratio
ratio of best kalman track GeV / slice GeV
Definition: SRCosRej.h:50
cosrej::CosRejObj::MinDist
float MinDist() const
Definition: CosRejObj.cxx:694
caf::SRContain::kalbakcell
int kalbakcell
Kalmantrack projected # cells from start point backwards to det edge.
Definition: SRContain.h:45
caf::SRCosRej::kalthetatrue
float kalthetatrue
cosine of angle between kalman trk dir and true dir of most contributing particle ...
Definition: SRCosRej.h:27
cosrej::CosRejObj::KalFwdCellND
int KalFwdCellND() const
Definition: CosRejObj.cxx:658
cosrej::CosRejObj::KFitSpeed
float KFitSpeed() const
Definition: CosRejObj.cxx:490
cosrej::CosRejObj::KalFwdAir
float KalFwdAir() const
Definition: CosRejObj.cxx:598
caf::SRCosRej::numuunconttunedpid
float numuunconttunedpid
tuned cosmic rejection PID for uncontained events - Jose&#39;s;
Definition: SRCosRej.h:36
cosrej::CosRejObj::CosFwdCellND
int CosFwdCellND() const
Definition: CosRejObj.cxx:634
cosrej::CosRejObj::ERatio
float ERatio() const
Definition: CosRejObj.cxx:700
caf::SRCosRej::numucontpid2019
float numucontpid2019
cosmic rejection PID for contained events; 2019 Analysis
Definition: SRCosRej.h:35
caf::SRContain::kalbakdist
float kalbakdist
Kalmantrack projected distance (cm) from start point backwards to det edge.
Definition: SRContain.h:44
caf::SRContain::cosbakair
float cosbakair
for cosmictrack projected distance backwards how much is through air (ND only, NYI) ...
Definition: SRContain.h:37
caf::SRCosRej::numucontpid2020
float numucontpid2020
cosmic rejection PID for contained events; 2020 Analysis
Definition: SRCosRej.h:34
caf::SRCosRej::kfitspeed
float kfitspeed
fit inverse speed (ns/cm) of Hough timing fit (TimingFit module) for best ReMId Kalman track ...
Definition: SRCosRej.h:55
caf::SRCosRej::anglecos
float anglecos
cos of angle of Cosmic Track
Definition: SRCosRej.h:24
caf::SRContain::cosbakcellnd
int cosbakcellnd
cosmictrack projected # cells from start point backwards to det edge, including muon catcher...
Definition: SRContain.h:36
cosrej::CosRejObj::KalChiDiff
float KalChiDiff() const
Definition: CosRejObj.cxx:743
caf::SRCosRej::kscorediff
float kscorediff
chisq score difference between best +c/-c assumption and free fit in Hough timing fit (TimingFit modu...
Definition: SRCosRej.h:59
caf::SRCosRej::fscattext
float fscattext
scattering variable (Fernanda)
Definition: SRCosRej.h:44
cosrej::CosRejObj::OldCosPID
float OldCosPID() const
Definition: CosRejObj.cxx:471
caf::SRCosRej::cdirscore
float cdirscore
chisq score difference between +c and -c assumptions in Hough timing fit (TimingFit module) for cosmi...
Definition: SRCosRej.h:58
caf::SRCosRej::fscattmax
float fscattmax
maximum scattering angle (Fernanda)
Definition: SRCosRej.h:42
cosrej::CosRejObj::AngleKal
float AngleKal() const
Definition: CosRejObj.cxx:453
caf::SRContain::kalbakair
float kalbakair
for Kalmantrack projected distance backwards how much is through air (ND only, NYI) ...
Definition: SRContain.h:47
void caf::FillCVNFeaturesVars ( const cvn::Features &  features,
caf::SRCVNFeatures &  srcvnfeatures,
int  maxcomponents 
)

Definition at line 699 of file FillPIDs.cxx.

References caf::SRCVNFeatures::components, cvn::Features::fOutput, cvn::Features::fPrincipalComponents, caf::SRCVNFeatures::ncomponents, caf::SRCVNFeatures::noutput, and caf::SRCVNFeatures::output.

702  {
703  srcvnfeatures.output = features.fOutput;
704  srcvnfeatures.noutput = features.fOutput.size();
705  srcvnfeatures.components = features.fPrincipalComponents;
706  if(int(srcvnfeatures.components.size()) > maxcomponents)
707  srcvnfeatures.components.resize(maxcomponents);
708  srcvnfeatures.ncomponents = srcvnfeatures.components.size();
709  }
caf::SRCVNFeatures::components
std::vector< float > components
Principal components of output.
Definition: SRCVNFeatures.h:18
cvn::Features::fOutput
std::vector< float > fOutput
Vector of outputs from neural net.
Definition: CVNFeatures.h:23
caf::SRCVNFeatures::output
std::vector< float > output
List of net output nodes.
Definition: SRCVNFeatures.h:15
caf::SRCVNFeatures::ncomponents
unsigned int ncomponents
Size of components vector.
Definition: SRCVNFeatures.h:19
caf::SRCVNFeatures::noutput
unsigned int noutput
Number of entries in output vector.
Definition: SRCVNFeatures.h:16
cvn::Features::fPrincipalComponents
std::vector< float > fPrincipalComponents
Linear combinations of fOutput.
Definition: CVNFeatures.h:25
void caf::FillCVNNeutronDaughterResultVars ( const std::vector< art::Ptr< rb::PID > > &  cvnneutrons,
caf::SRCVNNeutronDaughterResult &  cvnneutron 
)

Definition at line 95 of file FillPIDs.cxx.

References MECModelEnuComparisons::i, caf::SRCVNNeutronDaughterResult::otherid, caf::SRCVNNeutronDaughterResult::photonid, caf::SRCVNNeutronDaughterResult::protonid, and febshutoff_auto::val.

Referenced by caf::CAFMaker::produce().

97  {
98  for (unsigned int i=0; i<cvnneutrons.size(); ++i){
99  float val = cvnneutrons[i]->Value();
100  int id = cvnneutrons[i]->Pdg();
101  if(id==2212)
102  cvnneutron.protonid = val;
103  else if(id==22)
104  cvnneutron.photonid = val;
105  else
106  cvnneutron.otherid = val;
107  }
108 
109  }
febshutoff_auto.val
val
Definition: febshutoff_auto.py:277
void caf::FillCVNParticleResultVars ( const std::vector< art::Ptr< rb::PID > > &  cvnparts,
caf::SRCVNParticleResult &  cvnpart 
)

Definition at line 112 of file FillPIDs.cxx.

References caf::SRCVNParticleResult::electronid, caf::SRCVNParticleResult::emid, caf::SRCVNParticleResult::hadronid, MECModelEnuComparisons::i, caf::SRCVNParticleResult::maxval, caf::SRCVNParticleResult::muonid, caf::SRCVNParticleResult::neutronid, caf::SRCVNParticleResult::otherid, caf::SRCVNParticleResult::pdgmax, caf::SRCVNParticleResult::photonid, caf::SRCVNParticleResult::pionid, caf::SRCVNParticleResult::pizeroid, caf::SRCVNParticleResult::protonid, and febshutoff_auto::val.

Referenced by caf::CAFMaker::produce().

114  {
115 
116  float largestv = 0;
117  int largestid = 0;
118  for (unsigned int i=0; i<cvnparts.size(); ++i){
119  float val = cvnparts[i]->Value();
120  int id = cvnparts[i]->Pdg();
121  if (val > largestv){
122  largestv = val;
123  largestid = id;
124  }
125  switch (id) {
126  case 11:
127  cvnpart.electronid = val;
128  break;
129  case 13:
130  cvnpart.muonid = val;
131  break;
132  case 2212:
133  cvnpart.protonid = val;
134  break;
135  case 2112:
136  cvnpart.neutronid = val;
137  break;
138  case 211:
139  cvnpart.pionid = val;
140  break;
141  case 111:
142  cvnpart.pizeroid = val;
143  break;
144  case 22:
145  cvnpart.photonid = val;
146  break;
147  case 22 * 11:
148  cvnpart.emid = val;
149  break;
150  case 2212 * 211:
151  cvnpart.hadronid = val;
152  break;
153  default:
154  cvnpart.otherid = val;
155  break;
156  }
157  }
158  cvnpart.maxval = largestv;
159  cvnpart.pdgmax = largestid;
160 
161  if(cvnpart.emid < 0)
162  cvnpart.emid = cvnpart.electronid + cvnpart.photonid;
163  if(cvnpart.hadronid < 0)
164  cvnpart.hadronid = cvnpart.pionid + cvnpart.protonid;
165 
166  }
xml_schema::id
::xsd::cxx::tree::id< char, ncname > id
Definition: Database.h:165
caf::SRCVNParticleResult::neutronid
float neutronid
Likelihood neutron.
Definition: SRCVNParticleResult.h:22
caf::SRCVNParticleResult::pionid
float pionid
Likelihood pion.
Definition: SRCVNParticleResult.h:23
caf::SRCVNParticleResult::muonid
float muonid
Likelihood muon.
Definition: SRCVNParticleResult.h:19
caf::SRCVNParticleResult::photonid
float photonid
Likelihood photon.
Definition: SRCVNParticleResult.h:25
caf::SRCVNParticleResult::electronid
float electronid
Likelihood electron.
Definition: SRCVNParticleResult.h:20
caf::SRCVNParticleResult::pdgmax
unsigned int pdgmax
pdg of largest likelihood
Definition: SRCVNParticleResult.h:32
caf::SRCVNParticleResult::protonid
float protonid
Likelihood proton.
Definition: SRCVNParticleResult.h:21
caf::SRCVNParticleResult::emid
float emid
Likelihood electron or photon.
Definition: SRCVNParticleResult.h:28
caf::SRCVNParticleResult::maxval
float maxval
Maximum likelihood among net outputs.
Definition: SRCVNParticleResult.h:31
caf::SRCVNParticleResult::otherid
float otherid
Likelihood other.
Definition: SRCVNParticleResult.h:26
caf::SRCVNParticleResult::pizeroid
float pizeroid
Likelihood pizero.
Definition: SRCVNParticleResult.h:24
caf::SRCVNParticleResult::hadronid
float hadronid
Likelihood Proton or Pion.
Definition: SRCVNParticleResult.h:29
febshutoff_auto.val
val
Definition: febshutoff_auto.py:277
void caf::FillCVNPixelMaps ( const cvn::PixelMap &  map,
caf::SRPixelObjMap &  srcvnmap,
bool  useGeV 
)

Definition at line 617 of file FillPIDs.cxx.

References caf::SRPixelObjMap::cvnlabmap, caf::SRPixelObjMap::cvnmap, caf::SRPixelObjMap::cvnobjmap, cvn::PixelMap::fBound, cvn::PixelMap::fFracX, cvn::PixelMap::fFracY, cvn::Boundary::FirstCell(), caf::SRPixelObjMap::firstcellx, caf::SRPixelObjMap::firstcelly, caf::SRPixelObjMap::firstplane, cvn::Boundary::FirstPlane(), cvn::PixelMap::fNCell, cvn::PixelMap::fNPlane, caf::SRPixelObjMap::hitfracx, caf::SRPixelObjMap::hitfracy, cvn::Boundary::LastCell(), caf::SRPixelObjMap::lastcellx, caf::SRPixelObjMap::lastcelly, caf::SRPixelObjMap::lastplane, cvn::Boundary::LastPlane(), caf::SRPixelObjMap::ncells, caf::SRPixelObjMap::nchan, caf::SRPixelObjMap::nplanes, cvn::PixelMap::PixelMapToObjVector(), and cvn::PixelMap::PixelMapToVector().

Referenced by caf::CAFMaker::produce().

619  {
620  srcvnmap.nchan = 2; //two views to pixel map
621  srcvnmap.nplanes = map.fNPlane/2;
622  srcvnmap.ncells = map.fNCell;
623  srcvnmap.hitfracx = map.fFracX;
624  srcvnmap.hitfracy = map.fFracY;
625  srcvnmap.firstplane = map.fBound.FirstPlane();
626  srcvnmap.lastplane = map.fBound.LastPlane();
627  srcvnmap.firstcellx = map.fBound.FirstCell(0);
628  srcvnmap.lastcellx = map.fBound.LastCell(0);
629  srcvnmap.firstcelly = map.fBound.FirstCell(1);
630  srcvnmap.lastcelly = map.fBound.LastCell(1);
631 
632  std::vector<unsigned char> tval = map.PixelMapToVector(useGeV);
633  std::copy(tval.begin(), tval.end(), srcvnmap.cvnmap);
634  std::vector<unsigned char> tlab = map.PixelMapToObjVector(true);
635  std::copy(tlab.begin(), tlab.end(), srcvnmap.cvnlabmap);
636  std::vector<unsigned char> tobj = map.PixelMapToObjVector(false);
637  std::copy(tobj.begin(), tobj.end(), srcvnmap.cvnobjmap);
638  }
caf::SRPixelObjMap::cvnlabmap
unsigned char cvnlabmap[16000]
store 1D array with hit labels
Definition: SRPixelObjMap.h:33
cvn::Boundary::FirstCell
int FirstCell(const unsigned int &view) const
Definition: Boundary.h:44
caf::SRPixelObjMap::nplanes
unsigned int nplanes
number of planes in pixel map
Definition: SRPixelObjMap.h:17
cvn::PixelMap::PixelMapToObjVector
std::vector< unsigned char > PixelMapToObjVector(bool useLabels) const
Definition: PixelMap.cxx:417
caf::SRPixelObjMap::ncells
unsigned int ncells
number of cells in pixel map
Definition: SRPixelObjMap.h:21
cvn::PixelMap::fNPlane
unsigned int fNPlane
Number of planes, length of pixel map.
Definition: PixelMap.h:117
cvn::Boundary::LastPlane
int LastPlane() const
Definition: Boundary.h:43
cvn::Boundary::LastCell
int LastCell(const unsigned int &view) const
Definition: Boundary.h:45
caf::SRPixelObjMap::lastplane
unsigned int lastplane
global position of last plane in map
Definition: SRPixelObjMap.h:24
cvn::PixelMap::PixelMapToVector
std::vector< unsigned char > PixelMapToVector(bool useGeV) const
Definition: PixelMap.cxx:449
caf::SRPixelObjMap::cvnobjmap
unsigned char cvnobjmap[16000]
store 1D array with hit object id
Definition: SRPixelObjMap.h:34
cvn::PixelMap::fFracY
float fFracY
fraction of y view hits contained in map
Definition: PixelMap.h:120
cvn::PixelMap::fNCell
unsigned int fNCell
Number of cells, width of pixel map.
Definition: PixelMap.h:118
caf::SRPixelObjMap::cvnmap
unsigned char cvnmap[16000]
store 1D array the size of the pixel map, currently 80 cells x 100 planes x 2 views. Although sparse, store this way for easier conversion and compression in hdf5. First element is channel 0, plane 0, cell 0, then channel 0, plane 0, cell 1, etc. Channel 0 is x view, 1 is y view
Definition: SRPixelObjMap.h:32
cvn::PixelMap::fFracX
float fFracX
fraction of x view hits contained in map
Definition: PixelMap.h:119
caf::SRPixelObjMap::nchan
unsigned int nchan
number of channels in pixel map
Definition: SRPixelObjMap.h:22
caf::SRPixelObjMap::firstcellx
unsigned int firstcellx
global position of first x cell in map
Definition: SRPixelObjMap.h:25
caf::SRPixelObjMap::lastcelly
unsigned int lastcelly
global position of last y cell in map
Definition: SRPixelObjMap.h:28
caf::SRPixelObjMap::lastcellx
unsigned int lastcellx
global position of last x cell in map
Definition: SRPixelObjMap.h:26
cvn::PixelMap::fBound
Boundary fBound
Definition: PixelMap.h:153
caf::SRPixelObjMap::firstcelly
unsigned int firstcelly
global position of first y cell in map
Definition: SRPixelObjMap.h:27
caf::SRPixelObjMap::hitfracy
float hitfracy
fraction of y hits in slice in y view pixel map
Definition: SRPixelObjMap.h:30
cvn::Boundary::FirstPlane
int FirstPlane() const
Definition: Boundary.h:42
caf::SRPixelObjMap::hitfracx
float hitfracx
faction of x hits in slice in x view pixel map
Definition: SRPixelObjMap.h:29
caf::SRPixelObjMap::firstplane
unsigned int firstplane
global position of first plane in map
Definition: SRPixelObjMap.h:23
void caf::FillCVNPixelMaps ( const cvn::PixelMap &  map,
caf::SRPixelMap &  srcvnmap,
bool  useGeV 
)

Definition at line 641 of file FillPIDs.cxx.

References caf::SRPixelMap::cvnmap, cvn::PixelMap::fBound, cvn::PixelMap::fFracX, cvn::PixelMap::fFracY, cvn::Boundary::FirstCell(), caf::SRPixelMap::firstcellx, caf::SRPixelMap::firstcelly, caf::SRPixelMap::firstplane, cvn::Boundary::FirstPlane(), cvn::PixelMap::fNCell, cvn::PixelMap::fNPlane, caf::SRPixelMap::hitfracx, caf::SRPixelMap::hitfracy, cvn::Boundary::LastCell(), caf::SRPixelMap::lastcellx, caf::SRPixelMap::lastcelly, caf::SRPixelMap::lastplane, cvn::Boundary::LastPlane(), caf::SRPixelMap::ncells, caf::SRPixelMap::nchan, caf::SRPixelMap::nplanes, and cvn::PixelMap::PixelMapToVector().

643  {
644  srcvnmap.nchan = 2; //two views to pixel map
645  srcvnmap.nplanes = map.fNPlane/2;
646  srcvnmap.ncells = map.fNCell;
647  srcvnmap.hitfracx = map.fFracX;
648  srcvnmap.hitfracy = map.fFracY;
649  srcvnmap.firstplane = map.fBound.FirstPlane();
650  srcvnmap.lastplane = map.fBound.LastPlane();
651  srcvnmap.firstcellx = map.fBound.FirstCell(0);
652  srcvnmap.lastcellx = map.fBound.LastCell(0);
653  srcvnmap.firstcelly = map.fBound.FirstCell(1);
654  srcvnmap.lastcelly = map.fBound.LastCell(1);
655 
656  std::vector<unsigned char> tval = map.PixelMapToVector(useGeV);
657  std::copy(tval.begin(), tval.end(), srcvnmap.cvnmap);
658  }
cvn::Boundary::FirstCell
int FirstCell(const unsigned int &view) const
Definition: Boundary.h:44
caf::SRPixelMap::firstcellx
unsigned int firstcellx
global position of first x cell in map
Definition: SRPixelMap.h:26
caf::SRPixelMap::lastcelly
unsigned int lastcelly
global position of last y cell in map
Definition: SRPixelMap.h:29
caf::SRPixelMap::hitfracy
float hitfracy
fraction of y hits in slice in y view pixel map
Definition: SRPixelMap.h:31
cvn::PixelMap::fNPlane
unsigned int fNPlane
Number of planes, length of pixel map.
Definition: PixelMap.h:117
caf::SRPixelMap::ncells
unsigned int ncells
number of cells in pixel map
Definition: SRPixelMap.h:22
cvn::Boundary::LastPlane
int LastPlane() const
Definition: Boundary.h:43
cvn::Boundary::LastCell
int LastCell(const unsigned int &view) const
Definition: Boundary.h:45
caf::SRPixelMap::nplanes
unsigned int nplanes
number of planes in pixel map
Definition: SRPixelMap.h:18
cvn::PixelMap::PixelMapToVector
std::vector< unsigned char > PixelMapToVector(bool useGeV) const
Definition: PixelMap.cxx:449
caf::SRPixelMap::firstplane
unsigned int firstplane
global position of first plane in map
Definition: SRPixelMap.h:24
caf::SRPixelMap::lastcellx
unsigned int lastcellx
global position of last x cell in map
Definition: SRPixelMap.h:27
cvn::PixelMap::fFracY
float fFracY
fraction of y view hits contained in map
Definition: PixelMap.h:120
caf::SRPixelMap::nchan
unsigned int nchan
number of channels in pixel map
Definition: SRPixelMap.h:23
cvn::PixelMap::fNCell
unsigned int fNCell
Number of cells, width of pixel map.
Definition: PixelMap.h:118
caf::SRPixelMap::hitfracx
float hitfracx
faction of x hits in slice in x view pixel map
Definition: SRPixelMap.h:30
caf::SRPixelMap::lastplane
unsigned int lastplane
global position of last plane in map
Definition: SRPixelMap.h:25
caf::SRPixelMap::firstcelly
unsigned int firstcelly
global position of first y cell in map
Definition: SRPixelMap.h:28
cvn::PixelMap::fFracX
float fFracX
fraction of x view hits contained in map
Definition: PixelMap.h:119
caf::SRPixelMap::cvnmap
unsigned char cvnmap[16000]
store 1D array the size of the pixel map, currently 80 cells x 100 planes x 2 views. Although sparse, store this way for easier conversion and compression in hdf5. First element is channel 0, plane 0, cell 0, then channel 0, plane 0, cell 1, etc. Channel 0 is x view, 1 is y view
Definition: SRPixelMap.h:33
cvn::PixelMap::fBound
Boundary fBound
Definition: PixelMap.h:153
cvn::Boundary::FirstPlane
int FirstPlane() const
Definition: Boundary.h:42
void caf::FillCVNProngTrainingData ( const cvn::ProngTrainingData &  cvnpdata,
caf::SRProngTrainingData &  srpdata 
)

Definition at line 677 of file FillPIDs.cxx.

References cvn::ProngTrainingData::fIsPrimary, cvn::ProngTrainingData::fLabel3D, cvn::ProngTrainingData::fLabelX, cvn::ProngTrainingData::fLabelY, cvn::ProngTrainingData::fNCellX, cvn::ProngTrainingData::fNCellY, cvn::ProngTrainingData::fPAngle, cvn::ProngTrainingData::fPGap, cvn::ProngTrainingData::fPLength, cvn::ProngTrainingData::fPurity3D, cvn::ProngTrainingData::fPurityX, cvn::ProngTrainingData::fPurityY, cvn::ProngTrainingData::fRecE, cvn::ProngTrainingData::fVertX, cvn::ProngTrainingData::fVertY, cvn::ProngTrainingData::fVertZ, caf::SRProngTrainingData::isprimary, caf::SRProngTrainingData::label3d, caf::SRProngTrainingData::labelx, caf::SRProngTrainingData::labely, caf::SRProngTrainingData::ncellx, caf::SRProngTrainingData::ncelly, caf::SRProngTrainingData::pangle, caf::SRProngTrainingData::pgap, caf::SRProngTrainingData::plength, caf::SRProngTrainingData::purity3d, caf::SRProngTrainingData::purityx, caf::SRProngTrainingData::purityy, caf::SRProngTrainingData::rece, caf::SRProngTrainingData::vertx, caf::SRProngTrainingData::verty, and caf::SRProngTrainingData::vertz.

Referenced by caf::CAFMaker::produce().

679  {
680  srpdata.label3d = cvnpdata.fLabel3D;
681  srpdata.labelx = cvnpdata.fLabelX;
682  srpdata.labely = cvnpdata.fLabelY;
683  srpdata.isprimary = cvnpdata.fIsPrimary;
684  srpdata.purity3d = cvnpdata.fPurity3D;
685  srpdata.purityx = cvnpdata.fPurityX;
686  srpdata.purityy = cvnpdata.fPurityY;
687  srpdata.rece = cvnpdata.fRecE;
688  srpdata.ncellx = cvnpdata.fNCellX;
689  srpdata.ncelly = cvnpdata.fNCellY;
690  srpdata.vertx = cvnpdata.fVertX;
691  srpdata.verty = cvnpdata.fVertY;
692  srpdata.vertz = cvnpdata.fVertZ;
693  srpdata.plength = cvnpdata.fPLength;
694  srpdata.pangle = cvnpdata.fPAngle;
695  srpdata.pgap = cvnpdata.fPGap;
696  }
caf::SRProngTrainingData::labelx
int labelx
Class of the X view prong.
Definition: SRProngTrainingData.h:21
cvn::ProngTrainingData::fLabelY
ProngType fLabelY
Class of the Y view prong.
Definition: ProngTrainingData.h:35
cvn::ProngTrainingData::fLabelX
ProngType fLabelX
Class of the X view prong.
Definition: ProngTrainingData.h:34
caf::SRProngTrainingData::labely
int labely
Class of the Y view prong.
Definition: SRProngTrainingData.h:22
cvn::ProngTrainingData::fLabel3D
ProngType fLabel3D
Class of the 3D prong.
Definition: ProngTrainingData.h:33
caf::SRProngTrainingData::label3d
int label3d
Class of the 3D prong.
Definition: SRProngTrainingData.h:20
void caf::FillCVNResultVars ( const cvn::Result &  result,
caf::SRCVNResult &  srcvn,
bool  NoCosmics 
)

Definition at line 595 of file FillPIDs.cxx.

References caf::SRCVNResult::argmax, caf::SRCVNResult::cosmicid, distance(), cvn::Result::fOutput, caf::SRCVNResult::maxval, caf::SRCVNResult::ncid, caf::SRCVNResult::noutput, caf::SRCVNResult::nueid, caf::SRCVNResult::numuid, caf::SRCVNResult::nutauid, and caf::SRCVNResult::output.

Referenced by caf::CAFMaker::produce().

598  {
599  srcvn.output = result.fOutput;
600 
601  std::vector<float>::const_iterator maxIt =
602  std::max_element(result.fOutput.cbegin(), result.fOutput.cend());
603 
604  srcvn.maxval = *maxIt;
605  srcvn.argmax = std::distance(result.fOutput.cbegin(), maxIt);
606  srcvn.noutput = result.fOutput.size();
607 
608  srcvn.numuid = result.fOutput[0];
609  srcvn.nueid = result.fOutput[1];
610  srcvn.nutauid = -5;
611  srcvn.ncid = result.fOutput[2];
612  srcvn.cosmicid = result.fOutput[3];
613  if (NoCosmics) srcvn.cosmicid = -5;
614  }
caf::SRCVNResult::ncid
float ncid
Likelihood Neutral Current.
Definition: SRCVNResult.h:23
caf::SRCVNResult::maxval
float maxval
Maximum value among net outputs.
Definition: SRCVNResult.h:25
distance
unsigned distance(const T &t1, const T &t2)
Definition: Multiplet_module.cc:33
caf::SRCVNResult::nutauid
float nutauid
Likelihood Charge Current NuTau.
Definition: SRCVNResult.h:22
cvn::Result::fOutput
std::vector< float > fOutput
Vector of outputs from neural net.
Definition: Result.h:30
caf::SRCVNResult::nueid
float nueid
Likelihood Charge Current NuE.
Definition: SRCVNResult.h:21
caf::SRCVNResult::noutput
unsigned int noutput
Number of entries in output vector.
Definition: SRCVNResult.h:27
caf::SRCVNResult::cosmicid
float cosmicid
Likelihood Cosmic Interaction.
Definition: SRCVNResult.h:24
caf::SRCVNResult::output
std::vector< float > output
List of net output nodes.
Definition: SRCVNResult.h:19
caf::SRCVNResult::numuid
float numuid
Likelihood Charge Current NuMu.
Definition: SRCVNResult.h:20
caf::SRCVNResult::argmax
unsigned int argmax
Index of maximum value.
Definition: SRCVNResult.h:26
void caf::FillCVNTrainingData ( const cvn::TrainingData &  cvntdata,
caf::SRTrainingData &  srtdata 
)

Definition at line 661 of file FillPIDs.cxx.

References cvn::TrainingData::fFSlong, cvn::TrainingData::fFSp, cvn::TrainingData::fFSt, caf::SRTrainingData::finalstate, caf::SRTrainingData::finalstateprong, cvn::TrainingData::fInt, cvn::TrainingData::fLepEnergy, cvn::TrainingData::fNuEnergy, cvn::TrainingData::fPPt, cvn::TrainingData::fVtxX, cvn::TrainingData::fVtxY, cvn::TrainingData::fVtxZ, caf::SRTrainingData::interaction, caf::SRTrainingData::lepenergy, caf::SRTrainingData::nuenergy, caf::SRTrainingData::parent, caf::SRTrainingData::particles, caf::SRTrainingData::vtxx, caf::SRTrainingData::vtxy, and caf::SRTrainingData::vtxz.

Referenced by caf::CAFMaker::produce().

663  {
664  srtdata.parent = cvntdata.fPPt;
665  srtdata.interaction = cvntdata.fInt;
666  srtdata.finalstate = cvntdata.fFSt;
667  srtdata.finalstateprong = cvntdata.fFSp;
668  srtdata.particles = cvntdata.fFSlong;
669  srtdata.nuenergy = cvntdata.fNuEnergy;
670  srtdata.lepenergy = cvntdata.fLepEnergy;
671  srtdata.vtxx = cvntdata.fVtxX;
672  srtdata.vtxy = cvntdata.fVtxY;
673  srtdata.vtxz = cvntdata.fVtxZ;
674  }
caf::SRTrainingData::interaction
unsigned int interaction
Class of the event.
Definition: SRTrainingData.h:21
cvn::TrainingData::fInt
InteractionType fInt
Class of the event.
Definition: TrainingData.h:60
cvn::TrainingData::fFSlong
ParticlesType fFSlong
Class of the event by primary n prongs.
Definition: TrainingData.h:63
caf::SRTrainingData::parent
unsigned int parent
Class of the event by parent particle.
Definition: SRTrainingData.h:20
caf::SRTrainingData::particles
unsigned int particles
Class of the event by primary n prongs.
Definition: SRTrainingData.h:24
cvn::TrainingData::fVtxY
double fVtxY
True vertex of Y position.
Definition: TrainingData.h:67
caf::SRTrainingData::vtxx
double vtxx
True vertex of X position.
Definition: SRTrainingData.h:27
cvn::TrainingData::fLepEnergy
float fLepEnergy
True energy of outgoing lepton.
Definition: TrainingData.h:65
cvn::TrainingData::fNuEnergy
float fNuEnergy
True energy of neutrino event.
Definition: TrainingData.h:64
caf::SRTrainingData::vtxz
double vtxz
True vertex of Z position.
Definition: SRTrainingData.h:29
caf::SRTrainingData::vtxy
double vtxy
True vertex of Y position.
Definition: SRTrainingData.h:28
cvn::TrainingData::fVtxX
double fVtxX
True vertex of X position.
Definition: TrainingData.h:66
cvn::TrainingData::fFSp
FinalStateProngsType fFSp
Class of the event by primary n prongs.
Definition: TrainingData.h:62
caf::SRTrainingData::nuenergy
float nuenergy
True energy of neutrino event.
Definition: SRTrainingData.h:25
cvn::TrainingData::fVtxZ
double fVtxZ
True vertex of Z position.
Definition: TrainingData.h:68
caf::SRTrainingData::finalstate
unsigned int finalstate
Class of the event by final state.
Definition: SRTrainingData.h:22
cvn::TrainingData::fPPt
ParentParticleType fPPt
Class of the event by parent particle.
Definition: TrainingData.h:59
cvn::TrainingData::fFSt
FinalStateType fFSt
Class of the event by final state.
Definition: TrainingData.h:61
caf::SRTrainingData::finalstateprong
unsigned int finalstateprong
Definition: SRTrainingData.h:23
caf::SRTrainingData::lepenergy
float lepenergy
True energy of outgoing lepton.
Definition: SRTrainingData.h:26
void caf::FillDiFShowerVars ( const rb::Shower &  shw,
const std::vector< rb::Cluster > &  sliceList,
std::vector< cheat::TrackIDE > &  allTracks,
std::vector< cheat::TrackIDE > &  sliceTracks,
std::vector< cheat::TrackIDE > &  allTracksBirks,
std::vector< cheat::TrackIDE > &  sliceTracksBirks,
int  sliceIdx,
SRMRProperties &  srMR 
)

Definition at line 504 of file FillReco.cxx.

References rb::Cluster::AllCells(), caf::SRShower::calE, rb::Cluster::CalorimetricEnergy(), caf::SRShower::dir, rb::Prong::Dir(), rb::Cluster::ExtentPlane(), FillParticleTruth(), if(), caf::SRMRProperties::isShwDiF, IsTrueDiF(), kX, geo::kX, kXorY, geo::kXorY, kY, geo::kY, caf::SRShower::len, caf::SRMRProperties::lid, caf::SRShower::maxplanecont, caf::SRShower::maxplanegap, rb::Cluster::MostContiguousPlanes(), rb::Cluster::MostMissingPlanes(), rb::Cluster::NCell(), caf::SRShower::nhit, caf::SRShower::nhitx, caf::SRShower::nhity, caf::SRShower::nplane, caf::SRShower::nplanex, caf::SRShower::nplaney, rb::Cluster::NXCell(), rb::Cluster::NYCell(), caf::SRShower::start, rb::Prong::Start(), caf::SRShower::stop, rb::Shower::Stop(), rb::Shower::TotalLength(), caf::SRShower::truth, caf::SRShower::truthXView, caf::SRShower::truthYView, caf::SRShower::view, rb::Cluster::XCells(), and rb::Cluster::YCells().

Referenced by caf::CAFMaker::produce().

511  {
512  srMR.isShwDiF = IsTrueDiF(shw);
513 
514  srMR.lid.nhit = shw.NCell();
515  srMR.lid.nhitx = shw.NXCell();
516  srMR.lid.nhity = shw.NYCell();
517  srMR.lid.nplane = shw.ExtentPlane();
518  srMR.lid.nplanex = shw.ExtentPlane(geo::kX);
519  srMR.lid.nplaney = shw.ExtentPlane(geo::kY);
520  srMR.lid.maxplanecont = shw.MostContiguousPlanes(geo::kXorY);
521  srMR.lid.maxplanegap = shw.MostMissingPlanes(geo::kXorY);
522  srMR.lid.len = shw.TotalLength();
523  srMR.lid.calE = shw.CalorimetricEnergy();
524  srMR.lid.start = shw.Start();
525  srMR.lid.dir = shw.Dir();
526  srMR.lid.stop = shw.Stop();
527 
528  srMR.lid.truth = FillParticleTruth(sliceList,
529  shw.AllCells(),
530  allTracks,
531  sliceTracks,
532  allTracksBirks,
533  sliceTracksBirks,
534  sliceIdx);
535  srMR.lid.truthXView = FillParticleTruth(sliceList,
536  shw.XCells(),
537  allTracks,
538  sliceTracks,
539  allTracksBirks,
540  sliceTracksBirks,
541  sliceIdx);
542  srMR.lid.truthYView = FillParticleTruth(sliceList,
543  shw.YCells(),
544  allTracks,
545  sliceTracks,
546  allTracksBirks,
547  sliceTracksBirks,
548  sliceIdx);
549  if (shw.NXCell() == 0) srMR.lid.view = kY;
550  else if (shw.NYCell() == 0) srMR.lid.view = kX;
551  else srMR.lid.view = kXorY;
552  }
rb::Cluster::NCell
unsigned int NCell(geo::View_t view) const
Number of cells in view view.
Definition: Cluster.cxx:134
rb::Cluster::XCells
const art::PtrVector< rb::CellHit > & XCells() const
Get all cells from the x-view.
Definition: Cluster.h:124
geo::kXorY
X or Y views.
Definition: PlaneGeo.h:30
kXorY
Definition: geo.h:1
geo::kX
Vertical planes which measure X.
Definition: PlaneGeo.h:28
rb::Prong::Start
virtual TVector3 Start() const
Definition: Prong.h:73
geo::kY
Horizontal planes which measure Y.
Definition: PlaneGeo.h:29
kY
const Var kY([](const caf::SRProxy *sr){float tmp=0.f;if(sr->mc.nu.empty()) return tmp;tmp=sr->mc.nu[0].y;return tmp;})
rb::Cluster::AllCells
art::PtrVector< rb::CellHit > AllCells() const
Get all cells from both views.
Definition: Cluster.cxx:180
rb::Shower::TotalLength
virtual double TotalLength() const
Length (cm) of a shower.
Definition: Shower.cxx:43
caf::FillParticleTruth
SRParticleTruth FillParticleTruth(const std::vector< rb::Cluster > &sliceList, const art::PtrVector< rb::CellHit > &hits, std::vector< cheat::TrackIDE > &allTracks, std::vector< cheat::TrackIDE > &sliceTracks, std::vector< cheat::TrackIDE > &allTracksBirks, std::vector< cheat::TrackIDE > &sliceTracksBirks, int sliceIdx)
Function to fill the particle truth for a set of hits.
Definition: FillTruth.cxx:28
rb::Cluster::MostMissingPlanes
int MostMissingPlanes(geo::View_t view) const
Longest run of adjacent planes with no hits.
Definition: Cluster.cxx:668
rb::Cluster::CalorimetricEnergy
double CalorimetricEnergy(EEnergyCalcScheme escheme=kRecomputeEnergy) const
Simple estimate of neutrino energy.
Definition: Cluster.cxx:439
if
if(dump)
Definition: checkMeasured_Target.C:110
rb::Cluster::YCells
const art::PtrVector< rb::CellHit > & YCells() const
Get all cells from the x-view.
Definition: Cluster.h:126
rb::Prong::Dir
virtual TVector3 Dir() const
Unit vector describing prong direction.
Definition: Prong.h:77
rb::Shower::Stop
TVector3 Stop() const
Endpoint of the shower.
Definition: Shower.cxx:55
rb::Cluster::NYCell
unsigned int NYCell() const
Number of cells in the y-view.
Definition: Cluster.h:108
caf::IsTrueDiF
bool IsTrueDiF(const rb::Cluster &slice)
Helper for FillDiFVars and FillDiFShowerVars.
Definition: FillReco.cxx:67
rb::Cluster::NXCell
unsigned int NXCell() const
Number of cells in the x-view.
Definition: Cluster.h:106
rb::Cluster::MostContiguousPlanes
int MostContiguousPlanes(geo::View_t view) const
Longest run of adjacent planes with hits.
Definition: Cluster.cxx:635
kX
Definition: geo.h:1
rb::Cluster::ExtentPlane
unsigned int ExtentPlane(geo::View_t view=geo::kXorY) const
Definition: Cluster.h:250
void caf::FillDiFVars ( const rb::Cluster &  slice,
const std::vector< rb::Cluster > &  sliceList,
std::vector< cheat::TrackIDE > &  allTracks,
std::vector< cheat::TrackIDE > &  sliceTracks,
std::vector< cheat::TrackIDE > &  allTracksBirks,
std::vector< cheat::TrackIDE > &  sliceTracksBirks,
int  sliceIdx,
SRMRProperties &  srMR 
)

Definition at line 97 of file FillReco.cxx.

References rb::Cluster::AllCells(), caf::SRMRProperties::calE, rb::Cluster::CalorimetricEnergy(), rb::Cluster::ExtentPlane(), FillParticleTruth(), if(), caf::SRMRProperties::isDiF, IsTrueDiF(), kX, geo::kX, kXorY, geo::kXorY, kY, geo::kY, caf::SRMRProperties::maxplane, rb::Cluster::MaxPlane(), caf::SRMRProperties::maxplanecont, caf::SRMRProperties::maxplanegap, caf::SRMRProperties::minplane, rb::Cluster::MinPlane(), rb::Cluster::MostContiguousPlanes(), rb::Cluster::MostMissingPlanes(), rb::Cluster::NCell(), caf::SRMRProperties::nhit, caf::SRMRProperties::nhitx, caf::SRMRProperties::nhity, caf::SRMRProperties::nplane, caf::SRMRProperties::nplanex, caf::SRMRProperties::nplaney, rb::Cluster::NXCell(), rb::Cluster::NYCell(), caf::SRMRProperties::truth, caf::SRMRProperties::truthXView, caf::SRMRProperties::truthYView, caf::SRMRProperties::view, rb::Cluster::XCells(), and rb::Cluster::YCells().

Referenced by caf::CAFMaker::produce().

105  {
106  srMR.isDiF = IsTrueDiF(slice);
107 
108  srMR.nhit = slice.NCell();
109  srMR.nhitx = slice.NXCell();
110  srMR.nhity = slice.NYCell();
111  srMR.nplane = slice.ExtentPlane();
112  srMR.nplanex = slice.ExtentPlane(geo::kX);
113  srMR.nplaney = slice.ExtentPlane(geo::kY);
114  srMR.maxplanecont = slice.MostContiguousPlanes(geo::kXorY);
115  srMR.maxplanegap = slice.MostMissingPlanes(geo::kXorY);
116  srMR.minplane = slice.MinPlane(geo::kXorY);
117  srMR.maxplane = slice.MaxPlane(geo::kXorY);
118  srMR.calE = slice.CalorimetricEnergy();
119 
120  srMR.truth = FillParticleTruth(sliceList,
121  slice.AllCells(),
122  allTracks,
123  sliceTracks,
124  allTracksBirks,
125  sliceTracksBirks,
126  sliceIdx);
127  srMR.truthXView = FillParticleTruth(sliceList,
128  slice.XCells(),
129  allTracks,
130  sliceTracks,
131  allTracksBirks,
132  sliceTracksBirks,
133  sliceIdx);
134  srMR.truthYView = FillParticleTruth(sliceList,
135  slice.YCells(),
136  allTracks,
137  sliceTracks,
138  allTracksBirks,
139  sliceTracksBirks,
140  sliceIdx);
141  if (slice.NXCell() == 0) srMR.view = kY;
142  else if (slice.NYCell() == 0) srMR.view = kX;
143  else srMR.view = kXorY;
144  }
rb::Cluster::NCell
unsigned int NCell(geo::View_t view) const
Number of cells in view view.
Definition: Cluster.cxx:134
rb::Cluster::XCells
const art::PtrVector< rb::CellHit > & XCells() const
Get all cells from the x-view.
Definition: Cluster.h:124
geo::kXorY
X or Y views.
Definition: PlaneGeo.h:30
kXorY
Definition: geo.h:1
geo::kX
Vertical planes which measure X.
Definition: PlaneGeo.h:28
geo::kY
Horizontal planes which measure Y.
Definition: PlaneGeo.h:29
kY
const Var kY([](const caf::SRProxy *sr){float tmp=0.f;if(sr->mc.nu.empty()) return tmp;tmp=sr->mc.nu[0].y;return tmp;})
rb::Cluster::AllCells
art::PtrVector< rb::CellHit > AllCells() const
Get all cells from both views.
Definition: Cluster.cxx:180
caf::FillParticleTruth
SRParticleTruth FillParticleTruth(const std::vector< rb::Cluster > &sliceList, const art::PtrVector< rb::CellHit > &hits, std::vector< cheat::TrackIDE > &allTracks, std::vector< cheat::TrackIDE > &sliceTracks, std::vector< cheat::TrackIDE > &allTracksBirks, std::vector< cheat::TrackIDE > &sliceTracksBirks, int sliceIdx)
Function to fill the particle truth for a set of hits.
Definition: FillTruth.cxx:28
rb::Cluster::MostMissingPlanes
int MostMissingPlanes(geo::View_t view) const
Longest run of adjacent planes with no hits.
Definition: Cluster.cxx:668
rb::Cluster::CalorimetricEnergy
double CalorimetricEnergy(EEnergyCalcScheme escheme=kRecomputeEnergy) const
Simple estimate of neutrino energy.
Definition: Cluster.cxx:439
if
if(dump)
Definition: checkMeasured_Target.C:110
rb::Cluster::YCells
const art::PtrVector< rb::CellHit > & YCells() const
Get all cells from the x-view.
Definition: Cluster.h:126
rb::Cluster::NYCell
unsigned int NYCell() const
Number of cells in the y-view.
Definition: Cluster.h:108
rb::Cluster::MinPlane
unsigned int MinPlane(geo::View_t view=geo::kXorY) const
Definition: Cluster.cxx:462
caf::IsTrueDiF
bool IsTrueDiF(const rb::Cluster &slice)
Helper for FillDiFVars and FillDiFShowerVars.
Definition: FillReco.cxx:67
rb::Cluster::NXCell
unsigned int NXCell() const
Number of cells in the x-view.
Definition: Cluster.h:106
rb::Cluster::MostContiguousPlanes
int MostContiguousPlanes(geo::View_t view) const
Longest run of adjacent planes with hits.
Definition: Cluster.cxx:635
kX
Definition: geo.h:1
rb::Cluster::ExtentPlane
unsigned int ExtentPlane(geo::View_t view=geo::kXorY) const
Definition: Cluster.h:250
rb::Cluster::MaxPlane
unsigned int MaxPlane(geo::View_t view=geo::kXorY) const
Definition: Cluster.cxx:508
void caf::FillHadClustVars ( const rb::Cluster &  hadclust,
caf::SRHadClust &  srhadclust 
)

Definition at line 101 of file FillEnergies.cxx.

References caf::SRHadClust::boxmax, caf::SRHadClust::boxmin, caf::SRHadClust::calE, rb::Cluster::CalorimetricEnergy(), rb::CellHit::Cell(), rb::Cluster::Cell(), caf::SRHadClust::firstcell, caf::SRHadClust::firstplane, geom(), art::Ptr< T >::get(), rb::RecoHit::IsCalibrated(), geo::kX, geo::kXorY, caf::SRHadClust::lastcell, caf::SRHadClust::lastplane, calib::Calibrator::MakeRecoHit(), rb::Cluster::MaxCell(), rb::Cluster::MaxPlane(), caf::SRHadClust::meanpos, std::min(), rb::Cluster::MinCell(), rb::Cluster::MinMaxMeanXYZ(), rb::Cluster::MinPlane(), rb::Cluster::MostContiguousPlanes(), caf::SRHadClust::ncalhit, rb::Cluster::NCell(), geo::PlaneGeo::Ncells(), caf::SRHadClust::ncellsfromedge, caf::SRHadClust::ncontplanes, caf::SRHadClust::nhit, caf::SRHadClust::nmiphit, rb::Cluster::NXCell(), rb::Cluster::NYCell(), rb::RecoHit::PECorr(), rb::CellHit::Plane(), geo::GeometryBase::Plane(), caf::SRVector3D::SetX(), caf::SRVector3D::SetY(), rb::CellHit::View(), rb::Cluster::W(), rb::Cluster::XCell(), and rb::Cluster::YCell().

Referenced by FillNumuEnergyVars().

103  {
104  static art::ServiceHandle<calib::Calibrator> cal;
105  static art::ServiceHandle<geo::Geometry> geom;
106 
107  const double wx = (hadclust.NXCell() > 0) ? hadclust.W(hadclust.XCell(0).get()) : 0;
108  const double wy = (hadclust.NYCell() > 0) ? hadclust.W(hadclust.YCell(0).get()) : 0;
109 
110  // Variables needed:
111  // Holds contiguous planes info:
112  unsigned int cellsFromEdge = -1;
113 
114  // Hit counters:
115  int nCalHit = 0;
116  int nMipHit = 0;
117 
118  unsigned int minCellsFromEdge = 99999999;
119 
120  // Loop over hadclust in question here:
121  for(unsigned int hitIdx = 0; hitIdx < hadclust.NCell(); ++hitIdx)
122  {
123  const art::Ptr<rb::CellHit>& chit = hadclust.Cell(hitIdx);
124  // Get calbrated reco hit:
125  const rb::RecoHit rhit(cal->MakeRecoHit(*chit,
126  chit->View() == geo::kX ? wx : wy));
127 
128  const int planeNum = chit->Plane();
129  cellsFromEdge = std::min((unsigned int)chit->Cell(), geom->Plane(planeNum)->Ncells() - 1 - chit->Cell());
130 
131  if(cellsFromEdge < minCellsFromEdge)
132  {
133  minCellsFromEdge = cellsFromEdge;
134  }
135 
136  // Make sure things have been calibrated:
137  if(!rhit.IsCalibrated())
138  {
139  // MF_LOG_DEBUG("CAFMaker") << "Not calibrated?! "
140  // << chit->Plane() << " " << chit->Cell();
141  continue;
142  }
143 
144  // Count up both the total calibrated hits
145  // and the MIP hits from 100-245 PECorr
146  // Count up total hits:
147  nCalHit++;
148  // Count up MIP hits:
149  if ( (rhit.PECorr()>100.0) && (rhit.PECorr()<245.0) ) nMipHit++;
150  } // end of hadclust loop
151 
152  srhadclust.nhit = hadclust.NCell();
153  srhadclust.ncalhit = nCalHit;
154  srhadclust.nmiphit = nMipHit;
155  srhadclust.ncontplanes = hadclust.MostContiguousPlanes(geo::kXorY);
156  if (!(hadclust.NCell()==0)){
157 
158  srhadclust.firstplane = hadclust.MinPlane();
159  srhadclust.lastplane = hadclust.MaxPlane();
160  srhadclust.firstcell = hadclust.MinCell(geo::kXorY);
161  srhadclust.lastcell = hadclust.MaxCell(geo::kXorY);
162  srhadclust.ncellsfromedge = minCellsFromEdge;
163 
164  srhadclust.calE = hadclust.CalorimetricEnergy();
165 
166  TVector3 boxmin, boxmax, meanpos;
167  hadclust.MinMaxMeanXYZ(boxmin, boxmax, meanpos);
168  srhadclust.boxmin = boxmin;
169  srhadclust.boxmax = boxmax;
170  srhadclust.meanpos = meanpos;
171 
172  // TBrowser can't draw DBL_MAX
173  if (hadclust.NXCell()==0){
174  srhadclust.boxmin.SetX(-5555);
175  srhadclust.boxmax.SetX(5555);
176  }
177  if (hadclust.NYCell()==0){
178  srhadclust.boxmin.SetY(-5555);
179  srhadclust.boxmax.SetY(5555);
180  }
181 
182  } // end of loop over variables to only fill for hadclust if it has hits
183  else{
184  srhadclust.firstplane = 0;
185  srhadclust.lastplane = 0;
186  srhadclust.firstcell = 0;
187  srhadclust.lastcell = 0;
188  srhadclust.ncellsfromedge = 0;
189 
190  srhadclust.calE = -5;
191 
192  srhadclust.boxmin = TVector3(-5555,-5555,-5555);
193  srhadclust.boxmax = TVector3(-5555,-5555,-5555);
194  srhadclust.