Public Types | Public Member Functions | Static Public Member Functions | Protected Member Functions | Private Member Functions | Private Attributes | List of all members
skim::CosmicBeamComparison Class Reference
Inheritance diagram for skim::CosmicBeamComparison:
art::EDAnalyzer art::EventObserverBase art::Consumer art::EngineCreator

Public Types

using WorkerType = WorkerT< EDAnalyzer >
 
using ModuleType = EDAnalyzer
 

Public Member Functions

 CosmicBeamComparison (fhicl::ParameterSet const &p)
 
virtual ~CosmicBeamComparison ()
 
 CosmicBeamComparison (CosmicBeamComparison const &)=delete
 
 CosmicBeamComparison (CosmicBeamComparison &&)=delete
 
CosmicBeamComparisonoperator= (CosmicBeamComparison const &)=delete
 
CosmicBeamComparisonoperator= (CosmicBeamComparison &&)=delete
 
void analyze (art::Event const &e) override
 
void reconfigure (fhicl::ParameterSet const &p)
 
void beginJob () override
 
void beginRun (art::Run const &r) override
 
void beginSubRun (art::SubRun const &sr) override
 
std::string workerType () const
 
bool modifiesEvent () const
 
void registerProducts (MasterProductRegistry &, ProductDescriptions &, ModuleDescription const &)
 
std::string const & processName () const
 
bool wantAllEvents () const
 
bool wantEvent (Event const &e)
 
fhicl::ParameterSetID selectorConfig () const
 
art::Handle< art::TriggerResultsgetTriggerResults (Event const &e) const
 
template<typename T , BranchType = InEvent>
ProductToken< Tconsumes (InputTag const &)
 
template<typename T , art::BranchType BT>
art::ProductToken< Tconsumes (InputTag const &it)
 
template<typename T , BranchType = InEvent>
void consumesMany ()
 
template<typename Element , BranchType = InEvent>
ViewToken< Element > consumesView (InputTag const &)
 
template<typename T , art::BranchType BT>
art::ViewToken< TconsumesView (InputTag const &it)
 
template<typename T , BranchType = InEvent>
ProductToken< TmayConsume (InputTag const &)
 
template<typename T , art::BranchType BT>
art::ProductToken< TmayConsume (InputTag const &it)
 
template<typename T , BranchType = InEvent>
void mayConsumeMany ()
 
template<typename Element , BranchType = InEvent>
ViewToken< Element > mayConsumeView (InputTag const &)
 
template<typename T , art::BranchType BT>
art::ViewToken< TmayConsumeView (InputTag const &it)
 
base_engine_tcreateEngine (seed_t seed)
 
base_engine_tcreateEngine (seed_t seed, std::string const &kind_of_engine_to_make)
 
base_engine_tcreateEngine (seed_t seed, std::string const &kind_of_engine_to_make, label_t const &engine_label)
 
seed_t get_seed_value (fhicl::ParameterSet const &pset, char const key[]="seed", seed_t const implicit_seed=-1)
 

Static Public Member Functions

static cet::exempt_ptr< Consumernon_module_context ()
 

Protected Member Functions

CurrentProcessingContext const * currentContext () const
 
detail::CachedProducts & cachedProducts ()
 
void validateConsumedProduct (BranchType const bt, ProductInfo const &pi)
 
void prepareForJob (fhicl::ParameterSet const &pset)
 
void showMissingConsumes () const
 

Private Member Functions

void FillNueHistograms (art::Event const &e)
 
void FillNumuHistograms (art::Event const &e)
 

Private Attributes

std::string fSkimmedLabel
 label of module creating skimmed file More...
 
std::string fCosmicExposureLabel
 label of module creating cosmic exposure More...
 
std::string fPOTSumLabel
 label of module creating pot sum More...
 
std::string fNumuInstance
 label for numu products More...
 
std::string fNueInstance
 label for nue products More...
 
bool fCompareNumuCuts
 bool for comparing the numu cuts More...
 
bool fCompareNueCuts
 bool for comparing the nue cuts More...
 
novadaq::cnv::DetId fDetId
 the id of the detector we are looking at More...
 
art::ServiceHandle< ds::DetectorServicefDetService
 detector service More...
 
art::ServiceHandle< geo::GeometryfGeom
 
TTree * fNumuTree
 Tree to hold values of cut parameters for selected events. More...
 
TTree * fNueTree
 Tree to hold values of cut parameters for selected events. More...
 
TTree * fCosmicTree
 Tree filled directly from cosrej etc. objects (bypasses parameters struct) More...
 
TTree * fHitTree
 Tree to hold information about every hit in a kalman track. More...
 
TTree * fCosmicHitTree
 Tree to hold information about every hit in a kalman track. More...
 
NumuCutParameters fNumuCutParams
 from the ParametersNumu object More...
 
NueCutParameters fNueCutParams
 from the ParametersNumu object More...
 
NumuCutParameters fNumuSRParams
 from the ParametersNumu object More...
 
NueCutParameters fNueSRParams
 from the ParametersNumu object More...
 
EventInfo fEventInfo
 keep track of run and subrun, and event More...
 
TH1F * fPOTByRun
 
TH1F * fGoodPOTByRun
 
TH1F * fCosmicExposureByRun
 
TH1F * fNeutrinoE
 
TH1F * fReMId
 
TH1F * fSliceHits
 
TH1F * fSliceContigPlanes
 
TH1F * fCellsFromEdge
 
TH1F * fQePId1Track
 
TH1F * fQePId2Track
 
TH1F * fFDKalmanFwdCell
 
TH1F * fFDKalmanBakCell
 
TH1F * fFDCosFwdCell
 
TH1F * fFDCosBakCell
 
TH1F * fFDPlanesToFront
 
TH1F * fFDPlanesToBack
 
TH1F * fFDKalmanAngle
 
TH1F * fFDKalmanDirX
 
TH1F * fFDKalmanDirY
 
TH1F * fFDKalmanDirZ
 
TH1F * fFDKalSpeed
 
TH1F * fFDCosSpeed
 
TH1F * fFDDeltaT
 
TH1F * fFDDeltaX
 
TH1F * fFDDeltaY
 
TH1F * fFDDeltaZ
 
TH1F * fFDSpeed
 
TH1F * fFDTrackLength
 
TH2F * fFDDTvsDY
 
TH2F * fFDCosSpeedVsScore
 
TH1F * fFDContPID
 
TH1F * fFDSliceHits
 
TH1F * fFDDiffLastFirstLivePlanes
 
TH1F * fNDMaxPlaneHit
 
TH1F * fNDMinPlaneHit
 
TH1F * fNDTrackStartZ
 
TH1F * fNDTrackStopZ
 
TH1F * fNDKalmanFwdCell
 
TH1F * fNDKalmanBakCell
 
TH1F * fNDHadCalEnergySum
 
TH1F * fFDSliceHitsPerPlane
 
TH1F * fFDSliceHitsXView
 
TH1F * fFDSliceHitsYView
 
TH1F * fFDSliceHitsViewAsym
 
TH1F * fFDSliceCosAngleShw
 
TH1F * fFDSliceDistShwVtx
 
TH1F * fFDSliceFracShwHits
 
TH1F * fFDVtxDistToEastWall
 
TH1F * fFDVtxDistToWestWall
 
TH1F * fFDVtxDistToTop
 
TH1F * fFDVtxDistToBottom
 
TH1F * fFDVtxDistToFront
 
TH1F * fFDVtxDistToBack
 
TH1F * fFDPt
 
TH1F * fCalorimetricEnergy
 
TH1F * fLEM
 
TH1F * fCVN
 
TH1F * fLID
 
TH1F * fNDVtxX
 
TH1F * fNDVtxY
 
TH1F * fNDVtxZ
 
TH1F * fNDNShowers
 
TH1F * fNDShwVtxX
 
TH1F * fNDShwVtxY
 
TH1F * fNDShwVtxZ
 
TH1F * fNDShwEndX
 
TH1F * fNDShwEndY
 
TH1F * fNDShwEndZ
 
TH1F * fNDShowerCalE
 
TH1F * fNDShowerCalEfraction
 
TH1F * fNDShowerDistToVtx
 
TH1F * fNDShowerLength
 
TH1F * fNDCalorimetricE
 
TH1F * fNDCalEperHit
 
TH1F * fNDPlanesToFront
 
TH1F * fNDSliceHits
 
TH1F * fNDSliceNHitsPerPlane
 
TH1F * fNDSliceNoShowerCalE
 
float crNeutrinoE = 0
 
double crReMIdVal
 
float crSliceHits
 
float crSliceXViewHits
 
float crSliceYViewHits
 
float crSliceContigPlanes
 
float crSliceMaxPlane
 
float crSliceMinPlane
 
float crSliceTotalADC
 
float crSliceTotalGeV
 
float crSliceCalorimetricEnergy
 
TVector3 crSliceMinXYZ
 
TVector3 crSliceMaxXYZ
 
float crSliceMinTNS
 
float crSliceMaxTNS
 
float crCellsFromEdge
 
int crSliceClosestEdge
 
float crPlanesToFront
 
float crPlanesToBack
 
int crSliceNHitsNearEdge0
 
int crSliceNHitsNearEdge1
 
int crSliceNHitsNearEdge2
 
int crSliceNHitsNearEdge3
 
int crSliceNHitsNearEdge4
 
int crSliceNHitsNearEdge5
 
int crSliceNHitsNearEdge6
 
int crSliceNHitsNearEdge7
 
int crSliceNHitsNearEdge8
 
int crSliceNHitsNearEdge9
 
float crKalAngle
 
float crCosAngle
 
float crKalFwdCell
 
float crKalBakCell
 
float crCosFwdCell
 
float crCosBakCell
 
float crKalFwdDist
 
float crKalBakDist
 
float crCosFwdDist
 
float crCosBakDist
 
float crConCosPID
 
TVector3 crTrackStart
 
TVector3 crTrackStop
 
int crTrackHits
 
int crTrackNHitsNearEdge0
 
int crTrackNHitsNearEdge1
 
int crTrackNHitsNearEdge2
 
int crTrackNHitsNearEdge3
 
int crTrackNHitsNearEdge4
 
int crTrackNHitsNearEdge5
 
int crTrackNHitsNearEdge6
 
int crTrackNHitsNearEdge7
 
int crTrackNHitsNearEdge8
 
int crTrackNHitsNearEdge9
 
int crTrackClosestEdge
 
int crTrackCellsFromEdge
 
int crTrackXViewHits
 
int crTrackYViewHits
 
float crTrackTotalADC
 
float crTrackTotalGeV
 
float crTrackCalorimetricEnergy
 
float crTrackMinTNS
 
float crTrackMaxTNS
 
TVector3 crTrackMinXYZ
 
TVector3 crTrackMaxXYZ
 
TVector3 crTrackEarliestHitPos
 
TVector3 crTrackLatestHitPos
 
float crTrackMinPlane
 
float crTrackMaxPlane
 
float crTrackTotalLength
 
TVector3 crKalDir
 
float crTrackStartTime
 
float crTrackStopTime
 
TVector3 crCosTrackStart
 
TVector3 crCosTrackStop
 
TVector3 crCosTrackDir
 
int crCosTrackHits
 
int crCosTrackClosestEdge
 
int crCosTrackCellsFromEdge
 
int crCosTrackNHitsNearEdge0
 
int crCosTrackNHitsNearEdge1
 
int crCosTrackNHitsNearEdge2
 
int crCosTrackNHitsNearEdge3
 
int crCosTrackNHitsNearEdge4
 
int crCosTrackNHitsNearEdge5
 
int crCosTrackNHitsNearEdge6
 
int crCosTrackNHitsNearEdge7
 
int crCosTrackNHitsNearEdge8
 
int crCosTrackNHitsNearEdge9
 
int crCosTrackXViewHits
 
int crCosTrackYViewHits
 
float crCosTrackTotalADC
 
float crCosTrackTotalGeV
 
float crCosTrackCalorimetricEnergy
 
float crCosTrackMinTNS
 
float crCosTrackMaxTNS
 
TVector3 crCosTrackMinXYZ
 
TVector3 crCosTrackMaxXYZ
 
TVector3 crCosTrackEarliestHitPos
 
TVector3 crCosTrackLatestHitPos
 
float crCosTrackMinPlane
 
float crCosTrackMaxPlane
 
float crCosTrackTotalLength
 
float crCosTrackStartTime
 
float crCosTrackStopTime
 
float crKalSpeed
 
float crCosSpeed
 
float crCosDirScore
 
float crKalDirScore
 
float crCosScoreDiff
 
float crKalScoreDiff
 
float crKalSlope
 
float crCosSlope
 
float crKalChiSq
 
float crCosChiSq
 
float crKalChiDiff
 
float crCosChiDiff
 
float crNCloseLiveDiblocks
 
int crNTracks3D
 
float crScatt
 
float crQePIDVal
 
int crQePIDNTracks
 
float crXmin2
 
float crYmin2
 
float crZmin2
 
float crXmax2
 
float crYmax2
 
float crZmax2
 
float crscatllh
 
float crdedxllh
 
float crmeasfrac
 
int crTruthNNu
 
int crTruthNuNHitsSlc
 
int crTruthNuNHitsTot
 
bool crTruthNuIsCC
 
int crTruthNuPDG
 
int crTruthNuPDGOrig
 
TVector3 crTruthNuVertex
 
float crTruthNuTime
 
int crRun
 
int crSubrun
 
int crEvent
 
int htIndex
 
float htUncalibTNS
 
TVector3 htUncalibXYZ
 
double htUncalibW
 
float htCalibTNS
 
TVector3 htCalibXYZ
 
float htCalibPE
 
double htDistFromStart
 
double htDistFromStop
 
TVector3 htDir
 
TVector3 htTrackXYZ
 
int htPlane
 
int htCell
 
int htTrackIndex = 0
 
int chtIndex
 
float chtUncalibTNS
 
TVector3 chtUncalibXYZ
 
double chtUncalibW
 
double chtDistFromStart
 
double chtDistFromStop
 
TVector3 chtDir
 
TVector3 chtTrackXYZ
 
int chtPlane
 
int chtCell
 

Detailed Description

Definition at line 64 of file CosmicBeamComparison_module.cc.

Member Typedef Documentation

Definition at line 39 of file EDAnalyzer.h.

Definition at line 38 of file EDAnalyzer.h.

Constructor & Destructor Documentation

skim::CosmicBeamComparison::CosmicBeamComparison ( fhicl::ParameterSet const &  p)
explicit

Definition at line 375 of file CosmicBeamComparison_module.cc.

References reconfigure().

376 : EDAnalyzer(p)
377 {
378  this->reconfigure(p);
379 
380  return;
381 }
const char * p
Definition: xmltok.h:285
EDAnalyzer(Table< Config > const &config)
Definition: EDAnalyzer.h:100
void reconfigure(fhicl::ParameterSet const &p)
skim::CosmicBeamComparison::~CosmicBeamComparison ( )
virtual

Definition at line 384 of file CosmicBeamComparison_module.cc.

385 {
386 }
skim::CosmicBeamComparison::CosmicBeamComparison ( CosmicBeamComparison const &  )
delete
skim::CosmicBeamComparison::CosmicBeamComparison ( CosmicBeamComparison &&  )
delete

Member Function Documentation

void skim::CosmicBeamComparison::analyze ( art::Event const &  e)
overridevirtual

Implements art::EDAnalyzer.

Definition at line 752 of file CosmicBeamComparison_module.cc.

References art::EventID::event(), EventInfo::event, fEventInfo, FillNueHistograms(), FillNumuHistograms(), fNueInstance, fNumuInstance, art::Event::id(), LOG_DEBUG, art::Event::run(), EventInfo::run, art::Event::subRun(), and EventInfo::subrun.

753 {
754  fEventInfo.run = e.run();
755  fEventInfo.subrun = e.subRun();
756  fEventInfo.event = e.id().event();
757 
758  LOG_DEBUG("CosmicBeamComparison")
759  << e.run() << " " << e.subRun() << " " << e.id().event();
760 
761  if( !fNumuInstance.empty() ) this->FillNumuHistograms(e);
762  if( !fNueInstance.empty() ) this->FillNueHistograms(e);
763 
764  return;
765 }
std::string fNumuInstance
label for numu products
std::string fNueInstance
label for nue products
void FillNueHistograms(art::Event const &e)
#define LOG_DEBUG(stream)
Definition: Messenger.h:149
void FillNumuHistograms(art::Event const &e)
EventInfo fEventInfo
keep track of run and subrun, and event
Float_t e
Definition: plot.C:35
void skim::CosmicBeamComparison::beginJob ( )
overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 389 of file CosmicBeamComparison_module.cc.

