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

This module creates Common Analysis Files. More...

Classes

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::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::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::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::SRNumuEnergy >
 Proxy for caf::SRNumuEnergy. 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  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  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  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  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  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  SRNumuEnergy
 Numu energy estimator output. More...
 
class  SRNumuSandbox
 
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  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...
 

Typedefs

using SRProxy = caf::Proxy< caf::StandardRecord >
 
using SRSLidEnergyProxy = caf::Proxy< SRSLidEnergy >
 
using SRNueEnergyProxy = caf::Proxy< SRNueEnergy >
 
using SRBPFEnergyProxy = caf::Proxy< SRBPFEnergy >
 
using SRVector3DProxy = caf::Proxy< SRVector3D >
 
using SRHadClustProxy = caf::Proxy< SRHadClust >
 
using SRTrueNumuEnergyProxy = caf::Proxy< SRTrueNumuEnergy >
 
using SRNumuEnergyProxy = caf::Proxy< SRNumuEnergy >
 
using SREnergyBranchProxy = caf::Proxy< SREnergyBranch >
 
using SRHeaderProxy = caf::Proxy< SRHeader >
 
using SRLorentzVectorProxy = caf::Proxy< SRLorentzVector >
 
using SRTrueMichelEProxy = caf::Proxy< SRTrueMichelE >
 
using SRCosmicProxy = caf::Proxy< SRCosmic >
 
using SRBeamProxy = caf::Proxy< SRBeam >
 
using SRTrueParticleProxy = caf::Proxy< SRTrueParticle >
 
using SRGenieWeightsProxy = caf::Proxy< SRGenieWeights >
 
using SRFluxWeightsProxy = caf::Proxy< SRFluxWeights >
 
using SRMCReweightProxy = caf::Proxy< SRMCReweight >
 
using SRNeutrinoProxy = caf::Proxy< SRNeutrino >
 
using SRGlobalTruthProxy = caf::Proxy< SRGlobalTruth >
 
using SRTruthBranchProxy = caf::Proxy< SRTruthBranch >
 
using SRParticleTruthProxy = caf::Proxy< SRParticleTruth >
 
using SRSlcMEProxy = caf::Proxy< SRSlcME >
 
using SRTrkMEProxy = caf::Proxy< SRTrkME >
 
using SRMichelEProxy = caf::Proxy< SRMichelE >
 
using SRMRCCParentProxy = caf::Proxy< SRMRCCParent >
 
using SRParentBranchProxy = caf::Proxy< SRParentBranch >
 
using SRBpfIdProxy = caf::Proxy< SRBpfId >
 
using SRContainProxy = caf::Proxy< SRContain >
 
using SRCosRejProxy = caf::Proxy< SRCosRej >
 
using SRCVNResultProxy = caf::Proxy< SRCVNResult >
 
using SRLemProxy = caf::Proxy< SRLem >
 
using SRELidProxy = caf::Proxy< SRELid >
 
using SRNCCosRejProxy = caf::Proxy< SRNCCosRej >
 
using SRNCPi0BkgRejProxy = caf::Proxy< SRNCPi0BkgRej >
 
using SRNueCosRejProxy = caf::Proxy< SRNueCosRej >
 
using SRPreselProxy = caf::Proxy< SRPresel >
 
using SRRemidProxy = caf::Proxy< SRRemid >
 
using SRRvpProxy = caf::Proxy< SRRvp >
 
using SRSliceLIDProxy = caf::Proxy< SRSliceLID >
 
using SRVetoProxy = caf::Proxy< SRVeto >
 
using SRXnueProxy = caf::Proxy< SRXnue >
 
using SRIDBranchProxy = caf::Proxy< SRIDBranch >
 
using SRSliceProxy = caf::Proxy< SRSlice >
 
using SRCosmicCVNProxy = caf::Proxy< SRCosmicCVN >
 
using SRSpillProxy = caf::Proxy< SRSpill >
 
using SRCVNFeaturesProxy = caf::Proxy< SRCVNFeatures >
 
using SRPixelObjMapProxy = caf::Proxy< SRPixelObjMap >
 
using SRSliceMapProxy = caf::Proxy< SRSliceMap >
 
using SRTrainingDataProxy = caf::Proxy< SRTrainingData >
 
using SRTrainingBranchProxy = caf::Proxy< SRTrainingBranch >
 
using SRShowerProxy = caf::Proxy< SRShower >
 
using SRSLidProxy = caf::Proxy< SRSLid >
 
using SRShowerLIDProxy = caf::Proxy< SRShowerLID >
 
using SRMRPropertiesProxy = caf::Proxy< SRMRProperties >
 
using SRTrackProxy = caf::Proxy< SRTrack >
 
using SRTrackBaseProxy = caf::Proxy< SRTrackBase >
 
using SRKalmanTrackProxy = caf::Proxy< SRKalmanTrack >
 
using SRKalmanProxy = caf::Proxy< SRKalman >
 
using SRTrackBranchProxy = caf::Proxy< SRTrackBranch >
 
using SRPixelMapProxy = caf::Proxy< SRPixelMap >
 
using SRCVNParticleResultProxy = caf::Proxy< SRCVNParticleResult >
 
using SRProngTrainingDataProxy = caf::Proxy< SRProngTrainingData >
 
using SRProngProxy = caf::Proxy< SRProng >
 
using SRBpfTrackProxy = caf::Proxy< SRBpfTrack >
 
using SRBpfProxy = caf::Proxy< SRBpf >
 
using SRRegCVNResultProxy = caf::Proxy< SRRegCVNResult >
 
using SRFuzzyKProngProxy = caf::Proxy< SRFuzzyKProng >
 
using SRFuzzyKProxy = caf::Proxy< SRFuzzyK >
 
using SRElasticProxy = caf::Proxy< SRElastic >
 
using SRHoughVertexProxy = caf::Proxy< SRHoughVertex >
 
using SRVertexDTProxy = caf::Proxy< SRVertexDT >
 
using SRVertexBranchProxy = caf::Proxy< SRVertexBranch >
 
using StandardRecordProxy = caf::Proxy< StandardRecord >
 

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

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, caf::SRNumuEnergy &cafE)
 
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 FillJMEIDVars (const jmshower::EID &eid, caf::SRJMEid &sreid)
 
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 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 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 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)
 
double TrueNeutrinoDistance (novadaq::cnv::DetId det, const SRNeutrino &nu)
 
std::vector< SRGenieWeights > GenieReweightTable (const rwgt::GENIEReweightTable &table)
 
SRFluxWeights FluxReweights (const fxwgt::FluxWeights &flxwgts)
 
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)
 

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
 

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

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

Definition at line 86 of file FwdDeclare.h.

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

Definition at line 26 of file FwdDeclare.h.

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

Definition at line 170 of file FwdDeclare.h.

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

Definition at line 398 of file FwdDeclare.h.

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

Definition at line 392 of file FwdDeclare.h.

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

Definition at line 176 of file FwdDeclare.h.

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

Definition at line 272 of file FwdDeclare.h.

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

Definition at line 80 of file FwdDeclare.h.

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

Definition at line 182 of file FwdDeclare.h.

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

Definition at line 284 of file FwdDeclare.h.

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

Definition at line 374 of file FwdDeclare.h.

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

Definition at line 188 of file FwdDeclare.h.

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

Definition at line 422 of file FwdDeclare.h.

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

Definition at line 200 of file FwdDeclare.h.

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

Definition at line 56 of file FwdDeclare.h.

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

Definition at line 104 of file FwdDeclare.h.

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

Definition at line 410 of file FwdDeclare.h.

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

Definition at line 416 of file FwdDeclare.h.

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

Definition at line 98 of file FwdDeclare.h.

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

Definition at line 122 of file FwdDeclare.h.

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

Definition at line 38 of file FwdDeclare.h.

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

Definition at line 62 of file FwdDeclare.h.

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

Definition at line 428 of file FwdDeclare.h.

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

Definition at line 260 of file FwdDeclare.h.

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

Definition at line 356 of file FwdDeclare.h.

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

Definition at line 350 of file FwdDeclare.h.

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

Definition at line 194 of file FwdDeclare.h.

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

Definition at line 68 of file FwdDeclare.h.

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

Definition at line 110 of file FwdDeclare.h.

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

Definition at line 152 of file FwdDeclare.h.

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

Definition at line 158 of file FwdDeclare.h.

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

Definition at line 332 of file FwdDeclare.h.

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

Definition at line 206 of file FwdDeclare.h.

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

Definition at line 212 of file FwdDeclare.h.

using caf::SRNeutrinoProxy = typedef caf::Proxy<SRNeutrino>

Definition at line 116 of file FwdDeclare.h.

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

Definition at line 218 of file FwdDeclare.h.

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

Definition at line 20 of file FwdDeclare.h.

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

Definition at line 50 of file FwdDeclare.h.

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

Definition at line 164 of file FwdDeclare.h.

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

Definition at line 134 of file FwdDeclare.h.

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

Definition at line 368 of file FwdDeclare.h.

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

Definition at line 290 of file FwdDeclare.h.

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

Definition at line 224 of file FwdDeclare.h.

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

Definition at line 386 of file FwdDeclare.h.

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

Definition at line 380 of file FwdDeclare.h.

typedef caf::Proxy< caf::StandardRecord > caf::SRProxy

Definition at line 2 of file EpilogFwd.h.

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

Definition at line 404 of file FwdDeclare.h.

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

Definition at line 230 of file FwdDeclare.h.

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

Definition at line 236 of file FwdDeclare.h.

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

Definition at line 326 of file FwdDeclare.h.

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

Definition at line 314 of file FwdDeclare.h.

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

Definition at line 140 of file FwdDeclare.h.

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

Definition at line 242 of file FwdDeclare.h.

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

Definition at line 296 of file FwdDeclare.h.

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

Definition at line 266 of file FwdDeclare.h.

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

Definition at line 14 of file FwdDeclare.h.

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

Definition at line 320 of file FwdDeclare.h.

using caf::SRSpillProxy = typedef caf::Proxy<SRSpill>

Definition at line 278 of file FwdDeclare.h.

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

Definition at line 344 of file FwdDeclare.h.

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

Definition at line 362 of file FwdDeclare.h.

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

Definition at line 338 of file FwdDeclare.h.

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

Definition at line 308 of file FwdDeclare.h.

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

Definition at line 302 of file FwdDeclare.h.

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

Definition at line 146 of file FwdDeclare.h.

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

Definition at line 74 of file FwdDeclare.h.

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

Definition at line 44 of file FwdDeclare.h.

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

Definition at line 92 of file FwdDeclare.h.

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

Definition at line 128 of file FwdDeclare.h.

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

Definition at line 32 of file FwdDeclare.h.

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

Definition at line 440 of file FwdDeclare.h.

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

Definition at line 434 of file FwdDeclare.h.

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

Definition at line 248 of file FwdDeclare.h.

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

Definition at line 254 of file FwdDeclare.h.

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

Definition at line 446 of file FwdDeclare.h.

Enumeration Type Documentation

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 532 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, stan::math::fabs(), 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_flatrecord::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().

