Skip to content

Smanthey/solarnuana updates #112

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
mantheys wants to merge 7 commits into
base: develop
Choose a base branch
from
Open

Smanthey/solarnuana updates #112

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
7 commits
Select commit Hold shift + click to select a range
8555472
ADD WAVEFORM OUTPUT TO SOLARNUANA
mantheys Nov 18, 2025
bee2f85
Merge branch 'smanthey/solarnuana_updates' of https://github.com/DUNE…
mantheys Nov 18, 2025
d80c5dc
ADD NEW FLASH VARIABLES WITH OPHIT AND WAVEFORM BACKTRACKING
mantheys Nov 19, 2025
23a3315
ADD OPTION TO CONVERT UNIT OF OPHITS IF THEY ARE SAVED IN TICKS
mantheys Nov 21, 2025
62aa852
Merge branch 'smanthey/solar_waveform_analysis' of https://github.com…
mantheys Nov 21, 2025
5e71961
UPDATE DEFAULT
mantheys Nov 21, 2025
a3ea143
Merge branch 'smanthey/solar_waveform_analysis' of https://github.com…
mantheys Nov 21, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
254 changes: 175 additions & 79 deletions duneana/SolarNuAna/SolarNuAna_module.cc
View file
Open in desktop

Large diffs are not rendered by default.