References chtCell, chtDir, chtDistFromStart, chtDistFromStop, chtIndex, chtPlane, chtTrackXYZ, chtUncalibTNS, chtUncalibW, chtUncalibXYZ, crCellsFromEdge, crConCosPID, crCosAngle, crCosBakCell, crCosBakDist, crCosChiDiff, crCosChiSq, crCosDirScore, crCosFwdCell, crCosFwdDist, crCosScoreDiff, crCosSlope, crCosSpeed, crCosTrackCalorimetricEnergy, crCosTrackCellsFromEdge, crCosTrackClosestEdge, crCosTrackEarliestHitPos, crCosTrackHits, crCosTrackLatestHitPos, crCosTrackMaxPlane, crCosTrackMaxTNS, crCosTrackMaxXYZ, crCosTrackMinPlane, crCosTrackMinTNS, crCosTrackMinXYZ, crCosTrackNHitsNearEdge0, crCosTrackNHitsNearEdge1, crCosTrackNHitsNearEdge2, crCosTrackNHitsNearEdge3, crCosTrackNHitsNearEdge4, crCosTrackNHitsNearEdge5, crCosTrackNHitsNearEdge6, crCosTrackNHitsNearEdge7, crCosTrackNHitsNearEdge8, crCosTrackNHitsNearEdge9, crCosTrackStart, crCosTrackStartTime, crCosTrackStop, crCosTrackStopTime, crCosTrackTotalADC, crCosTrackTotalGeV, crCosTrackTotalLength, crCosTrackXViewHits, crCosTrackYViewHits, crdedxllh, crEvent, crKalAngle, crKalBakCell, crKalBakDist, crKalChiDiff, crKalChiSq, crKalDir, crKalDirScore, crKalFwdCell, crKalFwdDist, crKalScoreDiff, crKalSlope, crKalSpeed, crmeasfrac, crNCloseLiveDiblocks, crNeutrinoE, crNTracks3D, crPlanesToBack, crPlanesToFront, crQePIDNTracks, crQePIDVal, crReMIdVal, crRun, crscatllh, crScatt, crSliceCalorimetricEnergy, crSliceClosestEdge, crSliceContigPlanes, crSliceHits, crSliceMaxPlane, crSliceMaxTNS, crSliceMaxXYZ, crSliceMinPlane, crSliceMinTNS, crSliceMinXYZ, crSliceNHitsNearEdge0, crSliceNHitsNearEdge1, crSliceNHitsNearEdge2, crSliceNHitsNearEdge3, crSliceNHitsNearEdge4, crSliceNHitsNearEdge5, crSliceNHitsNearEdge6, crSliceNHitsNearEdge7, crSliceNHitsNearEdge8, crSliceNHitsNearEdge9, crSliceTotalADC, crSliceTotalGeV, crSliceXViewHits, crSliceYViewHits, crSubrun, crTrackCalorimetricEnergy, crTrackCellsFromEdge, crTrackClosestEdge, crTrackEarliestHitPos, crTrackHits, crTrackLatestHitPos, crTrackMaxPlane, crTrackMaxTNS, crTrackMaxXYZ, crTrackMinPlane, crTrackMinTNS, crTrackMinXYZ, crTrackNHitsNearEdge0, crTrackNHitsNearEdge1, crTrackNHitsNearEdge2, crTrackNHitsNearEdge3, crTrackNHitsNearEdge4, crTrackNHitsNearEdge5, crTrackNHitsNearEdge6, crTrackNHitsNearEdge7, crTrackNHitsNearEdge8, crTrackNHitsNearEdge9, crTrackStart, crTrackStartTime, crTrackStop, crTrackStopTime, crTrackTotalADC, crTrackTotalGeV, crTrackTotalLength, crTrackXViewHits, crTrackYViewHits, crTruthNNu, crTruthNuIsCC, crTruthNuNHitsSlc, crTruthNuNHitsTot, crTruthNuPDG, crTruthNuPDGOrig, crTruthNuTime, crTruthNuVertex, crXmax2, crXmin2, crYmax2, crYmin2, crZmax2, crZmin2, fCalorimetricEnergy, fCellsFromEdge, fCosmicExposureByRun, fCosmicHitTree, fCosmicTree, fCVN, fEventInfo, fFDContPID, fFDCosBakCell, fFDCosFwdCell, fFDCosSpeed, fFDCosSpeedVsScore, fFDDeltaT, fFDDeltaX, fFDDeltaY, fFDDeltaZ, fFDDiffLastFirstLivePlanes, fFDDTvsDY, fFDKalmanAngle, fFDKalmanBakCell, fFDKalmanDirX, fFDKalmanDirY, fFDKalmanDirZ, fFDKalmanFwdCell, fFDKalSpeed, fFDPlanesToBack, fFDPlanesToFront, fFDPt, fFDSliceCosAngleShw, fFDSliceDistShwVtx, fFDSliceFracShwHits, fFDSliceHits, fFDSliceHitsPerPlane, fFDSliceHitsViewAsym, fFDSliceHitsXView, fFDSliceHitsYView, fFDSpeed, fFDTrackLength, fFDVtxDistToBack, fFDVtxDistToBottom, fFDVtxDistToEastWall, fFDVtxDistToFront, fFDVtxDistToTop, fFDVtxDistToWestWall, fGoodPOTByRun, fHitTree, fLEM, fLID, fNDCalEperHit, fNDCalorimetricE, fNDHadCalEnergySum, fNDKalmanBakCell, fNDKalmanFwdCell, fNDMaxPlaneHit, fNDMinPlaneHit, fNDNShowers, fNDPlanesToFront, fNDShowerCalE, fNDShowerCalEfraction, fNDShowerDistToVtx, fNDShowerLength, fNDShwEndX, fNDShwEndY, fNDShwEndZ, fNDShwVtxX, fNDShwVtxY, fNDShwVtxZ, fNDSliceHits, fNDSliceNHitsPerPlane, fNDSliceNoShowerCalE, fNDTrackStartZ, fNDTrackStopZ, fNDVtxX, fNDVtxY, fNDVtxZ, fNeutrinoE, fNueCutParams, fNueInstance, fNueSRParams, fNueTree, fNumuCutParams, fNumuInstance, fNumuSRParams, fNumuTree, fPOTByRun, fQePId1Track, fQePId2Track, fReMId, fSliceContigPlanes, fSliceHits, htCalibPE, htCalibTNS, htCalibXYZ, htCell, htDir, htDistFromStart, htDistFromStop, htIndex, htPlane, htTrackIndex, htTrackXYZ, htUncalibTNS, htUncalibW, htUncalibXYZ, and art::TFileDirectory::make().

