Hello,
I have several G4Regions in a geometry and want to activate the DNA physics in one of them. As an example I use the microdosimetry example. Here I have added another volume and region:
G4Box* testSolid = new G4Box("Test",worldSizeX/2, worldSizeY/2,TargetSizeZ/2);
G4LogicalVolume* logicTest = new G4LogicalVolume(testSolid, fpWaterMaterial,"Test");
new G4PVPlacement(0, G4ThreeVector(0,0,TargetSizeZ*2), "Test", logicTest, physiWorld, false, 0);
testRegion = new G4Region("Test");
testRegion->SetProductionCuts(cuts);
testRegion->AddRootLogicalVolume(logicTest);
If I now run the example, I would expect the DNA processes to be active in the target region, but the standard processes to be active in the test region, since the DNA models were only assigned to the target region. However, I only observe the DNA processes in the test region and no standard process. The same happens when I use the G4EmDNAPhysicsActivator instead of the PhysicsList of the example.
G4VModularPhysicsList* physicsList = new G4VModularPhysicsList();
physicsList->RegisterPhysics( new G4EmLivermorePhysics() );
physicsList->RegisterPhysics(new G4EmDNAPhysicsActivator());
G4EmParameters* param = G4EmParameters::Instance();
param->AddDNA("Target", "DNA_Opt4");
What is the reason for this? How can the DNA models be deactivated in the test region?
Geant4 Version: 10.07
Background: I am actually simulating protons in a water phantom with a gold nanoparticle, using the DNA models for water in water and the DNA models for gold in gold for electrons and the Livermore processes for protons. The above example is just to demonstrate the problem.