Efficiency Correction of Particle-Averaged Quantities

Basic Information

• Paper ID: 2510.13838
• Title: Efficiency correction of particle-averaged quantities
• Authors: Masakiyo Kitazawa, ShinIchi Esumi, Takafumi Niida, Toshihiro Nonaka
• Classification: physics.data-an hep-ex hep-ph nucl-ex nucl-th
• Publication Date: October 11, 2025 (arXiv preprint)
• Paper Link: https://arxiv.org/abs/2510.13838v1

Abstract

This paper derives analytical formulas for reconstructing particle-averaged quantities from experimental results affected by particle measurement efficiency losses. These formulas are based on the assumption that the probabilities of observing individual particles are mutually independent. The study finds that these formulas are inconsistent with the intuitive formulas traditionally employed.

Research Background and Motivation

Core Problem

In relativistic heavy-ion collision experiments, researchers frequently need to measure the average values of physical quantities for specific particle species in each collision event, such as average transverse momentum p_T and flow anisotropy parameters (directed flow and elliptic flow). These quantities take the form:

$\frac{1}{N}\sum_{i=1}^{N}\xi_i$

where N is the number of particles in the event and ξ_i is the physical quantity carried by the i-th particle.

Significance

1. Experimental Universality: Such measurements are very common in relativistic heavy-ion collisions, including measurements of flow anisotropy parameters and average transverse momentum
2. Physical Significance: Event-by-event fluctuation studies are important for exploring physical phenomena such as the shape of colliding nuclei
3. Detector Effects: In actual experiments, detector performance is imperfect, and the probability (efficiency) of observing each particle is less than 1

Limitations of Existing Methods

1. Lack of Theoretical Foundation: The traditional efficiency correction formulas used in literature lack rigorous theoretical derivation
2. Problems with Intuitive Formulas: Traditional formulas take the form $\langle\langle\frac{\sum_i \xi_i/r_i}{\sum_i 1/r_i}\rangle\rangle$, but their correctness has not been verified
3. Limitations of Numerical Methods: While expansion methods are feasible, they typically require substantial numerical resources and lack transparency in the reconstruction process

Core Contributions

1. Derivation of Rigorous Analytical Formulas: Under the assumption of mutually independent particle observation probabilities, provides rigorous mathematical derivation for efficiency correction of particle-averaged quantities
2. Extension to Higher-Order Correlations: Extends the formulas to efficiency correction of second-order, third-order, and higher-order correlation functions
3. Treatment of Multiple Efficiency Scenarios: Extends from uniform efficiency to the general case of different efficiencies among particles
4. Verification of Formula Correctness: Validates the correctness of new formulas through simple models and demonstrates deviations in traditional formulas

Detailed Methodology

Task Definition

Input: Experimental observational data affected by efficiency losses (observed particle number n, physical quantities q) Output: True particle-averaged quantity $\langle Q/N \rangle$Constraints:

• Particle observation probabilities are mutually independent
• Efficiency values are known and determined for all observed particles
• Other detector effects (such as measurement errors, misidentification) are not considered

Core Methodological Framework

1. Uniform Efficiency Case

For the case where all particles have the same efficiency r, the derivation proceeds as follows:

Probability Distribution Relationship: $\tilde{P}(n,q) = \sum\int P \sum_{\{b_i\}} (1-r)^{N-n}r^n \delta_{n,\sum_i b_i} \delta(q-\sum_i b_i\xi_i)$

where b_i is a binary variable indicating whether the i-th particle is observed.

Generating Function Method: Introducing the generating function $\tilde{G}(s,\theta) = \sum\int \tilde{P} s^n e^{q\theta}$, by taking derivatives with respect to θ and integrating:

$\langle\frac{Q}{N}\rangle = \langle\langle\frac{q}{n}(1-\alpha^n)\rangle\rangle_{n\neq 0}$

where $\alpha = (r-1)/r$.

2. Multiple Efficiency Case

For the general case where particles have different efficiencies r_i, the final formula is:

$\langle\frac{Q}{N}\rangle = \langle\langle\sum_{i=1}^{n} \xi_i k_{1;i}\rangle\rangle_{n\neq 0}$

where the weight factor is: $k_{1;i} = \left(\prod_i r_i\right)^{-1} \int_0^1 d\sigma \prod_{j\neq i}(\sigma + r_j - 1)$

Higher-Order Correlation Functions

For the second-order correlation function $\frac{\sum_{i\neq j}\xi_i^{(1)}\xi_j^{(2)}}{N(N-1)}$, the correction formula is:

$\langle\frac{\{Q_1Q_2\}}{N(N-1)}\rangle = \langle\langle\sum_{i\neq j} q_{1,i}q_{2,j}k_{2;i,j}\rangle\rangle_{n\neq 0,1}$

where: $k_{2;i,j} = \left(\prod_i r_i\right)^{-1} \int_0^1 d\sigma' \int_0^{\sigma'} d\sigma \prod_{l\neq i,l\neq j}(\sigma + r_l - 1)$

Technical Innovations

