zoglauer / megalib

MEGAlib - the Medium-Energy Gamma-ray Astronomy library
http://megalibtoolkit.com
Other
41 stars 34 forks source link

COSIMA Activation - Skips Events and Freezes #98

Open joshhooker opened 6 days ago

joshhooker commented 6 days ago

On the 11.2 branch, the simulation will run until the event number is skipped and the events stop updating. All of the files do not continue to be updated.

Running COSIMA with 1 thread and the following options:

Physics parameters:

PhysicsListHD qgsp-bic-hp PhysicsListEM LivermorePol DecayMode ActivationBuildup

Storage parameters

StoreCalibrated true StoreSimulationInfo all DiscretizeHits true

Particle range cuts:

DefaultRangeCut 0.005

Simulate the primary protons and all secondaries:

Run Step1 Step1.FileName output.Phase1.BuildUp Step1.Events 200000000 Step1.IsotopeProductionFile output.Phase1.BuildUp.dat

Screenshot 2024-09-23 at 1 06 04 PM

Any ideas on what would cause this?

zoglauer commented 6 days ago

Thanks for reporting.

Which branch are you using? All branches using Geant4 11.2 are early development branches, which just have been tested with COSI. And I know there are still a few bugs left to fix.

Do you want to switch to the main branch which uses Geant4 10.2?

joshhooker commented 6 days ago

Thank you for the fast reply, Andreas!

I am currently on the feature/geant_11.2 branch.

I am unable to use MEGAlib with Geant4 10.2 because I am using materials with a Z > 92 and Geant4 10.2 appears to have issues with these heavy materials even with the AllowForHeavyElements flag turned on.

zoglauer commented 6 days ago

Is there a chance you could either email me or upload a geometry file which reproduces this issue?