390 {
391  // make histograms for each selection criteria
393 
394  std::string infost("run/I:subrun/I:event/I");
395 
396  fPOTByRun = tfs->make<TH1F>("potByRun", ";Run Number;POT (#times 10^{15})", 10000, 10000., 20000.);
397  fGoodPOTByRun = tfs->make<TH1F>("goodPOTByRun", ";Run Number;POT (#times 10^{15})", 10000, 10000., 20000.);
398  fCosmicExposureByRun = tfs->make<TH1F>("cosmicExposureByRun", ";Run Number;Exposure (s)", 10000, 10000., 20000.);
399 
400  if( !fNumuInstance.empty() ){
401  fNeutrinoE = tfs->make<TH1F>("neutrinoE", ";E_{#nu} (GeV);Events", 100, 0., 10. );
402  fReMId = tfs->make<TH1F>("remid", ";ReMId;Events", 200, -1., 1. );
403  fSliceHits = tfs->make<TH1F>("sliceHits", ";Hits In Slice;Events", 100, 0., 400. );
404  fSliceContigPlanes = tfs->make<TH1F>("sliceContigPlanes", ";Contiguous Planes in Slice;Events", 50, 0., 50. );
405  fCellsFromEdge = tfs->make<TH1F>("cellsFromEdge", ";Cells From Edge;Events", 50, 0., 50. );
406  fQePId1Track = tfs->make<TH1F>("qepid1Track", ";QePId_{1 track};Events", 200, -1., 1. );
407  fQePId2Track = tfs->make<TH1F>("qepid2Track", ";QePId_{2 tracks};Events", 200, -1., 1. );
408  fFDKalmanFwdCell = tfs->make<TH1F>("fdKalmanFwdCell", ";Kalman Foward Cell;Events", 105, -5., 100. );
409  fFDKalmanBakCell = tfs->make<TH1F>("fdKalmanBakCell", ";Kalman Backward Cell;Events", 105, -5., 100. );
410  fFDCosFwdCell = tfs->make<TH1F>("fdCosFwdCell", ";cos Foward Cell;Events", 105, -5., 100. );
411  fFDCosBakCell = tfs->make<TH1F>("fdCosBakCell", ";cos Backward Cell;Events", 105, -5., 100. );
412  fFDPlanesToFront = tfs->make<TH1F>("fdPlanesToFront", ";Planes to Front;Events", 100, 0., 100. );
413  fFDPlanesToBack = tfs->make<TH1F>("fdPlanesToBack", ";Planes to Back;Events", 100, 0., 100. );
414  fFDKalmanAngle = tfs->make<TH1F>("fdKalmanAngle", ";Kalman Angle;Events", 200, -4., 4. );
415  fFDContPID = tfs->make<TH1F>("fdContPID", ";ContPID;Events", 200, -1., 1. );
416  fFDSliceHits = tfs->make<TH1F>("fdSliceHits", ";Slice Hits;Events", 500, 0., 500. );
417  fFDKalmanDirX = tfs->make<TH1F>("fdKalmanDirX", ";X component of Kalman Direction;Events", 200, -1.1,1.1 );
418  fFDKalmanDirY = tfs->make<TH1F>("fdKalmanDirY", ";Y component of Kalman Direction;Events", 200, -1.1,1.1 );
419  fFDKalmanDirZ = tfs->make<TH1F>("fdKalmanDirZ", ";Z component of Kalman Direction;Events", 200, -1.1,1.1 );
420  fFDKalSpeed = tfs->make<TH1F>("fdKalSpeed", ";Inverse speed from timing fit of Kalman track (ns/cm);Events", 200, -1.1, 1.1);
421  fFDCosSpeed = tfs->make<TH1F>("fdCosSpeed", ";Inverse speed from timing fit of Cosmic track (ns/cm);Events", 200, -1.1, 1.1);
422  fFDDeltaT = tfs->make<TH1F>("fdDeltaT", ";Diff between max and min hit time (ns);Events", 200, -1000,1000 );
423  fFDDeltaX = tfs->make<TH1F>("fdDeltaX", ";Diff between start and end X pos (cm);Events", 400, -2000,2000 );
424  fFDDeltaY = tfs->make<TH1F>("fdDeltaY", ";Diff between start and end Y pos (cm);Events", 400, -2000,2000 );
425  fFDDeltaZ = tfs->make<TH1F>("fdDeltaZ", ";Diff between start and end Z pos (cm);Events", 700, -1000,6000 );
426  fFDSpeed = tfs->make<TH1F>("fdSpeed", ";Speed (m/s);Events", 1000, 0, 10000000 );
427  fFDTrackLength = tfs->make<TH1F>("fdTracklength", ";Tracklength (cm);Events", 500, 0, 5000 );
428  fFDDTvsDY = tfs->make<TH2F>("fdDtVsDY", ";Diff between max and min hit time (ns);Diff between start and end Y pos (cm)", 200, 0, 1000, 400, -2000, 2000);
429  fFDCosSpeedVsScore = tfs->make<TH2F>("fdCosSpeedVsScore", ";Inverse speed of Cosmic track (ns/cm);Hough fit score in forward dir - reverse dir", 200, -0.1, 0.1, 200, -1, 1);
430  fFDDiffLastFirstLivePlanes = tfs->make<TH1F>("fddifflastfirstliveplane", ";Last good diblock in track's section of mask - first", 16, -1, 15);
431 
432  fNDMaxPlaneHit = tfs->make<TH1F>("ndMaxPlaneHit", ";Max Hit Plane;Events", 220, 0., 220. );
433  fNDMinPlaneHit = tfs->make<TH1F>("ndMinPlaneHit", ";Min Hit Plane;Events", 220, 0., 220. );
434  fNDTrackStartZ = tfs->make<TH1F>("ndTrackStartZ", ";Track Starting Z Pos (cm);Events", 200, 0., 2000.);
435  fNDTrackStopZ = tfs->make<TH1F>("ndTrackStopZ", ";Track Ending Z Pos (cm);Events", 200, 0., 2000.);
436  fNDKalmanFwdCell = tfs->make<TH1F>("ndKalmanFwdCell", ";Kalman Foward Cell;Events", 100, 0., 100. );
437  fNDKalmanBakCell = tfs->make<TH1F>("ndKalmanBakCell", ";Kalman Backward Cell;Events", 100, 0., 100. );
438  fNDHadCalEnergySum = tfs->make<TH1F>("ndHadCalEnergySum", ";Hadronic Energy (GeV);Events", 100, 0., 10. );
439 
440  std::string paramst("fBestPIDTrack/F:fNeutrinoE/F:fNDHadronicCal/F:");
441  paramst.append("fReMIdVal/D:fQePIDVal/F:fQePIDNTracks/F:");
442  paramst.append("fSliceHits/F:fSliceContigPlanes/F:fSliceMaxPlane/F:fSliceMinPlane/F:fSliceCellsFromEdge/F:");
443  paramst.append("fPlanesToFront/F:fPlanesToBack/F:fCosRejKalFwdCell/F:fCosRejKalBakCell/F:fCosRejAngleKal/F:");
444  paramst.append("fCosRejCosBakCell/F:fCosRejCosFwdCell/F:fCosRejConCosPID/F:fCosRejKalBakCellND/F:fCosRejKalFwdCellND/F:");
445 
446  std::string srst = paramst;
447  srst.append("fCosRejKalYPosAtTrans/F:fTrackStartZ/F:fTrackStopZ/F:fTrackCCE/F");
448 
449  paramst.append("fCosRejKalYPosAtTrans/F:fTrackStart:fTrackStop:fTrackCCE/F:");
450  paramst.append("fTrackStartTime/F:fTrackStopTime/F:");
451  paramst.append("fCosRejKalSpeed/F:fCosRejCosSpeed/F:fCosRejCosDirScore/F:");
452  paramst.append("fDiffLastFirstLiveDiblock/F");
453  fNumuTree = tfs->make<TTree>("Numu", "Skimmed Numu Interaction Cut Parameters");
454  fNumuTree->Branch("info", &fEventInfo, infost .c_str());
455  fNumuTree->Branch("params", &fNumuCutParams, paramst.c_str());
456  fNumuTree->Branch("sr", &fNumuSRParams, srst.c_str());//paramst.c_str());
457  }
458 
459  fCosmicTree = tfs->make<TTree>("FDCosmics", "Information for cosmic rej");
460  fCosmicTree->Branch("NeutrinoE", &crNeutrinoE);
461  fCosmicTree->Branch("ReMId", &crReMIdVal);
462  fCosmicTree->Branch("SliceHits", &crSliceHits);
463  fCosmicTree->Branch("SliceXViewHits", &crSliceXViewHits);
464  fCosmicTree->Branch("SliceYViewHits", &crSliceYViewHits);
465  fCosmicTree->Branch("SliceContigPlanes", &crSliceContigPlanes);
466  fCosmicTree->Branch("SliceMaxPlane", &crSliceMaxPlane);
467  fCosmicTree->Branch("SliceMinPlane", &crSliceMinPlane);
468  fCosmicTree->Branch("SliceTotalADC", &crSliceTotalADC);
469  fCosmicTree->Branch("SliceTotalGeV", &crSliceTotalGeV);
470  fCosmicTree->Branch("SliceCalorimetricEnergy", &crSliceCalorimetricEnergy);
471  fCosmicTree->Branch("SliceMinXYZ", &crSliceMinXYZ);
472  fCosmicTree->Branch("SliceMaxXYZ", &crSliceMaxXYZ);
473  fCosmicTree->Branch("SliceMinTNS", &crSliceMinTNS);
474  fCosmicTree->Branch("SliceMaxTNS", &crSliceMaxTNS);
475  fCosmicTree->Branch("CellsFromEdge", &crCellsFromEdge);
476  fCosmicTree->Branch("SliceClosestEdge", &crSliceClosestEdge);
477  fCosmicTree->Branch("PlanesToFront", &crPlanesToFront);
478  fCosmicTree->Branch("PlanesToBack", &crPlanesToBack);
479  fCosmicTree->Branch("SliceNHitsNearEdge0", &crSliceNHitsNearEdge0);
480  fCosmicTree->Branch("SliceNHitsNearEdge1", &crSliceNHitsNearEdge1);
481  fCosmicTree->Branch("SliceNHitsNearEdge2", &crSliceNHitsNearEdge2);
482  fCosmicTree->Branch("SliceNHitsNearEdge3", &crSliceNHitsNearEdge3);
483  fCosmicTree->Branch("SliceNHitsNearEdge4", &crSliceNHitsNearEdge4);
484  fCosmicTree->Branch("SliceNHitsNearEdge5", &crSliceNHitsNearEdge5);
485  fCosmicTree->Branch("SliceNHitsNearEdge6", &crSliceNHitsNearEdge6);
486  fCosmicTree->Branch("SliceNHitsNearEdge7", &crSliceNHitsNearEdge7);
487  fCosmicTree->Branch("SliceNHitsNearEdge8", &crSliceNHitsNearEdge8);
488  fCosmicTree->Branch("SliceNHitsNearEdge9", &crSliceNHitsNearEdge9);
489  fCosmicTree->Branch("KalAngle", &crKalAngle);
490  fCosmicTree->Branch("CosAngle", &crCosAngle);
491  fCosmicTree->Branch("KalFwdCell", &crKalFwdCell);
492  fCosmicTree->Branch("KalBakCell", &crKalBakCell);
493  fCosmicTree->Branch("CosFwdCell", &crCosFwdCell);
494  fCosmicTree->Branch("CosBakCell", &crCosBakCell);
495  fCosmicTree->Branch("KalFwdDist", &crKalFwdDist);
496  fCosmicTree->Branch("KalBakDist", &crKalBakDist);
497  fCosmicTree->Branch("CosFwdDist", &crCosFwdDist);
498  fCosmicTree->Branch("CosBakDist", &crCosBakDist);
499  fCosmicTree->Branch("ConCosPID", &crConCosPID);
500  fCosmicTree->Branch("KalDir", &crKalDir);
501  fCosmicTree->Branch("NTracks3D", &crNTracks3D);
502  fCosmicTree->Branch("TrackStart", &crTrackStart);
503  fCosmicTree->Branch("TrackStop", &crTrackStop);
504  fCosmicTree->Branch("TrackStartTime", &crTrackStartTime);
505  fCosmicTree->Branch("TrackStopTime", &crTrackStopTime);
506  fCosmicTree->Branch("TrackHits", &crTrackHits);
507  fCosmicTree->Branch("TrackXViewHits", &crTrackXViewHits);
508  fCosmicTree->Branch("TrackYViewHits", &crTrackYViewHits);
509  fCosmicTree->Branch("TrackCellsFromEdge", &crTrackCellsFromEdge);
510  fCosmicTree->Branch("TrackClosestEdge", &crTrackClosestEdge);
511  fCosmicTree->Branch("TrackNHitsNearEdge0", &crTrackNHitsNearEdge0);
512  fCosmicTree->Branch("TrackNHitsNearEdge1", &crTrackNHitsNearEdge1);
513  fCosmicTree->Branch("TrackNHitsNearEdge2", &crTrackNHitsNearEdge2);
514  fCosmicTree->Branch("TrackNHitsNearEdge3", &crTrackNHitsNearEdge3);
515  fCosmicTree->Branch("TrackNHitsNearEdge4", &crTrackNHitsNearEdge4);
516  fCosmicTree->Branch("TrackNHitsNearEdge5", &crTrackNHitsNearEdge5);
517  fCosmicTree->Branch("TrackNHitsNearEdge6", &crTrackNHitsNearEdge6);
518  fCosmicTree->Branch("TrackNHitsNearEdge7", &crTrackNHitsNearEdge7);
519  fCosmicTree->Branch("TrackNHitsNearEdge8", &crTrackNHitsNearEdge8);
520  fCosmicTree->Branch("TrackNHitsNearEdge9", &crTrackNHitsNearEdge9);
521  fCosmicTree->Branch("TrackTotalADC", &crTrackTotalADC);
522  fCosmicTree->Branch("TrackTotalGeV", &crTrackTotalGeV);
523  fCosmicTree->Branch("TrackCalorimetricEnergy", &crTrackCalorimetricEnergy);
524  fCosmicTree->Branch("TrackMinXYZ", &crTrackMinXYZ);
525  fCosmicTree->Branch("TrackMaxXYZ", &crTrackMaxXYZ);
526  fCosmicTree->Branch("TrackEarliestHitPos", &crTrackEarliestHitPos);
527  fCosmicTree->Branch("TrackLatestHitPos", &crTrackLatestHitPos);
528  fCosmicTree->Branch("TrackMinTNS", &crTrackMinTNS);
529  fCosmicTree->Branch("TrackMaxTNS", &crTrackMaxTNS);
530  fCosmicTree->Branch("TrackMaxPlane", &crTrackMaxPlane);
531  fCosmicTree->Branch("TrackMinPlane", &crTrackMinPlane);
532  fCosmicTree->Branch("TrackTotalLength", &crTrackTotalLength);
533  fCosmicTree->Branch("CosTrackStart", &crCosTrackStart);
534  fCosmicTree->Branch("CosTrackStop", &crCosTrackStop);
535  fCosmicTree->Branch("CosTrackStartTime", &crCosTrackStartTime);
536  fCosmicTree->Branch("CosTrackStopTime", &crCosTrackStopTime);
537  fCosmicTree->Branch("CosTrackHits", &crCosTrackHits);
538  fCosmicTree->Branch("CosTrackXViewHits", &crCosTrackXViewHits);
539  fCosmicTree->Branch("CosTrackYViewHits", &crCosTrackYViewHits);
540  fCosmicTree->Branch("CosTrackCellsFromEdge", &crCosTrackCellsFromEdge);
541  fCosmicTree->Branch("CosTrackClosestEdge", &crCosTrackClosestEdge);
542  fCosmicTree->Branch("CosTrackNHitsNearEdge0", &crCosTrackNHitsNearEdge0);
543  fCosmicTree->Branch("CosTrackNHitsNearEdge1", &crCosTrackNHitsNearEdge1);
544  fCosmicTree->Branch("CosTrackNHitsNearEdge2", &crCosTrackNHitsNearEdge2);
545  fCosmicTree->Branch("CosTrackNHitsNearEdge3", &crCosTrackNHitsNearEdge3);
546  fCosmicTree->Branch("CosTrackNHitsNearEdge4", &crCosTrackNHitsNearEdge4);
547  fCosmicTree->Branch("CosTrackNHitsNearEdge5", &crCosTrackNHitsNearEdge5);
548  fCosmicTree->Branch("CosTrackNHitsNearEdge6", &crCosTrackNHitsNearEdge6);
549  fCosmicTree->Branch("CosTrackNHitsNearEdge7", &crCosTrackNHitsNearEdge7);
550  fCosmicTree->Branch("CosTrackNHitsNearEdge8", &crCosTrackNHitsNearEdge8);
551  fCosmicTree->Branch("CosTrackNHitsNearEdge9", &crCosTrackNHitsNearEdge9);
552  fCosmicTree->Branch("CosTrackTotalADC", &crCosTrackTotalADC);
553  fCosmicTree->Branch("CosTrackTotalGeV", &crCosTrackTotalGeV);
554  fCosmicTree->Branch("CosTrackCalorimetricEnergy", &crCosTrackCalorimetricEnergy);
555  fCosmicTree->Branch("CosTrackMinXYZ", &crCosTrackMinXYZ);
556  fCosmicTree->Branch("CosTrackMaxXYZ", &crCosTrackMaxXYZ);
557  fCosmicTree->Branch("CosTrackEarliestHitPos", &crCosTrackEarliestHitPos);
558  fCosmicTree->Branch("CosTrackLatestHitPos", &crCosTrackLatestHitPos);
559  fCosmicTree->Branch("CosTrackMinTNS", &crCosTrackMinTNS);
560  fCosmicTree->Branch("CosTrackMaxTNS", &crCosTrackMaxTNS);
561  fCosmicTree->Branch("CosTrackMaxPlane", &crCosTrackMaxPlane);
562  fCosmicTree->Branch("CosTrackMinPlane", &crCosTrackMinPlane);
563  fCosmicTree->Branch("CosTrackTotalLength", &crCosTrackTotalLength);
564  fCosmicTree->Branch("Scatt", &crScatt);
565  fCosmicTree->Branch("QePIDVal", &crQePIDVal);
566  fCosmicTree->Branch("QePIDNTracks", &crQePIDNTracks);
567  fCosmicTree->Branch("KalSpeed", &crKalSpeed);
568  fCosmicTree->Branch("CosSpeed", &crCosSpeed);
569  fCosmicTree->Branch("KalDirScore", &crKalDirScore);
570  fCosmicTree->Branch("CosDirScore", &crCosDirScore);
571  fCosmicTree->Branch("KalScoreDiff", &crKalScoreDiff);
572  fCosmicTree->Branch("CosScoreDiff", &crCosScoreDiff);
573  fCosmicTree->Branch("KalSlope", &crKalSlope);
574  fCosmicTree->Branch("CosSlope", &crCosSlope);
575  fCosmicTree->Branch("KalChiSq", &crKalChiSq);
576  fCosmicTree->Branch("CosChiSq", &crCosChiSq);
577  fCosmicTree->Branch("KalChiDiff", &crKalChiDiff);
578  fCosmicTree->Branch("CosChiDiff", &crCosChiDiff);
579  fCosmicTree->Branch("NCloseLiveDiblocks", &crNCloseLiveDiblocks);
580  fCosmicTree->Branch("Xmin2", &crXmin2);
581  fCosmicTree->Branch("Ymin2", &crYmin2);
582  fCosmicTree->Branch("Zmin2", &crZmin2);
583  fCosmicTree->Branch("Xmax2", &crXmax2);
584  fCosmicTree->Branch("Ymax2", &crYmax2);
585  fCosmicTree->Branch("Zmax2", &crZmax2);
586  fCosmicTree->Branch("scatllh", &crscatllh);
587  fCosmicTree->Branch("dedxllh", &crdedxllh);
588  fCosmicTree->Branch("measfrac", &crmeasfrac);
589  fCosmicTree->Branch("TruthNNu", &crTruthNNu);
590  fCosmicTree->Branch("TruthNuNHitsSlc", &crTruthNuNHitsSlc);
591  fCosmicTree->Branch("TruthNuNHitsTot", &crTruthNuNHitsTot);
592  fCosmicTree->Branch("TruthNuIsCC", &crTruthNuIsCC);
593  fCosmicTree->Branch("TruthNuVertex", &crTruthNuVertex);
594  fCosmicTree->Branch("TruthNuPDG", &crTruthNuPDG);
595  fCosmicTree->Branch("TruthNuPDGOrig", &crTruthNuPDGOrig);
596  fCosmicTree->Branch("TruthNuVertex", &crTruthNuVertex);
597  fCosmicTree->Branch("TruthNuTime", &crTruthNuTime);
598  fCosmicTree->Branch("Run", &crRun);
599  fCosmicTree->Branch("Subrun", &crSubrun);
600  fCosmicTree->Branch("Event", &crEvent);
601 
602  fHitTree = tfs->make<TTree>("HitTree", "Entry per hit in Kalman track");
603  fHitTree->Branch("Index", &htIndex);
604  fHitTree->Branch("TrackIndex", &htTrackIndex);
605  fHitTree->Branch("UncalibTNS", &htUncalibTNS);
606  fHitTree->Branch("UncalibXYZ", &htUncalibXYZ);
607  fHitTree->Branch("CalibTNS", &htCalibTNS);
608  fHitTree->Branch("CalibXYZ", &htCalibXYZ);
609  fHitTree->Branch("CalibPE", &htCalibPE);
610  fHitTree->Branch("W", &htUncalibW);
611  fHitTree->Branch("DistFromStart", &htDistFromStart);
612  fHitTree->Branch("DistFromStop", &htDistFromStop);
613  fHitTree->Branch("Dir", &htDir);
614  fHitTree->Branch("TrackXYZ", &htTrackXYZ);
615  fHitTree->Branch("Plane", &htPlane);
616  fHitTree->Branch("Cell", &htCell);
617 
618  fCosmicHitTree = tfs->make<TTree>("CosmicHitTree", "Entry per hit in Cosmic track");
619  fCosmicHitTree->Branch("Index", &chtIndex);
620  fCosmicHitTree->Branch("TNS", &chtUncalibTNS);
621  fCosmicHitTree->Branch("XYZ", &chtUncalibXYZ);
622  fCosmicHitTree->Branch("W", &chtUncalibW);
623  fCosmicHitTree->Branch("DistFromStart", &chtDistFromStart);
624  fCosmicHitTree->Branch("DistFromStop", &chtDistFromStop);
625  fCosmicHitTree->Branch("Dir", &chtDir);
626  fCosmicHitTree->Branch("TrackXYZ", &chtTrackXYZ);
627  fCosmicHitTree->Branch("Plane", &chtPlane);
628  fCosmicHitTree->Branch("Cell", &chtCell);
629 
630  if( !fNueInstance.empty() ){
631  fLEM = tfs->make<TH1F>("lem", ";LEM;Slices", 100, 0., 1. );
632  fCVN = tfs->make<TH1F>("cvn", ";CVN;Slices", 100, 0., 1. );
633  fLID = tfs->make<TH1F>("lid", ";LID;Slices", 240, -0.1, 1.1);
634  fCalorimetricEnergy = tfs->make<TH1F>("calorimetricE", ";Calorimetric E (GeV)", 100, 0., 50.);
635 
636  fFDSliceHitsPerPlane = tfs->make<TH1F>("fdSliceHitsPerPlane", ";Slice Hits/Plane;Slices", 100, 0., 50.);
637  fFDSliceHitsXView = tfs->make<TH1F>("fdSliceHitsXView", ";Slice Hits (X view);Slices", 100, 0., 100.);
638  fFDSliceHitsYView = tfs->make<TH1F>("fdSliceHitsYView", ";Slice Hits (Y view);Slices", 100, 0., 100.);
639  fFDSliceHitsViewAsym = tfs->make<TH1F>("fdSliceHitsViewAsym", ";Slice Hits View Asymmetry;Slices", 200, -10., 10.);
640  fFDSliceCosAngleShw = tfs->make<TH1F>("fdSliceCosAngleShw", ";cos(#theta_{showers});Slices", 100, -1., 1.);
641  fFDSliceDistShwVtx = tfs->make<TH1F>("fdSliceDistShwVtx", ";#DeltaL_{shower,vertex};Slices", 100, 0., 400.);
642  fFDSliceFracShwHits = tfs->make<TH1F>("fdSliceFracShwHits", ";Hits_{shower}/Hits_{slice};Slices", 50, 0., 2.);
643  fFDVtxDistToEastWall = tfs->make<TH1F>("fdVtxDistToEastWall", ";Distance to East Wall;Showers", 100, 0., 100.);
644  fFDVtxDistToWestWall = tfs->make<TH1F>("fdVtxDistToWestWall", ";Distance to West Wall;Showers", 100, 0., 100.);
645  fFDVtxDistToTop = tfs->make<TH1F>("fdVtxDistToTop", ";Distance to Top;Showers", 100, 0., 400.);
646  fFDVtxDistToBottom = tfs->make<TH1F>("fdVtxDistToBottom", ";Distance to Bottom;Showers", 100, 0., 400.);
647  fFDVtxDistToFront = tfs->make<TH1F>("fdVtxDistToFront", ";Distance to Front;Showers", 100, 0., 100.);
648  fFDVtxDistToBack = tfs->make<TH1F>("fdVtxDistToBack", ";Distance to Back;Showers", 100, 0., 400.);
649  fFDPt = tfs->make<TH1F>("fdPt", ";p_{t};Showers", 50, 0., 2.);
650 
651  // now the ND histograms
652  fNDNShowers = tfs->make<TH1F>("ndNshowers", ";Nnumber of Showers;Slices", 10, 0., 10.);
653  fNDShowerCalE = tfs->make<TH1F>("ndshowercale", ";Shower E_{calorimetric};Slices", 120, 0., 6.);
654  fNDShowerCalEfraction = tfs->make<TH1F>("ndshowercalefrac", ";Fraction of Slice E in Shower;Slices", 110, 0., 1.1);
655  fNDSliceNoShowerCalE = tfs->make<TH1F>("ndnoshowercale", ";(Slice -Shower) E_{calorimetric};Slices",120, 0., 6.);
656  fNDShowerDistToVtx = tfs->make<TH1F>("ndNshowergap", ";Shower Distance to Vertex;Slices", 25, 0., 200.);
657  fNDVtxX = tfs->make<TH1F>("ndVtxX", ";x_{vertex} (cm);Showers", 200, -200., 200.);
658  fNDVtxY = tfs->make<TH1F>("ndVtxY", ";y_{vertex} (cm);Showers", 200, -200., 200.);
659  fNDVtxZ = tfs->make<TH1F>("ndVtxZ", ";z_{vertex} (cm);Showers", 700, 0., 1400.);
660  fNDShwVtxX = tfs->make<TH1F>("ndShwVtxX", ";x_{vertex} (cm);Showers", 200, -200., 200.);
661  fNDShwVtxY = tfs->make<TH1F>("ndShwVtxY", ";y_{vertex} (cm);Showers", 200, -200., 200.);
662  fNDShwVtxZ = tfs->make<TH1F>("ndShwVtxZ", ";z_{vertex} (cm);Showers", 700, 0., 1400.);
663  fNDShwEndX = tfs->make<TH1F>("ndShwEndX", ";x_{vertex} (cm);Showers", 200, -200., 200.);
664  fNDShwEndY = tfs->make<TH1F>("ndShwEndY", ";y_{vertex} (cm);Showers", 200, -200., 200.);
665  fNDShwEndZ = tfs->make<TH1F>("ndShwEndZ", ";z_{vertex} (cm);Showers", 700, 0., 1400.);
666  fNDSliceHits = tfs->make<TH1F>("ndSliceHits", ";Hits in Slice;Slices", 500, 0., 500.);
667  fNDCalorimetricE = tfs->make<TH1F>("ndCalorimetricE", ";E_{calorimetric} (GeV);Showers", 120, 0., 6.);
668  fNDCalEperHit = tfs->make<TH1F>("ndCalEperHit", ";E_{calorimetric}/ hit ;Slices", 100, 0., 0.1);
669  fNDShowerLength = tfs->make<TH1F>("ndShowerLength", ";L_{shower} (cm);Showers", 700, 0., 1400.);
670  fNDSliceNHitsPerPlane = tfs->make<TH1F>("ndSliceHitsPerPlane", ";Slice Hits/Plane;Slices", 100, 0., 50.);
671  fNDPlanesToFront = tfs->make<TH1F>("ndPlanesToFront", ";Planes to Front;Events", 600, 0., 600.);
672 
673  std::string paramst("fNumPlane/F:fPlanesToFront/F:fCellsPerPlane/F:fNumXCell/F:fNumYCell/F:");
674  paramst.append("fSliceHits/F:fHitAsymmetry/F:fNumShowers/F:fShowerNumXCell/F:fShowerNumYCell/F:");
675  paramst.append("fCosShowers/F:fShowerDirX/F:fShowerDirY/F:fShowerDirZ/F:");
676  paramst.append("fShowerCalE/F:fShowerLength/F:fFracShowerHits/F:");
677  paramst.append("fShowerVertexX/F:fShowerVertexY/F:fShowerVertexZ/F:");
678  paramst.append("fShowerEndX/F:fShowerEndY/F:fShowerEndZ/F:");
679  paramst.append("fShowerMinX/F:fShowerMinY/F:fShowerMaxX/F:fShowerMaxY/F:fShowerMinVtxZ/F:");
680  paramst.append("fShowerMinEndZ/F:fShowerMaxVtxZ/F:fShowerMaxEndZ/F:fShowerVtxDist/F:fShowerPt/F:");
681  paramst.append("fSliceMinPlane/F:fVertexX/F:fVertexY/F:fVertexZ/F:fVertexMaxW/F:fCalorimetricE/F:fNueEnergy/F:");
682  paramst.append("fMinWestDist/F:fMinEastDist/F:fMinBotDist/F:fMinTopDist/F:fMinFrontDist/F:fMinBackDist:");
683  paramst.append("fLIDVal/F:fLEMVal/F:fCVNVal/F");
684  fNueTree = tfs->make<TTree>("Nue", "Skimmed Nue Interaction Cut Parameters");
685  fNueTree->Branch("info", &fEventInfo, infost .c_str());
686  fNueTree->Branch("params", &fNueCutParams, paramst.c_str());
687  fNueTree->Branch("sr", &fNueSRParams, paramst.c_str());
688 
689  }
690 
691  // scan all the parameters in an interactive ROOT session
692  // Numu->Scan("info.run:info.event:params.fBestPIDTrack:params.fNeutrinoE:params.fNDHadronicCal:params.fReMIdVal:params.fQePIDVal:params.fQePIDNTracks:params.fSliceHits:params.fSliceContigPlanes:params.fSliceMaxPlane:params.fSliceMinPlane:params.fSliceCellsFromEdge:params.fPlanesToFront:params.fPlanesToBack:params.fCosRejKalFwdCell:params.fCosRejKalBakCell:params.fCosRejAngleKal:params.fCosRejCosBakCell:params.fCosRejCosFwdCell:params.fCosRejConCosPID:params.fCosRejKalBakCellND:params.fCosRejKalFwdCellND:params.fCosRejKalYPosAtTrans:params.fTrackStartZ:params.fTrackStopZ:params.fTrackCCE");
693 
694  // compare the values from skimming vs CAF
695  // Numu->Scan("info.run:info.event:params.fBestPIDTrack/sr.fBestPIDTrack:params.fNeutrinoE/sr.fNeutrinoE:params.fNDHadronicCal/sr.fNDHadronicCal:params.fReMIdVal/sr.fReMIdVal:params.fQePIDVal/sr.fQePIDVal:params.fQePIDNTracks/sr.fQePIDNTracks:params.fSliceHits/sr.fSliceHits:params.fSliceContigPlanes/sr.fSliceContigPlanes:params.fSliceMaxPlane/sr.fSliceMaxPlane:params.fSliceMinPlane/sr.fSliceMinPlane:params.fSliceCellsFromEdge/sr.fSliceCellsFromEdge:params.fPlanesToFront/sr.fPlanesToFront:params.fPlanesToBack/sr.fPlanesToBack:params.fCosRejKalFwdCell/sr.fCosRejKalFwdCell:params.fCosRejKalBakCell/sr.fCosRejKalBakCell:params.fCosRejAngleKal/sr.fCosRejAngleKal:params.fCosRejCosBakCell/sr.fCosRejCosBakCell:params.fCosRejCosFwdCell/sr.fCosRejCosFwdCell:params.fCosRejConCosPID/sr.fCosRejConCosPID:params.fCosRejKalBakCellND/sr.fCosRejKalBakCellND:params.fCosRejKalFwdCellND/sr.fCosRejKalFwdCellND:params.fCosRejKalYPosAtTrans/sr.fCosRejKalYPosAtTrans:params.fTrackStartZ/sr.fTrackStartZ:params.fTrackStopZ/sr.fTrackStopZ:params.fTrackCCE/sr.fTrackCCE");
696 
697 // Nue->Scan("info.run:info.event:params.fNumPlane:params.fCellsPerPlane:fNumXCell:params.fNumYCell:params.fSliceHits:params.fHitAsymmetry:params.fCosShowers:params.fShowerDirX:params.fShowerDirY:params.fShowerDirZ:params.fShowerLength:params.fFracShowerHits:params.fShowerVertexX:params.fShowerVertexY:params.fShowerVertexZ:params.fShowerEndX:params.fShowerEndY:params.fShowerEndZ:params.fShowerMaxW:params.fShowerMaxW:params.fShowerMinVtxZ:params.fShowerMinEndZ:params.fShowerMaxVtxZ:params.fShowerMaxEndZ:params.fShowerVtxDist:params.fShowerPt:params.fVertexX:params.fVertexY:params.fVertexZ:params.fVertexMaxW:params.fCalorimetricE:params.fMinWestDist:params.fMinEastDist:params.fMinBotDist:params.fMinTopDist:params.fMinFrontDist:params.fMinBackDist");
698 
699 
700  return;
701 }
std::string fNumuInstance
label for numu products
std::string fNueInstance
label for nue products
NumuCutParameters fNumuSRParams
from the ParametersNumu object
TTree * fCosmicTree
Tree filled directly from cosrej etc. objects (bypasses parameters struct)
TTree * fNumuTree
Tree to hold values of cut parameters for selected events.
NueCutParameters fNueCutParams
from the ParametersNumu object
TTree * fHitTree
Tree to hold information about every hit in a kalman track.
::xsd::cxx::tree::string< char, simple_type > string
Definition: Database.h:154
T * make(ARGS...args) const
NueCutParameters fNueSRParams
from the ParametersNumu object
EventInfo fEventInfo
keep track of run and subrun, and event
NumuCutParameters fNumuCutParams
from the ParametersNumu object
TTree * fCosmicHitTree
Tree to hold information about every hit in a kalman track.
TTree * fNueTree
Tree to hold values of cut parameters for selected events.
void skim::CosmicBeamComparison::beginRun ( art::Run const &  r)
overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 704 of file CosmicBeamComparison_module.cc.