538  {
539  static art::ServiceHandle<cheat::BackTracker> bt;
540  static art::ServiceHandle<geo::Geometry> geom;
541 
542  // Get information about true cosmic
543  vec->push_back(SRCosmic());
544  SRCosmic& srTruth = vec->back();
545 
546  const art::Ptr<simb::MCTruth> &tru = effPur->neutrinoInt;
547 
548  // BUT GetParticle(int i) requires looping over list of particles in
549  // single MCTruth object. So we have made an MCTruth object for each
550  // "cosmic" and GetParticle(0) corresponds to the only particle then
551  srTruth.pdg = tru->GetParticle(0).PdgCode();
552 
553  srTruth.E = tru->GetParticle(0).Momentum().E();
554  srTruth.visE = effPur->energyTotal;
555  srTruth.visEinslc = effPur->energySlice;
556  srTruth.eff = effPur->efficiency;
557  srTruth.pur = effPur->purity;
558  srTruth.nhitslc = effPur->nSliceHits;
559  srTruth.nhittot = effPur->nTotalHits;
560 
561  srTruth.time = tru->GetParticle(0).Position().T();
562  srTruth.vtx = tru->GetParticle(0).Position().Vect();
563  srTruth.p = tru->GetParticle(0).Momentum();
564 
565 
566  srTruth.azimuth = atan2(tru->GetParticle(0).Momentum().Pz(),
567  tru->GetParticle(0).Momentum().Px());
568  srTruth.zenith = -(tru->GetParticle(0).Momentum().Py())/tru->GetParticle(0).Momentum().P();
569 
570  std::vector<const sim::Particle*> particles = bt->MCTruthToParticles(tru);
571 
572 
573  const simb::MCTrajectory& traj = particles[0]->Trajectory();
574 
575 
576  bool entered = 0;
577  bool exited = 0;
578 
579  for(size_t pt = 0; pt < traj.size(); ++pt){
580  TVector3 point = traj.Position(pt).Vect();
581 
582  // Check to see if trajectory point is in the detector
583  bool inDet = fabs(point.X()) <= geom->DetHalfWidth() // Check X inside detector
584  && fabs(point.Y()) <= geom->DetHalfHeight() // then check Y
585  && point.Z() >= 0 && point.Z() <= geom->DetLength(); // then check Z
586 
587  if(inDet && !entered){
588  // Then we just entered the detector
589  entered = true;
590 
591  if(pt > 0){
592  // Put xyz from last point into an array for geo function call
593  TVector3 lastPoint = traj.Position(pt - 1).Vect();
594  double pointArray[3] = {point.X(), point.Y(), point.Z()};
595 
596  // Put xyz from direction vector into an array.
597  TVector3 dir = (lastPoint - point).Unit();
598  double dirArray[3] = {dir.X(), dir.Y(), dir.Z()};
599 
600  double enterArray[3] = {0,0,0};
601 
602  geo::ProjectToBoxEdge(pointArray, dirArray,
603  - geom->DetHalfWidth() , geom->DetHalfWidth(),
604  - geom->DetHalfHeight() , geom->DetHalfHeight(),
605  0 , geom->DetLength(),
606  enterArray);
607 
608  srTruth.enter.SetXYZ(enterArray[0], enterArray[1], enterArray[2]);
609  // No good way to interpolate between momentum at pt-1 and pt.
610  srTruth.penter = traj.Momentum(pt-1);
611  }
612  else{ // If this is the case, we started inside the detector.
613  srTruth.enter = point;
614  srTruth.penter = traj.Momentum(pt);
615  }
616  }
617 
618  if(!inDet && entered && !exited){
619  // Then we just exited the detector
620  exited= true;
621 
622  // Put xyz from last point into an array for geo function call
623  TVector3 lastPoint = traj.Position(pt - 1).Vect();
624  double lastPointArray[3] = {lastPoint.X(), lastPoint.Y(), lastPoint.Z()};
625 
626  // Put xyz from direction vector into an array.
627  TVector3 dir = (point - lastPoint).Unit();
628  double dirArray[3] = {dir.X(), dir.Y(), dir.Z()};
629 
630  double exitArray[3] = {0,0,0};
631 
632  geo::ProjectToBoxEdge(lastPointArray, dirArray,
633  - geom->DetHalfWidth() , geom->DetHalfWidth(),
634  - geom->DetHalfHeight() , geom->DetHalfHeight(),
635  0 , geom->DetLength(),
636  exitArray);
637 
638  srTruth.exit.SetXYZ(exitArray[0], exitArray[1], exitArray[2]);
639  }
640  }
641  srTruth.stop = traj.Position(traj.size() - 1).Vect();
642 
643  // Make sure the particle actually entered the detector
644  if(!entered){
645  // If the primary never entered the detector, set the
646  // coordinates to -5*10^9
647 
648  srTruth.enter.SetXYZ(-5e9, -5e9, -5e9);
649  srTruth.exit.SetXYZ(-5e9, -5e9, -5e9);
650  }
651 
652  // If the particle entered and did not exit, set the exit position to stop
653  if(entered && !exited){
654  srTruth.exit = srTruth.stop;
655  }
656  // Find the Michel electrons
657  FindAndAddMichels(particles, allTracks, &srTruth.michel);
658  }
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
stan::math::fabs
fvar< T > fabs(const fvar< T > &x)
Definition: fabs.hpp:15
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:54
cheat::NeutrinoEffPur::energyTotal
double energyTotal
Sum of FLS hits from the neutrino contributing to any hit in the event.
Definition: BackTracker.h:52
cheat::NeutrinoEffPur::efficiency
double efficiency
Efficiency (based on FLS energy) of neutrino interaction relative to slice.
Definition: BackTracker.h:48
caf::FindAndAddMichels
void FindAndAddMichels(std::vector< const sim::Particle * > particles, std::vector< cheat::TrackIDE > &allTracks, std::vector< SRTrueMichelE > *michelVec)
Definition: FillTruth.cxx:661
vec
Eigen::VectorXd vec
Definition: rvalue_return_test.cpp:23
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
simb::MCTruth::GetParticle
const simb::MCParticle & GetParticle(int i) const
Definition: MCTruth.h:75
cheat::NeutrinoEffPur::nSliceHits
int nSliceHits
Number of hits from this neutrino in this slice.
Definition: BackTracker.h:53
gen_flatrecord.pt
pt
Definition: gen_flatrecord.py:160
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:51
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:47
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:527
cheat::NeutrinoEffPur::purity
double purity
Purity (based on FLS energy) of neutrino interaction relative to slice.
Definition: BackTracker.h:49
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 255 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().

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

153  {
154  vec->push_back(SRSlcME());
155  SRSlcME& srME = vec->back();
156 
157  srME.mid = michel.MID();
158  srME.nhitx = michel.NXCell();
159  srME.nhity = michel.NYCell();
160  srME.adc = michel.TotalADC();
161  srME.calE = michel.CalorimetricEnergy();
162  srME.deltat = michel.DeltaT();
163  srME.disttoslc = michel.DistToSlc();
164  srME.meanpos = TVector3(michel.MeanX(),michel.MeanY(),michel.MeanZ());
165 
166  srME.truth = FillParticleTruth(sliceList,
167  michel.AllCells(),
168  allTracks,
169  sliceTracks,
170  allTracksBirks,
171  sliceTracksBirks,
172  sliceIdx);
173  }
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:25
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
vec
Eigen::VectorXd vec
Definition: rvalue_return_test.cpp:23
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 31 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().

40  {
41  vec->push_back(SRTrkME());
42  SRTrkME& srME = vec->back();
43 
44  srME.nhitx = michel.NXCell();
45  srME.nhity = michel.NYCell();
46  srME.mid = michel.MID();
47  srME.calE = michel.CalorimetricEnergy();
48  srME.adc = michel.TotalADC();
49  srME.deltat = michel.DeltaT();
50  srME.disttotrack = michel.DistToTrk();
51  srME.meanpos = TVector3(michel.MeanX(),michel.MeanY(),michel.MeanZ());
52  srME.truth = FillParticleTruth(sliceList,
53  michel.AllCells(),
54  allTracks,
55  sliceTracks,
56  allTracksBirks,
57  sliceTracksBirks,
58  sliceIdx);
59 
60  srTrk.me.push_back(srME); // Also copy this TrkME to parent track
61  }
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:25
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
vec
Eigen::VectorXd vec
Definition: rvalue_return_test.cpp:23
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 695 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().

699  {
700  vec->push_back(SRTrueMichelE());
701  SRTrueMichelE& srMichel = vec->back();
702 
703  // Find the visible energy from the FLS hits
704  float visE = 0;
705  for (std::vector<cheat::TrackIDE>::iterator it = allTracks.begin();
706  it != allTracks.end(); ++it)
707  {
708  if(michel.TrackId() == it->trackID)
709  {
710  visE = it->energy;
711  break;
712  }
713  }
714 
715  srMichel.E = michel.E();
716  srMichel.visE = visE;
717  srMichel.time = michel.T();
718  srMichel.mustop = motherMuon.EndPosition();
719  srMichel.p = michel.Momentum();
720  }
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
vec
Eigen::VectorXd vec
Definition: rvalue_return_test.cpp:23
simb::MCParticle::Momentum
const TLorentzVector & Momentum(const int i=0) const
Definition: MCParticle.h:219
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:21
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:53
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:44
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:52
caf::SRIDBranch::cvnloosepreselptp
SRCVNResult cvnloosepreselptp
Output from CVN - Loose Presel plus PtP cut (many-class PID)
Definition: SRIDBranch.h:51
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:43
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:57
caf::SRTrackBranch::kalman
SRKalman kalman
Tracks produced by KalmanTrack.
Definition: SRTrackBranch.h:24
novadaq::cnv::kFARDET
Far Detector at Ash River, MN.
Definition: DAQConventions.h:31
caf::SRIDBranch::remid
SRRemid remid
Output from RecoMuonID (ReMId) package.
Definition: SRIDBranch.h:37
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:38
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:39
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:40
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, cet::sqlite::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  }
caf::blind_val
static const int blind_val
Definition: Blinding.cxx:16
submit_hadd.u
u
Definition: submit_hadd.py:83
caf::blind_val_unsigned
static const unsigned int blind_val_unsigned
Definition: Blinding.cxx:17
cet::sqlite::max
T max(sqlite3 *const db, std::string const &table_name, std::string const &column_name)
Definition: statistics.h:68
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 786 of file FillTruth.cxx.

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

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

787  {
788  switch (simbGeneratorEnum)
789  {
790  case simb::Generator_t::kGENIE:
791  return kGENIE;
792 
793  case simb::Generator_t::kGIBUU:
794  return kGIBUU;
795 
796  default:
797  return kUnknownGenerator;
798 
799  }
800  }
void caf::CopyRemidVars ( const SRKalmanTrack &  srTrk,
caf::SRRemid &  remid 
)

Definition at line 659 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().

660  {
661  remid.pid = srTrk.rempid;
662  remid.scatllh = srTrk.scatllh;
663  remid.dedxllh = srTrk.dedxllh;
664  remid.len = srTrk.len;
665  remid.measfrac = srTrk.measfrac;
666  remid.cont = srTrk.remcont;
667  }
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 803 of file FillTruth.cxx.

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

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

804  {
805  std::vector<unsigned int> ret;
806  std::regex pattern("(\\d+)");
807  auto vals_begin = std::sregex_iterator(versionString.begin(), versionString.end(), pattern);
808  auto vals_end = std::sregex_iterator();
809 
810  for (auto it = vals_begin; it != vals_end; it++)
811  ret.push_back(std::stoi(it->str()));
812 
813  return ret;
814  }
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:283
void caf::FillCosRejVars ( const cosrej::CosRejObj &  cosrej,
caf::SRCosRej &  srcosrej,
caf::SRContain &  srcontain 
)

Definition at line 360 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().

