Direct light from muon

The following example demonstrates how to calculate the number of photoelectrons for a given environment and PMT.

using LumenManufaktur

Let's create the parameters with the defaults but override the dispersion model with the one used in KM3NeT ORCA detectors in the Mediterranean:

params = LMParameters(dispersion_model=DispersionORCA)
LMParameters:
  minimum distance = 0.1 m
  module_radius = 0.25 m
  lambda min / max = 300.0 nm / 700.0 nm
  number of integrations points = 20
  degree of Legendre polynomials = 5
  dispersion model = BaileyDispersion(240.0, 1.3201, 1.4e-5, 16.2566, -4383.0, 1.1455e6)
  scattering model = LumenManufaktur.Kopelevich()
  scattering probability model = LumenManufaktur.Scatteringp00075()
  absorption model = LumenManufaktur.DefaultAbsorption()

To calculate the number of photoelectrons from direct light induced by a muon at the closest distance of R=23.5, with a standard KM3NeT PMT pointing in the direction of θ=π and ϕ=π/2, we call directlightfrommuon():

directlightfrommuon(params, LumenManufaktur.KM3NeTPMT, 23.5, π, π/2)
0.03345844066791786