References ds::DetectorService::DetId(), fDetId, and fDetService.

705 {
706  fDetId = fDetService->DetId();
707 
708  return;
709 }
novadaq::cnv::DetId DetId() const
What detector are we in?
art::ServiceHandle< ds::DetectorService > fDetService
detector service
novadaq::cnv::DetId fDetId
the id of the detector we are looking at
void skim::CosmicBeamComparison::beginSubRun ( art::SubRun const &  sr)
overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 712 of file CosmicBeamComparison_module.cc.

References fCosmicExposureByRun, fGoodPOTByRun, fPOTByRun, art::DataViewImpl::getManyByType(), LOG_DEBUG, pot, and art::SubRun::run().

713 {
714  // get the CosmicExposure and POTSum out of the subrun
715  //art::Handle<sumdata::CosmicExposure> ce;
716  //art::Handle<sumdata::POTSum> pot;
717  //sr.getByLabel(fCosmicExposureLabel, ce );
718  //sr.getByLabel(fPOTSumLabel, pot);
719 
720 // if( ce.isValid() ){
721 // LOG_DEBUG("CosmicBeamComparison") << sr.run() << " " << ce->totlivetime;
722 // fCosmicExposureByRun->Fill(sr.run(), ce->totlivetime);
723 // }
724 // else LOG_DEBUG("CosmicBeamComparison") << "No cosmic exposure found";
725 // if( pot.isValid() ){
726 // LOG_DEBUG("CosmicBeamComparison") << sr.run() << " " << pot->totgoodpot;
727 // fPOTByRun->Fill(sr.run(), 1.e-16 * pot->totgoodpot);
728 // }
729 // else LOG_DEBUG("CosmicBeamComparison") << "No pot summary found";
730 
731  std::vector<art::Handle<sumdata::CosmicExposure>> ce;
732  std::vector<art::Handle<sumdata::POTSum>> pot;
733  sr.getManyByType(ce );
734  sr.getManyByType(pot);
735  if( ce.size() > 0 ){
736  LOG_DEBUG("CosmicBeamComparison") << sr.run() << " " << ce[0]->totlivetime;
737  fCosmicExposureByRun->Fill(sr.run(), ce[0]->totlivetime);
738  }
739  else LOG_DEBUG("CosmicBeamComparison") << "No cosmic exposure found";
740  if( pot.size() > 0 ){
741  LOG_DEBUG("CosmicBeamComparison") << sr.run() << " " << pot[0]->totgoodpot;
742  fPOTByRun->Fill(sr.run(), 1.e-3 * pot[0]->totpot);
743  fGoodPOTByRun->Fill(sr.run(), 1.e-3 * pot[0]->totgoodpot);
744  }
745  else LOG_DEBUG("CosmicBeamComparison") << "No pot summary found";
746 
747 
748  return;
749 }
#define LOG_DEBUG(stream)
Definition: Messenger.h:149
#define pot
static constexpr Double_t sr
Definition: Munits.h:164
detail::CachedProducts& art::EventObserverBase::cachedProducts ( )
inlineprotectedinherited

Definition at line 79 of file EventObserverBase.h.

References art::EventObserverBase::selectors_.

80  {
81  return selectors_;
82  }
detail::CachedProducts selectors_
template<typename T , BranchType = InEvent>
ProductToken<T> art::Consumer::consumes ( InputTag const &  )
inherited
template<typename T , art::BranchType BT>
art::ProductToken<T> art::Consumer::consumes ( InputTag const &  it)
inherited

Definition at line 146 of file Consumer.h.

References art::InputTag::instance(), PandAna.reco_validation.prod5_pid_validation::invalid, art::InputTag::label(), art::InputTag::process(), and T.

147 {
148  if (!moduleContext_)
149  return ProductToken<T>::invalid();
150 
151  consumables_[BT].emplace_back(ConsumableType::Product,
152  TypeID{typeid(T)},
153  it.label(),
154  it.instance(),
155  it.process());
156  return ProductToken<T>{it};
157 }
set< int >::iterator it
static ProductToken< T > invalid()
Definition: ProductToken.h:47
ConsumableProducts consumables_
Definition: Consumer.h:137
double T
Definition: Xdiff_gwt.C:5
bool moduleContext_
Definition: Consumer.h:135
template<typename T , art::BranchType BT>
void art::Consumer::consumesMany ( )
inherited

Definition at line 161 of file Consumer.h.

References T.

162 {
163  if (!moduleContext_)
164  return;
165 
166  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
167 }
ConsumableProducts consumables_
Definition: Consumer.h:137
double T
Definition: Xdiff_gwt.C:5
bool moduleContext_
Definition: Consumer.h:135
template<typename Element , BranchType = InEvent>
ViewToken<Element> art::Consumer::consumesView ( InputTag const &  )
inherited
template<typename T , art::BranchType BT>
art::ViewToken<T> art::Consumer::consumesView ( InputTag const &  it)
inherited

Definition at line 171 of file Consumer.h.

References art::InputTag::instance(), PandAna.reco_validation.prod5_pid_validation::invalid, art::InputTag::label(), art::InputTag::process(), and T.

172 {
173  if (!moduleContext_)
174  return ViewToken<T>::invalid();
175 
176  consumables_[BT].emplace_back(ConsumableType::ViewElement,
177  TypeID{typeid(T)},
178  it.label(),
179  it.instance(),
180  it.process());
181  return ViewToken<T>{it};
182 }
set< int >::iterator it
static ViewToken< Element > invalid()
Definition: ProductToken.h:75
ConsumableProducts consumables_
Definition: Consumer.h:137
double T
Definition: Xdiff_gwt.C:5
bool moduleContext_
Definition: Consumer.h:135
base_engine_t& art::EngineCreator::createEngine ( seed_t  seed)
inherited
base_engine_t& art::EngineCreator::createEngine ( seed_t  seed,
std::string const &  kind_of_engine_to_make 
)
inherited
base_engine_t& art::EngineCreator::createEngine ( seed_t  seed,
std::string const &  kind_of_engine_to_make,
label_t const &  engine_label 
)
inherited
CurrentProcessingContext const* art::EDAnalyzer::currentContext ( ) const
protectedinherited
void skim::CosmicBeamComparison::FillNueHistograms ( art::Event const &  e)
private

Definition at line 787 of file CosmicBeamComparison_module.cc.

References skim::ParametersNue::CalorimetricE(), skim::ParametersNue::CellsPerPlane(), skim::ParametersNue::CosShowers(), skim::ParametersNue::CVNVal(), art::Handle< T >::failedToGet(), fCalorimetricEnergy, fCVN, fDetId, fFDPt, fFDSliceCosAngleShw, fFDSliceDistShwVtx, fFDSliceFracShwHits, fFDSliceHitsPerPlane, fFDSliceHitsViewAsym, fFDSliceHitsXView, fFDSliceHitsYView, fFDVtxDistToBack, fFDVtxDistToBottom, fFDVtxDistToEastWall, fFDVtxDistToFront, fFDVtxDistToTop, fFDVtxDistToWestWall, skim::SkimmingUtils::FillSliceVector(), fLEM, fLID, fmc, fNDCalEperHit, fNDCalorimetricE, fNDNShowers, fNDPlanesToFront, fNDShowerCalE, fNDShowerCalEfraction, fNDShowerDistToVtx, fNDShowerLength, fNDShwEndX, fNDShwEndY, fNDShwEndZ, fNDShwVtxX, fNDShwVtxY, fNDShwVtxZ, fNDSliceHits, fNDSliceNHitsPerPlane, fNDSliceNoShowerCalE, fNDVtxX, fNDVtxY, fNDVtxZ, fNueCutParams, fNueInstance, fNueTree, skim::ParametersNue::FracShowerHits(), fSkimmedLabel, art::DataViewImpl::getByLabel(), skim::ParametersNue::HitAsymmetry(), novadaq::cnv::kFARDET, novadaq::cnv::kNEARDET, skim::ParametersNue::LEMVal(), skim::ParametersNue::LIDVal(), LOG_DEBUG, LOG_WARNING, skim::ParametersNue::MinBackDist(), skim::ParametersNue::MinBotDist(), skim::ParametersNue::MinEastDist(), skim::ParametersNue::MinFrontDist(), skim::ParametersNue::MinTopDist(), skim::ParametersNue::MinWestDist(), skim::ParametersNue::NumShowers(), skim::ParametersNue::ParametersStruct(), skim::ParametersNue::PlanesToFront(), skim::ParametersNue::ProngLength(), skim::ParametersNue::ShowerCalE(), skim::ParametersNue::ShowerEnd(), skim::ParametersNue::ShowerNumXCell(), skim::ParametersNue::ShowerNumYCell(), skim::ParametersNue::ShowerPt(), skim::ParametersNue::ShowerVertex(), skim::ParametersNue::ShowerVtxDist(), skim::ParametersNue::SliceHits(), getGoodRuns4SAM::tag, and skim::ParametersNue::Vertex().

Referenced by analyze().

