diff --git a/ubana/CombinedReco/run_combinedrecotree_overlay_nuwro.fcl b/ubana/CombinedReco/run_combinedrecotree_overlay_nuwro.fcl new file mode 100644 index 000000000..48569f6d9 --- /dev/null +++ b/ubana/CombinedReco/run_combinedrecotree_overlay_nuwro.fcl @@ -0,0 +1,4 @@ +#include "run_combinedrecotree_overlay.fcl" +physics.filters.nuselection.AnalysisTools.default.ignoreMCFlux: true +physics.filters.nuselection.AnalysisTools.default.ignoreGTruth: true +physics.filters.singlephotonana.isTextGen: true diff --git a/ubana/CombinedReco/run_combinedrecotree_overlay_sce_nuwro.fcl b/ubana/CombinedReco/run_combinedrecotree_overlay_sce_nuwro.fcl new file mode 100644 index 000000000..07584eca9 --- /dev/null +++ b/ubana/CombinedReco/run_combinedrecotree_overlay_sce_nuwro.fcl @@ -0,0 +1,4 @@ +#include "run_combinedrecotree_overlay_sce.fcl" +physics.filters.nuselection.AnalysisTools.default.ignoreMCFlux: true +physics.filters.nuselection.AnalysisTools.default.ignoreGTruth: true +physics.filters.singlephotonana.isTextGen: true diff --git a/ubana/CombinedReco/run_combinedrecotree_run1_overlay_numi_nuwro.fcl b/ubana/CombinedReco/run_combinedrecotree_run1_overlay_numi_nuwro.fcl new file mode 100644 index 000000000..2cd49e115 --- /dev/null +++ b/ubana/CombinedReco/run_combinedrecotree_run1_overlay_numi_nuwro.fcl @@ -0,0 +1,4 @@ +#include "run_combinedrecotree_run1_overlay_numi.fcl" +physics.filters.nuselection.AnalysisTools.default.ignoreMCFlux: true +physics.filters.nuselection.AnalysisTools.default.ignoreGTruth: true +physics.filters.singlephotonana.isTextGen: true diff --git a/ubana/CombinedReco/run_combinedrecotree_run1_overlay_numi_sce_nuwro.fcl b/ubana/CombinedReco/run_combinedrecotree_run1_overlay_numi_sce_nuwro.fcl new file mode 100644 index 000000000..73f08685a --- /dev/null +++ b/ubana/CombinedReco/run_combinedrecotree_run1_overlay_numi_sce_nuwro.fcl @@ -0,0 +1,4 @@ +#include "run_combinedrecotree_run1_overlay_numi_sce.fcl" +physics.filters.nuselection.AnalysisTools.default.ignoreMCFlux: true +physics.filters.nuselection.AnalysisTools.default.ignoreGTruth: true +physics.filters.singlephotonana.isTextGen: true diff --git a/ubana/CombinedReco/run_combinedrecotree_run3_altLArPIDWeights_overlay_nuwro.fcl b/ubana/CombinedReco/run_combinedrecotree_run3_altLArPIDWeights_overlay_nuwro.fcl new file mode 100644 index 000000000..4baa0b7a9 --- /dev/null +++ b/ubana/CombinedReco/run_combinedrecotree_run3_altLArPIDWeights_overlay_nuwro.fcl @@ -0,0 +1,4 @@ +#include "run_combinedrecotree_run3_altLArPIDWeights_overlay.fcl" +physics.filters.nuselection.AnalysisTools.default.ignoreMCFlux: true +physics.filters.nuselection.AnalysisTools.default.ignoreGTruth: true +physics.filters.singlephotonana.isTextGen: true diff --git a/ubana/CombinedReco/run_combinedrecotree_run3_altLArPIDWeights_overlay_sce_nuwro.fcl b/ubana/CombinedReco/run_combinedrecotree_run3_altLArPIDWeights_overlay_sce_nuwro.fcl new file mode 100644 index 000000000..dd2a4e2a8 --- /dev/null +++ b/ubana/CombinedReco/run_combinedrecotree_run3_altLArPIDWeights_overlay_sce_nuwro.fcl @@ -0,0 +1,4 @@ +#include "run_combinedrecotree_run3_altLArPIDWeights_overlay_sce.fcl" +physics.filters.nuselection.AnalysisTools.default.ignoreMCFlux: true +physics.filters.nuselection.AnalysisTools.default.ignoreGTruth: true +physics.filters.singlephotonana.isTextGen: true diff --git a/ubana/CombinedReco/run_combinedrecotree_run3_normLArPIDWeights_overlay_nuwro.fcl b/ubana/CombinedReco/run_combinedrecotree_run3_normLArPIDWeights_overlay_nuwro.fcl new file mode 100644 index 000000000..bbdf2b60d --- /dev/null +++ b/ubana/CombinedReco/run_combinedrecotree_run3_normLArPIDWeights_overlay_nuwro.fcl @@ -0,0 +1,4 @@ +#include "run_combinedrecotree_run3_normLArPIDWeights_overlay.fcl" +physics.filters.nuselection.AnalysisTools.default.ignoreMCFlux: true +physics.filters.nuselection.AnalysisTools.default.ignoreGTruth: true +physics.filters.singlephotonana.isTextGen: true diff --git a/ubana/CombinedReco/run_combinedrecotree_run3_normLArPIDWeights_overlay_sce_nuwro.fcl b/ubana/CombinedReco/run_combinedrecotree_run3_normLArPIDWeights_overlay_sce_nuwro.fcl new file mode 100644 index 000000000..c890a3da1 --- /dev/null +++ b/ubana/CombinedReco/run_combinedrecotree_run3_normLArPIDWeights_overlay_sce_nuwro.fcl @@ -0,0 +1,4 @@ +#include "run_combinedrecotree_run3_normLArPIDWeights_overlay_sce.fcl" +physics.filters.nuselection.AnalysisTools.default.ignoreMCFlux: true +physics.filters.nuselection.AnalysisTools.default.ignoreGTruth: true +physics.filters.singlephotonana.isTextGen: true diff --git a/ubana/CombinedReco/run_combinedrecotree_run3_overlay_numi_nuwro.fcl b/ubana/CombinedReco/run_combinedrecotree_run3_overlay_numi_nuwro.fcl new file mode 100644 index 000000000..42ad7f4fe --- /dev/null +++ b/ubana/CombinedReco/run_combinedrecotree_run3_overlay_numi_nuwro.fcl @@ -0,0 +1,4 @@ +#include "run_combinedrecotree_run3_overlay_numi.fcl" +physics.filters.nuselection.AnalysisTools.default.ignoreMCFlux: true +physics.filters.nuselection.AnalysisTools.default.ignoreGTruth: true +physics.filters.singlephotonana.isTextGen: true diff --git a/ubana/CombinedReco/run_combinedrecotree_run3_overlay_numi_sce_nuwro.fcl b/ubana/CombinedReco/run_combinedrecotree_run3_overlay_numi_sce_nuwro.fcl new file mode 100644 index 000000000..5a0349944 --- /dev/null +++ b/ubana/CombinedReco/run_combinedrecotree_run3_overlay_numi_sce_nuwro.fcl @@ -0,0 +1,4 @@ +#include "run_combinedrecotree_run3_overlay_numi_sce.fcl" +physics.filters.nuselection.AnalysisTools.default.ignoreMCFlux: true +physics.filters.nuselection.AnalysisTools.default.ignoreGTruth: true +physics.filters.singlephotonana.isTextGen: true diff --git a/ubana/CombinedReco/run_combinedrecotree_run3_overlay_nuwro.fcl b/ubana/CombinedReco/run_combinedrecotree_run3_overlay_nuwro.fcl new file mode 100644 index 000000000..4af1a85e9 --- /dev/null +++ b/ubana/CombinedReco/run_combinedrecotree_run3_overlay_nuwro.fcl @@ -0,0 +1,4 @@ +#include "run_combinedrecotree_run3_overlay.fcl" +physics.filters.nuselection.AnalysisTools.default.ignoreMCFlux: true +physics.filters.nuselection.AnalysisTools.default.ignoreGTruth: true +physics.filters.singlephotonana.isTextGen: true diff --git a/ubana/CombinedReco/run_combinedrecotree_run3_overlay_sce_nuwro.fcl b/ubana/CombinedReco/run_combinedrecotree_run3_overlay_sce_nuwro.fcl new file mode 100644 index 000000000..29d61b50e --- /dev/null +++ b/ubana/CombinedReco/run_combinedrecotree_run3_overlay_sce_nuwro.fcl @@ -0,0 +1,4 @@ +#include "run_combinedrecotree_run3_overlay_sce.fcl" +physics.filters.nuselection.AnalysisTools.default.ignoreMCFlux: true +physics.filters.nuselection.AnalysisTools.default.ignoreGTruth: true +physics.filters.singlephotonana.isTextGen: true diff --git a/ubana/CombinedReco/run_combinedrecotree_run4_overlay_numi_nuwro.fcl b/ubana/CombinedReco/run_combinedrecotree_run4_overlay_numi_nuwro.fcl new file mode 100644 index 000000000..516fa2d16 --- /dev/null +++ b/ubana/CombinedReco/run_combinedrecotree_run4_overlay_numi_nuwro.fcl @@ -0,0 +1,4 @@ +#include "run_combinedrecotree_run4_overlay_numi.fcl" +physics.filters.nuselection.AnalysisTools.default.ignoreMCFlux: true +physics.filters.nuselection.AnalysisTools.default.ignoreGTruth: true +physics.filters.singlephotonana.isTextGen: true diff --git a/ubana/CombinedReco/run_combinedrecotree_run4_overlay_nuwro.fcl b/ubana/CombinedReco/run_combinedrecotree_run4_overlay_nuwro.fcl new file mode 100644 index 000000000..1676cf3aa --- /dev/null +++ b/ubana/CombinedReco/run_combinedrecotree_run4_overlay_nuwro.fcl @@ -0,0 +1,4 @@ +#include "run_combinedrecotree_run4_overlay.fcl" +physics.filters.nuselection.AnalysisTools.default.ignoreMCFlux: true +physics.filters.nuselection.AnalysisTools.default.ignoreGTruth: true +physics.filters.singlephotonana.isTextGen: true diff --git a/ubana/MicroBooNEWireCell/WireCellAnaTree_module.cc b/ubana/MicroBooNEWireCell/WireCellAnaTree_module.cc index 2b9e35b2c..c0ab929af 100644 --- a/ubana/MicroBooNEWireCell/WireCellAnaTree_module.cc +++ b/ubana/MicroBooNEWireCell/WireCellAnaTree_module.cc @@ -259,6 +259,7 @@ class WireCellAnaTree : public art::EDAnalyzer { bool f_get_spill_time; bool f_get_redk2nu_time; bool f_get_reboone_time; + bool f_no_mcflux; // LArPID and MC backtracking (from LANTERN) options: bool fRunLArPID; //run LArPID network over reco particles @@ -2066,6 +2067,7 @@ void WireCellAnaTree::reconfigure(fhicl::ParameterSet const& pset) f_get_spill_time = pset.get("get_spill_time", false); f_get_reboone_time = pset.get("get_reboone_time", false); f_get_redk2nu_time = pset.get("get_redk2nu_time", false); + f_no_mcflux = pset.get("no_mcflux", false); fTimeBetweenBuckets = pset.get("TimeBetweenBuckets",1e9/53.103e6); fBucketTimeSigma = pset.get("BucketTimeSigma",0.750); fNBucketsPerBatch = pset.get("NBucketsPerBatch",84); @@ -4998,13 +5000,13 @@ void WireCellAnaTree::analyze(art::Event const& e) } }//f_get_reboone_time - else if(f_get_spill_time){ + else if(f_get_spill_time && !f_no_mcflux){ propegation_time = (f_mcflux_dk2gen+f_mcflux_gen2vtx+f_mcflux_vx/100)*100*0.033356; } f_cor_nu_time = propegation_time + decay_time + spill_time; - f_cor_nu_time_nospill = propegation_time + decay_time; - f_cor_nu_time_spill = spill_time; + f_cor_nu_time_nospill = propegation_time + decay_time; + f_cor_nu_time_spill = spill_time; f_cor_nu_deltatime = f_cor_nu_time - f_truth_nu_pos[3]; std::cout<<"Origional nu time "<("makeNuMINtuple", false); fIgnoreMCFlux = p.get("ignoreMCFlux", false); + fIgnoreGTruth = p.get("ignoreGTruth", false); NuMIOpFilterProd = p.get("NuMIOpFiltProcName",""); NuMISWTrigProd = p.get("NuMISWTriggerProcName","" ); } @@ -1470,7 +1472,11 @@ void DefaultAnalysis::SaveTruth(art::Event const &e) { // load GTruth - auto const >_h = e.getValidHandle>(fMCTproducer); + if(!fIgnoreGTruth){ + auto const >_h = e.getValidHandle>(fMCTproducer); + auto gt = gt_h->at(0); + _resid = gt.fResNum; + } // load MCTruth auto const &mct_h = e.getValidHandle>(fMCTproducer); @@ -1529,9 +1535,6 @@ void DefaultAnalysis::SaveTruth(art::Event const &e) }// if flux handle is valid }// if ignore MCFlux is false - auto gt = gt_h->at(0); - _resid = gt.fResNum; - auto mct = mct_h->at(0); if (mct.NeutrinoSet()){ // NeutrinoSet(): whether the neutrino information has been set auto neutrino = mct.GetNeutrino();