ACEM: Approximate Cherenkov Emission Model

acem is a Python package designed for generation of Cherenkov light yield profiles from particle showers. One of the main goals is to more accurately include shower-to-shower fluctuations while preserving model simplicity.

Installation

You can install acem directly from PyPI.

pip install acem

Alternatively, clone this repository and from the directory run pip install . or to work with scripts in misc/ use pip install .[misc]. The latter will allow you to run checks including plots like this

cd misc/scripts/
./check.py 3 42 (energy at 1TeV, initial seed of 42)

It's also possible to install directly from github with

pip install git+https://github.com/pathtorepo/acem

Usage

Here is an example of how to import the models and use the package:

from acem import model, media

par = model.Parametrization1D(media.ICE)

# Sample 100 showers profiles initiated by a 1 TeV pi+
shos = par.sample(211, 1e3, 100)

shos will be a list of Shower1D objects, which is a simple container that consists of the amplitude and 1D shape of the Cherenkov-weighted track length (units are in cm). Its .dldx(xs) method can be used to evaluate the shower profile at some given distance(s) xs from the start position.

Example

The 1D models are based on FLUKA simulations. As an example, when used in combination with a hadronization library such as PYTHIA8, it can be used to generate neutrino DIS shower profiles.