788 {
789  // get the relevant data products out of the event record
790  // need the slice, the vertex for the slice, and the showers in the slice
792 
794  e.getByLabel(tag, sliceHandle);
795  if( sliceHandle.failedToGet() ){
796  LOG_WARNING("CosmicBeamComparison") << "No slices passing nue selection in this event";
797  return;
798  }
799 
800  // Get the standard record to compare the values used for the cuts
801  //art::FindOne<caf::StandardRecord> fosr(sliceHandle, e, tag);
802 
803  art::FindMany<rb::Shower> fms (sliceHandle, e, tag);
804  art::FindMany<rb::Prong> fmp (sliceHandle, e, tag);
805  art::FindMany<rb::Vertex> fmv (sliceHandle, e, tag);
806  art::FindMany<cosrej::NueCosRej> fmc (sliceHandle, e, tag);
807  art::FindMany<slid::EventLID> fmlid (sliceHandle, e, tag);
808  art::FindMany<lem::PIDDetails> fmlem (sliceHandle, e, tag);
809  art::FindMany<cvn::Result> fmcvn (sliceHandle, e, tag);
810  art::FindMany<cvn::Result> fmocvn (sliceHandle, e, tag);
811  art::FindMany<presel::Veto> fmnueveto(sliceHandle, e, tag);
812 
813 
815 
816  for(size_t s = 0; s < sliceHandle->size(); ++s){
817 
818  std::vector<rb::Shower const*> showers;
819  std::vector<rb::Prong const*> prongs;
820  std::vector<rb::Vertex const*> vertices;
821  std::vector<slid::ShowerLID const*> slids;
822  std::vector<slid::EventLID const*> elids;
823  std::vector<lem::PIDDetails const*> lems;
824  std::vector<cvn::Result const*> cvns;
825  std::vector<presel::Veto const*> nuevetos;
826  std::vector<cosrej::NueCosRej const*> cosrejs;
827 
828  su.FillSliceVector(s, fmv, vertices);
829  su.FillSliceVector(s, fmc, cosrejs);
830  su.FillSliceVector(s, fms, showers);
831  su.FillSliceVector(s, fmp, prongs);
832  su.FillSliceVector(s, fmlid, elids);
833  su.FillSliceVector(s, fmlem, lems);
834  su.FillSliceVector(s, fmcvn, cvns);
835  su.FillSliceVector(s, fmnueveto, nuevetos);
836  su.FillSliceVector(e, tag, showers, slids);
837 
838  if(vertices.size() != 1 ||
839  cosrejs .size() != 1 ||
840  lems .size() != 1 ||
841  cvns .size() != 1 ||
842  elids .size() != 1 ||
843  showers .size() < 1 ||
844  prongs .size() < 1 ||
845  slids .size() != showers.size()){
846  LOG_DEBUG("CosmicBeamComparison")
847  << "Wrong number of vertices: " << vertices.size()
848  << " : 1 or cosrejs: " << cosrejs .size()
849  << " : 1 or lids: " << elids .size()
850  << " : 1 or lems: " << lems .size()
851  << " : 1 or cvns: " << cvns .size()
852  << " : 1 or showers: " << showers .size()
853  << " : or prongs: " << prongs .size()
854  << " or slids: " << slids .size()
855  << " for the skimmed slice, that should never happen";
856  continue;
857  }
858 
859  // make the parameters object
860  skim::ParametersNue params((*sliceHandle)[s],
861  *(vertices[0]),
862  *(cosrejs[0]),
863  showers,
864  prongs,
865  slids,
866  *(elids[0]),
867  *(lems[0]),
868  *(cvns[0]),
869  *(nuevetos[0]));
870 
871 
872  //skim::ParametersNue sr(fosr.at(s).ref());
873 
874  fNueCutParams = params.ParametersStruct();
875  //fNueSRParams = sr .ParametersStruct();
876 
877  //if (fCompareNueCuts) params.Compare(fNueSRParams);
878 
879  fLEM ->Fill(params.LEMVal());
880  fCVN ->Fill(params.CVNVal());
881  fLID ->Fill(params.LIDVal());
882  fCalorimetricEnergy->Fill(params.CalorimetricE());
883 
885  fFDSliceHitsPerPlane->Fill(params.CellsPerPlane());
886  fFDSliceHitsXView ->Fill(params.ShowerNumXCell());
887  fFDSliceHitsYView ->Fill(params.ShowerNumYCell());
888  fFDSliceHitsViewAsym->Fill(params.HitAsymmetry());
889  fFDSliceCosAngleShw ->Fill(params.CosShowers());
890  fFDSliceDistShwVtx ->Fill(params.ShowerVtxDist());
891  fFDSliceFracShwHits ->Fill(params.FracShowerHits());
892  fFDVtxDistToEastWall->Fill(params.MinEastDist());
893  fFDVtxDistToWestWall->Fill(params.MinWestDist());
894  fFDVtxDistToTop ->Fill(params.MinTopDist());
895  fFDVtxDistToBottom ->Fill(params.MinBotDist());
896  fFDVtxDistToFront ->Fill(params.MinFrontDist());
897  fFDVtxDistToBack ->Fill(params.MinBackDist());
898  fFDPt ->Fill(params.ShowerPt());
899  }
900 
902  fNDVtxX ->Fill(params.Vertex().X() );
903  fNDVtxY ->Fill(params.Vertex().Y() );
904  fNDVtxZ ->Fill(params.Vertex().Z() );
905  fNDShwVtxX ->Fill(params.ShowerVertex().X());
906  fNDShwVtxY ->Fill(params.ShowerVertex().Y());
907  fNDShwVtxZ ->Fill(params.ShowerVertex().Z());
908  fNDShwEndX ->Fill(params.ShowerEnd().X() );
909  fNDShwEndY ->Fill(params.ShowerEnd().Y() );
910  fNDShwEndZ ->Fill(params.ShowerEnd().Z() );
911  fNDSliceHits ->Fill(params.SliceHits() );
912  fNDCalorimetricE ->Fill(params.CalorimetricE() );
913  fNDShowerLength ->Fill(params.ProngLength() );
914  fNDSliceNHitsPerPlane->Fill(params.CellsPerPlane() );
915  fNDPlanesToFront ->Fill(params.PlanesToFront() );
916  fNDNShowers ->Fill(params.NumShowers() );
917  fNDShowerDistToVtx ->Fill(params.ShowerVtxDist() );
918  fNDShowerCalE ->Fill(params.ShowerCalE() );
919  fNDShowerCalEfraction->Fill(params.ShowerCalE() / params.CalorimetricE());
920  fNDSliceNoShowerCalE ->Fill(params.CalorimetricE() - params.ShowerCalE() );
921  fNDCalEperHit ->Fill(params.CalorimetricE() / params.SliceHits() );
922  }
923 
924  fNueTree->Fill();
925 
926  } // end loop over skimmed slices
927 
928  return;
929 }
std::string fNueInstance
label for nue products
#define LOG_DEBUG(stream)
Definition: Messenger.h:149
TFile * fmc
Definition: bdt_com.C:14
void FillSliceVector(size_t const &s, art::FindMany< T > &fmp, std::vector< T const * > &ptrMap)
std::string fSkimmedLabel
label of module creating skimmed file
const XML_Char * s
Definition: expat.h:262
Far Detector at Ash River, MN.
novadaq::cnv::DetId fDetId
the id of the detector we are looking at
NueCutParameters fNueCutParams
from the ParametersNumu object
Near Detector in the NuMI cavern.
#define LOG_WARNING(category)
Float_t e
Definition: plot.C:35
TTree * fNueTree
Tree to hold values of cut parameters for selected events.
bool failedToGet() const
Definition: Handle.h:196
void skim::CosmicBeamComparison::FillNumuHistograms ( art::Event const &  e)
private

2017

2018

Definition at line 953 of file CosmicBeamComparison_module.cc.

References calcFirstLastLivePlane(), geo::GeometryBase::CellInfo(), chtCell, chtDir, chtDistFromStart, chtDistFromStop, chtIndex, chtPlane, chtTrackXYZ, chtUncalibTNS, chtUncalibW, chtUncalibXYZ, crCellsFromEdge, crConCosPID, crCosAngle, crCosBakCell, crCosBakDist, crCosChiDiff, crCosChiSq, crCosDirScore, crCosFwdCell, crCosFwdDist, crCosScoreDiff, crCosSlope, crCosSpeed, crCosTrackCalorimetricEnergy, crCosTrackCellsFromEdge, crCosTrackClosestEdge, crCosTrackDir, crCosTrackEarliestHitPos, crCosTrackHits, crCosTrackLatestHitPos, crCosTrackMaxPlane, crCosTrackMaxTNS, crCosTrackMaxXYZ, crCosTrackMinPlane, crCosTrackMinTNS, crCosTrackMinXYZ, crCosTrackNHitsNearEdge0, crCosTrackNHitsNearEdge1, crCosTrackNHitsNearEdge2, crCosTrackNHitsNearEdge3, crCosTrackNHitsNearEdge4, crCosTrackNHitsNearEdge5, crCosTrackNHitsNearEdge6, crCosTrackNHitsNearEdge7, crCosTrackNHitsNearEdge8, crCosTrackNHitsNearEdge9, crCosTrackStart, crCosTrackStartTime, crCosTrackStop, crCosTrackStopTime, crCosTrackTotalADC, crCosTrackTotalGeV, crCosTrackTotalLength, crCosTrackXViewHits, crCosTrackYViewHits, crEvent, crKalAngle, crKalBakCell, crKalBakDist, crKalChiDiff, crKalChiSq, crKalDir, crKalDirScore, crKalFwdCell, crKalFwdDist, crKalScoreDiff, crKalSlope, crKalSpeed, crmeasfrac, crNCloseLiveDiblocks, crNeutrinoE, crNTracks3D, crReMIdVal, crRun, crScatt, crSliceCalorimetricEnergy, crSliceClosestEdge, crSliceContigPlanes, crSliceHits, crSliceMaxPlane, crSliceMaxTNS, crSliceMaxXYZ, crSliceMinPlane, crSliceMinTNS, crSliceMinXYZ, crSliceNHitsNearEdge0, crSliceNHitsNearEdge1, crSliceNHitsNearEdge2, crSliceNHitsNearEdge3, crSliceNHitsNearEdge4, crSliceNHitsNearEdge5, crSliceNHitsNearEdge6, crSliceNHitsNearEdge7, crSliceNHitsNearEdge8, crSliceNHitsNearEdge9, crSliceTotalADC, crSliceTotalGeV, crSliceXViewHits, crSliceYViewHits, crSubrun, crTrackCalorimetricEnergy, crTrackCellsFromEdge, crTrackClosestEdge, crTrackEarliestHitPos, crTrackHits, crTrackLatestHitPos, crTrackMaxPlane, crTrackMaxTNS, crTrackMaxXYZ, crTrackMinPlane, crTrackMinTNS, crTrackMinXYZ, crTrackNHitsNearEdge0, crTrackNHitsNearEdge1, crTrackNHitsNearEdge2, crTrackNHitsNearEdge3, crTrackNHitsNearEdge4, crTrackNHitsNearEdge5, crTrackNHitsNearEdge6, crTrackNHitsNearEdge7, crTrackNHitsNearEdge8, crTrackNHitsNearEdge9, crTrackStart, crTrackStartTime, crTrackStop, crTrackStopTime, crTrackTotalADC, crTrackTotalGeV, crTrackTotalLength, crTrackXViewHits, crTrackYViewHits, DEFINE_ART_MODULE(), art::EventID::event(), fCellsFromEdge, fCosmicHitTree, fCosmicTree, fDetId, fFDContPID, fFDCosBakCell, fFDCosFwdCell, fFDCosSpeed, fFDCosSpeedVsScore, fFDDeltaT, fFDDeltaX, fFDDeltaY, fFDDeltaZ, fFDDTvsDY, fFDKalmanAngle, fFDKalmanBakCell, fFDKalmanDirX, fFDKalmanDirY, fFDKalmanDirZ, fFDKalmanFwdCell, fFDKalSpeed, fFDPlanesToBack, fFDPlanesToFront, fFDSliceHits, fFDSpeed, fFDTrackLength, fGeom, fHitTree, skim::SkimmingUtils::FillSliceVector(), PandAna.Demos.pi0_spectra::fmt, fNDHadCalEnergySum, fNDKalmanBakCell, fNDKalmanFwdCell, fNDMaxPlaneHit, fNDMinPlaneHit, fNDTrackStartZ, fNDTrackStopZ, fNeutrinoE, fNumuCutParams, fNumuInstance, fNumuTree, fQePId1Track, fQePId2Track, fReMId, fSkimmedLabel, fSliceContigPlanes, fSliceHits, art::FindMany< ProdB, Data >::get(), nova::dbi::RunHistory::GetConfiguration(), art::DataViewImpl::getValidHandle(), nova::dbi::RunHistory::GoodDiBlockMask(), remid::HighestPIDTrack(), htCalibPE, htCalibTNS, htCalibXYZ, htCell, htDir, htDistFromStart, htDistFromStop, htIndex, htPlane, htTrackIndex, htTrackXYZ, htUncalibTNS, htUncalibW, htUncalibXYZ, art::Event::id(), art::Event::isRealData(), novadaq::cnv::kFARDET, novadaq::cnv::kNEARDET, geo::kXorY, LOG_DEBUG, LOG_VERBATIM, LOG_WARNING, Mag(), cet::sqlite::max(), rb::RecoHit::PECorr(), rh, art::Event::run(), art::Event::subRun(), rb::RecoHit::T(), POTSpillRate::view, W, rb::RecoHit::X(), rb::RecoHit::Y(), and rb::RecoHit::Z().

Referenced by analyze().

