i

UNIVERSIDAD NACIONAL DE SAN AGUSTÍN DE AREQUIPA

FACULTAD DE CIENCIAS NATURALES Y FORMALES
ESCUELA PROFESIONAL DE FÍSICA

MODELADO DEL CABEZAL DE UN ACELERADOR LINEAL CLÍNICO

PARA LA SIMULACIÓN DEL TRATAMIENTO DE CÁNCER DE
PRÓSTATA LOCALIZADO CON FOTONES DE 6 MEV

Tesis presentada por el Bachiller:

Wilson Sullca Chambi

Para optar el título profesional de:

Licenciado en Física

Asesor:
Dr. Vega Ramírez, José Luis Javier

AREQUIPA – PERÚ
2022
 ii

UNIVERSIDAD NACIONAL DE SAN AGUSTÍN DE AREQUIPA

FACULTAD DE CIENCIAS NATURALES Y FORMALES
ESCUELA PROFESIONAL DE FÍSICA

MODELADO DEL CABEZAL DE UN ACELERADOR LINEAL CLÍNICO

PARA LA SIMULACIÓN DEL TRATAMIENTO DE CÁNCER DE
PRÓSTATA LOCALIZADO CON FOTONES DE 6 MEV

Wilson Sullca Chambi

Tesis presentada a la Facultad de Ciencias

Naturales y Formales, de la Universidad
Nacional de San Agustín, para optar el
título profesional de Licenciado en Física.

ASESOR: Dr. José Luis Javier Vega Ramírez.

Arequipa

(2022)
 iii

FICHA CATALOGRÁFICA

Sullca Chambi, Wilson

Modelado del cabezal de un acelerador lineal clínico para la Simulación del

tratamiento de cáncer de próstata localizado con fotones de 6 MeV

Tesis presentada a la Facultad de Ciencias Naturales y Formales, de

la Universidad Nacional de San Agustín, para optar el título profesional de
Licenciado en Física.

ASESOR: Dr. José Luis Javier Vega Ramírez.

1. Cáncer de Próstata 2. Acelerador lineal (Linac) 3. Dosimetría 4.

Porcentaje de dosis en profundidad (PDD) 5. TPR 6. Código PENELOPE
 iv

FICHA DE APROBACIÓN

Modelado del cabezal de un acelerador lineal clínico para la Simulación del tratamiento de cáncer de próstata
localizado con fotones de 6 MeV

Wilson Sullca chambi

Tesis presentada a la Facultad de Ciencias Naturales y Formales, de la UniversidadNacional de San Agustín,

para optar el título profesional de Licenciado en Física.

Fecha de aprobación:

Banca Examinadora

DRA. ROXANA FILOMENA TORRES GUILÉN

Presidente
Facultad de Ciencias Naturales y Formales – Universidad Nacional de San Agustín.

MG. LUDOLFO FREDY CAYLLAHUA QUISPE

Integrante.
Facultad de Ciencias Naturales y Formales – Universidad Nacional de San Agustín.

MG. ELFER VENANCIO ARENAS HERRERA

Accesitario
Facultad de Ciencias Naturales y Formales – Universidad Nacional de San Agustín.
 v

DEDICATORIA

Dedico este trabajo a Dios, mis padres y a mis hermanos por el gran apoyo que
me dan.
 vi

AGRADECIMIENTOS

Primeramente, quisiera agradecer a Dios por la salud.

A mis padres Nicolás, Marcelina y mis hermanos: Quienes siempre me guiaron por el
camino del estudio y me dieron valores que siempre los tendré presente.

A mi asesor Dr. José Luis Javier Vega Ramírez, por sus enseñanzas, dedicatoria y apoyo
incondicional en la realización de este trabajo de tesis.
 vii

ÍNDICE GENERAL

ÍNDICE DE FIGURAS ...................................................................................................................... 9

ÍNDICE DE TABLAS ..................................................................................................................... 10
RESUMEN ...................................................................................................................................... 11
ABSTRACT ..................................................................................................................................... 13
1. INTRODUCCIÓN ....................................................................................................................... 15
2. FUNDAMENTOS TEÓRICOS .................................................................................................. 18
2.1. Cáncer de próstata ..................................................................................................................... 18
2.2. Aceleradores Lineales de Uso Médico (Linac). ........................................................................ 20
2.3. Interacción de los fotones con la materia. ................................................................................. 22
2.3.1. El efecto fotoeléctrico. ........................................................................................................... 22
2.3.2. El efecto Compton.................................................................................................................. 23
2.3.3. Producción de pares ............................................................................................................... 24
2.3.4. Scattering coherente o Rayleigh ............................................................................................ 25
2.3.5. Interacciones Fotonucleares ................................................................................................... 26
2.4. Interacción de la radiación con tejidos biológicos .................................................................... 28
2.5. Unidades dosimétricas básicas. ................................................................................................. 29
2.5.1. Dosis absorbida ...................................................................................................................... 30
2.5.2. Kerma ..................................................................................................................................... 30
2.6. Especiación de la calidad del haz de fotones. ........................................................................... 31
2.6.1. Rendimiento en profundidad porcentual. .............................................................................. 33
2.7. Análisis de histograma dosis-volumen. .................................................................................... 34
2.8. Descripción de un haz de fotones clínico.................................................................................. 35
2.8.1. Región de Build up. ............................................................................................................... 38
2.9. Método Montecarlo. .................................................................................................................. 39
2.9.1. Probabilidad y sorteo de números aleatorios. ........................................................................ 40
2.9.2. Código PENELOPE ............................................................................................................... 41
2.9.1. La estructura de PENELOPE. ................................................................................................ 42
2.10. Sistemas de Planificación de Tratamientos. ............................................................................ 42
3. MATERIALES Y MÉTODOS .................................................................................................... 44
3.1. Conceptos básicos. .................................................................................................................... 44
 viii
3.2. Instrumentación en dosimetría en teleterapia. ......................................................................... 45
3.2.1. Cámaras de ionización gaseosas. ........................................................................................... 46
3.2.2. Electrómetro. .......................................................................................................................... 48
3.2.3. Maniquí o cuba de agua. ........................................................................................................ 48
3.3. Acelerador Lineal Varian. ......................................................................................................... 49
3.4. Geometría del cabezal. .............................................................................................................. 50
3.4.1. Blanco. ................................................................................................................................... 50
3.4.2. Filtro aplanador. ..................................................................................................................... 51
3.5. Estructura del código PENELOPE ........................................................................................... 52
3.6. El archivo de geometría del cabezal Clinac. ............................................................................. 54
3.7. Creación de la geómetra del cabezal. ........................................................................................ 54
3.7.1. Creación de materiales en PENELOPE. ................................................................................ 55
3.7.2. Archivo de entrada. ................................................................................................................ 56
3.8. Verificación de datos experimentales. ...................................................................................... 58
4. RESULTADOS Y DISCUSIONES ............................................................................................. 60
4.1. Obtención del espectro de Fluencia por electrón incidente ...................................................... 60
4.2. Calculo del factor que determina la calidad del haz 𝑻𝑷𝑹𝟐𝟎; 𝟏𝟎 ............................................. 61
4.3. Medición experimental del PDD. ............................................................................................. 62
4.3.1. comparación de TPR .............................................................................................................. 64
4.4. Comparación de los resultados. ................................................................................................ 65
4.4.1. Porcentaje de dosis en profundidad. ...................................................................................... 65
4.5. Simulación de 4 campos para el cáncer de próstata.................................................................. 67
5.CONCLUSIONES ........................................................................................................................ 75
REFERENCIAS BIBLIOGRÁFICAS ............................................................................................. 77

APENDICES…………………………………………………………………………………………………81
 ix

ÍNDICE DE FIGURAS
Figura 2.1: 19
Figura 2.2: 20
Figura 2.3: 22
Figura 2.4: 24
Figura 2.5: 25
Figura 2.6: 26
Figura 2.7: 27
Figura 2.8: 27
Figura 2.9: 28
Figura 2.10: 34
Figura 2.11 35
Figura 2.12: 37
Figura 2.13: 38
Figura 2.14: 40
Figura 2.15: 43
Figura 3.1: 47
Figura 3.2: 49
Figura 3.3: 50
Figura 3.4: 51
Figura 3.5: 52
Figura 3.6: 53
Figura 3.7: 56
Figura 3.8: 57
Figura 3.9: 60
Figura 4.1: 61
Figura 4.2: 62
Figura 4.3: 63
Figura 4.4: 63
Figura 4.5: 66
Figura 4.6: 67
Figura 4.7: 69
Figura 4.8: 69
Figura 4.9: 70
Figura 4.10: 70
Figura 4.11: 72
Figura 4.12: 72
Figura 4.13: 73
Figura 4.14: 73
Figura 4.15: 74
 x

ÍNDICE DE TABLAS

Tabla 3.1: 51
Tabla 4.1: 63
Tabla 4.2: 64
Tabla 4.3: 71
Tabla 4.4: 73
Tabla 4.5: 74
Tabla 4.6: 74
 xi

RESUMEN
El cáncer de próstata, es el cáncer no cutáneo más común y la segunda causa
principal de muerte por cáncer en los hombres en el Perú (INEN, 2021). Uno de los
tratamientos más aplicados en la actualidad para esta enfermedad es la radioterapia. Sin
embargo, la efectividad de este tipo de tratamientos que involucran dosis elevadas de
radiación viene asociada la posibilidad de sufrir los efectos secundarios de la radiación
sobre los tejidos sanos circundantes. Este tipo riesgo de lesión de las células sanas depende
fundamentalmente de la orientación del haz emitido por la unidad de radioterapia y de la
intensidad de la radiación recibida hacia el paciente.

En este trabajo se utiliza el método Monte Carlo con el código PENELOPE para la
simulación de la planificación de los tratamientos en radioterapia, ya que es una poderosa
herramienta para mejorar la eficacia en el cálculo de la distribución de dosis en un
determinado medio frente a los sistemas actuales de tratamiento. Una simulación de este
tipo implica crear con mucho realismo la geometría del cabezal del acelerador, así como
definir los parámetros físicos que rigen el transporte de las partículas. Además, es
imprescindible el conocimiento detallado del espectro de energía de radiación emitido por
el acelerador lineal modelizado, pues de la energía del haz dependen los factores
dosimétricos a cuantificar.

Con el espectro de energía obtenido se procede a hacer una simulación de un fantoma

de un paciente con cáncer de próstata localizado con dosis de 50 𝐺𝑦 en 25 sesiones, a
200 𝑐𝐺𝑦 por día, creando un fantoma con estructuras como la próstata de volumen
3,33 𝑐𝑚3 también se crean los órganos de riesgo que en este caso son el recto, la vejiga
y las cabezas femorales para saber la dosis real que le puede llegar con más realismo, para
la comparación se utiliza un sistema de planificación como es el eclipse sotfwere del
acelerador lineal Varian, utilizando el mismo fantoma y los mismos parámetros como son
los cuatro campos 04 (lateral derecho, anterior, lateral izquierdo y posterior), con tamaño
de campo de 5 𝑐𝑚 𝑥 5 𝑐𝑚 y con pesos iguales en cada campo. Se compararon los
resultados notando que hay más dosis que el que se puede esperar en la próstata con un
margen de 1,8 % de más, al recto encontramos una dosis máxima de variación de 4,13 %,
 xii
la vejiga recibe una dosis de más de 3,32 % y las cabezas femorales reciben una dosis
de 1,35 % de más, esto nos indica que se debe tener en cuenta estos valores que no son
tan elevados, quero si se trabajan con dosis de más de 70 𝐺𝑦 estos órganos de riesgo
pueden ser dañados por la misma radiación y lo que se espera en la radioterapia es poder
curar al paciente y no generar efectos secundarios.
 xiii

ABSTRACT

Prostate cancer is the most common non-skin cancer and the second leading cause
of cancer death in men in Peru (INEN, 2021). One of the most applied treatments today
for this disease is radiotherapy. However, the efficacy of this type of treatment that
involves high doses of radiation is associated with the possibility of suffering the
secondary effects of radiation on surrounding healthy tissues. This type of risk of injury
to healthy cells depends mainly on the orientation of the beam received by the
radiotherapy unit and the intensity of the radiation towards the patient.

In this work, Monte Carlo is used by the PENELOPE code for the simulation of the
planning of radiotherapy treatments, since it is a powerful tool to improve the efficiency
in the calculation of the dose distribution in a certain medium compared to conventional
systems. current treatment. Such a simulation involves creating the geometry of the
accelerator head with great realism, as well as defining the physical parameters that
govern the transport of the particles. In addition, detailed knowledge of the radiation
energy spectrum emitted by the modeled linear accelerator is essential, since the
dosimetric factors to be quantified depend on the energy of the beam.

With the energy spectrum obtained, a simulation of a phantom of a patient with

localized prostate cancer is made with a dose of 50 𝐺𝑦 in 25 sections, at 200 𝑐𝐺𝑦 per
day, creating a phantom with structures such as the prostate of volume 3,33 𝑐𝑚3 the
organs at risk are also created, which in this case are the rectum, the bladder and the
femoral heads to know the real dose that can reach it more realistically, for comparison a
planning system such as is the sotfwere eclipse of the Varian linear accelerator, using the
same phantom and the same parameters as the four fields 04 (right side, front, left side
and back), with a field size of 5 𝑐𝑚 𝑥 5 𝑐𝑚 and with equal weights in each field. The
results were compared, noting that there is more dose than can be expected in the prostate
with a margin of 1,8% more, in the rectum we found a maximum dose variation of 4,13%,
the bladder receives a dose of more of 3,32% and the femoral heads receive a dose of
1,35% more, this tells us that these values should be taken into account, which are not so
high, but if you work with doses of more than 70 Gy, these organs at risk can be damaged
 xiv
by the same radiation and what is expected in radiotherapy is to be able to cure to the
patient and not generate secondary effects.
 15

1. INTRODUCCIÓN

El cáncer desde el punto de vista biológico, es un trastorno, que se caracteriza por

la alteración del equilibrio entre la proliferación y los mecanismos normales de muerte
celular; su consecuencia es el desarrollo de una clona (células) que puede invadir y
destruir los tejidos adyacentes, y diseminarse hacia sitios distantes en los que se forman
nuevas colonias u ocurre propagación metastásica. En la práctica, cáncer es un término
genérico empleado para referirse a más de un centenar de enfermedades distintas, con
epidemiología, origen, factores de riesgo, patrones de diseminación, respuesta al
tratamiento y de diversos pronósticos.

A pesar de los grandes avances en los últimos años, tanto en prevención, como en
detección y tratamiento del cáncer, cada año se diagnostica en hombres 4,000 nuevos
casos de cáncer de próstata (INEN, 2021) y la mayoría de casos de cáncer son detectados
en etapas cuando ya en cáncer está conformado, ya que nuestra población no está
acostumbrada a realizarse chequeos continuos.

La radiación tiene la capacidad para penetrar en la materia e interaccionar con

moléculas, átomos, partículas, en estas interacciones, la radiación pierde parte o la
totalidad de su energía, siendo depositada en el medio material que está atravesando. Esta
transferencia de energía es la causa de los distintos efectos producidos por las radiaciones,
como pueden ser los efectos biológicos sobre la materia viva (Attix F,.2004).

En radioterapia se trata a pacientes generalmente oncológicos, por ello la cantidad

de radiación que va a recibir un paciente debe ser muy cuidadoso en la cantidad de dosis
de radiación administrada. Cuando no se realiza una estimación dosimétrica previa de la
irradiación que se va a aplicar, resulta complicado evaluar con rigor los resultados
obtenidos, como el alcance completo al tumor y el porcentaje de irradiación de tejido
sano. Por ello, la calidad del tratamiento aplicado y la probabilidad de control local
dependen, en gran medida, en la experiencia del especialista. (Pascau et al., 2012).
 16

Hoy en día, el cálculo más preciso de la dosis aplicada mediante radiaciones

ionizantes en el tratamiento médico se basa en la simulación detallada del transporte de
la radiación en medios materiales mediante métodos Monte Carlo (MC) (Krieger et al.,
2005). El empleo de métodos MC para el cálculo de dosis, permite considerar efectos
como la interacción de la radiación con estructuras heterogéneas, así como la
consideración de los efectos de materiales empleados para proteger órganos sanos
(Oshima et al., 2009).

