First Result of WireCell Signal Processing in Proto

  • Slides: 14
Download presentation
First Result of Wire-Cell Signal Processing in Proto. DUNE Wenqiang Gu on behalf of

First Result of Wire-Cell Signal Processing in Proto. DUNE Wenqiang Gu on behalf of the Wire-Cell team BNL Proto. DUNE Sim/Reco Meeting, 11/28/2018

Outline • Signal Processing in Wire-Cell toolkit • 2 D deconvolution • Region of

Outline • Signal Processing in Wire-Cell toolkit • 2 D deconvolution • Region of interest (ROI) Ionization Electron Signal Processing in Single Phase LAr. TPCs I. JINST 13 P 07006 (2018) • Software integration in LAr. Soft • Performance of signal processing • Full TPC simulation sample • Data sample 2

Signal processing (SP): deconvolution & filter • Fourier transform Deconvolution + Filter Inverse Fourier

Signal processing (SP): deconvolution & filter • Fourier transform Deconvolution + Filter Inverse Fourier transform Liquid Argon TPC Signal Formation, Signal Processing and Hit Reconstruction Bruce Baller, JINST 12 (2017) no. 07, P 07010 3

Long-range induction 2 D deconvolution • However, the induction from neighboring ionization electrons has

Long-range induction 2 D deconvolution • However, the induction from neighboring ionization electrons has to been considered The inversion of matrix R can again be done with deconvolution through 2 -D FFT 2 D: both time and wires dimensions 4

Just 2 D deconvolution will not be enough ROI + Adaptive Baseline • The

Just 2 D deconvolution will not be enough ROI + Adaptive Baseline • The bi-polar nature of induction signal amplifies low-frequency noise during deconvolution • Improved through region of interest (ROI) and the adaptive baseline technique Only for illustration, not a proto. DUNE version Given N time bins with 2 MHz digitization frequency, • The highest freq is 1 MHz • The lowest freq (above 0) is 2/N MHz e. g. , 200 bins 10 k. Hz • Obviously not sensitive to noise < 2/N MHz • Adaptive baseline linear baseline correction instead of flat baseline correction 5

Software integration in LAr. Soft framwork Raw Data Raw Decoder ADC mitigation (module: Data.

Software integration in LAr. Soft framwork Raw Data Raw Decoder ADC mitigation (module: Data. Prep) • sticky code • FEMB 302 • undershoot raw: : Raw. Digit / recob: : Wire Reco Chain … Downstream analysis recob: : Wire (2 D decon. ) Larwirecell • consumes raw: : Raw. Digit, or recob: : Wire Noise Filter / ADC mitigation (in development) Wire. Cell Toolkit • • • coherent noise ADC nonlinearity etc. Signal Processing • 2 D deconvolution • ROI Imaging TPC drift simulation … 6

Software integration in LAr. Soft (cont’) • Wire-Cell Toolkit • Repository https: //github. com/Wire.

Software integration in LAr. Soft (cont’) • Wire-Cell Toolkit • Repository https: //github. com/Wire. Cell • Document https: //wirecell. github. io/ • larwirecell (https: //cdcvs. fnal. gov/redmine/projects/larwirecell) --- usage example $ lar -n 1 -c Run. Raw. Decoder. fcl np 04_raw_run 005141_0017_dl 1. root $ lar -n 1 -c nfsp. fcl np 04_raw_run 005141_0017_dl 1_decode. root $ lar -n 1 -c wcls-nf-sp. fcl np 04_raw_run 005141_0017_dl 1_decode_reco. root # get output. root $ lar -n 1 -c eventdump. fcl output. root Upstream noise filtered raw waveforms from Data. Prep 7

SP performance test in a full TPC simulation APA#3 A MIP (~5000 e/mm) track

SP performance test in a full TPC simulation APA#3 A MIP (~5000 e/mm) track from bottom to top across the TPC APA#2 APA#6 APA#4 Full TPC includes: • Ionized electron absorption, diffusion, fluctuation • Field response, electronics response, etc. • Noise Clear tracks from SP Consistent with the channel map 8

SP Performance in proto. DUNE beam data 1 D deconvolution Run 5141, Event 23865

SP Performance in proto. DUNE beam data 1 D deconvolution Run 5141, Event 23865 Threshold: 5 From the offline reco chain (proto. DUNE_reco_data. fcl) 2 D deconvolution* *: There is still room for improving the software filter and some thresholds, etc. **: Noise filtering has not been applied here for both 1 D & 2 D. 9

Detailed example 1: U plane After Noise Filtering 1 -D Deconvolution • Re-normalize 1

Detailed example 1: U plane After Noise Filtering 1 -D Deconvolution • Re-normalize 1 D & 2 D to the same scale • No significant negative component after 2 D deconvolution • Long tracks (in time) are more visible in the 2 D deconvolution Wire-Cell 2 -D Deconvolution Ch 545 10

Example 2: V plane After Noise Filtering 1 -D Deconvolution 2 -D Deconvolution 11

Example 2: V plane After Noise Filtering 1 -D Deconvolution 2 -D Deconvolution 11

Example 3: W plane After Noise Filtering 1 -D Deconvolution 2 -D Deconvolution •

Example 3: W plane After Noise Filtering 1 -D Deconvolution 2 -D Deconvolution • 1 D & 2 D deconvolution are consistent in collection plane 12

Other efforts from Wire-Cell team TPC signal/noise simulation Or Data Existing Noise filtering Module

Other efforts from Wire-Cell team TPC signal/noise simulation Or Data Existing Noise filtering Module ADC sticky code mitigation ADC nonlinearity correction Coherent noise removal Ledge and dead channel identification More work needed Already have good progresses in • Ledge identification by Zeyuan Yu • Noisy channel by Carlos Sarasty Electronics response calibration as part of ADC nonlinearity calibration Signal Processing Module FFT to Frequency domain: i) Misconfigured channels ii) Timing issue for FEMB 302 iii) Baseline undershoot 2 D deconvolution + ROI for general signal processing High-level reconstruction modules 13

Summary • With 2 D deconvolution + ROI, Wire-Cell toolkit has successfully achieved the

Summary • With 2 D deconvolution + ROI, Wire-Cell toolkit has successfully achieved the signal processing in proto. DUNE • Still have some room for improving software filters, thresholds, etc. • Wire-Cell toolkit has been integrated in the LAr. Soft via an interface module larwirecell • Consumes the existing ADC mitigation in the reco chain for the December production • More efforts will be made to improve the noise filtering and ADC problems in proto. DUNE 14