Skip to content

Compute ceneter time of particle transit in laser beam (tao) using max likelihood#82

Merged
jtgasparik merged 9 commits intoARM-DOE:mainfrom
jtgasparik:tao_max_likelihood
Apr 6, 2026
Merged

Compute ceneter time of particle transit in laser beam (tao) using max likelihood#82
jtgasparik merged 9 commits intoARM-DOE:mainfrom
jtgasparik:tao_max_likelihood

Conversation

@jtgasparik
Copy link
Copy Markdown
Collaborator

Estimate tau_hat using the Moteki & Kondo (2008) grid-search MLE.

This implements the normalized-derivative framework from the paper:
  - S'(t)/S(t) = I'(t)/I(t) + C_sca'(t)/C_sca(t)  [Eq. (4)]
  - For the Gaussian beam, I'(t)/I(t) = -(t - tau)/sigma^2  [Eq. (5)]
  - Likelihood of a p-point subset is multivariate normal [Eq. (A.9)]
  - Mean vector uses Eq. (A.4)
  - Covariance/variance use Eqs. (A.10a), (A.10b)
  - Signal-noise model uses Eq. (A.6)
  - Error propagation to the derivative uses Eq. (A.7)
  - tau is found by grid search over a 1D parameter space [Appendix A.5]

Computed values are tested for non-evaporative particles (events 499 and 1040).

@jtgasparik jtgasparik requested a review from rcjackson March 20, 2026 21:02


@dataclass(frozen=True)
class MLEConfig:
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I would add a docstring to the MLEConfig class so that users can know what these settings are.

@jtgasparik jtgasparik requested a review from rcjackson April 6, 2026 14:03
@jtgasparik jtgasparik merged commit 45ba380 into ARM-DOE:main Apr 6, 2026
20 checks passed
@jtgasparik jtgasparik deleted the tao_max_likelihood branch April 6, 2026 14:10
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants