Skip to content

Commit f5903f8

Browse files
authored
Merge pull request #98 from DUNE/mantheys/duneana_vd
SolarNuAna_module VD Flash Match Studies
2 parents 6ca410b + cfe1ada commit f5903f8

9 files changed

+201
-80
lines changed

duneana/SolarNuAna/SolarNuAna_module.cc

Lines changed: 109 additions & 27 deletions
Large diffs are not rendered by default.

duneana/SolarNuAna/fcl/SolarNuAna.fcl

Lines changed: 4 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -23,11 +23,11 @@ BEGIN_PROLOG
2323
DetectorSizeY: 600 # 1 TPC Y length in [cm].
2424
DetectorSizeZ: 1400 # 1 TPC Z length in [cm].
2525
DetectorDriftTime: 4500 # Don't change this value for HD unless you know what you are doing!!!
26+
ClusterChargeVariable: "Integral" # Charge variable to use for cluster matching: "Integral" or "SummedADC".
2627
XACathodeX: 0 # X position of the VD cathode XAs in [cm].
2728
XAMembraneY: 0 # Y position of the VD membrane XAs in [cm].
2829
XAFinalCapZ: 0 # Z position of the VD final cap XAs in [cm].
2930
XAStartCapZ: 0 # Z position of the VD start cap XAs in [cm].
30-
3131
ClusterAlgoTime: 25 # Time window to look for plane clusters in [tick] units.
3232
ClusterAlgoAdjChannel: 3 # Number of adjacent channels to look for plane clusters.
3333

@@ -57,6 +57,7 @@ BEGIN_PROLOG
5757
# OpFlashAlgoCentroid: false # Use the centroid computation of the flash based on a likelihood estimation.
5858

5959
AdjOpFlashTime: 4500 # Time window to look for adjacent flashes in [tick] units should be one full drift length.
60+
AdjOpFlashX: 100. # X distance to search for adj. OpFlashes reconstructed in [cm] units.
6061
AdjOpFlashY: 100. # Y distance to search for adj. OpFlashes reconstructed in [cm] units.
6162
AdjOpFlashZ: 100. # Z distance to search for adj. OpFlashes reconstructed in [cm] units.
6263
AdjOpFlashMaxPERatioCut: 1.00 # Cut on the maximum OpHit PE contribution to the total OpFlash PE.
@@ -91,10 +92,10 @@ BEGIN_PROLOG
9192
solar_nu_ana_vd.BackgroundLabelVector: @local::generator_dunevd10kt_1x8x14_3view_30deg
9293