363  {
364  srcosrej.anglekal = cosrej.AngleKal();
365  srcosrej.anglecos = cosrej.AngleCos();
366 
367  srcosrej.numucontpid2020 = cosrej.ConCosPID();
368  srcosrej.numucontpid2019 = cosrej.OldCosPID();
369  srcosrej.numuunconttunedpid = cosrej.UnconTunedCosPID();
370 
371  srcosrej.nkal3d = cosrej.NTracks3D();
372  srcosrej.mindist = cosrej.MinDist();
373  srcosrej.mincell = cosrej.MinCell();
374  srcosrej.eratio = cosrej.ERatio();
375  srcosrej.scatt = cosrej.Scatt();
376  srcosrej.hadE = cosrej.StartAct();
377  srcosrej.costhetatrue = cosrej.CosThetaTrue();
378  srcosrej.kalthetatrue = cosrej.KalThetaTrue();
379  srcosrej.pdgbest = cosrej.PDGBest();
380  srcosrej.fscattmax = cosrej.FScattMax();
381  srcosrej.fscattsum = cosrej.FScattSum();
382  srcosrej.fscattext = cosrej.FScattExt();
383  srcosrej.fscattsig = cosrej.FScattSig();
384  srcosrej.kfitspeed = cosrej.KFitSpeed();
385  srcosrej.cfitspeed = cosrej.CFitSpeed();
386  srcosrej.kdirscore = cosrej.KDirScore();
387  srcosrej.cdirscore = cosrej.CDirScore();
388  srcosrej.kscorediff = cosrej.KScoreDiff();
389  srcosrej.cscorediff = cosrej.CScoreDiff();
390  srcosrej.kalchisq = cosrej.KalChisq();
391  srcosrej.kalslope = cosrej.KalSlope();
392  srcosrej.kalchidiff = cosrej.KalChiDiff();
393  srcosrej.coschisq = cosrej.CosChisq();
394  srcosrej.cosslope = cosrej.CosSlope();
395  srcosrej.coschidiff = cosrej.CosChiDiff();
396 
397  srcontain.vtxdist = cosrej.VtxDist();
398  srcontain.enddist = cosrej.EndDist();
399  srcontain.cosfwddist = cosrej.CosFwdDist();
400  srcontain.cosbakdist = cosrej.CosBakDist();
401  srcontain.cosfwdcell = cosrej.CosFwdCell();
402  srcontain.cosfwdcellnd = cosrej.CosFwdCellND();
403  srcontain.cosbakcell = cosrej.CosBakCell();
404  srcontain.cosbakcellnd = cosrej.CosBakCellND();
405  srcontain.cosfwdair = cosrej.CosFwdAir();
406  srcontain.cosfwdsteel = cosrej.CosFwdSteel();
407  srcontain.cosbakair = cosrej.CosBakAir();
408  srcontain.cosbaksteel = cosrej.CosBakSteel();
409  srcontain.kalfwdair = cosrej.KalFwdAir();
410  srcontain.kalfwdsteel = cosrej.KalFwdSteel();
411  srcontain.kalbakair = cosrej.KalBakAir();
412  srcontain.kalbaksteel = cosrej.KalBakSteel();
413  srcontain.kalfwddist = cosrej.KalFwdDist();
414  srcontain.kalbakdist = cosrej.KalBakDist();
415  srcontain.kalfwdcell = cosrej.KalFwdCell();
416  srcontain.kalfwdcellnd = cosrej.KalFwdCellND();
417  srcontain.kalbakcell = cosrej.KalBakCell();
418  srcontain.kalbakcellnd = cosrej.KalBakCellND();
419  srcontain.cosyposattrans = cosrej.CosYPosAtTrans();
420  srcontain.kalyposattrans = cosrej.KalYPosAtTrans();
421  }
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 646 of file FillPIDs.cxx.

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

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

649  {
650  srcvnfeatures.output = features.fOutput;
651  srcvnfeatures.noutput = features.fOutput.size();
652  srcvnfeatures.components = features.fPrincipalComponents;
653  if(int(srcvnfeatures.components.size()) > maxcomponents)
654  srcvnfeatures.components.resize(maxcomponents);
655  srcvnfeatures.ncomponents = srcvnfeatures.components.size();
656  }
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::FillCVNParticleResultVars ( const std::vector< art::Ptr< rb::PID > > &  cvnparts,
caf::SRCVNParticleResult &  cvnpart 
)

Definition at line 89 of file FillPIDs.cxx.

References caf::SRCVNParticleResult::electronid, 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().

91  {
92 
93  float largestv = 0;
94  int largestid = 0;
95  for (unsigned int i=0; i<cvnparts.size(); ++i){
96  float val = cvnparts[i]->Value();
97  int id = cvnparts[i]->Pdg();
98  if (val > largestv){
99  largestv = val;
100  largestid = id;
101  }
102  switch (id) {
103  case 11:
104  cvnpart.electronid = val;
105  break;
106  case 13:
107  cvnpart.muonid = val;
108  break;
109  case 2212:
110  cvnpart.protonid = val;
111  break;
112  case 2112:
113  cvnpart.neutronid = val;
114  break;
115  case 211:
116  cvnpart.pionid = val;
117  break;
118  case 111:
119  cvnpart.pizeroid = val;
120  break;
121  case 22:
122  cvnpart.photonid = val;
123  break;
124  default:
125  cvnpart.otherid = val;
126  break;
127  }
128  }
129  cvnpart.maxval = largestv;
130  cvnpart.pdgmax = largestid;
131 
132  }
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:28
caf::SRCVNParticleResult::protonid
float protonid
Likelihood proton.
Definition: SRCVNParticleResult.h:21
caf::SRCVNParticleResult::maxval
float maxval
Maximum likelihood among net outputs.
Definition: SRCVNParticleResult.h:27
caf::SRCVNParticleResult::otherid
float otherid
Likelihood other.
Definition: SRCVNParticleResult.h:26
caf::SRCVNParticleResult::pizeroid
float pizeroid
Likelihood pizero.
Definition: SRCVNParticleResult.h:24
febshutoff_auto.val
val
Definition: febshutoff_auto.py:277
void caf::FillCVNPixelMaps ( const cvn::PixelMap &  map,
caf::SRPixelObjMap &  srcvnmap,
bool  useGeV 
)

Definition at line 564 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().

566  {
567  srcvnmap.nchan = 2; //two views to pixel map
568  srcvnmap.nplanes = map.fNPlane/2;
569  srcvnmap.ncells = map.fNCell;
570  srcvnmap.hitfracx = map.fFracX;
571  srcvnmap.hitfracy = map.fFracY;
572  srcvnmap.firstplane = map.fBound.FirstPlane();
573  srcvnmap.lastplane = map.fBound.LastPlane();
574  srcvnmap.firstcellx = map.fBound.FirstCell(0);
575  srcvnmap.lastcellx = map.fBound.LastCell(0);
576  srcvnmap.firstcelly = map.fBound.FirstCell(1);
577  srcvnmap.lastcelly = map.fBound.LastCell(1);
578 
579  std::vector<unsigned char> tval = map.PixelMapToVector(useGeV);
580  std::copy(tval.begin(), tval.end(), srcvnmap.cvnmap);
581  std::vector<unsigned char> tlab = map.PixelMapToObjVector(true);
582  std::copy(tlab.begin(), tlab.end(), srcvnmap.cvnlabmap);
583  std::vector<unsigned char> tobj = map.PixelMapToObjVector(false);
584  std::copy(tobj.begin(), tobj.end(), srcvnmap.cvnobjmap);
585  }
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 588 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().

590  {
591  srcvnmap.nchan = 2; //two views to pixel map
592  srcvnmap.nplanes = map.fNPlane/2;
593  srcvnmap.ncells = map.fNCell;
594  srcvnmap.hitfracx = map.fFracX;
595  srcvnmap.hitfracy = map.fFracY;
596  srcvnmap.firstplane = map.fBound.FirstPlane();
597  srcvnmap.lastplane = map.fBound.LastPlane();
598  srcvnmap.firstcellx = map.fBound.FirstCell(0);
599  srcvnmap.lastcellx = map.fBound.LastCell(0);
600  srcvnmap.firstcelly = map.fBound.FirstCell(1);
601  srcvnmap.lastcelly = map.fBound.LastCell(1);
602 
603  std::vector<unsigned char> tval = map.PixelMapToVector(useGeV);
604  std::copy(tval.begin(), tval.end(), srcvnmap.cvnmap);
605  }
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 624 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().

626  {
627  srpdata.label3d = cvnpdata.fLabel3D;
628  srpdata.labelx = cvnpdata.fLabelX;
629  srpdata.labely = cvnpdata.fLabelY;
630  srpdata.isprimary = cvnpdata.fIsPrimary;
631  srpdata.purity3d = cvnpdata.fPurity3D;
632  srpdata.purityx = cvnpdata.fPurityX;
633  srpdata.purityy = cvnpdata.fPurityY;
634  srpdata.rece = cvnpdata.fRecE;
635  srpdata.ncellx = cvnpdata.fNCellX;
636  srpdata.ncelly = cvnpdata.fNCellY;
637  srpdata.vertx = cvnpdata.fVertX;
638  srpdata.verty = cvnpdata.fVertY;
639  srpdata.vertz = cvnpdata.fVertZ;
640  srpdata.plength = cvnpdata.fPLength;
641  srpdata.pangle = cvnpdata.fPAngle;
642  srpdata.pgap = cvnpdata.fPGap;
643  }
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 542 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().

545  {
546  srcvn.output = result.fOutput;
547 
548  std::vector<float>::const_iterator maxIt =
549  std::max_element(result.fOutput.cbegin(), result.fOutput.cend());
550 
551  srcvn.maxval = *maxIt;
552  srcvn.argmax = std::distance(result.fOutput.cbegin(), maxIt);
553  srcvn.noutput = result.fOutput.size();
554 
555  srcvn.numuid = result.fOutput[0];
556  srcvn.nueid = result.fOutput[1];
557  srcvn.nutauid = -5;
558  srcvn.ncid = result.fOutput[2];
559  srcvn.cosmicid = result.fOutput[3];
560  if (NoCosmics) srcvn.cosmicid = -5;
561  }
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 608 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().

610  {
611  srtdata.parent = cvntdata.fPPt;
612  srtdata.interaction = cvntdata.fInt;
613  srtdata.finalstate = cvntdata.fFSt;
614  srtdata.finalstateprong = cvntdata.fFSp;
615  srtdata.particles = cvntdata.fFSlong;
616  srtdata.nuenergy = cvntdata.fNuEnergy;
617  srtdata.lepenergy = cvntdata.fLepEnergy;
618  srtdata.vtxx = cvntdata.fVtxX;
619  srtdata.vtxy = cvntdata.fVtxY;
620  srtdata.vtxz = cvntdata.fVtxZ;
621  }
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 501 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().

508  {
509  srMR.isShwDiF = IsTrueDiF(shw);
510 
511  srMR.lid.nhit = shw.NCell();
512  srMR.lid.nhitx = shw.NXCell();
513  srMR.lid.nhity = shw.NYCell();
514  srMR.lid.nplane = shw.ExtentPlane();
515  srMR.lid.nplanex = shw.ExtentPlane(geo::kX);
516  srMR.lid.nplaney = shw.ExtentPlane(geo::kY);
517  srMR.lid.maxplanecont = shw.MostContiguousPlanes(geo::kXorY);
518  srMR.lid.maxplanegap = shw.MostMissingPlanes(geo::kXorY);
519  srMR.lid.len = shw.TotalLength();
520  srMR.lid.calE = shw.CalorimetricEnergy();
521  srMR.lid.start = shw.Start();
522  srMR.lid.dir = shw.Dir();
523  srMR.lid.stop = shw.Stop();
524 
525  srMR.lid.truth = FillParticleTruth(sliceList,
526  shw.AllCells(),
527  allTracks,
528  sliceTracks,
529  allTracksBirks,
530  sliceTracksBirks,
531  sliceIdx);
532  srMR.lid.truthXView = FillParticleTruth(sliceList,
533  shw.XCells(),
534  allTracks,
535  sliceTracks,
536  allTracksBirks,
537  sliceTracksBirks,
538  sliceIdx);
539  srMR.lid.truthYView = FillParticleTruth(sliceList,
540  shw.YCells(),
541  allTracks,
542  sliceTracks,
543  allTracksBirks,
544  sliceTracksBirks,
545  sliceIdx);
546  if (shw.NXCell() == 0) srMR.lid.view = kY;
547  else if (shw.NYCell() == 0) srMR.lid.view = kX;
548  else srMR.lid.view = kXorY;
549  }
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:25
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:64
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 94 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().