954 {
956  art::InputTag cosinputTag(fSkimmedLabel, fNumuInstance+"Cosmics");
957 
958  // get the relevant data products out of the event record
959  // need the slice, rb::Energy, remid::ReMId, cosrej::CosRejObj,
960  // qeef::QePId, numue::NumuE, rb::Track
961  auto sliceHandle = e.getValidHandle< std::vector<rb::Cluster> >(inputTag);
962  if( sliceHandle.failedToGet() ){
963  LOG_WARNING("CosmicBeamComparison") << "No slices passing numu selection in this event";
964  return;
965  }
966 
967  // Get the standard record to compare the values used for the cuts
968  //art::FindOne<caf::StandardRecord> fosr (sliceHandle, e, inputTag);
969 
970  art::FindManyP<rb::Track> fmt (sliceHandle, e, inputTag);
971  art::FindManyP<rb::Track> fmct (sliceHandle, e, cosinputTag);
972  art::FindMany<rb::Energy> fme (sliceHandle, e, inputTag);
973  art::FindMany<cosrej::CosRejObj> fmcro (sliceHandle, e, inputTag);
974  art::FindMany<qeef::QePId> fmqe (sliceHandle, e, inputTag);
975  art::FindMany<numue::NumuE> fmnue (sliceHandle, e, inputTag);
976  art::FindMany<cvn::Result> fmcvn (sliceHandle, e, inputTag);
977  art::FindMany<cvn::Result> fmocvn(sliceHandle, e, inputTag);
978 
979  art::FindMany<cosrej::NueCosRej> fmnuecro (sliceHandle, e, inputTag); // Nue CosRej added for Analysis 2017
980  art::FindMany<rb::Track> fmtks (sliceHandle, e, inputTag);
981  art::FindMany<rb::Shower> fms (sliceHandle, e, inputTag);
982  art::FindMany<rb::Prong> fmp (sliceHandle, e, inputTag);
983  art::FindMany<rb::Vertex> fmv (sliceHandle, e, inputTag);
984  art::FindMany<slid::EventLID> fmlid (sliceHandle, e, inputTag);
985 
986  /*
987  std::vector< art::Ptr<rb::Track> > trackPtrs;
988  std::vector< art::Ptr<rb::Track> > costrackPtrs;
989  std::vector<cosrej::CosRejObj const*> cosrej;
990  std::vector<qeef::QePId const*> qepid;
991  std::vector<numue::NumuE const*> numue;
992  */
993  /*
994  // std::vector<cosrej::NueCosRej const*> nuecosrej; // Nue CosRej added for Analysis 2017
995  std::vector<cvn::Result const*> cvns; /// 2017
996  std::vector<rb::Shower const*> showers;
997  std::vector<rb::Prong const*> prongs;
998  std::vector<rb::Vertex const*> vertices;
999  std::vector<slid::ShowerLID const*> slids;
1000  */
1001  //==> size_t bestPIDTrack = std::numeric_limits<size_t>::max();
1002 
1004 
1005  for(size_t s = 0; s < sliceHandle->size(); ++s){
1006 
1007 
1008  std::vector< art::Ptr<rb::Track> > trackPtrs;
1009  std::vector< art::Ptr<rb::Track> > costrackPtrs;
1010  std::vector<cosrej::CosRejObj const*> cosrej;
1011  std::vector<qeef::QePId const*> qepid;
1012  std::vector<numue::NumuE const*> numue;
1013 
1014  std::vector<cosrej::NueCosRej const*> nuecosrej; // Nue CosRej added for Analysis 2017
1015  std::vector<cvn::Result const*> cvns; /// 2017
1016  std::vector<cvn::Result const*> cvn2017s; /// 2018
1017  std::vector<rb::Track const*> tracks;
1018  std::vector<rb::Shower const*> showers;
1019  std::vector<rb::Prong const*> prongs;
1020  std::vector<rb::Vertex const*> vertices;
1021  std::vector<slid::ShowerLID const*> slids;
1022 
1023  size_t bestPIDTrack = std::numeric_limits<size_t>::max();
1024 
1025 // if(fosr.at(s).ref().hdr.ismc)//only fill truth information if it's there
1026 // {
1027 // crTruthNNu = fosr.at(s).ref().mc.allnus.size();
1028 // if(fosr.at(s).ref().mc.nu.size() > 0)
1029 // {
1030 // crTruthNuNHitsSlc = fosr.at(s).ref().mc.nu[0].nhitslc;
1031 // crTruthNuNHitsTot = fosr.at(s).ref().mc.nu[0].nhittot;
1032 // crTruthNuIsCC = fosr.at(s).ref().mc.nu[0].iscc;
1033 // crTruthNuPDG = fosr.at(s).ref().mc.nu[0].pdg;
1034 // crTruthNuPDGOrig = fosr.at(s).ref().mc.nu[0].pdgorig;
1035 // crTruthNuVertex = fosr.at(s).ref().mc.nu[0].vtx;
1036 // crTruthNuTime = fosr.at(s).ref().mc.nu[0].time;
1037 // }
1038 // }
1039 
1040 
1041 // fmt.get(s, trackPtrs);
1042 // fmct.get(s, costrackPtrs);
1043 
1044  // better have at least one track
1045  if(trackPtrs.size() < 1) continue;
1046 
1047  // get the best pid for the slice - this index will not be the
1048  // same as in the preskimmed slice
1049  bestPIDTrack = remid::HighestPIDTrack(trackPtrs, inputTag, e);
1050 
1051  fmcro.get(s, cosrej);
1052  fmnue.get(s, numue);
1053  fmnuecro.get(s, nuecosrej); // Nue CosRej added for Analysis 2017
1054  fmcvn.get( s, cvns ); // 2017
1055  fmocvn.get(s, cvn2017s); // 2018
1056 
1057  su.FillSliceVector(s, fmtks, tracks);
1058  su.FillSliceVector(s, fmv, vertices);
1059  su.FillSliceVector(s, fms, showers);
1060  su.FillSliceVector(s, fmp, prongs);
1061  su.FillSliceVector(e, inputTag, showers, slids);
1062 
1063 
1064  //fmqe .get(s, qepid);
1065  //qeef::QePId qep;
1066  //if( qepid.size() > 0 ) qep = *(qepid.front());
1067 
1068 
1069  if(vertices.size() != 1 ||
1070  cvns .size() != 1 ||
1071  cvn2017s.size() != 1 ||
1072  showers .size() < 1 ||
1073  prongs .size() < 1 ||
1074  tracks .size() < 1 ||
1075  slids .size() != showers.size()){
1076  LOG_DEBUG("CosmicBeamComparison")
1077  << "Wrong number of vertices: " << vertices.size()
1078  << " : 1 or cvns: " << cvns .size()
1079  << " : 1 or cvn2017s: " << cvn2017s.size()
1080  << " : 1 or showers: " << showers .size()
1081  << " : or prongs: " << prongs .size()
1082  << " : or tracks: " << tracks .size()
1083  << " or slids: " << slids .size()
1084  << " for the skimmed slice, that should never happen";
1085  continue;
1086  }
1087 
1088  // this is not terribly efficient, as we are making the FindOneP
1089  // for every slice in the event, but I don't see a good way around
1090  // that while still getting the bookkeeping right
1091  art::FindOne<rb::Energy> foe (trackPtrs, e, inputTag);
1092  art::FindOne<remid::ReMId> foid(trackPtrs, e, inputTag);
1093 
1094  //LOG_VERBATIM("CosmicBeamComparison") << "Filling params";
1095  skim::ParametersNumu params(bestPIDTrack,
1096  foe.at(bestPIDTrack).ref(),
1097  foid.at(bestPIDTrack).ref(),
1098  *(cosrej.front()),
1099  *(nuecosrej.front()),
1100  //qep,
1101  *(vertices[0]),
1102  tracks,
1103  showers,
1104  prongs,
1105  slids,
1106  *(costrackPtrs[0]), //Not sure what else to do here...
1107  *(numue.front()),
1108  sliceHandle->at(s),
1109  *(trackPtrs[bestPIDTrack]),
1110  *(cvns.front()),
1111  *(cvn2017s.front()),
1112  e.isRealData());
1113 
1114 // skim::ParametersNumu sr(fosr.at(s).ref(),
1115 // e.isRealData());
1116 
1117  fNumuCutParams = params.ParametersStruct();
1118  //LOG_VERBATIM("CosmicBeamComparison") << "fNumuCutParams.fReMIdVal = " << fNumuCutParams.fReMIdVal;
1119 // fNumuSRParams = sr .ParametersStruct();
1120 //
1121 // if (fCompareNumuCuts) params.Compare(fNumuSRParams);
1122 
1123  // Fill sanity check histograms
1124  fNeutrinoE ->Fill(params.NeutrinoE());
1125  fReMId ->Fill(params.ReMIdValue());
1126  fSliceHits ->Fill(params.SliceHits());
1127  fSliceContigPlanes->Fill(params.SliceContigPlanes());
1128  fCellsFromEdge ->Fill(params.SliceCellsFromEdge());
1129  if(params.QePIDNTracks() == 1) fQePId1Track->Fill(params.QePIDValue());
1130  if(params.QePIDNTracks() == 2) fQePId2Track->Fill(params.QePIDValue());
1131 
1133  fFDKalmanFwdCell->Fill(params.CosRejKalFwdCell());
1134  fFDKalmanBakCell->Fill(params.CosRejKalBakCell());
1135  fFDKalmanAngle ->Fill(params.CosRejAngleKal());
1136  fFDCosBakCell ->Fill(params.CosRejCosBakCell());
1137  fFDCosFwdCell ->Fill(params.CosRejCosFwdCell());
1138  fFDContPID ->Fill(params.CosRejConCosPID());
1139  fFDPlanesToBack ->Fill(params.PlanesToBack());
1140  fFDPlanesToFront->Fill(params.PlanesToFront());
1141  fFDSliceHits ->Fill(params.SliceHits());
1142  fFDKalmanDirX ->Fill(params.TrackDir().X());
1143  fFDKalmanDirY ->Fill(params.TrackDir().Y());
1144  fFDKalmanDirZ ->Fill(params.TrackDir().Z());
1145  fFDKalSpeed ->Fill(params.CosRejKalSpeed());
1146  fFDCosSpeed ->Fill(params.CosRejCosSpeed());
1147  fFDDeltaT ->Fill(params.TrackStopTime()-params.TrackStartTime());
1148  fFDDeltaX ->Fill(params.TrackStop().X()-params.TrackStart().X());
1149  fFDDeltaY ->Fill(params.TrackStop().Y()-params.TrackStart().Y());
1150  fFDDeltaZ ->Fill(params.TrackStop().Z()-params.TrackStart().Z());
1151  fFDDTvsDY ->Fill(params.TrackStopTime()-params.TrackStartTime(), params.TrackStop().Y()-params.TrackStart().Y());
1152  fFDCosSpeedVsScore->Fill(params.CosRejCosSpeed(),params.CosRejCosDirScore());
1153  fFDSpeed ->Fill(((params.TrackStop()-params.TrackStart()).Mag() * 1e-2) / ((params.TrackStopTime()-params.TrackStartTime()) * 1e-9));//Convert cm/ns to m/s
1154  fFDTrackLength ->Fill((params.TrackStop()-params.TrackStart()).Mag());
1155 
1156  //Fill CosmicTree - mostly the same as params, but ReMId works
1157  crNeutrinoE = params.NeutrinoE();
1158  crReMIdVal = params.ReMIdValue();
1159  crSliceHits = sliceHandle->at(s).NCell();
1160  crSliceContigPlanes = sliceHandle->at(s).MostContiguousPlanes(geo::kXorY);
1161  crSliceMaxPlane = sliceHandle->at(s).MaxPlane();
1162  crSliceMinPlane = sliceHandle->at(s).MinPlane();
1163  crSliceXViewHits = sliceHandle->at(s).NXCell();
1164  crSliceYViewHits = sliceHandle->at(s).NYCell();
1165  crSliceTotalADC = sliceHandle->at(s).TotalADC();
1166  crSliceTotalGeV = sliceHandle->at(s).TotalGeV();
1167  crSliceCalorimetricEnergy = sliceHandle->at(s).CalorimetricEnergy();
1168  crSliceMinXYZ = sliceHandle->at(s).MinXYZ();
1169  crSliceMaxXYZ = sliceHandle->at(s).MaxXYZ();
1170  crSliceMinTNS = sliceHandle->at(s).MinTNS();
1171  crSliceMaxTNS = sliceHandle->at(s).MaxTNS();
1172  crCellsFromEdge = params.SliceCellsFromEdge();
1173  crSliceClosestEdge = params.SliceClosestEdge();
1174  crSliceNHitsNearEdge0 = params.SliceNHitsNearEdge(0);
1175  crSliceNHitsNearEdge1 = params.SliceNHitsNearEdge(1);
1176  crSliceNHitsNearEdge2 = params.SliceNHitsNearEdge(2);
1177  crSliceNHitsNearEdge3 = params.SliceNHitsNearEdge(3);
1178  crSliceNHitsNearEdge4 = params.SliceNHitsNearEdge(4);
1179  crSliceNHitsNearEdge5 = params.SliceNHitsNearEdge(5);
1180  crSliceNHitsNearEdge6 = params.SliceNHitsNearEdge(6);
1181  crSliceNHitsNearEdge7 = params.SliceNHitsNearEdge(7);
1182  crSliceNHitsNearEdge8 = params.SliceNHitsNearEdge(8);
1183  crSliceNHitsNearEdge9 = params.SliceNHitsNearEdge(9);
1184  //crPlanesToFront = sr.PlanesToFront();//params.PlanesToFront()
1185  //crPlanesToBack = sr.PlanesToBack();//params.PlanesToBack() returns a lot of negative numbers, I'm not entirely convinced by it.
1186  crKalAngle = params.CosRejAngleKal();
1187  crCosAngle = cosrej.front()->AngleCos();
1188  crKalFwdCell = params.CosRejKalFwdCell();
1189  crKalBakCell = params.CosRejKalBakCell();
1190  crCosFwdCell = params.CosRejCosFwdCell();
1191  crCosBakCell = params.CosRejCosBakCell();
1192  crKalFwdDist = cosrej.front()->KalFwdDist();
1193  crKalBakDist = cosrej.front()->KalBakDist();
1194  crCosFwdDist = cosrej.front()->CosFwdDist();
1195  crCosBakDist = cosrej.front()->CosBakDist();
1196  crConCosPID = params.CosRejConCosPID();
1197  crKalDir = params.TrackDir();
1198  crTrackStart = params.TrackStart();
1199  crTrackStop = params.TrackStop();
1200  crTrackStartTime = params.TrackStartTime();
1201  crTrackStopTime = params.TrackStopTime();
1202  crTrackHits = (trackPtrs[bestPIDTrack])->NCell();
1203  crTrackCellsFromEdge = params.TrackNCellsFromEdge();
1204  crTrackClosestEdge = params.TrackClosestEdge();
1205  crTrackNHitsNearEdge0 = params.TrackNHitsNearEdge(0);
1206  crTrackNHitsNearEdge1 = params.TrackNHitsNearEdge(1);
1207  crTrackNHitsNearEdge2 = params.TrackNHitsNearEdge(2);
1208  crTrackNHitsNearEdge3 = params.TrackNHitsNearEdge(3);
1209  crTrackNHitsNearEdge4 = params.TrackNHitsNearEdge(4);
1210  crTrackNHitsNearEdge5 = params.TrackNHitsNearEdge(5);
1211  crTrackNHitsNearEdge6 = params.TrackNHitsNearEdge(6);
1212  crTrackNHitsNearEdge7 = params.TrackNHitsNearEdge(7);
1213  crTrackNHitsNearEdge8 = params.TrackNHitsNearEdge(8);
1214  crTrackNHitsNearEdge9 = params.TrackNHitsNearEdge(9);
1215  crTrackXViewHits = (trackPtrs[bestPIDTrack])->NXCell();
1216  crTrackYViewHits = (trackPtrs[bestPIDTrack])->NYCell();
1217  crTrackEarliestHitPos = params.TrackEarliestHitPos();
1218  crTrackLatestHitPos = params.TrackLatestHitPos();
1219  crTrackTotalADC = (trackPtrs[bestPIDTrack])->TotalADC();//sum of adc of all hits
1220  crTrackTotalGeV = (trackPtrs[bestPIDTrack])->TotalGeV();//simple sum of estimated gev of all hits
1221  crTrackCalorimetricEnergy = (trackPtrs[bestPIDTrack])->CalorimetricEnergy();//simple estimate of neutrino energy
1222  crTrackMinTNS = (trackPtrs[bestPIDTrack])->MinTNS();
1223  crTrackMaxTNS = (trackPtrs[bestPIDTrack])->MaxTNS();
1224  crTrackMinXYZ = (trackPtrs[bestPIDTrack])->MinXYZ();
1225  crTrackMaxXYZ = (trackPtrs[bestPIDTrack])->MaxXYZ();
1226  crTrackMinPlane = (trackPtrs[bestPIDTrack])->MinPlane();
1227  crTrackMaxPlane = (trackPtrs[bestPIDTrack])->MaxPlane();
1228  crTrackTotalLength = (trackPtrs[bestPIDTrack])->TotalLength();
1229  crCosTrackStart = params.CosmicTrackStart();
1230  crCosTrackStop = params.CosmicTrackStop();
1231  crCosTrackStartTime = params.CosTrackStartTime();
1232  crCosTrackStopTime = params.CosTrackStopTime();
1233  crCosTrackDir = params.CosmicTrackDir();
1234  crCosTrackHits = params.CosTrackHits();
1235  crCosTrackCellsFromEdge = params.CosTrackNCellsFromEdge();
1236  crCosTrackClosestEdge = params.CosTrackClosestEdge();
1237  crCosTrackNHitsNearEdge0 = params.CosTrackNHitsNearEdge(0);
1238  crCosTrackNHitsNearEdge1 = params.CosTrackNHitsNearEdge(1);
1239  crCosTrackNHitsNearEdge2 = params.CosTrackNHitsNearEdge(2);
1240  crCosTrackNHitsNearEdge3 = params.CosTrackNHitsNearEdge(3);
1241  crCosTrackNHitsNearEdge4 = params.CosTrackNHitsNearEdge(4);
1242  crCosTrackNHitsNearEdge5 = params.CosTrackNHitsNearEdge(5);
1243  crCosTrackNHitsNearEdge6 = params.CosTrackNHitsNearEdge(6);
1244  crCosTrackNHitsNearEdge7 = params.CosTrackNHitsNearEdge(7);
1245  crCosTrackNHitsNearEdge8 = params.CosTrackNHitsNearEdge(8);
1246  crCosTrackNHitsNearEdge9 = params.CosTrackNHitsNearEdge(9);
1247  crCosTrackXViewHits = params.CosTrackXViewHits();
1248  crCosTrackYViewHits = params.CosTrackXViewHits();
1249  crCosTrackTotalADC = params.CosTrackTotalADC();
1250  crCosTrackTotalGeV = params.CosTrackTotalGeV();
1251  crCosTrackCalorimetricEnergy = params.CosTrackTotalCalorimetricEnergy();
1252  crCosTrackMinTNS = params.CosTrackMinTNS();
1253  crCosTrackMaxTNS = params.CosTrackMaxTNS();
1254  crCosTrackMinXYZ = params.CosTrackMinXYZ();
1255  crCosTrackMaxXYZ = params.CosTrackMaxXYZ();
1256  crCosTrackEarliestHitPos = params.CosmicTrackEarliestHitPos();
1257  crCosTrackLatestHitPos = params.CosmicTrackLatestHitPos();
1258  crCosTrackMinPlane = params.CosTrackMinPlane();
1259  crCosTrackMaxPlane = params.CosTrackMaxPlane();
1260  crCosTrackTotalLength = (costrackPtrs[0])->TotalLength();
1261  crKalSpeed = params.CosRejKalSpeed();
1262  crCosSpeed = params.CosRejCosSpeed();
1263  crCosDirScore = params.CosRejCosDirScore();
1264  crKalDirScore = cosrej.front()->KDirScore();
1265  crCosScoreDiff = cosrej.front()->CScoreDiff();
1266  crKalScoreDiff = cosrej.front()->KScoreDiff();
1267  crKalSlope = cosrej.front()->KalSlope();
1268  crCosSlope = cosrej.front()->CosSlope();
1269  crKalChiSq = cosrej.front()->KalChisq();
1270  crCosChiSq = cosrej.front()->CosChisq();
1271  crKalChiDiff = cosrej.front()->KalChiDiff();
1272  crCosChiDiff = cosrej.front()->CosChiDiff();
1273  crNTracks3D = cosrej.front()->NTracks3D();
1274  crScatt = cosrej.front()->Scatt();
1275  //crQePIDVal = qep.Value();
1276  //crQePIDNTracks = qep.Ntrk();
1277 // crXmin2 = fosr.at(s).ref().sel.contain.xmin2;
1278 // crYmin2 = fosr.at(s).ref().sel.contain.ymin2;
1279 // crZmin2 = fosr.at(s).ref().sel.contain.zmin2;
1280 // crXmax2 = fosr.at(s).ref().sel.contain.xmax2;
1281 // crYmax2 = fosr.at(s).ref().sel.contain.ymax2;
1282 // crZmax2 = fosr.at(s).ref().sel.contain.zmax2;
1283 // crscatllh = fosr.at(s).ref().sel.remid.scatllh;
1284 // crdedxllh = fosr.at(s).ref().sel.remid.dedxllh;
1285  crmeasfrac = foid.at(bestPIDTrack).ref().MeasurementFraction();
1286 
1288  unsigned int dibmask = (rh->GoodDiBlockMask()&rh->GetConfiguration());
1289  std::bitset<14> binary(dibmask);
1290  std::pair<int,int> planesA = calcFirstLastLivePlane(params.SliceMinPlane(), binary, fDetId);
1291  std::pair<int,int> planesB = calcFirstLastLivePlane(params.SliceMaxPlane(), binary, fDetId);
1292  if ((planesA.first != planesB.first) || (planesA.second != planesB.second)) crNCloseLiveDiblocks = -1000;
1293  else crNCloseLiveDiblocks = (planesA.second - planesA.first + 1) * 0.015625;
1294 
1295  crRun = e.run();
1296  crSubrun = e.subRun();
1297  crEvent = e.id().event();
1298 
1299  fCosmicTree->Fill();
1300  }
1301 
1302  LOG_VERBATIM("CosmicBeamComparison") << "About to fill hittree";
1303  for(unsigned int hitnumber = 0; hitnumber < (trackPtrs[bestPIDTrack])->NCell(); hitnumber++){
1304  rb::RecoHit calibhit = (trackPtrs[bestPIDTrack])->RecoHit((trackPtrs[bestPIDTrack])->Cell(hitnumber));//Calibrated hit
1305  htIndex = hitnumber;
1306  if(htIndex==0) htTrackIndex++;
1307  htUncalibTNS = (trackPtrs[bestPIDTrack])->Cell(hitnumber)->TNS();
1308  htPlane = (trackPtrs[bestPIDTrack])->Cell(hitnumber)->Plane();
1309  htCell = (trackPtrs[bestPIDTrack])->Cell(hitnumber)->Cell();
1310  geo::View_t* view = 0;
1311  double hitpos[3];
1312  double dpos[3];//half dimensions of cell in xyz - need it for the cellinfo function
1313  //const geo::OfflineChan offline_channel(htPlane, htCell);
1314  fGeom->CellInfo(htPlane, htCell, view, hitpos, dpos);
1315  //offline_channel.GetCenter
1316  htUncalibXYZ.SetXYZ(hitpos[0], hitpos[1], hitpos[2]);
1317 
1318  htCalibTNS = calibhit.T();
1319  htCalibXYZ.SetX(calibhit.X());
1320  htCalibXYZ.SetY(calibhit.Y());
1321  htCalibXYZ.SetZ(calibhit.Z());
1322  htCalibPE = calibhit.PECorr();
1323 
1324  //*view = (trackPtrs[bestPIDTrack])->Cell(hitnumber)->View();
1325  //if(!view) LOG_VERBATIM("CosmicBeamComparison") << "View not assigned";
1326  //if(htPlane %2 == 0) htXYZ.SetY((trackPtrs[bestPIDTrack])->W((trackPtrs[bestPIDTrack]->Cell(hitnumber)).get()));//was if(*view == geo::kX)
1327  //else htXYZ.SetX((trackPtrs[bestPIDTrack])->W((trackPtrs[bestPIDTrack]->Cell(hitnumber)).get()));
1328  htUncalibW = (trackPtrs[bestPIDTrack])->W((trackPtrs[bestPIDTrack]->Cell(hitnumber)).get());
1329  htDistFromStart = (trackPtrs[bestPIDTrack])->DistanceFromStart(htCalibXYZ.Z());
1330  htDistFromStop = (trackPtrs[bestPIDTrack])->DistanceFromEnd(htCalibXYZ.Z());
1331  htDir = (trackPtrs[bestPIDTrack])->InterpolateDir(htCalibXYZ.Z());
1332  double closestx;
1333  double closesty;
1334  (trackPtrs[bestPIDTrack])->InterpolateXY(htCalibXYZ.Z(), closestx, closesty);
1335  htTrackXYZ.SetXYZ(closestx, closesty, htCalibXYZ.Z());
1336 
1337  fHitTree->Fill();
1338  }
1339 
1340  LOG_VERBATIM("CosmicBeamComparison") << "About to fill cosmic hittree";
1341  for(unsigned int hitnumber = 0; hitnumber < (costrackPtrs[0])->NCell(); hitnumber++){
1342  chtIndex = hitnumber;
1343  chtUncalibTNS = (costrackPtrs[0])->Cell(hitnumber)->TNS();
1344  chtPlane = (costrackPtrs[0])->Cell(hitnumber)->Plane();
1345  chtCell = (costrackPtrs[0])->Cell(hitnumber)->Cell();
1346  geo::View_t* view = 0;
1347  double hitpos[3];
1348  double dpos[3];//half dimensions of cell in xyz - need it for the cellinfo function
1349  //const geo::OfflineChan offline_channel(chtPlane, chtCell);
1350  fGeom->CellInfo(chtPlane, chtCell, view, hitpos, dpos);
1351  //offline_channel.GetCenter
1352  chtUncalibXYZ.SetXYZ(hitpos[0], hitpos[1], hitpos[2]);
1353  //*view = (costrackPtrs[0])->Cell(hitnumber)->View();
1354  //if(!view) LOG_VERBATIM("CosmicBeamComparison") << "View not assigned";
1355  //if(chtPlane %2 == 0) chtUncalibXYZ.SetY((costrackPtrs[0])->W((costrackPtrs[0]->Cell(hitnumber)).get()));//was if(*view == geo::kX)
1356  //else chtUncalibXYZ.SetX((costrackPtrs[0])->W((costrackPtrs[0]->Cell(hitnumber)).get()));
1357  chtUncalibW = (costrackPtrs[0])->W((costrackPtrs[0]->Cell(hitnumber)).get());
1358  chtDistFromStart = (costrackPtrs[0])->DistanceFromStart(chtUncalibXYZ.Z());
1359  chtDistFromStop = (costrackPtrs[0])->DistanceFromEnd(chtUncalibXYZ.Z());
1360  chtDir = (costrackPtrs[0])->InterpolateDir(chtUncalibXYZ.Z());
1361  double closestx;
1362  double closesty;
1363  (costrackPtrs[0])->InterpolateXY(chtUncalibXYZ.Z(), closestx, closesty);
1364  chtTrackXYZ.SetXYZ(closestx, closesty, chtUncalibXYZ.Z());
1365 
1366  fCosmicHitTree->Fill();
1367  }
1368 
1370  fNDMaxPlaneHit ->Fill(params.SliceMaxPlane());
1371  fNDMinPlaneHit ->Fill(params.SliceMinPlane());
1372  fNDTrackStartZ ->Fill(params.TrackStart().Z());
1373  fNDTrackStopZ ->Fill(params.TrackStop().Z());
1374  fNDKalmanBakCell ->Fill(params.CosRejKalBakCellND());
1375  fNDKalmanFwdCell ->Fill(params.CosRejKalFwdCellND());
1376  fNDHadCalEnergySum->Fill(params.NDHadronicCal());
1377  }
1378 
1379  fNumuTree->Fill();
1380 
1381  }// end loop over slices
1382 
1383  return;
1384 }
std::string fNumuInstance
label for numu products
float T() const
Definition: RecoHit.h:39
#define LOG_DEBUG(stream)
Definition: Messenger.h:149
Energy estimators for CC events.
Definition: FillEnergies.h:7
TH2 * rh
Definition: drawXsec.C:5
X or Y views.
Definition: PlaneGeo.h:30
enum geo::_plane_proj View_t
Enumerate the possible plane projections.
float Z() const
Definition: RecoHit.h:38
std::pair< int, int > calcFirstLastLivePlane(int plane, std::bitset< 14 > binary, novadaq::cnv::DetId det)
void CellInfo(unsigned int ip, unsigned int ic, View_t *view=0, double *pos=0, double *dpos=0) const
void FillSliceVector(size_t const &s, art::FindMany< T > &fmp, std::vector< T const * > &ptrMap)
Calibrated quantities relying on position in the orthogonal view. To generate a rb::CellHit from a rb...
Definition: RecoHit.h:19
std::string fSkimmedLabel
label of module creating skimmed file
const XML_Char * s
Definition: expat.h:262
TTree * fCosmicTree
Tree filled directly from cosrej etc. objects (bypasses parameters struct)
Far Detector at Ash River, MN.
novadaq::cnv::DetId fDetId
the id of the detector we are looking at
art::ServiceHandle< geo::Geometry > fGeom
TTree * fNumuTree
Tree to hold values of cut parameters for selected events.
Near Detector in the NuMI cavern.
TTree * fHitTree
Tree to hold information about every hit in a kalman track.
#define LOG_WARNING(category)
Cosmic Rejection PIDs for Numu analysis.
Definition: FillParentInfo.h:9
float X() const
Definition: RecoHit.h:36
float Mag() const
float Y() const
Definition: RecoHit.h:37
#define LOG_VERBATIM(category)
float PECorr() const
Definition: RecoHit.cxx:47
int GoodDiBlockMask(int subrun=-1, bool reload=false)
Float_t e
Definition: plot.C:35
T max(sqlite3 *const db, std::string const &table_name, std::string const &column_name)
Definition: statistics.h:68
NumuCutParameters fNumuCutParams
from the ParametersNumu object
#define W(x)
TTree * fCosmicHitTree
Tree to hold information about every hit in a kalman track.
unsigned int HighestPIDTrack(const std::vector< art::Ptr< rb::Track > > &sliceTracks, const std::string &remidModuleLabel, const art::Event &e)
Definition: ReMId.cxx:249
seed_t art::EngineCreator::get_seed_value ( fhicl::ParameterSet const &  pset,
char const  key[] = "seed",
seed_t const  implicit_seed = -1 
)
inherited
art::Handle<art::TriggerResults> art::EventObserverBase::getTriggerResults ( Event const &  e) const
inlineinherited

