refactor: Stable RNG for smeared digitization in Examples

[andiwand]

Our smeared digitization will create very different random numbers if the simulations changes. This can be unexpected for pile-up studies where you would think your measurements from the hard scatter stay the same and you just get more pile-up hits.

I am trying to improve the situation by creating stable random number streams for each particle on a sensitive. To avoid correlations for multiple hits of a single particle on a single surface (should only happen for loopers) I introduced an index which is also part of the seed for the rng.

[coderabbitai]

Important

Review skipped

Draft detected.

Please check the settings in the CodeRabbit UI or the .coderabbit.yaml file in this repository. To trigger a single review, invoke the @coderabbitai review command.

You can disable this status message by setting the reviews.review_status to false in the CodeRabbit configuration file.

---

Thank you for using CodeRabbit. We offer it for free to the OSS community and would appreciate your support in helping us grow. If you find it useful, would you consider giving us a shout-out on your favorite social media?

❤️ Share

• X
• Mastodon
• Reddit
• LinkedIn

🪧 Tips

Chat

There are 3 ways to chat with CodeRabbit:

• Review comments: Directly reply to a review comment made by CodeRabbit. Example:
  • I pushed a fix in commit <commit_id>, please review it.
  • Generate unit testing code for this file.
  • Open a follow-up GitHub issue for this discussion.
• Files and specific lines of code (under the "Files changed" tab): Tag @coderabbitai in a new review comment at the desired location with your query. Examples:
  • @coderabbitai generate unit testing code for this file.
  • @coderabbitai modularize this function.
• PR comments: Tag @coderabbitai in a new PR comment to ask questions about the PR branch. For the best results, please provide a very specific query, as very limited context is provided in this mode. Examples:
  • @coderabbitai gather interesting stats about this repository and render them as a table. Additionally, render a pie chart showing the language distribution in the codebase.
  • @coderabbitai read src/utils.ts and generate unit testing code.
  • @coderabbitai read the files in the src/scheduler package and generate a class diagram using mermaid and a README in the markdown format.
  • @coderabbitai help me debug CodeRabbit configuration file.

Note: Be mindful of the bot's finite context window. It's strongly recommended to break down tasks such as reading entire modules into smaller chunks. For a focused discussion, use review comments to chat about specific files and their changes, instead of using the PR comments.

CodeRabbit Commands (Invoked using PR comments)

• @coderabbitai pause to pause the reviews on a PR.
• @coderabbitai resume to resume the paused reviews.
• @coderabbitai review to trigger an incremental review. This is useful when automatic reviews are disabled for the repository.
• @coderabbitai full review to do a full review from scratch and review all the files again.
• @coderabbitai summary to regenerate the summary of the PR.
• @coderabbitai resolve resolve all the CodeRabbit review comments.
• @coderabbitai configuration to show the current CodeRabbit configuration for the repository.
• @coderabbitai help to get help.

Other keywords and placeholders

• Add @coderabbitai ignore anywhere in the PR description to prevent this PR from being reviewed.
• Add @coderabbitai summary to generate the high-level summary at a specific location in the PR description.
• Add @coderabbitai anywhere in the PR title to generate the title automatically.

CodeRabbit Configuration File (.coderabbit.yaml)

• You can programmatically configure CodeRabbit by adding a .coderabbit.yaml file to the root of your repository.
• Please see the configuration documentation for more information.
• If your editor has YAML language server enabled, you can add the path at the top of this file to enable auto-completion and validation: # yaml-language-server: $schema=https://coderabbit.ai/integrations/schema.v2.json

Documentation and Community

• Visit our Documentation for detailed information on how to use CodeRabbit.
• Join our Discord Community to get help, request features, and share feedback.
• Follow us on X/Twitter for updates and announcements.

[github-actions]

📊: Physics performance monitoring for a3036b6

Full contents

physmon summary

