SBNSoftware · gputnam · Jul 9, 2024 · Jul 9, 2024 · Jul 9, 2024 · Jul 9, 2024
diff --git a/fcl/detsim/detsim_2d_icarus_SConly.fcl b/fcl/detsim/detsim_2d_icarus_SConly.fcl
@@ -0,0 +1,9 @@
+#include "detsim_2d_icarus.fcl"
+
+physics.producers.tpcsim: @local::icarus_simwire_wirecell_SConly
+
+physics.simulate: ["rns", "tpcsim"]
+physics.stream: []
+outputs: {}
+
+process_name: TPCSim
diff --git a/fcl/detsim/detsim_2d_icarus_fitFR.fcl b/fcl/detsim/detsim_2d_icarus_fitFR.fcl
@@ -1,4 +1,4 @@
 #include "detsim_2d_icarus.fcl"
 
-physics.producers.daq: @local::icarus_simwire_wirecell_fitSR
 physics.producers.daq.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-multitpc-sim-drift-simchannel.jsonnet"]
+physics.producers.daq: @local::icarus_simwire_wirecell_fitSR_P0nom
diff --git a/fcl/reco/Definitions/sp_filter_parameters.fcl b/fcl/reco/Definitions/sp_filter_parameters.fcl
@@ -0,0 +1,15 @@
+BEGIN_PROLOG
+Wiener_tight_U_sigma: 0.065
+Wiener_tight_U_power: 4.4
+Wiener_tight_V_sigma: 0.065
+Wiener_tight_V_power: 2.6
+Wiener_tight_W_sigma: 0.07
+Wiener_tight_W_power: 3.4
+
+Wire_ind_sigma: 0.4
+Wire_col_sigma: 2.2
+
+Gaus_wide_sigma: 0.06
+
+gain_ADC_per_e: 0.01214 # ADC / e-. For frame_scale
+END_PROLOG
diff --git a/fcl/reco/Definitions/stage0_icarus_defs.fcl b/fcl/reco/Definitions/stage0_icarus_defs.fcl
diff --git a/fcl/reco/Stage0/Run2/stage0E_run2_wc_icarus.fcl b/fcl/reco/Stage0/Run2/stage0E_run2_wc_icarus.fcl
@@ -0,0 +1,10 @@
+#include "stage0_run2_wc_icarus.fcl"
+
+physics.path: [
+    filterdataintegrity,
+    @sequence::icarus_stage0_PMT,
+    @sequence::icarus_stage0_CRT,
+    daqTPCROI,
+    @sequence::icarus_stage0_multiTPC_2d_TPC_E,
+    @sequence::icarus_stage0_EastHits2d_TPC
+]
diff --git a/fcl/reco/Stage0/Run2/stage0E_run2_wc_icarus_mc.fcl b/fcl/reco/Stage0/Run2/stage0E_run2_wc_icarus_mc.fcl
@@ -0,0 +1,10 @@
+#include "stage0_run2_wc_icarus_mc.fcl"
+
+# Run the PMT+CRT stage0, initial TPC and East 2D deconv
+physics.path: [   @sequence::icarus_stage0_mc_PMT,
+                  MCDecodeTPCROI,
+                  @sequence::icarus_stage0_2d_multiTPC_E,
+                  @sequence::icarus_stage0_mc_crt
+              ]
+
+process_name: MCstage0E
diff --git a/fcl/reco/Stage0/Run2/stage0E_run2_wcdnn_icarus.fcl b/fcl/reco/Stage0/Run2/stage0E_run2_wcdnn_icarus.fcl
@@ -0,0 +1,14 @@
+#include "stage0E_run2_wc_icarus.fcl"
+#
+# Turn on PCA noise suppression
+# physics.producers.MCDecodeTPCROI.DecoderTool.DenoiserType: "pca"
+
+physics.producers.decon2droiEE.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+physics.producers.decon2droiEW.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+physics.producers.decon2droiWE.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+physics.producers.decon2droiWW.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+
+physics.producers.gaushit2dTPCEE.CalDataModuleLabel: "decon2droiEE:dnnsp"
+physics.producers.gaushit2dTPCEW.CalDataModuleLabel: "decon2droiEW:dnnsp"
+physics.producers.gaushit2dTPCWE.CalDataModuleLabel: "decon2droiWE:dnnsp"
+physics.producers.gaushit2dTPCWW.CalDataModuleLabel: "decon2droiWW:dnnsp"
diff --git a/fcl/reco/Stage0/Run2/stage0E_run2_wcdnn_icarus_mc.fcl b/fcl/reco/Stage0/Run2/stage0E_run2_wcdnn_icarus_mc.fcl
@@ -0,0 +1,14 @@
+#include "stage0E_run2_wc_icarus_mc.fcl"
+#
+# Turn on PCA noise suppression
+# physics.producers.MCDecodeTPCROI.DecoderTool.DenoiserType: "pca"
+
+physics.producers.decon2droiEE.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+physics.producers.decon2droiEW.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+physics.producers.decon2droiWE.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+physics.producers.decon2droiWW.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+
+physics.producers.gaushit2dTPCEE.CalDataModuleLabel: "decon2droiEE:dnnsp"
+physics.producers.gaushit2dTPCEW.CalDataModuleLabel: "decon2droiEW:dnnsp"
+physics.producers.gaushit2dTPCWE.CalDataModuleLabel: "decon2droiWE:dnnsp"
+physics.producers.gaushit2dTPCWW.CalDataModuleLabel: "decon2droiWW:dnnsp"
diff --git a/fcl/reco/Stage0/Run2/stage0W_run2_wc_icarus.fcl b/fcl/reco/Stage0/Run2/stage0W_run2_wc_icarus.fcl
@@ -0,0 +1,8 @@
+#include "stage0_run2_wc_icarus.fcl"
+
+physics.path: [
+    @sequence::icarus_stage0_multiTPC_2d_TPC_W,
+    @sequence::icarus_stage0_WestHits2d_TPC
+]
+
+process_name: stage0W
diff --git a/fcl/reco/Stage0/Run2/stage0W_run2_wc_icarus_mc.fcl b/fcl/reco/Stage0/Run2/stage0W_run2_wc_icarus_mc.fcl
@@ -0,0 +1,6 @@
+#include "stage0_run2_wc_icarus_mc.fcl"
+
+# Run the West 2D deconv and finish TPC stage0
+physics.path: [@sequence::icarus_stage0_2d_multiTPC_W]
+
+process_name: MCstage0W
diff --git a/fcl/reco/Stage0/Run2/stage0W_run2_wcdnn_icarus.fcl b/fcl/reco/Stage0/Run2/stage0W_run2_wcdnn_icarus.fcl
@@ -0,0 +1,14 @@
+#include "stage0W_run2_wc_icarus.fcl"
+
+# Turn on PCA noise suppression
+# physics.producers.MCDecodeTPCROI.DecoderTool.DenoiserType: "pca"
+
+physics.producers.decon2droiEE.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+physics.producers.decon2droiEW.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+physics.producers.decon2droiWE.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+physics.producers.decon2droiWW.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+
+physics.producers.gaushit2dTPCEE.CalDataModuleLabel: "decon2droiEE:dnnsp"
+physics.producers.gaushit2dTPCEW.CalDataModuleLabel: "decon2droiEW:dnnsp"
+physics.producers.gaushit2dTPCWE.CalDataModuleLabel: "decon2droiWE:dnnsp"
+physics.producers.gaushit2dTPCWW.CalDataModuleLabel: "decon2droiWW:dnnsp"
diff --git a/fcl/reco/Stage0/Run2/stage0W_run2_wcdnn_icarus_mc.fcl b/fcl/reco/Stage0/Run2/stage0W_run2_wcdnn_icarus_mc.fcl
@@ -0,0 +1,14 @@
+#include "stage0W_run2_wc_icarus_mc.fcl"
+
+# Turn on PCA noise suppression
+# physics.producers.MCDecodeTPCROI.DecoderTool.DenoiserType: "pca"
+
+physics.producers.decon2droiEE.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+physics.producers.decon2droiEW.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+physics.producers.decon2droiWE.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+physics.producers.decon2droiWW.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+
+physics.producers.gaushit2dTPCEE.CalDataModuleLabel: "decon2droiEE:dnnsp"
+physics.producers.gaushit2dTPCEW.CalDataModuleLabel: "decon2droiEW:dnnsp"
+physics.producers.gaushit2dTPCWE.CalDataModuleLabel: "decon2droiWE:dnnsp"
+physics.producers.gaushit2dTPCWW.CalDataModuleLabel: "decon2droiWW:dnnsp"
diff --git a/fcl/reco/Stage0/Run2/stage0_run2_wc_icarus.fcl b/fcl/reco/Stage0/Run2/stage0_run2_wc_icarus.fcl
@@ -32,39 +32,6 @@ outputs.rootOutput.outputCommands:         ["keep *_*_*_*",
                                             "keep *_daq_ICARUSTriggerV*_*"
                                            ]
 
