Core-Removal¶
Original Author: Doug Rutledge
This module splits a given pulse-series (InputRecoPulseSeries
) into
two parts: The one close around a given vertex (VertexName
) (it is
called the core pulses and stored in CorePulsesName
) and the part
far out (corona pulses stored in OutputRecoPulseSeries
).
The pulses of an OM are part of the core pulses if the OM distance
to the vertex position is smaller than a radius \(R=k R_{SPE}\),
where \(k\) is a scaling factor given by
SPEFraction
and the SPE radius \(R_{SPE}\) (stored as
SPERadiusName
) is calculated with the formula
where \(\lambda_{Att}\) (LambdaAttn
), \(C_\lambda\),
\(E_c\) (CriticalEnergy
) and \(R_{min}\)
(MinimumSPERadius
) are given as parameters. The spline parameter
is calculated as \(c_{spline}=\frac{\lambda_{Att}}{2\ln E_c}\)
and constant is calculated as
\(c_\lambda=R_{min} - \frac{\lambda_{Att}}{2}\ln E_c\).
The energy \(E\) is taken from the vertex particle or if
NChCalib
is set estimated from the numbers of channels with
different algorithms depending on DetectorConf
(0: Amanda, 1:
IceCube, 2: Combined).