Modelling Algorithms Flashcards

1
Q

Measurement vs modelling approach

A

Measurement: measure absorbed dose in water and correct for inhomogeneities (still used in RadCalc etc)
Modelling: model physics of what is going on, can use pencil beam and collapsed cone to aggregate large no’s of particles, don’t need water dose distribution as pre-requisite

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
2
Q

What is a history?

A

The transport of a single primary particle and all the subsequent secondary particle

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
3
Q

What is a track?

A

The path of an individual particle until it is absorbed. Collection of tracks of primary and secondaries makes history.
Comprised of many steps

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
4
Q

What is a simulation?

A

All of the histories in a given geometry

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
5
Q

What 4 elements characterise the history?

A

Free path between events
Interaction
Energy loss and angular deflection
New particles

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
6
Q

What random number is generated?

A

Between 0 and 1

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
7
Q

How do we get a random number to correspond to a quantity?

A

Sampling using Cumulative probability distribution method

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
8
Q

Steps to get random number

A

Integrate probability distribution to find cumulative probability function
Invert cumulative probability function

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
9
Q

What is the geometry and what defines it?

A

The description of the surfaces and bodies through which the particles are transported
Defined by linac and patient
Patient derived from CT scan

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
10
Q

What is scoring

A

Accumulation of dose in each voxel for each interaction for a single particle

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
11
Q

How is dose dist. built up?

A

Summing contributions in each voxel from large number of particle histories

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
12
Q

What is statistics and how does increasing the number of histories affect it?

A

The dose in each voxel has an uncertainty which is referred to as the statistics, increasing no. of histories reduces the uncertainty but increases computation time

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
13
Q

What models are used to model particles starting from linac target?

A

Phase space model
Patient specific model

(Phase space often divided into:
Beam phase space
Patient modifiers phase space)

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
14
Q

What is phase space?

A

A record consisting of each particles position, energy, direction as it exits treatment head

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
15
Q

What is input for patient specific model?

A

Phase space data, used as basis of dose calculation in specific patient

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
16
Q

What is a condensed electron history?

A

Condense many small interactions into one large virtual reaction, as nearly all electron interactions involve small energy losses and scattering angles.
Each step transfers the same amount of energy as the large number of small losses, and scatters the electron by an angle equal to a large number of scattering angles.
10^5-10^6 electron interactions per electron history - impossible to model all of them.

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
17
Q

How can calculation time be cut down?

A

Variance reduction techniques:
Transport cutoffs
Zonal discard

Decrease calculation time to achieve a given statistical variance which is still clinically acceptable

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
18
Q

What are transport cutoffs?

A

Stop tracking an energy below a threshold and deposit all energy at last position. Photon cut off lower than electron due to range

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
19
Q

What is zonal discard?

A

If electron does not have enough energy to escape the voxel, deposit all energy in that voxel.

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
20
Q

What parameters can you have in raystation?

A

Number of histories per cm^2
Calculation resolution

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
21
Q

When is plan approved?

A

When mean relative statistical uncertainty is less than 2%

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
22
Q

How does number of histories affect uncertainty?

A

More histories, lower the uncertainty, more accurate and less noisy doses.

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
23
Q

How does resolution affect uncertainty?

A

Increasing dose grid voxel size means number of histories per voxel increases, uncertainty decreases but resolution worse.

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
24
Q

What does uncertainty ultimately depend on?

A

Number of histories per voxel

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
25
Q

When is MC modelling crucial?

A

In MR-linacs.
The B field makes the movement of secondary electrons non-isotropic. Electron return effect.

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
26
Q

What is simulation approach

A

Simulates man individual particles and their history with probabilities and interaction processes, most accurate but computationally intensive.

Can use MC to help with model based approaches

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
27
Q

Idea behind Digital Signal Processing

A

Take (usually) analogue data and convert it into digital form.

Uses signals that can originate as sensory data and convert to digital form.

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
28
Q

Signal and system in DSP

A

Signal: real world data or parameter that relates to input signal
System: process that produces and output signal in response to an input signal

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
29
Q

What is superposition and what does it require?

A

Can break down signal into number of smaller simpler components that are processed individually
Requires homogeneity and additivity

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
30
Q