-# Add below as per Gray Putnam (should we move to jsonnet files?)
-physics.producers.decon2droiEE.wcls_main.structs.gain0: 17.05212 
-physics.producers.decon2droiEW.wcls_main.structs.gain0: 17.05212 
-physics.producers.decon2droiWE.wcls_main.structs.gain0: 17.05212 
-physics.producers.decon2droiWW.wcls_main.structs.gain0: 17.05212 
-
-
-physics.producers.decon2droiEE.wcls_main.structs.gain1: 12.1420344 
-physics.producers.decon2droiEW.wcls_main.structs.gain1: 12.1420344 
-physics.producers.decon2droiWE.wcls_main.structs.gain1: 12.1420344 
-physics.producers.decon2droiWW.wcls_main.structs.gain1: 12.1420344 
-
-physics.producers.decon2droiEE.wcls_main.structs.gain2: 13.0261362 
-physics.producers.decon2droiEW.wcls_main.structs.gain2: 13.0261362 
-physics.producers.decon2droiWE.wcls_main.structs.gain2: 13.0261362 
-physics.producers.decon2droiWW.wcls_main.structs.gain2: 13.0261362 
-
-
-physics.producers.decon2droiEE.wcls_main.structs.shaping0: 1.3  
-physics.producers.decon2droiEW.wcls_main.structs.shaping0: 1.3  
-physics.producers.decon2droiWE.wcls_main.structs.shaping0: 1.3  
-physics.producers.decon2droiWW.wcls_main.structs.shaping0: 1.3  
-
-physics.producers.decon2droiEE.wcls_main.structs.shaping1: 1.45  
-physics.producers.decon2droiEW.wcls_main.structs.shaping1: 1.45  
-physics.producers.decon2droiWE.wcls_main.structs.shaping1: 1.45  
-physics.producers.decon2droiWW.wcls_main.structs.shaping1: 1.45  
-
-physics.producers.decon2droiEE.wcls_main.structs.shaping2: 1.3  
-physics.producers.decon2droiEW.wcls_main.structs.shaping2: 1.3  
-physics.producers.decon2droiWE.wcls_main.structs.shaping2: 1.3  
-physics.producers.decon2droiWW.wcls_main.structs.shaping2: 1.3  
-
 services.message.destinations :
 {
   STDCOUT:

diff --git a/fcl/reco/Stage0/Run2/stage0_run2_wc_icarus_mc.fcl b/fcl/reco/Stage0/Run2/stage0_run2_wc_icarus_mc.fcl
@@ -40,45 +40,10 @@ physics.producers.ophit.InputModule: "opdaq"
 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.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"
 
-# As per Gray Putname...
-physics.producers.decon2droiEE.wcls_main.structs.gain0: 17.05212 
-physics.producers.decon2droiEW.wcls_main.structs.gain0: 17.05212 
-physics.producers.decon2droiWE.wcls_main.structs.gain0: 17.05212 
-physics.producers.decon2droiWW.wcls_main.structs.gain0: 17.05212 
-
-
-physics.producers.decon2droiEE.wcls_main.structs.gain1: 12.1420344 
-physics.producers.decon2droiEW.wcls_main.structs.gain1: 12.1420344 
-physics.producers.decon2droiWE.wcls_main.structs.gain1: 12.1420344 
-physics.producers.decon2droiWW.wcls_main.structs.gain1: 12.1420344 
-
-physics.producers.decon2droiEE.wcls_main.structs.gain2: 13.0261362 
-physics.producers.decon2droiEW.wcls_main.structs.gain2: 13.0261362 
-physics.producers.decon2droiWE.wcls_main.structs.gain2: 13.0261362 
-physics.producers.decon2droiWW.wcls_main.structs.gain2: 13.0261362 
-
-
-physics.producers.decon2droiEE.wcls_main.structs.shaping0: 1.3  
-physics.producers.decon2droiEW.wcls_main.structs.shaping0: 1.3  
-physics.producers.decon2droiWE.wcls_main.structs.shaping0: 1.3  
-physics.producers.decon2droiWW.wcls_main.structs.shaping0: 1.3  
-
-physics.producers.decon2droiEE.wcls_main.structs.shaping1: 1.45  
-physics.producers.decon2droiEW.wcls_main.structs.shaping1: 1.45  
-physics.producers.decon2droiWE.wcls_main.structs.shaping1: 1.45  
-physics.producers.decon2droiWW.wcls_main.structs.shaping1: 1.45  
-
-physics.producers.decon2droiEE.wcls_main.structs.shaping2: 1.3  
-physics.producers.decon2droiEW.wcls_main.structs.shaping2: 1.3  
-physics.producers.decon2droiWE.wcls_main.structs.shaping2: 1.3  
-physics.producers.decon2droiWW.wcls_main.structs.shaping2: 1.3  
-
 # restore legacy G4 labels
 physics.producers.mcophit.SimPhotonsProducer: "largeant"
diff --git a/fcl/reco/Stage0/Run2/stage0_run2_wcdnn_icarus.fcl b/fcl/reco/Stage0/Run2/stage0_run2_wcdnn_icarus.fcl
@@ -0,0 +1,14 @@
+#include "stage0_run2_wc_icarus.fcl"
+
+# Turn on PCA noise suppression
+# physics.producers.MCDecodeTPCROI.DecoderTool.DenoiserType: "pca"
+
+physics.producers.decon2droiEE.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+physics.producers.decon2droiEW.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+physics.producers.decon2droiWE.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+physics.producers.decon2droiWW.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+
+physics.producers.gaushit2dTPCEE.CalDataModuleLabel: "decon2droiEE:dnnsp"
+physics.producers.gaushit2dTPCEW.CalDataModuleLabel: "decon2droiEW:dnnsp"
+physics.producers.gaushit2dTPCWE.CalDataModuleLabel: "decon2droiWE:dnnsp"
+physics.producers.gaushit2dTPCWW.CalDataModuleLabel: "decon2droiWW:dnnsp"
diff --git a/fcl/reco/Stage0/Run2/stage0_run2_wcdnn_icarus_mc.fcl b/fcl/reco/Stage0/Run2/stage0_run2_wcdnn_icarus_mc.fcl
@@ -0,0 +1,40 @@
+#include "stage0_run2_wc_icarus_mc.fcl"
+
+# Turn on PCA noise suppression
+# physics.producers.MCDecodeTPCROI.DecoderTool.DenoiserType: "pca"
+
+physics.producers.decon2droiEE.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+physics.producers.decon2droiEW.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+physics.producers.decon2droiWE.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+physics.producers.decon2droiWW.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-decode-to-sig-twofaced-dnnroi.jsonnet"]
+
+physics.producers.gaushit2dTPCEE.CalDataModuleLabel: "decon2droiEE:dnnsp"
+physics.producers.gaushit2dTPCEW.CalDataModuleLabel: "decon2droiEW:dnnsp"
+physics.producers.gaushit2dTPCWE.CalDataModuleLabel: "decon2droiWE:dnnsp"
+physics.producers.gaushit2dTPCWW.CalDataModuleLabel: "decon2droiWW:dnnsp"
+
+physics.producers.gaushit2dTPCEE.HitFinderToolVec.CandidateHitsPlane0.RoiThreshold:      2.5
+physics.producers.gaushit2dTPCEE.HitFinderToolVec.CandidateHitsPlane1.RoiThreshold:      2.5
+physics.producers.gaushit2dTPCEE.HitFinderToolVec.CandidateHitsPlane2.RoiThreshold:      2.5
+physics.producers.gaushit2dTPCEE.HitFilterAlg.MinPulseHeight:                            [2., 2., 2.]
+
+physics.producers.gaushit2dTPCEW.HitFinderToolVec.CandidateHitsPlane0.RoiThreshold:      2.5
+physics.producers.gaushit2dTPCEW.HitFinderToolVec.CandidateHitsPlane1.RoiThreshold:      2.5
+physics.producers.gaushit2dTPCEW.HitFinderToolVec.CandidateHitsPlane2.RoiThreshold:      2.5
+physics.producers.gaushit2dTPCEW.HitFilterAlg.MinPulseHeight:                            [2., 2., 2.]
+
+physics.producers.gaushit2dTPCWE.HitFinderToolVec.CandidateHitsPlane0.RoiThreshold:      2.5
+physics.producers.gaushit2dTPCWE.HitFinderToolVec.CandidateHitsPlane1.RoiThreshold:      2.5
+physics.producers.gaushit2dTPCWE.HitFinderToolVec.CandidateHitsPlane2.RoiThreshold:      2.5
+physics.producers.gaushit2dTPCWE.HitFilterAlg.MinPulseHeight:                            [2., 2., 2.]
+
+physics.producers.gaushit2dTPCWW.HitFinderToolVec.CandidateHitsPlane0.RoiThreshold:      2.5
+physics.producers.gaushit2dTPCWW.HitFinderToolVec.CandidateHitsPlane1.RoiThreshold:      2.5
+physics.producers.gaushit2dTPCWW.HitFinderToolVec.CandidateHitsPlane2.RoiThreshold:      2.5
+physics.producers.gaushit2dTPCWW.HitFilterAlg.MinPulseHeight:                            [2., 2., 2.]
+
+physics.analyzers.spanaE.HitProducers: ["gaushit2dTPCEE", "gaushit2dTPCEW"]
+physics.analyzers.spanaE.WireProducers: ["decon2droiEE:dnnsp", "decon2droiEW:dnnsp"]
+
+physics.analyzers.spanaW.HitProducers: ["gaushit2dTPCWE", "gaushit2dTPCWW"]
+physics.analyzers.spanaW.WireProducers: ["decon2droiWE:dnnsp", "decon2droiWW:dnnsp"]
diff --git a/icaruscode/TPC/CMakeLists.txt b/icaruscode/TPC/CMakeLists.txt
@@ -8,3 +8,4 @@ add_subdirectory(SignalProcessing)
 add_subdirectory(Simulation)
 add_subdirectory(Tracking)
 add_subdirectory(ICARUSWireCell)
+add_subdirectory(SPAna)
diff --git a/icaruscode/TPC/ICARUSWireCell/detsimmodules_wirecell_ICARUS.fcl b/icaruscode/TPC/ICARUSWireCell/detsimmodules_wirecell_ICARUS.fcl
@@ -12,9 +12,10 @@ icarus_simwire_wirecell:
     {
         tool_type: WCLS
         apps: ["Pgrapher"]
-        // logsinks: ["stdout"]
+        logsinks: ["stdout"]
         // loglevels: ["magnify:debug"]
-        plugins: ["WireCellPgraph", "WireCellGen","WireCellSio","WireCellRoot","WireCellLarsoft"]
+        loglevels: ["debug", "pgraph:info"]
+        plugins: ["WireCellPgraph", "WireCellGen","WireCellSio","WireCellRoot","WireCellLarsoft", "WireCellHio"]
         // needs to be found via your WIRECELL_PATH 
         configs: ["pgrapher/experiment/icarus/wcls-multitpc-sim-drift-simchannel-refactored.jsonnet"]
         // Contract note: these exact "type:name" must be used to identify
@@ -75,6 +76,13 @@ icarus_simwire_wirecell:
     }
 }
 