102  {
103  srMR.isDiF = IsTrueDiF(slice);
104 
105  srMR.nhit = slice.NCell();
106  srMR.nhitx = slice.NXCell();
107  srMR.nhity = slice.NYCell();
108  srMR.nplane = slice.ExtentPlane();
109  srMR.nplanex = slice.ExtentPlane(geo::kX);
110  srMR.nplaney = slice.ExtentPlane(geo::kY);
111  srMR.maxplanecont = slice.MostContiguousPlanes(geo::kXorY);
112  srMR.maxplanegap = slice.MostMissingPlanes(geo::kXorY);
113  srMR.minplane = slice.MinPlane(geo::kXorY);
114  srMR.maxplane = slice.MaxPlane(geo::kXorY);
115  srMR.calE = slice.CalorimetricEnergy();
116 
117  srMR.truth = FillParticleTruth(sliceList,
118  slice.AllCells(),
119  allTracks,
120  sliceTracks,
121  allTracksBirks,
122  sliceTracksBirks,
123  sliceIdx);
124  srMR.truthXView = FillParticleTruth(sliceList,
125  slice.XCells(),
126  allTracks,
127  sliceTracks,
128  allTracksBirks,
129  sliceTracksBirks,
130  sliceIdx);
131  srMR.truthYView = FillParticleTruth(sliceList,
132  slice.YCells(),
133  allTracks,
134  sliceTracks,
135  allTracksBirks,
136  sliceTracksBirks,
137  sliceIdx);
138  if (slice.NXCell() == 0) srMR.view = kY;
139  else if (slice.NYCell() == 0) srMR.view = kX;
140  else srMR.view = kXorY;
141  }
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:25
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:64
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 103 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, 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().

105  {
106  static art::ServiceHandle<calib::Calibrator> cal;
107  static art::ServiceHandle<geo::Geometry> geom;
108 
109  const double wx = (hadclust.NXCell() > 0) ? hadclust.W(hadclust.XCell(0).get()) : 0;
110  const double wy = (hadclust.NYCell() > 0) ? hadclust.W(hadclust.YCell(0).get()) : 0;
111 
112  // Variables needed:
113  // Holds contiguous planes info:
114  unsigned int cellsFromEdge = -1;
115 
116  // Hit counters:
117  int nCalHit = 0;
118  int nMipHit = 0;
119 
120  unsigned int minCellsFromEdge = 99999999;
121 
122  // Loop over hadclust in question here:
123  for(unsigned int hitIdx = 0; hitIdx < hadclust.NCell(); ++hitIdx)
124  {
125  const art::Ptr<rb::CellHit>& chit = hadclust.Cell(hitIdx);
126  // Get calbrated reco hit:
127  const rb::RecoHit rhit(cal->MakeRecoHit(*chit,
128  chit->View() == geo::kX ? wx : wy));
129 
130  const int planeNum = chit->Plane();
131  cellsFromEdge = std::min((unsigned int)chit->Cell(), geom->Plane(planeNum)->Ncells() - 1 - chit->Cell());
132 
133  if(cellsFromEdge < minCellsFromEdge)
134  {
135  minCellsFromEdge = cellsFromEdge;
136  }
137 
138  // Make sure things have been calibrated:
139  if(!rhit.IsCalibrated())
140  {
141  // LOG_DEBUG("CAFMaker") << "Not calibrated?! "
142  // << chit->Plane() << " " << chit->Cell();
143  continue;
144  }
145 
146  // Count up both the total calibrated hits
147  // and the MIP hits from 100-245 PECorr
148  // Count up total hits:
149  nCalHit++;
150  // Count up MIP hits:
151  if ( (rhit.PECorr()>100.0) && (rhit.PECorr()<245.0) ) nMipHit++;
152  } // end of hadclust loop
153 
154  srhadclust.nhit = hadclust.NCell();
155  srhadclust.ncalhit = nCalHit;
156  srhadclust.nmiphit = nMipHit;
157  srhadclust.ncontplanes = hadclust.MostContiguousPlanes(geo::kXorY);
158  if (!(hadclust.NCell()==0)){
159 
160  srhadclust.firstplane = hadclust.MinPlane();
161  srhadclust.lastplane = hadclust.MaxPlane();
162  srhadclust.firstcell = hadclust.MinCell(geo::kXorY);
163  srhadclust.lastcell = hadclust.MaxCell(geo::kXorY);
164  srhadclust.ncellsfromedge = minCellsFromEdge;
165 
166  srhadclust.calE = hadclust.CalorimetricEnergy();
167 
168  TVector3 boxmin, boxmax, meanpos;
169  hadclust.MinMaxMeanXYZ(boxmin, boxmax, meanpos);
170  srhadclust.boxmin = boxmin;
171  srhadclust.boxmax = boxmax;
172  srhadclust.meanpos = meanpos;
173 
174  // TBrowser can't draw DBL_MAX
175  if (hadclust.NXCell()==0){
176  srhadclust.boxmin.SetX(-5555);
177  srhadclust.boxmax.SetX(5555);
178  }
179  if (hadclust.NYCell()==0){
180  srhadclust.boxmin.SetY(-5555);
181  srhadclust.boxmax.SetY(5555);
182  }
183 
184  } // end of loop over variables to only fill for hadclust if it has hits
185  else{
186  srhadclust.firstplane = 0;
187  srhadclust.lastplane = 0;
188  srhadclust.firstcell = 0;
189  srhadclust.lastcell = 0;
190  srhadclust.ncellsfromedge = 0;
191 
192  srhadclust.calE = -5;
193 
194  srhadclust.boxmin = TVector3(-5555,-5555,-5555);
195  srhadclust.boxmax = TVector3(-5555,-5555,-5555);
196  srhadclust.meanpos = TVector3(-5555,-5555,-5555);
197  }
198  }
rb::Cluster::W
virtual double W(const rb::CellHit *chit) const
Estimate the unmeasured coordinate of chit.
Definition: Cluster.cxx:121
rb::Cluster::NCell
unsigned int NCell(geo::View_t view) const
Number of cells in view view.
Definition: Cluster.cxx:134
rb::Cluster::XCell
art::Ptr< rb::CellHit > XCell(unsigned int xIdx) const
Get the ith cell in the x-view.
Definition: Cluster.cxx:157
caf::SRHadClust::ncellsfromedge
unsigned int ncellsfromedge
minimum number of cells to edge of detector
Definition: SRHadClust.h:30
caf::SRVector3D::SetX
void SetX(float _x)
Definition: SRVector3D.h:28
geo::kXorY
X or Y views.
Definition: PlaneGeo.h:30
rb::CellHit::Plane
unsigned short Plane() const
Definition: CellHit.h:39
rb::CellHit::View
geo::View_t View() const
Definition: CellHit.h:41
geo::kX
Vertical planes which measure X.
Definition: PlaneGeo.h:28
geo::PlaneGeo::Ncells
unsigned int Ncells() const
Number of cells in this plane.
Definition: PlaneGeo.h:43
calib::Calibrator::MakeRecoHit
rb::RecoHit MakeRecoHit(rb::CellHit const &cellhit, double w)
Definition: Calibrator_service.cc:575
caf::SRHadClust::ncontplanes
unsigned int ncontplanes
number of continuous planes
Definition: SRHadClust.h:25
rb::Cluster::MaxCell
unsigned int MaxCell(geo::View_t view) const
Definition: Cluster.cxx:518
caf::SRHadClust::meanpos
SRVector3D meanpos
Mean position of hits in cluster, weighted by charge [cm].
Definition: SRHadClust.h:34
geo::GeometryBase::Plane
const PlaneGeo * Plane(unsigned int i) const
Definition: GeometryBase.cxx:749
caf::SRHadClust::nhit
unsigned int nhit
number of hits
Definition: SRHadClust.h:22
rb::RecoHit
Calibrated quantities relying on position in the orthogonal view. To generate a rb::CellHit from a rb...
Definition: RecoHit.h:19
rb::CellHit::Cell
unsigned short Cell() const
Definition: CellHit.h:40
rb::Cluster::CalorimetricEnergy
double CalorimetricEnergy(EEnergyCalcScheme escheme=kRecomputeEnergy) const
Simple estimate of neutrino energy.
Definition: Cluster.cxx:439
rb::Cluster::YCell
art::Ptr< rb::CellHit > YCell(unsigned int yIdx) const
Get the ith cell in the y-view.
Definition: Cluster.cxx:165
caf::SRHadClust::nmiphit
unsigned int nmiphit
number of minimum ionizing hits
Definition: SRHadClust.h:24
min
static float min(const float a, const float b, const float c)
Definition: absgeo.cxx:45
rb::Cluster::NYCell
unsigned int NYCell() const
Number of cells in the y-view.
Definition: Cluster.h:108
caf::SRHadClust::calE
float calE
Calorimetric energy of the cluster [GeV].
Definition: SRHadClust.h:31
rb::Cluster::MinPlane
unsigned int MinPlane(geo::View_t view=geo::kXorY) const
Definition: Cluster.cxx:462
rb::Cluster::Cell
art::Ptr< rb::CellHit > Cell(geo::View_t view, unsigned int viewIdx) const
Get the ith cell from view view.
Definition: Cluster.cxx:145
art::Ptr::get
T const * get() const
Definition: Ptr.h:321
rb::Cluster::NXCell
unsigned int NXCell() const
Number of cells in the x-view.
Definition: Cluster.h:106
rb::Cluster::MinCell
unsigned int MinCell(geo::View_t view) const
Definition: Cluster.cxx:472
rb::Cluster::MostContiguousPlanes
int MostContiguousPlanes(geo::View_t view) const
Longest run of adjacent planes with hits.
Definition: Cluster.cxx:635
geom
void geom(int which=0)
Definition: geom.C:163
caf::SRHadClust::lastcell
unsigned int lastcell
last cell
Definition: SRHadClust.h:29
caf::SRHadClust::ncalhit
unsigned int ncalhit
number of hits with calibration
Definition: SRHadClust.h:23
rb::Cluster::MaxPlane
unsigned int MaxPlane(geo::View_t view=geo::kXorY) const
Definition: Cluster.cxx:508
rb::Cluster::MinMaxMeanXYZ
void MinMaxMeanXYZ(TVector3 &lo, TVector3 &hi, TVector3 &mean, rb::AveragingScheme scheme=kDefaultScheme) const
Gets the min/max/mean all at once, called by the functions above.
Definition: Cluster.cxx:582
caf::SRHadClust::boxmin
SRVector3D boxmin
Minimum coordinates box containing all the hits [cm].
Definition: SRHadClust.h:32
caf::SRHadClust::firstcell
unsigned int firstcell
first cell
Definition: SRHadClust.h:28
caf::SRVector3D::SetY
void SetY(float _y)
Definition: SRVector3D.h:29
caf::SRHadClust::firstplane
unsigned int firstplane
first plane
Definition: SRHadClust.h:26
caf::SRHadClust::lastplane
unsigned int lastplane
last plane
Definition: SRHadClust.h:27
caf::SRHadClust::boxmax
SRVector3D boxmax
Maximum coordinates box containing all the hits [cm].
Definition: SRHadClust.h:33
void caf::FillJMEIDVars ( const jmshower::EID &  eid,
caf::SRJMEid &  sreid 
)

Definition at line 47 of file FillPIDs.cxx.

