Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

fix: FPE in FATRAS hit direction #2189

Merged
merged 1 commit into from
Jun 9, 2023

Conversation

paulgessinger
Copy link
Member

@paulgessinger paulgessinger commented Jun 9, 2023

I've encountered FPEs in the calculation of the average direction of FATRAS hits.

Currently we do this:

  • Divide four-momenta before and after divide by 2x the norm of the three-momentum vector
  • Add both four-vectors, then take only the three-momentum out, then normalize thos

This fails if the norm of either of the three momenta is invalid or zero. I think Eigen has checks against this but we do the division manually.

I think this is equivalent to:

  • Getting the normalized three-momenta
  • Adding them and dividing by 2, then normalize again

Thoughts?

@paulgessinger paulgessinger added this to the next milestone Jun 9, 2023
@andiwand andiwand added Component - Fatras Affects the Fatras module automerge labels Jun 9, 2023
@codecov
Copy link

codecov bot commented Jun 9, 2023

Codecov Report

Merging #2189 (a92bb9d) into main (078a01e) will not change coverage.
The diff coverage is n/a.

@@           Coverage Diff           @@
##             main    #2189   +/-   ##
=======================================
  Coverage   49.41%   49.41%           
=======================================
  Files         441      441           
  Lines       25176    25176           
  Branches    11617    11617           
=======================================
  Hits        12441    12441           
  Misses       4481     4481           
  Partials     8254     8254           

📣 We’re building smart automated test selection to slash your CI/CD build times. Learn more

@kodiakhq kodiakhq bot merged commit 4083f82 into acts-project:main Jun 9, 2023
53 checks passed
@github-actions github-actions bot removed the automerge label Jun 9, 2023
@github-actions
Copy link

github-actions bot commented Jun 9, 2023

📊 Physics performance monitoring for a92bb9d

Summary
Full report
Seeding: seeded, truth estimated, orthogonal
CKF: seeded, truth smeared, truth estimated, orthogonal
IVF: seeded, truth smeared, truth estimated, orthogonal
AMVF: seeded, truth smeared, truth estimated, orthogonal
Ambiguity resolution: seeded, orthogonal
Truth tracking
Truth tracking (GSF)

Vertexing

Vertexing vs. mu
IVF seeded

IVF truth_smeared

IVF truth_estimated

IVF orthogonal

AMVF seeded

AMVF truth_smeared

AMVF truth_estimated

AMVF orthogonal

Seeding

Seeding seeded

Seeding truth_estimated

Seeding orthogonal

CKF

CKF seeded

CKF truth_smeared

CKF truth_estimated

CKF orthogonal

Ambiguity resolution

seeded

Truth tracking (Kalman Filter)

Truth tracking

Truth tracking (GSF)

Truth tracking

@@ -88,9 +88,9 @@ class Hit {
}
/// Average normalized particle direction vector through the surface.
Vector3 unitDirection() const {
Copy link
Contributor

Choose a reason for hiding this comment

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

just for reference - it might be better to not average the direction but stick to either before or after. I think for digitization we should use before?

@paulgessinger paulgessinger deleted the fix/fatras-hit-dir branch June 9, 2023 13:44
@paulgessinger paulgessinger modified the milestones: next, v27.0.0 Jun 20, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Component - Fatras Affects the Fatras module
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

2 participants