+# SimChannel only (SConly) CONFIG
+icarus_simwire_wirecell_SConly: @local::icarus_simwire_wirecell
+icarus_simwire_wirecell_SConly.wcls_main.configs: ["pgrapher/experiment/icarus/wcls-multitpc-sim-drift-simchannel-true.jsonnet"]
+icarus_simwire_wirecell_SConly.wcls_main.outputers:[
+            "wclsDepoFluxWriter:postdrift"
+]
+
 # TUNED FR CONFIG
 icarus_simwire_wirecell_fitSR: @local::icarus_simwire_wirecell
 # Add in the ER tail

diff --git a/icaruscode/TPC/ICARUSWireCell/icarus/chndb-base.jsonnet b/icaruscode/TPC/ICARUSWireCell/icarus/chndb-base.jsonnet
@@ -28,7 +28,7 @@ function(params, anode, field, n, rms_cuts=[])
 
 
     // Externally determined "bad" channels.
-    bad: [],
+    bad: [167, 1795, 1796, 1824, 1826, 1828, 1830, 1860, 1888, 1889, 1890, 1891, 1922, 2013, 2806, 3725, 3743, 4262, 5572, 8523, 8524, 8525, 8526, 8527, 8945, 9180, 9387, 9391, 9664, 10138, 12051, 14028, 14403, 15307, 15737, 16017, 16018, 16019, 16020, 16021, 16022, 16023, 16024, 16025, 16026, 16027, 16028, 16029, 16030, 16031, 16937, 17201, 19580, 23225, 23457, 23984, 23985, 25662, 27219, 27231, 27335, 27967, 28640, 28804, 28805, 28806, 28807, 28808, 28809, 28810, 30496, 30497, 30498, 30499, 30500, 30501, 30502, 30503, 30504, 30505, 30506, 30507, 30508, 30509, 30510, 30511, 30512, 30513, 30514, 30515, 30516, 30517, 30518, 30519, 30520, 30521, 30522, 30523, 30524, 30525, 30526, 30527, 30705, 31316, 31317, 31726, 34220, 38294, 41501, 41975, 42851, 43231, 43518, 44051, 44410, 49216, 49792, 49793, 49794, 49795, 49796, 49797, 49798, 49799, 49800, 49801, 49802, 49803, 49804, 49805, 49806, 49807, 49808, 49809, 49810, 49811, 49812, 49813, 49814, 49815, 49816, 49817, 49818, 49819, 49820, 49821, 49822, 49823, 50080, 50081, 50082, 50083, 50084, 50085, 50086, 50087, 50088, 50089, 50090, 50091, 50092, 50093, 50094, 50095, 50096, 50097, 50098, 50099, 50100, 50101, 50102, 50103, 50104, 50105, 50106, 50107, 50108, 50109, 50110, 50111, 50673, 51900],
 
     // Overide defaults for specific channels.  If an info is
     // mentioned for a particular channel in multiple objects in this
@@ -41,7 +41,7 @@ function(params, anode, field, n, rms_cuts=[])
       // repeat values found here in subsequent entries unless you
       // wish to change them.
       {
-        channels: util.anode_channels(n),
+        channels: util.anode_channels_twofaced(n),
         nominal_baseline: 2048.0,  // adc count
         gain_correction: 1.0,  // unitless
         response_offset: 0.0,  // ticks?