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

Feedback Re: Assigning regions in a nested parameterisation 

Forum: Geometry
Re: Question Assigning regions in a nested parameterisation (Jake Forster)
Date: 19 Aug, 2016
From: Gabriele Cosmo <Gabriele Cosmo>

Regions can only be defined and assigned in your detector construction.

Cheers, Gabriele

On Fri, 22 Jul 2016 06:57:52 GMT, Jake Forster wrote:
> 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 &#65533;Target&#65533; (as in the microdosimetry example).
> In detector construction I define the region &#65533;Target&#65533; and apply cuts:
> G4Region* fpRegion = new G4Region(&#65533;Target&#65533;);
> G4ProductionCuts* cuts = new G4ProductionCuts();
> G4Double defCut = 1*nanometer;
> cuts->SetProductionCut(defCut,&#65533;gamma&#65533;);
> cuts->SetProductionCut(defCut,&#65533;e-&#65533;);
> cuts->SetProductionCut(defCut,&#65533;e+&#65533;);
> cuts->SetProductionCut(defCut,&#65533;proton&#65533;);
> 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(&#65533;Target&#65533;);
> 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: "Re: 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 ]