Merge pull request #606 from SBNSoftware/feature/lynnt_wirecell_teardrops

WireCell NF+SP updates (teardrop mitigation)

Author: bear-is-asleep

sbndcode/WireCell/cfg/pgrapher/experiment/sbnd/chndb-base.jsonnet

@@ -20,11 +20,18 @@ function(params, anode, field, n, rms_cuts=[])
     // For MicroBooNE, channel groups is a 2D list.  Each element is
     // one group of channels which should be considered together for
     // coherent noise filtering.
-    // 5638 is the number of channels in a single APA, (1984*2 + 1670), including the 6 channel gap.
-    // 4806 is the number of channels for both induction + half of real collection wires + 6 channel gap
-    groups: [std.range(   0 + n * 5638 + g*32,    0 + n * 5638 + (g+1)*32 - 1) for g in std.range(0,150)] +
-            [std.range(4806 + n * 5638 + g*32, 4806 + n * 5638 + (g+1)*32 - 1) for g in std.range(0,26)] ,
+    // 5638 is the number of channels in a single APA, (1984*2 + 1670), 
+    // including the 6 channel gap.
+    // The induction planes have to types of ch grouping (due to FEMB), 
+    // one where the grouping is 32 channels wide and one where it is 128 channels wide.
+    // The collection planes are grouped by 64. 
 
+    groups: [std.range(   0 + n * 5638 + g*32 ,    0 + n * 5638 + (g+1)*32  - 1) for g in std.range(0,26)] + # first section of u 
+            [std.range( 832 + n * 5638 + g*128,  832 + n * 5638 + (g+1)*128 - 1) for g in std.range(0,9)] +  # second section of u
+            [std.range(1984 + n * 5638 + g*128, 1984 + n * 5638 + (g+1)*128 - 1) for g in std.range(0,9)] +  # first section of v
+            [std.range(3136 + n * 5638 + g*32 , 3136 + n * 5638 + (g+1)*32  - 1) for g in std.range(0,26)] + # second section of v 
+            [std.range(3974 + n * 5638 + g*64 , 3974 + n * 5638 + (g+1)*64  - 1) for g in std.range(0,13)] + # first half of w
+            [std.range(4806 + n * 5638 + g*64 , 4806 + n * 5638 + (g+1)*64  - 1) for g in std.range(0,13)] , # second half of w
 
     // Externally determined "bad" channels.
     //

sbndcode/WireCell/cfg/pgrapher/experiment/sbnd/chndb-perfect.jsonnet

@@ -20,9 +20,18 @@ function(params, anode, field, n, rms_cuts=[])
     // For MicroBooNE, channel groups is a 2D list.  Each element is
     // one group of channels which should be considered together for
     // coherent noise filtering.
-    groups: [std.range(   0 + n * 5638 + g*32,    0 + n * 5638 + (g+1)*32 - 1) for g in std.range(0,149)] +
-            [std.range(4806 + n * 5638 + g*32, 4806 + n * 5638 + (g+1)*32 - 1) for g in std.range(0,25)] ,
+    // 5638 is the number of channels in a single APA, (1984*2 + 1670), 
+    // including the 6 channel gap.
+    // The induction planes have to types of ch grouping (due to FEMB), 
+    // one where the grouping is 32 channels wide and one where it is 128 channels wide.
+    // The collection planes are grouped by 64. 
 
+    groups: [std.range(   0 + n * 5638 + g*32 ,    0 + n * 5638 + (g+1)*32  - 1) for g in std.range(0,26)] + # first section of u 
+            [std.range( 832 + n * 5638 + g*128,  832 + n * 5638 + (g+1)*128 - 1) for g in std.range(0,9)] +  # second section of u
+            [std.range(1984 + n * 5638 + g*128, 1984 + n * 5638 + (g+1)*128 - 1) for g in std.range(0,9)] +  # first section of v
+            [std.range(3136 + n * 5638 + g*32 , 3136 + n * 5638 + (g+1)*32  - 1) for g in std.range(0,26)] + # second section of v 
+            [std.range(3974 + n * 5638 + g*64 , 3974 + n * 5638 + (g+1)*64  - 1) for g in std.range(0,13)] + # first half of w
+            [std.range(4806 + n * 5638 + g*64 , 4806 + n * 5638 + (g+1)*64  - 1) for g in std.range(0,13)] , # second half of w
 
     // Externally determined "bad" channels.
     //

sbndcode/WireCell/cfg/pgrapher/experiment/sbnd/nf.jsonnet

@@ -61,7 +61,7 @@ function(params, anode, chndbobj, n, name='', dft=default_dft)
           wc.tn(single),
         ],
         grouped_filters: [
-          // wc.tn(grouped),
+          wc.tn(grouped),
         ],
         channel_status_filters: [
         ],

sbndcode/WireCell/cfg/pgrapher/experiment/sbnd/params.jsonnet

@@ -123,7 +123,7 @@ base {
         // The "absolute" time (ie, relative to trigger time?) that the lower edge
         // of final readout tick #0 should correspond to.  This is a
         // "fixed" notion.
-        local tick0_time = -200*wc.us,
+        local tick0_time = -205*wc.us, // updated in January 2025, front porch in data = 410 ticks
 
         // Open the ductor's gate a bit early.
         local response_time_offset = $.det.response_plane / $.lar.drift_speed,

sbndcode/WireCell/cfg/pgrapher/experiment/sbnd/sp-filters.jsonnet

@@ -75,7 +75,7 @@ local wf(name, data={}) = {
 [
   lf('ROI_tight_lf', { tau: 0.0185 * wc.megahertz }),  // 0.02 
   lf('ROI_tighter_lf', { tau: 0.145 * wc.megahertz }),  // 0.1 
-  lf('ROI_loose_lf', { tau: 0.00175 * wc.megahertz }),  // 0.0025 
+  lf('ROI_loose_lf', { tau: 0.006 * wc.megahertz }),  // updated from 0.00175 for data teardrop removal
 
   hf('Gaus_tight'),
   hf('Gaus_wide', { sigma: 0.10 * wc.megahertz }), 

sbndcode/WireCell/cfg/pgrapher/experiment/sbnd/wcls-sim-drift-depoflux-nf-sp.jsonnet

@@ -116,7 +116,7 @@ local wcls_depoflux_writer = g.pnode({
     field_response: wc.tn(tools.field),
     tick: 0.5 * wc.us,
     window_start: 0.0 * wc.ms,
-    window_duration: self.tick * 3400,
+    window_duration: self.tick * params.daq.nticks,
     nsigma: 3.0,
 
     reference_time: -1700 * wc.us,