En este trabajo se utilizó el código PENELOPE (PENetration and Energy Loss of

Positrons and Electrons) (Salvat et al., 2003). La forma más habitual de obtener la
información sobre la radiación emitida por el acelerador lineal, y que se codifica en
ficheros denominados espacios de fase, es mediante una simulación MC realista del
acelerador. Ello requiere un conocimiento detallado del cabezal del mismo (Ma and
Jiang., 1999; Capote et al., 2006), lo que no siempre está disponible. Este método lo que
se pretende es simular de forma completa el transporte del haz de radiación a través de
los componentes del cabezal del LINAC y su interacción con un fantoma de agua en el
que se determinarían las distribuciones de dosis a comparar con los datos experimentales
del acelerador., Sin embargo, este procedimiento requiere información detallada acerca
de la geometría y de la composición de los materiales que conforman el cabezal del
LINAC que no siempre está disponible, por ser información reservada que no suelen
proporcionar los fabricantes. Además, el objetivo es determinar las propiedades del haz
de electrones que impacta el blanco, lo cual requiere bastante tiempo de cálculo (Sheikh
et al., 2002).

El objetivo de este trabajo primero es caracterizar el haz de fotones de tal manera

que se pueda tener conocimiento completo de este, y obtener una mejor distribución de
la dosis, así saber exactamente cuanta dosis va recibir un paciente ya que los efectos son
determinísticos (Finlay H et al., 2016). Para la segunda parte con el espectro que se tiene
se hará una simulación en un fantoma que simulara el tratamiento de un cáncer de próstata
 17

localizado (Bainbridge et al., 2017) y así comparar cuanta dosis va recibir el tumor y los
órganos circundantes como la vejiga. En este caso se simulará con 4 campos (lateral
derecho, lateral izquierdo, anterior y posterior) y se hará una comparación con los
planificadores que como eclipse y Mónaco. Para esto en el primer capítulo se trata de la
parte teórica el cual abarca conceptos de fundamentales como el cáncer de próstata y su
tratamiento en radioterapia, también la interacción de la radiación con la materia y
conceptos que se usan con el código PENELOPE. Para el segundo capítulo se muestra
los materiales y métodos usados más que todo para la simulación del cabezal y en la parte
práctica se usa el acelerador lineal de marca VARIAN que se procede a hacer un control
según el TRS 398, y para el último capítulo se muestran los resultados obtenidos de la
simulación y se realiza una simulación de un paciente con cáncer de próstata localizado
en un fantoma creado a dimensión real el cual se le incide cuatro (04) campos de
tratamiento a dosis de 50 Gy en 25 sesiones.
 18

2. FUNDAMENTOS TEÓRICOS

2.1. Cáncer de próstata

La glándula prostática es un órgano reproductor masculino cuya función
principal es secretar líquido prostático, uno de los componentes del semen. Rodea la
uretra y se localiza posterior a la cara inferior de la sínfisis del pubis, superior al
diafragma urogenital y anterior al recto (McNeal JE, 1981). Se encuentra debajo de
la vejiga urinaria y frente al recto. Es un órgano en forma de pirámide con la base
ubicada en la parte superior contigua a la vejiga urinaria y el ápice apuntando hacia
abajo. Mide aproximadamente 5 cm × 4 cm × 3 cm y pesa 20 g entre los 20 y 50 años
y luego aumenta a 30 g entre los 60 y 80 años.

El cáncer de próstata se produce cuando algunas células prostáticas mutan y

comienzan a multiplicarse descontroladamente. Estas también podrían propagarse
desde la próstata a otras partes del cuerpo, especialmente los huesos y los ganglios
linfáticos originando una metástasis. Esta afección puede causar dolor, micción
dificultosa, disfunción eréctil, entre otros síntomas, como se muestra en la figura 2.1.

Figura 2.1: Esquema de la prostata, arriba prostata normal, abajo prostata anormal (Michael J., 2019)
 19

Esta enfermedad se desarrolla frecuentemente en individuos mayores de 50 años.

Es el segundo tipo de cáncer más común en hombres. Sin embargo, muchos hombres que
desarrollan cáncer de próstata generalmente nunca tienen síntomas, ni son sometidos a
terapia. Diversos factores, como la genética y la dieta, han sido implicados en su
desarrollo, pero hasta la fecha, las modalidades de prevención primaria conocidas son
insuficientes para eliminar el riesgo de contraer la enfermedad. La detección se lleva a
cabo principalmente por la prueba en sangre del antígeno prostático específico, llamado
PSA (acrónimo inglés de prostate-specific antigen) o por exploración física de la glándula
prostática (tacto rectal). Los resultados sospechosos típicamente dan lugar a la posterior
extracción de una muestra tisular de la próstata (biopsia), que es examinada en
microscopio. El tratamiento puede incluir cirugía, radioterapia, quimioterapia, o bien una
combinación de todas. La edad y el estado de salud general del afectado, tanto como el
alcance de la diseminación, la apariencia de los tejidos examinados al microscopio y la
respuesta del cáncer al tratamiento inicial, son vitales en la determinación del resultado
terapéutico.

Figura 2.2: Contorno de un cáncer de próstata localizado en un planificador eclipse, se distinguen

los órganos de riesgo (fuente propia de un paciente con cáncer de próstata localizado).
 20

2.2. Aceleradores Lineales de Uso Médico (Linac).

La radioterapia es una técnica que consiste en el uso de radiación de alta energía

para el tratamiento de males como el cáncer. En esta técnica es posible utilizar Rayos-
X, Rayos Gamma y partículas cargadas. La radiación es administrada por medio de un
dispositivo externo al cuerpo (radioterapia externa) o insertando material radioactivo
en el cuerpo, de manera intersticial (braquiterapia).

En nuestro país se están implementando estos equipos para el tratamiento de

estos males, las sesiones de radioterapia externa se efectúan utilizando un acelerador
lineal (figura 2.3), el cual produce Rayos-X de alta energía. Esta radiación se
administra al paciente, en la región del tumor, de manera que las células tumorales
sean destruidas y los tejidos sanos alrededor del tumor se salven.

La producción de radiación por un acelerador lineal está basada en emisión

termo-iónica, de manera similar a lo que ocurre en Rayos-X convencionales. Así, un
haz de electrones se emite de un cañón de electrones y es acelerado en un tubo debido
a microondas de alta frecuencia generadas por un Magnetrón o un Klystron (Pope J.,
1999). Estos electrones pueden ser utilizados directamente para tratamientos
superficiales o bien pueden interaccionar con un blanco de metal de número atómico
elevado para ser transformados en fotones para usarlos en tratamientos de mayor
profundidad en el cuerpo. Por esto, los aceleradores lineales pueden disponer de
distintas energías de fotones y electrones.
 21

Figura 2.3: Esquema general de un acelerador lineal de uso médico (Podgorsak E. B., 2005).

La modulación del haz y su alineación dependen de colimadores de haz de

radiación dentro del cabezal del acelerador lineal. Todos los aceleradores lineales
poseen colimadores, los cuales se pueden dividir entre colimadores primarios y
secundarios. Los colimadores primarios están hechos de bloques de plomo (o de una
aleación de metales con número atómico alto), con el fin de producir un haz cónico
transversal.
 22

2.3. Interacción de los fotones con la materia.

Los fotones cuando interaccionan con algún medio, en principio, pueden
interactuar con los electrones orbitales, el núcleo, o el campo electromagnético de los
átomos que conforman la materia, perdiendo en el proceso toda, parte o nada de su
energía cinética. De los variados modos de interacción posibles, los preponderantes en
física de radiaciones son:

2.3.1. El efecto fotoeléctrico.

Es la interacción entre el fotón y un electrón orbital del átomo, en la cual el

electrón absorbe por completo la energía del fotón incidente, denominándose
fotoelectrón, como se observa en la figura 2.4. Los fotoelectrones no pueden ser
eyectados de un determinado nivel energético del átomo, a excepción de que la energía
del fotón incidente exceda la energía de ligadura del mismo. La expulsión de un
electrón de la capa K es la que más probabilidad de ocurrencia posee, siendo de cuatro
a siete veces más probable que la de la capa L, asumiendo que se satisfacen los
requerimientos energéticos (Cherry et al., 2003). La energía cinética del electrón
eyectado corresponde entonces, a la predicción de la ecuación 2.1 donde 𝐾𝐵 indica la
energía de ligadura correspondiente.

El efecto fotoeléctrico genera, de este modo, una vacancia en el nivel de energía

del electrón eyectado, lo cual lleva a la emisión de un rayo X característico o un
electrón Auger, por lo cual la energía cinética del fotoelectrón es depositada,
típicamente, cerca del sitio en donde ocurrió la interacción.

𝐸𝑓𝑜𝑡𝑜 𝑒𝑙𝑒𝑐𝑡𝑟𝑜𝑛 = ℎ𝜈 − 𝐾𝐵 (2.1)

 23

Figura 2.4: Representación del efecto fotoeléctrico. (Gunderson and Tepper, 2012)

2.3.2. El efecto Compton

Esta la interacción entre el fotón incidente y un electrón orbital de la capa

externa, débilmente ligado al átomo. La energía del fotón (ℎ𝜈) es mucho mayor que
la energía de ligadura del electrón orbital, por lo cual éste cede parte de la misma al
electrón liberado, llamado electrón Compton o de retroceso, y posteriormente es
dispersado un ángulo θ con una energía ℎ𝑣´ (Podgorsak, 2005) como puede
apreciarse en la figura 2.4.

Al plantear la conservación de la energía y del momento se obtiene una

expresión que relaciona la energía del fotón dispersado con el ángulo de dispersión,
como se muestra en la ecuación 2.3, siendo entonces 𝐸𝑘 (ecuación 2.2) la energía
transferida al electrón Compton. Por otro lado, el ángulo de dispersión θ se relaciona
con el ángulo φ del electrón de retroceso mediante la expresión 2.4, siendo así
evidente que el rango de φ es de 0 (para θ = π, retro-dispersión del fotón) a 180 grados
(para θ = 0, dispersión hacia adelante del fotón) dado un valor arbitrario de energía
del fotón.
 24

𝐸𝑘 = ℎ𝜈 − ℎ𝑣´ (2.2)
1
ℎ𝑣´ = ℎ𝜈 ℎ𝜈 (2.3)
1+ (1−cos 𝜗)
𝑚𝑒 𝑐2

ℎ𝜈 𝜗
cot 𝜑 = (1 + 𝑚 2
) tg ( 2 ) (2.4)
𝑒𝑐

Figura 2.5: Representación del efecto Compton (Gunderson and Tepper, 2012).

2.3.3. Producción de pares

Esta interacción se da cuando los fotones interactúan con el campo de fuerza
Coulombiano del átomo. El mismo absorbe por completo al fotón y utiliza su energía
para la creación de un par electrón-positrón: 𝑒 − − 𝑒 + . Dado que la masa en reposo
de un electrón es de 0,511 𝑀𝑒𝑉, este proceso no es posible salvo que la energía del
fotón sea mayor o igual a la masa en reposo de dos electrones, es decir, de 2𝑚𝑐 2 =
1,02 𝑀𝑒𝑉 (Greening, 1985).

La diferencia de energía entre ℎ𝜈 y 1.02 𝑀𝑒𝑉 se reparte como energía cinética

entre las dos partículas originadas, acorde a la expresión 2.5. Tanto el 𝑒 − como el 𝑒 +
disipan su energía mediante radiaciones directamente ionizantes. Debido a que el
positrón pierde toda su energía cinética se aniquilan mutuamente con un electrón
presente en el medio generando un par de fotones, denominados fotones de
aniquilación, que viajan en direcciones opuestas con energía de 0,511 𝑀𝑒𝑉.
 25

𝐸𝑒 + + 𝐸𝑒 − = ℎ𝜈 + 1.02 𝑀𝑒𝑉 (2.5)

Otra opción en este caso es la producción de tripletes, en lugar de pares, es un

proceso posible; sin violar ningún principio físico; pero comparativamente mucho
menos probable a igualdad de condiciones (Valente, 2020).

Figura 2.6: Representación de la interacción producción de pares. (Gunderson and Tepper, 2012)

2.3.4. Scattering coherente o Rayleigh

Es la interacción entre el fotón incidente y el átomo como un todo. Dado que la masa
de este último es naturalmente mayor que la del fotón, luego de la colisión entre ambos, el
átomo absorbe muy poca energía induciendo un leve retroceso dando como resultado un
fotón dispersado sin pérdidas significativas de energía, tal como puede observarse en la
figura 2.6.

Debido a que este tipo de interacción tiene relevancia para energías menores a 50 𝑘𝑒𝑉
y que la transferencia de energía a la materia es muy poca, la cantidad de dispersión coherente
es insignificante en los rangos de energía utilizados.
 26

Figura 2.7: Representación de la interacción scattering de Rayleigh (Gunderson and Tepper, 2012).

2.3.5. Interacciones Fotonucleares

Esta interacción se da cuando fotones muy energéticos (de algunos MeV)

interactúan con el núcleo del átomo, excitándolo y dando como resultado un protón
(𝑝) o un neutrón (𝑛). Las interacciones (𝑛) tienen mucha importancia debido a que los
neutrones pueden ocasionar problemas de radioprotección (energías mayores a
10 𝑀𝑒𝑉). Reacciones fotonucleares de muy alta energía, fuera del rango de interés
para la física médica y para este trabajo, pueden también inducir estados nucleares que
se desexcitan emitiendo partículas más masivas aún, como deuterio, tritio, partículas,
etc. (Hussein, 2007).

Figura 2.8: Representación de interacciones fotonucleares. (Gunderson and Tepper, 2012).

 27

De todas estas interacciones descritas de la radiación ionizante con la materia

expuestos, los primeros tres (efecto fotoeléctrico, Compton y producción de pares) son
los más relevantes en el rango energético de interés para la física médica, ya que
resultan en transferencia de energía a los electrones, los cuales posteriormente la
depositan en el medio mediante interacciones Coulombianas a lo largo de su
trayectoria.

Si analizamos la importancia relativa de estas interacciones en función del número

atómico y de la energía de los fotones incidentes, se obtienen los resultados presentes
en la figura 2.9. El efecto fotoeléctrico predomina a bajas energías para la mayoría de
los materiales, pero disminuye rápidamente a medida que ésta aumenta, en
comparación con el efecto Compton que se vuelve posteriormente predominante. Para
un rango de energías entre 0.1 − 10 𝑀𝑒𝑉 es el mecanismo de interacción que posee
mayor probabilidad de ocurrencia, hasta que los fotones logran energías de algunas
decenas de 𝑀𝑒𝑉 , siendo así la producción de pares el proceso que contribuye
mayoritariamente en las interacciones de los fotones con el medio

Figura 2.9: representación de la dependencia de la energía y del número atómico de los efectos
Compton, fotoeléctrico y producción de pares. (Greening, 1985)
 28

2.4. Interacción de la radiación con tejidos biológicos

La interacción de la radiación con la materia produce excitaciones y ionizaciones

entre los átomos y moléculas, así como también una gran cantidad de electrones
secundarios que, eventualmente, podrán producir ionizaciones adicionales (efecto
cascada). La radiación ionizante es considerada como un agente, que puede llegar a
dar en un daño biológico. Si bien los efectos que ésta produce sobre el organismo son
evidentes, los mecanismos que generan estos daños no lo son.

Como los sistemas biológicos están constituidos mayoritariamente por agua, la

mayor parte de la ionización producida por irradiación ocurre en las moléculas de este
medio. La molécula de agua que ha perdido un electrón se convierte en un ion positivo
altamente reactivo, de corta vida media, que se desexcita rápidamente para producir
un ion de hidrógeno (𝐻 + ) y un radical 𝑂𝐻 sin carga

𝐻2 𝑂 → 𝐻 + + 𝑂𝐻 (2.6)

El radical OH, que tiene un electrón desapareado, es altamente reactivo, por lo

cual tendrá un rango de unos pocos milímetros. Los radicales libres, que son simples
fragmentos de moléculas rotas, son especies extremadamente inestables, por lo cual
tendrá una vida media corta (Wright HA et al.,1985).

Sobre los efectos biológicos, se ha determinado que el principal objetivo de la

radiación es el ácido desoxirribunocleico (ADN), el cual se organiza en cromosomas
dentro del núcleo de la célula (Wright HA et al.,1985). La reacción más importante de
estos radicales es la que se produce con la molécula de ADN. El daño producido puede
ser directo (cuando la radiación ioniza directamente a la molécula de ADN) e indirecto
(debido al efecto de los radicales libres) y dependerá del tipo de radiación incidente y
la estructura y organización de los cromosomas.
 29

Los daños producidos por la radiación incidente y/o los radicales tienen varios
grados de severidad, dependiendo de la capacidad de reparación de la célula. Un daño
de gran importancia es la ruptura de la doble hélice del ADN que resulta muy difícil
de reparar para la célula, los daños de esta magnitud tendrán consecuencias biológicas
significativas, tal como la aberración del cromosoma, mutaciones y muerte celular. En
caso de que se produzcan mutaciones, la información genética será transmitida a las
células descendientes de la célula dañada pudiendo, eventualmente, terminar en cáncer
(Olive PL, 1998)