References caf::SRJMEid::ann, jmshower::EID::ANN(), caf::SRJMEid::annE, jmshower::EID::ANNE(), caf::SRJMEid::annEL, jmshower::EID::ANNEL(), caf::SRJMEid::annENoCos, jmshower::EID::ANNENoCos(), caf::SRJMEid::annEPi0, jmshower::EID::ANNEPi0(), caf::SRJMEid::annEPi0EL, jmshower::EID::ANNEPi0EL(), caf::SRJMEid::annNoCos, jmshower::EID::ANNNoCos(), caf::SRJMEid::dedx0, jmshower::EID::Dedx0(), caf::SRJMEid::dedx1, jmshower::EID::Dedx1(), caf::SRJMEid::dedx2, jmshower::EID::Dedx2(), caf::SRJMEid::dedx3, jmshower::EID::Dedx3(), caf::SRJMEid::dedx4, jmshower::EID::Dedx4(), caf::SRJMEid::dedx5, jmshower::EID::Dedx5(), caf::SRJMEid::eellll, jmshower::EID::EelLLL(), caf::SRJMEid::eelllt, jmshower::EID::EelLLT(), caf::SRJMEid::eglll, jmshower::EID::EGLLL(), caf::SRJMEid::egllt, jmshower::EID::EGLLT(), caf::SRJMEid::elll, jmshower::EID::ELLL(), caf::SRJMEid::ellt, jmshower::EID::ELLT(), caf::SRJMEid::emulll, jmshower::EID::EMuLLL(), caf::SRJMEid::emullt, jmshower::EID::EMuLLT(), caf::SRJMEid::enlll, jmshower::EID::ENLLL(), caf::SRJMEid::enllt, jmshower::EID::ENLLT(), caf::SRJMEid::epi0lll, jmshower::EID::EPi0LLL(), caf::SRJMEid::epi0llt, jmshower::EID::EPi0LLT(), caf::SRJMEid::epilll, jmshower::EID::EPiLLL(), caf::SRJMEid::epillt, jmshower::EID::EPiLLT(), caf::SRJMEid::eplll, jmshower::EID::EPLLL(), caf::SRJMEid::epllt, jmshower::EID::EPLLT(), caf::SRJMEid::gap, jmshower::EID::Gap(), caf::SRJMEid::ismuon, jmshower::EID::IsMuon(), caf::SRJMEid::mulll, jmshower::EID::MuLLL(), caf::SRJMEid::mullt, jmshower::EID::MuLLT(), caf::SRJMEid::pi0mass, jmshower::EID::Pi0Mass(), jmshower::EID::ShE(), caf::SRJMEid::shwEFrac, caf::SRJMEid::vtxgev, and jmshower::EID::VtxGeV().

49  {
50  sreid.ann = eid.ANN();
51  sreid.annE = eid.ANNE();
52  sreid.annNoCos = eid.ANNNoCos();
53  sreid.annENoCos = eid.ANNENoCos();
54  sreid.annEL = eid.ANNEL();
55  sreid.annEPi0 = eid.ANNEPi0();
56  sreid.annEPi0EL = eid.ANNEPi0EL();
57  sreid.ismuon = eid.IsMuon();
58  sreid.dedx0 = eid.Dedx0();
59  sreid.dedx1 = eid.Dedx1();
60  sreid.dedx2 = eid.Dedx2();
61  sreid.dedx3 = eid.Dedx3();
62  sreid.dedx4 = eid.Dedx4();
63  sreid.dedx5 = eid.Dedx5();
64  sreid.eglll = eid.EGLLL();
65  sreid.egllt = eid.EGLLT();
66  sreid.emulll = eid.EMuLLL();
67  sreid.emullt = eid.EMuLLT();
68  sreid.epi0lll = eid.EPi0LLL();
69  sreid.epi0llt = eid.EPi0LLT();
70  sreid.eplll = eid.EPLLL();
71  sreid.epllt = eid.EPLLT();
72  sreid.enlll = eid.ENLLL();
73  sreid.enllt = eid.ENLLT();
74  sreid.epilll = eid.EPiLLL();
75  sreid.epillt = eid.EPiLLT();
76  sreid.vtxgev = eid.VtxGeV();
77  sreid.pi0mass = eid.Pi0Mass();
78  sreid.shwEFrac = eid.ShE();
79  sreid.gap = eid.Gap();
80  sreid.elll = eid.ELLL();
81  sreid.ellt = eid.ELLT();
82  sreid.eellll = eid.EelLLL();
83  sreid.eelllt = eid.EelLLT();
84  sreid.mulll = eid.MuLLL();
85  sreid.mullt = eid.MuLLT();
86  }
jmshower::EID::ANNEL
double ANNEL() const
Definition: EID.cxx:34
jmshower::EID::ANNNoCos
double ANNNoCos() const
Definition: EID.cxx:61
caf::SRJMEid::annENoCos
float annENoCos
ann output, with E without cosTheta
Definition: SRJMEid.h:27
jmshower::EID::ELLL
double ELLL() const
Definition: EID.cxx:798
caf::SRJMEid::gap
float gap
gap from vertex to start of shower
Definition: SRJMEid.h:57
jmshower::EID::EPiLLT
double EPiLLT() const
Definition: EID.cxx:748
caf::SRJMEid::dedx4
float dedx4
Plane Dedx.
Definition: SRJMEid.h:40
caf::SRJMEid::ann
float ann
ann output
Definition: SRJMEid.h:21
jmshower::EID::MuLLT
double MuLLT() const
Definition: EID.cxx:850
jmshower::EID::Dedx4
double Dedx4() const
Definition: EID.cxx:622
caf::SRJMEid::epi0llt
float epi0llt
Electron - Pi0 ll for ltransverse shower.
Definition: SRJMEid.h:47
caf::SRJMEid::elll
float elll
electron ll for longitudinal shower
Definition: SRJMEid.h:61
jmshower::EID::ANNEPi0EL
double ANNEPi0EL() const
Definition: EID.cxx:52
jmshower::EID::EGLLL
double EGLLL() const
Definition: EID.cxx:638
jmshower::EID::EPiLLL
double EPiLLL() const
Definition: EID.cxx:738
caf::SRJMEid::emullt
float emullt
Electron - muon ll for ltransverse shower.
Definition: SRJMEid.h:45
caf::SRJMEid::eglll
float eglll
Plane Dedx.
Definition: SRJMEid.h:42
caf::SRJMEid::vtxgev
float vtxgev
Energy of slice in vertex region.
Definition: SRJMEid.h:54
jmshower::EID::IsMuon
bool IsMuon() const
Definition: EID.cxx:547
caf::SRJMEid::epillt
float epillt
Electron - pion ll for ltransverse shower.
Definition: SRJMEid.h:53
caf::SRJMEid::eplll
float eplll
Electron - proton ll for longitudinal shower.
Definition: SRJMEid.h:48
jmshower::EID::EelLLT
double EelLLT() const
Definition: EID.cxx:828
jmshower::EID::ENLLL
double ENLLL() const
Definition: EID.cxx:718
caf::SRJMEid::epllt
float epllt
Electron - proton ll for ltransverse shower.
Definition: SRJMEid.h:49
jmshower::EID::Dedx3
double Dedx3() const
Definition: EID.cxx:615
caf::SRJMEid::epilll
float epilll
Electron - pion ll for longitudinal shower.
Definition: SRJMEid.h:52
jmshower::EID::EMuLLL
double EMuLLL() const
Definition: EID.cxx:658
jmshower::EID::Dedx1
double Dedx1() const
Definition: EID.cxx:601
caf::SRJMEid::dedx3
float dedx3
Plane Dedx.
Definition: SRJMEid.h:39
jmshower::EID::Gap
double Gap() const
Definition: EID.cxx:788
caf::SRJMEid::annNoCos
float annNoCos
ann output, without cosTheta
Definition: SRJMEid.h:26
jmshower::EID::Dedx2
double Dedx2() const
Definition: EID.cxx:608
jmshower::EID::ShE
double ShE() const
Definition: EID.cxx:778
caf::SRJMEid::ellt
float ellt
electron ll for transverse shower
Definition: SRJMEid.h:62
jmshower::EID::ANN
double ANN() const
Definition: EID.cxx:16
jmshower::EID::EMuLLT
double EMuLLT() const
Definition: EID.cxx:668
jmshower::EID::ANNENoCos
double ANNENoCos() const
Definition: EID.cxx:70
caf::SRJMEid::ismuon
int ismuon
Boolean to distinguish muons from electrons.
Definition: SRJMEid.h:28
caf::SRJMEid::annEPi0EL
float annEPi0EL
ann output, EPi0 Elastic scattering
Definition: SRJMEid.h:25
caf::SRJMEid::annEL
float annEL
ann output,
Definition: SRJMEid.h:23
jmshower::EID::ENLLT
double ENLLT() const
Definition: EID.cxx:728
caf::SRJMEid::annE
float annE
ann output, with E
Definition: SRJMEid.h:22
caf::SRJMEid::dedx1
float dedx1
Plane Dedx.
Definition: SRJMEid.h:37
jmshower::EID::EPLLL
double EPLLL() const
Definition: EID.cxx:698
caf::SRJMEid::shwEFrac
float shwEFrac
fraction of energy of leading shower out of total energy of slice
Definition: SRJMEid.h:56
caf::SRJMEid::emulll
float emulll
Electron - muon ll for longitudinal shower.
Definition: SRJMEid.h:44
caf::SRJMEid::enllt
float enllt
Electron - neutron ll for ltransverse shower.
Definition: SRJMEid.h:51
jmshower::EID::ANNE
double ANNE() const
Definition: EID.cxx:25
caf::SRJMEid::enlll
float enlll
Electron - neutron ll for longitudinal shower.
Definition: SRJMEid.h:50
jmshower::EID::MuLLL
double MuLLL() const
Definition: EID.cxx:840
caf::SRJMEid::pi0mass
float pi0mass
best pi0 mass hypothesis coming from combinations of JMShowers
Definition: SRJMEid.h:55
jmshower::EID::EPi0LLT
double EPi0LLT() const
Definition: EID.cxx:688
jmshower::EID::EelLLL
double EelLLL() const
Definition: EID.cxx:818
jmshower::EID::Dedx0
double Dedx0() const
Definition: EID.cxx:594
jmshower::EID::ELLT
double ELLT() const
Definition: EID.cxx:808
jmshower::EID::EPLLT
double EPLLT() const
Definition: EID.cxx:708
caf::SRJMEid::epi0lll
float epi0lll
Electron - Pi0 ll for longitudinal shower.
Definition: SRJMEid.h:46
caf::SRJMEid::eelllt
float eelllt
electron el ll for transverse shower
Definition: SRJMEid.h:64
caf::SRJMEid::dedx0
float dedx0
Definition: SRJMEid.h:36
jmshower::EID::EPi0LLL
double EPi0LLL() const
Definition: EID.cxx:678
caf::SRJMEid::dedx2
float dedx2
Plane Dedx.
Definition: SRJMEid.h:38
jmshower::EID::Dedx5
double Dedx5() const
Definition: EID.cxx:628
caf::SRJMEid::egllt
float egllt
Electron - gamma ll for ltransverse shower.
Definition: SRJMEid.h:43
caf::SRJMEid::mullt
float mullt
muon ll for transverse shower
Definition: SRJMEid.h:66
caf::SRJMEid::dedx5
float dedx5
Plane Dedx.
Definition: SRJMEid.h:41
jmshower::EID::Pi0Mass
double Pi0Mass() const
Definition: EID.cxx:768
caf::SRJMEid::annEPi0
float annEPi0
ann output, EPi0
Definition: SRJMEid.h:24
jmshower::EID::VtxGeV
double VtxGeV() const
Definition: EID.cxx:758
jmshower::EID::ANNEPi0
double ANNEPi0() const
Definition: EID.cxx:43
caf::SRJMEid::mulll
float mulll
muon ll for longitudinal shower
Definition: SRJMEid.h:65
jmshower::EID::EGLLT
double EGLLT() const
Definition: EID.cxx:648
caf::SRJMEid::eellll
float eellll
electron el ll for longitudinal shower
Definition: SRJMEid.h:63
void caf::FillLEMVars ( const lem::PIDDetails &  lem,
caf::SRLem &  srlem 
)

