Applying Quenching Corrections in CalDigi
- For each MCIntegratingHit including particles with Z > 1, store:
- Flag indicating Z > 1 present
- Ein for each particle with Z > 1
- Eout for each particle with Z > 1
- Z for each particle with Z > 1
- If, in a single MCIntegratingHit (i.e. for one "cell" in a single xtal), there are N particles with Z > 1:
E effective = [E IHit – (? i (E in i – E out i ))] +
? i [(E in i – E out i ) × Q((E in i + E out i )/2 , Z)] [1]- Q(E,Z) is the quenching correction as a function of E and Z, assumed here = 1 for Z = 1
- This formulation assumes that all energy from ionization is deposited locally i.e. that the contribution from longer range delta electrons is negligible
- This assumption has to be tested!
- Note that violation of this assumption causes errors in quenching correction, not errors in deposited energy directly
The first term in [1] is the (approximate) contribution of Z = 1 particles
The second term in [1] is the the sum of contributions from all Z > 1 particles contributing to this integrating hit (i.e. in this cell), corrected for quenching
- Note that, due to delta electrons