Repository navigation
Conversation
Adds two entries to AngularDistributionFormulas describing the cosmic muon flux surviving the 2450 m.w.e. rock overburden of the Canfranc Underground Laboratory. The slant rock thickness grows as h0/cos(theta), so the intensity per unit solid angle is much steeper than at sea level, I(theta) ~ cos^6(theta). Reference: Trzaska et al., Eur. Phys. J. C 79, 721 (2019). - LscMuons: sin(theta) * cos^6(theta), for the 'cosmic' spatial generator, whose projected sampling area does not depend on the direction. - LscMuonsWall: sin(theta) * cos^7(theta), for a horizontal 'wall' or 'circle' surface generator, where position and direction are sampled independently and the cos(theta) projection of the plane must be included explicitly. The sin(theta) factor is the jacobian of dOmega = sin(theta) dtheta dphi, needed because the formula is sampled with TF1::GetRandom on theta. This is the convention of SinCos2, not the one of Cos2 / Cos3. The behaviour of the existing formulas is unchanged, and a note documenting that both conventions coexist has been added to the enum. GetCosmicFluxInCountsPerCm2PerSecond() now also supports these formulas, so GetEquivalentSimulatedTime() returns the live time for this generator. Geant4 does not provide an absolute rate, so the reference measured flux through a horizontal plane is used, J = Phi_horiz * (n + 2) / (n + 1) for LscMuons and J = Phi_horiz for LscMuonsWall. Co-Authored-By: Claude Opus 5 <noreply@anthropic.com>
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
1 participant
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
Adds two entries to AngularDistributionFormulas describing the cosmic muon flux surviving the 2450 m.w.e. rock overburden of the Canfranc Underground Laboratory. The slant rock thickness grows as h0/cos(theta), so the intensity per unit solid angle is much steeper than at sea level, I(theta) ~ cos^6(theta). Reference: Trzaska et al., Eur. Phys. J. C 79, 721 (2019).
The sin(theta) factor is the jacobian of dOmega = sin(theta) dtheta dphi, needed because the formula is sampled with TF1::GetRandom on theta. This is the convention of SinCos2, not the one of Cos2 / Cos3. The behaviour of the existing formulas is unchanged, and a note documenting that both conventions coexist has been added to the enum.
GetCosmicFluxInCountsPerCm2PerSecond() now also supports these formulas, so GetEquivalentSimulatedTime() returns the live time for this generator. Geant4 does not provide an absolute rate, so the reference measured flux through a horizontal plane is used, J = Phi_horiz * (n + 2) / (n + 1) for LscMuons and J = Phi_horiz for LscMuonsWall.