Definition at line 282 of file FillPIDs.cxx.

References caf::SRLem::avgexpE, caf::SRLem::avginvE, caf::SRLem::avgsigE, caf::SRLem::chibkg, caf::SRLem::chibkgexp, caf::SRLem::chisig, caf::SRLem::chisigexp, caf::SRLem::dectree, caf::SRLem::energydiff, caf::SRLem::energydiffexp, caf::SRLem::enrichfrac, caf::SRLem::enrichfracexp, caf::SRLem::fannid, caf::SRLem::fannidenrich, lem::PIDExtraVars::fAvgExpE, lem::PIDExtraVars::fAvgInvE, lem::PIDExtraVars::fAvgSigE, lem::PIDExtraVars::fChiBkg, lem::PIDExtraVars::fChiBkgExp, lem::PIDExtraVars::fChiSig, lem::PIDExtraVars::fChiSigExp, lem::PIDExtraVars::fDecTreePID, lem::PIDExtraVars::fEnergyDiff, lem::PIDExtraVars::fEnergyDiffExp, lem::PIDExtraVars::fEnrichFrac, lem::PIDExtraVars::fEnrichFracExp, lem::PIDExtraVars::fFannId, lem::PIDExtraVars::fFannIdEnrich, lem::PIDExtraVars::fMeanQFrac, lem::PIDExtraVars::fMeanQFracBkg, lem::PIDExtraVars::fMeanQFracBkgExp, lem::PIDExtraVars::fMeanQFracExp, lem::PIDExtraVars::fMeanQFracFit, lem::PIDExtraVars::fMeanQFracSig, lem::PIDExtraVars::fMeanQFracSigExp, lem::PIDExtraVars::fMeanY, lem::PIDExtraVars::fMeanYExp, lem::PIDExtraVars::fMeanYFit, lem::PIDExtraVars::fMeanYSig, lem::PIDExtraVars::fMeanYSigExp, lem::PIDExtraVars::fPIDExp, lem::PIDExtraVars::fPIDExpFit, lem::PIDExtraVars::fPIDFit, lem::PIDExtraVars::fTMVABDT, lem::PIDExtraVars::fTMVABDTD, lem::PIDExtraVars::fTMVABDTG, lem::PIDExtraVars::fTMVABDTMF, lem::PIDExtraVars::fTMVAMLP, caf::SRLem::meanqfrac, caf::SRLem::meanqfracbkg, caf::SRLem::meanqfracbkgexp, caf::SRLem::meanqfracexp, caf::SRLem::meanqfracfit, caf::SRLem::meanqfracsig, caf::SRLem::meanqfracsigexp, caf::SRLem::meany, caf::SRLem::meanyexp, caf::SRLem::meanyfit, caf::SRLem::meanysig, caf::SRLem::meanysigexp, caf::SRLem::pid, caf::SRLem::pidexp, caf::SRLem::pidexpfit, caf::SRLem::pidfit, caf::SRLem::tmvabdt, caf::SRLem::tmvabdtd, caf::SRLem::tmvabdtg, caf::SRLem::tmvabdtmf, caf::SRLem::tmvamlp, rb::PID::Value(), and lem::PIDDetails::vars.

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

284  {
285  srlem.pid = lem.Value();
286  srlem.pidfit = lem.vars.fPIDFit;
287  srlem.pidexp = lem.vars.fPIDExp;
288  srlem.pidexpfit = lem.vars.fPIDExpFit;
289  srlem.meany = lem.vars.fMeanY;
290  srlem.meanyfit = lem.vars.fMeanYFit;
291  srlem.meanyexp = lem.vars.fMeanYExp;
292  srlem.meanysig = lem.vars.fMeanYSig;
293  srlem.meanysigexp = lem.vars.fMeanYSigExp;
294  srlem.meanqfrac = lem.vars.fMeanQFrac;
295  srlem.meanqfracfit = lem.vars.fMeanQFracFit;
296  srlem.meanqfracexp = lem.vars.fMeanQFracExp;
297  srlem.meanqfracsig = lem.vars.fMeanQFracSig;
298  srlem.meanqfracsigexp = lem.vars.fMeanQFracSigExp;
299  srlem.meanqfracbkg = lem.vars.fMeanQFracBkg;
300  srlem.meanqfracbkgexp = lem.vars.fMeanQFracBkgExp;
301  srlem.energydiff = lem.vars.fEnergyDiff;
302  srlem.energydiffexp = lem.vars.fEnergyDiffExp;
303  srlem.enrichfrac = lem.vars.fEnrichFrac;
304  srlem.enrichfracexp = lem.vars.fEnrichFracExp;
305  srlem.fannid = lem.vars.fFannId;
306  srlem.fannidenrich = lem.vars.fFannIdEnrich;
307  srlem.chisig = lem.vars.fChiSig;
308  srlem.chibkg = lem.vars.fChiBkg;
309  srlem.chisigexp = lem.vars.fChiSigExp;
310  srlem.chibkgexp = lem.vars.fChiBkgExp;
311  srlem.avginvE = lem.vars.fAvgInvE;
312  srlem.avgexpE = lem.vars.fAvgExpE;
313  srlem.avgsigE = lem.vars.fAvgSigE;
314  srlem.tmvabdtg = lem.vars.fTMVABDTG;
315  srlem.tmvabdt = lem.vars.fTMVABDT;
316  srlem.tmvabdtd = lem.vars.fTMVABDTD;
317  srlem.tmvabdtmf = lem.vars.fTMVABDTMF;
318  srlem.tmvamlp = lem.vars.fTMVAMLP;
319  srlem.dectree = lem.vars.fDecTreePID;
320  }
caf::SRLem::meanyexp
float meanyexp
Hadronic y of matches, "exp".
Definition: SRLem.h:32
caf::SRLem::pid
float pid
The default PID value - normally use this one.
Definition: SRLem.h:24
caf::SRLem::tmvabdtmf
float tmvabdtmf
TMVA boosted decision tree, Fisher discriminant.
Definition: SRLem.h:65
caf::SRLem::tmvabdt
float tmvabdt
TMVA boosted decision tree, adaptive boost.
Definition: SRLem.h:63
caf::SRLem::energydiffexp
float energydiffexp
Pot. diff between sig and bkg matches, "exp".
Definition: SRLem.h:45
caf::SRLem::meanyfit
float meanyfit
Hadronic y of matches, "fit".
Definition: SRLem.h:31
caf::SRLem::meanqfracbkg
float meanqfracbkg
Mean fraction matched in background matches.
Definition: SRLem.h:41
lem::PIDExtraVars::fMeanQFracBkgExp
double fMeanQFracBkgExp
Definition: PIDExtraVars.h:47
caf::SRLem::enrichfrac
float enrichfrac
Fraction of matches that are "enriched".
Definition: SRLem.h:47
caf::SRLem::chibkgexp
float chibkgexp
chisq of a fit assuming best match is bkg, "exp"
Definition: SRLem.h:56
rb::PID::Value
double Value() const
Definition: PID.h:22
caf::SRLem::tmvabdtd
float tmvabdtd
TMVA boosted decision tree, decorrelation.
Definition: SRLem.h:64
caf::SRLem::meanqfracbkgexp
float meanqfracbkgexp
Fraction matched in background matches, "exp".
Definition: SRLem.h:42
caf::SRLem::meanqfracfit
float meanqfracfit
Fraction of charge matched, "fit".
Definition: SRLem.h:37
caf::SRLem::meanysig
float meanysig
Mean hadronic y of matches that are signal.
Definition: SRLem.h:33
caf::SRLem::fannid
float fannid
ANN output, "fann" package.
Definition: SRLem.h:50
caf::SRLem::meanqfracsig
float meanqfracsig
Fraction matched in signal matches.
Definition: SRLem.h:39
lem::PIDExtraVars::fMeanQFracSigExp
double fMeanQFracSigExp
Definition: PIDExtraVars.h:44
caf::SRLem::pidexp
float pidexp
Fraction of matches that are signal, "exp".
Definition: SRLem.h:27
caf::SRLem::energydiff
float energydiff
Potential diff between sig and bkg matches.
Definition: SRLem.h:44
caf::SRLem::meanqfracsigexp
float meanqfracsigexp
Fraction matched in signal matches "exp".
Definition: SRLem.h:40
caf::SRLem::avgsigE
float avgsigE
Fraction of signal matches, sigmoid weighting.
Definition: SRLem.h:60
caf::SRLem::chisig
float chisig
chi-square of a fit assuming best match is signal
Definition: SRLem.h:53
caf::SRLem::meanysigexp
float meanysigexp
Hadronic y of matches that are signal, "exp".
Definition: SRLem.h:34
caf::SRLem::enrichfracexp
float enrichfracexp
Fraction of matches that are "enriched", "exp".
Definition: SRLem.h:48
caf::SRLem::avgexpE
float avgexpE
Fraction of signal matches, "exp" weighting.
Definition: SRLem.h:59
lem::PIDDetails::vars
PIDExtraVars vars
Definition: PIDDetails.h:32
caf::SRLem::meanqfracexp
float meanqfracexp
Fraction of charge matched, "exp".
Definition: SRLem.h:38
caf::SRLem::fannidenrich
float fannidenrich
ANN output, including enrichfracexp.
Definition: SRLem.h:51
caf::SRLem::pidexpfit
float pidexpfit
Fraction of matches that are signal, "expfit".
Definition: SRLem.h:28
caf::SRLem::meany
float meany
Mean hadronic y of matches.
Definition: SRLem.h:30
caf::SRLem::tmvamlp
float tmvamlp
TMVA ANN ("multi layer perceptron")
Definition: SRLem.h:66
caf::SRLem::meanqfrac
float meanqfrac
Mean fraction of charge matched.
Definition: SRLem.h:36
caf::SRLem::avginvE
float avginvE
Fraction of signal matches, 1/E weighting.
Definition: SRLem.h:58
caf::SRLem::chisigexp
float chisigexp
chisq of a fit assuming best match is signal, "exp"
Definition: SRLem.h:55
caf::SRLem::pidfit
float pidfit
Fraction of matches that are signal, "fit".
Definition: SRLem.h:26
caf::SRLem::chibkg
float chibkg
chisq of a fit assuming best match is background
Definition: SRLem.h:54
caf::SRLem::dectree
float dectree
Decision tree PID.
Definition: SRLem.h:68
caf::SRLem::tmvabdtg
float tmvabdtg
TMVA boosted decision tree, gradient boost.
Definition: SRLem.h:62
void caf::FillLIDEnergyVars ( const slid::ShowerLID &  slid,
caf::SRFuzzyKProng &  png 
)

Definition at line 272 of file FillPIDs.cxx.

References caf::SRSLidEnergy::depE, caf::SRSLidEnergy::E, caf::SRSLidEnergy::hadE, slid::ShowerLID::HadronicEnergy(), caf::SRShowerLID::lidE, slid::ShowerLID::NueEnergy(), slid::ShowerLID::ShowerDepEnergy(), slid::ShowerLID::ShowerEnergy(), caf::SRSLidEnergy::shwE, and caf::SRFuzzyKProng::shwlid.

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

