From 88f04650039b5ec65bb9f77fa0c10e767c871fca Mon Sep 17 00:00:00 2001 From: Yoann Kermaidic Date: Tue, 16 Dec 2025 17:06:01 +0100 Subject: [PATCH 1/2] PDVD beam data reprocessing --- .../standard_reco_stage0_protodunevd_beam.fcl | 51 +++++++++++++++++++ .../standard_reco_stage1_protodunevd_beam.fcl | 31 +++++++++++ 2 files changed, 82 insertions(+) create mode 100644 fcl/protodunevd/reco/standard_reco_stage0_protodunevd_beam.fcl create mode 100644 fcl/protodunevd/reco/standard_reco_stage1_protodunevd_beam.fcl diff --git a/fcl/protodunevd/reco/standard_reco_stage0_protodunevd_beam.fcl b/fcl/protodunevd/reco/standard_reco_stage0_protodunevd_beam.fcl new file mode 100644 index 00000000..9502f89e --- /dev/null +++ b/fcl/protodunevd/reco/standard_reco_stage0_protodunevd_beam.fcl @@ -0,0 +1,51 @@ +#include "PDVDTriggerReader4.fcl" +#include "PDVDTimingRawDecoder.fcl" +#include "messageservices_dune.fcl" +#include "BeamEvent.fcl" +#include "PDHDCTBRawDecoder.fcl" +#include "standard_reco_protodunevd_keepup.fcl" + +services.IFBeam: {} + +physics.producers.ctbrawdecoder: @local::PDHDCTBRawDecoderDefaults +physics.producers.timingrawdecoder: @local::PDVDTimingRawDecoder +physics.producers.beamevent: @local::pdvd_beamevent + +#Put all of the raw decoders here since they're quick +physics.produce: [ + tpcrawdecoder, + triggerrawdecoder, + timingrawdecoder, + ctbrawdecoder, + beamevent, + pdvddaphne, + ophit, + wclsdatavd, + gaushit, + nhitsfilter, + reco3d, + hitpdune +] + +#don't need a timestamp here -- only intermittent +services.TFileService.fileName: "%ifb_reco_stage0_%tc_beam_hists.root" + +services.DetectorClocksService.DefaultTrigTime: 2500 +services.DetectorClocksService.TriggerOffsetTPC: 2500 +#services.DetectorClocksService.ClockSpeedTPC: 1.953125 +services.DetectorClocksService.ClockSpeedTPC: 2 + +services.DetectorPropertiesService.Efield: [0.450,3.125, 5e-1, 3.125] +services.DetectorPropertiesService.NumberTimeSamples: 10000 +services.DetectorPropertiesService.ReadOutWindowSize: 10000 + +outputs.out1.fileName: "%ifb_reco_stage0_%tc_beam.root" +outputs.out1.outputCommands: [ "keep *"] +process_name: pdvdbeamstage0 + +physics.producers.tpcrawdecoder.CrateList: [-1] +physics.producers.tpcrawdecoder.DecoderToolParams: @local::PDVDDataInterfaceWIBEthDefaults + +physics.producers.beamevent.ForceRead: false +physics.producers.beamevent.MatchTime: true + diff --git a/fcl/protodunevd/reco/standard_reco_stage1_protodunevd_beam.fcl b/fcl/protodunevd/reco/standard_reco_stage1_protodunevd_beam.fcl new file mode 100644 index 00000000..b77ddf9a --- /dev/null +++ b/fcl/protodunevd/reco/standard_reco_stage1_protodunevd_beam.fcl @@ -0,0 +1,31 @@ +#include "standard_reco_protodunevd_keepup.fcl" + +#Put all of the raw decoders here since they're quick +physics.produce: [ + opflash, + wclsdatavd, + gaushit, + nhitsfilter, + reco3d, + hitpdune, + pandora, + pandoraTrack, + pandoraShower, +## pandoraStdcalo, + pandoraGnocalo #, +] + +#don't need a timestamp here -- only intermittent +services.TFileService.fileName: "%ifb_reco_stage0_%tc_beam_hists.root" + +services.DetectorClocksService.DefaultTrigTime: 2500 +services.DetectorClocksService.TriggerOffsetTPC: 2500 +services.DetectorClocksService.ClockSpeedTPC: 2 + +services.DetectorPropertiesService.Efield: [0.450,3.125, 5e-1, 3.125] +services.DetectorPropertiesService.NumberTimeSamples: 10000 +services.DetectorPropertiesService.ReadOutWindowSize: 10000 + +outputs.out1.fileName: "%ifb_reco_stage0_%tc_beam.root" +outputs.out1.outputCommands: [ "keep *"] +process_name: pdvdbeamstage1 From 7a91df179031ac8cd45b91b62c0c68332588c6b3 Mon Sep 17 00:00:00 2001 From: Yoann Kermaidic Date: Tue, 16 Dec 2025 21:13:09 +0100 Subject: [PATCH 2/2] Remove duplicated producers --- .../standard_reco_stage1_protodunevd_beam.fcl | 17 ++++++++--------- 1 file changed, 8 insertions(+), 9 deletions(-) diff --git a/fcl/protodunevd/reco/standard_reco_stage1_protodunevd_beam.fcl b/fcl/protodunevd/reco/standard_reco_stage1_protodunevd_beam.fcl index b77ddf9a..e78c4211 100644 --- a/fcl/protodunevd/reco/standard_reco_stage1_protodunevd_beam.fcl +++ b/fcl/protodunevd/reco/standard_reco_stage1_protodunevd_beam.fcl @@ -1,13 +1,12 @@ #include "standard_reco_protodunevd_keepup.fcl" -#Put all of the raw decoders here since they're quick +source:{ + module_type: RootInput + maxEvents: 100 +} + physics.produce: [ opflash, - wclsdatavd, - gaushit, - nhitsfilter, - reco3d, - hitpdune, pandora, pandoraTrack, pandoraShower, @@ -16,7 +15,7 @@ physics.produce: [ ] #don't need a timestamp here -- only intermittent -services.TFileService.fileName: "%ifb_reco_stage0_%tc_beam_hists.root" +services.TFileService.fileName: "%ifb_reco_stage1_%tc_beam_hists.root" services.DetectorClocksService.DefaultTrigTime: 2500 services.DetectorClocksService.TriggerOffsetTPC: 2500 @@ -26,6 +25,6 @@ services.DetectorPropertiesService.Efield: [0.450,3.125, 5e-1, 3.125] services.DetectorPropertiesService.NumberTimeSamples: 10000 services.DetectorPropertiesService.ReadOutWindowSize: 10000 -outputs.out1.fileName: "%ifb_reco_stage0_%tc_beam.root" -outputs.out1.outputCommands: [ "keep *"] +outputs.out1.fileName: "%ifb_reco_stage1_%tc_beam.root" +outputs.out1.outputCommands: [ "keep *", "drop raw::RawDigits_*_*_*", "drop *_wclsdatavd_wiener_*", "drop *_wclsdatavd_gaus_*"] process_name: pdvdbeamstage1