La resistencia a la radiación de las células depende fuertemente de la etapa del

ciclo celular en la que se encuentre. Cuando ésta se está dividiendo (mitosis), todas las
moléculas de ADN se encuentran desplegadas en el núcleo, por lo cual la sección
eficaz será mucho mayor, aumentando significativamente la sensibilidad a la
radiación. Debido a que las células cancerígenas poseen una gran actividad
reproductiva, la mayoría de ellas son más sensibles a la radiación que el tejido sano.
Esta ventaja terapéutica constituye la razón por la cual la radioterapia representa, en la
mayoría de los casos, una parte fundamental en el tratamiento contra el cáncer.

2.5. Unidades dosimétricas básicas.

La primera unidad utilizada para la medición de dosis fue el Roentgen (R),

establecida en 1928, que es una unidad de exposición (𝑋). La exposición es una medida
de la ionización producida por los fotones en el aire y que es definida como el cociente
entre el valor absoluto de la carga total de iones del mismo signo que son producidos
en el aire (𝑑𝑄) cuando los electrones liberados por ionizaciones en una masa de aire
𝑑𝑚 son detenidos completamente por el medio. La misma es una unidad de exposición.

𝑑𝑄
𝑋 = 𝑑𝑚 (2.7)
 30

En las unidades del sistema internacional (SI), el Roentgen se expresa de la

siguiente manera:

1 𝑅 = 2,58 × 104 𝐶/𝐾𝑔 (2.8)

Esta unidad sólo podía definirse para rayos X, pero se necesita hacerle
correcciones dependientes de la energía y del medio a irradiar para fotones de alta
energía, electrones y otras partículas, por lo cual fue sustituida rápidamente
(Podgorsak E. B., 2005).

2.5.1. Dosis absorbida

Es la cantidad dosis absorbida sirve para medir los efectos biológicos producidos
por todos los tipos de radiación ionizante, incluidas partículas cargadas, no cargadas,
en todos los materiales y energías. Se define como el cociente entre la energía media
𝛥𝐸 depositada por la radiación ionizante en un material de masa 𝜌𝛥𝑉 , donde 𝛥𝑉 es
muy pequeño, es decir:
𝛥𝐸
𝐷 = lim (2.9)
𝛥𝑉→0 𝜌𝛥𝑉

Antiguamente se utilizaba el 𝑟𝑎𝑑 como unidad de dosis, y representaba la

absorción de 100 𝑒𝑟𝑔 de energía por gramo. En el SI la unidad utilizada es el Gray
(Gy) y se define como 1 𝐺𝑦 = 1 𝐽/𝑘𝑔, por lo que 1 𝐺𝑦 = 100 𝑟𝑎𝑑.

2.5.2. Kerma

La cantidad denominada kerma (K), acrónimo de kinetic energy released in

matter, esta definida como el cociente entre la suma de la energía cinética inicial de
todas las partículas ionizantes cargadas 𝑑𝐸𝑡𝑟 liberadas por partículas sin carga en un
material de masa 𝑑𝑚.
 31

𝑑𝐸𝑡𝑟
𝐾 = (2.10)
𝑑𝑚

Su unidad es la misma que para la dosis (𝐽/𝑘𝑔) y también recibe el nombre de

Gy en el SI.
El kerma para un haz de fotones que atraviesa un medio, en un punto es
directamente proporcional a la fluencia de energía de los fotones 𝛹 , y está dada por:

𝜇𝑡𝑟
𝐾 = 𝛹( ) (2.11)
𝜌
𝜇𝑡𝑟
donde ⁄𝜌 es el coeficiente másico de energía trasferida en un medio, promediado
sobre el espectro de fluencia de energía.

2.6. Especiación de la calidad del haz de fotones.

La mayoría de los protocolos de dosimetría, se basan en patrones de kerma en

aire y en patrones de dosis absorbida en agua, han recomendado la razón tejido–
maniquí, 𝑇𝑃𝑅20;10 , es el que define y da característica de la calidad de un haz de
fotones de alta energía (TRS 398, 2006). Para fotones de alta energía producidos por
aceleradores de uso clínico, la calidad 𝑄 del haz se especifica mediante la razón tejido
maniquí, 𝑇𝑃𝑅20;10 es el cociente entre las dosis absorbidas en las profundidades de
20 𝑐𝑚 y 10 𝑐𝑚 en el maniquí de agua, medidas con una SSD constante de 100 𝑐𝑚 y
un tamaño de campo de 10𝑐𝑚 𝑥 10𝑐𝑚.

La característica más importante del índice de calidad del haz, 𝑇𝑃𝑅20;10 , es su

independencia de la contaminación electrónica en el haz incidente. Es también una
medida del coeficiente efectivo de atenuación que describe la disminución
aproximadamente exponencial de una curva de dosis en profundidad de fotones más
allá de la profundidad del máximo de dosis (Greening J., 2005). Como la 𝑇𝑃𝑅20;10 se
obtiene mediante el cociente de las dosis, no se necesita el uso del factor de corrección
 32

por desplazamiento en las dos profundidades cuando se utilizan cámaras cilíndricas.

Además, la 𝑇𝑃𝑅20;10 , en la mayor parte de las configuraciones clínicas, no se ve
afectada por los pequeños errores sistemáticos en el posicionamiento de la cámara en
cada profundidad debido a que la colocación en las dos posiciones se verá afectada de
la misma manera.

En la literatura se han propuesto otros descriptores de la calidad del haz, como

el porcentaje de dosis en profundidad a 10 𝑐𝑚 y la profundidad del 80 % de la dosis
en profundidad (Andreo P., 2000) en donde se hace una revisión de los descriptores de
la calidad de los haces de fotones, basada en una descripción proporcionada por ICRU
(ICRU, 2001). Se debe recalcar, sin embargo, que no hay un índice de calidad del haz
que satisfaga todos los posibles requisitos de ser un índice único para todo el rango de
energías de este Código de Practica y todos los posibles aceleradores usados en
hospitales y laboratorios de calibración. Esto es muy importante ya que los haces
producidos por los aceleradores no clínicos presentes en algunos laboratorios de
calibración, en general no serán idénticos a los de los aceleradores de uso clínico.

El 𝑇𝑃𝑅20;10 también puede obtenerse mediante la relación sencilla (TRS 398, 2006).

𝑇𝑃𝑅20;10 = 1,2661𝑃𝐷𝐷20;10 − 0,0595 (2.12)

Donde 𝑃𝐷𝐷20;10 es el cociente de los porcentajes de dosis en profundidad en las

profundidades de 20 𝑐𝑚 𝑦 10 𝑐𝑚 para un tamaño de campo del 10 𝑐𝑚 × 10 𝑐𝑚,
definido en la superficie del maniquí a una SSD de 100 cm.
 33

Figura 2.10: Configuración experimental para la determinación del índice de calidad del haz
Q (𝑇𝑃𝑅20;10 ). La distancia fuente-cámara (SCD) se mantiene constante a 100 cm y las medidas
se hacen con 10𝑔/𝑐𝑚2 𝑦 20𝑔/𝑐𝑚2 de agua sobre la cámara. El tamaño del campo en la
posición del punto de referencia de la cámara es de 10 𝑐𝑚 × 10 𝑐𝑚. Se puede usar una camara
cilíndrica o una plano-paralela (TRS 398, 2006)

2.6.1. Rendimiento en profundidad porcentual.

Una forma de caracterizar la distribución de dosis a lo largo del eje central del
haz, es normalizar la dosis de cada punto respecto a un valor medido a una dada
profundidad de referencia. Por lo general dicha normalización se hace cuando se
irradia un fantoma de agua. La cantidad rendimiento en profundidad porcentual (PDD
por sus siglas en inglés), se define como el cociente, expresado en porcentaje, de la
dosis a cualquier profundidad 𝑑, respecto de la dosis a una profundidad de referencia
𝑑0 (en general, el punto donde la dosis es máxima) a lo largo del eje central del haz
(figura 2.10). Es decir:

𝐷
𝑃𝐷𝐷 [ %] = 𝐷 𝑑 × 100 (2.13)
𝑑0
 34

Figura 2.11: Rendimiento en profundidad porcentual. 𝑑 es una profundidad cualquiera y 𝑑0 una de

referencia (fuente propia).

2.7. Análisis de histograma dosis-volumen.

Con el incremento de la tecnología en radioterapia, la interpretación de las

curvas de isodosis para el gran número de planos que involucran las placas
tomográficas, se volvió cada vez más compleja debido a la gran cantidad de
información requerida para estos análisis. Una solución a este inconveniente fue
desarrollada por Austin-Seymour et al. en 1986 (Austin S et al., 1986). Esta propuesta
consistía en agrupar la distribución de dosis tridimensional en histogramas de dosis-
volumen (DVH, por sus siglas en inglés). En este trabajo como se simulará la
irradiación de un tumor localizado de próstata, se hará un histograma con el volumen
que recibía una dosis perteneciente a un intervalo específico, lo cual se conoce como
histograma de dosis-volumen diferencial (dDVH).
 35

Otra forma de representar el volumen irradiado es el histograma de dosis-

volumen acumulado, o simplemente histograma dosis-volumen (Drzymala R. et al.,
1991), en el cual se suman los histogramas diferenciales partiendo de una dosis de
interés hasta una dosis máxima, tal como se expresa en la ecuación (2.14).

𝐷𝑉𝐻(𝐷) = ∑ 𝑑𝐷𝑉𝐻(𝑥)𝐷𝑚𝑎𝑥
𝑥=𝐷 (2.14)

Mientras que el histograma diferencial da una idea de cuánto volumen recibe

una determinada dosis, el acumulado muestra el porcentaje de volumen que recibe al
menos una dosis D. En un tratamiento ideal, el histograma acumulado será una función
escalón, en el cual el 100% del volumen recibirá una dosis menor o igual a la máxima
𝐷𝑚𝑎𝑥 y el 0% del volumen recibirá dosis mayores a ésta.

2.8. Descripción de un haz de fotones clínico.

La dosimetría de radiación trata con dos entidades claramente diferentes: una

describe el haz de radiación de fotones en sí mismo en términos del número y energía
de las partículas que constituyen el haz y la otra describe la cantidad de energía que el
haz de fotones puede depositar en un medio determinado como aire, agua o material
biológico.

Un haz de fotones propagándose en vacío está gobernado por la ley de cuadrado

