LARSOFT REFACTORISATION AND FUTURE FAR DETECTOR PRODUCTION PLANS - DOM BRAILSFORD HIGH ENERGY WORKING GROUP 15TH OCTOBER 2021
←
→
Page content transcription
If your browser does not render page correctly, please read the page content below
LArSoft refactorisation
and future far detector
production plans
Dom Brailsford
High energy working group
15th October 2021
MCC11 and the current state of
FD simulation/reconstruction
• The last production that involved simulating LBNF neutrino interactions in the DUNE FD was
MCC11, circa 2018
• MCC11 has formed the benchmark that all ongoing FD sim/reco developments compare to
• Some of those benchmarks are no longer relevant
• MCC11 les do not run out of the box with the latest dunetpc
• Conversely, some areas of the current DUNE FD work ow still rely on ancient services/
tunings/con gurations from the MCC11 era
• Many developments featured in the ProtoDUNE work ow have not yet been ported to
DUNE FD
• We are now updating the dunetpc sim/reco work ow and will then run a new FD production
• For this production spin, we’re only expecting to simulate LBNF interactions
2
fi
fi
fl
fl
fl
Updates to the sim/reco
that are HE-related
• A non-exhaustive list of improvements that might be relevant to the low energy
group are:
• Refactored LArG4 (H. Wenzel)
• 2D signal simulation/processing via wirecell (BNL team)
• Wire length-dependent noise as used in ProtoDUNE and MicroBooNE (various)
• ProtoDUNE-style hit disambiguation using spsolve (C. Backhouse and R. Sulej)
• A reorganisation/restructure of the ART fcl services and fcl les for the dune far
detector
• Updated lifetime and di usion constants based on ProtoDUNE studies and
measurements (various)
3
ff
fi
Refactored LArG4
• The ‘legacy’ larsim has been used extensively in DUNE FD, but has a few limitations
• An ‘all in one’ simulation
• A physics list that is hard to change/extend
• Di cult to retrieve the energy deposits in the active detector volume
• Many others
• The refactored larg4 xes all of those problems.
• For anyone who is active in simulation, this is a great addition to the work ow
• If you’re removed from the simulation, you shouldn’t notice a di erence in your analysis
• One key bene t is the particle tracking and electron drift simulation are completely separate
entities -> you are able to use completely di erent drift simulation models without needing
to edit larsim
4
ffi
fi
fi
ff
ff
fl
Wirecell for DUNE FD provided by H. Yu et. al. at BNL Plots by W. Gu, A. Paudel
Signal simulation and Reconstructed cosmic dQ/dx MPVs in protodune
processing using wirecell
• TDR-era drift simulation simulates the e ect an electron has on a wire
as it marches towards it
• In truth, the electron induces a signal on neighbouring wires too
• The wirecell package from BNL provides
• a 2D drift simulation, encapsulating this e ect
• A 2D wire signal deconvolution, disentangling this correlated e ect
• Extensively tested in ProtoDUNE and MicroBooNE
5
ff
ff
ff
Wire-length dependent
noise simulation
• TDR-era simulation uses an
exponential function to model
frequency
H. Yu
• Noise studied extensively in
• MicroBooNE
• ProtoDUNE
• ICEBERG
• New noise model introduces a
frequency spectrum that is wire-
length dependent (RHS plot)
6
Using 35t-era
C. Backhouse and R. Sulej disambiguation
Hit disambiguation
• Hit disambiguation for DUNE FD static since the TDR
• The in-use algorithm was developed for the 35t
and matches hit triplets using timing information
• Top image show pandora event display of di-muon
event in FD using the 35t-era disambiguation
ProtoDUNE-era
• Ghost tracks near the anode gap in the induction disambiguation
views caused by disambiguation
• Bottom display is the same event using the
ProtoDUNE-era disambiguation
• Protodune disambiguation is Space Point Solver
(C. Backhouse) + SP->Hit disambiguation (R. Sulej)
• Ghost tracks near the anode plane are gone
7Services/fcl reorganisation
• A dry subject but unfortunately one you do
need to hear about
• Dunetpc fcl services restructured to make
the refactored larg4 the standard work ow
• Legacy service blocks have been created
to make maintain the older g4 work ow
• Service fcls have been generally
reorganised and grouped by detector
• The ‘standard’ job fcls run the refactored
larg4, wirecell and other improvements
listed in this talk
• ‘Legacy’ fcl les are provided to run the
older work ow
8
fl
fi
fl
flProduction status
• All core features are in place and we’re entering the validation stage using
LBNF simulations
• All changes lives on feature/dbrailsf_refactor
• Once validations are nished we will start a new production to produce new
LBNF neutrino events
• The service restructure and legacy fcls are already on dunetpc develop with
some caveats
• The refactor -> standard and standard -> legacy restructure is in place
but some additional modi cations completely undo those changes
• The ‘legacy’ job fcls are already provided but the ‘standard’ still run the
legacy work ow
9
fl
fi
fiHow will the changes affect
the HE group
• The standard software is going to change so your HE sim/
reco work ows will pick up these changes
• We are already providing a 1D signal simulation/signal
processing work ow using the refactored LArG4
• This will include the protodune noise model (J. Shen
and M. Wang)
• If you want to completely reject the changes, legacy fcl
les are also provided to run the software as it is now
10
fi
fl
flFuture high energy
reconstruction in pandora
• A few analyses are now starting to make use of the Pandora
pattern recognition
• A longstanding feature in pandora is that the vertex
reconstruction assumes the topology is accelerator driven
• Reconstructed clusters each provide a reconstructed
vertex candidate
• An MVA is used to select the ‘best’ vertex candidate
• Some of the discriminators in the MVA prefer the vertex
to be more upstream in Z
• Further pattern recognition algs use the vertex as an
anchor point -> the pattern recognition really relies on the Reconstructed vertex candidates
vertex being in the right place Provided by H. Lay
• There will be a more detailed talk on atmospheric
reconstruction in a few weeks
• In the interim, we suggest cheating the vertex reconstruction
11M. He
Ongoing
reconstruction
validation
• Initial atmospheric neutrino analysis using pandora
• Top RHS plot reconstructed and true vertical direction
components of the leading lepton
• Some events are aimed the right way
• Some are aimed exactly the wrong way
• Bottom plot shows separation in X position between the true
and reconstructed vertex
• Bi-modal structure expected to be due to a non-zero T0
• Pandora’s neutrino reconstruction treats the entire TPC
as one global volume to cluster across
• An uncorrected T0 is problematic when clustering
across the central CPA
• We have very recently requested a production respin of the
atmospheric generation that xes the T0, to be used for
reconstruction studies
12
fiV. Di. Benedetto, R. Cross
Continuous Integration (CI)/
Validation system
• dunetpc makes heavy use of the CI system as part of the
LArSoft pull request system
• Developed by V. Di Benedetto
• DUNE FD more recently makes use of an extended high
statistics validation system
• Developed by Vi. Di Benedetto and lled in by R. Cross
• The validation system provides
• an easy means of quickly making high statistics plots
after making software changes
• Comparing the output to a set of reference les
• Allows you to make simple comparisons (track
completeness) OR analysis-level comparisons
(reconstructed atmospheric oscillation pro le)
• Group support would be needed to ll in the physics details
but we can help with the technical implementation
13
fi
fi
fi
fiSummary
• FD sim/reco is updating the sim/reco work ow to bring in the
latest developments from various LArTPCs
• We are ultimately looking to upgrade the MCC11 LBNF
simulation that people still use for comparisons
• The changes will propagate to the high energy working group
but it will still be possible to run the old work ows
• We will be back to give a more detailed talk on atmospheric
reconstruction in a few weeks
• Please consider using the validation system
14
fl
flYou can also read
