flamedisx 0.1.1

Fast likelihood analysis in many dimensions for xenon TPCs

Flamedisx

Fast likelihood analysis in many dimensions for xenon TPCs.

By Jelle Aalbers, Bart Pelssers, and Christian Antiochi

Description

Flamedisx aims to increase the practical number of dimensions (e.g. s1, s2, x, y, z and time) and parameters (g1, g2, recombination model coefficients, electron lifetime, ...) in LXe TPC likelihoods.

Traditionally, we evaluate (the probability density functions used in) our likelihoods using histograms created from high-statistics MC simulations. We precompute these histograms for several parameter combinations, then interpolate between them during inference ("verical template morphing" in collider physics jargon). The precomputation time is exponential in the number of likelihood/histogram dimensions and the number of parameters used.

Flamedisx instead computes the probability density directly at each observed event, without using MC integration (or approximating the model). The commonly used LXe emission model is simple enough that the integral equivalent to an MC simulation can be computed with a few matrix multiplications, at a speed of a few ms -- instead of a high-statistics MC simulation that takes O(minute) or more.

This has several advantages:

• Each event has its "private" detector model computation at the observed (x, y, z, time), so making the likelihood time- and position dependent incurs no additional computational burden.
• The likelihood for a dataset takes O(seconds) to compute, so we can do this at each of optimizer's proposed points during inference. We thus remove the precomputaion step exponential in the number of parameters -- and can thus fit a great deal more parameters.
• Since the likelihood consists of deterministic matrix multiplications, it can be implemented in tensorflow/pytorch. This enables automatic differentiation, which unlocks the gradient during minimizing, drastically reducing the number of needed interations for a fit or profile likelihood.

Note this is under construction, so it probably has numerous bugs and little or no documentation.

0.1.1 / 2019-07-21

• 5x speedup for Hessian (#9)
• Fix pip install

0.1.0 / 2019-07-16

• Batching (#7)
• Inference (#6)
• Ported to tensorflow / GPU support (#1, #2, #3, #5)

0.0.1 / 2019-03-17

• Initial numpy-based version