Message: Assigning regions in a nested parameterisation Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Assigning regions in a nested parameterisation 

Forum: Physics List
Date: 18 Jul, 2016
From: Jake Forster <Jake Forster>

Hello,
 
I would like to apply different physics to different voxels in a nested parameterisation. The nested parameterisation is done following the RE02 example.

In the physics list I have defined default physics as well as physics for the region “Target” (as in the microdosimetry example).

In detector construction I define the region “Target” and apply cuts:
G4Region* fpRegion = new G4Region(“Target”);
G4ProductionCuts* cuts = new G4ProductionCuts();
G4Double defCut = 1*nanometer;
cuts->SetProductionCut(defCut,”gamma”);
cuts->SetProductionCut(defCut,”e-”);
cuts->SetProductionCut(defCut,”e+”);
cuts->SetProductionCut(defCut,”proton”);
fpRegion->SetProductionCuts(cuts);
 
Then in nested phantom parameterisation, in the ComputeMaterial function I add a voxel to the Target region by:
G4Region* aRegion = G4RegionStore::GetInstance()->GetRegion(“Target”);
G4LogicalVolume* logVolume = physVolume->GetLogicalVolume();
aRegion->AddRootLogicalVolume(logVolume);
 
The code compiles and begins to run, but as soon as the particle enters the phantom it segmentation faults.
 
Is there something wrong with this approach? Is there another way to have different physics in different voxels?
 
Thank you.
Kind regards,
Jake

 Add Message Add Message
to: "Assigning regions in a nested parameterisation"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]