Is dose distribution a linear system?

A

Yes: homogeneity, greater number of photons means greater dose
Additivity: dose from two beams can be added to give the total dose, or dose from individual sites can be added

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
31
Q

Synthesis and decomposition

A

Synthesis is adding things up
Decomposition is breaking things down into constituent parts, can be done an infinite number of ways

32
Q

How can we decompose dose distribution problem?

A

Split into primary fluence and head scatter fluence so these can be adjusted for a variety of factors

33
Q

Input signal and impulse response in dose distribution

A

Input signal is fluence map, modelling primary beam photons and photons scattered in head
Impulse response is scatter kernels

34
Q

What is fluence map?

A

Map of photons passing through medium from head of linac
Shows how beam exits linac and diverges and is attenuated
In inhomogeneous medium this is modified to account for differing densities and effect of attenuation

35
Q

What is scatter kernel and how is it created?

A

Describes how does is deposited around an interaction site by electrons released at the site and photons scattered
Known as PSF/impulse response/ filter or convolution kernel
Built up from modelling many interactions at same site with MC

36
Q

What is the total dose to a point and how are the fluence map and scatter kernal related to it?

A

Total dose is sum of contributions to that point from interactions in all other points
Fluence map dictates how many photons interact at given point and scatter kernel describes how dose will be deposited

37
Q

Convolution for dose

A

Fluence * K = D

38
Q

Convolution theorum

A

Taking the inverse FT of the product of the FTs of fluence and dose kernel gives the convolution
Using fast fourier transform does this quickly

39
Q

Problems with correction based approach?

A

Radiation scatters (photon interacts again, scatter depends on patient internal)
Patient more complex than water tank (inhomogeneities and irregularities)
Fields are the same as measured (can’t measure all irregular fields, treatment fields not likely to be measured)

40
Q

Advantage of model based approach

A

need less measured data (but maybe more non-standard conditions data), good models can predict dose far better in complex situations

41
Q

What are two steps of model based dose calculation algorithm?

A

Model head of linac and what comes out
Model patient and how beam is deposited as energy

42
Q

What is the difference between type A and type B algorithms?

A

Type A does not model the variation of penumbra with density
Type B does

43
Q

What is input for dose calculations?

A

Energy fluence

Defined as 2D array describing how energy of beam is distributed across plane

44
Q

Fluence vs energy fluence

A

Fluence is particles passing per unit ae = number/area

Energy fluence = particles x energy / beam cross section

45
Q

What are 2 energy fluence maps used?

A

Primary energy fluence - direct photons
Head scatter energy fluence - photons that have been scattered at least once in head

Sum arrays for total if needed

46
Q

Direct energy fluence

A

Usually determined at isocentre distance
Source is origin focus
Divergence considered by ISL

47
Q

Head scatter energy fluence

A

Modelled as if created at FF (more complex than this but this is source of most)
Scaled with ISL

48
Q

Model of source / target

A

Oval shape
Fixed in space, does not rotate
Gaussian distribution
Discretised and represented as 2D array of numbers in TPS
Focal spot of electron beam on W target

49
Q

How is direct fluence obtained?

A

Modulating open beam fluence with attenuation from elements of treatment head
Ray trace performed

50
Q

What are large and small field penumbra most effected by?

A

Large field penumbra more affected by FF size

Small field penumbra affected by primary source

Projecting small point like primary, resulting fluence has sharp edges

51
Q

How are HU values converted to mass density?

A

Linear interpolation in a CT to density table

52
Q

Density scaling vs electron density scaling

A

Density scaling: all material assumed to be water, path length differs in water and other materials, use scaling based on mass density, but large discrepancies found for some materials

Electron density scaling more relevant because CS dominates, error much lower

53
Q

Effective density scaling

A

Takes into account pair production, which becomes an issue at higher energies
Modification of electron density
Lower error over larger energy range

54
Q

What is TERMA

A

Total energy released per unit mass
(by photons to secondary particles)
(energy taken away from beam at a point)

Total energy of electrons from PE effect, electrons from CS, photons participating in CS, electrons and positrons from pair production

55
Q

