diff --git a/sbndcode/JobConfigurations/standard/detsim/detector_variations/detsim_detvar_PDSonly_sbnd.fcl b/sbndcode/JobConfigurations/standard/detsim/detector_variations/detsim_detvar_PDSonly_sbnd.fcl index 82cf34122..98221a6fd 100644 --- a/sbndcode/JobConfigurations/standard/detsim/detector_variations/detsim_detvar_PDSonly_sbnd.fcl +++ b/sbndcode/JobConfigurations/standard/detsim/detector_variations/detsim_detvar_PDSonly_sbnd.fcl @@ -9,10 +9,12 @@ producers: { rns: { module_type: "RandomNumberSaver" } opdaq: @local::sbnd_opdetdigitizer +pmtpulseoscillation: @local::sbnd_pmtpulseoscillation +pmtmetricmc: @local::sbnd_pmtmcmetricproducer } #define the producer and filter modules for this path, order matters, -simulate: [rns, opdaq] +simulate: [rns, opdaq, pmtpulseoscillation, pmtmetricmc] #define the output stream, there could be more than one if using filters stream1: [ out1 ] diff --git a/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTGainFluct.fcl b/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTGainFluct.fcl deleted file mode 100644 index 8db6961e9..000000000 --- a/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTGainFluct.fcl +++ /dev/null @@ -1,2 +0,0 @@ -#include "standard_detsim_sbnd.fcl" -physics.producers.opdaq.GainFluctuationsParams.DynodeK: 0.25 \ No newline at end of file diff --git a/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTGainFluct_HighNoise.fcl b/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTGainFluct_HighNoise.fcl deleted file mode 100644 index efa572693..000000000 --- a/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTGainFluct_HighNoise.fcl +++ /dev/null @@ -1,3 +0,0 @@ -#include "standard_detsim_sbnd.fcl" -physics.producers.opdaq.GainFluctuationsParams.DynodeK: 0.25 -physics.producers.opdaq.PMTBaselineRMS: 3.5 \ No newline at end of file diff --git a/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTGainFluct_HighNoise_LowEff.fcl b/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTGainFluct_HighNoise_LowEff.fcl deleted file mode 100644 index a0e579154..000000000 --- a/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTGainFluct_HighNoise_LowEff.fcl +++ /dev/null @@ -1,9 +0,0 @@ -#include "standard_detsim_sbnd.fcl" -physics.producers.opdaq.GainFluctuationsParams.DynodeK: 0.25 -physics.producers.opdaq.PMTBaselineRMS: 3.5 -physics.producers.opdaq.PMTCoatedVUVEff_tpc0: 0.021 -physics.producers.opdaq.PMTCoatedVISEff_tpc0: 0.02329 -physics.producers.opdaq.PMTUncoatedEff_tpc0: 0.0225 -physics.producers.opdaq.PMTCoatedVUVEff_tpc1: 0.018 -physics.producers.opdaq.PMTCoatedVISEff_tpc1: 0.02329 -physics.producers.opdaq.PMTUncoatedEff_tpc1: 0.0225 \ No newline at end of file diff --git a/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTGainFluct_LowEff.fcl b/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTGainFluct_LowEff.fcl deleted file mode 100644 index 85241bc89..000000000 --- a/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTGainFluct_LowEff.fcl +++ /dev/null @@ -1,8 +0,0 @@ -#include "standard_detsim_sbnd.fcl" -physics.producers.opdaq.GainFluctuationsParams.DynodeK: 0.25 -physics.producers.opdaq.PMTCoatedVUVEff_tpc0: 0.021 -physics.producers.opdaq.PMTCoatedVISEff_tpc0: 0.02329 -physics.producers.opdaq.PMTUncoatedEff_tpc0: 0.0225 -physics.producers.opdaq.PMTCoatedVUVEff_tpc1: 0.018 -physics.producers.opdaq.PMTCoatedVISEff_tpc1: 0.02329 -physics.producers.opdaq.PMTUncoatedEff_tpc1: 0.0225 \ No newline at end of file diff --git a/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTHighEff.fcl b/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTHighEff.fcl new file mode 100644 index 000000000..984f8e352 --- /dev/null +++ b/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTHighEff.fcl @@ -0,0 +1,8 @@ +#include "detsim_detvar_PDSonly_sbnd.fcl" +#factor of 1.19 +physics.producers.opdaq.PMTCoatedVUVEff_tpc0: 0.04665 +physics.producers.opdaq.PMTCoatedVISEff_tpc0: 0.03094 +physics.producers.opdaq.PMTUncoatedEff_tpc0: 0.04248 +physics.producers.opdaq.PMTCoatedVUVEff_tpc1: 0.04665 +physics.producers.opdaq.PMTCoatedVISEff_tpc1: 0.03094 +physics.producers.opdaq.PMTUncoatedEff_tpc1: 0.04248 diff --git a/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTHighNoise.fcl b/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTHighNoise.fcl deleted file mode 100644 index 83eb95ee3..000000000 --- a/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTHighNoise.fcl +++ /dev/null @@ -1,2 +0,0 @@ -#include "standard_detsim_sbnd.fcl" -physics.producers.opdaq.PMTBaselineRMS: 3.5 \ No newline at end of file diff --git a/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTHighNoise_LowEff.fcl b/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTHighNoise_LowEff.fcl deleted file mode 100644 index 047e7d9cf..000000000 --- a/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTHighNoise_LowEff.fcl +++ /dev/null @@ -1,8 +0,0 @@ -#include "standard_detsim_sbnd.fcl" -physics.producers.opdaq.PMTBaselineRMS: 3.5 -physics.producers.opdaq.PMTCoatedVUVEff_tpc0: 0.021 -physics.producers.opdaq.PMTCoatedVISEff_tpc0: 0.02329 -physics.producers.opdaq.PMTUncoatedEff_tpc0: 0.0225 -physics.producers.opdaq.PMTCoatedVUVEff_tpc1: 0.018 -physics.producers.opdaq.PMTCoatedVISEff_tpc1: 0.02329 -physics.producers.opdaq.PMTUncoatedEff_tpc1: 0.0225 \ No newline at end of file diff --git a/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTLowEff.fcl b/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTLowEff.fcl index ddb999453..fea4ed43b 100644 --- a/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTLowEff.fcl +++ b/sbndcode/JobConfigurations/standard/detsim/detector_variations/pmt_variations/standard_detsim_sbnd_PMTLowEff.fcl @@ -1,7 +1,8 @@ -#include "standard_detsim_sbnd.fcl" -physics.producers.opdaq.PMTCoatedVUVEff_tpc0: 0.021 -physics.producers.opdaq.PMTCoatedVISEff_tpc0: 0.02329 -physics.producers.opdaq.PMTUncoatedEff_tpc0: 0.0225 -physics.producers.opdaq.PMTCoatedVUVEff_tpc1: 0.018 -physics.producers.opdaq.PMTCoatedVISEff_tpc1: 0.02329 -physics.producers.opdaq.PMTUncoatedEff_tpc1: 0.0225 \ No newline at end of file +#include "detsim_detvar_PDSonly_sbnd.fcl" +#factor of 0.94 +physics.producers.opdaq.PMTCoatedVUVEff_tpc0: 0.03685 +physics.producers.opdaq.PMTCoatedVISEff_tpc0: 0.02444 +physics.producers.opdaq.PMTUncoatedEff_tpc0: 0.03356 +physics.producers.opdaq.PMTCoatedVUVEff_tpc1: 0.03685 +physics.producers.opdaq.PMTCoatedVISEff_tpc1: 0.02444 +physics.producers.opdaq.PMTUncoatedEff_tpc1: 0.03356 diff --git a/sbndcode/JobConfigurations/standard/reco/detector_variations/reco1_detvar_sbnd.fcl b/sbndcode/JobConfigurations/standard/reco/detector_variations/reco1_detvar_sbnd.fcl index 305bcf680..9984a4d1e 100644 --- a/sbndcode/JobConfigurations/standard/reco/detector_variations/reco1_detvar_sbnd.fcl +++ b/sbndcode/JobConfigurations/standard/reco/detector_variations/reco1_detvar_sbnd.fcl @@ -1,3 +1,32 @@ #include "standard_reco1_sbnd.fcl" process_name: Reco1Var + +physics.producers.specialblipgaushit.CalDataModuleLabel: "simtpc2d:dnnsp" +physics.producers.simplemerge: @local::sbnd_simplemerge +physics.producers.simplemerge.InputSourcesLabels: ["largeant","largeantdropped"] + +sbnd_reco1_producer_sequence_DetVar: [ + rns + , simplemerge + , opdecopmt +#, opdecoxarapuca + , ophitpmt +#, ophitxarapuca + , opflashtpc0 + , opflashtpc1 +#, opflashtpc0xarapuca +#, opflashtpc1xarapuca + , gaushit + , gaushitTruthMatch + , specialblipgaushit + , blipgaushitTruthMatch + , crtstrips + , crtclustering + , crtspacepoints + , crttracks + , cluster3d +] + +physics.reco1: @local::sbnd_reco1_producer_sequence_DetVar +physics.analyzers.supera.supera_params: "sbnd_supera_parameters.fcl" diff --git a/sbndcode/JobConfigurations/standard/reco/detector_variations/sbnd_supera_parameters.fcl b/sbndcode/JobConfigurations/standard/reco/detector_variations/sbnd_supera_parameters.fcl new file mode 100644 index 000000000..ad891af97 --- /dev/null +++ b/sbndcode/JobConfigurations/standard/reco/detector_variations/sbnd_supera_parameters.fcl @@ -0,0 +1,183 @@ +#B. Carlson +#bcarlson1@ufl.edu + +ProcessDriver: { + + Verbosity: 2 + EnableFilter: true + RandomAccess: false + ProcessType: ["SuperaMCTruth","SuperaMCTruth","SuperaBBoxInteraction","SuperaMCParticleCluster","SuperaSpacePoint","SuperaOptical","SuperaCRT","Tensor3DFromCluster3D","CombineTensor3D","ThresholdTensor3D","CombineTensor3D","ParticleCorrector","EmptyTensorFilter","RescaleChargeTensor3D"] + ProcessName: ["MultiPartVrtx","MultiPartRain","SuperaBBoxInteraction","SuperaMCParticleCluster","SuperaSpacePoint","SuperaOptical","SuperaCRT","Tensor3DFromCluster3D","CombineTensor3DGhost","ThresholdTensor3D","CombineTensor3D","ParticleCorrector","EmptyTensorFilter","RescaleChargeTensor3D"] + + IOManager: { + Verbosity: 2 + Name: "IOManager" + IOMode: 1 + OutFileName: "out_test.root" + InputFiles: [] + InputDirs: [] + StoreOnlyType: [] + StoreOnlyName: [] + } + + ProcessList: { + EmptyTensorFilter: { + Tensor3DProducerList: ["pcluster_semantics_ghost"] + MinVoxel3DCountList: [1] + } + + RescaleChargeTensor3D: { + HitKeyProducerList: ["reco_hit_key0","reco_hit_key1","reco_hit_key2"] + HitChargeProducerList: ["reco_hit_charge0","reco_hit_charge1","reco_hit_charge2"] + OutputProducer: "reco_rescaled" + ReferenceProducer: "pcluster" + } + + ThresholdTensor3D: { # fill with ghost value (5) + TargetProducer: "pcluster_semantics_ghost" + OutputProducer: "pcluster_semantics_ghost" + PaintValue: 5 + } + + CombineTensor3DGhost: { # Combine voxels + OutputProducer: "pcluster_semantics_ghost" + Tensor3DProducers: ["reco"] + PoolType: 0 + } + + CombineTensor3D: { + Tensor3DProducers: ["pcluster_semantics_ghost","pcluster_semantics"] + OutputProducer: "pcluster_semantics_ghost" + PoolType: 0 + } + + SuperaMCParticleCluster: { + OutputLabel: "pcluster" + LArMCParticleProducer: "simplemerge" # produced when merging from previous stage + LArMCShowerProducer: "mcreco" + LArMCTrackProducer: "mcreco" + DeltaSize: 10 + LArSimEnergyDepositLiteProducer: "sedlite" + Meta3DFromCluster3D: "mcst" + Meta2DFromTensor2D: "" + Verbosity: 2 + UseSimEnergyDeposit: false #AnalyzeSimEnergyDeposit if true, else AnalyzeFirstLastStep - default false + UseSimEnergyDepositLite_for_FirstLastStep: false + UseSimEnergyDepositLite: false + UseSimEnergyDepositPoints: false + UseOrigTrackID: true + CryostatList: [0,0] + TPCList: [0,1] + PlaneList: [] + SemanticPriority: [1,2,0,3,4] # 0-4 for shower track michel delta LE-scattering + CheckParticleValidity: false #set false, but this needs to be checked later + AssertParentTrackID: true + SuperaTrue2RecoVoxel3D: { + UseOrigTrackID: true + DebugMode: true + Profile: true + Verbosity: 2 + Meta3DFromCluster3D: "pcluster" + LArSimChProducer: "simtpc2d simpleSC" + LArSpacePointProducers: ["cluster3d"] + OutputTensor3D: "masked_true" + OutputCluster3D: "masked_true2reco" + TwofoldMatching: true + UseTruePosition: true + HitThresholdNe: 75 # Updated SBND Fall 2025 Production + HitWindowTicks: 10.0 # Updated SBND Fall 2025 Production + HitPeakFinding: false + PostAveraging: true + PostAveragingThreshold_cm: 0.50 # Updated SBND Fall 2025 Production + DumpToCSV: false + RecoChargeRange: [-1000,50000] + VoxelDistanceThreshold: 2.0 # Updated SBND Fall 2025 Production + } + } + + MultiPartVrtx: { + Verbosity: 2 + LArMCTruthProducer: "generator" + OutParticleLabel: "mpv" + Origin: 0 + } + MultiPartRain: { + Verbosity: 2 + LArMCTruthProducer: "corsika" + OutParticleLabel: "mpr" + Origin: 0 + } + + SuperaBBoxInteraction: { + Verbosity: 2 + LArMCTruthProducer: "generator" + #LArSimEnergyDepositProducer: "largeant TPCActive" + LArSimEnergyDepositLiteProducer: "sedlite" + UseSEDLite: true + Origin: 0 + Cluster3DLabels: ["mcst","pcluster","sed","masked_true2reco"] + Tensor3DLabels: ["reco","pcluster_index","masked_true"] + BBoxSize: [614.4, 614.4, 614.4] # 2048 pixels of 0.3 mm + BBoxBottom: [-307.2, -307.2, -57.4] # centered on the detector center with some padding + VoxelSize: [0.3,0.3,0.3] + UseFixedBBox: true + CryostatList: [0,0] + TPCList: [0,1] + } + + SuperaSimEnergyDeposit: { + Verbosity: 2 + #LArSimEnergyDepositProducer: "largeant TPCActive" + LArSimEnergyDepositLiteProducer: "sedlite" + LArMCShowerProducer: "mcreco" + UseSEDLite: true + ParticleProducer: "pcluster" + OutCluster3DLabel: "sed" + StoreLength: false + StoreCharge: false + StorePhoton: false + StoreDiffTime: false + StoreAbsTime: true + StoreDEDX: false + TPCList: [0,1] + CryostatList: [0,0] + } + + ParticleCorrector: { + Verbosity: 2 + Cluster3DProducer: "pcluster_highE" + ParticleProducer: "pcluster" + OutputProducer: "corrected" + VoxelMinValue: -1000 + } + + + Tensor3DFromCluster3D: { + Verbosity: 2 + Cluster3DProducerList: ["pcluster","sed"] + OutputProducerList: ["pcluster","sed"] + PITypeList: [1,1] + FixedPIList: [0.,0.] + } + + SuperaSpacePoint: { + Verbosity: 2 + SpacePointProducers: ["cluster3d"] + OutputLabel: "reco" + DropOutput: ["hit_amp","hit_rms","hit_mult","hit_time"] + StoreWireInfo: true + RecoChargeRange: [-1000, 50000] + } + SuperaOptical: { + OpFlashProducers: ["opflashtpc0","opflashtpc1"] # Add XA products when valid + OpFlashOutputs: ["tpc0","tpc1"] # Add XA products when valid + } + SuperaCRT: { + CRTHitProducers: ["crtspacepoints"] + CRTHitOutputs: ["crthit"] + UseCRTSpacePoint: true + } + + } +} +