| Message: Model for reaction Li7(d,p)Li8 | Not Logged In (login) |
|
Hello,
I am new to geant4 and try to simulate a Li7 beam with a kinetic energy of about 25 MeV which should interact with a deuterium gas target. Thus Li8 should be produced by the reaction Li7(d,p)Li8. The detector design und electromagnetic interaction (I studied the energy deposition in the target first) seems to work fine. But now I am trying to implement a hadronic inelastic model. Concerning the following page it seems that the only model I can use for so low energies is the Wilson abrasion model. http://geant4.cern.ch/support/proc_mod_catalog/models/hadronic/G4WilsonAbrasionModel.html I have added this model to my physics list and simulated 100,000,000 events to see if there is any Li8 produced. But it is not the case! Possibly I am using the wrong cross sections. Here a snippet of my physics list: ################################################################################ G4HadronInelasticProcess* IonInelasticProc = new G4HadronInelasticProcess("IonInelastic", G4GenericIon::GenericIon()); // applying a model // Wilson Abrasion Model G4WilsonAbrasionModel* wilsonAbrasionModel = new G4WilsonAbrasionModel(); wilsonAbrasionModel->SetMinEnergy(0*MeV); IonInelasticProc->RegisterMe(wilsonAbrasionModel); // applying cross section G4TripathiCrossSection* TripathiCS = new G4TripathiCrossSection; G4IonsShenCrossSection* ShenCS = new G4IonsShenCrossSection; IonInelasticProc->AddDataSet(ShenCS); IonInelasticProc->AddDataSet(TripathiCS); processManager->AddDiscreteProcess(IonInelasticProc); ################################################################################ Possibly you have some hints for me how to manage this. Many thanks and best regards, Jakob Wehner
|
| Inline Depth: | Outline Depth: | Add message: |
|
to: |