274  {
275  png.shwlid.lidE.E = slid.NueEnergy();
276  png.shwlid.lidE.depE = slid.ShowerDepEnergy();
277  png.shwlid.lidE.shwE = slid.ShowerEnergy();
278  png.shwlid.lidE.hadE = slid.HadronicEnergy();
279  }
caf::SRFuzzyKProng::shwlid
SRShowerLID shwlid
Shower information.
Definition: SRFuzzyKProng.h:18
caf::SRSLidEnergy::shwE
float shwE
Energy of shower [GeV].
Definition: SRSLidEnergy.h:21
caf::SRSLidEnergy::E
float E
Energy [GeV].
Definition: SRSLidEnergy.h:18
caf::SRSLidEnergy::hadE
float hadE
Hadronic energy [GeV].
Definition: SRSLidEnergy.h:22
slid::ShowerLID::HadronicEnergy
float HadronicEnergy() const
Definition: ShowerLID.h:159
caf::SRSLidEnergy::depE
float depE
Total energy deposited in shower [GeV].
Definition: SRSLidEnergy.h:20
slid::ShowerLID::NueEnergy
float NueEnergy() const
Definition: ShowerLID.h:158
slid::ShowerLID::ShowerDepEnergy
float ShowerDepEnergy() const
Definition: ShowerLID.h:156
slid::ShowerLID::ShowerEnergy
float ShowerEnergy() const
Definition: ShowerLID.h:157
caf::SRShowerLID::lidE
SRSLidEnergy lidE
Definition: SRShowerLID.h:34
void caf::FillMRCCParentInfo ( const murem::MRCCParent &  parent,
caf::SRMRCCParent &  srparent,
const art::Event &  evt 
)

Definition at line 23 of file FillParentInfo.cxx.

References cosrejlabel, cvnlabel, caf::SRMRCCParent::eff, make_root_from_grid_output::evt, FillNumuInfo(), FillSliceInfo(), art::DataViewImpl::getByLabel(), numuelabel, murem::MRCCParent::ParentEff(), murem::MRCCParent::ParentIndex(), murem::MRCCParent::ParentPur(), caf::SRMRCCParent::pur, caf::SRMRCCParent::slcidx, slicelabel, trklabelbpf, trklabelkalman, and whattracker.

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

26  {
27 
28  const int parentSliceIndex = parent.ParentIndex();
29  srparent.slcidx = parentSliceIndex;
30  srparent.eff = parent.ParentEff();
31  srparent.pur = parent.ParentPur();
32 
33  art::Handle<std::vector<rb::Cluster> > sliceHandle;
34  evt.getByLabel(slicelabel, sliceHandle);
35 
36  art::Ptr<rb::Cluster> parentSlice(sliceHandle,
37  parentSliceIndex);
38 
39  // find numu objects associated with parentslice
40  art::FindOneP<cvn::Result> cvnassn({parentSlice},
41  evt,
42  cvnlabel);
43  std::string loadtrack = "";
44  if ( whattracker == "kalman" ) loadtrack = trklabelkalman;
45  if ( whattracker == "bpf" ) loadtrack = trklabelbpf;
46  art::FindManyP<rb::Track> trackassn({parentSlice},
47  evt,
48  loadtrack);
49  art::FindOneP<numue::NumuE> numueassn({parentSlice},
50  evt,
51  numuelabel);
52  art::FindOneP<cosrej::CosRejObj> cosrejassn({parentSlice},
53  evt,
54  cosrejlabel);
55  FillSliceInfo( parentSlice, srparent);
56  if(cvnassn.isValid() &&
57  trackassn.isValid() &&
58  cosrejassn.isValid() &&
59  numueassn.isValid()){
60  FillNumuInfo( cvnassn.at(0),
61  trackassn.at(0),
62  numueassn.at(0),
63  cosrejassn.at(0),
64  evt,
65  srparent,
66  parentSlice );
67  }
68  }// end of FillParentInfo
caf::slicelabel
const std::string slicelabel
Definition: FillParentInfo.h:19
caf::cvnlabel
const std::string cvnlabel
Definition: FillParentInfo.h:25
caf::SRMRCCParent::pur
float pur
Purity with which this event matched to the parent slice.
Definition: SRMRCCParent.h:30
make_root_from_grid_output.evt
evt
Definition: make_root_from_grid_output.py:82
murem::MRCCParent::ParentEff
float ParentEff() const
Returns the efficiency of MRCC slice wrt the parent.
Definition: MRCCParent.h:41
caf::FillNumuInfo
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)
Definition: FillParentInfo.cxx:73
caf::trklabelkalman
const std::string trklabelkalman
Definition: FillParentInfo.h:20
caf::FillSliceInfo
void FillSliceInfo(const art::Ptr< rb::Cluster > &slice, SRMRCCParent &parent)
Definition: FillParentInfo.cxx:190
caf::trklabelbpf
const std::string trklabelbpf
Definition: FillParentInfo.h:21
caf::SRMRCCParent::eff
float eff
Efficiency with which this event matched to the parent slice.
Definition: SRMRCCParent.h:28
caf::whattracker
const std::string whattracker
Definition: FillParentInfo.h:18
xml_schema::string
::xsd::cxx::tree::string< char, simple_type > string
Definition: Database.h:154
art::DataViewImpl::getByLabel
bool getByLabel(std::string const &label, std::string const &productInstanceName, Handle< PROD > &result) const
Definition: DataViewImpl.h:344
caf::numuelabel
const std::string numuelabel
Definition: FillParentInfo.h:22
murem::MRCCParent::ParentIndex
int ParentIndex() const
Returns the index of the parent slice.
Definition: MRCCParent.h:35
caf::SRMRCCParent::slcidx
int slcidx
Index of the parent slice.
Definition: SRMRCCParent.h:25
murem::MRCCParent::ParentPur
float ParentPur() const
Returns the purity of MRCC slice wrt the parent.
Definition: MRCCParent.h:47
caf::cosrejlabel
const std::string cosrejlabel
Definition: FillParentInfo.h:23
void caf::FillMuIdVars ( const remid::ReMId &  remid,
SRKalmanTrack &  srTrk 
)

Definition at line 36 of file FillPIDs.cxx.

References remid::ReMId::AvgDedx(), remid::ReMId::Contained(), caf::SRKalmanTrack::dedxllh, remid::ReMId::DedxSeparation(), caf::SRKalmanTrack::measfrac, remid::ReMId::MeasurementFraction(), caf::SRKalmanTrack::remavededx, caf::SRKalmanTrack::remcont, caf::SRKalmanTrack::rempid, caf::SRKalmanTrack::scatllh, remid::ReMId::ScatSeparation(), and rb::PID::Value().

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

37  {
38  srTrk.remavededx = remid.AvgDedx();
39  srTrk.rempid = remid.Value();
40  srTrk.scatllh = remid.ScatSeparation();
41  srTrk.dedxllh = remid.DedxSeparation();
42  srTrk.measfrac = remid.MeasurementFraction();
43  srTrk.remcont = remid.Contained();
44  }
remid::ReMId::ScatSeparation
double ScatSeparation() const
Return the scattering separation variable used as an input to the kNN that determines the pid value...
Definition: ReMId.cxx:206
remid::ReMId::AvgDedx
double AvgDedx() const
Definition: ReMId.cxx:68
rb::PID::Value
double Value() const
Definition: PID.h:22
remid::ReMId::Contained
bool Contained() const
Definition: ReMId.cxx:243
remid::ReMId::MeasurementFraction
double MeasurementFraction() const
Return the measurement fraction variable used as an input to the kNN that determines the pid value...
Definition: ReMId.cxx:213
remid::ReMId::DedxSeparation
double DedxSeparation() const
Return the dE/dx separation variable used as an input to the kNN that determines the pid value...
Definition: ReMId.cxx:199
void caf::FillNCCosRejVars ( const ncid::NCCosRej &  nccosrej,
caf::SRNCCosRej &  srnccosrej 
)

Definition at line 487 of file FillPIDs.cxx.

References caf::SRNCCosRej::cospidbdt, caf::SRNCCosRej::cospidbdtg, ncid::NCCosRej::CosPIDDT(), ncid::NCCosRej::CosPIDDTG(), caf::SRNCCosRej::pngptp, and ncid::NCCosRej::ProngTransMom().

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

489  {
490  srnccosrej.pngptp = nccosrej.ProngTransMom();
491  srnccosrej.cospidbdt = nccosrej.CosPIDDT();
492  srnccosrej.cospidbdtg = nccosrej.CosPIDDTG();
493  }
caf::SRNCCosRej::cospidbdtg
float cospidbdtg
NC cosrej PID. A Keras using 15 variables.
Definition: SRNCCosRej.h:23
ncid::NCCosRej::CosPIDDTG
double CosPIDDTG() const
Gradient BDT trained, based on simple containment variables.
Definition: NCCosRej.h:26
caf::SRNCCosRej::cospidbdt
float cospidbdt
NC Cosrej PID. A real adaptive BDT using 13 variables.
Definition: SRNCCosRej.h:20
caf::SRNCCosRej::pngptp
float pngptp
Event transverse momentum fraction, based on prongs. See ncid::NCCosRej::ProngTransMom.
Definition: SRNCCosRej.h:17
ncid::NCCosRej::CosPIDDT
double CosPIDDT() const
Real Adaptive BDT trained, based on simple containment variables.
Definition: NCCosRej.h:23
ncid::NCCosRej::ProngTransMom
double ProngTransMom() const
Transverse component of the energy-weighted average of prong directions.
Definition: NCCosRej.h:20
void caf::FillNCPi0BkgRejVars ( const ncpi0::NCPi0BkgRej &  ncpi0bkgrej,
caf::SRNCPi0BkgRej &  srncpi0bkgrej 
)

Definition at line 496 of file FillPIDs.cxx.

References caf::SRNCPi0BkgRej::ncpidbdtg, ncpi0::NCPi0BkgRej::NCPIDBDTG(), caf::SRNCPi0BkgRej::ncpidbdtglt, and ncpi0::NCPi0BkgRej::NCPIDBDTGLT().

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

498  {
499  srncpi0bkgrej.ncpidbdtg = ncpi0bkgrej.NCPIDBDTG();
500  srncpi0bkgrej.ncpidbdtglt = ncpi0bkgrej.NCPIDBDTGLT();
501  }
ncpi0::NCPi0BkgRej::NCPIDBDTGLT
double NCPIDBDTGLT() const
Definition: NCPi0BkgRej.h:22
ncpi0::NCPi0BkgRej::NCPIDBDTG
double NCPIDBDTG() const
Gradient BDT trained, based on prong/shower variables.
Definition: NCPi0BkgRej.h:19
void caf::FillNueCosRejVars ( const cosrej::NueCosRej &  nuecosrej,
caf::SRNueCosRej &  srnuecosrej 
)

Definition at line 424 of file FillPIDs.cxx.