9394
solar_nu_ana_vd.Geometry: "VD"
94-
solar_nu_ana_vd.DetectorSizeX: 340
95+
solar_nu_ana_vd.DetectorSizeX: 330
9596
solar_nu_ana_vd.DetectorSizeY: 680
9697
solar_nu_ana_vd.DetectorSizeZ: 2100
97-
solar_nu_ana_vd.DetectorDriftTIme: 8500
98+
solar_nu_ana_vd.DetectorDriftTime: 8500
9899
solar_nu_ana_vd.XACathodeX: -327.5
99100
solar_nu_ana_vd.XAMembraneY: 743.302
100101
solar_nu_ana_vd.XAFinalCapZ: 2188.38
Lines changed: 53 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,53 @@
1+
# solar_ana_dunevd10kt_1x8x14_3view_30deg.fcl
2+
#
3+
# Run SolarNuAna on the output of a standard reco workflow for DUNE FD
4+
5+
#include "services_dune.fcl"
6+
#include "SolarNuAna.fcl"
7+
8+
process_name: SolarNuAna
9+
10+
services:
11+
{
12+
@table::dunefd_services
13+
TFileService: { fileName: "solar_ana_dunevd10kt_1x8x14_3view_30deg_hist.root" }
14+
TimeTracker: {}
15+
MemoryTracker: {} # default is one
16+
RandomNumberGenerator: {} # ART native random number generator
17+
FileCatalogMetadata: @local::art_file_catalog_mc
18+
@table::dunefdvd_1x8x14_3view_30deg_services
19+
# @table::dunefdvd_1x8x6_3view_simulation_services
20+
}
21+
22+
# a temporary fix until the photon back tracker gets fixed...
23+
#services.PhotonBackTrackerService: @erase
24+
25+
physics:
26+
{
27+
analyzers:
28+
{
29+
solarnuana: @local::solar_nu_ana_vd
30+
}
31+
ana: [ solarnuana ]
32+
end_paths: [ ana ]
33+
}
34+
35+
source:
36+
{
37+
module_type: RootInput
38+
maxEvents: -1 # Number of events to create
39+
}
40+
41+
outputs:
42+
{
43+
out1:
44+
{
45+
module_type: RootOutput
46+
fileName: "%ifb_solar.root"
47+
dataTier: "analized"
48+
compressionLevel: 1
49+
}
50+
}
51+
52+
services.message.destinations.LogStandardOut.threshold: "INFO"
53+
services.message.destinations.LogStandardOut.type: "cout"
Lines changed: 7 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,7 @@
1+
# Specific configuration for signal events with coustom flash generation
2+
# Run SolarNuAna on the output of a standard reco workflow for DUNE FD-VD
3+
4+
#include "solar_ana_dunevd10kt_1x8x14_3view_30deg.fcl"
5+
6+
physics.analyzers.solarnuana.GenerateAdjOpFlash: true
7+
physics.analyzers.solarnuana.SaveOpFlashInfo: true
Lines changed: 6 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,6 @@
1+
# Specific configuration for investigating energy depositions of signal events.
2+
# Run SolarNuAna on the output of a standard reco workflow for DUNE FD
3+
4+
#include "solar_ana_dunevd10kt_1x8x14_3view_30deg.fcl"
5+
6+
physics.analyzers.solarnuana.SaveSignalDaughters: true
Lines changed: 6 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,6 @@
1+
# Specific configuration for investigating energy depositions of signal events.
2+
# Run SolarNuAna on the output of a standard reco workflow for DUNE FD-VD
3+
4+
#include "solar_ana_flash_dunevd10kt_1x8x14_3view_30deg.fcl
5+
6+
physics.analyzers.solarnuana.SaveSignalDaughters: true
Lines changed: 6 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,6 @@
1+
# Specific configuration for investigating energy depositions of signal events.
2+
# Run SolarNuAna on the output of a standard reco workflow for DUNE FD-VD
3+
4+
#include "solar_ana_flash_radiological_decay0_dunevd10kt_1x8x14_3view_30deg.fcl"
5+
6+
physics.analyzers.solarnuana.SaveSignalDaughters: true
Lines changed: 8 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,8 @@
1+
# Specific configuration for signal+bkg events with coustom flash generation
2+
# Run SolarNuAna on the output of a standard reco workflow for DUNE FD
3+
4+
#include "solar_ana_dunevd10kt_1x8x14_3view_30deg.fcl"
5+
6+
physics.analyzers.solarnuana.SaveOpFlashInfo: true
7+
physics.analyzers.solarnuana.SaveSignalOpHits: true
8+
physics.analyzers.solarnuana.GenerateAdjOpFlash: true
Lines changed: 2 additions & 50 deletions
Original file line numberDiff line numberDiff line change
@@ -1,53 +1,5 @@
11
# solar_ana_dunevd10kt_1x8x14_3view_30deg.fcl
22
#
3-
# Run SolarNuAna on the output of a standard reco workflow for DUNE FD
3+
# Run SolarNuAna on the output of a standard reco workflow for DUNE FD-VD
44

5-
#include "services_dune.fcl"
6-
#include "SolarNuAna.fcl"
7-
8-
process_name: SolarNuAna
9-
10-
services:
11-
{
12-
@table::dunefd_services
13-
TFileService: { fileName: "solar_ana_dunevd10kt_1x8x14_3view_30deg_hist.root" }
14-
TimeTracker: {}
15-
MemoryTracker: {} # default is one
16-
RandomNumberGenerator: {} # ART native random number generator
17-
FileCatalogMetadata: @local::art_file_catalog_mc
18-
@table::dunefdvd_1x8x14_3view_30deg_services
19-
# @table::dunefdvd_1x8x6_3view_simulation_services
20-
}
21-
22-
# a temporary fix until the photon back tracker gets fixed...
23-
#services.PhotonBackTrackerService: @erase
24-
25-
physics:
26-
{
27-
analyzers:
28-
{
29-
solarnuana: @local::solar_nu_ana_vd
30-
}
31-
ana: [ solarnuana ]
32-
end_paths: [ ana ]
33-
}
34-
35-
source:
36-
{
37-
module_type: RootInput
38-
maxEvents: -1 # Number of events to create
39-
}
40-
41-
outputs:
42-
{
43-
out1:
44-
{
45-
module_type: RootOutput
46-
fileName: "%ifb_solar.root"
47-
dataTier: "analized"
48-
compressionLevel: 1
49-
}
50-
}
51-
52-
services.message.destinations.LogStandardOut.threshold: "INFO"
53-
services.message.destinations.LogStandardOut.type: "cout"
5+
#include "solar_ana_dunevd10kt_1x8x14_3view_30deg.fcl"

0 commit comments

Comments
 (0)