Problem 2369

Created attachment 692 [details]
SecondSetup

Hi, 

This bug is reported again after what we believe was an premature resolution to problem 2359. Previously, only the second comment, which was on excessive warning reports were remarked upon. However, problem 2359 had a more pressing issue that revolved around many killed particles when using G4PhantomParameterisation. This error is produced when utilizing G4PhantomParameterisation: 

*** Event Must Be Aborted ***

G4WT1 > G4Track (0x7ff04e07c9b0) - track ID = 4, parent ID = 1

G4WT1 >  Particle type : e- - creator process : RadioactiveDecayBase, creator model : gammaPIXE

G4WT1 >  Kinetic energy : 961.646 keV - Momentum direction : (0.650302,-0.000215437,0.759676)

G4WT1 >  Step length : 1 Ang - total energy deposit : 0 eV 

G4WT1 >  Pre-step point : (0.662544,-0.5,1.47516) - Physical volume : PatientRealBoxPV (Water__1.000 g/cm3)

G4WT1 >  - defined by : ScoringWorld - step status : 3

G4WT1 >  Post-step point : (0.662544,-0.5,1.47516) - Physical volume : PatientRealBoxPV (Water__1.000 g/cm3)

G4WT1 >  - defined by : msc - step status : 7

G4WT1 >  *** Note: Step information might not be properly updated.

G4WT1 > 

-------- EEEE -------- G4Exception-END --------- EEEE -------

To quickly recap, our simulation had two setups, which were both based upon S-value or Dose-Point-Kernel Medical Physics simulations. 

The first setup had a 5 m x 5m x 5m world filled with air, which contains a cubed detector geometry filled with water that had the dimensions 81mm x 81mm x 81mm. This cubed geometry was voxelized by a voxel volume of 1 mm3 in x,y,z with G4PhantomParameterisation. The scoring was achieved by creating a parallelworld with no mass or radiation physics turned on. This world  created the same exact voxelized cubed geometry (utilized G4PhantomParameterisation) in the parallelworld as was used as the scoring volume. 
The primitive scorer G4VPSDoseDeposit was utilized and the radiation physics list used was Penelope with G4DecayPhysics and G4RadioactiveDecayPhysics. A source was placed at the center of the cube with GPS. All deposited dose was scored. This simulation results are indicated in the G4Phant.png that showed results that had many killed particles. 

The second setup was the exact same setup as the first simulation except the cubed detector geometry was not voxelized with G4PhantomParameterisation and was filled entirely with water. A parallel world was again created and voxelized with G4PhantomParameterisation. The result of this simulation is seen in noG4Phant.png and seems to produce expected results. 

Best, 
Peter