Download Cross-gray

Reconstructing Neutrino Interactions in
Liquid Argon TPCs
Ben Newell
Steve Dennis
1
Outline

LAr-TPCs

Automation desirable

Algorithmic recognition
2
Cellular Automata

Conway's 'Game of Life'

Local rules

Cell states update simultaneously
3
CATS – The Cellular Automaton at HERA-B

HERA-B experiment uses eight 'superlayers'

Create 'track segments' between layers

Cellular automaton on track segments
4
Cellular Automata for Track Reconstruction

Cells have an index – initially 1

Local – only neighbours


Common endpoint

'Breaking angle'
Principal Direction
5
The CA Algorithm – Forward Pass

For each cell:



Look for leftward neighbours
Check if any have same index
Mark index to update

Update all indices

Repeat
6
The CA Algorithm – Reverse Pass

Start at highest index cell

Run to cell of index 1 using steps of 1

Mark cells used

Repeat with unused cells

Build all possible paths

Result: List of track candidates
7
CARLA – A Cellular Automaton for Reconstruction in
Liquid Argon

Implemented in Python

Extra steps required to suit our needs
8
Clustering the Data

LAr-TPC resolution ~mm

Thousands of voxels in principal direction

Performance problems

Clustering



Voxel size
Clustering orthogonal to principal direction
Reject 'lone' points
9
Post-Processing

Problems:



Breaking
Kinks
Clones

Filtering by shared points

Track cleaning


Breaker
Merger
10
Generalisation to 3D

Simple to work in higher dimensions

Directionality

May miss tracks

Solution: permute the axes and run on each

Recombine the results
11
CARLA in 3D
12
Parameters for reconstruction

Voxel size

Clustering radius

Cell tolerance

Filtering tolerance

Breaking Angle

Merger


Direction Tolerance

Distance Tolerance
Breaker

Correlation Tolerance

Segment length
13
Early results
14
15
16
17
18
19
20
21
22
Efficiency of CARLA
23
2D: Efficiency of reconstructing correct 2 tracks
24
2D: Variation of efficiency with breaking angle
25
2D: Variation of efficiency with voxel size
26
3D: Opening Angle Variance
27
3D: Opening Angle Variance
28
CARLA in 3D
29
Future Developments

Improvements to filtering

Documentation

User Interface
30