References cosrej::NueCosRej::CoreBDT(), cosrej::NueCosRej::CosAngleToNextProng(), caf::SRNueCosRej::cosdang, caf::SRNueCosRej::cospidcontain, cosrej::NueCosRej::CosPIDContain(), caf::SRNueCosRej::cospidcontainxy, cosrej::NueCosRej::CosPIDContainXY(), caf::SRNueCosRej::cospidcorebdt, caf::SRNueCosRej::cospidlight, cosrej::NueCosRej::CosPIDLight(), caf::SRNueCosRej::cospidperibdt, caf::SRNueCosRej::distallpngback, caf::SRNueCosRej::distallpngbottom, caf::SRNueCosRej::distallpngeast, caf::SRNueCosRej::distallpngfront, caf::SRNueCosRej::distallpngtop, caf::SRNueCosRej::distallpngwest, caf::SRNueCosRej::hitsperplane, cosrej::NueCosRej::HitsPerPlane(), caf::SRNueCosRej::hitsperplaneasymm, cosrej::NueCosRej::HitsPerPlaneAsymm(), caf::SRNueCosRej::hitsperplaneasymmslice, cosrej::NueCosRej::HitsPerPlaneAsymmSlice(), min(), caf::SRNueCosRej::muanglediffbydist, cosrej::NueCosRej::MuAngleDiffByDist(), caf::SRNueCosRej::muanglediffbytime, cosrej::NueCosRej::MuAngleDiffByTime(), caf::SRNueCosRej::muclosestapproachbydist, cosrej::NueCosRej::MuClosestApproachByDist(), caf::SRNueCosRej::muclosestapproachbytime, cosrej::NueCosRej::MuClosestApproachByTime(), caf::SRNueCosRej::musliceidxbydist, cosrej::NueCosRej::MuSliceIdxByDist(), caf::SRNueCosRej::musliceidxbytime, cosrej::NueCosRej::MuSliceIdxByTime(), caf::SRNueCosRej::mutimediffbydist, cosrej::NueCosRej::MuTimeDiffByDist(), caf::SRNueCosRej::mutimediffbytime, cosrej::NueCosRej::MuTimeDiffByTime(), cosrej::NueCosRej::ParticleShowerTransMom(), caf::SRNueCosRej::partptp, cosrej::NueCosRej::PeriBDT(), cosrej::NueCosRej::PhotonShowerMomX(), cosrej::NueCosRej::PhotonShowerMomY(), cosrej::NueCosRej::PhotonShowerTransMom(), caf::SRNueCosRej::photptp, caf::SRNueCosRej::photpxp, caf::SRNueCosRej::photpyp, caf::SRNueCosRej::pngptp, cosrej::NueCosRej::ProngDistToVtx(), caf::SRNueCosRej::prongmaxx, cosrej::NueCosRej::ProngMaxX(), caf::SRNueCosRej::prongmaxy, cosrej::NueCosRej::ProngMaxY(), caf::SRNueCosRej::prongmaxz, cosrej::NueCosRej::ProngMaxZ(), caf::SRNueCosRej::prongminx, cosrej::NueCosRej::ProngMinX(), caf::SRNueCosRej::prongminy, cosrej::NueCosRej::ProngMinY(), caf::SRNueCosRej::prongminz, cosrej::NueCosRej::ProngMinZ(), cosrej::NueCosRej::ProngTransMom(), caf::SRNueCosRej::sparsenessasymm, cosrej::NueCosRej::SparsenessAsymm(), caf::SRNueCosRej::sparsenessasymmslice, cosrej::NueCosRej::SparsenessAsymmSlice(), caf::SRNueCosRej::startback, caf::SRNueCosRej::startbottom, cosrej::NueCosRej::StartDistToBack(), cosrej::NueCosRej::StartDistToBottom(), cosrej::NueCosRej::StartDistToEast(), cosrej::NueCosRej::StartDistToFront(), cosrej::NueCosRej::StartDistToTop(), cosrej::NueCosRej::StartDistToWest(), caf::SRNueCosRej::starteast, caf::SRNueCosRej::startfront, cosrej::NueCosRej::StartMinDistToBack(), cosrej::NueCosRej::StartMinDistToBottom(), cosrej::NueCosRej::StartMinDistToEast(), cosrej::NueCosRej::StartMinDistToFront(), cosrej::NueCosRej::StartMinDistToTop(), cosrej::NueCosRej::StartMinDistToWest(), caf::SRNueCosRej::starttop, caf::SRNueCosRej::startwest, caf::SRNueCosRej::stopback, caf::SRNueCosRej::stopbottom, cosrej::NueCosRej::StopDistToBack(), cosrej::NueCosRej::StopDistToBottom(), cosrej::NueCosRej::StopDistToEast(), cosrej::NueCosRej::StopDistToFront(), cosrej::NueCosRej::StopDistToTop(), cosrej::NueCosRej::StopDistToWest(), caf::SRNueCosRej::stopeast, caf::SRNueCosRej::stopfront, cosrej::NueCosRej::StopMinDistToBack(), cosrej::NueCosRej::StopMinDistToBottom(), cosrej::NueCosRej::StopMinDistToEast(), cosrej::NueCosRej::StopMinDistToFront(), cosrej::NueCosRej::StopMinDistToTop(), cosrej::NueCosRej::StopMinDistToWest(), caf::SRNueCosRej::stoptop, caf::SRNueCosRej::stopwest, and caf::SRNueCosRej::vtxdoca.

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

426  {
427  srnuecosrej.hitsperplane = nuecosrej.HitsPerPlane();
428  srnuecosrej.pngptp = nuecosrej.ProngTransMom();
429  srnuecosrej.partptp = nuecosrej.ParticleShowerTransMom();
430  srnuecosrej.photptp = nuecosrej.PhotonShowerTransMom();
431  srnuecosrej.photpxp = nuecosrej.PhotonShowerMomX();
432  srnuecosrej.photpyp = nuecosrej.PhotonShowerMomY();
433  srnuecosrej.starttop = nuecosrej.StartDistToTop();
434  srnuecosrej.startbottom = nuecosrej.StartDistToBottom();
435  srnuecosrej.startfront = nuecosrej.StartDistToFront();
436  srnuecosrej.startback = nuecosrej.StartDistToBack();
437  srnuecosrej.startwest = nuecosrej.StartDistToWest();
438  srnuecosrej.starteast = nuecosrej.StartDistToEast();
439  srnuecosrej.stoptop = nuecosrej.StopDistToTop();
440  srnuecosrej.stopbottom = nuecosrej.StopDistToBottom();
441  srnuecosrej.stopfront = nuecosrej.StopDistToFront();
442  srnuecosrej.stopback = nuecosrej.StopDistToBack();
443  srnuecosrej.stopwest = nuecosrej.StopDistToWest();
444  srnuecosrej.stopeast = nuecosrej.StopDistToEast();
445  srnuecosrej.cospidcontain= nuecosrej.CosPIDContain();
446  srnuecosrej.cospidcontainxy= nuecosrej.CosPIDContainXY();
447  srnuecosrej.cospidlight = nuecosrej.CosPIDLight();
448  srnuecosrej.cospidperibdt = nuecosrej.PeriBDT();
449  srnuecosrej.cospidcorebdt = nuecosrej.CoreBDT();
450  srnuecosrej.cosdang = nuecosrej.CosAngleToNextProng();
451  srnuecosrej.vtxdoca = nuecosrej.ProngDistToVtx();
452  srnuecosrej.prongmaxx = nuecosrej.ProngMaxX();
453  srnuecosrej.prongmaxy = nuecosrej.ProngMaxY();
454  srnuecosrej.prongmaxz = nuecosrej.ProngMaxZ();
455  srnuecosrej.prongminx = nuecosrej.ProngMinX();
456  srnuecosrej.prongminy = nuecosrej.ProngMinY();
457  srnuecosrej.prongminz = nuecosrej.ProngMinZ();
458  srnuecosrej.distallpngtop = std::min(nuecosrej.StartMinDistToTop(),
459  nuecosrej.StopMinDistToTop());
460  srnuecosrej.distallpngbottom = std::min(nuecosrej.StartMinDistToBottom(),
461  nuecosrej.StopMinDistToBottom());
462  srnuecosrej.distallpngfront = std::min(nuecosrej.StartMinDistToFront(),
463  nuecosrej.StopMinDistToFront());
464  srnuecosrej.distallpngback = std::min(nuecosrej.StartMinDistToBack(),
465  nuecosrej.StopMinDistToBack());
466  srnuecosrej.distallpngwest = std::min(nuecosrej.StartMinDistToWest(),
467  nuecosrej.StopMinDistToWest());
468  srnuecosrej.distallpngeast = std::min(nuecosrej.StartMinDistToEast(),
469  nuecosrej.StopMinDistToEast());
470 
471  srnuecosrej.sparsenessasymm = nuecosrej.SparsenessAsymm();
472  srnuecosrej.hitsperplaneasymm = nuecosrej.HitsPerPlaneAsymm();
473  srnuecosrej.sparsenessasymmslice = nuecosrej.SparsenessAsymmSlice();
474  srnuecosrej.hitsperplaneasymmslice = nuecosrej.HitsPerPlaneAsymmSlice();
475  srnuecosrej.musliceidxbydist = nuecosrej.MuSliceIdxByDist();
476  srnuecosrej.muanglediffbydist = nuecosrej.MuAngleDiffByDist();
477  srnuecosrej.mutimediffbydist = nuecosrej.MuTimeDiffByDist();
478  srnuecosrej.muclosestapproachbydist = nuecosrej.MuClosestApproachByDist();
479  srnuecosrej.musliceidxbytime = nuecosrej.MuSliceIdxByTime();
480  srnuecosrej.muanglediffbytime = nuecosrej.MuAngleDiffByTime();
481  srnuecosrej.mutimediffbytime = nuecosrej.MuTimeDiffByTime();
482  srnuecosrej.muclosestapproachbytime = nuecosrej.MuClosestApproachByTime();
483 
484  }
cosrej::NueCosRej::CosAngleToNextProng
double CosAngleToNextProng() const
Cosine of Angle between the leading and the next leading prong.
Definition: NueCosRej.h:164
cosrej::NueCosRej::MuClosestApproachByDist
double MuClosestApproachByDist() const
Distance of closest approach between this slice and the possible muon parent (selected by closest app...
Definition: NueCosRej.h:224
cosrej::NueCosRej::ProngMaxY
double ProngMaxY() const
The maximum Y between the start and stop points of the leading prong in the slice.
Definition: NueCosRej.h:176
cosrej::NueCosRej::MuClosestApproachByTime
double MuClosestApproachByTime() const
Distance of closest approach between this slice and the possible muon parent (selected by difference ...
Definition: NueCosRej.h:241
cosrej::NueCosRej::StopMinDistToBottom
double StopMinDistToBottom() const
Minimum perpendicular distance of all prongs stop points to the bottom edge of the detector...
Definition: NueCosRej.h:124
caf::SRNueCosRej::distallpngbottom
float distallpngbottom
Definition: SRNueCosRej.h:95
caf::SRNueCosRej::mutimediffbydist
float mutimediffbydist
Definition: SRNueCosRej.h:193
caf::SRNueCosRej::muanglediffbytime
float muanglediffbytime
Definition: SRNueCosRej.h:206
cosrej::NueCosRej::MuAngleDiffByDist
double MuAngleDiffByDist() const
Minimum angle difference between this slice and the possible muon parent (selected by closest approac...
Definition: NueCosRej.h:220
cosrej::NueCosRej::MuAngleDiffByTime
double MuAngleDiffByTime() const
Minimum angle difference between this slice and the possible muon parent (selected by difference in t...
Definition: NueCosRej.h:237
cosrej::NueCosRej::StopDistToTop
double StopDistToTop() const
Perpendicular distance of stop point of the leading prong to the top edge of the detector.
Definition: NueCosRej.h:72
cosrej::NueCosRej::CoreBDT
double CoreBDT() const
Definition: NueCosRej.h:158
cosrej::NueCosRej::StopMinDistToWest
double StopMinDistToWest() const
Minimum perpendicular distance of all prongs stop points to the west edge of the detector (west is po...
Definition: NueCosRej.h:136
cosrej::NueCosRej::StartMinDistToTop
double StartMinDistToTop() const
Minimum perpendicular distance of all prongs start points.
Definition: NueCosRej.h:96
caf::SRNueCosRej::muclosestapproachbydist
float muclosestapproachbydist
Definition: SRNueCosRej.h:197
cosrej::NueCosRej::PeriBDT
double PeriBDT() const
Definition: NueCosRej.h:157
cosrej::NueCosRej::StopDistToWest
double StopDistToWest() const
Perpendicular distance of stop point of the leading prong to the west edge of the detector (west is p...
Definition: NueCosRej.h:88
cosrej::NueCosRej::StopDistToFront
double StopDistToFront() const
Perpendicular distance of stop point of the leading prong to the front edge of the detector...
Definition: NueCosRej.h:80
cosrej::NueCosRej::HitsPerPlaneAsymmSlice
double HitsPerPlaneAsymmSlice() const
Asymmetry in the hitsperplane at the start and end of the slice.
Definition: NueCosRej.h:207
cosrej::NueCosRej::PhotonShowerTransMom
double PhotonShowerTransMom() const
Transverse component of the event momentum. Particle masses in the momentum calculation are assumed t...
Definition: NueCosRej.h:36
caf::SRNueCosRej::muanglediffbydist
float muanglediffbydist
Definition: SRNueCosRej.h:189