Created attachment 614 [details] out_In - I tryed geant4.10.5, 4.10.6 and 4.10.3 and I used QGSP_BIC* pl = new QGSP_BIC(); pl->RegisterPhysics(new G4RadioactiveDecayPhysics()); So I made very simple geometry: initial beam and In target, if I use neutrons (for example 1 eV energy) there are excited states of In for example: In116[289.660] In116[127.267] In114[501.948] In114[190.268] But in case of using gamma with energy 14 MeV or 8 MeV or 20 MeV (doesn’t matter which energy), excited states are not appeared at all. But there are should be these states (channels are on the attached image "list NR" in archive NR.rar). I used more than 100 000000 gammas. Also I installed indium plate and used neutron and gamma beam with energies as on the image. All reactions with neutrons work fine (marked yellow color) but with gamma reactions doesn’t work at all (marked red color). (10 000 000 events each run). I also did simulation in Hadr03 example with using QGSP_BIC physic list and G4RadioactiveDecayPhysics. I simulated 100 000 000 gammas with 14 MeV energy into In plate 0.25 mm. And I didn’t see any exited states of In (out_In in archive NR.rar). And I also tryed it with other material (about 10 different materials) but the result is the same. Steps to Reproduce: Just run Hadr03 example using In.mac in archive NR.rar and paste: QGSP_BIC* pl = new QGSP_BIC; pl->RegisterPhysics(new G4RadioactiveDecayPhysics()); runManager->SetUserInitialization(pl); instead of Hadr03 PhysicsList.
Created attachment 615 [details] NR
Hi Aleksei, The gamma-nuclear models in QGSP_BIC and related physics lists do not have this capability. The HP models do not handle incident gammas either. For this, you must use a physics list with the LEND models. These models use the GND database which can be obtained at ftp://gdo-nuclear.ucllnl.org/ Note that these models are still under development and may have some problems.
(In reply to dennis.herbert.wright from comment #2) > Hi Aleksei, > > The gamma-nuclear models in QGSP_BIC and related physics lists do not have > this capability. The HP models do not handle incident gammas either. > For this, you must use a physics list with the LEND models. These models > use the GND database which can be obtained at > > ftp://gdo-nuclear.ucllnl.org/ > > Note that these models are still under development and may have some > problems. Hello! Thank you very much for your answer. Actually I have already tried to use physics list with LEND model (Shielding LEND), but still saw nothing. Also I tryed to use LEND model by means: G4ProcessManager* pManager = G4Gamma::Gamma()->GetProcessManager(); G4PhotoNuclearProcess* process = new G4PhotoNuclearProcess(); G4LENDorBERTModel* lend = new G4LENDorBERTModel(G4Gamma::Gamma()); lend->SetMaxEnergy(80*MeV); process->RegisterMe(lend); G4LENDCombinedCrossSection* lendXS = new G4LENDCombinedCrossSection(G4Gamma::Gamma()); process->AddDataSet(lendXS); pManager->AddDiscreteProcess(process); but I didn't see any metastable states from photons. Also I tryed new model (patch02, geant4.10.06): G4ProcessManager* pManager = G4Gamma::Gamma()->GetProcessManager(); G4PhotoNuclearProcess* process = new G4PhotoNuclearProcess(); G4LENDGammaModel* lend = new G4LENDGammaModel(G4Gamma::Gamma()); lend->SetMaxEnergy(80*MeV); process->RegisterMe(lend); G4LENDGammaCrossSection* lendXS = new G4LENDGammaCrossSection(G4Gamma::Gamma()); process->AddDataSet(lendXS); pManager->AddDiscreteProcess(process); and the situation is the same. What did I miss? Should I use a specific physics list? Could you give an advice about that?