Definition at line 61 of file EventObserverBase.h.

References art::detail::CachedProducts::getOneTriggerResults(), and art::EventObserverBase::selectors_.

62  {
64  }
detail::CachedProducts selectors_
art::Handle< art::TriggerResults > getOneTriggerResults(Event const &) const
Float_t e
Definition: plot.C:35
template<typename T , BranchType = InEvent>
ProductToken<T> art::Consumer::mayConsume ( InputTag const &  )
inherited
template<typename T , art::BranchType BT>
art::ProductToken<T> art::Consumer::mayConsume ( InputTag const &  it)
inherited

Definition at line 189 of file Consumer.h.

References art::InputTag::instance(), PandAna.reco_validation.prod5_pid_validation::invalid, art::InputTag::label(), art::InputTag::process(), and T.

190 {
191  if (!moduleContext_)
192  return ProductToken<T>::invalid();
193 
194  consumables_[BT].emplace_back(ConsumableType::Product,
195  TypeID{typeid(T)},
196  it.label(),
197  it.instance(),
198  it.process());
199  return ProductToken<T>{it};
200 }
set< int >::iterator it
static ProductToken< T > invalid()
Definition: ProductToken.h:47
ConsumableProducts consumables_
Definition: Consumer.h:137
double T
Definition: Xdiff_gwt.C:5
bool moduleContext_
Definition: Consumer.h:135
template<typename T , art::BranchType BT>
void art::Consumer::mayConsumeMany ( )
inherited

Definition at line 204 of file Consumer.h.

References T.

205 {
206  if (!moduleContext_)
207  return;
208 
209  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
210 }
ConsumableProducts consumables_
Definition: Consumer.h:137
double T
Definition: Xdiff_gwt.C:5
bool moduleContext_
Definition: Consumer.h:135
template<typename Element , BranchType = InEvent>
ViewToken<Element> art::Consumer::mayConsumeView ( InputTag const &  )
inherited
template<typename T , art::BranchType BT>
art::ViewToken<T> art::Consumer::mayConsumeView ( InputTag const &  it)
inherited

Definition at line 214 of file Consumer.h.

References art::InputTag::instance(), PandAna.reco_validation.prod5_pid_validation::invalid, art::InputTag::label(), art::InputTag::process(), and T.

215 {
216  if (!moduleContext_)
217  return ViewToken<T>::invalid();
218 
219  consumables_[BT].emplace_back(ConsumableType::ViewElement,
220  TypeID{typeid(T)},
221  it.label(),
222  it.instance(),
223  it.process());
224  return ViewToken<T>{it};
225 }
set< int >::iterator it
static ViewToken< Element > invalid()
Definition: ProductToken.h:75
ConsumableProducts consumables_
Definition: Consumer.h:137
double T
Definition: Xdiff_gwt.C:5
bool moduleContext_
Definition: Consumer.h:135
bool art::EventObserverBase::modifiesEvent ( ) const
inlineinherited

Definition at line 25 of file EventObserverBase.h.

26  {
27  return false;
28  }
static cet::exempt_ptr<Consumer> art::Consumer::non_module_context ( )
staticinherited
CosmicBeamComparison& skim::CosmicBeamComparison::operator= ( CosmicBeamComparison const &  )
delete
CosmicBeamComparison& skim::CosmicBeamComparison::operator= ( CosmicBeamComparison &&  )
delete
void art::Consumer::prepareForJob ( fhicl::ParameterSet const &  pset)
protectedinherited
std::string const& art::EventObserverBase::processName ( ) const
inlineinherited
void skim::CosmicBeamComparison::reconfigure ( fhicl::ParameterSet const &  p)

Definition at line 768 of file CosmicBeamComparison_module.cc.

References fCompareNueCuts, fCompareNumuCuts, fCosmicExposureLabel, fNueInstance, fNumuInstance, fPOTSumLabel, fSkimmedLabel, and fhicl::ParameterSet::get().

Referenced by CosmicBeamComparison().

769 {
770 
771  fSkimmedLabel = p.get<std::string >("SkimmedLabel");
772  fCosmicExposureLabel = p.get<std::string >("CosmicExposureLabel", "ifdbspillinfo");
773  fPOTSumLabel = p.get<std::string >("POTSumLabel", "exposure" );
774  fCompareNumuCuts = p.get<bool >("CompareNumuCuts", false);
775  fCompareNueCuts = p.get<bool >("CompareNueCuts", false);
776  auto instanceLabels = p.get<std::vector<std::string> >("Instances");
777 
778  for(auto itr : instanceLabels){
779  if(itr.find("numu") != std::string::npos) fNumuInstance = itr;
780  if(itr.find("nue") != std::string::npos) fNueInstance = itr;
781  }
782 
783  return;
784 }
std::string fNumuInstance
label for numu products
std::string fNueInstance
label for nue products
bool fCompareNueCuts
bool for comparing the nue cuts
std::string fCosmicExposureLabel
label of module creating cosmic exposure
const char * p
Definition: xmltok.h:285
std::string fPOTSumLabel
label of module creating pot sum
std::string fSkimmedLabel
label of module creating skimmed file
bool fCompareNumuCuts
bool for comparing the numu cuts
::xsd::cxx::tree::string< char, simple_type > string
Definition: Database.h:154
void art::EventObserverBase::registerProducts ( MasterProductRegistry ,
ProductDescriptions ,
ModuleDescription const &   
)
inlineinherited

Definition at line 33 of file EventObserverBase.h.

36  {}
fhicl::ParameterSetID art::EventObserverBase::selectorConfig ( ) const
inlineinherited

Definition at line 56 of file EventObserverBase.h.

References art::EventObserverBase::selector_config_id_.

57  {
58  return selector_config_id_;
59  }
fhicl::ParameterSetID selector_config_id_
void art::Consumer::showMissingConsumes ( ) const
protectedinherited

Referenced by art::RootOutput::endJob().

void art::Consumer::validateConsumedProduct ( BranchType const  bt,
ProductInfo const &  pi 
)
protectedinherited
bool art::EventObserverBase::wantAllEvents ( ) const
inlineinherited

Definition at line 46 of file EventObserverBase.h.

References art::EventObserverBase::wantAllEvents_.

Referenced by art::RootOutput::RootOutput().

47  {
48  return wantAllEvents_;
49  }
bool art::EventObserverBase::wantEvent ( Event const &  e)
inlineinherited

Definition at line 51 of file EventObserverBase.h.

References art::EventObserverBase::selectors_, and art::detail::CachedProducts::wantEvent().

52  {
53  return selectors_.wantEvent(e);
54  }
detail::CachedProducts selectors_
Float_t e
Definition: plot.C:35
bool wantEvent(Event const &)
std::string art::EDAnalyzer::workerType ( ) const
inlineinherited

Definition at line 109 of file EDAnalyzer.h.

References art::EDAnalyzer::currentContext().

110  {
111  return "WorkerT<EDAnalyzer>";
112  }

Member Data Documentation

int skim::CosmicBeamComparison::chtCell
private

Definition at line 370 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::chtDir
private

Definition at line 367 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

double skim::CosmicBeamComparison::chtDistFromStart
private

Definition at line 365 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

double skim::CosmicBeamComparison::chtDistFromStop
private

Definition at line 366 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::chtIndex
private

Definition at line 361 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::chtPlane
private

Definition at line 369 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::chtTrackXYZ
private

Definition at line 368 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::chtUncalibTNS
private

Definition at line 362 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

double skim::CosmicBeamComparison::chtUncalibW
private

Definition at line 364 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::chtUncalibXYZ
private

Definition at line 363 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crCellsFromEdge
private

Definition at line 218 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crConCosPID
private

Definition at line 242 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crCosAngle
private

Definition at line 233 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crCosBakCell
private

Definition at line 237 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crCosBakDist
private

Definition at line 241 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crCosChiDiff
private

Definition at line 318 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crCosChiSq
private

Definition at line 316 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crCosDirScore
private

Definition at line 309 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crCosFwdCell
private

Definition at line 236 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crCosFwdDist
private

Definition at line 240 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crCosScoreDiff
private

Definition at line 311 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crCosSlope
private

Definition at line 314 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crCosSpeed
private

