Skip to content

Conversation

@Aske-Rosted
Copy link
Collaborator

This PR introduces two new features to the I3HighestEparticleExtractor.

  • primary_type - Records the PID of the primary of the particle that produces the final particle used as base for the Highest Energy Particle
  • primary_is_nu - Checks whether the primary of the Highest Energy Particle is a neutrino.

primary_is_nu can be used to check whether the particle found was produced by a neutrino or is from a coincident muon in the case of the nugen files.

The changes introduced by #822 included the following function call.

frame[self.mctree].is_in_subtree(
                        p.id, track.GetI3Particle().id
                    )

Which is very computationally expensive especially for high energy CORSIKA files. Furthermore the check is irrelevant for CORSIKA where everything is considered to be in the sub-tree of the atmospheric shower.

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.

1 participant