|Message: Reversing microdosimetry example||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)|
Dear Geant4 Users,
I would like to adapt the microdosimetry example: instead of having standard physics in the world, and DNA in a region, I need to have DNA in the world and Livermore physics in a region.
However, when I try to set Livermore models for my region and deactivate DNA in this region, Livermore gets used everywhere.
Here is what I wrote for proton ionisation:
ph->RegisterProcess(new G4DNAIonisation("proton_G4DNAIonisation"), particle);
G4hIonisation* hIoni = new G4hIonisation(); hIoni->SetEmModel(new G4DummyModel(),1); hIoni->SetEmModel(new G4DummyModel(), 2); hIoni->SetStepFunction(0.2, 50 * um); ph->RegisterProcess(hIoni, particle);
//For region mod = new G4DNARuddIonisationModel(); mod->SetActivationLowEnergyLimit(100 * MeV); em_config->SetExtraEmModel("proton", "proton_G4DNAIonisation", mod, reg, 0.0, 0.5*MeV);
mod = new G4DNABornIonisationModel(); mod->SetActivationLowEnergyLimit(100 * MeV); em_config->SetExtraEmModel("proton", "proton_G4DNAIonisation", mod, reg, 0.5*MeV, 100 * MeV);
mod = new G4BraggModel(); em_config->SetExtraEmModel("proton", "hIoni", mod, reg, 0.0, 2 * MeV, new G4IonFluctuations());
mod = new G4BetheBlochModel(); em_config->SetExtraEmModel("proton", "hIoni", mod, reg, 2 * MeV, 10 * TeV, new G4UniversalFluctuation());
I can't figure out why it's not working. Does someone have any idea? Thanks !
|Inline Depth:||Outline Depth:||Add message:|