1. Generating Function Technique: Cleverly employs the generating function method to handle complex probability distribution transformations
2. Integral Transformation: Connects true and observed distributions through appropriate integral transformations
3. Weight Factor Design: Derives precise weight factors that account for the mutual influence of all particle efficiencies
4. Systematic Extension: Provides systematic treatment from first-order to higher-order correlation functions

Experimental Setup

Verification Model

A simple binary model is employed for verification:

• Particle physical quantity $\xi_i = \pm 1$
• Fixed particle number $N_+ = N_- = N/2$ (producing equal numbers of positive and negative particles)
• Different efficiencies $r_+, r_-$

Evaluation Methods

1. Analytical Calculation: Direct computation of theoretical expectation values
2. Numerical Verification: Calculation of reconstructed values using derived formulas
3. Comparative Analysis: Comparison with results from traditional formulas

Experimental Results

Main Results

1. Verification of New Formulas

In the simple model, the new formula perfectly reconstructs the true values:

• Theoretical value: $\langle Q/N \rangle = 0$ (symmetric case)
• Reconstructed value: Exact zero value obtained using formula (41)

2. Deviations in Traditional Formulas

The traditional formula $\langle\langle\frac{n_+/r_+ - n_-/r_-}{n_+/r_+ + n_-/r_-}\rangle\rangle$ exhibits systematic deviations:

• Non-zero results for all efficiency values except $r_- = 1$
• Deviation increases with decreasing efficiency
• Approaches correct values only in the limits $N \to \infty$ or $r \to 1$

3. Convergence Analysis

For cases with $r < 1/2$, since $|\alpha| \geq 1$, the term $\alpha^n$ may cause numerical instability, requiring special attention.

Important Findings

1. Physical Meaning of Correction Terms: The additional term $-\langle\langle(q/n)\alpha^n\rangle\rangle$ represents the contribution from n=0 events
2. Efficiency Threshold Effect: When $r \leq 1/2$, reconstruction may become difficult
3. Applicability Range of Traditional Formulas: Traditional formulas are valid only under specific limiting conditions

Related Work

Efficiency Correction Research

1. Particle Number Correction: Existing literature has established efficiency correction methods for particle number cumulants
2. Higher-Order Correlations: Analytical methods exist for handling efficiency correction of particle number higher-order correlations
3. Expansion Methods: Serve as numerical alternatives but with high computational cost and low transparency

Innovations in This Work

Compared to existing work, this paper provides for the first time a rigorous analytical treatment for particle-averaged quantities, filling a theoretical gap in this field.

Conclusions and Discussion

Main Conclusions

1. Rigorous Formulas: Derives rigorous efficiency correction formulas for particle-averaged quantities and their higher-order correlations
2. Defects in Traditional Formulas: Demonstrates the incorrectness of traditional intuitive formulas
3. Practical Value: Provides theoretical foundation and practical tools for experimental analysis

Limitations

1. Independence Assumption: Requires mutually independent particle observation probabilities, which may not be fully satisfied in actual experiments
2. Numerical Complexity: Formulas involve multiple integrals, which may be computationally complex
3. Efficiency Range: For very low efficiency cases ($r < 1/2$), reconstruction may be unstable

Future Directions

1. Numerical Optimization: Development of efficient numerical computation methods, such as utilizing relationships of elementary symmetric polynomials
2. Self-Correlation Extension: Treatment of higher-order correlation functions containing self-correlation terms
3. Different Acceptance Regions: Extension to correlation functions between completely different acceptance regions

In-Depth Evaluation

Strengths

1. Theoretical Rigor: Provides complete mathematical derivation, filling theoretical gaps
2. Methodological Systematicity: Systematic treatment from simple to complex, from first-order to higher-order
3. Practical Value: Directly applicable to heavy-ion collision experiment analysis
4. Sufficient Verification: Formula correctness verified through specific models

Weaknesses

1. Computational Complexity: Numerical computation in practical applications may be relatively complex
2. Assumption Limitations: Independence assumption may not be fully satisfied in actual experiments
3. Lack of Intuitiveness: Complex formula form lacks intuitive physical interpretation

Impact

1. Theoretical Contribution: Provides important theoretical tools for particle physics experimental data analysis
2. Experimental Application: Directly applicable to data analysis in RHIC, LHC, and other large-scale experiments
3. Methodological Value: The generating function method may inspire solutions to other similar problems

Applicable Scenarios

1. Heavy-Ion Collision Experiments: Flow anisotropy analysis in RHIC, LHC, and other experiments
2. Particle Physics Experiments: Any experiments involving particle-averaged quantity measurements
3. Statistical Physics: Statistical analysis of particle systems with efficiency losses

References

The paper cites 39 related references, covering multiple aspects including heavy-ion collision physics, efficiency correction methods, and flow anisotropy measurements, demonstrating the comprehensiveness and depth of the research. Important references include pioneering work on efficiency correction methods 19-28 and key results from heavy-ion collision experiments 1-6.

Overall Assessment: This is a high-quality theoretical physics paper that addresses an important problem in experimental data analysis. Although the mathematical derivation is complex, it provides important theoretical tools for experimental physicists and possesses significant academic value and practical significance.