Dear all,
I implemented two cases, one where I tried adding QMD to a modular physics list and another adding BIC using the following lines respectively:
physicslist->RegisterPhysics(new G4IonQMDPhysics())
and then in another simulation
physicslist->RegisterPhysics(new G4IonBinaryCascadePhysics())
When I simulate a 64Ni beam at 250Mev/u incident on a 1975 mg/cm2 aluminium target in geant4:
The fraction of 64Ni that survive the target is 0.901,
64Co are produced at a rate of 2e-5.,
62Ni are produced at a rate of 1.06e-4,
When there numbers are compared to the rates I get in LISE++ (which I’m using as a control to validate the geant4 simulation) http://lise.nscl.msu.edu/lise.html:
0.89 64Ni should survive the target,
64Co should be produced at a rate of 6.87e-4
62Ni should be produced at a rate of 2.35e-3
With the exception of the 64Ni, it seems the reaction rates are off by an order of magnitude.
I notice that both QMD and BIC produce reaction products at identical rates, is there something wrong with the way I’ve added them to the physics list? Otherwise I can only conclude the models currently available aren’t accurate enough for the nuclei I’m trying to study.
Best,
Michael