0-Lifetime Physics Function
Last Update: 08-Dec-05
Implementation in Code
The 0-lifetime physics function is approximately a delta-function
in Lxy(true) or VPDL(true).
This is implemented in code as:
|
ncc(xtr,tag;δ)
= Ptag / δ
|
for 0 < xtr < δ
|
|
ncc(xtr,tag;δ)
= 0
|
for xtr >= δ
|
This is correctly normalized to Ptag if:
∑i=bins
ncc(xi,tag;δ)
Δxbin
with
Problems with this Implementation
-
When looping over xtr (to fill histo's or
calculate pdf normalizations, for example) have to be sure that
x < δ is a point in the loop.
-
Normalization integrals have to use
Δxbin = δ
which generally means a lot of steps since δ should be small
to avoid the problem below.
-
δ needs to be small or the 0-lifetime PDF will get too much
weight (will be non-zero for too many events).
Solution to Problem 2
The issue is evaluating the sum giving the Phys-Func * Eff
normalization:
1/A =
∑i Δx
ncc(xi)
ε(xi)
This can simply be set to
Then, as long as histo's are filled with an x-loop using
Δx = δ,
the full PDF will be correctly normalized:
PDFcc =
H(y)
∑ ΔK H(K)
ncc(xtr)
ε(xtr) A
where y = Ptm for Lxy fits or
σ(vpdl) for VPDL fits.