1 Practice GEANT4 - ML2 Geant4 School, 12-14 October 2009, Catania Italy Dr. C. Andenna DIPIA – ISPESL e Gruppo Collegato Sanità INFN (Roma) claudio.andenna@ispesl.it.

1

Practice GEANT4 - ML2

Geant4 School, 12-14 October 2009, Catania Italy

Dr. C. Andenna

DIPIA – ISPESL e Gruppo Collegato Sanità INFN (Roma)

claudio.andenna@ispesl.itclaudio.andenna@iss.infn.it

Geant4 Low Energy Electromagnetic Physics Working Group

Some aspects of the program

/primaryParticleData/idParticleSource 1

switch (this->idCurrentParticleSource){case idParticleSource_Random:

this->GenerateFromRandom();break;

case idParticleSource_CalculatedPhaseSpace:this->GenerateFromCalculatedPhaseSpace();break;

}

2Geant4 School, 12-14 October 2009, Catania Italy

Geant4 Low Energy Electromagnetic Physics Working Group

Some aspects of the program

/primaryParticleData/idParticleSource 2switch (this->idCurrentParticleSource)

{case idParticleSource_Random:this->GenerateFromRandom();break;case idParticleSource_CalculatedPhaseSpace:this->GenerateFromCalculatedPhaseSpace();break;

}

3Geant4 School, 12-14 October 2009, Catania Italy

Geant4 Low Energy Electromagnetic Physics Working Group

Some aspects of the program Dose calculation and its standard deviationIn each voxel:

4

111

2

2222

2

NNN

d

d

NN

dNd

NN

dd

N

dd

N

iiN

ii

N

iiN

i

id

N

ii

11

22

2

2

22

N

DND

N

ddN

N

dDdD

N

ii

N

ii

dD

N

ii

N

ii

Geant4 School, 12-14 October 2009, Catania Italy

Geant4 Low Energy Electromagnetic Physics Working Group

Some aspects of the program

Dose calculation:

G4double depEnergy= aStep->GetTotalEnergyDeposit(); voxelMass=density*voxels.halfSize.X*voxels.halfSize.Y*voxels.h

alfSize.Z; dose=depEnergy/voxelMass; Voxels.nEvents++; voxels.depEnergy+=dose; voxels.depEnergy2+=dose*dose;

5Geant4 School, 12-14 October 2009, Catania Italy

Geant4 Low Energy Electromagnetic Physics Working Group

Content

6

An overview of the results

Phase space data Phantom voxelization Comparison between experimental and calculated

doses Graphical view

Geant4 School, 12-14 October 2009, Catania Italy

Geant4 Low Energy Electromagnetic Physics Working Group

Phase space data

Type of data saved

7

n X[mm]

Y[mm]

Z[mm]

dx dy dz kinEnergy

[MeV]

Part Type

Prim part type

Prim part

0 18.57 45.06 -831 0.03934 0.095413 0.99455 1.41684 22 11 33

1 57.61 27.74 -831 0.12265 0.058820 0.99004 1.35906 22 11 71

Geant4 School, 12-14 October 2009, Catania Italy

Geant4 Low Energy Electromagnetic Physics Working Group

Phantom voxelization

Type of data saved

8

Phys Volume X[mm]

Y[mm]

Z[mm]

ix iy iz Dose[Gy]

Dose^2[Gy^2]

nEvents

fullWaterPhantomPV

-297 -297 -177 0 0 27 1.78e-09 1.36e-18 9

.... .. .. ...

Geant4 School, 12-14 October 2009, Catania Italy

Geant4 Low Energy Electromagnetic Physics Working Group

Experimental vs Calculated data

Type of data saved

9

% 6 MeV Tutte le misure gaf 10 luglio 2009 campo 3x3n208=2;fh208=[0 10 3.10291e+0710 758 2.73489e+07];% x [mm], y [mm], z [mm], expDose [Gy], Calculated dose [Gy], Calculated dose2 [Gy^2], nEvents, normDose [Gy], normDoseError [Gy]d208=[0 0 -300 0.6644 5.93877e-09 1.01606e-18 490.184275 0.01706470 0 -290 4.3128 1.44254e-07 1.05554e-15 2744.47607 0.9729410 0 -258 2.7057 7.59039e-08 2.80751e-16 652.35523 0.433426….];

Geant4 School, 12-14 October 2009, Catania Italy

Geant4 Low Energy Electromagnetic Physics Working Group

Graphical view

Dedicated programs have to be used to generate the wanted graphs

We use matlab programs

10Geant4 School, 12-14 October 2009, Catania Italy

Geant4 Low Energy Electromagnetic Physics Working Group

Graphical view: examples

Comparison between calculated and measured data

11Geant4 School, 12-14 October 2009, Catania Italy

Geant4 Low Energy Electromagnetic Physics Working Group

Graphical view: 3D

Two examples of a section of the ROG

12Geant4 School, 12-14 October 2009, Catania Italy

13

Thank you !

Should you have any question,

do not hesitate to contact us