|Message: Re: Is there any Hadronic Processes||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
This looks fine as far as EM physics goes. In your energy range you will likely need a hadronic physics process for ions.
You can define a hadronic inelastic process as follows:
G4HadronInelasticProcess* ionInelasticProcess = new G4HadronInelasticProcess("IonInelasticProcess", G4GenericIon::GenericIon());
ionInelasticProcess->RegisterMe(binaryCascade); ionInelasticProcess->RegisterMe(qmd); ionInelasticProcess->AddDataSet(shenXS); ionInelasticProcess->AddDataSet(tripXS); ionInelasticProcess->AddDataSet(tripLightXS);
where binaryCascade and qmd are, respectively, the low energy and high energy ion hadronic models. They can be defined as follows:
binaryCascade = new G4BinaryLightIonReaction; binaryCascade->SetMinEnergy(0.0); binaryCascade->SetMaxEnergy(110*MeV);
qmd = new G4QMDReaction; qmd->SetMinEnergy(100*MeV); qmd->SetMaxEnergy(10*GeV);
The data sets are the hadronic ion-ion cross sections:
G4IonsShenCrossSection* shenXS = new G4IonsShenCrossSection; G4TripathiCrossSection* tripXS = new G4TripathiCrossSection; G4TripathiLightCrossSection* tripLightXS = new G4TripathiLightCrossSection;
|Inline Depth:||Outline Depth:||Add message:|