Definition at line 308 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crCosTrackCalorimetricEnergy
private

Definition at line 295 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crCosTrackCellsFromEdge
private

Definition at line 280 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crCosTrackClosestEdge
private

Definition at line 279 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::crCosTrackDir
private

Definition at line 277 of file CosmicBeamComparison_module.cc.

Referenced by FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::crCosTrackEarliestHitPos
private

Definition at line 300 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crCosTrackHits
private

Definition at line 278 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::crCosTrackLatestHitPos
private

Definition at line 301 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crCosTrackMaxPlane
private

Definition at line 303 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crCosTrackMaxTNS
private

Definition at line 297 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::crCosTrackMaxXYZ
private

Definition at line 299 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crCosTrackMinPlane
private

Definition at line 302 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crCosTrackMinTNS
private

Definition at line 296 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::crCosTrackMinXYZ
private

Definition at line 298 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crCosTrackNHitsNearEdge0
private

Definition at line 281 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crCosTrackNHitsNearEdge1
private

Definition at line 282 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crCosTrackNHitsNearEdge2
private

Definition at line 283 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crCosTrackNHitsNearEdge3
private

Definition at line 284 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crCosTrackNHitsNearEdge4
private

Definition at line 285 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crCosTrackNHitsNearEdge5
private

Definition at line 286 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crCosTrackNHitsNearEdge6
private

Definition at line 287 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crCosTrackNHitsNearEdge7
private

Definition at line 288 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crCosTrackNHitsNearEdge8
private

Definition at line 289 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crCosTrackNHitsNearEdge9
private

Definition at line 290 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::crCosTrackStart
private

Definition at line 275 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crCosTrackStartTime
private

Definition at line 305 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::crCosTrackStop
private

Definition at line 276 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crCosTrackStopTime
private

Definition at line 306 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crCosTrackTotalADC
private

Definition at line 293 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crCosTrackTotalGeV
private

Definition at line 294 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crCosTrackTotalLength
private

Definition at line 304 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crCosTrackXViewHits
private

Definition at line 291 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crCosTrackYViewHits
private

Definition at line 292 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crdedxllh
private

Definition at line 331 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

int skim::CosmicBeamComparison::crEvent
private

Definition at line 343 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crKalAngle
private

Definition at line 232 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crKalBakCell
private

Definition at line 235 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crKalBakDist
private

Definition at line 239 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crKalChiDiff
private

Definition at line 317 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crKalChiSq
private

Definition at line 315 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::crKalDir
private

Definition at line 272 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crKalDirScore
private

Definition at line 310 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crKalFwdCell
private

Definition at line 234 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crKalFwdDist
private

Definition at line 238 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crKalScoreDiff
private

Definition at line 312 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crKalSlope
private

Definition at line 313 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crKalSpeed
private

Definition at line 307 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crmeasfrac
private

Definition at line 332 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crNCloseLiveDiblocks
private

Definition at line 319 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crNeutrinoE = 0
private

Definition at line 203 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crNTracks3D
private

Definition at line 320 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crPlanesToBack
private

Definition at line 221 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

float skim::CosmicBeamComparison::crPlanesToFront
private

Definition at line 220 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

int skim::CosmicBeamComparison::crQePIDNTracks
private

Definition at line 323 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

float skim::CosmicBeamComparison::crQePIDVal
private

Definition at line 322 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

double skim::CosmicBeamComparison::crReMIdVal
private

Definition at line 204 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crRun
private

Definition at line 341 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crscatllh
private

Definition at line 330 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

float skim::CosmicBeamComparison::crScatt
private

Definition at line 321 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crSliceCalorimetricEnergy
private

Definition at line 213 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crSliceClosestEdge
private

Definition at line 219 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crSliceContigPlanes
private

Definition at line 208 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crSliceHits
private

Definition at line 205 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crSliceMaxPlane
private

Definition at line 209 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crSliceMaxTNS
private

Definition at line 217 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::crSliceMaxXYZ
private

Definition at line 215 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crSliceMinPlane
private

Definition at line 210 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crSliceMinTNS
private

Definition at line 216 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::crSliceMinXYZ
private

Definition at line 214 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crSliceNHitsNearEdge0
private

Definition at line 222 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crSliceNHitsNearEdge1
private

Definition at line 223 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crSliceNHitsNearEdge2
private

Definition at line 224 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crSliceNHitsNearEdge3
private

Definition at line 225 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crSliceNHitsNearEdge4
private

Definition at line 226 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crSliceNHitsNearEdge5
private

Definition at line 227 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crSliceNHitsNearEdge6
private

Definition at line 228 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crSliceNHitsNearEdge7
private

Definition at line 229 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crSliceNHitsNearEdge8
private

Definition at line 230 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crSliceNHitsNearEdge9
private

Definition at line 231 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crSliceTotalADC
private

Definition at line 211 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crSliceTotalGeV
private

Definition at line 212 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crSliceXViewHits
private

Definition at line 206 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crSliceYViewHits
private

Definition at line 207 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crSubrun
private

Definition at line 342 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crTrackCalorimetricEnergy
private

Definition at line 262 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crTrackCellsFromEdge
private

Definition at line 257 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crTrackClosestEdge
private

Definition at line 256 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::crTrackEarliestHitPos
private

Definition at line 267 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crTrackHits
private

Definition at line 245 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::crTrackLatestHitPos
private

Definition at line 268 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crTrackMaxPlane
private

Definition at line 270 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crTrackMaxTNS
private

Definition at line 264 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::crTrackMaxXYZ
private

Definition at line 266 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crTrackMinPlane
private

Definition at line 269 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crTrackMinTNS
private

Definition at line 263 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::crTrackMinXYZ
private

Definition at line 265 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crTrackNHitsNearEdge0
private

Definition at line 246 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crTrackNHitsNearEdge1
private

Definition at line 247 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crTrackNHitsNearEdge2
private

Definition at line 248 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crTrackNHitsNearEdge3
private

Definition at line 249 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crTrackNHitsNearEdge4
private

Definition at line 250 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crTrackNHitsNearEdge5
private

Definition at line 251 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crTrackNHitsNearEdge6
private

Definition at line 252 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crTrackNHitsNearEdge7
private

Definition at line 253 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crTrackNHitsNearEdge8
private

Definition at line 254 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crTrackNHitsNearEdge9
private

Definition at line 255 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::crTrackStart
private

Definition at line 243 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crTrackStartTime
private

Definition at line 273 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::crTrackStop
private

Definition at line 244 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crTrackStopTime
private

Definition at line 274 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crTrackTotalADC
private

Definition at line 260 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crTrackTotalGeV
private

Definition at line 261 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::crTrackTotalLength
private

Definition at line 271 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crTrackXViewHits
private

Definition at line 258 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crTrackYViewHits
private

Definition at line 259 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::crTruthNNu
private

Definition at line 333 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

bool skim::CosmicBeamComparison::crTruthNuIsCC
private

Definition at line 336 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

int skim::CosmicBeamComparison::crTruthNuNHitsSlc
private

Definition at line 334 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

int skim::CosmicBeamComparison::crTruthNuNHitsTot
private

Definition at line 335 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

int skim::CosmicBeamComparison::crTruthNuPDG
private

Definition at line 337 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

int skim::CosmicBeamComparison::crTruthNuPDGOrig
private

Definition at line 338 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

float skim::CosmicBeamComparison::crTruthNuTime
private

Definition at line 340 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

TVector3 skim::CosmicBeamComparison::crTruthNuVertex
private

Definition at line 339 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

float skim::CosmicBeamComparison::crXmax2
private

Definition at line 327 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

float skim::CosmicBeamComparison::crXmin2
private

Definition at line 324 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

float skim::CosmicBeamComparison::crYmax2
private

Definition at line 328 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

float skim::CosmicBeamComparison::crYmin2
private

Definition at line 325 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

float skim::CosmicBeamComparison::crZmax2
private

Definition at line 329 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

float skim::CosmicBeamComparison::crZmin2
private

Definition at line 326 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

TH1F* skim::CosmicBeamComparison::fCalorimetricEnergy
private

Definition at line 174 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fCellsFromEdge
private

Definition at line 123 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

bool skim::CosmicBeamComparison::fCompareNueCuts
private

bool for comparing the nue cuts

Definition at line 96 of file CosmicBeamComparison_module.cc.

Referenced by reconfigure().

bool skim::CosmicBeamComparison::fCompareNumuCuts
private

bool for comparing the numu cuts

Definition at line 95 of file CosmicBeamComparison_module.cc.

Referenced by reconfigure().

TH1F* skim::CosmicBeamComparison::fCosmicExposureByRun
private

Definition at line 116 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and beginSubRun().

std::string skim::CosmicBeamComparison::fCosmicExposureLabel
private

label of module creating cosmic exposure

Definition at line 91 of file CosmicBeamComparison_module.cc.

Referenced by reconfigure().

TTree* skim::CosmicBeamComparison::fCosmicHitTree
private

Tree to hold information about every hit in a kalman track.

Definition at line 105 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TTree* skim::CosmicBeamComparison::fCosmicTree
private

Tree filled directly from cosrej etc. objects (bypasses parameters struct)

Definition at line 103 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fCVN
private

Definition at line 176 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

novadaq::cnv::DetId skim::CosmicBeamComparison::fDetId
private

the id of the detector we are looking at

Definition at line 97 of file CosmicBeamComparison_module.cc.

Referenced by beginRun(), FillNueHistograms(), and FillNumuHistograms().

art::ServiceHandle<ds::DetectorService> skim::CosmicBeamComparison::fDetService
private

detector service

Definition at line 98 of file CosmicBeamComparison_module.cc.

Referenced by beginRun().

EventInfo skim::CosmicBeamComparison::fEventInfo
private

keep track of run and subrun, and event

Definition at line 111 of file CosmicBeamComparison_module.cc.

Referenced by analyze(), and beginJob().

TH1F* skim::CosmicBeamComparison::fFDContPID
private

Definition at line 146 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fFDCosBakCell
private

Definition at line 129 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fFDCosFwdCell
private

Definition at line 128 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fFDCosSpeed
private

Definition at line 137 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH2F* skim::CosmicBeamComparison::fFDCosSpeedVsScore
private

Definition at line 145 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fFDDeltaT
private

Definition at line 138 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fFDDeltaX
private

Definition at line 139 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fFDDeltaY
private

Definition at line 140 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fFDDeltaZ
private

Definition at line 141 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fFDDiffLastFirstLivePlanes
private

Definition at line 148 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

TH2F* skim::CosmicBeamComparison::fFDDTvsDY
private

Definition at line 144 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fFDKalmanAngle
private

Definition at line 132 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fFDKalmanBakCell
private

Definition at line 127 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fFDKalmanDirX
private

Definition at line 133 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fFDKalmanDirY
private

Definition at line 134 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fFDKalmanDirZ
private

Definition at line 135 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fFDKalmanFwdCell
private

Definition at line 126 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fFDKalSpeed
private

Definition at line 136 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fFDPlanesToBack
private

Definition at line 131 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fFDPlanesToFront
private

Definition at line 130 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fFDPt
private

Definition at line 172 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fFDSliceCosAngleShw
private

Definition at line 163 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fFDSliceDistShwVtx
private

Definition at line 164 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fFDSliceFracShwHits
private

Definition at line 165 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fFDSliceHits
private

Definition at line 147 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fFDSliceHitsPerPlane
private

Definition at line 159 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fFDSliceHitsViewAsym
private

Definition at line 162 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fFDSliceHitsXView
private

Definition at line 160 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fFDSliceHitsYView
private

Definition at line 161 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fFDSpeed
private

Definition at line 142 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fFDTrackLength
private

Definition at line 143 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fFDVtxDistToBack
private

Definition at line 171 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fFDVtxDistToBottom
private

Definition at line 169 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fFDVtxDistToEastWall
private

Definition at line 166 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fFDVtxDistToFront
private

Definition at line 170 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fFDVtxDistToTop
private

Definition at line 168 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fFDVtxDistToWestWall
private

Definition at line 167 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

art::ServiceHandle<geo::Geometry> skim::CosmicBeamComparison::fGeom
private

Definition at line 99 of file CosmicBeamComparison_module.cc.

Referenced by FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fGoodPOTByRun
private

Definition at line 115 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and beginSubRun().

TTree* skim::CosmicBeamComparison::fHitTree
private

Tree to hold information about every hit in a kalman track.

Definition at line 104 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fLEM
private

Definition at line 175 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fLID
private

Definition at line 177 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fNDCalEperHit
private

Definition at line 193 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fNDCalorimetricE
private

Definition at line 192 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fNDHadCalEnergySum
private

Definition at line 156 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fNDKalmanBakCell
private

Definition at line 155 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fNDKalmanFwdCell
private

Definition at line 154 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fNDMaxPlaneHit
private

Definition at line 150 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fNDMinPlaneHit
private

Definition at line 151 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fNDNShowers
private

Definition at line 181 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fNDPlanesToFront
private

Definition at line 194 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fNDShowerCalE
private

Definition at line 188 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fNDShowerCalEfraction
private

Definition at line 189 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fNDShowerDistToVtx
private

Definition at line 190 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fNDShowerLength
private

Definition at line 191 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fNDShwEndX
private

Definition at line 185 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fNDShwEndY
private

Definition at line 186 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fNDShwEndZ
private

Definition at line 187 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fNDShwVtxX
private

Definition at line 182 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fNDShwVtxY
private

Definition at line 183 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fNDShwVtxZ
private

Definition at line 184 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fNDSliceHits
private

Definition at line 195 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fNDSliceNHitsPerPlane
private

Definition at line 196 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fNDSliceNoShowerCalE
private

Definition at line 197 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fNDTrackStartZ
private

Definition at line 152 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fNDTrackStopZ
private

Definition at line 153 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fNDVtxX
private

Definition at line 178 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fNDVtxY
private

Definition at line 179 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fNDVtxZ
private

Definition at line 180 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

TH1F* skim::CosmicBeamComparison::fNeutrinoE
private

Definition at line 119 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

NueCutParameters skim::CosmicBeamComparison::fNueCutParams
private

from the ParametersNumu object

Definition at line 108 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

std::string skim::CosmicBeamComparison::fNueInstance
private

label for nue products

Definition at line 94 of file CosmicBeamComparison_module.cc.

Referenced by analyze(), beginJob(), FillNueHistograms(), and reconfigure().

NueCutParameters skim::CosmicBeamComparison::fNueSRParams
private

from the ParametersNumu object

Definition at line 110 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

TTree* skim::CosmicBeamComparison::fNueTree
private

Tree to hold values of cut parameters for selected events.

Definition at line 102 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNueHistograms().

NumuCutParameters skim::CosmicBeamComparison::fNumuCutParams
private

from the ParametersNumu object

Definition at line 107 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

std::string skim::CosmicBeamComparison::fNumuInstance
private

label for numu products

Definition at line 93 of file CosmicBeamComparison_module.cc.

Referenced by analyze(), beginJob(), FillNumuHistograms(), and reconfigure().

NumuCutParameters skim::CosmicBeamComparison::fNumuSRParams
private

from the ParametersNumu object

Definition at line 109 of file CosmicBeamComparison_module.cc.

Referenced by beginJob().

TTree* skim::CosmicBeamComparison::fNumuTree
private

Tree to hold values of cut parameters for selected events.

Definition at line 101 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fPOTByRun
private

Definition at line 114 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and beginSubRun().

std::string skim::CosmicBeamComparison::fPOTSumLabel
private

label of module creating pot sum

Definition at line 92 of file CosmicBeamComparison_module.cc.

Referenced by reconfigure().

TH1F* skim::CosmicBeamComparison::fQePId1Track
private

Definition at line 124 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fQePId2Track
private

Definition at line 125 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fReMId
private

Definition at line 120 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

std::string skim::CosmicBeamComparison::fSkimmedLabel
private

label of module creating skimmed file

Definition at line 90 of file CosmicBeamComparison_module.cc.

Referenced by FillNueHistograms(), FillNumuHistograms(), and reconfigure().

TH1F* skim::CosmicBeamComparison::fSliceContigPlanes
private

Definition at line 122 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TH1F* skim::CosmicBeamComparison::fSliceHits
private

Definition at line 121 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::htCalibPE
private

Definition at line 352 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::htCalibTNS
private

Definition at line 350 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::htCalibXYZ
private

Definition at line 351 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::htCell
private

Definition at line 358 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::htDir
private

Definition at line 355 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

double skim::CosmicBeamComparison::htDistFromStart
private

Definition at line 353 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

double skim::CosmicBeamComparison::htDistFromStop
private

Definition at line 354 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::htIndex
private

Definition at line 346 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::htPlane
private

Definition at line 357 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

int skim::CosmicBeamComparison::htTrackIndex = 0
private

Definition at line 359 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::htTrackXYZ
private

Definition at line 356 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

float skim::CosmicBeamComparison::htUncalibTNS
private

Definition at line 347 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

double skim::CosmicBeamComparison::htUncalibW
private

Definition at line 349 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().

TVector3 skim::CosmicBeamComparison::htUncalibXYZ
private

Definition at line 348 of file CosmicBeamComparison_module.cc.

Referenced by beginJob(), and FillNumuHistograms().


The documentation for this class was generated from the following file: