diff --git a/fcl/caf/cafmakerjob_icarus_detsim2d_noyzsim.fcl b/fcl/caf/cafmakerjob_icarus_detsim2d_noyzsim.fcl new file mode 100644 index 000000000..cbcfaf1e1 --- /dev/null +++ b/fcl/caf/cafmakerjob_icarus_detsim2d_noyzsim.fcl @@ -0,0 +1,11 @@ +#include "cafmakerjob_icarus.fcl" +#include "cafmaker_add_detsim2d_icarus.fcl" + +# Fix for CAFMaker simchannel -- using correct (?) +physics.producers.cafmaker.SimChannelLabel: "daq:simpleSC" + +# Fix for BackTracker simchannel +services.BackTrackerService.BackTracker.SimChannelModuleLabel: "daq:simpleSC" + +# Fix for MCReco simchannel +physics.producers.mcreco.SimChannelLabel: "daq:simpleSC" \ No newline at end of file diff --git a/fcl/detsim/detsim_2d_icarus_refactored.fcl b/fcl/detsim/detsim_2d_icarus_refactored.fcl index 993ce8586..799170d15 100644 --- a/fcl/detsim/detsim_2d_icarus_refactored.fcl +++ b/fcl/detsim/detsim_2d_icarus_refactored.fcl @@ -19,7 +19,7 @@ physics: { crtdaq: @local::icarus_crtsim opdaq: @local::icarus_simpmt - daq: @local::icarus_simwire_wirecell_filtersed + daq: @local::icarus_simwire_wirecell_shifted rns: { module_type: "RandomNumberSaver" } } # producers diff --git a/fcl/detsim/detsim_2d_icarus_refactored_keepSimChannel.fcl b/fcl/detsim/detsim_2d_icarus_refactored_keepSimChannel.fcl new file mode 100644 index 000000000..799170d15 --- /dev/null +++ b/fcl/detsim/detsim_2d_icarus_refactored_keepSimChannel.fcl @@ -0,0 +1,57 @@ +#include "services_icarus_simulation.fcl" +#include "detsimmodules_wirecell_ICARUS.fcl" +#include "opdetsim_pmt_icarus.fcl" +#include "crtsimmodules_icarus.fcl" +#include "trigger_emulation_icarus.fcl" +#include "rootoutput_icarus.fcl" + +process_name: DetSim + +services: { + @table::icarus_detsim_services + #FileCatalogMetadata: @local::art_file_catalog_mc +} # services + +physics: { + + producers: { + @table::icarus_shifting_triggersim.producers + + crtdaq: @local::icarus_crtsim + opdaq: @local::icarus_simpmt + daq: @local::icarus_simwire_wirecell_shifted + + rns: { module_type: "RandomNumberSaver" } + } # producers + + simulate: [ rns, opdaq, @sequence::icarus_shifting_triggersim.path, daq, crtdaq ] + + # define the output stream, there could be more than one if using filters + stream: [ rootoutput ] + +} # physics + +outputs: { + rootoutput: { + @table::icarus_rootoutput + outputCommands: [ + "keep *" + , "drop *_ionization_*_*" + , "drop *_simdrift_*_*" + , "drop *_pdfastsim_*_*" + , "drop raw::OpDetWaveform*_opdaq_*_*" + , "drop sim::SimEnergyDeposits_largeant_*_*" + , "drop sim::SimEnergyDepositLites_sedlite_*_*" + , "drop sim::SimEnergyDeposits_shifted_*_*" + , "drop *_pmtlvdsgatesinit_*_*" + , "drop *_pmtfixedthrinit_*_*" + , "drop *_pmttriggerwindowsinit_*_*" + , "drop *_triggersimgatesinit_*_*" + ] + } +} + +physics.producers.crtdaq.G4ModuleLabel: "shifted" +physics.producers.opdaq.InputModule: "pdfastsim" +physics.producers.shifted.InitAuxDetSimChannelLabel: "genericcrt" +physics.producers.shifted.InitSimPhotonsLabel: "pdfastsim" diff --git a/fcl/detsim/detsim_2d_icarus_refactored_yzsim_allSimChannels.fcl b/fcl/detsim/detsim_2d_icarus_refactored_yzsim_allSimChannels.fcl new file mode 100644 index 000000000..481fde4bb --- /dev/null +++ b/fcl/detsim/detsim_2d_icarus_refactored_yzsim_allSimChannels.fcl @@ -0,0 +1,69 @@ +#include "services_icarus_simulation.fcl" +#include "larg4_services_icarus.fcl" +#include "detsimmodules_wirecell_ICARUS.fcl" +#include "opdetsim_pmt_icarus.fcl" +#include "crtsimmodules_icarus.fcl" +#include "trigger_emulation_icarus.fcl" +#include "rootoutput_icarus.fcl" +#include "mergesimsources_icarus.fcl" +#include "icarus_simedepfilter.fcl" + +process_name: DetSim + +services: { + @table::icarus_detsim_services + @table::icarus_g4_services + @table::icarus_larg4_services + #FileCatalogMetadata: @local::art_file_catalog_mc +} # services + +physics: { + + producers: { + @table::icarus_shifting_triggersim.producers + + crtdaq: @local::icarus_crtsim + opdaq: @local::icarus_simpmt + filtersed: @local::simedepfilter_ind1gap + daq: @local::icarus_simwire_wirecell_yz + merge: @local::icarus_merge_yz_wirecell_simchans + + rns: { module_type: "RandomNumberSaver" } + } # producers + + simulate: [ rns, opdaq, @sequence::icarus_shifting_triggersim.path, filtersed, daq, merge, crtdaq ] + + # define the output stream, there could be more than one if using filters + stream: [ rootoutput ] + +} # physics + +outputs: { + rootoutput: { + @table::icarus_rootoutput + outputCommands: [ + "keep *" + , "drop *_ionization_*_*" + , "drop *_simdrift_*_*" + , "drop *_pdfastsim_*_*" + , "drop raw::OpDetWaveform*_opdaq_*_*" + , "drop sim::SimEnergyDeposits_largeant_*_*" + , "drop sim::SimEnergyDepositLites_sedlite_*_*" + , "drop sim::SimEnergyDeposits_shifted_*_*" + , "drop sim::SimEnergyDepositLites_shifted_*_*" +# , "drop sim::SimChannels_daq_*_*" +# , "drop *_pmtlvdsgatesinit_*_*" +# , "drop *_pmtfixedthrinit_*_*" +# , "drop *_pmttriggerwindowsinit_*_*" +# , "drop *_triggersimgatesinit_*_*" + ] + } +} + +physics.producers.crtdaq.G4ModuleLabel: "shifted" +physics.producers.opdaq.InputModule: "pdfastsim" +physics.producers.shifted.InitAuxDetSimChannelLabel: "genericcrt" +physics.producers.shifted.InitSimPhotonsLabel: "pdfastsim" +physics.producers.filtersed.InitSimEnergyDepositLabel: "shifted" +physics.producers.filtersed.InitSimEnergyDepositLiteLabel: "shifted" + diff --git a/fcl/gen/genie/simulation_genie_icarus_bnb_volDetEnclosure_nue.fcl b/fcl/gen/genie/simulation_genie_icarus_bnb_volDetEnclosure_nue.fcl new file mode 100644 index 000000000..d5e902d4e --- /dev/null +++ b/fcl/gen/genie/simulation_genie_icarus_bnb_volDetEnclosure_nue.fcl @@ -0,0 +1,13 @@ +# File: simulation_genie_icarus_bnb_volDetEnclosure_nue.fcl +# Purpose: generate neutrino interactions from BNB in the whole DetEnclosure volume -- only nue +# Date: November 26, 2025 +# Author: Mattia Sotgia (msotgia@ge.infn.it) +# +# This configuration is based on `simulation_genie_icarus_bnb_volDetEnclosure.fcl`, +# the only difference being that the generator flavour are only nue. +# + +#include "simulation_genie_icarus_bnb_volDetEnclosure.fcl" + +physics.producers.generator.GenFlavors: [ 12, -12 ] + diff --git a/fcl/gen/single/prod_muon_1.5GeV_isotropic_icarus_active.fcl b/fcl/gen/single/prod_muon_1.5GeV_isotropic_icarus_active.fcl new file mode 100644 index 000000000..65a86c311 --- /dev/null +++ b/fcl/gen/single/prod_muon_1.5GeV_isotropic_icarus_active.fcl @@ -0,0 +1,15 @@ +#include "prod_muon_0-1GeV_isotropic_icarus.fcl" + +process_name: MuonGen + +outputs.out1.fileName: "prod_muon_1.5GeV_isotropic_icarus_active_%tc_gen.root" + +physics.producers.generator.X0: [ -220.215 ] +physics.producers.generator.Y0: [ -23.5 ] +physics.producers.generator.Z0: [ 0 ] +physics.producers.generator.SigmaX: [ 148.275 ] +physics.producers.generator.SigmaY: [ 158.46 ] +physics.producers.generator.SigmaZ: [ 894.951 ] + +physics.producers.generator.P0: [ 1.5 ] +physics.producers.generator.SigmaP: [ 0. ] \ No newline at end of file diff --git a/fcl/reco/Stage0/mc/stage0_run2_icarus_noyzsim_mc.fcl b/fcl/reco/Stage0/mc/stage0_run2_icarus_noyzsim_mc.fcl new file mode 100644 index 000000000..6c6674d8b --- /dev/null +++ b/fcl/reco/Stage0/mc/stage0_run2_icarus_noyzsim_mc.fcl @@ -0,0 +1,68 @@ +## Author(s): Riccardo Triozzi (rtriozzi@phd.unipd.it), Mattia Sotgia (msotgia@ge.infn.it) +## Date: November 27, 2025 +## This fhicl file is used to run "stage0" processing specifically for the case where all +## TPC data is included in an artdaq data product from a single instance +## Without any simulation of the YZ distortion +## +#include "stage0_icarus_mc_defs.fcl" +#include "stage0_icarus_driver_common.fcl" + +process_name: MCstage0 + +## Add the MC module to the list of producers +physics.producers: { @table::icarus_stage0_producers + @table::icarus_stage0_mc_producers + } + +## Use the following to run the full defined stage0 set of modules +physics.path: [ @sequence::icarus_stage0_mc_PMT, + MCDecodeTPCROI, + @sequence::icarus_stage0_multiTPC, + simChannelROI, + @sequence::icarus_stage0_mc_crthit, + @sequence::icarus_stage0_mc_crtreco + ] + +## boiler plate... +physics.outana: [ purityinfoana0, purityinfoana1 ] +physics.trigger_paths: [ path ] +physics.end_paths: [ outana, streamROOT ] + +# Drop data products that are no longer needed, but make sure to keep important items! +# For example, we need to drop the RawDigits from the detector simulation stage but want to keep the SimChannel info from WireCell... +outputs.rootOutput.outputCommands: ["keep *_*_*_*", + "drop *_daq*_*_*", + "drop raw::RawDigit*_*_*_*", + "drop raw::OpDetWaveform*_opdaq_*_*", + "drop *_MCDecodeTPCROI_*_*", + "drop *_decon1droi_*_*", + "drop recob::Wire*_roifinder*_*_*", + "keep *_daq_simpleSC*_*"] + +# Set the expected input for ophit +physics.producers.ophit.InputModule: "shifted" + +# Note the default assumption is that our detector simulation input will come from the WireCell 2D drift simulation, a la 'daq' +# If we are running the 1D drift simulation we need to switch to using: +# `physics.producers.MCDecodeTPCROI.FragmentsLabelVec: ["daq3:PHYSCRATEDATATPCWW","daq2:PHYSCRATEDATATPCWE","daq1:PHYSCRATEDATATPCEW","daq0:PHYSCRATEDATATPCEE"]` +# +physics.producers.MCDecodeTPCROI.FragmentsLabelVec: ["daq:TPCWW","daq:TPCWE","daq:TPCEW","daq:TPCEE"] +physics.producers.MCDecodeTPCROI.OutInstanceLabelVec: ["PHYSCRATEDATATPCWW", "PHYSCRATEDATATPCWE", "PHYSCRATEDATATPCEW", "PHYSCRATEDATATPCEE"] +physics.producers.decon1droi.RawDigitLabelVec: ["MCDecodeTPCROI:PHYSCRATEDATATPCWW","MCDecodeTPCROI:PHYSCRATEDATATPCWE","MCDecodeTPCROI:PHYSCRATEDATATPCEW","MCDecodeTPCROI:PHYSCRATEDATATPCEE"] + +physics.producers.simChannelROI.SimChannelLabelVec: ["daq:simpleSC"] +physics.producers.simChannelROI.OutInstanceLabelVec: ["All"] + +physics.producers.decon2droiEE.wcls_main.params.raw_input_label: "MCDecodeTPCROI:PHYSCRATEDATATPCEE" +physics.producers.decon2droiEW.wcls_main.params.raw_input_label: "MCDecodeTPCROI:PHYSCRATEDATATPCEW" +physics.producers.decon2droiWE.wcls_main.params.raw_input_label: "MCDecodeTPCROI:PHYSCRATEDATATPCWE" +physics.producers.decon2droiWW.wcls_main.params.raw_input_label: "MCDecodeTPCROI:PHYSCRATEDATATPCWW" + +## Need overrides for the purity monitor +physics.analyzers.purityinfoana0.SelectEvents: [ path ] +physics.analyzers.purityinfoana1.SelectEvents: [ path ] +physics.producers.purityana0.RawModuleLabel: ["MCDecodeTPCROI:PHYSCRATEDATATPCWW","MCDecodeTPCROI:PHYSCRATEDATATPCWE","MCDecodeTPCROI:PHYSCRATEDATATPCEW","MCDecodeTPCROI:PHYSCRATEDATATPCEE"] +physics.producers.purityana1.RawModuleLabel: ["MCDecodeTPCROI:PHYSCRATEDATATPCWW","MCDecodeTPCROI:PHYSCRATEDATATPCWE","MCDecodeTPCROI:PHYSCRATEDATATPCEW","MCDecodeTPCROI:PHYSCRATEDATATPCEE"] + +# restore legacy G4 labels +physics.producers.mcophit.SimPhotonsProducer: "shifted" \ No newline at end of file diff --git a/fcl/reco/Stage1/mc/stage1_run2_1d_icarus_MC.fcl b/fcl/reco/Stage1/mc/stage1_run2_1d_icarus_MC.fcl index 85b155397..8f5af8791 100644 --- a/fcl/reco/Stage1/mc/stage1_run2_1d_icarus_MC.fcl +++ b/fcl/reco/Stage1/mc/stage1_run2_1d_icarus_MC.fcl @@ -14,4 +14,4 @@ physics.reco: [ ] physics.producers.cluster3DCryoW.Hit3DBuilderAlg.HitFinderTagVec: ["gaushit1dTPCWW", "gaushit1dTPCWE"] -physics.producers.cluster3DCryoE.Hit3DBuilderAlg.HitFinderTagVec: ["gaushit1dTPCEW", "gaushit1dTPCEE"] +physics.producers.cluster3DCryoE.Hit3DBuilderAlg.HitFinderTagVec: ["gaushit1dTPCEW", "gaushit1dTPCEE"] \ No newline at end of file diff --git a/fcl/reco/Stage1/mc/stage1_run2_1d_icarus_noyzsim_MC.fcl b/fcl/reco/Stage1/mc/stage1_run2_1d_icarus_noyzsim_MC.fcl new file mode 100644 index 000000000..b95d6307f --- /dev/null +++ b/fcl/reco/Stage1/mc/stage1_run2_1d_icarus_noyzsim_MC.fcl @@ -0,0 +1,12 @@ + +#include "stage1_run2_1d_icarus_MC.fcl" + +# Tentative fix for mcreco +physics.producers.mcreco.SimChannelLabel: "shifted" + +# Tentative fix for caloskim +physics.analyzers.caloskimW.SimChannelproducer: "daq:simpleSC" +physics.analyzers.caloskimE.SimChannelproducer: "daq:simpleSC" + +# Tentative fix for backtracker +services.BackTrackerService.BackTracker.SimChannelModuleLabel: "daq:simpleSC" \ No newline at end of file diff --git a/fcl/reco/Stage1/mc/stage1_run2_icarus_MC_pandoraMonitoring.fcl b/fcl/reco/Stage1/mc/stage1_run2_icarus_MC_pandoraMonitoring.fcl new file mode 100644 index 000000000..3c7b54be8 --- /dev/null +++ b/fcl/reco/Stage1/mc/stage1_run2_icarus_MC_pandoraMonitoring.fcl @@ -0,0 +1,10 @@ + +#include "stage1_run2_icarus_MC.fcl" + +physics.producers.pandoraGausCryoE.ConfigFile: "PandoraSettings_Master_ICARUS_PandoraMonitoring.xml" +physics.producers.pandoraGausCryoE.EnableMCParticles: true +physics.producers.pandoraGausCryoE.SimChannelModuleLabel: "daq:simpleSC" + +physics.producers.pandoraGausCryoW.ConfigFile: "PandoraSettings_Master_ICARUS_PandoraMonitoring.xml" +physics.producers.pandoraGausCryoW.EnableMCParticles: true +physics.producers.pandoraGausCryoW.SimChannelModuleLabel: "daq:simpleSC" \ No newline at end of file diff --git a/fcl/reco/Stage1/mc/stage1_run2_icarus_newhitfinder_width_MC.fcl b/fcl/reco/Stage1/mc/stage1_run2_icarus_newhitfinder_width_MC.fcl new file mode 100644 index 000000000..5bc3b7c8b --- /dev/null +++ b/fcl/reco/Stage1/mc/stage1_run2_icarus_newhitfinder_width_MC.fcl @@ -0,0 +1,41 @@ +## File: stage1_run2_icarus_newhitfinder_width_MC.fcl +## Author: Mattia Sotgia (msotgia@ge.infn.it) +## +## This FHiCL introduces the re-tuned hit-finder values and uses the hit widths in Pandora algorithms +## Pulse trains are off. +## + + +#include "stage1_run2_icarus_MC.fcl" + + +# Disable pulse trains (setting them all to 1) +physics.producers.gaushit2dTPCEE.LongMaxHits: [1, 1, 1] +physics.producers.gaushit2dTPCEW.LongMaxHits: [1, 1, 1] +physics.producers.gaushit2dTPCWE.LongMaxHits: [1, 1, 1] +physics.producers.gaushit2dTPCWW.LongMaxHits: [1, 1, 1] + +# Allow for more Gaussians +physics.producers.gaushit2dTPCEE.MaxMultiHit: 30 +physics.producers.gaushit2dTPCEW.MaxMultiHit: 30 +physics.producers.gaushit2dTPCWE.MaxMultiHit: 30 +physics.producers.gaushit2dTPCWW.MaxMultiHit: 30 + +# Way looser Chi2 requirement +physics.producers.gaushit2dTPCEE.Chi2NDF: 2000 +physics.producers.gaushit2dTPCEW.Chi2NDF: 2000 +physics.producers.gaushit2dTPCWE.Chi2NDF: 2000 +physics.producers.gaushit2dTPCWW.Chi2NDF: 2000 + +# Re-fitting strategy +# Anecdotally, this did not help +# physics.producers.gaushit2dTPCEE.TryNplus1Fits: true +# physics.producers.gaushit2dTPCEW.TryNplus1Fits: true +# physics.producers.gaushit2dTPCWE.TryNplus1Fits: true +# physics.producers.gaushit2dTPCWW.TryNplus1Fits: true + +# Pandora settings new +physics.producers.pandoraGausCryoE.ConfigFile: "PandoraSettings_Master_ICARUS_NeutrinoHitWidth.xml" +physics.producers.pandoraGausCryoE.UseHitWidths: true +physics.producers.pandoraGausCryoW.ConfigFile: "PandoraSettings_Master_ICARUS_NeutrinoHitWidth.xml" +physics.producers.pandoraGausCryoW.UseHitWidths: true \ No newline at end of file diff --git a/fcl/reco/Stage1/mc/stage1_run2_icarus_newhitfinder_width_intermediatept_MC.fcl b/fcl/reco/Stage1/mc/stage1_run2_icarus_newhitfinder_width_intermediatept_MC.fcl new file mode 100644 index 000000000..3502f9635 --- /dev/null +++ b/fcl/reco/Stage1/mc/stage1_run2_icarus_newhitfinder_width_intermediatept_MC.fcl @@ -0,0 +1,16 @@ +## File: stage1_run2_icarus_newhitfinder_width_intermediatept_MC.fcl +## Author: Mattia Sotgia (msotgia@ge.infn.it) +## +## This FHiCL introduces the re-tuned hit-finder values and uses the hit widths in Pandora algorithms +## Pulse trains are intermediate = 10. +## + + +#include "stage1_run2_icarus_newhitfinder_width_MC.fcl" + + +# Intermediate pulse train settings +physics.producers.gaushit2dTPCEE.LongMaxHits: [10, 10, 10] +physics.producers.gaushit2dTPCEW.LongMaxHits: [10, 10, 10] +physics.producers.gaushit2dTPCWE.LongMaxHits: [10, 10, 10] +physics.producers.gaushit2dTPCWW.LongMaxHits: [10, 10, 10] \ No newline at end of file diff --git a/fcl/reco/Stage1/mc/stage1_run2_icarus_noyzsim_MC.fcl b/fcl/reco/Stage1/mc/stage1_run2_icarus_noyzsim_MC.fcl new file mode 100644 index 000000000..6f7e39ffc --- /dev/null +++ b/fcl/reco/Stage1/mc/stage1_run2_icarus_noyzsim_MC.fcl @@ -0,0 +1,12 @@ + +#include "stage1_run2_icarus_MC.fcl" + +# Tentative fix for mcreco +physics.producers.mcreco.SimChannelLabel: "shifted" + +# Tentative fix for caloskim +physics.analyzers.caloskimW.SimChannelproducer: "daq:simpleSC" +physics.analyzers.caloskimE.SimChannelproducer: "daq:simpleSC" + +# Tentative fix for backtracker +services.BackTrackerService.BackTracker.SimChannelModuleLabel: "daq:simpleSC" \ No newline at end of file diff --git a/fcl/reco/Stage1/mc/stage1_run2_icarus_noyzsim_MC_pandoraMonitoring.fcl b/fcl/reco/Stage1/mc/stage1_run2_icarus_noyzsim_MC_pandoraMonitoring.fcl new file mode 100644 index 000000000..224457e6e --- /dev/null +++ b/fcl/reco/Stage1/mc/stage1_run2_icarus_noyzsim_MC_pandoraMonitoring.fcl @@ -0,0 +1,10 @@ + +#include "stage1_run2_icarus_noyzsim_MC.fcl" + +physics.producers.pandoraGausCryoE.ConfigFile: "PandoraSettings_Master_ICARUS_PandoraMonitoring.xml" +physics.producers.pandoraGausCryoE.EnableMCParticles: true +physics.producers.pandoraGausCryoE.SimChannelModuleLabel: "daq:simpleSC" + +physics.producers.pandoraGausCryoW.ConfigFile: "PandoraSettings_Master_ICARUS_PandoraMonitoring.xml" +physics.producers.pandoraGausCryoW.EnableMCParticles: true +physics.producers.pandoraGausCryoW.SimChannelModuleLabel: "daq:simpleSC" \ No newline at end of file diff --git a/fcl/reco/Stage1/mc/stage1_run2_icarus_noyzsim_newhitfinder_width_MC.fcl b/fcl/reco/Stage1/mc/stage1_run2_icarus_noyzsim_newhitfinder_width_MC.fcl new file mode 100644 index 000000000..7f6db883d --- /dev/null +++ b/fcl/reco/Stage1/mc/stage1_run2_icarus_noyzsim_newhitfinder_width_MC.fcl @@ -0,0 +1,41 @@ +## File: stage1_run2_icarus_newhitfinder_width_MC.fcl +## Author: Mattia Sotgia (msotgia@ge.infn.it) +## +## This FHiCL introduces the re-tuned hit-finder values and uses the hit widths in Pandora algorithms +## Pulse trains are off. +## + + +#include "stage1_run2_icarus_noyzsim_MC.fcl" + + +# Disable pulse trains (setting them all to 1) +physics.producers.gaushit2dTPCEE.LongMaxHits: [1, 1, 1] +physics.producers.gaushit2dTPCEW.LongMaxHits: [1, 1, 1] +physics.producers.gaushit2dTPCWE.LongMaxHits: [1, 1, 1] +physics.producers.gaushit2dTPCWW.LongMaxHits: [1, 1, 1] + +# Allow for more Gaussians +physics.producers.gaushit2dTPCEE.MaxMultiHit: 30 +physics.producers.gaushit2dTPCEW.MaxMultiHit: 30 +physics.producers.gaushit2dTPCWE.MaxMultiHit: 30 +physics.producers.gaushit2dTPCWW.MaxMultiHit: 30 + +# Way looser Chi2 requirement +physics.producers.gaushit2dTPCEE.Chi2NDF: 2000 +physics.producers.gaushit2dTPCEW.Chi2NDF: 2000 +physics.producers.gaushit2dTPCWE.Chi2NDF: 2000 +physics.producers.gaushit2dTPCWW.Chi2NDF: 2000 + +# Re-fitting strategy +# Anecdotally, this did not help +# physics.producers.gaushit2dTPCEE.TryNplus1Fits: true +# physics.producers.gaushit2dTPCEW.TryNplus1Fits: true +# physics.producers.gaushit2dTPCWE.TryNplus1Fits: true +# physics.producers.gaushit2dTPCWW.TryNplus1Fits: true + +# Pandora settings new +physics.producers.pandoraGausCryoE.ConfigFile: "PandoraSettings_Master_ICARUS_NeutrinoHitWidth.xml" +physics.producers.pandoraGausCryoE.UseHitWidths: true +physics.producers.pandoraGausCryoW.ConfigFile: "PandoraSettings_Master_ICARUS_NeutrinoHitWidth.xml" +physics.producers.pandoraGausCryoW.UseHitWidths: true \ No newline at end of file diff --git a/fcl/reco/Stage1/mc/stage1_run2_icarus_noyzsim_newhitfinder_width_intermediatept_MC.fcl b/fcl/reco/Stage1/mc/stage1_run2_icarus_noyzsim_newhitfinder_width_intermediatept_MC.fcl new file mode 100644 index 000000000..093b2dc22 --- /dev/null +++ b/fcl/reco/Stage1/mc/stage1_run2_icarus_noyzsim_newhitfinder_width_intermediatept_MC.fcl @@ -0,0 +1,16 @@ +## File: stage1_run2_icarus_newhitfinder_width_intermediatept_MC.fcl +## Author: Mattia Sotgia (msotgia@ge.infn.it) +## +## This FHiCL introduces the re-tuned hit-finder values and uses the hit widths in Pandora algorithms +## Pulse trains are intermediate = 10. +## + + +#include "stage1_run2_icarus_noyzsim_newhitfinder_width_MC.fcl" + + +# Intermediate pulse train settings +physics.producers.gaushit2dTPCEE.LongMaxHits: [10, 10, 10] +physics.producers.gaushit2dTPCEW.LongMaxHits: [10, 10, 10] +physics.producers.gaushit2dTPCWE.LongMaxHits: [10, 10, 10] +physics.producers.gaushit2dTPCWW.LongMaxHits: [10, 10, 10] \ No newline at end of file diff --git a/fcl/reco/Stage1/mc/stage1_run2_icarus_noyzsim_newhitfinder_width_slicing_MC.fcl b/fcl/reco/Stage1/mc/stage1_run2_icarus_noyzsim_newhitfinder_width_slicing_MC.fcl new file mode 100644 index 000000000..25e44b62f --- /dev/null +++ b/fcl/reco/Stage1/mc/stage1_run2_icarus_noyzsim_newhitfinder_width_slicing_MC.fcl @@ -0,0 +1,13 @@ +## File: stage1_run2_icarus_newhitfinder_width_slicing_MC.fcl +## Author: Mattia Sotgia (msotgia@ge.infn.it) +## +## This FHiCL introduces the re-tuned hit-finder values and uses the hit widths in Pandora algorithms (also in slice creation) +## Pulse trains are off. +## + + +#include "stage1_run2_icarus_noyzsim_newhitfinder_width_MC.fcl" + +# Pandora settings new +physics.producers.pandoraGausCryoE.ConfigFile: "PandoraSettings_Master_ICARUS_NeutrinoHitWidth_SlicingHitWidth.xml" +physics.producers.pandoraGausCryoW.ConfigFile: "PandoraSettings_Master_ICARUS_NeutrinoHitWidth_SlicingHitWidth.xml" \ No newline at end of file diff --git a/fcl/reco/Stage1/mc/stage1_run2_icarus_noyzsim_newhitfinder_width_slicing_intermediatept_MC.fcl b/fcl/reco/Stage1/mc/stage1_run2_icarus_noyzsim_newhitfinder_width_slicing_intermediatept_MC.fcl new file mode 100644 index 000000000..ab4c8f937 --- /dev/null +++ b/fcl/reco/Stage1/mc/stage1_run2_icarus_noyzsim_newhitfinder_width_slicing_intermediatept_MC.fcl @@ -0,0 +1,13 @@ +## File: stage1_run2_icarus_newhitfinder_width_slicing_intermediatept_MC.fcl +## Author: Mattia Sotgia (msotgia@ge.infn.it) +## +## This FHiCL introduces the re-tuned hit-finder values and uses the hit widths in Pandora algorithms (also in slice creation) +## Pulse trains are intermediate = 10. +## + + +#include "stage1_run2_icarus_noyzsim_newhitfinder_width_intermediatept_MC.fcl" + +# Pandora settings new +physics.producers.pandoraGausCryoE.ConfigFile: "PandoraSettings_Master_ICARUS_NeutrinoHitWidth_SlicingHitWidth.xml" +physics.producers.pandoraGausCryoW.ConfigFile: "PandoraSettings_Master_ICARUS_NeutrinoHitWidth_SlicingHitWidth.xml" \ No newline at end of file diff --git a/icaruscode/TPC/Tracking/ICARUSPandora/scripts/PandoraSettings_Master_ICARUS_NeutrinoHitWidth.xml b/icaruscode/TPC/Tracking/ICARUSPandora/scripts/PandoraSettings_Master_ICARUS_NeutrinoHitWidth.xml new file mode 100644 index 000000000..d7c645a93 --- /dev/null +++ b/icaruscode/TPC/Tracking/ICARUSPandora/scripts/PandoraSettings_Master_ICARUS_NeutrinoHitWidth.xml @@ -0,0 +1,61 @@ + + + false + false + true + + + 1000. + + + + + CaloHitListU + CaloHitListV + CaloHitListW + CaloHitList2D + CaloHitList2D + + + CaloHitListU CaloHitListV CaloHitListW + true + + + + PandoraSettings_Cosmic_Standard.xml + PandoraSettings_Neutrino_ICARUS_HitWidth.xml + PandoraSettings_Slicing_Standard.xml + + true + false + + + + 2.0 + cautious + + + + + PandoraBdt_v08_33_00_SBND.xml + NeutrinoId + 0 + 999 + true + + + Input + Input + false + RecreatedPfos + RecreatedClusters + RecreatedVertices + false + 2.0 + + + + true + true + + diff --git a/icaruscode/TPC/Tracking/ICARUSPandora/scripts/PandoraSettings_Master_ICARUS_NeutrinoHitWidth_SlicingHitWidth.xml b/icaruscode/TPC/Tracking/ICARUSPandora/scripts/PandoraSettings_Master_ICARUS_NeutrinoHitWidth_SlicingHitWidth.xml new file mode 100644 index 000000000..6eb8faeee --- /dev/null +++ b/icaruscode/TPC/Tracking/ICARUSPandora/scripts/PandoraSettings_Master_ICARUS_NeutrinoHitWidth_SlicingHitWidth.xml @@ -0,0 +1,61 @@ + + + false + false + true + + + 1000. + + + + + CaloHitListU + CaloHitListV + CaloHitListW + CaloHitList2D + CaloHitList2D + + + CaloHitListU CaloHitListV CaloHitListW + true + + + + PandoraSettings_Cosmic_Standard.xml + PandoraSettings_Neutrino_ICARUS_HitWidth.xml + PandoraSettings_Slicing_HitWidth.xml + + true + false + + + + 2.0 + cautious + + + + + PandoraBdt_v08_33_00_SBND.xml + NeutrinoId + 0 + 999 + true + + + Input + Input + false + RecreatedPfos + RecreatedClusters + RecreatedVertices + false + 2.0 + + + + true + true + + diff --git a/icaruscode/TPC/Tracking/ICARUSPandora/scripts/PandoraSettings_Master_ICARUS_PandoraMonitoring.xml b/icaruscode/TPC/Tracking/ICARUSPandora/scripts/PandoraSettings_Master_ICARUS_PandoraMonitoring.xml new file mode 100644 index 000000000..bd2171da3 --- /dev/null +++ b/icaruscode/TPC/Tracking/ICARUSPandora/scripts/PandoraSettings_Master_ICARUS_PandoraMonitoring.xml @@ -0,0 +1,61 @@ + + + true + true + true + + + 1000. + + + + + CaloHitListU + CaloHitListV + CaloHitListW + CaloHitList2D + CaloHitList2D + + + CaloHitListU CaloHitListV CaloHitListW + false + + + + PandoraSettings_Cosmic_Standard.xml + PandoraSettings_Neutrino_ICARUS_PandoraMonitoring.xml + PandoraSettings_Slicing_Standard.xml + + true + false + + + + 2.0 + cautious + + + + + PandoraBdt_v08_33_00_SBND.xml + NeutrinoId + 0 + 999 + true + + + Input + Input + false + RecreatedPfos + RecreatedClusters + RecreatedVertices + false + 2.0 + + + + true + true + + diff --git a/icaruscode/TPC/Tracking/ICARUSPandora/scripts/PandoraSettings_Neutrino_ICARUS_HitWidth.xml b/icaruscode/TPC/Tracking/ICARUSPandora/scripts/PandoraSettings_Neutrino_ICARUS_HitWidth.xml new file mode 100644 index 000000000..d56452c93 --- /dev/null +++ b/icaruscode/TPC/Tracking/ICARUSPandora/scripts/PandoraSettings_Neutrino_ICARUS_HitWidth.xml @@ -0,0 +1,374 @@ + + + + true + false + true + + + + CaloHitListU + CaloHitListV + CaloHitListW + CaloHitList2D + CaloHitList2D + + + + + + CaloHitListU + ClustersU + true + true + + + + + + + + + + + + + + + + + + + CaloHitListV + ClustersV + true + true + + + + + + + + + + + + + + + + + + + CaloHitListW + ClustersW + true + true + + + + + + + + + + + + + + + + + + + ClustersU ClustersV ClustersW + + + ClustersU ClustersV ClustersW + CandidateVertices3D + true + false + + + CaloHitListU CaloHitListV CaloHitListW + ClustersU ClustersV ClustersW + NeutrinoVertices3D + true + PandoraBdt_v09_14_00_ICARUS_vertex_bnb.xml + VertexBDTRegion + VertexBDTVertex + + + + + + + + + + ClustersU ClustersV ClustersW + true + + + ClustersU ClustersV ClustersW + + + + + ClustersU + ClustersV + ClustersW + TrackParticles3D + + + + true + true + false + false + + + 0.750.75 + 0.750.75 + + + + + ClustersU + ClustersV + ClustersW + TrackParticles3D + + + + + + + 5. + ClustersU + ClustersV + ClustersW + TrackParticles3D + + + + + + + + + TrackParticles3D + ShowerParticles3D + false + + + ShowerParticles3D + + + ClustersU ClustersV ClustersW + true + + + ClustersU ClustersV ClustersW + + + ClustersU + ClustersV + ClustersW + ShowerParticles3D + + + + + + + + + + ClustersU + ClustersV + ClustersW + TrackParticles3D + + + + true + true + false + false + + + 0.750.75 + 0.750.75 + + + + + ClustersU + ClustersV + ClustersW + TrackParticles3D + + + + + + + 5. + ClustersU + ClustersV + ClustersW + TrackParticles3D + + + + + + + + + ClustersU ClustersV ClustersW + TrackParticles3D + + + ClustersU ClustersV ClustersW + TrackParticles3D + + + ClustersU ClustersV ClustersW + TrackParticles3D + true + 0.5 + 5 + 1. + + + + + ShowerParticles3D + ClustersU ClustersV ClustersW + + + ShowerParticles3D + ClustersU ClustersV ClustersW + + + ShowerParticles3D + ClustersU ClustersV ClustersW + + + + + TrackParticles3D + ShowerParticles3D + true + false + + + TrackParticles3D + TrackCaloHits3D + TrackClusters3D + + 3 + 3 + 3 + 3 + 2 + 2 + 2 + + + + ShowerParticles3D + ShowerCaloHits3D + ShowerClusters3D + + + + + + + + + + TrackParticles3D ShowerParticles3D + ClustersU ClustersV ClustersW TrackClusters3D ShowerClusters3D + + + TrackParticles3D ShowerParticles3D + ClustersU ClustersV ClustersW + + + TrackParticles3D ShowerParticles3D + ClustersU ClustersV ClustersW + true + + + + + NeutrinoVertices3D + NeutrinoParticles3D + + + NeutrinoParticles3D + TrackParticles3D ShowerParticles3D + false + + + + + + + + NeutrinoParticles3D + DaughterVertices3D + + + false + CaloHitList2D + Input + BNBPfoCharTraining + TrackParticles3D + ShowerParticles3D + true + true + PandoraBdt_v09_72_01_ICARUS_pfochar_bnb.xml + PfoCharBDT + PandoraBdt_v09_72_01_ICARUS_pfochar_bnb.xml + PfoCharBDTNoChargeInfo + + + + + + + + + + + + + + + + + NeutrinoParticles3D + + + + + TrackParticles3D + DaughterVertices3D + + + + + NeutrinoParticles3D TrackParticles3D ShowerParticles3D + NeutrinoVertices3D DaughterVertices3D CandidateVertices3D + ClustersU ClustersV ClustersW TrackClusters3D ShowerClusters3D + CaloHitListU CaloHitListV CaloHitListW CaloHitList2D + NeutrinoParticles3D + + diff --git a/icaruscode/TPC/Tracking/ICARUSPandora/scripts/PandoraSettings_Neutrino_ICARUS_PandoraMonitoring.xml b/icaruscode/TPC/Tracking/ICARUSPandora/scripts/PandoraSettings_Neutrino_ICARUS_PandoraMonitoring.xml new file mode 100644 index 000000000..12b9cbf38 --- /dev/null +++ b/icaruscode/TPC/Tracking/ICARUSPandora/scripts/PandoraSettings_Neutrino_ICARUS_PandoraMonitoring.xml @@ -0,0 +1,386 @@ + + + + true + false + true + + + + CaloHitListU + CaloHitListV + CaloHitListW + CaloHitList2D + CaloHitList2D + + + + CaloHitListU CaloHitListV CaloHitListW + false + + + + + + CaloHitListU + ClustersU + true + true + + + + + + + + + + + + + + + + + + CaloHitListV + ClustersV + true + true + + + + + + + + + + + + + + + + + + CaloHitListW + ClustersW + true + true + + + + + + + + + + + + + + + + + ClustersU ClustersV ClustersW + false + + + + + ClustersU ClustersV ClustersW + + + ClustersU ClustersV ClustersW + CandidateVertices3D + true + false + + + CaloHitListU CaloHitListV CaloHitListW + ClustersU ClustersV ClustersW + NeutrinoVertices3D + true + PandoraBdt_v09_14_00_ICARUS_vertex_bnb.xml + VertexBDTRegion + VertexBDTVertex + + + + + + + + + + ClustersU ClustersV ClustersW + true + + + ClustersU ClustersV ClustersW + + + + + ClustersU + ClustersV + ClustersW + TrackParticles3D + + + + true + true + false + false + + + 0.750.75 + 0.750.75 + + + + + ClustersU + ClustersV + ClustersW + TrackParticles3D + + + + + + + 5. + ClustersU + ClustersV + ClustersW + TrackParticles3D + + + + + + + + + TrackParticles3D + ShowerParticles3D + false + + + ShowerParticles3D + + + ClustersU ClustersV ClustersW + true + + + ClustersU ClustersV ClustersW + + + ClustersU + ClustersV + ClustersW + ShowerParticles3D + + + + + + + + + + ClustersU + ClustersV + ClustersW + TrackParticles3D + + + + true + true + false + false + + + 0.750.75 + 0.750.75 + + + + + ClustersU + ClustersV + ClustersW + TrackParticles3D + + + + + + + 5. + ClustersU + ClustersV + ClustersW + TrackParticles3D + + + + + + + + + ClustersU ClustersV ClustersW + TrackParticles3D + + + ClustersU ClustersV ClustersW + TrackParticles3D + + + ClustersU ClustersV ClustersW + TrackParticles3D + true + 0.5 + 5 + 1. + + + + + ShowerParticles3D + ClustersU ClustersV ClustersW + + + ShowerParticles3D + ClustersU ClustersV ClustersW + + + ShowerParticles3D + ClustersU ClustersV ClustersW + + + + ClustersU ClustersV ClustersW + false + + + + + TrackParticles3D + ShowerParticles3D + true + false + + + TrackParticles3D + TrackCaloHits3D + TrackClusters3D + + 3 + 3 + 3 + 3 + 2 + 2 + 2 + + + + ShowerParticles3D + ShowerCaloHits3D + ShowerClusters3D + + + + + + + + + + TrackParticles3D ShowerParticles3D + ClustersU ClustersV ClustersW TrackClusters3D ShowerClusters3D + + + TrackParticles3D ShowerParticles3D + ClustersU ClustersV ClustersW + + + TrackParticles3D ShowerParticles3D + ClustersU ClustersV ClustersW + true + + + + + NeutrinoVertices3D + NeutrinoParticles3D + + + NeutrinoParticles3D + TrackParticles3D ShowerParticles3D + false + + + + + + + + NeutrinoParticles3D + DaughterVertices3D + + + false + CaloHitList2D + Input + BNBPfoCharTraining + TrackParticles3D + ShowerParticles3D + true + true + PandoraBdt_v09_72_01_ICARUS_pfochar_bnb.xml + PfoCharBDT + PandoraBdt_v09_72_01_ICARUS_pfochar_bnb.xml + PfoCharBDTNoChargeInfo + + + + + + + + + + + + + + + + + NeutrinoParticles3D + + + + + TrackParticles3D + DaughterVertices3D + + + + + NeutrinoParticles3D TrackParticles3D ShowerParticles3D + NeutrinoVertices3D DaughterVertices3D CandidateVertices3D + ClustersU ClustersV ClustersW TrackClusters3D ShowerClusters3D + CaloHitListU CaloHitListV CaloHitListW CaloHitList2D + NeutrinoParticles3D + + diff --git a/icaruscode/TPC/Tracking/ICARUSPandora/scripts/PandoraSettings_Slicing_HitWidth.xml b/icaruscode/TPC/Tracking/ICARUSPandora/scripts/PandoraSettings_Slicing_HitWidth.xml new file mode 100644 index 000000000..1fcd721c7 --- /dev/null +++ b/icaruscode/TPC/Tracking/ICARUSPandora/scripts/PandoraSettings_Slicing_HitWidth.xml @@ -0,0 +1,267 @@ + + + true + false + true + + + + CaloHitListU + CaloHitListV + CaloHitListW + CaloHitList2D + CaloHitList2D + + + + + + CaloHitListU + ClustersU + true + true + + + + + + + + + + + + + + + + + + + CaloHitListV + ClustersV + true + true + + + + + + + + + + + + + + + + + + + CaloHitListW + ClustersW + true + true + + + + + + + + + + + + + + + + + + + ClustersU + ClustersV + ClustersW + TrackParticles3D + + + + true + true + false + false + + + 0.750.75 + 0.750.75 + + + + + ClustersU + ClustersV + ClustersW + TrackParticles3D + + + + + + + 5. + ClustersU + ClustersV + ClustersW + TrackParticles3D + + + + + + + + + TrackParticles3D + ShowerParticles3D + false + + + ShowerParticles3D + + + ClustersU ClustersV ClustersW + true + + + ClustersU ClustersV ClustersW + + + ClustersU + ClustersV + ClustersW + ShowerParticles3D + + + + + + + + + + ClustersU + ClustersV + ClustersW + TrackParticles3D + + + + true + true + false + false + + + 0.750.75 + 0.750.75 + + + + + ClustersU + ClustersV + ClustersW + TrackParticles3D + + + + + + + 5. + ClustersU + ClustersV + ClustersW + TrackParticles3D + + + + + + + + + ClustersU ClustersV ClustersW + TrackParticles3D + + + ClustersU ClustersV ClustersW + TrackParticles3D + + + ClustersU ClustersV ClustersW + TrackParticles3D + true + 0.5 + 5 + 1. + + + + + TrackParticles3D + ShowerParticles3D + false + true + + + TrackParticles3D + TrackCaloHits3D + TrackClusters3D + + 3 + 3 + 3 + 3 + 2 + 2 + 2 + + + + ShowerParticles3D + ShowerCaloHits3D + ShowerClusters3D + + + + + + + + + + CaloHitListU + CaloHitListV + CaloHitListW + ClustersU + ClustersV + ClustersW + SliceClusters + SliceParticles + + TrackParticles3D + ShowerParticles3D + + + TrackParticles3D ShowerParticles3D + ClustersU ClustersV ClustersW TrackClusters3D ShowerClusters3D + + + + SliceParticles + +