inverso (la fluencia de fotones es inversamente proporcional al cuadrado de la distancia
desde la fuente; un haz de fotones propagándose a través de un fantoma (Núñez G.,
2019) o un paciente, por otro lado, está afectado no solo por la ley de cuadrado inverso
sino también por la atenuación del haz de fotones dentro de dichos medios. Estos
efectos hacen que la deposición de dosis en un fantoma o paciente sea un proceso
complicado y su determinación una tarea difícil.
 36

La medición directa de la distribución de dosis dentro del paciente es

básicamente imposible, sin embargo, para un resultado exitoso del tratamiento del
paciente es requerimiento que dicha distribución de dosis en el volumen irradiado sea
conocida con exactitud y precisión. Esto se logra mediante el uso de varias funciones
que vinculan la dosis en cualquier punto arbitrario dentro del paciente con la dosis
conocida en el punto de calibración o referencia en un fantoma. Estos valores de dosis
son generalmente medidos con detectores de radiación adecuados. En una distribución
de dosis típica pueden identificarse varios puntos o regiones importantes. El haz
ingresa en el paciente en la superficie donde libera cierta dosis superficial 𝐷𝑠 . Debajo
de la superficie la dosis aumenta rápidamente, alcanza un valor máximo en la
profundidad 𝑧𝑚𝑎𝑥 y luego disminuye exponencialmente hasta alcanzar un valor 𝐷𝑒𝑥
en el punto de salida del paciente, denominada dosis de salida.

Figura 2.12: Curva típica de distribución de dosis en profundidad en un paciente, derecha grafico
teórico, izquierda curva de dosis en profundidad. (Podgorsak E., 2005).

La dosis superficial representa contribuciones a la dosis a partir de: fotones

dispersados de colimadores, filtro aplanador y aire, fotones retro-dispersados desde el
paciente y por electrones de alta energía producidos por interacciones de fotones en el
aire y cualquier estructura de protección en las cercanías del paciente.
 37

La zona entre la superficie (𝑧 = 0) y 𝑧𝑚𝑎𝑥 se conoce como región de build-up

o zona de acumulación y resulta del rango relativamente largo de partículas cargadas
secundarias (electrones y positrones) que primero son liberadas en el paciente por
interacciones de fotones y luego depositan su energía cinética en él. Es la región
inmediatamente debajo de la superficie del paciente, la condición de equilibrio de
partículas cargadas (CPE) no existe, sin embargo, conforme la profundidad aumenta,
la condición de equilibrio electrónico llega a su alcance en 𝑧𝑚𝑎𝑥 . Luego de 𝑧𝑚𝑎𝑥 la
dosis disminuye debido a la atenuación del fotón en el paciente, lo que da como
resultado un CPE transitorio en lugar de verdadero condiciones de equilibrio de
partículas cargadas transitorias (TCPE). La profundidad de la dosis máxima 𝑧𝑚𝑎𝑥
debajo de la superficie del paciente depende de la energía del haz y del tamaño del
campo. La segunda a menudo se ignora pues representa solo un efecto menor. La
distribución de dosis a lo largo del eje central del haz proporciona solo parte de la
información requerida para una descripción de dosis precisa dentro del paciente; es
por ello que se trabaja con dichos datos en conjunto con los llamados perfiles de dosis
fuera del eje. Estos perfiles presentan tres regiones a distinguir: central, penumbra y
umbra.

Figura 2.13: Curva típica de distribución de dosis lateral (Podgorsak E., 2005).
 38

La región central representa la porción central del perfil que se extiende desde
el eje central del haz hasta 1.0 − 1,5 𝑐𝑚 (para fotones de 6 MV) de los bordes del
campo geométrico del haz. El tamaño de campo geométrico, indicado por el campo de
luz óptico, usualmente se define como la separación entre los puntos de nivel de dosis
del 50% en el perfil del haz. Para LINACS, la región central del perfil del haz se ve
afectada por la energía de los electrones que impactan el blanco, por el número atómico
del blanco y por el número atómico y forma geométrica del filtro aplanador.

En la penumbra del perfil, la dosis cambia rápidamente y depende del campo

que definen los colimadores, del tamaño finito del punto focal (tamaño de la fuente) y
del desequilibrio electrónico lateral. La caída de dosis alrededor del borde geométrico
del haz tiene forma sigmoidea y se extiende debajo de las mordazas colimadoras hacia
la región de la cola prenumbral, donde hay una pequeña componente de dosis debido
a la transmisión a través de las mordazas.

2.8.1. Región de Build up.

La dosis en profundidad relativa de un haz de fotones de alta energía aumenta

desde la superficie hasta un máximo, a una profundidad que depende de su misma
energía, antes de disminuir de manera casi exponencial. La region entre la superficie
y la profundidad de la dosis máxima es conocida como la región de Build up o región
de acumulación, y el conocimiento de sus características es muy importante a fin de
estimar el riesgo de una dosificación inferior en una lesión superficial. Esta zona en la
curva de dosis en profundidad, comprendida entre la superficie (𝑧 = 0) y la
profundidad de máxima dosis 𝑧𝑚𝑎𝑥 es consecuencia de los relativamente largos
recorridos de las partículas secundarias (cargadas) con alta energía (electrones y
positrones). Estas partículas secundarias originadas dentro del material irradiado por
interacciones de los fotones primarios, viajan una cierta distancia y a lo largo de esta
trayectoria depositan su energía en el material. En la región inmediata a la superficie
 39

de ingreso dentro del material irradiado, no se satisfacen las condiciones de equilibrio

electrónico, y consecuentemente la dosis absorbida es mucho menor al Kerma de
colisión. Sin embargo, a medida que la profundidad 𝑧 aumenta, las condiciones de
equilibrio electrónico son alcanzadas para un cierto valor de profundidad, 𝑧 = 𝑧𝑚𝑎𝑥 .
Más allá de 𝑧𝑚𝑎𝑥 , tanto la dosis como el Kerma de colisión disminuyen debido al
efecto de atenuación de los fotones incidentes por parte del material irradiado,
resultando en una especie de pseudo-equilibrio electrónico.

Figura 2.14: Build up de la curva de dosis en profundidad (Podgorsak E., 2005).

2.9. Método Montecarlo.

El método Monte Carlo (MC) se define como una técnica numérica no determinista
basada en el muestreo aleatorio repetido para obtener el resultado, en esta medida se
constituye como un método extensamente usado para simular sistemas físicos cuando es
imposible tener una solución exacta con un algoritmo determinista. A diferencia de otras
técnicas numéricas, que producen soluciones aproximadas, la ventaja del método MC es
que tiene un error absoluto de la estimación que decrece como ~1/√𝑁, donde N es el
número de historias simuladas, independientemente de la dimensión del problema a
estudiar.
 40

2.9.1. Probabilidad y sorteo de números aleatorios.

La probabilidad que una variable aleatoria este entre un valor 𝑥1 y 𝑥1 + 𝑑𝑥 viene

determinada por:
𝑃 {𝑥1 < 𝑥 < 𝑥1 + 𝑑𝑥 } = 𝑝(𝑥1 )𝑑𝑥 (2,15)

Donde 𝑝(𝑥) es la función de distribución de probabilidad. La función de

distribución de probabilidad acumulada está dada como:

𝑥
𝑃(𝑥) = ∫𝑥 𝑝(𝑥 ′ ) 𝑑𝑥 ′ (2.16)
𝑚𝑖𝑛

Donde 𝑃(𝑥 ′ ) es una función monótona creciente que varía de 0 a 1. Por lo tanto, si
hacemos la transformación 𝜉 = 𝑃(𝑥), la nueva variable aleatoria 𝜉 toma valores en el
intervalo (0,1). La función de distribución de probabilidad de la variable 𝜉 está
relacionada con la de x de la siguiente forma

𝑑𝜉 −1 𝑑𝑃
𝑝𝜉 (𝜉) = 𝑝(𝑥) ( ) = 𝑝(𝑥) =1 (2.17)
𝑑𝑥 𝑑𝑥

De donde se deduce que 𝜉 está distribuida uniformemente en el intervalo (0,1). Así,

sorteando valores uniformemente distribuidos entre 0 y 1, podemos obtener 𝑥 que siguen
la distribución de probabilidad 𝑝(𝑥) invirtiendo la transformación 𝑥 = 𝑃 −1 (𝜉) (Salvat
et al., 2011).

Esto se llama método de la transformada inversa, pero existen otros métodos para
obtener variables aleatorias que sigan cierta distribución de probabilidad a partir de la
generación de números aleatorios (Martínez I., 2008):
 41

2.9.2. Código PENELOPE

PENELOPE (Penetration and Energy Loss of Positrons and Electrons) es un

paquete de subrutinas escritas en fortran 77 que realizan simulaciones en Montecarlo
(MC) de transporte acoplado de electrones y fotones en sistemas de materiales arbitrarios
que constan de una serie de regiones homogéneas (cuerpos) limitadas por interfaces
(Salvat et al., 2011), para un amplio rango de energías (50 𝑒𝑉 a 1 𝐺𝑒𝑉). Inicialmente, y
como su nombre lo indica, el código se utilizaba para el transporte de electrones y
positrones mediante una simulación mixta. La misma consiste en una combinación entre
simulaciones detalladas (computan evento por evento) y técnicas de simulación
condensada (acumulan una serie de eventos en un único evento equivalente
estadísticamente) (Pérez and Valente, 2018).

Dado que PENELOPE se trata de una serie de subrutinas que obedecen a una
secuencia, no puede operar por sí mismo, es el usuario el que debe proveer un programa
principal, main, que direccione y detalle la situación particular a simular. El mismo debe
contener información que controle la geometría, la evolución de las trayectorias, y guarde
registro de las cantidades relevantes. Por este motivo también el usuario se debe encargar
de proporcionar el archivo general de input y los archivos (típicamente de extensión . 𝑔𝑒𝑜
y . 𝑚𝑎𝑡) que se encarguen de dar detalles de la geometría deseada y los materiales a
utilizar, respectivamente.

Una vez finalizada la simulación se generan archivos de salida denominados

outputs, que guardan los resultados generales, estimaciones de las incertezas y tablas que
resumen el proceso de simulación (Pérez and Valente, 2018). Mediante la incorporación
de rutinas definidas por el usuario, PENELOPE permite incluir la presencia de campos
electromagnéticos. Por lo tanto, es posible calcular el transporte de radiación y las
distribuciones de dosis tomando en cuenta los efectos producidos por los campos.
 42

Por último, PENELOPE es propiedad de la NEA (Nuclear Energy Agency),

organismo que lo gestiona y coordina la distribución del paquete en cumplimiento de las
normativas vigentes. A los fines del presente trabajo se utilizó la licencia oficial de
PENELOPE v. 2008 a nombre de Mauro Valente, en calidad de Liaison officer de NEA
para Argentina y la región.

2.9.1. La estructura de PENELOPE.

Estructura general PENELOPE es un paquete de subrutinas; por lo tanto, requiere

de un programa principal para funcionar, además de ciertos datos de entrada. En la Figura
2.14 se ilustra la estructura general de PENELOPE.

Figura 2.15: Esquema de la operación del paquete PENELOPE (Salvat F., 2011).

2.10. Sistemas de Planificación de Tratamientos.

Los sistemas de planificación de tratamientos (TPS) se utilizan en radioterapia

externa para simular la configuración geométrica de los haces de radiación y su
distribución dosimétrica. Con esto se busca maximizar el control tumoral y minimizar
daños a los órganos sanos circundantes. La utilización de sistemas complejos de cálculo
 43

es necesaria para la planificación de tratamientos (Sharpe, M. B., 2006).

La fase previa a la planificación del tratamiento consiste en la adquisición de
imágenes CT, tomando especial cuidado en la reproducibilidad y comodidad del paciente,
ya que esta será la misma durante todo el tratamiento. En algunas ocasiones cuando es
difícil localizar el volumen tumoral, se emplean otras modalidades de imágenes como la
tomografía por emisión de positrones (PET) o imágenes de resonancia magnética (MRI,
por sus siglas en ingles). La siguiente fase consiste en realizar los trazos para delimitar
los blancos anatómicos a tratar por el medico radio-oncólogo (Jones, D., 1994).

GTV (volumen de tumor macroscópico): es el tumor en sí, con su crecimiento

maligno y puede consistir de tumor primario, linfodenopatıas metástasis y metástasis a
distancia.

CTV (Volumen blanco clínico): Es un volumen de tejido que contiene al GTV, más
un margen de seguridad que incluye enfermedad maligna microscópica sub-clínica que
debe ser eliminada (es decir, tejido con sospecha de tumor que no puede demostrarse
clínicamente).

PTV (Volumen blanco de tratamiento): es el volumen que Incluye al CTV más el

margen interno y el margen de preparación, y este ligado al marco de referencia del
acelerador, a este volumen usualmente se da tratamiento.

OAR (órganos de riesgo): son los Tejidos sanos, órganos cerca del tumor, cuya
tolerancia de dosis puede influenciar el volumen a tratar y/o la prescripción de la dosis.
 44

3. MATERIALES Y MÉTODOS

3.1. Conceptos básicos.

A continuación, se da los nombres técnicos que se usa en este trabajo.

 PDD (Percentage Depth Dose): Porcentaje de dosis absorbida a una profundidad

del eje central, normalizada respecto a la dosis a una profundidad de referencia
de dicho eje, manteniendo constante la distancia fuente-superficie y el tamaño
de campo (Pinza C. et al., 2009).

 Perfiles de dosis absorbida: Representan la variación relativa de dosis absorbida

en función de la distancia al eje del haz, normalizados al valor de dosis absorbida
en dicho eje (Pinza C. et al., 2009).

 Penumbra: Dosimétricamente se define como “la distancia lateral existente entre

los puntos que reciben el 80 % y 20 % de la dosis absorbida respecto al eje del
haz” (Pinza C. et al., 2009). Es una región en la que el gradiente de dosis varia
rápidamente con la distancia al eje, por lo que es un factor importante a tener en
cuenta a la hora de planificar un tratamiento.

 Campo irradiado: Es el área situada en el plano perpendicular al eje del haz de

radiación que contiene al isocentro y que delimita la isodosis del 50 % del
máximo de dosis en el plano. El tamaño de campo de los campos regulares
irradiados corresponde a la anchura del perfil de dosis al 50 % (TRS 398; 2006).
En el apartado 3.1 destacaremos la relevancia de los tamaños de campo
pequeños.

 La razón tejido-Fantona TPR (Tissue Phantom Ratio), es el cociente de la dosis

absorbida en agua a una profundidad del eje central, normalizada respecto a la
dosis a una profundidad de referencia en dicho eje, manteniendo constante la
 45

distancia fuente detector y el tamaño de campo (Pinza C. et al., 2009).

 𝑇𝑃𝑅20,10: Es la razón tejido-fantoma en agua a las profundidades de 20 y 10

centímetros obtenidas con 𝑆𝐷𝐷 = 100 𝑐𝑚 y tamaño de campo de 10 𝑥 10 𝑐𝑚2.
Por razones prácticas a la hora de hacer las medidas, se utiliza el 𝑃𝐷𝐷20,10 y
posteriormente la conversión a 𝑇𝑃𝑅20,10 mediante la relación:

𝑇𝑃𝑅20,10 = 1. 2661 𝑃𝐷𝐷20,10 − 0.0595 (3.1)

 Output factors: Factores obtenidos para haces colimados respecto a un haz de

referencia obtenidos como el cociente de las dosis en un punto del eje central del
haz, medidas en las mismas condiciones de SDD y SSD. Estos factores muestran
la variación de la dosis en el punto de referencia cuando cambia el tamaño de
campo.

3.2. Instrumentación en dosimetría en teleterapia.

La magnitud dosimétrica que mejor se correlaciona con los efectos biológicos

producidos por la radiación ionizante es la dosis absorbida (Brosed A. et al., 2011).
Para medir la dosis absorbida que producen los haces de radiación se usa como
material de referencia el agua y diferentes instrumentos, como dosímetros, que
cuantifican los efectos producidos por la radiación ionizante en el agua. Algunos de
los instrumentos utilizados en este trabajo son:
 46

3.2.1. Cámaras de ionización gaseosas.

Son detectores, formados por una cavidad (en este caso cilíndrica) rellena de
gas, recubierta por una pared de grafito o plástico que hace de cátodo, y por un
electrodo de aluminio o grafito, que hace de ánodo. Entre ambos elementos se
genera una tensión que permite recoger las cargas formadas dentro de la cavidad
como consecuencia de la radiación ionizante. El ánodo y el cátodo deben estar
constituidos por materiales de bajo número atómico, equivalente al aire o al tejido,
para para que el proceso de medida interfiera lo mínimo posible en el depósito de
la energía. (Pinza C. et al., 2009)

Figura 3.1: Esquema de una cámara de ionización con sus componentes (Pinza C. et al., 2009).

La dosis absorbida en agua se determina a partir de medidas de ionización en

aire a las que se aplican factores de corrección para tener en cuenta la variación de
las siguientes magnitudes de influencia respecto de las condiciones de calibración
del detector:

 Corrección por temperatura y presión: Debido a que la mayoría de cámaras

están abiertas al aire, es necesario conocer las condiciones atmosféricas en
las que se realizó la medida ya que cambian la masa contenida en el volumen
de ionización. Por ello se utiliza el factor de calibración 𝑘𝑇,𝑃 .
 47

 Efecto de polaridad: El factor de corrección 𝑘𝑝𝑜𝑙 tiene en cuenta la

diferencia de carga recogida por el detector al cambiar la polarización de la
tensión eléctrica.

 Recombinación de iones: El factor de corrección por recombinación 𝑘𝑠 tiene

en cuenta la disminución en la recolección de cargas por la recombinación
entre ellas en el medio ionizante.

Teniendo en cuenta todos estos factores, la dosis absorbida en agua

depositada por un haz de fotones de calidad Q, y medida por una cámara de
ionización calibrada en agua (w) con un haz de calidad de referencia 𝑄0 , vendrá
dada por: (TRS 398; 2006).

𝐷𝑤,𝑄 = 𝑀𝑄 𝑁𝐷,𝑤,𝑄0 𝑘𝑄,𝑄0 (3.2)

Donde:
 𝑀𝑄 es la lectura del dosímetro corregida por las magnitudes de influencia.
Suponiendo que son independientes, se pueden expresar como producto de
los factores de corrección: 𝑀𝑄 = 𝑀𝑄∗ 𝑘𝑇,𝑃 𝑘𝑠 𝑘𝑝𝑜𝑙 , donde 𝑀𝑄∗ es la lectura del
dosímetro sin corregir.
 𝑁𝐷,𝑤,𝑄0 es el coeficiente de calibración del conjunto cámara - electrómetro.
 𝑘𝑄,𝑄0 es un factor que tiene en cuenta la diferencia entre la calidad del haz
medido Q, y la calidad del haz de referencia con el que se calibro la cámara
60
𝑄0 . Normalmente la calidad de referencia es la radiación gamma del 𝐶𝑜.
 48

3.2.2. Electrómetro.

Son un complemento necesario para las cámaras de ionización. Están

formados por una fuente de alimentación que proporciona el voltaje necesario para
generar el campo eléctrico que colecta las cargas, y por un amplificador operacional
de alta ganancia que permite aumentar la amplitud de la señal que se genera en los
dosímetros. Además, permite realizar el cambio de polaridad y medir su efecto.
(Brosed A, 2012).

Figura 3.2: Electrómetro usado para este trabajo, es para la toma de datos, donde se puede
observar que mide solo carga (fuente propia).

3.2.3. Maniquí o cuba de agua.

El agua es el medio en el que se realizan las medidas por su gran

disponibilidad y porque es el principal componente del cuerpo humano. La cuba de
agua es el maniquí más usado, se trata de un recipiente lleno de agua destilada de
grandes dimensiones con guías que permiten mover los detectores de forma remota
en las tres direcciones del espacio. Además, dispone de un software que permite
automatizar la adquisición de datos de forma indexada en el espacio.
 49

Figura 3.3: Cuba de agua usado para este trabajo en la toma de datos, donde se puede observar
la cámara de ionización a una altura de 10 cm sobre la superficie del agua (fuente propia).

3.3. Acelerador Lineal Varian.

Como ya se tiene una descripción general de los componentes de un acelerador

lineal en general, a continuación, nos centraremos en la descripción general del
acelerador lineal, Clinac 2100C que se usó para este trabajo, está situado en el
instituto de oncología y radioterapia de la clínica Ricardo Palma.

Inicialmente tenemos que el haz de electrones golpea en el blanco,

generalmente compuesto de dos capas (algunos son cobre únicamente). La primera
capa está hecha de tungsteno (la mayor parte de la producción de fotones
bremsstrahlung se produce aquí), el segundo es de cobre (para una rápida disipación
de calor).

El tamaño del campo de los haces de rayos X producidos por bremsstrahlung

en el blanco están limitados por el colimador primario el cual está hecho de
tungsteno, luego debajo de este colimador primario se encuentra el filtro aplanador,
el cual va a depender de la energía, este filtro tiene más volumen en medio que en los
bordes exteriores para producir un campo de radiación de intensidad uniforme en la
profundidad. El filtro de 6 MV está hecho de cobre.
 50

3.4. Geometría del cabezal.

En este trabajo se estudia el cabezal del acelerador lineal en modo fotones del
equipo Varian Clinac 2100 (Varian Medical System, Palo Alto, EE.UU.). La
construcción de la geometría de simulación se ha llevado a cabo teniendo en cuenta
la información suministrada por los fabricantes. En la figura 3.4 se esquematizan la
estructura geométrica y los materiales de los componentes que más contribuyen a
la dosis en el maniquí en la configuración de fotones.

Figura 3.4: Componentes de cabezal del acelerador lineal clinac 2100C

3.4.1. Blanco.

El blanco del acelerador Clinac está compuesto de tungsteno, el material

envolvente del blanco es cobre en ambos aceleradores y para las dos
configuraciones de energía. los espesores de la lámina de W son 0.89 mm y 0.64
mm y de la base de Cu 1.57 mm.
 51

Tabla 3.1. Materiales de los elementos que conforman las diferentes

configuraciones del cabezal LINAC estudiado y sus respectivas composiciones y
densidades.

Blanco Materiales W
composición [%] 100
densidad (𝑔/𝑐𝑚3 ) 19,3
Envolvente del blanco Materiales Cu
composición [%] 100
densidad (𝑔/𝑐𝑚3 ) 8,96
Colimadores Materiales W
composición [%] 100
densidad (𝑔/𝑐𝑚3 ) 19,3
Filtro aplanador Materiales Cu
composición [%] 100
densidad (𝑔/𝑐𝑚3 ) 8,96

Figura 3.5: Esquema del colimador primario (Las medidas están en mm) (fuente propia)

3.4.2. Filtro aplanador.

Como se ha indicado anteriormente, la distribución de fotones producidos por

el blanco presenta un pico pronunciado en la dirección adelante. Sin embargo, en
las aplicaciones clínicas se requiere un haz que tenga un frente plano y que sea lo
suficientemente extenso para poder irradiar todo el volumen requerido. Por eso el
haz se hace pasar por un filtro aplanador que cumple esta función.
En nuestro caso para la energía de 6 MV se utiliza un solo filtro aplanador,
 52

donde se aprecia en la figura 3.6 es esquema del filtro aplanador usado para este
trabajo.

Figura 3.6: Filtro aplanador 6 MV. (obtenida de código PENELOPE 2008)

3.5. Estructura del código PENELOPE

Las fuentes usadas básicas del código PENELOPE son las siguientes:

a) penelope.f. Es el paquete de subrutinas para la simulación MC del

