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

Question/Issue with Status ID #817

Open
Meatbally2 opened this issue Dec 11, 2024 · 1 comment
Open

Question/Issue with Status ID #817

Meatbally2 opened this issue Dec 11, 2024 · 1 comment

Comments

@Meatbally2
Copy link
Contributor

Hi,

I generated some eHe3 scattering events with BeAGLE, and used eicsmear and hepmc3ascii2root to convert the output to hepmc3 and hepmc3.tree.root formats. The status IDs from "particle.status" in the file looked like what I would expect for the hepmc3 status code for my sample. However, when I used the files (I tried both) as an input to run epic simulation, the IDs from "MCParticle.generatorStatus" in the simulation output were mostly 0. Below are the screenshots from the hepmc3.tree.root and the simulation output file:

particle.status from hepmc3.tree.root
Screen Shot 2024-12-10 at 19 40 43

MCParticle.generatorStatus from edm4hep.root (In this figure, I only analysed 10 events from the above file)
Screen Shot 2024-12-10 at 19 37 21

I am not sure if this is a bug for running simulation on BeAGLE events, or if I am doing anything wrong. Any suggestions will be appreciated! Thank you so much!

Environment: (where does this bug occur, have you tried other environments)

  • Which branch (main for latest released): main (from the most up to date eic-shell)
  • Which revision (HEAD for the most recent):
  • Any specific OS or system where the issue occurs? I ran the analysis on BNL SDCC
  • Any special versions of ROOT or Geant4?

Steps to reproduce: (give a step by step account of how to trigger the bug)

Expected Result: (what do you expect when you execute the steps above)

I expect the MCParticle.generatorStatus to have similar result as particle.status in the hepmc3.tree.root.

Actual Result: (what do you get when you execute the steps above)

@veprbl
Copy link
Member

veprbl commented Dec 11, 2024

You need to filter out some entries by MCParticles.simulatorStatus's bit 30
https://github.com/key4hep/EDM4hep/blob/main/edm4hep.yaml#L303 to remove Geant tracks

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

No branches or pull requests

2 participants