the 35ton data David Adams BNL June 8 2016 DUNE 35t simreco Updated 08jun2016 1130 EDT Introduction G oals I would like to study detector response in the 35ton TPC We see clear tracks in 35ton data ID: 798630
Download The PPT/PDF document "Measuring track response in" is the property of its rightful owner. Permission is granted to download and print the materials on this web site for personal, non-commercial use only, and to display it on your personal computer provided you do not modify the materials and that you retain all copyright notices contained in the materials. By downloading content from our website, you accept the terms of this agreement.
Slide1
Measuring track response inthe 35-ton data
David AdamsBNLJune 8, 2016
DUNE 35t sim/reco
Updated 08jun2016 11:30 EDT
Slide2IntroductionGoalsI would like to study detector response in the 35-ton TPC
We see clear tracks in 35-ton dataClearly evident in collection viewOften (but not always) there are corresponding tracks in induction planesSee my talk at the the May 11 35-ton sim/reco meetingCan we compare signal shapes with MC?Bad channels, stuck bits and high noise make this challenging
D. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 20162
Slide3IntroductionProcedureMostly working in my own analysis module:
DXDisplay Repository: https://github.com/dladams/dune_extensionsDevelop and ask for tools that can be used in this and other environmentsTSI tools in dunetpc
Using sliced samples from recent reprocessing/pnfs/lbne/scratch/users/php13tkw/v05_10_00/slice/run*/*/lbne*_split.rootThanks to Karl for processing theseOnly looking at run 13843 right nowD. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 20163
Slide4Bad channelsI tried using the bad channel list in dunetpcI.e. the service with channels listed in FCL file
I am working with regularly-updated head of dunetpcIt removes some bad channels but not all and it removes channels that appear to be goodSee following slidesShould I make my own bad channel list?
D. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 20164
Slide5TPC 1 spectrum keeping all channels
D. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 2016
5
Slide6TPC 1 spectrum excluding “bad” channels
D. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 2016
6
Slide7Sticky bitsSticky bits are commonAssume sticking any time the low six bits are all on or all off
I.e. ADC count is 63, 64, 127, 128, 191, 192, …, 4031, 40320 is underflow4095 is overflowFollowing pages show this is pretty commonI have modified my code so bins with sticky bits can be ignoredADC spectra have holes but look much nicerNew service StandardRawDigitExtractService
extracts data from larsoft RawDigit,subtracts pedestal andsets flag indicating underflow, overflow or sticky bitUsed here to identify sticky bitsD. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 20167
Slide8TPC 1 sticky bits
D. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 2016
8
Slide9Channel zoom of preceding page
D. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 2016
9
Slide10NoiseWell known that much of noise is coherentI.e. strong correlation with nearby channels
But correlation is complex and strongest for groups of 64 or 32With many exceptions within these groupsFollowing page show overlay of a signal in a channel with the spectra from many other channels in the same group of 32Bins with sticky bits are removedSignal from track is evident in the spectraNoise shows strong albeit not perfect correlation
Assuming this spectrum is typical, it should be possible in some cases to use neighbors to estimate the background under the signalExcluding nearest neighbors polluted by signalSelect tracks where neighbors have consistent spectraD. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 201610
Slide11ADC spectra from run 13843 event 1
D. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 2016
11Pedestal-subtractedADC countTickChannel 938Neighboring channelsin same group of 32
Slide12Plan of actionIdentify bad channels by eyeLook at 2D ADC spectra and sticky-bit plots
Ignore bins with apparent sticky bitsDevelop procedure for coherent BG estimationFind sets of channels with coherent BG signalsStart with groups of 32Can we include channels with different wire orientations?Must exclude channels with tracksSerious issue because our goal is to study the response to a track
I.e. neighbors are guaranteed to be polluted by track signalWe must avoid the nearest neighborsSurvey dataFind track peaks with BG that can be estimatedAccumulate BG-subtracted track spectraLike to have track angles—use reco (3D) tracks?D. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 201612
Slide13SW organizationDeveloping tools for “data preparation”Pedestal subtraction
Identification of bad channelsIdentification of bins with underflow, overflow or stuck bitsCoherent noise subtractionPackaging these tools as TSI services in dunetpcEasy to mix and match different service implementationsProvide high-level steeringHigh-level service that starts from larsoft RawDigit
And art module wrapper for this serviceBoth will work for all Dune detectors (35t, FD, protoDUNE)Use existing code where availableRepackage as servicesMore at next FD sim/reco meetingD. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 201613
Slide14Summary/conclusionsI continue to look at 35t dataBut not a full-time activity
Goal to study the response to charged tracksCompare with and eventually tune MCThere are challenges to using the 35t data in this mannerMany bin with stuck bits—ignore theseHigh noise levelChannel-to-channel coherence provides an avenue to remove much of this noiseBut coherence is complicated and extends to a limited # channels
And there will likely be signal pollution in the channels used to estimate the noisePlan outlined for extracting response spectraDo we want to carry this out?Generally used code being packaged as dunetpc servicesD. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 201614
Slide15ExtrasOnline-ordered channel-tick spectra w/o and w/ bad channelsExample track
spectrumCorrelation plotsD. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 2016
15
Slide16Online-ordered spectrum for all channelsD. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 2016
16
Slide17Online-ordered spectrum w/o “bad” channels”D. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 2016
17
Slide18Event 1, track 1
D. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 201618
Slide19Track in collection planeD. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 2016
19
Slide20Track in induction planeD. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 2016
20
Slide21Correlation plots
D. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 201621
Slide22APA 3 (1/2)
D. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 201622
Slide23APA 3 (2/2)
D. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 201623
Slide24APA 2 (1/2)
D. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 201624
Slide25APA 2 (2/2)
D. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 201625
Slide26APA 0 (1/2)
D. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 201626
Slide27APA 0 (2/2)
D. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 201627
Slide28APA 1 (1/2)
D. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 201628
Slide29APA 1 (2/2)
D. Adams, BNL DUNE 35t sim/reco Track response in 35-ton June 8, 201629