transporte de partículas acoplado de fotones y electrones en un medio
homogéneo. Está compuesto de las siguientes rutinas:

 PEINIT. Inicializa el paquete de simulación.

 CLEANS. Inicia una pila donde se guarda el estado de las partículas
generadas por una misma historia.
 START. Inicia el transporte de una partícula.
 SECPAR. Lee el estado de una nueva partícula de la pila. Nos retorna LEFT,
el número de partículas secundarias generadas que todavía se guardan en la
pila.
 JUMP. Determina la distancia hasta la siguiente interacción.
 KNOCK. Simula una interacción de partículas.
 RAND. Genera números aleatorios uniformemente distribuidos entre 0 y 1.
 53

b) pengeom.f. Permite la simulación en cualquier sistema de materiales

formado por diferentes cuerpos homogéneos limitados por superficies
cuadrigas.

 GEOMIN. Lee las variables de entrada del archivo de geometría e

inicializa el paquete de geometría.
 LOCATE. Determina el cuerpo (IBODY) y material (MAT) donde se
mueve la partícula en cada momento. Si 𝑀𝐴𝑇 = 0, la partícula ha
escapado del sistema.
 STEP. Mueve la partícula una cierta distancia DS. Si al mover la partícula
una cierta distancia se encuentra una interfaz (NCROSS), se para la
partícula en la interfaz y se vuelve a empezar la simulación de la partícula
con las condiciones de antes de cruzar la frontera.

c) material.f. Permite crear el archivo de materiales que contiene la

información física de cada material que usaremos en nuestra simulación.
d) penvared.f. Permite aplicar técnicas de reducción de varianza (técnicas
matemáticas para reducir el tiempo en conseguir una cierta incertidumbre
en la simulación).

También se dispone de otras subrutinas/programas en general:

 GVIEW. Programa para visualizar las geometrías.
 SHOWER. Programa para visualizar la trayectoria de las partículas.
 EMFIELDS. Paquete de subrutinas para simular el transporte de partículas
en campos electromagnéticos estáticos.
 TABLES. Genera tablas con los datos físicos del transporte de radiación.
 54

3.6. El archivo de geometría del cabezal Clinac.

Para crear el archivo de geometría, uno de los conceptos básicos a tener en cuenta
es el que una superficie se define por una función 𝛷(𝑥; 𝑦; 𝑧) = 0 (𝛷(𝑥; 𝑦; 𝑧) > 0
para los puntos de fuera de la superficie y 𝛷(𝑥; 𝑦; 𝑧) < 0 para los de dentro).
Generaremos nuestros objetos a partir de las llamadas cuadricas (planos, esferas,
cilindros, paraboloides) que se pueden expresar muy fácilmente en su forma reducida:

𝛷(𝑥; 𝑦; 𝑧) = 𝐼1 𝑥 2 + 𝐼2 𝑦 2 + 𝐼3 𝑧 2 + 𝐼4 𝑧 + 𝐼5 = 0 (3.3)

Donde los coeficientes 𝐼𝑖 solo pueden tomar valores −1, 0 y 1. Aplicando

transformaciones de escalado (X-scale, Y-scale, Z-scale), de rotación (Omega, Theta) y
de traslación (X-shift, Y-shift, Z-shift) podemos obtener cualquier cuadrica en cualquier
posición del espacio que se va a ingresar (Salvat F., 2011).

3.7. Creación de la geómetra del cabezal.

Inicialmente se arma la geometría teniendo en cuenta la ecuación 3.3. En un archivo

con extensión . 𝒈𝒆𝒐, según la información suministrada por el fabricante. La geometría
del cabezal de acelerador lineal se armó tomando varias figuras geométricas y
superponiéndolas, entre estas figuras tenemos los planos, cilindros, conos y esferas. A
continuación, se muestra el archivo 𝒂𝒄𝒆𝒍. 𝒈𝒆𝒐 del acelerador lineal de energía 6 MV.

En este archivo, podemos notar que el acelerador fue encapsulado en un módulo de

aire para hacer mucho más real la simulación, además fue introducido un módulo de agua
a 100 cm de la fuente para simular el fantoma y obtener la curva de dosis en profundidad
en la segunda parte de la simulación. El ángulo dado a las mordazas está dado para
proyectar un campo 10𝑐𝑚 × 10𝑐𝑚. Se tuvieron en cuenta la mayor parte de la geometría
suministrada por el fabricante, en donde se no se tuvo en cuenta los aplicadores ni las
multilaminas ya que en este trabajo eran irrelevantes.
 55

Figura 3.7: Vista en 3D cabezal del acelerador Varian de 6 MV (obtenida del código PENELOPE 2008).

3.7.1. Creación de materiales en PENELOPE.

Siguiendo las especificaciones del paquete Montecarlo entregado por la empresa

Varian, se generaron cada uno de los materiales que componen el cabezal del acelerador
en 𝑝𝑒𝑛𝑑𝑏𝑎𝑠𝑒 con el ejecutable 𝑚𝑎𝑡𝑒𝑟𝑖𝑎𝑙. 𝑒𝑥𝑒.

Figura 3.8: Ejemplo de la generación de materiales con material.exe, en este caso específicamente el
tungsteno.
 56

3.7.2. Archivo de entrada.

En el archivo de entrada con extensión . 𝒊𝒏, en nuestro caso lo llamamos accel1.in

para la primera parte de la simulación en donde interaccionan electrones con el blanco y
en la segunda parte de la simulación llamamos accel2.in la cual ya entraría a la interacción
de fotones con la materia en donde nuestro blanco se convertirá en fuente de fotones. A
continuación, el archivo de entrada de datos iniciales para el acelerador Varian de 6 MV.

.
 57

Podemos notar que en este archivo de entrada encontramos información de la

creación del archivo de espacio de fase, que es fundamental para la segunda etapa de la
simulación, además de la energía inicial de los electrones, posición de la fuente, lectura
de los materiales, y otros parámetros que nos permitirán generar también el archivo de
fluencia.

En la segunda etapa de la simulación encontramos que nuestro blanco será nuestra

fuente de fotones. Este nuevo archivo de entrada accel2.in permitirá generar información
para obtener la curva de dosis en profundidad, la cual es fundamental para los cálculos
del 𝑇𝑃𝑅20;10

Se establecieron los materiales para la representación de la geometría modelada,

según el manual de fabricación del LINAC. El colimador primario del cabezal en su
interior posee el target de tungsteno-cobre.

Se modificó el archivo input original de PENELOPE, adaptando algunos datos para

este trabajo. Este archivo cuenta con distintas secciones a modificar como la definición
de la fuente, información sobre los materiales y parámetros de simulación, distribución
de dosis, etc. Para la definición de la fuente es necesario ingresar el tipo de partículas
primarias, su energía, coordenadas y forma del haz a simular. Esta información debe
representar la medición de la manera más realista posible.
 58

En el acelerador CLINAC 2100C los electrones se aceleran a un potencial nominal

de 6 𝑀𝑉; sin embargo, la energía media del espectro emitido desde la guía hacia el
blanco no necesariamente será 6 MeV, por esta razón en lugar de ingresar una energía
fija de partículas primarias se introdujo una función continua de distribución normal con
energía media de 𝐸̅ = 6,56 𝑀𝑒𝑉 y una desviación estándar de 𝜎 = 0,50 𝑀𝑒𝑉 en 20
divisiones de 0,1 𝑀𝑒𝑉 que va desde 5,5 𝑀𝑒𝑉 hasta 7,5 𝑀𝑒𝑉. Los valores de energía
media y desviación estándar se extrajeron de un estudio previo que consistió en verificar
que la energía de los electrones emitidos desde la guía hacia el blanco no es fija (Salvat
F., 2011).

3.8. Verificación de datos experimentales.

 Antes de la obtención de los datos en el acelerador, se debe verificar que la

posición del gantry fuera de cero grados, también el sistema de colimadores se
debe colocar en la posición de cero grados.

 Los ejes de la cuba se nivelaron para que la cámara de ionización se encontrara

perpendicular al haz de radiación.

 Para el campo de referencia se utilizó un campo de 10 𝑐𝑚 × 10 𝑐𝑚. Para

todos los casos, se verifico que el tamaño de campo fuera el adecuado a una
SSD de 100 𝑐𝑚 sobre la superficie de la cuba de agua.

 La cámara de ionización se colocó en su punto efectivo (1,5 𝑐𝑚) debajo de la

superficie del agua. La verificación de la colocación se realizó por medio de
los láseres laterales y frontal, así como por la proyección del telemetro sobre la
superficie del agua.

 Esta se proyecta cuando se tiene el campo luminoso encendido y el centro de

esta proyección corresponde al eje central del haz de radiación y el campo
 59

luminoso. El barrido en profundidad en el eje central de la cámara de ionización

se realiza colocando la cámara en cada punto de medición escogido
previamente. En este punto

Figura 3.9: Diagrama de cómo se colocó el detector en las distancias de 20 cm y 10 cm

respectivamente en nuestra simulación en PENELOPE (TRS 398, 2006).
 RESULTADOS Y DISCUSIONES
60

4. RESULTADOS Y DISCUSIONES

4.1. Obtención del espectro de Fluencia por electrón incidente

Todas las simulaciones fueron ejecutadas en una computadora instalada,

especialmente para este trabajo. El procesador del computador fue de cuatro núcleos, que
permitió realizar 04 simulaciones en simultáneo, para luego poder hacer un promedio de
las mismas. Es recomendación importante realizar varias simulaciones con distintas
semillas para la obtención de mejores resultados.

En la primera parte de la simulación se obtuvo el espectros de fluencia de fotones

por electrón incidente, uno a la entrada del filtro aplanador y otro a la salida, para un
campo de 10𝑐𝑚 × 10𝑐𝑚, nuestro detector esta hecho específicamente de aire y es un
cilindro de 1,5 cm de radio por 0,1 de espesor. a continuación, se muestra el espectro de
nuestro cabezal VARIAN que se simulo.

Figura 4.1: Espectros de fluencia de fotones por electrón incidente, Energía 6 MV VARIAN (Clinac
2100C). (la tabla se muestra en el apéndice A)
 RESULTADOS Y DISCUSIONES
61

4.2. Calculo del factor que determina la calidad del haz 𝑻𝑷𝑹𝟐𝟎;𝟏𝟎

Con el código PENELOPE 2008 se empieza con la simulación, utilizamos el haz

de 6 𝑀𝑉 con un campo 10𝑐𝑚 × 10𝑐𝑚, como se muestra en la figura 4.3 y a una
distancia fuente superficie de 100 𝑐𝑚 como es recomendada por el Protocolo TRS 398.
Con un fantoma de agua de dimensiones 50𝑥50𝑥50 𝑐𝑚3 Se utilizó una cámara de
ionización de 0,6𝑐𝑚 de diámetro (cilindro) colocados a 1,50 𝑐𝑚, 10,00 𝑐𝑚 y a
20,00 𝑐𝑚, se halló la dosis referencial en cada punto (66,9835 𝐺𝑦 y 39,5872𝑔𝑦).
Utilizamos la ecuación dada en la teoría (3.1), en donde solo nos queda tomar el
𝑃𝐷𝐷20;10 que es la razón de la dosis a 20 𝑐𝑚 𝑦 10 𝑐𝑚 de profundidad en la curva de
dosis en profundidad. Por lo que se interpolaron ciertos datos para encontrar la dosis a
esas profundidades cuyo valor nos da 𝑃𝐷𝐷20;10 = 0,591 tenemos que:

𝑇𝑃𝑅20;10 = 1,2661. 𝑃𝐷𝐷20;10 − 0,0595 = 0, 688

Figura 4.2: Fantoma de agua vista 3D y 2D usado en la simulación del código PENELOPE (fantoma
de dimensiones 50𝑥50𝑥50𝑐𝑚3 ).
 RESULTADOS Y DISCUSIONES
62

Figura 4.3: Espectro de distribución de dosis en el fantoma de agua, el tamaño de campo es de

10𝑐𝑚 × 10𝑐𝑚.

4.3. Medición experimental del PDD.

Las mediciones fueron realizadas en el instituto peruano de oncología y

radioterapia, usando un acelerador lineal Varian 2100 (como se muestra en la figura
4.2).

Figura 4.4: Acelerador lineal Clinac usado para este trabajo de marca Varian.
 RESULTADOS Y DISCUSIONES
63

La figura 4.4 muestra el acelerador lineal Clinac 2100c de la marca VARIAN,

donde puede emitir fotones y electrones respectivamente, aquí se realiza mensualmente
el control de la calidad del haz tanto para electrones y fotones, esta parte es fundamental
para los procesos clínicos ya que se debe entregar las dosis precisas a cada plan de
tratamiento y se espera que cada plan de tratamiento sea igual al simulado en el
planificador y al que se administra a cada paciente.

Para la parte experimental los datos del PDD se obtuvieron utilizando una cámara
ionización marca PTW-freiburg modelo TN 30013 con número de serie 01340, tipo
cilíndrica para mediciones de fotones. las lecturas de los datos se tomaron con el
Electrómetro marca PTW-freiburg serie 00242 con registro de calibración número R339-
19, con factor de calibración de 95,16 𝑚𝐺𝑦/𝑛𝐶. La cámara de ionización que se colocó
en el eje central del campo, de la cuba lleno de agua. Los movimientos en los ejes se
manipularon mediante el controlador que tiene la cuba SNC 1D Scanner, de la marca
SUN NUCLEAR. La cámara de ionización fue calibrada en el Laboratorio Secundario de
Calibración Dosimétrica del Instituto Peruano de Energía Nuclear (IPEN). Las
mediciones se realizaron siguiendo las recomendaciones del código de practica TRS 398,
las cuales se enlistan en la Tabla 4.1.

Tabla 4.1: Condiciones de referencia para la determinación de la dosis absorbida

en agua para haces de fotones de alta energía (TRS 398, 2005).

Magnitud de influencia Valor de referencia

Material del maniquí Agua
Tipo de cámara Cilíndrica
Profundidad de medida 𝑇𝑃𝑅20;10 𝑎 10 𝑔/𝑐𝑚2
𝑧𝑟𝑒𝑓 𝑇𝑃𝑅20;10 𝑎 20 𝑔/𝑐𝑚2
Punto de referencia de la Eje central en el centro del
cámara volumen de la cavidad
Posición del punto de En la profundidad de medida 𝑧𝑟𝑒𝑓
referencia de la cámara
SSD 100 𝑐𝑚
Tamaño de campo 10 𝑐𝑚 × 10 𝑐𝑚
 RESULTADOS Y DISCUSIONES
64

Los valores obtenidos en el electrómetro fueron los siguientes:

Tabla 4.2: Medidas experimentales tomadas del electrómetro PTV de las medidas
de la carga en los puntos de 1.5 cm, 10.0 cm y 20.0 cm respectivamente.

Distancia (𝒄𝒎) Medida (𝒏𝑪)

1,50 19,76
10,00 13,39
20,00 7,81

 Entonces el PDD que es la relación de las dosis en profundidad de 10 𝑦 20 𝑐𝑚

respectivamente sería un valor de 0,583 (dato obtenido de la tabla 4.2 relación
entre las cargas a 10 cm y 20 cm respectivamente).

 El TPR saldría (según la ecuación 3.1) un valor de 0,679.

4.3.1. comparación de TPR

Dados los resultados obtenidos se tiene el error:

𝐸𝑡𝑒𝑜𝑟𝑖𝑐𝑜 − 𝐸𝑚𝑒𝑑𝑖𝑑𝑜
𝐸𝑟 =
𝐸𝑡𝑒𝑜𝑟𝑖𝑐𝑜

0,679 − 0,688
𝐸𝑟 = 𝑥 100%
0,679

