Validation of electron multiple scattering 13 to 20 MeV Daren - - PowerPoint PPT Presentation

validation of electron multiple scattering
SMART_READER_LITE
LIVE PREVIEW

Validation of electron multiple scattering 13 to 20 MeV Daren - - PowerPoint PPT Presentation

Validation of electron multiple scattering 13 to 20 MeV Daren Sawkey September 29, 2015 Geant4 collaboration meeting Experiment: Ross et al., Med. Phys. 35 , 4121 (2008) 7 materials, various thickness, 13 and 20 MeV Simulation:


slide-1
SLIDE 1

Validation of electron multiple scattering

13 to 20 MeV

Daren Sawkey

September 29, 2015 Geant4 collaboration meeting

slide-2
SLIDE 2

Experiment: Ross et al., Med. Phys. 35, 4121 (2008)

7 materials, various thickness, 13 and 20 MeV

Simulation: electronScattering2 example

slide-3
SLIDE 3
  • 1. Width of Gaussian fit to central part
  • 2. Overall goodness of fit
slide-4
SLIDE 4

Improvements:

  • Use physics lists from Geant4 source rather than local
  • same # of incident particles for each physics list (previously SS was 10%)
  • calculation of p-value
  • Include simulation, measurement uncertainties in chi-squared

Simulation uncertainty: Measurement uncertainty: 1% in width—what is it for each data point?

determine empirically

slide-5
SLIDE 5

Width good, chi-sq good

Examples

slide-6
SLIDE 6

Width good, chi-sq poor

Examples

slide-7
SLIDE 7

Width poor, chi-sq poor

Examples

slide-8
SLIDE 8

What is the uncertainty of the measured data points?

slide-9
SLIDE 9

Assume 0% measurement uncertainty

slide-10
SLIDE 10

Assume 1% measurement uncertainty

slide-11
SLIDE 11

Assume 2% measurement uncertainty

slide-12
SLIDE 12

Summarize results

  • Calculate average chi-square
  • Root-mean-square (ratio of widths – 1)
  • Express as %
  • Function of Geant4 version
slide-13
SLIDE 13

20 MeV

ratio of widths

slide-14
SLIDE 14

Chi-squared

20 MeV 20 MeV

slide-15
SLIDE 15

20 MeV

p-value

slide-16
SLIDE 16

IDEAL

All data vs. version: Widths

1-

slide-17
SLIDE 17

IDEAL

All data vs. version: chi-square

slide-18
SLIDE 18

All data vs. version: p-value

slide-19
SLIDE 19

Examples of two problems found

slide-20
SLIDE 20

LateralDisplacementBeyondSafety

  • n
  • ff
slide-21
SLIDE 21

Step Limit

Boundary Safety SafetyPlus

slide-22
SLIDE 22

TestEm5, 10 MeV e- through 10 um Cu

Distance to boundary: msc limits step; strict near boundary SafetyPlus: msc doesn’t limit step in absorber.

slide-23
SLIDE 23

Conclusions:

Rigorous test of scattering algorithms! Good calculation of chi-square to test tails Two (3?) metrics per physics list per Geant4 version Should Chi-square, p-value, or something else be used?