What is needed to calculate TERMA?

A

Attenuation coefficients
Energy spectrum
Gaussians describing penumbra

56
Q

KERMA and SCERMA

A

KERMA: energy transferred to collisions at point
SCERMA: energy transferred to photon scatter at point

57
Q

Arguments for dose to water?

A

Consistent with absorbed dose to water CoP
Caution option for OARs: Dw usually greater than Dm
Clinical experience with Dw

58
Q

Argument for dose to medium

A

More accurate representation of what patient receives
Calculating Dm then converting to Dw involves additional assumptions increasing uncertainty

59
Q

How is dose deposition estimated?

A

Using predetermined energy deposition kernels to describe the energy spread around a primary interaction

60
Q

Types of kernels

A

Pencil kernels: point kernel pre-convoluted over depth dimension
Planar kernels: point kernel pre-convoluted over 2D (slab dose)
Broad beam kernels: point kernel pre-convoluted over 3D (3D dose distribution)

61
Q

Why use different kernel?

A

Speed advantages to using planar kernel if incident radiation only changing in 1D

If beam intensity intentionally varied in 2D, pencil beam more appropriate

If beam fluence changing in more complex way then point kernel should be used and full 3D integration necessary

62
Q

What are the 2 main types of photon dose algorithms

A

Superposition / convolution: collapsed cone based on point kernels

Simple convolution (utilises FFT): pencil beam, pre convolved point dose kernels

63
Q

Type A algorithm

A

Pencil beam

Pre convolve in depth direction
Only need to convolve in 1D, saves computational time
Sacrifice ability to perform density scaling for lateral heterogeneities
Problem with divergence too

64
Q

Absorbed dose from TERMA type A

A

D = T * K
(TERMA convoluted with dose deposition kernel)

65
Q

Type B Algorithm

A

Collapsed cone: convolution / superposition

Composition of patient affects spread of dose around point, can’t be accurately represented by one invariant pre-calculated point kernel

66
Q

How do collapsed cone point kernels work?

A

PSF discretised into cones travelling out from interaction site
Energy emitted in solid angle cone assumed to be transported along cone axis
Therefore adapts to environment
Directions more closely distributed in forward direction where most of energy flows
Balance between accuracy and speed, TPS use hundreds of directions
Can no longer use FFTs to speed up calculation, not true convolution

67
Q

Kernel depth dependence

A

Needed to account for depth hardening and off-axis softening

Mono-energetic PS kernels simulated and build poly-energetic kernels for each depth

68
Q

Type A vs Type B

A

Type A: invariant kernels use fast FFT convolution techniques
Type B: adapt kernels to local environment and are therefore no longer true convolution and cannot use FFT techniques

69
Q

What 5 components are needed for MC method?

A

Random numbers
Sampling
Photon and electron interaction models and probabilities
Geometry
Scoring and statistics

Uses random number generators and the probability distributions of all the interaction processes.

70
Q

How are the following modelled in fluence:
no. particles
position
direction
energy

A

No. Particles: matrix element
Position: matrix element position
Direction: as if they’re coming from respective source to matrix element
Energy: given by beam spectrum common to all elements

71
Q

What two things combine to give us the energy fluence?

A

Beam model
Treatment plan

72
Q

What question determines type A vs type B algorithms?

A

Does the algorithm model the variation of penumbra with density?

73
Q

What are overall steps to get dose?

A

Work out fluence
Use fluence along with attenuation information to get TERMA
Use TERMA and PSF to get dose

74
Q

How is a pencil kernel created?

A

Pre-convolution over the depth direction

75
Q

How do we get TERMA from CT information?

A

HU converted to mass density through linear interpolation of CT to density
Work out effective density from densities and then electron densities
Discretise effective densities into voxels
Use effective density to get average radiological depth
Calculate TERMA using average radiological depth in each voxel and energy fluence

76
Q

How is effective density used to get radiological depth?

A

Effective density scaling
What path length in water results in same attenuation as corresponding path length in material i
l_w = l_i . rho_i/rho_w

77
Q

What are the two obvious problems with pencil beam?

A

Lack of divergence
Heterogeneities: can’t be scaled laterally