𝐸𝑟 = 1,32%

Una discrepancia del 1,32 % que se da, la concordancia entre los resultados experimentales y
simulados.
 RESULTADOS Y DISCUSIONES
65

4.4. Comparación de los resultados.

4.4.1. Porcentaje de dosis en profundidad.

Se observa que para el PDD calculado por simulación, el máximo porcentaje de

dosis absorbida se encuentra a 1,48 𝑐𝑚 de profundidad y a 10,00 𝑐𝑚 le cae el 67,9 %
de la dosis absorbida (profundidad a la cual el protocolo TRS 398 recomienda medir).

De la misma manera, realizando los cálculos con la cámara de ionización, para la

parte experimenta, el perfil PDD es el máximo porcentaje de dosis absorbida se localiza
a 1,50 𝑐𝑚 de profundidad y el 67,6% de dicha dosis se encuentra a 10,00 𝑐𝑚. La
comparación gráfica entre ambos perfiles PDD se muestra en la figura 4.5.

Figura 4.5: Comparación entre el porcentaje de dosis en profundidad experimental y simulado, para un
campo de 10𝑥10 𝑐𝑚2 .
 RESULTADOS Y DISCUSIONES
66

Puede observarse la correspondencia entre los valores simulados por código

PENELOPE y valores experimentales, estos valores se obtuvieron tomando como la
superficie como valor inicial y a partir de este se tomó los datos cada 0,5 cm hasta llegar
a una profundidad de 36 cm, esto comprueba que el cabezal de nuestro acelerador
simulado da valores semejantes a los del acelerador lineal que se tomaron los datos, y
también comprueba que los datos dosimétricos con el planificador eclipse dan los valores
correctos en la planificación de los tratamientos.

También se puede simular los perfiles del campo de 10𝑥10 𝑐𝑚2 a medidos a
diferentes profundidades como se muestra en la figura siguiente:

Figura 4.6: Perfiles de la dosis de entrada de un campo de 10𝑥10 𝑐𝑚2 dónde: a) representa la
medida de la dosis a 1.5 cm de profundidad, b) es la medida a 10 cm de profundidad y c) es a 20 cm
de profundidad.

Los perfiles de campo nos muestras la dosis que se entrega el cabezal siendo un
tamaño de campo constante pero las medidas de la dosis varían según la profundidad
medida. Se aprecia en la figura 4.6 en la parte a) que la dosis es homogénea ya que se
trabaja con una energía de 6 𝑀𝑉 y el dmax para esta energía en agua es a 1,5 cm de
profundidad, esto indica que esta energía entrega su dosis optima a esta profundidad. En
la parte b) que es la dosis a 10 cm de profundidad, se aprecia que la curvatura empieza a
pronunciarse, esto debido al build up y a los electrones secundarios que se produce. Para
la parte c) es el perfil de dosis medida a 20 cm de profundidad, se aprecia que hay más
radiación dispersa y más build up en zonas fuera del campo de tratamiento.
RESULTADOS Y DISCUSIONES
67

Según esto se puede ver que a profundidades altas la energía de 6 MV no será muy
conveniente en su uso, se puede usar otra de mayor energía. Esto puede relacionarse
cuando se trata a pacientes con más masa corporal ya que se tiene mucha más dispersión
y creación de electrones secundarios que pueden generar daño cuando se trate a
pacientes.

4.5. Simulación de 4 campos para el cáncer de próstata.

Una vez que tenemos la geometría de nuestro acelerador lineal y el espectro

de nuestro cabezal, se procede a hacer una simulación de un paciente que es
diagnosticado con cáncer de próstata localizado cuya dosis convencional será de
50 𝐺𝑦 en 25 secciones, en cada sección de tratamiento se dará 200 𝑐𝐺𝑦 de dosis
(Nachbar, M. et al. 2020) para ello se creara un fantoma en forma de elipse que
simulara a la sección pélvica de un paciente varón, con lados en 𝑦 = 20 𝑐𝑚 y en
𝑥 = 36 𝑐𝑚, en este caso se hará un contorno de una próstata en forma esférica de
3,33 𝑐𝑚3 , la vejiga en forma elíptica de volumen estándar de 4,74 𝑐𝑚3 , el recto que
es de forma cilíndrica con un volumen de 5,45 𝑐𝑚3, también se simulo las cabezas
femorales en formas esféricas de volumen estándar de 8,43 𝑐𝑚3 , estos son valores
estándar de un paciente hombre de 60 años de edad con diagnóstico de cáncer de
próstata localizado donde se tratará directamente a la próstata que es en forma
esférica para nuestros cálculos.

Una vez tenido el fantoma con los órganos de riego se precede a simular cuatro
(04) campos de tratamiento (lateral derecho, lateral izquierdo, anterior y posterior),
con tamaño de campo de 5 𝑥 5 𝑐𝑚 en todos los campos como se muestra en la figura
4.5. (el algoritmo de la entrada está dada en el apéndice).
RESULTADOS Y DISCUSIONES
68

Figura 4.7: Vista Transversal 2D del fantoma usado para la simulación en PENELOPE.

Figura 4.8: Vista Sagital 2D del fantoma usado para la simulación en PENELOPE.
RESULTADOS Y DISCUSIONES
69

Figura 4.9: Vista 3D del fantoma usado con loa campos de tratamiento.

En la figura de 4.7 a figura 4.9 se muestra como ingresa los haces en nuestra
simulación de en nuestro fantoma, los cuatro campos de tratamientos. Según nuestros
datos ingresados cada campo tiene el mismo peso y el mismo tamaño de campo.

Figura 4.10: Porcentaje de dosis en profundidad de los campos de tratamiento.

 RESULTADOS Y DISCUSIONES
70

Se puede apreciar como son los PDDs de cada campo de tratamiento, la dosis
máxima claramente llega al centro de nuestro fantoma, también se aprecia la dispersión
que da cada campo de tratamiento.

Los valores obtenidos de nuestra simulación de la dosis absorbida y de cuanto

le llega a cada órgano se muestra en la tabla 4.3 siguiente:

Tabla 4.3: valores calculados en con la simulación en código PENELOPE.

Estructura Volumen Dosis mínima Dosis Máxima Dosis media
(cm3) (cGy) (cGy) (cGy)
Vejiga 47,495 1774,0 4804,3 2837,0
Fantoma 30336,015 0,5 5175,0 124,4
Próstata 34,069 4829,9 5208,3 5146,2
Recto 54,565 156,0 3911,3 1895,0
Cab fem der. 84,378 650,7 3294,1 2676,3
Cab fem izq. 84,696 658,2 3303,3 2691,4

Estos datos se obtuvieron interceptando la dosis de cada campo y en este caso solo
nos interesa la dosis máxima que va a recibir cada órgano de riesgo, ya que se puede
obtener una curva de cada campo más detallado, pero para este trabajo solo se contó con
hacer la intersección de cada campo y ver cuanta dosis le llega a nuestro fantoma con los
órganos de riesgo, se puede simular más órganos de riesgo según el diagnóstico de cada
paciente y que estadio estaría en tumor, pero en este caso solo es un cáncer localizado que
no se expandió aun y que puede ser controlado por radiación.

Para comprobar los datos simulados en código PENELOPE, se simulo en el

planificador usado para la referencia es eclipse con Modelo de cálculo AAA_13335,
donde se usa para tratamientos que se dan a pacientes en la clínica. Con las mismas
características dadas en la simulación anterior, se toma el mismo fantoma con la próstata
y los órganos de riego, ya que puede sacar por un archivo en PENELOPE 2008. La dosis
es similar 50 𝐺𝑦 en 25 sesiones a 200 𝑐𝐺𝑦 por sesión, se usa cuatro campos de
 RESULTADOS Y DISCUSIONES
71

tratamiento lateral derecho (270°), anterior (0°), lateral izquierdo (90°) y posterior (180°)
el tamaño de campo 5 𝑐𝑚 𝑥 5 𝑐𝑚, los pesos de los campos iguales, la prescripción de
dosis es sobre la próstata, los valores obtenidos son dados en la tabla 4.4.

Figura 4.11: Vista de la entrada de la dosis en el fantoma (Planificador eclipse).

Figura 4.12: Vista de la entrada de los campos de tratamiento en el fantoma (Planificador eclipse).
RESULTADOS Y DISCUSIONES
72

Figura 4.13: vista transversal del fantoma y de la planificación en eclipse.

De la figura 4.11 a la figura 4.13 se muestra la entrada de la dosis de cada campo

de tratamiento dados en el planificador eclipse.

Figura 4.14: porcentaje de la distribución de los campos de tratamiento en la planificación en

eclipse.

El eclipse al momento de planificar el tratamiento, nos da diferentes valores de

cuanta dosis por volumen puede recibir una región en específico, por la figura 4.14 se
tiene que la dosis que va a recibir en el centro de la próstata se muestra en la tabla 4.4.
RESULTADOS Y DISCUSIONES
73

Tabla 4.4: Valores calculados con el planificador eclipse de la próstata y órganos

de riesgo.

El histograma de volumen dosis que nos muestra el planificador se muestra en

la figura 4.6 con una dosis de 50,0 𝐺𝑦 que va a recibir el paciente.

Figura 4.15: Histograma volumen dosis de la irradiación sobre la próstata (imagen obtenida del
planificador eclipse).

El histograma de dosis volumen nos ayuda a ver cuánta dosis va recibiendo cada
estructura da planificación, en ella se muestra como la radiación ingresa y como va
cayendo durante el tratamiento, depende también de cuanta dosis se va a entregar al
tratamiento del paciente. En este caso son campos directos y sin uso de multilamina, solo
se usó las mordazas primarias para dar forma de nuestros campos.
RESULTADOS Y DISCUSIONES
74

Comparando los resultados los simulados con los obtenidos en el planificador se

tiene la siguiente tabla de comparación.

Tabla 4.5: comparación de las dosis recibidas tanto la próstata como los órganos de
riego.
SIMULADO PLANIFICADO
Dosis Dosis Dosis Dosis
ESTRUCTURA Máxima media Máxima media
(cGy) (cGy) (cGy) (cGy)
Vejiga 4804,3 2837,0 4963,9 2819,3
Fantoma 5175,0 124,4 5170,0 109,4
Próstata 5208,3 5146,2 5112,0 5070,3
Recto 3911,3 1895,0 3749,7 1883,9
Cab fem der. 3294,1 2676,3 3249,7 2727,3
Cab fem izq. 3303,3 2691,4 3260,3 2741,0

A continuacion se muesta el porcentaje de error que tiene cada estructura

tomando como referencia a la simulacion como dato teorico y a lo planificado en el
planificador eclipse.

Tabla 4.6: Tabla de la deviación de la comparación de la próstata y los órganos

de riesgo.

DESVIACION
Dosis Máxima Dosis media
ESTRUCTURA (%) (%)
Vejiga 3,32 0,62
Próstata 1,85 1,47
Recto 4,13 0,59
Cab fem der. 1,35 1,91
Cab fem izq. 1,30 1,84
Se puede observar un ligero desviacion de las dosis que entrega nuestro
cabezal que se simulo con las mismas condiciones que el planificador. Con estas
dosis se puede concluir las dosis dadas son de gran importancia debido a que pueden
probocar daños a los organos de riesgo.
CONCLUSIONES
75

5. CONCLUSIONES

Se presenta la simulación del cabezal del acelerador lineal Clinac 2100c por el
método Monte Carlo en código PENELOPE, este acelerador es utilizado en el instituto
peruano de oncología y radioterapia de lima para dar tratamientos de radioterapia a
pacientes oncológicos. Se ha desarrollado la geometría del cabezal del acelerador,
configurado para un tratamiento de radioterapia con fotones de energía 6 MV, se ha
simulado con el mayor nivel de detalle posible desde el blanco que es el generador de
fotones por radiación de frenado, el colimador primario junto con el filtro aplanador,
la cámara de monitorización, las mordazas y los colimadores. Para esto, se obtuvo la
caracterización del haz de fotones con su respectivo espectro, donde se nota un mayor
número de fluencia en el espectro de salida, esto físicamente es debido a la generación
de partículas segundarias en el filtro aplanador. Al tomar la curva de dosis en
profundidad para 6 MV los simulados y compararla con los datos obtenidos en una
medición con el acelerador, notamos que son similares y tienen pocas discrepancias.

Como el 𝑇𝑃𝑅20;10 es considerado uno de los factores más importantes que

determina la calidad del haz de un acelerador lineal (TRS 398, 2006), al hacer los
cálculos de este factor sobre la curva de dosis en profundidad PDD encontramos que
el valor obtenido fue de 0.591, el cual es un valor aceptable dentro de los protocolos
establecidos para este factor. El 𝑇𝑃𝑅20;10 simulado en código PENELOPE versión
2008, nos dio un valor de 0,688 y el valor medico experimentalmente con la cámara
de ionización nos da un valor de 0,679 el valor del error en estos valores fue de
1,32 %, que es un valor muy aceptable debido a las fluencia de fotones que presenta
el acelerador y está dentro de las tolerancias permitidas.

Con el espectro de nuestro acelerador lineal se procedió a simular un tratamiento

de un paciente con cáncer de próstata localizado, utilizando cuatro (04) campos de
tratamiento, y con sus órganos de riesgo los cuales se encontró las dosis que le llega a
la próstata (5208,3 𝑐𝐺𝑦) y a los órganos de riesgo como el recto (3911,3 𝑐𝐺𝑦), la
CONCLUSIONES
76

vejiga (4804,3 𝑐𝐺𝑦), cabeza femoral derecha (3294,1 𝑐𝐺𝑦) y cabeza femoral
izquierda (3303,3 𝑐𝐺𝑦). Ahora con el planificador eclipse que usa el modelo de
cálculo AAA 13435, y usando el mismo fantoma con la próstata con sus órganos de
riesgo, se puede notar que hay ligeros variaciones como la dosis máxima que recibe la
próstata difiere en 1,85% que está dentro de la tolerancia, eso indica que nuestro
cabezal esta aceptable. Sobre los órganos de riesgo como el recto que tiene un margen
de 4,13 % nos indica que recibe una dosis más elevada con entrega el planificador. La
vejiga tiene un 3,33 % de margen de error en dosis máxima y la dosis máxima a las
cabezas femorales el margen fue de 1.35% en promedio. Como se tiene entendido
algunos planificadores no toman en cuenta las densidades de algunos órganos como el
hueso que es denso que el agua y también por hacer rápida una simulación por el coste
computacional que trata de ser lo más rápido en tiempo posible, esto hace que no se
sepa con precisión la dosis que le va a llegar a los diferentes órganos, ya que es crucial
cuando se trabaja con dosis muy elevadas y las repercusiones que podía llegar al
paciente en algunos casos sería fatal.

Finalmente, se concluye que el uso del código PENELOPE, permite evaluar

diversas condiciones de irradiación de un paciente y poder medir la distribución de las
dosis que puede recibir, y así se puede saber con más precisión cuanta dosis puede
recibir cada órgano que pudiese ser tratado.
 REFERENCIA BIBLIOGRÁFICA 77

REFERENCIAS BIBLIOGRÁFICAS

Attix, F. H. (2004): Introduction to Radiological Physics ad Rdiation Dosimetry.

Madison, wisconsin: WILEY-VCH Verlag GmbH.
Andreo, P. (2000): On the beam quality specification of high-energy photons for
radiotherapy dosimetry. 27, p. 441–444
Azcona, J. (2007): Desarrollo de un sistema de verificación de tratamientos en
radioterapia con modulación de intensidad. Navarra, España, Universidad de
Navarra, Tesis de Doctorado.
Bainbridge, H., Menten, M., Fast, M., Nill, S., Oelfke, U., and McDonald, F. (2017).
Treating locally advanced lung cancer with a 1.5 T MR-Linac - Effects os the
magnetic field and irradiation geometry on conventionally fractionated and isotoxic
dose-escalated radiotherapy. Radiotherapy and Oncology, 129:280–285.
Brahme, H. (1979): Radiation beam characteristics of a 22 MeV microtron. 18, p. 244–
272
Brahme, P. (1986): Dosimetry and quality specification of high energy photon beams. 25,
p. 213–223
Brosed A. y Guirado D. (2011): Fundamentos de física medica vol.8: Radiobiología y
principios de Oncología. ADI, isbn: 978-84-944186-3-1.
Brosed A. (2012): Fundamentos de física medica vol.1: Medida de la radiación. ADI,
isbn: 978-84-938016-1-8.
Capote R, Jeraj R, Ma C.-M., Rogers D.W.O., Sánchez-Doblado F, Sempau J., Seuntjens
J. and Siebers J. V. (2006): Phase-Space Database for External Beam Radiotherapy:
Summary Report of a Consultants‟ Meeting Technical Report INDC(NDS)-0484,
International Nuclear Data Committee, INDC, IAEA, Vienna, Austria.
Pinza C. y Lliso F. (2009): Control de calidad en aceleradores de electrones para uso
médico. ADI, isbn: 978-84-692-3138-8.
Cherry, S., Sorenson, J., and Phelps, M. (2003). Physics in nuclear medicine. Saunders.
 REFERENCIA BIBLIOGRÁFICA 78