• ✅ Particles fatras
• ✅ Particles geant4
• ✅ Particles ttbar
• ✅ Vertices ttbar
• 🔴 Truth tracking (KF)
• 🔴 Truth tracking (GSF)
• 🔴 Truth tracking (GX2F)
• 🔴 Truth tracking (KF refit)
• 🔴 Truth tracking (GSF refit)
• 🔴 CKF finding performance | trackfinding | single muon | truth smeared seeding
• 🔴 CKF fitting performance | trackfinding | single muon | truth smeared seeding
• 🔴 CKF track summary | trackfinding | single muon | truth smeared seeding
• ✅ Seeding trackfinding | single muon | truth estimated seeding
• 🔴 CKF finding performance | trackfinding | single muon | truth estimated seeding
• 🔴 CKF fitting performance | trackfinding | single muon | truth estimated seeding
• 🔴 CKF track summary | trackfinding | single muon | truth estimated seeding
• ✅ Seeding trackfinding | single muon | default seeding
• 🔴 CKF finding performance | trackfinding | single muon | default seeding
• 🔴 CKF fitting performance | trackfinding | single muon | default seeding
• 🔴 CKF track summary | trackfinding | single muon | default seeding
• ✅ Seeding trackfinding | single muon | orthogonal seeding
• 🔴 CKF finding performance | trackfinding | single muon | orthogonal seeding
• 🔴 CKF fitting performance | trackfinding | single muon | orthogonal seeding
• 🔴 CKF track summary | trackfinding | single muon | orthogonal seeding
• ✅ Seeding trackfinding | 4 muon x 50 vertices | default seeding
• 🔴 CKF finding performance | trackfinding | 4 muon x 50 vertices | default seeding
• 🔴 CKF fitting performance | trackfinding | 4 muon x 50 vertices | default seeding
• 🔴 CKF track summary | trackfinding | 4 muon x 50 vertices | default seeding
• 🔴 Ambisolver finding performance | trackfinding | 4 muon x 50 vertices | default seeding
• 🔴 IVF notime | trackfinding | 4 muon x 50 vertices | default seeding
• 🔴 AMVF gauss notime | trackfinding | 4 muon x 50 vertices | default seeding
• 🔴 AMVF grid time | trackfinding | 4 muon x 50 vertices | default seeding
• 🔴 Seeding trackfinding | ttbar with 200 pileup | default seeding
• 🔴 CKF finding performance | trackfinding | ttbar with 200 pileup | default seeding
• 🔴 CKF fitting performance | trackfinding | ttbar with 200 pileup | default seeding
• 🔴 CKF track summary | trackfinding | ttbar with 200 pileup | default seeding
• 🔴 Ambisolver finding performance | trackfinding | ttbar with 200 pileup | default seeding
• 🔴 AMVF gauss notime | trackfinding | ttbar with 200 pileup | default seeding
• 🔴 AMVF grid time | trackfinding | ttbar with 200 pileup | default seeding

[andiwand]

Tested this locally and works fine

[andiwand]

cc @benjaminhuth @paulgessinger @asalzburger

I don't understand how this changes physmon so much? Did we have some weird correlations before? Or do we now?

[benjaminhuth]

One thing I'm not sure is if the the properties of any / this particular RNG are similar in situations like

rng = rng_from_seed(seed)
seq1 = [ rng() for  i in range(N) ]
# vs
seq2 = [ rng_from_seed(i)() for i in range(N) ]

Because effectively we have a measurement unique RNG, that is only called once or twice, right?

Another thing one could try to improve the quality of the seed is maybe using std::seed_seq, instead of just summing up the ids.

[andiwand]

looks like std::seed_seq only does 32 bits

[sonarqubecloud]

Quality Gate passed

Issues
1 New issue
0 Accepted issues

Measures
0 Security Hotspots
0.0% Coverage on New Code
0.0% Duplication on New Code

See analysis details on SonarQube Cloud