search

celeritas-project / celeritas

Celeritas is a new Monte Carlo transport code designed to accelerate scientific discovery in high energy physics by improving detector simulation throughput and energy efficiency using GPUs.
https://celeritas-project.github.io/celeritas/user/index.html
Other
63 stars 34 forks source link

LEGEND integration #1347

Open sethrj opened 3 months ago

sethrj commented 3 months ago

This is the primary tracking issue for integrating Celeritas with the simulation workflow for the LEGEND experiment. The primary Celeritas contact points are the "assignees" to the right.

Core capabilities

Notes

drbenmorgan commented 3 months ago

I'm happy to be added as a preliminary contact point if you want here as I'm officially in LEGEND.

sethrj commented 2 months ago

From Matthew Green:

To start, we simulate scintillation photons produced in liquid argon, in which we approximate the blue curve in the figure as a gaussian distribution centered at 128nm with sigma = 2.929nm, emitted isotropically. scintillation-probability-density

We voxelize the detector in 5mmx5mmx5mm voxels, generate scintillation photons, track. In addition to the scintillation, we also include wavelength shifting in TPB and wavelength-shifting optical fibers and scintillation in PEN plastic.

Everthing inside the double-walled steel cyrostat is argon, but it is partitioned by that vertical copper can into 2 regions. In the outer region, where things are essentially axially symmetric, we use the axial symmetry to effectively generate a 2D probability map in R, Z. In the inner region, where the axial symmetry is broken, we generate the map fully in 3D. In both cases we generate random points within the voxel boundary, sample whether that point is in LAr, and if so, throw a photon. If not, resample. legend-flask

Maps look something like this: optical-map Inner region (left-most) is a 2d slice from a 3d map, middle/outer regions are the 2D map in R,z.