Cunningham jhon Robert, J. H. (1983). The Physics of Radiology. U.S.A: CHARLES

THOMAS- PUBLISHER.
Fano, U. (1956): Differential Inelastic Scattering of Relativistic Charged Particles. En:
Phys. Rev. 102, p. 385–387
Followill, R. Tello V. Hanson W. (1998): An empirical relationship for determining
photon beam quality in TG-21 from a ratio of percent depth doses. 25, p. 1202– 205
Greening, J. (2005): Fundamentals of Radiation Dosimetry. New York, United States:
Taylor and Francis group.
Gunderson, L. L. and Tepper, J. E. (2021). Clinical Radiation Oncology. ELSEVIER. 5
To Edicion.
Krieger T. (2005): Sauer O.A., Monte Carlo- versus pencil-beam-/collapsed-cone-dose
calculation in a heterogeneous multilayer phantom. Phys. Med. Biol. 50 859-868.
Khan, F. (2003): The Physics of Radiation Therapy. Minneapolis, Minnesota: Lippincott
Williams and Wilkins.
King B W, L. (2011): X-ray coherent scattering form factors of tissues, water and plastics
using energy dispersion. Physics in Medicine and Biology, 4377-4397.
Henry Finlay Godson, M. Ravikumar, K.M. Ganesh, S. Sathiyan, Y. Retna Ponmalar,
(2016): Small field output factors: Comparison of measurements with various
detectors and effects of detector orientation with primary jaw setting”, Radiation
Measurements 85, 99-110.
Hubbell. (1982). Photon Mass Attenuation and Energy-absorption Coefficients from
1KeV to 20MeV. Int. J. Appl. Radiat. Isot., 1269 - 1290.
Hussein, E. M. A. (2007). Radiation Mechanics Principles and Practice. ELSEVIER.
Iaea. (2017):Dosimetry of small static fields used in external beam radiotherapy: An
international code of practice for reference and relative dose determination.,
Tecnical Report Series No. 483.
INEN (2021): https://portal.inen.sld.pe/indicadores-anuales-de-gestion-produccion-
hospitalaria.
Nachbar, M., Monnich, D., Boeke, S., Gani, C., Weidner, N., Heinrich, V., lo Russo, M.,
Winter, L. L. L. J., Tsitsekidis, S., Dohm, O., Thorwarth, D., Zips, D., and De-
Colle, C. (2020). Partial breast irradiation with the 1.5T MR-Linac: First patient
 REFERENCIA BIBLIOGRÁFICA 79

treatment and analysis of electron return and stream effects. Radiotherapy and
Oncology, 145:30–35.
Núñez, G. (2019.): Determinación de fluencia de electrones de 6 MeV en terapia
superficial. Córdoba, Argentina: s.n.
Ma CM, Mok E, Kapur A, Pawlicki T, Findley D, et al. (1999): Clinical implementation
of a Monte Carlo treatment planning system. Medical Physics 26 2133-43
Martínez, I. (2008): Diseño y optimización del aplicador de un acelerador tipo microtron
para radioterapia intraoperatoria, Universidad de Barcelona, Tesis de Grado.
Medina, A. (2005): Influencia de la contaminación de electrones en la dosimetría de haces
de fotones de uso clínico. Granada, España, Universidad de Granada, Tesis de
Doctorado.
Oshima T, Aoyama Y, Shimozato T, Sawaki M, Imai T, Ito Y, et al. (2009): An
experimental attenuation plate to improve the dose distribution in intraoperative
electron beam radiotherapy for breast cancer. Phys Med Biol; 54:3491
Pacheco, M. (2009): Efectos de la dispersión múltiple en la simulación realista de
aceleradores lineales de uso clínico. Navarra, España, Universidad de Navarra,
Tesis de Doctorado.
Salvat F., J. Sempau y J. M. Fernández Varea. (2011): PENELOPE, a code system for
Montecarlo simulation of electron and photon transport. España: Facultad de física
(ECM and ICC), Universidad de Barcelona.
Sheikh-Bagheri D. y Rogers D. W. (2002): A Monte Carlo calculation of nine
megavoltage photon beam spectra using the BEAM code Med. Phys. 29 391-402
Jones, D. (1994). “ICRU Report 50, Prescribing, Recording and Reporting Photon Beam
Therapy”. Med. Phys., 21: 833-834.
Juste, B. (2011): Aplicacion del Metodo de Monte Carlo a la Planificacion en
Radioterapia y a la Reconstruccion de Espectros de Fotones de Aceleradores
Lineales de Particulas (LinAc). Valencia, España, Universidad Politécnica de
Valencia, Tesis de Doctorado.
Sharpe, M. B. (2006). “IAEA Technical Reports Series No. 430: Commissioning And
Quality Assurance Of Computerized Planning Systems For Radiation Treatment Of
Cancer”. Med. Phys., 33: 561-561.
 REFERENCIA BIBLIOGRÁFICA 80

TRS 398, (2006): IAEA, Absorbed dose determination in external beam radiotherapy: An
international code of practice for dosimetry based on standards of absorbed dose to
water, Tecnical Report Series No. 398.
Valente, M. (2020). Fundamentos de física médica. Facultad de Matemática, Astronomía,
Física y Computación, UNC.
Wright H. A., Turner J. E., (1985): Hamm R., A Monte Carlo Primer For Health
Physicists.
APÉNDICE 81

APÉNDICES
 APÉNDICE 82

Apéndice A

Se tomó el valor de 50 Gy como dosis referencial ya que la dosis para un cáncer localizada puede
llegar a 70Gy. A continuación, se toma el Quantec que son valores de las tolerancias de los distintos
órganos de riesgo en este caso se toma solo para la próstata (Gunderson, L. L. and Tepper, J. E. ,2021).
 APÉNDICE 83

APENDICE B
Valor del espectro de energía encontrada en la simulación del cabezal.

Elow 1.75 3.850E-04 3.60 1.222E-04

(MeV) MV 1.80 3.734E-04 3.65 1.180E-04
0.00 9.195E-07 1.85 3.591E-04 3.70 1.153E-04
0.05 1.590E-05 1.90 3.506E-04 3.75 1.118E-04
0.10 8.232E-05 1.95 3.377E-04 3.80 1.072E-04
0.15 2.950E-04 2.00 3.249E-04 3.85 1.042E-04
0.20 5.962E-04 2.05 3.158E-04 3.90 1.008E-04
0.25 8.113E-04 2.10 3.058E-04 3.95 9.805E-05
0.30 9.220E-04 2.15 2.974E-04 4.00 9.493E-05
0.35 9.639E-04 2.20 2.867E-04 4.05 9.166E-05
0.40 9.731E-04 2.25 2.776E-04 4.10 8.842E-05
0.45 9.533E-04 2.30 2.686E-04 4.15 8.535E-05
0.50 9.814E-04 2.35 2.592E-04 4.20 8.318E-05
0.55 8.996E-04 2.40 2.519E-04 4.25 7.986E-05
0.60 8.702E-04 2.45 2.438E-04 4.30 7.744E-05
0.65 8.383E-04 2.50 2.371E-04 4.35 7.325E-05
0.70 8.084E-04 2.55 2.286E-04 4.40 7.100E-05
0.75 7.799E-04 2.60 2.233E-04 4.45 6.847E-05
0.80 7.553E-04 2.65 2.154E-04 4.50 6.620E-05
0.85 7.270E-04 2.70 2.086E-04 4.55 6.267E-05
0.90 6.972E-04 2.75 2.037E-04 4.60 5.988E-05
0.95 6.727E-04 2.80 1.972E-04 4.65 5.706E-05
1.00 6.513E-04 2.85 1.919E-04 4.70 5.494E-05
1.05 6.297E-04 2.90 1.845E-04 4.75 5.085E-05
1.10 6.053E-04 2.95 1.797E-04 4.80 4.896E-05
1.15 5.843E-04 3.00 1.743E-04 4.85 4.627E-05
1.20 5.630E-04 3.05 1.699E-04 4.90 4.306E-05
1.25 5.452E-04 3.10 1.646E-04 4.95 4.035E-05
1.30 5.254E-04 3.15 1.599E-04 5.00 3.728E-05
1.35 5.106E-04 3.20 1.561E-04 5.05 3.462E-05
1.40 4.901E-04 3.25 1.492E-04 5.10 3.116E-05
1.45 4.735E-04 3.30 1.461E-04 5.15 2.776E-05
1.50 4.597E-04 3.35 1.414E-04 5.20 2.452E-05
1.55 4.457E-04 3.40 1.377E-04 5.25 2.087E-05
1.60 4.283E-04 3.45 1.329E-04 5.30 1.606E-05
1.65 4.129E-04 3.50 1.292E-04 5.35 9.496E-06
1.70 4.001E-04 3.55 1.255E-04 5.40 3.426E-06
 APÉNDICE 84

APENDICE C

Algoritmo de entrada de haz de energía de cada campo de tratamiento.

ARCHIVO entrada.in en Código PENELOPE

SPOSIT 0 0 100.90 [Coordinates of the source]

FIELD 0 0 0.270 [Conical beam; angles in deg]
.
>>>>>>>> Material data and simulation parameters.
MFNAME prostate.mat [Material file, up to 20 chars]
MSIMPA 1e5 1e5 1e5 0.05 0.05 5e4 5e3 [EABS(1:4),C1,C2,WCC,WCR]
MFNAME CF.mat [Material file, up to 20 chars]
MSIMPA 1e5 1e5 1e5 0.05 0.05 5e4 5e3 [EABS(1:4),C1,C2,WCC,WCR]
MFNAME water.mat [Material file, up to 20 chars]
MSIMPA 1e5 1e5 1e5 0.05 0.05 5e4 5e3 [EABS(1:4),C1,C2,WCC,WCR]
MFNAME aire.mat [Material file, up to 20 chars]
MSIMPA 1e5 1e5 1e5 0.05 0.05 5e4 5e3 [EABS(1:4),C1,C2,WCC,WCR]
.
SPOSIT 0 0 100.90 [Coordinates of the source]
FIELD 0 0 0.0 [Conical beam; angles in deg]
.
>>>>>>>> Material data and simulation parameters.
MFNAME prostate.mat [Material file, up to 20 chars]
MSIMPA 1e5 1e5 1e5 0.05 0.05 5e4 5e3 [EABS(1:4),C1,C2,WCC,WCR]
MFNAME vegiga.mat [Material file, up to 20 chars]
MSIMPA 1e5 1e5 1e5 0.05 0.05 5e4 5e3 [EABS(1:4),C1,C2,WCC,WCR]
MFNAME water.mat [Material file, up to 20 chars]
MSIMPA 1e5 1e5 1e5 0.05 0.05 5e4 5e3 [EABS(1:4),C1,C2,WCC,WCR]
MFNAME aire.mat [Material file, up to 20 chars]
MSIMPA 1e5 1e5 1e5 0.05 0.05 5e4 5e3 [EABS(1:4),C1,C2,WCC,WCR]

SPOSIT 0 0 100.90 [Coordinates of the source]

FIELD 0 0 0.90 [Conical beam; angles in deg]
.
>>>>>>>> Material data and simulation parameters.
MFNAME prostate.mat [Material file, up to 20 chars]
MSIMPA 1e5 1e5 1e5 0.05 0.05 5e4 5e3 [EABS(1:4),C1,C2,WCC,WCR]
MFNAME CF.mat [Material file, up to 20 chars]
MSIMPA 1e5 1e5 1e5 0.05 0.05 5e4 5e3 [EABS(1:4),C1,C2,WCC,WCR]
MFNAME water.mat [Material file, up to 20 chars]
MSIMPA 1e5 1e5 1e5 0.05 0.05 5e4 5e3 [EABS(1:4),C1,C2,WCC,WCR]
MFNAME aire.mat [Material file, up to 20 chars]
MSIMPA 1e5 1e5 1e5 0.05 0.05 5e4 5e3 [EABS(1:4),C1,C2,WCC,WCR]

SPOSIT 0 0 100.90 [Coordinates of the source]

FIELD 0 0 0.180 [Conical beam; angles in deg]
.
>>>>>>>> Material data and simulation parameters.
MFNAME prostate.mat [Material file, up to 20 chars]
MSIMPA 1e5 1e5 1e5 0.05 0.05 5e4 5e3 [EABS(1:4),C1,C2,WCC,WCR]
MFNAME rectum.mat [Material file, up to 20 chars]
MSIMPA 1e5 1e5 1e5 0.05 0.05 5e4 5e3 [EABS(1:4),C1,C2,WCC,WCR]
MFNAME water.mat [Material file, up to 20 chars]
MSIMPA 1e5 1e5 1e5 0.05 0.05 5e4 5e3 [EABS(1:4),C1,C2,WCC,WCR]
MFNAME aire.mat [Material file, up to 20 chars]
 APÉNDICE 85

MSIMPA 1e5 1e5 1e5 0.05 0.05 5e4 5e3 [EABS(1:4),C1,C2,WCC,WCR]

>>>>>>>> Geometry definition file.

GEOMFN Pelvis.geo [Geometry definition file, 20 chars]
.
>>>>>>>> Impact detectors (up to 20 caracteres).
IMPDET 1e5 25e6 100 0 1 [E-window, no. of bins, IPSF, IDCUT] &
IDSPC impdet-prostate.dat [Spectrum file name, 20 chars] &
IDBODY 1 [Active body; one line for each body] &*
IDKPAR 2 [Kind of detected particles, one line each] &*
IMPDET 1e5 25e6 100 0 1 [E-window, no. of bins, IPSF, IDCUT] &
IDSPC impdet-water.dat [Spectrum file name, 20 chars] &
IDBODY 2 [Active body; one line for each body] &*
IDKPAR 2 [Kind of detected particles, one line each] &*
IMPDET 1e5 25e6 100 0 1 [E-window, no. of bins, IPSF, IDCUT] &
IDSPC impdet-CF.dat [Spectrum file name, 20 chars] &
IDBODY 3 [Active body; one line for each body] &*
IDKPAR 2 [Kind of detected particles, one line each] &*
IMPDET 1e5 25e6 100 0 1 [E-window, no. of bins, IPSF, IDCUT] &
IDSPC impdet-rectum.dat [Spectrum file name, 20 chars] &
IDBODY 3 [Active body; one line for each body] &*
IDKPAR 2 [Kind of detected particles, one line each] &*
IMPDET 1e5 25e6 100 0 1 [E-window, no. of bins, IPSF, IDCUT] &
IDSPC impdet-vejiga.dat [Spectrum file name, 20 chars] &
IDBODY 3 [Active body; one line for each body] &*
IDKPAR 2 [Kind of detected particles, one line each] &*
IMPDET 1e5 25e6 100 0 1 [E-window, no. of bins, IPSF, IDCUT] &
IDSPC impdet-skin.dat [Spectrum file name, 20 chars] &
IDBODY 4 [Active body; one line for each body] &*
IDKPAR 2 [Kind of detected particles, one line each] &*
.
>>>>>>>> Distribution of dose
GRIDX -5 5 [X coordinates of the enclosure vertices]
GRIDY -5 5 [Y coordinates of the enclosure vertices]
GRIDZ 0 16.0 [Z coordinates of the enclosure vertices]

GRIDBN 45 45 101 1
.
>>>>>>>> Job properties
RESUME dump2.dmp [Resume from this dump file, 20 chars]
DUMPTO dump2.dmp [Generate this dump file, 20 chars]
DUMPP 180 [Dumping period, in sec]
.
NSIMSH 3e8 [Desired number of simulated showers]
TIME 1e8 [Allotted simulation time, in sec]
.
END [Ends the reading of input data]
 APÉNDICE 86

APENDICE D

Archivo de la creación del cabezal del acelerador lineal Clinac 2100C, paso a paso.

ARCHIVO acel.geo en Código PENELOPE

ACELERADOR LINEAL CLINAC 2100C- 6 MV

