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:

with

Problems with this Implementation

  1. 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.
  2. Normalization integrals have to use Δxbin = δ which generally means a lot of steps since δ should be small to avoid the problem below.
  3. δ 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:

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:

where y = Ptm for Lxy fits or σ(vpdl) for VPDL fits.