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

Thin layer optics has to be tweaked #7

Open
3 tasks
martinunland opened this issue Aug 15, 2024 · 0 comments
Open
3 tasks

Thin layer optics has to be tweaked #7

martinunland opened this issue Aug 15, 2024 · 0 comments
Assignees

Comments

@martinunland
Copy link
Collaborator

martinunland commented Aug 15, 2024

There are two issues with the newly implemented thin-layer optics physics:

  1. Unexpected Absorption/Transmission/Reflection Coefficients: The observed coefficients do not match the expected values, though the discrepancy is minor. The issue likely arises because the reflection at the glass boundary does not account for the imaginary part of the refractive index. Although Geant4 has a function (G4OpBoundaryProcess::GetReflectivity) that correctly calculates the Fresnel coefficients, it seems this function is not being invoked in the current setup.

  2. Transmission Angle Calculation: The transmission angle of photons is being calculated using Geant4's default Snell's law, which does not include the imaginary part of the refractive index. The important code can be found starting at line 1777 of OMSimOpBoundaryProcess.cc. The difference introduced by this omission is expected to be very small.

Action Items:

  • Identify which property needs to be defined to ensure G4OpBoundaryProcess::GetReflectivity is called at the glass boundary.
  • Simulate the four-layer system to verify if this resolves the issue.
  • Adjust the costheta variable used to calculate the new photon momentum and confirm that the simulation behaves as expected.
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