0000000000000000000000000000000000000000000000000000000000000000
SURFACE (1) Cilindro del Blanco 1 W
INDICES=(1, 1, 0, 0,-1)
X-SCALE=(9.000000000000000E-01, 0)
Y-SCALE=(9.000000000000000E-01, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE (2) par de planos para el cilindro del blanco de W
INDICES=(0, 0, 1, 0,-1)
Z-SCALE=(4.445000000000000E-02, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE (1) Cilindro Blanco 1.
MATERIAL (1)
SURFACE (1), SIDE POINTER=(-1)
SURFACE (2), SIDE POINTER=(-1)
1111111111111111111111111111111111111111111111111111111111111111
Z-SHIFT=(4.445000000000000E-02, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 3) par de planos para el cilindro del blanco 2 de Cu
INDICES=( 0, 0, 1, 0,-1)
Z-SCALE=( 7.875000000000000E-02, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE ( 2) Cilindro Blanco 2 Cu.
MATERIAL( 2)
SURFACE ( 1), SIDE POINTER=(-1)
SURFACE ( 3), SIDE POINTER=(-1)
1111111111111111111111111111111111111111111111111111111111111111
Z-SHIFT=( 1.676500000000000E-01, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 4) plano para el cilindro de Aire de los blancos
INDICES=( 0, 0, 0, 1, 0)
Z-SHIFT=(-2.000000000000000E+00, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 5) Cicindro de aire de los blancos
INDICES=( 1, 1, 0, 0,-1)
X-SCALE=( 3.000000000000000E-01, 0)
Y-SCALE=( 3.000000000000000E-01, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 6) Plano Z=0. para el cono1 del colimador primario
INDICES=( 0, 0, 0, 1, 0)
Z-SHIFT=( 0.000000000000000E+00, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE ( 3) Cilindro Aire.
MATERIAL( 0)
 APÉNDICE 87

SURFACE ( 4), SIDE POINTER=( 1)

SURFACE ( 5), SIDE POINTER=(-1)
SURFACE ( 6), SIDE POINTER=(-1)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 7) sphere diameter=1.1 cm
INDICES=( 1, 1, 1, 0,-1)
X-SCALE=( 1.500000000000000E+00, 0)
Y-SCALE=( 1.500000000000000E+00, 0)
Z-SCALE=( 1.500000000000000E+00, 0)
Z-SHIFT=( 3.700000000000000E-01, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE ( 4) Esfera Aire.
MATERIAL( 0)
SURFACE ( 7), SIDE POINTER=(-1)
MODULE ( 1)
MODULE ( 2)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 8) Plano para el cono1 del colimador primario
INDICES=( 0, 0, 0, 1, 0)
Z-SHIFT=( 5.500000000000000E+00, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 9) cone 1
INDICES=( 1, 1,-1, 0, 0)
Z-SCALE=( 4.049904764000000E+00, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE (5) cono1 del colimador primario de Aire
MATERIAL ( 0)
SURFACE ( 8), SIDE POINTER=(-1)
SURFACE ( 9), SIDE POINTER=(-1)
SURFACE ( 6), SIDE POINTER=( 1)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 10) Cicindro 2 exterior. Colimador primario
INDICES=( 1, 1, 0, 0,-1)
X-SCALE=( 5.000000000000000E+00, 0)
Y-SCALE=( 5.000000000000000E+00, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE ( 6) Cilindro W.
MATERIAL( 1)
SURFACE ( 8), SIDE POINTER=(-1)
SURFACE ( 10), SIDE POINTER=(-1)
SURFACE ( 4), SIDE POINTER=( 1)
MODULE ( 3)
MODULE ( 5)
MODULE ( 4)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 11) Plano para el cono2 del colimador primario
INDICES=( 0, 0, 0, 1, 0)
Z-SHIFT=( 8.464300956000000E+00, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE ( 7) cono2 del colimador primario de Aire
 APÉNDICE 88

MATERIAL( 0)
SURFACE ( 11), SIDE POINTER=(-1)
SURFACE ( 9), SIDE POINTER=(-1)
SURFACE ( 8), SIDE POINTER=( 1)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 12) Cilindro, colimador primario
INDICES=( 1, 1, 0, 0,-1)
X-SCALE=( 3.800000000000000E+00, 0)
Y-SCALE=( 3.800000000000000E+00, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE ( 8) Cilindro salida del colimador primario
MATERIAL ( 1)
SURFACE ( 8), SIDE POINTER=( 1)
SURFACE ( 12), SIDE POINTER=(-1)
SURFACE ( 11), SIDE POINTER=(-1)
MODULE ( 7)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 13) Ventana de Vacío
INDICES=( 1, 1, 0, 0,-1)
X-SCALE=( 3.900000000000000E+00, 0)
Y-SCALE=( 3.900000000000000E+00, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 14) par de planos para la ventana de vacío
INDICES=( 0, 0, 1, 0,-1)
Z-SCALE=( 8.000000000000000E-01, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE ( 9) Ventana de Vacío.
MATERIAL( 0)
SURFACE ( 13), SIDE POINTER=(-1)
SURFACE ( 14), SIDE POINTER=(-1)
1111111111111111111111111111111111111111111111111111111111111111
Z-SHIFT=( 9.264300956000000E+00, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 15) Cilindro, impact detector
INDICES=( 1, 1, 0, 0,-1)
X-SCALE=( 1.500000000000000E+00, 0)
Y-SCALE=( 1.500000000000000E+00, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 16) Plane Z=-100
INDICES=( 0, 0, 1, 0,-1)
Z-SCALE=( 5.000000000000000E-02, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE ( 10) Detector de impacto 1
MATERIAL( 3)
SURFACE ( 15), SIDE POINTER=(-1)
SURFACE ( 16), SIDE POINTER=(-1)
1111111111111111111111111111111111111111111111111111111111111111
Z-SHIFT=( 1.262300000000000E+01, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 17) Plano inferior para el cono1 del filtro aplanador
 APÉNDICE 89

INDICES=( 0, 0, 0, 1, 0)
Z-SHIFT=( 0.000000000000000E+00, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 18) Plano superior para el cono1 del filtro aplanador
INDICES=( 0, 0, 0, 1, 0)
Z-SHIFT=( 5.300000000000000E-02, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 19) cono 1 del filtro aplanador
INDICES=( 1, 1,-1, 0, 0)
Z-SCALE=( 5.300000000000000E-01, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE ( 11) Filtro Aplanador cono 1. Cu
MATERIAL( 2)
SURFACE ( 17), SIDE POINTER=( 1)
SURFACE ( 19), SIDE POINTER=(-1)
SURFACE ( 18), SIDE POINTER=(-1)
1111111111111111111111111111111111111111111111111111111111111111
Z-SHIFT=( 1.170000000000000E+01, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 20) Plano superior para el cono2 del filtro aplanador
INDICES=( 0, 0, 0, 1, 0)
Z-SHIFT=( 2.083333330000000E-01, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 21) Plano inferior para el cono2 del filtro aplanador
INDICES=( 0, 0, 0, 1, 0)
Z-SHIFT=( 8.333333300000000E-02, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 22) cono 2 del filtro aplanador
INDICES=( 1, 1,-1, 0, 0)
Z-SCALE=( 8.333333330000000E-01, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE ( 12) Filtro Aplanador cono 2 . Cu
MATERIAL( 2)
SURFACE ( 20), SIDE POINTER=(-1)
SURFACE ( 22), SIDE POINTER=(-1)
SURFACE ( 21), SIDE POINTER=( 1)
1111111111111111111111111111111111111111111111111111111111111111
Z-SHIFT=( 1.166966667000000E+01, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 23) Plano superior para el cono3 del filtro aplanador
INDICES=( 0, 0, 0, 1, 0)
Z-SHIFT=( 7.661176470000000E-01, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 24) Plano inferior para el cono3 del filtro aplanador
INDICES=( 0, 0, 0, 1, 0)
Z-SHIFT=( 1.741176470000000E-01, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 25) cone 3 del filtro aplanador
INDICES=( 1, 1,-1, 0, 0)
Z-SCALE=( 6.964705880000000E-01, 0)
 APÉNDICE 90

0000000000000000000000000000000000000000000000000000000000000000
MODULE ( 13) Filtro aplanador cono 3. Cu
MATERIAL( 2)
SURFACE ( 23), SIDE POINTER=(-1)
SURFACE ( 25), SIDE POINTER=(-1)
SURFACE ( 24), SIDE POINTER=( 1)
1111111111111111111111111111111111111111111111111111111111111111
Z-SHIFT=( 1.170388235000000E+01, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 26) Cilindro base del filtro aplanador Cu
INDICES=( 1, 1, 0, 0,-1)
X-SCALE=( 1.500000000000000E+00, 0)
Y-SCALE=( 1.500000000000000E+00, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 27) par de plano en Z para el cilindro
INDICES=( 0, 0, 1, 0,-1)
Z-SCALE=( 5.150000000000000E-02, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE ( 14) cilindro
MATERIAL( 2)
SURFACE ( 26), SIDE POINTER=(-1)
SURFACE ( 27), SIDE POINTER=(-1)
1111111111111111111111111111111111111111111111111111111111111111
Z-SHIFT=( 1.252150000000000E+01, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 28) Cilindro cámara de ionización
INDICES=( 1, 1, 0, 0,-1)
X-SCALE=( 4.762500000000000E+00, 0)
Y-SCALE=( 4.762500000000000E+00, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 29) par de planos cámara de ionización
INDICES=( 0, 0, 1, 0,-1)
Z-SCALE=( 6.350000000000000E-03, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 30) par de planos cámara de ionización
INDICES=( 0, 0, 1, 0,-1)
Z-SCALE=( 6.350000000000000E-03, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 31) par de planos cámara de ionización
INDICES=( 0, 0, 1, 0,-1)
Z-SCALE=( 6.350000000000000E-03, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE ( 15) cilindro 1 cámara
MATERIAL( 4)
SURFACE ( 28), SIDE POINTER=(-1)
SURFACE ( 29), SIDE POINTER=(-1)
1111111111111111111111111111111111111111111111111111111111111111
Z-SHIFT=( 1.483500000000000E+01, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE ( 16) cilindro 1 cámara
 APÉNDICE 91

MATERIAL (4)
SURFACE (28), SIDE POINTER=(-1)
SURFACE (30), SIDE POINTER=(-1)
1111111111111111111111111111111111111111111111111111111111111111
Z-SHIFT=(1.555000000000000E+01, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE (17) cilindro 1 cámara
MATERIAL( 4)
SURFACE ( 28), SIDE POINTER=(-1)
SURFACE ( 31), SIDE POINTER=(-1)
1111111111111111111111111111111111111111111111111111111111111111
Z-SHIFT=( 1.626500000000000E+01, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 32) Mirror
INDICES=( 1, 1, 0, 0,-1)
X-SCALE=( 8.890000000000001E+00, 0)
Y-SCALE=( 8.890000000000001E+00, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 33) par de planos Mirror
INDICES=( 0, 0, 1, 0,-1)
Z-SCALE=( 1.000000000000000E-03, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE ( 18) cilindro Mirror
MATERIAL( 5)
SURFACE ( 32), SIDE POINTER=(-1)
SURFACE ( 33), SIDE POINTER=(-1)
1111111111111111111111111111111111111111111111111111111111111111
THETA=( 3.500000000000000E+01, 0) DEG
PHI=(-9.000000000000000E+01, 0) DEG
Z-SHIFT=( 2.199900000000000E+01, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 34) planos, Y colimador secundario.
INDICES=( 0, 0, 1, 0,-1)
Z-SCALE=( 3.900000000000000E+00, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 35) right-left planes, Y upper right jaw.
INDICES=( 0, 1, 0, 0,-1)
Y-SCALE=( 5.000000000000000E+00, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 36) front-back planes, Y upper right jaw.
INDICES=( 1, 0, 0, 0,-1)
X-SCALE=( 1.000000000000000E+01, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE ( 19) Y RIGHT UPPER JAW.
MATERIAL( 1)
SURFACE ( 34), SIDE POINTER=(-1)
SURFACE ( 35), SIDE POINTER=(-1)
SURFACE ( 36), SIDE POINTER=(-1)
1111111111111111111111111111111111111111111111111111111111111111
OMEGA=( 9.000000000000000E+01, 0) DEG
 APÉNDICE 92

THETA=( 2.862405226000000E+00, 0) DEG

PHI=(-9.000000000000000E+01, 0) DEG
Y-SHIFT=(-6.593253275000000E+00, 0)
Z-SHIFT=( 3.186506549000000E+01, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 37) planos,Y jaw.
INDICES=( 0, 0, 1, 0,-1)
Z-SCALE=( 3.900000000000000E+00, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 38) right-left planes, Y upper right jaw.
INDICES=( 0, 1, 0, 0,-1)
Y-SCALE=( 5.000000000000000E+00, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE ( 20) Y RIGHT UPPER JAW.
MATERIAL( 1)
SURFACE ( 37), SIDE POINTER=(-1)
SURFACE ( 38), SIDE POINTER=(-1)
SURFACE ( 36), SIDE POINTER=(-1)
1111111111111111111111111111111111111111111111111111111111111111
OMEGA=( 9.000000000000000E+01, 0) DEG
THETA=( 2.862405226000000E+00, 0) DEG
PHI=( 9.000000000000000E+01, 0) DEG
Y-SHIFT=( 6.593253275000000E+00, 0)
Z-SHIFT=( 3.186506549000000E+01, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 39) planos,X jaw.
INDICES=( 0, 0, 1, 0,-1)
Z-SCALE=( 3.900000000000000E+00, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 40) right-left planes, X upper right jaw.
INDICES=( 0, 1, 0, 0,-1)
Y-SCALE=( 1.000000000000000E+01, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 41) front-back planes, X upper right jaw.
INDICES=( 1, 0, 0, 0,-1)
X-SCALE=( 5.000000000000000E+00, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE ( 21) X RIGHT UPPER JAW.
MATERIAL( 1)
SURFACE ( 39), SIDE POINTER=(-1)
SURFACE ( 40), SIDE POINTER=(-1)
SURFACE ( 41), SIDE POINTER=(-1)
1111111111111111111111111111111111111111111111111111111111111111
THETA=(-2.862405226000000E+00, 0) DEG
X-SHIFT=(-7.027710542000000E+00, 0)
Z-SHIFT=( 4.055421084000000E+01, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 42) planos,X jaw.
INDICES=( 0, 0, 1, 0,-1)
Z-SCALE=( 3.900000000000000E+00, 0)
 APÉNDICE 93

0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 43) front-back planes, X upper right jaw.
INDICES=( 1, 0, 0, 0,-1)
X-SCALE=( 5.000000000000000E+00, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE ( 22) X RIGHT UPPER JAW.
MATERIAL( 1)
SURFACE ( 42), SIDE POINTER=(-1)
SURFACE ( 40), SIDE POINTER=(-1)
SURFACE ( 43), SIDE POINTER=(-1)
1111111111111111111111111111111111111111111111111111111111111111
THETA=( 2.862405226000000E+00, 0) DEG
X-SHIFT=( 7.027710542000000E+00, 0)
Z-SHIFT=( 4.055421084000000E+01, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 44) Pair of planes, x=+-20
INDICES=( 1, 0, 0, 0,-1)
X-SCALE=( 2.000000000000000E+01, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 45) Pair of planes, y=+-20
INDICES=( 0, 1, 0, 0,-1)
Y-SCALE=( 2.000000000000000E+01, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 46) Plane Z=-100
INDICES=( 0, 0, 1, 0,-1)
Z-SCALE=( 2.000000000000000E+01, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE ( 23) Water phantom
MATERIAL( 6)
SURFACE ( 46), SIDE POINTER=(-1)
SURFACE ( 44), SIDE POINTER=(-1)
SURFACE ( 45), SIDE POINTER=(-1)
1111111111111111111111111111111111111111111111111111111111111111
Z-SHIFT=( 1.200000000000000E+02, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( 47) Plano superior para el Modulo de Aire
INDICES=( 0, 0, 0, 1, 0)
Z-SHIFT=( 1.000000000000000E+02, 0)
0000000000000000000000000000000000000000000000000000000000000000
MODULE ( 24) Acelerador en Aire
MATERIAL( 3)
SURFACE ( 4), SIDE POINTER=( 1)
SURFACE ( 47), SIDE POINTER=(-1)
SURFACE ( 44), SIDE POINTER=(-1)
SURFACE ( 45), SIDE POINTER=(-1)
MODULE (6)
MODULE (10)
MODULE (8)
MODULE (11)
MODULE (12)
 APÉNDICE 94

MODULE (13)
MODULE (14)
MODULE (19)
MODULE (20)
MODULE (21)
MODULE (22)
MODULE (15)
MODULE (16)
MODULE (17)
MODULE (18)
MODULE (9)
0000000000000000000000000000000000000000000000000000000000000000
END 0000000000000000000000000000000000000000000000000000000