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

Implementing EmPhysicsList Livermore-5D #87

Open
wants to merge 10 commits into
base: feature/compair2
Choose a base branch
from

Conversation

AdrienLaviron
Copy link

@AdrienLaviron AdrienLaviron commented Jan 24, 2024

Implemented the new physics list as

  • Livermore-5D in MCParameterFile
  • EmPhysicsList number 6 in MCPhysicsList

Uses G4BetheHeitler5DModel from Geant4 v. 10.5 (and up) as pair production model, the rest of the EM processes are from the Livermore-Pol physics list.
The triplet-to-nuclear ratio is implemented as 1/Z for all energies.

Adrien Laviron added 10 commits January 24, 2024 18:33
LivPolwBH5D is a new EM physics list (EmPhysicsList zoglauer#6 in MCPhysicsList.*) based on Livermore-Pol, with the pair production model switched to the new G4BetheHeitler5DModel from Geant4 v10.5 and superior.
It is a polarized pair production model that implements 4-momentum conservation at the vertex and a fraction 1/(Z+1) of the pair production are modelled as triplet pair production.

Files MCParameterFile.cc and MCPhysicsList.* have been modified to declare that physics list as the EM physics list zoglauer#6 and to construct that physics list upon request.

- New file G4BetheHeitler5DModel.cc comes from Geant4 v11.1, has been modified to fit in Geant4 v10. List of changes:
TheElectron -> theElectron (it changed in G4PairProductionRelModel)
fTheGamma -> theGamma (idem)
Removed all references to G4EmParameters and theManager (instance of G4EmParameters). It means that we cannot shunt the physics (impose raw conversion instead of the default atom conversion, impose either nuclear or triplet pair conversion instead of the default "choose randomly between" behavior) from commands anymore.
SelectTargetAtom -> SelectRandomAtom (simple name change, also removes the need for GetLogKineticEnergy)
G4Element::GetZasInt -> lrint(G4Element::GetZ)

- New file G4BetheHeitler5DModel.hh : No changes from Geant4 v11.1

- New files LivermorePolarizedBH5DPhysics.* : Same as the Geant4 v10.02 except it uses G4BetheHeitler5DModel for pair production. No changes to the .hh.
…nd file names to Livermore-5D / Livermore5D*
FIX: *.sim not appearing in revan
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