20 changes: 12 additions & 8 deletions duneana/SolarNuAna/fcl/SolarNuAna.fcl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -7,12 +7,14 @@ BEGIN_PROLOG
{
module_type: "SolarNuAna"
GEANT4Label: "largeant" # The label for the process which ran GEANT4
OpWaveformLabel: "opdec" # The label for the process which ran the OpDetWaveform
HitLabel: "hitfd" # String for the process that made the reco hits
ClusterLabel: "planecluster" # The label for the process which ran the "perplane" hit clustering
SolarClusterLabel: "solarcluster" # The label for the process which ran the solar clustering
TrackLabel: "pmtracktc" # The label for the process which ran the PMTrack
OpHitLabel: "ophitspe" # The label for the process which ran the OpHit
OpFlashLabel: "solarflash" # The label for the process which ran the OpFlash
OpHitTime2us: false # If true, the OpHit time is in [ticks] and will be converted to [us].
OpFlashTime2us: false # If true, the OpFlash time is in [ticks] and will be converted to [us]. For other producers than "solarflash".
OpFlashTimeOffset: 18.1 # Time offset to be applied to the OpFlash time in [us] units.

Expand Down Expand Up @@ -53,18 +55,19 @@ BEGIN_PROLOG

GenerateAdjOpFlash: true # Generate OpFlashes.
OpFlashAlgoNHit: 0 # Min number of hits to consider a flash. Change to 3 for bkg run to avoid huge output.
OpFlashAlgoMinTime: 0.008 # Negative time window to look for adj. OpHits in [us] units.
OpFlashAlgoMaxTime: 0.016 # Positive time window to look for adj. OpHits in [us] units.
OpFlashAlgoRad: 300 # Distance to look for adj. OpHits in [cm] units.
OpFlashAlgoPE: 1.5 # PE threshold to look for adj. OpHits.
OpFlashAlgoTriggerPE: 1.5 # PE threshold to trigger an OpFlash.
OpFlashAlgoMinTime: 0.32 # Negative time window to look for adj. OpHits in [us] units.
OpFlashAlgoMaxTime: 0.96 # Positive time window to look for adj. OpHits in [us] units.
OpFlashAlgoWeightedTime: true # If true, use weighted time for OpFlash time calculation, otherwise use flash trigger time.
OpFlashAlgoRad: 600 # Distance to look for adj. OpHits in [cm] units.
OpFlashAlgoPE: 0.0 # PE threshold to look for adj. OpHits.
OpFlashAlgoTriggerPE: 0.0 # PE threshold to trigger an OpFlash.
OpFlashAlgoHotVertexThld: 0.3 # Relative threshold to consider a hit as hot for opflash vertex determination [0-1].
OpFlashAlgoHitDuplicates: true # If true, allow hits to be used in multiple flashes.

AdjOpFlashX: 140. # X distance to search for adj. OpFlashes reconstructed in [cm] units.
AdjOpFlashY: 140. # Y distance to search for adj. OpFlashes reconstructed in [cm] units.
AdjOpFlashZ: 140. # Z distance to search for adj. OpFlashes reconstructed in [cm] units.
AdjOpFlashMinNHitCut: 3 # Cut on the minimum number of OpHits in the OpFlash.
AdjOpFlashMinNHitCut: 0 # Cut on the minimum number of OpHits in the OpFlash.
AdjOpFlashMembraneProjection: true # If true, the OpFlash matching is projected on the membrane planes for VD.
AdjOpFlashEndCapProjection: false # If true, the OpFlash matching is projected on the end cap planes for VD.
AdjOpFlashMaxPERatioCut: 1.00 # Cut on the maximum OpHit PE contribution to the total OpFlash PE.
Expand Down Expand Up @@ -103,6 +106,7 @@ BEGIN_PROLOG

solar_nu_ana_vd: @local::solar_nu_ana_hd_centralAPA
solar_nu_ana_vd.HitLabel: "gaushit"
solar_nu_ana_vd.OpWaveformLabel: "opdigi10ppm"
solar_nu_ana_vd.OpHitLabel: "ophit10ppm"
solar_nu_ana_vd.OpFlashLabel: "solarflash"
solar_nu_ana_vd.OpFlashTimeOffset: 0
Expand All @@ -112,8 +116,8 @@ BEGIN_PROLOG
solar_nu_ana_vd.XAMembraneY: 743.302
solar_nu_ana_vd.XAFinalCapZ: 2188.38
solar_nu_ana_vd.XAStartCapZ: -96.5
solar_nu_ana_vd.OpFlashAlgoMinTime: 0.03
solar_nu_ana_vd.OpFlashAlgoMaxTime: 0.06
solar_nu_ana_vd.OpFlashAlgoMinTime: 0.96
solar_nu_ana_vd.OpFlashAlgoMaxTime: 3.20
solar_nu_ana_vd.OpFlashAlgoRad: 600.
solar_nu_ana_vd.AdjOpFlashMinPELightMap: [
["Amplitude", [1.7578e-7, -4.9589e-4, 1.2519e0]],
Expand Down
3 changes: 3 additions & 0 deletions ...a/SolarNuAna/fcl/dunefd/solar_ana_flash_radiological_decay0_dune10kt_1x2x6_centralAPA.fcl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6,6 +6,9 @@
physics.analyzers.solarnuana.GenerateAdjOpFlash: true
physics.analyzers.solarnuana.SaveOpFlashInfo: true
physics.analyzers.solarnuana.SaveAdjOpFlashInfo: true
physics.analyzers.solarnuana.OpFlashAlgoPE: 1.5
physics.analyzers.solarnuana.OpFlashAlgoTriggerPE: 1.5
physics.analyzers.solarnuana.AdjOpFlashMinNHitCut: 3
physics.analyzers.solarnuana.AdjOpFlashX: 100.
physics.analyzers.solarnuana.AdjOpFlashY: 100.
physics.analyzers.solarnuana.AdjOpFlashZ: 100.
3 changes: 3 additions & 0 deletions ...a/SolarNuAna/fcl/dunefd/solar_ana_flash_radiological_decay0_dune10kt_1x2x6_lateralAPA.fcl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6,6 +6,9 @@
physics.analyzers.solarnuana.GenerateAdjOpFlash: true
physics.analyzers.solarnuana.SaveOpFlashInfo: true
physics.analyzers.solarnuana.SaveAdjOpFlashInfo: true
physics.analyzers.solarnuana.OpFlashAlgoPE: 1.5
physics.analyzers.solarnuana.OpFlashAlgoTriggerPE: 1.5
physics.analyzers.solarnuana.AdjOpFlashMinNHitCut: 3
physics.analyzers.solarnuana.AdjOpFlashX: 100.
physics.analyzers.solarnuana.AdjOpFlashY: 100.
physics.analyzers.solarnuana.AdjOpFlashZ: 100.
10 changes: 3 additions & 7 deletions ...na/fcl/dunefd/solar_ana_marley_adjflash_radiological_decay0_dune10kt_1x2x6_centralAPA.fcl
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,11 +1,7 @@
# Specific configuration for investigating energy depositions of signal events.
# Run SolarNuAna on the output of a standard reco workflow for DUNE FD

#include "solar_ana_radiological_decay0_dune10kt_1x2x6_centralAPA.fcl"
#include "solar_ana_marley_flash_radiological_decay0_dune10kt_1x2x6_centralAPA.fcl"

physics.analyzers.solarnuana.SaveOpFlashInfo: false

physics.analyzers.solarnuana.SaveSignalDaughters: true
physics.analyzers.solarnuana.GenerateAdjOpFlash: true
physics.analyzers.solarnuana.SaveAdjOpFlashInfo: true
physics.analyzers.solarnuana.AdjOpFlashX: 100.
physics.analyzers.solarnuana.AdjOpFlashY: 100.
physics.analyzers.solarnuana.AdjOpFlashZ: 100.
9 changes: 2 additions & 7 deletions ...na/fcl/dunefd/solar_ana_marley_adjflash_radiological_decay0_dune10kt_1x2x6_lateralAPA.fcl
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,11 +1,6 @@
# Specific configuration for investigating energy depositions of signal events.
# Run SolarNuAna on the output of a standard reco workflow for DUNE FD

#include "solar_ana_radiological_decay0_dune10kt_1x2x6_lateralAPA.fcl"
#include "solar_ana_marley_flash_radiological_decay0_dune10kt_1x2x6_lateralAPA.fcl"

physics.analyzers.solarnuana.SaveSignalDaughters: true
physics.analyzers.solarnuana.GenerateAdjOpFlash: true
physics.analyzers.solarnuana.SaveAdjOpFlashInfo: true
physics.analyzers.solarnuana.AdjOpFlashX: 100.
physics.analyzers.solarnuana.AdjOpFlashY: 100.
physics.analyzers.solarnuana.AdjOpFlashZ: 100.
physics.analyzers.solarnuana.SaveOpFlashInfo: false
4 changes: 3 additions & 1 deletion duneana/SolarNuAna/fcl/dunefdvd/solar_ana_flash_dunevd10kt_1x8x14_3view_30deg.fcl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -5,4 +5,6 @@

physics.analyzers.solarnuana.GenerateAdjOpFlash: true
physics.analyzers.solarnuana.SaveOpFlashInfo: true
physics.analyzers.solarnuana.SaveAdjOpFlashInfo: true
physics.analyzers.solarnuana.SaveAdjOpFlashInfo: true
physics.analyzers.solarnuana.OpFlashAlgoPE: 1.5
physics.analyzers.solarnuana.OpFlashAlgoTriggerPE: 1.5
4 changes: 3 additions & 1 deletion ...dunefdvd/solar_ana_flash_radiological_decay0_dunevd10kt_1x8x14_3view_30deg_optimistic.fcl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6,8 +6,10 @@
physics.analyzers.solarnuana.GenerateAdjOpFlash: true
physics.analyzers.solarnuana.SaveOpFlashInfo: true
physics.analyzers.solarnuana.SaveAdjOpFlashInfo: true
physics.analyzers.solarnuana.OpFlashAlgoPE: 1.5
physics.analyzers.solarnuana.OpFlashAlgoTriggerPE: 1.5
physics.analyzers.solarnuana.AdjOpFlashMembraneProjection: false
physics.analyzers.solarnuana.AdjOpFlashMinNHitCut: 6
physics.analyzers.solarnuana.AdjOpFlashMinNHitCut: 3
physics.analyzers.solarnuana.AdjOpFlashX: 100.
physics.analyzers.solarnuana.AdjOpFlashY: 100.
physics.analyzers.solarnuana.AdjOpFlashZ: 100.
2 changes: 2 additions & 0 deletions ...l/dunefdvd/solar_ana_flash_radiological_decay0_dunevd10kt_1x8x14_3view_30deg_shielded.fcl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6,6 +6,8 @@
physics.analyzers.solarnuana.GenerateAdjOpFlash: true
physics.analyzers.solarnuana.SaveOpFlashInfo: true
physics.analyzers.solarnuana.SaveAdjOpFlashInfo: true
physics.analyzers.solarnuana.OpFlashAlgoPE: 1.5
physics.analyzers.solarnuana.OpFlashAlgoTriggerPE: 1.5
physics.analyzers.solarnuana.AdjOpFlashMembraneProjection: false
physics.analyzers.solarnuana.AdjOpFlashMinNHitCut: 3
physics.analyzers.solarnuana.AdjOpFlashX: 100.
Expand Down
6 changes: 6 additions & 0 deletions ...olar_ana_marley_adjflash_yprojected_radiological_decay0_dunevd10kt_1x8x14_3view_30deg.fcl
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,6 @@
# Specific configuration for investigating energy depositions of signal events.
# Run SolarNuAna on the output of a standard reco workflow for DUNE FD-VD

#include "solar_ana_adjflash_yprojected_radiological_decay0_dunevd10kt_1x8x14_3view_30deg.fcl"

physics.analyzers.solarnuana.SaveSignalDaughters: true
6 changes: 6 additions & 0 deletions ...marley_adjflash_yprojected_radiological_decay0_dunevd10kt_1x8x14_3view_30deg_shielded.fcl
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,6 @@
# Specific configuration for investigating energy depositions of signal events.
# Run SolarNuAna on the output of a standard reco workflow for DUNE FD-VD

#include "solar_ana_adjflash_yprojected_radiological_decay0_dunevd10kt_1x8x14_3view_30deg_shielded.fcl"

physics.analyzers.solarnuana.SaveSignalDaughters: true