Proton radiotherapy is a cancer treatment which has the potential to offer greater cure rates and/or fewer serious side effects than conventional radiotherapy. Its availability in the UK is currently limited to a single low-energy fixed beamline for the treatment of ocular tumours, but a number of facilities designed to treat deep-seated tumours are in development. This thesis focusses on the quantitative use of x-ray computed tomography (CT) images in planning proton radiotherapy treatments. It arrives at several recommendations that can be used to inform clinical protocols for the acquisition of planning scans, and their subsequent use in treatment planning systems. The primary tool developed is a software CT scanner, which simulates images of an anthropomorphic virtual phantom, informed by measurements taken on a clinical scanner. The software is used to investigate the accuracy of the stoichiometric method for calibrating CT image pixel values to proton stopping power, with particular attention paid to the impact of beam hardening and photon starvation artefacts. The strength of the method adopted is in allowing comparison between CT-estimated and exactly-calculated proton stopping powers derived from the same physical data (specified in the phantom), leading to results that are difficult to obtain otherwise. A number of variations of the stoichiometric method are examined, identifying the best-performing calibration phantom and CT tube voltage (kVp). Improvements in accuracy are observed when using a second-pass beam hardening correction algorithm. Also presented is a method for identifying the proton paths where stopping power uncertainties are likely to be greatest. Estimates of the proton range uncertainties caused by CT artefacts and calibration errors are obtained for a range of realistic clinical scenarios. The current practice of including planning margins equivalent to 3.5% of the range is found to ensure coverage in all but the very worst of cases. Results herein suggest margins could be reduced to <2% if the best-performing protocol is followed; however, an analysis specific to the CT scanner and treatment site in question should be carried out before such a change is made in the clinic.

Au cours du traitement d'une tumeur avec des ions légers, la position du pic de Bragg doit être connue avec précision. Une fraction importante des ions incidents subissent des collisions nucléaires avec les noyaux cibles conduisant à l'émission de particules promptes qui peuvent être porteuses d'informations sur le parcours des ions. Ce travail, qui se concentre sur les gamma prompts, montre que le rendement en profondeur de ces émissions est fortement corrélé au parcours des ions et que les taux de comptage mesurés permettent d'envisager un système d'imagerie réaliste, fonctionnant en temps réel. Dans ce but, nous avons réalisé des expériences au GANIL et au GSI avec un détecteur collimaté placé perpendiculairement à l'axe du faisceau et la technique du temps de vol a été utilisée pour réduire le bruit de fond induit par les neutrons et les particules chargées. Des simulations Geant4 ont été réalisées pour concevoir le dispositif expérimental et interpréter les données. Un accord qualitatif entre les simulations et l'expérience est observé pour la quantité d'énergie déposée dans le détecteur et pour la forme du spectre de temps de vol. Cependant, des divergences apparaissent pour le rendement de gamma prompts et la distribution en profondeur des gamma détectés. Ces divergences sont discutées, principalement en termes de modèles de physique nucléaire qui doivent être améliorés. Après avoir sélectionné les modèles physiques offrant les simulations les plus en accord avec les mesures, des études concernant les lieux d'émissions des gamma prompts et l'influence de la diffusion dans la cible ont été réalisés afin de déterminer l'impact sur la corrélation avec le parcours des ions / During the treatment of a tumor with light ions, the Bragg peak location must be accurately known. A significant fraction of the incident ions undergo nuclear collisions with the target nuclei leading to the prompt emission of particles which may carry information on the ion path. This work, which focuses on prompt gamma, shows that the depth profile of these emissions is highly correlated to the ions path and the measured counting rates allow to consider a realistic imaging system, operating in real time. For that purpose, we performed experiments at GANIL and at GSI with a collimated detector placed perpendicular to the beam axis and the time of flight technique was used in order to reduce the noise induced by neutrons and charged particles. Geant4 simulations were performed for the experimental design and data interpretation. A qualitative agreement between simulations and experiment is observed for the amount of energy deposited in the detector and the shape of the time of flight spectrum. However, discrepancies appear for the prompt gamma yield and the depth distribution of gamma detected. These discrepancies are discussed, mainly in terms of nuclear physics models that must be improved. After selecting the physical models which lead to the best agreement between simulations and measurements, studies on the location of prompt gamma emission and on the influence of diffusion in the target were performed to determine the impact on the correlation with the ion path

Ion beam therapy (IBT) is a promising treatment option in radiotherapy. The characteristic physical and biological properties of light ion beams allow for the delivery of highly tumour conformal dose distributions. Related to the sparing of surrounding healthy tissue and nearby organs at risk, it is feasible to escalate the dose in the tumour volume to reach higher tumour control and survival rates. Remarkable clinical outcome was achieved with IBT for radio-resistant, deep-seated, static and well fixated tumour entities. Presumably, more patients could benefit from the advantages of IBT if it would be available for more frequent tumour sites. Those located in the thorax and upper abdominal region are commonly subjected to intra-fractional, respiration related motion. Different motion compensated dose delivery techniques have been developed for active field shaping with scanned pencil beams and are at least available under experimental conditions at the GSI Helmholtzzentrum für Schwerionenforschung (GSI) in Darmstadt, Germany. High standards for quality assurance are required in IBT to ensure a safe and precise dose application. Both underdosage in the tumour and overdosage in the normal tissue might endanger the treatment success. Since minor unexpected anatomical changes e.g. related to patient mispositioning, tumour shrinkage or tissue swelling could already lead to remarkable deviations between planned and delivered dose distribution, a valuable dose monitoring system is desired for IBT. So far, positron emission tomography (PET) is the only in vivo, in situ and non-invasive qualitative dose monitoring method applied under clinical conditions. It makes use of the tissue autoactivation by nuclear fragmentation reactions occurring along the beam path. Among others, β+-emitting nuclides are generated and decay according to their half-life under the emission of a positron. The subsequent positron-electron annihilation creates two 511 keV photons which are emitted in opposite direction and can be detected as coincidence event by a dedicated PET scanner. The induced three-dimensional (3D) β+-activity distribution in the patient can be reconstructed from the measured coincidences. Conclusions about the delivered dose distribution can be drawn indirectly from a comparison between two β+-activity distributions: the measured one and an expected one generated by a Monte-Carlo simulation. This workflow has been proven to be valuable for the dose monitoring in IBT when it was applied for about 440 patients, mainly suffering from deep-seated head and neck tumours that have been treated with 12C ions at GSI. In the presence of intra-fractional target motion, the conventional 3D PET data processing will result in an inaccurate representation of the β+-activity distribution in the patient. Four-dimensional, time-resolved (4D) reconstruction algorithms adapted to the special geometry of in-beam PET scanners allow to compensate for the motion related blurring artefacts. Within this thesis, a 4D maximum likelihood expectation maximization (MLEM) reconstruction algorithm has been implemented for the double-head scanner Bastei installed at GSI. The proper functionality of the algorithm and its superior performance in terms of suppressing motion related blurring artefacts compared to an already applied co-registration approach has been demonstrated by a comparative simulation study and by dedicated measurements with moving radioactive sources and irradiated targets. Dedicated phantoms mainly made up of polymethyl methacrylate (PMMA) and a motion table for regular one-dimensional (1D) motion patterns have been designed and manufactured for the experiments. Furthermore, the general applicability of the 4D MLEM algorithm for more complex motion patterns has been demonstrated by the successful reduction of motion artefacts from a measurement with rotating (two-dimensional moving) radioactive sources. For 1D cos^2 and cos^4 motion, it has been clearly illustrated by systematic point source measurements that the motion influence can be better compensated with the same number of motion phases if amplitude-sorted instead of time-sorted phases are utilized. In any case, with an appropriate parameter selection to obtain a mean residual motion per phase of about half of the size of a PET crystal size, acceptable results have been achieved. Additionally, it has been validated that the 4D MLEM algorithm allows to reliably access the relevant parameters (particle range and lateral field position and gradients) for a dose verification in intra-fractionally moving targets even from the intrinsically low counting statistics of IBT-PET data. To evaluate the measured β+-activity distribution, it should be compared to a simulated one that is expected from the moving target irradiation. Thus, a 4D version of the simulation software is required. It has to emulate the generation of β+-emitters under consideration of the intra-fractional motion, their decay at motion state dependent coordinates and to create listmode data streams from the simulated coincidences. Such a revised and extended version that has been compiled for the special geometry of the Bastei PET scanner is presented within this thesis. The therapy control system provides information about the exact progress of the motion compensated dose delivery. This information and the intra-fractional target motion needs to be taken into account for simulating realistic β+-activity distributions. A dedicated preclinical phantom simulation study has been performed to demonstrate the correct functionality of the 4D simulation program and the necessity of the additional, motion-related input parameters. Different to the data evaluation for static targets, additional effort is required to avoid a potential misleading interpretation of the 4D measured and simulated β+-activity distribu- tions in the presence of deficient motion mitigation or data processing. It is presented that in the presence of treatment errors the results from the simulation might be in accordance to the measurement although the planned and delivered dose distribution are different. In contrast to that, deviations may occur between both distributions which are not related to anatomical changes but to deficient 4D data processing. Recommendations are given in this thesis to optimize the 4D IBT-PET workflow and to prevent the observer from a mis-interpretation of the dose monitoring data. In summary, the thesis contributes on a large scale to a potential future application of the IBT-PET monitoring for intra-fractionally moving target volumes by providing the required reconstruction and simulation algorithms. Systematic examinations with more realistic, multi-directional and irregular motion patterns are required for further improvements. For a final rating of the expectable benefit from a 4D IBT-PET dose monitoring, future investigations should include real treatment plans, breathing curves and 4D patient CT images. / Die Ionenstrahltherapie (englisch: ion beam therapy, IBT) ist eine vielversprechende Behandlungsoption im Bereich der Strahlentherapie. Die charakteristischen physikalischen und biologischen Eigenschaften der Ionenstrahlen werden genutzt, um tumorkonformale Dosisverteilungen zu erzeugen. Die verbesserte Schonung des an den Tumor angrenzenden Normalgewebes und eventuell naheliegender Risikoorgane ermöglicht eine Dosissteigerung im Zielgebiet und somit potentiell höhere Tumorkontroll- und Überlebensraten. Für tiefliegende, gegenüber konventioneller Strahlung resistente, statische und gut fixierte Tumore wurden bereits beachtliche klinische Resultate erzielt. Wahrscheinlich könnten noch mehr Patienten von den Vorteilen der IBT profitieren, wenn diese auch für häufiger auftretende und intrafraktionell bewegliche Tumore uneingeschränkt nutzbar wäre. Verschiedene bewegungskompensierte Bestrahlungsmethoden wurden entwickelt und stehen zumindest unter experimentellen Bedingungen für weitere Untersuchungen am GSI Helmholtzzentrum für Schwerionenforschung (GSI) in Darmstadt zur Verfügung. Um eine sichere und präzise Dosisapplikation in der IBT zu ermöglichen, werden hohe Anforderungen an die Qualitätssicherung gesetzt. Sowohl auftretende Überdosierungen im Normalgewebe als auch Unterdosierungen im Tumor können den Therapieerfolg gefährden. Da bereits kleine, unerwartete anatomische Veränderungen, zum Beispiel durch Fehlpositionierung des Patienten, Schrumpfung des Tumors oder Schwellungen, zu erheblichen Abweichungen zwischen geplanter und applizierter Dosisverteilung führen können, gibt es Bestrebungen, die applizierte Dosis zumindest qualitativ zu verifizieren. Die Positronen-Emissions-Tomografie (PET) ist derzeit die einzige, bereits klinisch erprobte Methode für ein in vivo, in situ und nicht-invasives qualitatives Dosismonitoring. Diese Methode ist im Stande, die Autoaktivierung des bestrahlten Gewebes zu erfassen, welche aufgrund von Kernfragmentierungsprozessen entlang des Strahlweges erzeugt wird. Unter anderem werden in diesen Reaktionen instabile Nuklide erzeugt, die entsprechend ihrer Halbwertszeit unter Emission eines Positrons zerfallen. Bei der anschließenden Positron-Elektron-Annihilation werden zwei 511keV Photonen in entgegengesetzter Richtung emittiert und können mittels eines geeigneten PET-Scanners als Koinzidenzereignis detektiert werden. Die im Patienten induzierte dreidimensionale (3D) β+-Aktivitätsverteilung kann aus den gemessenen Koinzidenzen rekonstruiert werden. Ein Vergleich der gemessenen mit einer erwarteten, mittels Monte-Carlo Simulation erzeugten β+-Aktivitätsverteilung erlaubt es, Schlussfolgerungen über die tatsächlich im Patienten deponierte 3D Dosisverteilung zu ziehen. Diese Art der Datenauswertung wurde erfolgreich für die qualitative Dosisverifikation von über 440 Patienten eingesetzt, deren Tumore (vorwiegend im Kopf- und Halsbereich) an der GSI mit 12C-Ionen bestrahlt wurden. Bei der konventionellen 3D IBT-PET-Datenverarbeitung wird eine mögliche intrafraktionelle Bewegung des Zielgebietes nicht berücksichtigt und fehlerhaft rekonstruierte β+-Aktivitätsverteilungen sind die Folge. Daher werden vierdimensionale, zeitaufgelöste (4D) Rekonstruktionsalgorithmen benötigt, die für die spezielle Geometrie eines in-beam PET-Scanner adaptiert wurden und eine Kompensation der bewegungsinduzierten Artefakte ermöglichen. Im Rahmen der vorliegenden Arbeit wurde für den an der GSI installierten Doppelkopf-PET-Scanner Bastei ein 4D Maximum-Likelihood-Expectation-Maximization (MLEM) Algorithmus implementiert. Die Funktionsfähigkeit des Algorithmus sowie dessen verbesserte Reduktion von Bewegungsartefakten im Vergleich zu einem bereits vorhandenen Koregistrierungsansatz wurde anhand verschiedener Messungen mit bewegten radioaktiven Quellen und bestrahlten Phantomen sowie einer vergleichenden Simulationsstudie dargelegt. Für die Experimente wurden entsprechende Phantomgeometrien (zumeist aus Polymethylmethacrylat (PMMA)) sowie ein Bewegungstisch für reguläre eindimensionale (1D) Bewegungsmuster entworfen und gefertigt. Zudem wurde durch die erfolgreiche, quasi-statische und nahezu artefaktfreie Rekonstruktion einer rotierenden und sich damit zweidimensional bewegenden Aktivitätsverteilung die prinzipielle Anwendbarkeit des 4D MLEM Algorithmus für komplexere Bewegungsmuster gezeigt. Systematische Punktquellenmessungen mit 1D cos^2- und cos^4-förmigen Bewegungsmustern haben deutlich gemacht, dass der Bewegungseinfluss mit der gleichen Anzahl an Bewegungsphasen besser kompensiert werden kann, wenn die Bewegungsphasen entsprechend der Bewegungsamplitude anstelle der -phase unterteilt sind. In jedem Fall können aber zufriedenstellende Rekonstruktionsergebnisse erzielt werden, wenn durch geeignete Parameterwahl eine mittlere Restbewegung pro Bewegungsphase von maximal etwa der halben Größe eines Detektorkristalls eingestellt wird. Durch weitere Experimente konnte gezeigt werden, dass nach der Rekonstruktion mit dem 4D MLEM Algorithmus die relevanten Parameter für die qualitative Dosisverifikation (Teilchenreichweite, laterale Feldposition und -gradienten) zuverlässig erfasst werden können. Dies ist auch dann der Fall, wenn nur eine verminderte Anzahl an Koinzidenzereignissen, so wie sie unter klinischen Bedingungen zu erwarten ist, für die Auswertung verwendet wird. Um die gemessene β+-Aktivitätsverteilung besser zu beurteilen, sollte sie mit einer simulierten, für die bewegungskompensierte Bestrahlung erwarteten Verteilung verglichen werden und es bedarf deshalb einer 4D Version der Simulationssoftware. Diese muss die Erzeugung sowie den Zerfall der Positronenemitter unter Berücksichtigung der intrafraktionellen Bewegung simulieren und aus den gültigen Koinzidenzereignissen Listmode-Datensätze erstellen. Eine derart überarbeitet Version des Simulationsprogramms wurde für den Bastei PET-Scanner erstellt und wird in dieser Arbeit vorgestellt. Informationen über den exakten Verlauf der bewegungskompensierten Bestrahlung werden durch das Therapiekontrollsystem geliefert. Diese Informationen sowie die intrafraktionelle Bewegung werden in die Simulation realistischer β+-Aktivitätsverteilungen bzw. der zugehörigen Listmode-Datensätze einbezogen. Anhand einer präklinischen Phantom-Simulationsstudie wurde die korrekte Funktionsweise des Simulationsprogramms sowie die Notwendigkeit der zusätzlichen Parameter gezeigt. Im Gegensatz zur Datenauswertung für statische Zielvolumina bedarf es bei intrafraktioneller Bewegung gegebenenfalls zusätzlichen Aufwand, um eine Fehlinterpretation aus dem Vergleich der gemessenen und simulierten β+-Aktivitätsverteilung zu vermeiden. In der vorliegenden Arbeit wird beispielhaft gezeigt, dass sich bei fehlerhafter Bewegungskompensation die gemessene und simulierte β+-Aktivitätsverteilung einander ähneln können, obwohl die applizierte Dosisverteilung deutlich von der geplanten abweicht. Im Gegensatz dazu können auch Abweichungen zwischen Messung und Simulation auftreten, die nicht auf anatomische Veränderungen, sondern auf eine ungenaue 4D Datenverarbeitung zurückzuführen sind. Es werden Vorschläge unterbreitet, um den Prozess der 4D IBT-PET Datenauswertung zu optimieren und somit Fehlinterpretationen zu vermeiden. Die vorliegende Dissertationsschrift enthält durch die Bereitstellung der benötigten 4D Rekonstruktions- und Simulationsprogramme grundlegende Arbeiten für eine mögliche zukünftige Anwendung der 4D IBT-PET als qualitatives Dosismonitoring bei intrafraktionell bewegten Zielvolumina. Für weitere Verbesserungen des Verfahrens sind zusätzliche systematische Betrachtungen mit realistischeren, mehrdimensionalen und unregelmäßigen Bewegungsmustern notwendig. Zukünftige Untersuchungen sollten außerdem echte Bestrahlungspläne, Atemkurven sowie 4D Patienten-CT-Daten einschließen, um den erwartbaren Nutzen eines 4D IBT-PET Dosismonitorings besser abschätzen zu können.

Ion beam therapy (IBT) is a promising treatment option in radiotherapy. The characteristic physical and biological properties of light ion beams allow for the delivery of highly tumour conformal dose distributions. Related to the sparing of surrounding healthy tissue and nearby organs at risk, it is feasible to escalate the dose in the tumour volume to reach higher tumour control and survival rates. Remarkable clinical outcome was achieved with IBT for radio-resistant, deep-seated, static and well fixated tumour entities. Presumably, more patients could benefit from the advantages of IBT if it would be available for more frequent tumour sites. Those located in the thorax and upper abdominal region are commonly subjected to intra-fractional, respiration related motion. Different motion compensated dose delivery techniques have been developed for active field shaping with scanned pencil beams and are at least available under experimental conditions at the GSI Helmholtzzentrum für Schwerionenforschung (GSI) in Darmstadt, Germany. High standards for quality assurance are required in IBT to ensure a safe and precise dose application. Both underdosage in the tumour and overdosage in the normal tissue might endanger the treatment success. Since minor unexpected anatomical changes e.g. related to patient mispositioning, tumour shrinkage or tissue swelling could already lead to remarkable deviations between planned and delivered dose distribution, a valuable dose monitoring system is desired for IBT. So far, positron emission tomography (PET) is the only in vivo, in situ and non-invasive qualitative dose monitoring method applied under clinical conditions. It makes use of the tissue autoactivation by nuclear fragmentation reactions occurring along the beam path. Among others, β+-emitting nuclides are generated and decay according to their half-life under the emission of a positron. The subsequent positron-electron annihilation creates two 511 keV photons which are emitted in opposite direction and can be detected as coincidence event by a dedicated PET scanner. The induced three-dimensional (3D) β+-activity distribution in the patient can be reconstructed from the measured coincidences. Conclusions about the delivered dose distribution can be drawn indirectly from a comparison between two β+-activity distributions: the measured one and an expected one generated by a Monte-Carlo simulation. This workflow has been proven to be valuable for the dose monitoring in IBT when it was applied for about 440 patients, mainly suffering from deep-seated head and neck tumours that have been treated with 12C ions at GSI. In the presence of intra-fractional target motion, the conventional 3D PET data processing will result in an inaccurate representation of the β+-activity distribution in the patient. Four-dimensional, time-resolved (4D) reconstruction algorithms adapted to the special geometry of in-beam PET scanners allow to compensate for the motion related blurring artefacts. Within this thesis, a 4D maximum likelihood expectation maximization (MLEM) reconstruction algorithm has been implemented for the double-head scanner Bastei installed at GSI. The proper functionality of the algorithm and its superior performance in terms of suppressing motion related blurring artefacts compared to an already applied co-registration approach has been demonstrated by a comparative simulation study and by dedicated measurements with moving radioactive sources and irradiated targets. Dedicated phantoms mainly made up of polymethyl methacrylate (PMMA) and a motion table for regular one-dimensional (1D) motion patterns have been designed and manufactured for the experiments. Furthermore, the general applicability of the 4D MLEM algorithm for more complex motion patterns has been demonstrated by the successful reduction of motion artefacts from a measurement with rotating (two-dimensional moving) radioactive sources. For 1D cos^2 and cos^4 motion, it has been clearly illustrated by systematic point source measurements that the motion influence can be better compensated with the same number of motion phases if amplitude-sorted instead of time-sorted phases are utilized. In any case, with an appropriate parameter selection to obtain a mean residual motion per phase of about half of the size of a PET crystal size, acceptable results have been achieved. Additionally, it has been validated that the 4D MLEM algorithm allows to reliably access the relevant parameters (particle range and lateral field position and gradients) for a dose verification in intra-fractionally moving targets even from the intrinsically low counting statistics of IBT-PET data. To evaluate the measured β+-activity distribution, it should be compared to a simulated one that is expected from the moving target irradiation. Thus, a 4D version of the simulation software is required. It has to emulate the generation of β+-emitters under consideration of the intra-fractional motion, their decay at motion state dependent coordinates and to create listmode data streams from the simulated coincidences. Such a revised and extended version that has been compiled for the special geometry of the Bastei PET scanner is presented within this thesis. The therapy control system provides information about the exact progress of the motion compensated dose delivery. This information and the intra-fractional target motion needs to be taken into account for simulating realistic β+-activity distributions. A dedicated preclinical phantom simulation study has been performed to demonstrate the correct functionality of the 4D simulation program and the necessity of the additional, motion-related input parameters. Different to the data evaluation for static targets, additional effort is required to avoid a potential misleading interpretation of the 4D measured and simulated β+-activity distribu- tions in the presence of deficient motion mitigation or data processing. It is presented that in the presence of treatment errors the results from the simulation might be in accordance to the measurement although the planned and delivered dose distribution are different. In contrast to that, deviations may occur between both distributions which are not related to anatomical changes but to deficient 4D data processing. Recommendations are given in this thesis to optimize the 4D IBT-PET workflow and to prevent the observer from a mis-interpretation of the dose monitoring data. In summary, the thesis contributes on a large scale to a potential future application of the IBT-PET monitoring for intra-fractionally moving target volumes by providing the required reconstruction and simulation algorithms. Systematic examinations with more realistic, multi-directional and irregular motion patterns are required for further improvements. For a final rating of the expectable benefit from a 4D IBT-PET dose monitoring, future investigations should include real treatment plans, breathing curves and 4D patient CT images.:1 Motivation 1.1 Potential and obstacles of ion beam therapy 1.2 Objectives of the thesis 2 Ion beam therapy and moving targets 2.1 Physical and biological properties of ion beams 2.1.1 Dose deposition 2.1.2 Biological effectivity 2.2 Technical aspects of ion beam delivery 2.2.1 Active and passive beam delivery technique 2.2.2 Beam monitoring for pencil beam scanning 2.2.3 Considerations in treatment planning related to patient CT image 2.3 Organ motion in ion beam therapy 2.3.1 Types of organ motion 2.3.2 Detection of intra-fractional motion 2.3.3 Motion compensated ion beam therapy 2.4 Dose monitoring by means of positron emission tomography 2.4.1 Principle of PET imaging in ion beam therapy 2.4.2 In-beam PET at GSI 3 Reconstruction of in-beam PET data taken from moving targets 3.1 Reconstruction algorithm 3.1.1 3D MLEM reconstruction applied at GSI 3.1.2 4D in-beam PET reconstruction methods 3.1.3 Comparison of gated co-registration and 4D MLEM 3.2 Experiments with moving radioactive sources 3.2.1 Rotation of radioactive sources 3.2.2 One-dimensional point source motion 3.3 In-beam PET measurements with moving targets 3.3.1 Verification of lateral field position and gradients 3.3.2 Verification of particle range 3.4 Summary and discussion 4 Simulation of phase-sorted in-beam PET data for moving targets 4.1 Upgrading the IBT-PET simulation from 3D to 4D 4.1.1 General and motion-related simulation demands 4.1.2 Input parameters for the 4D simulation program 4.1.3 Workflow of the 4D simulation program 4.2 Verification of the 4D simulation code by means of a preclinical phantom study 4.2.1 Experiment design 4.2.2 4D in-beam PET data simulation 4.2.3 Comparison with 3D simulation 4.3 Summary and discussion 5 Interpretation of 4D IBT-PET data with respect to deficient motion mitigation or data processing 5.1 Detectability of failed motion mitigation 5.1.1 Failure in gated beam delivery 5.1.2 Failure in lateral target tracking 5.2 Deficient correlation between motion and PET data 5.3 Recommendations for the 4D IBT-PET workflow 6 Summary and outlook 7 Appendix A Transformation matrices A.1 Composition of transformation matrices A.2 Storage of transformation matrices A.3 Transformation matrices for rotation B Noise reduction in analogue signals by FFT-based filtering C Motion tables and corresponding motion patterns C.1 Rotational motion C.2 Motion table with stepping motor for precise 1D motion patterns C.3 Motion table enabling relative target movement D Synchronisation of PET, motion and beam monitoring data E Sorting PET data by time or amplitude and calculating corresponding mean offsets Bibliography / Die Ionenstrahltherapie (englisch: ion beam therapy, IBT) ist eine vielversprechende Behandlungsoption im Bereich der Strahlentherapie. Die charakteristischen physikalischen und biologischen Eigenschaften der Ionenstrahlen werden genutzt, um tumorkonformale Dosisverteilungen zu erzeugen. Die verbesserte Schonung des an den Tumor angrenzenden Normalgewebes und eventuell naheliegender Risikoorgane ermöglicht eine Dosissteigerung im Zielgebiet und somit potentiell höhere Tumorkontroll- und Überlebensraten. Für tiefliegende, gegenüber konventioneller Strahlung resistente, statische und gut fixierte Tumore wurden bereits beachtliche klinische Resultate erzielt. Wahrscheinlich könnten noch mehr Patienten von den Vorteilen der IBT profitieren, wenn diese auch für häufiger auftretende und intrafraktionell bewegliche Tumore uneingeschränkt nutzbar wäre. Verschiedene bewegungskompensierte Bestrahlungsmethoden wurden entwickelt und stehen zumindest unter experimentellen Bedingungen für weitere Untersuchungen am GSI Helmholtzzentrum für Schwerionenforschung (GSI) in Darmstadt zur Verfügung. Um eine sichere und präzise Dosisapplikation in der IBT zu ermöglichen, werden hohe Anforderungen an die Qualitätssicherung gesetzt. Sowohl auftretende Überdosierungen im Normalgewebe als auch Unterdosierungen im Tumor können den Therapieerfolg gefährden. Da bereits kleine, unerwartete anatomische Veränderungen, zum Beispiel durch Fehlpositionierung des Patienten, Schrumpfung des Tumors oder Schwellungen, zu erheblichen Abweichungen zwischen geplanter und applizierter Dosisverteilung führen können, gibt es Bestrebungen, die applizierte Dosis zumindest qualitativ zu verifizieren. Die Positronen-Emissions-Tomografie (PET) ist derzeit die einzige, bereits klinisch erprobte Methode für ein in vivo, in situ und nicht-invasives qualitatives Dosismonitoring. Diese Methode ist im Stande, die Autoaktivierung des bestrahlten Gewebes zu erfassen, welche aufgrund von Kernfragmentierungsprozessen entlang des Strahlweges erzeugt wird. Unter anderem werden in diesen Reaktionen instabile Nuklide erzeugt, die entsprechend ihrer Halbwertszeit unter Emission eines Positrons zerfallen. Bei der anschließenden Positron-Elektron-Annihilation werden zwei 511keV Photonen in entgegengesetzter Richtung emittiert und können mittels eines geeigneten PET-Scanners als Koinzidenzereignis detektiert werden. Die im Patienten induzierte dreidimensionale (3D) β+-Aktivitätsverteilung kann aus den gemessenen Koinzidenzen rekonstruiert werden. Ein Vergleich der gemessenen mit einer erwarteten, mittels Monte-Carlo Simulation erzeugten β+-Aktivitätsverteilung erlaubt es, Schlussfolgerungen über die tatsächlich im Patienten deponierte 3D Dosisverteilung zu ziehen. Diese Art der Datenauswertung wurde erfolgreich für die qualitative Dosisverifikation von über 440 Patienten eingesetzt, deren Tumore (vorwiegend im Kopf- und Halsbereich) an der GSI mit 12C-Ionen bestrahlt wurden. Bei der konventionellen 3D IBT-PET-Datenverarbeitung wird eine mögliche intrafraktionelle Bewegung des Zielgebietes nicht berücksichtigt und fehlerhaft rekonstruierte β+-Aktivitätsverteilungen sind die Folge. Daher werden vierdimensionale, zeitaufgelöste (4D) Rekonstruktionsalgorithmen benötigt, die für die spezielle Geometrie eines in-beam PET-Scanner adaptiert wurden und eine Kompensation der bewegungsinduzierten Artefakte ermöglichen. Im Rahmen der vorliegenden Arbeit wurde für den an der GSI installierten Doppelkopf-PET-Scanner Bastei ein 4D Maximum-Likelihood-Expectation-Maximization (MLEM) Algorithmus implementiert. Die Funktionsfähigkeit des Algorithmus sowie dessen verbesserte Reduktion von Bewegungsartefakten im Vergleich zu einem bereits vorhandenen Koregistrierungsansatz wurde anhand verschiedener Messungen mit bewegten radioaktiven Quellen und bestrahlten Phantomen sowie einer vergleichenden Simulationsstudie dargelegt. Für die Experimente wurden entsprechende Phantomgeometrien (zumeist aus Polymethylmethacrylat (PMMA)) sowie ein Bewegungstisch für reguläre eindimensionale (1D) Bewegungsmuster entworfen und gefertigt. Zudem wurde durch die erfolgreiche, quasi-statische und nahezu artefaktfreie Rekonstruktion einer rotierenden und sich damit zweidimensional bewegenden Aktivitätsverteilung die prinzipielle Anwendbarkeit des 4D MLEM Algorithmus für komplexere Bewegungsmuster gezeigt. Systematische Punktquellenmessungen mit 1D cos^2- und cos^4-förmigen Bewegungsmustern haben deutlich gemacht, dass der Bewegungseinfluss mit der gleichen Anzahl an Bewegungsphasen besser kompensiert werden kann, wenn die Bewegungsphasen entsprechend der Bewegungsamplitude anstelle der -phase unterteilt sind. In jedem Fall können aber zufriedenstellende Rekonstruktionsergebnisse erzielt werden, wenn durch geeignete Parameterwahl eine mittlere Restbewegung pro Bewegungsphase von maximal etwa der halben Größe eines Detektorkristalls eingestellt wird. Durch weitere Experimente konnte gezeigt werden, dass nach der Rekonstruktion mit dem 4D MLEM Algorithmus die relevanten Parameter für die qualitative Dosisverifikation (Teilchenreichweite, laterale Feldposition und -gradienten) zuverlässig erfasst werden können. Dies ist auch dann der Fall, wenn nur eine verminderte Anzahl an Koinzidenzereignissen, so wie sie unter klinischen Bedingungen zu erwarten ist, für die Auswertung verwendet wird. Um die gemessene β+-Aktivitätsverteilung besser zu beurteilen, sollte sie mit einer simulierten, für die bewegungskompensierte Bestrahlung erwarteten Verteilung verglichen werden und es bedarf deshalb einer 4D Version der Simulationssoftware. Diese muss die Erzeugung sowie den Zerfall der Positronenemitter unter Berücksichtigung der intrafraktionellen Bewegung simulieren und aus den gültigen Koinzidenzereignissen Listmode-Datensätze erstellen. Eine derart überarbeitet Version des Simulationsprogramms wurde für den Bastei PET-Scanner erstellt und wird in dieser Arbeit vorgestellt. Informationen über den exakten Verlauf der bewegungskompensierten Bestrahlung werden durch das Therapiekontrollsystem geliefert. Diese Informationen sowie die intrafraktionelle Bewegung werden in die Simulation realistischer β+-Aktivitätsverteilungen bzw. der zugehörigen Listmode-Datensätze einbezogen. Anhand einer präklinischen Phantom-Simulationsstudie wurde die korrekte Funktionsweise des Simulationsprogramms sowie die Notwendigkeit der zusätzlichen Parameter gezeigt. Im Gegensatz zur Datenauswertung für statische Zielvolumina bedarf es bei intrafraktioneller Bewegung gegebenenfalls zusätzlichen Aufwand, um eine Fehlinterpretation aus dem Vergleich der gemessenen und simulierten β+-Aktivitätsverteilung zu vermeiden. In der vorliegenden Arbeit wird beispielhaft gezeigt, dass sich bei fehlerhafter Bewegungskompensation die gemessene und simulierte β+-Aktivitätsverteilung einander ähneln können, obwohl die applizierte Dosisverteilung deutlich von der geplanten abweicht. Im Gegensatz dazu können auch Abweichungen zwischen Messung und Simulation auftreten, die nicht auf anatomische Veränderungen, sondern auf eine ungenaue 4D Datenverarbeitung zurückzuführen sind. Es werden Vorschläge unterbreitet, um den Prozess der 4D IBT-PET Datenauswertung zu optimieren und somit Fehlinterpretationen zu vermeiden. Die vorliegende Dissertationsschrift enthält durch die Bereitstellung der benötigten 4D Rekonstruktions- und Simulationsprogramme grundlegende Arbeiten für eine mögliche zukünftige Anwendung der 4D IBT-PET als qualitatives Dosismonitoring bei intrafraktionell bewegten Zielvolumina. Für weitere Verbesserungen des Verfahrens sind zusätzliche systematische Betrachtungen mit realistischeren, mehrdimensionalen und unregelmäßigen Bewegungsmustern notwendig. Zukünftige Untersuchungen sollten außerdem echte Bestrahlungspläne, Atemkurven sowie 4D Patienten-CT-Daten einschließen, um den erwartbaren Nutzen eines 4D IBT-PET Dosismonitorings besser abschätzen zu können.:1 Motivation 1.1 Potential and obstacles of ion beam therapy 1.2 Objectives of the thesis 2 Ion beam therapy and moving targets 2.1 Physical and biological properties of ion beams 2.1.1 Dose deposition 2.1.2 Biological effectivity 2.2 Technical aspects of ion beam delivery 2.2.1 Active and passive beam delivery technique 2.2.2 Beam monitoring for pencil beam scanning 2.2.3 Considerations in treatment planning related to patient CT image 2.3 Organ motion in ion beam therapy 2.3.1 Types of organ motion 2.3.2 Detection of intra-fractional motion 2.3.3 Motion compensated ion beam therapy 2.4 Dose monitoring by means of positron emission tomography 2.4.1 Principle of PET imaging in ion beam therapy 2.4.2 In-beam PET at GSI 3 Reconstruction of in-beam PET data taken from moving targets 3.1 Reconstruction algorithm 3.1.1 3D MLEM reconstruction applied at GSI 3.1.2 4D in-beam PET reconstruction methods 3.1.3 Comparison of gated co-registration and 4D MLEM 3.2 Experiments with moving radioactive sources 3.2.1 Rotation of radioactive sources 3.2.2 One-dimensional point source motion 3.3 In-beam PET measurements with moving targets 3.3.1 Verification of lateral field position and gradients 3.3.2 Verification of particle range 3.4 Summary and discussion 4 Simulation of phase-sorted in-beam PET data for moving targets 4.1 Upgrading the IBT-PET simulation from 3D to 4D 4.1.1 General and motion-related simulation demands 4.1.2 Input parameters for the 4D simulation program 4.1.3 Workflow of the 4D simulation program 4.2 Verification of the 4D simulation code by means of a preclinical phantom study 4.2.1 Experiment design 4.2.2 4D in-beam PET data simulation 4.2.3 Comparison with 3D simulation 4.3 Summary and discussion 5 Interpretation of 4D IBT-PET data with respect to deficient motion mitigation or data processing 5.1 Detectability of failed motion mitigation 5.1.1 Failure in gated beam delivery 5.1.2 Failure in lateral target tracking 5.2 Deficient correlation between motion and PET data 5.3 Recommendations for the 4D IBT-PET workflow 6 Summary and outlook 7 Appendix A Transformation matrices A.1 Composition of transformation matrices A.2 Storage of transformation matrices A.3 Transformation matrices for rotation B Noise reduction in analogue signals by FFT-based filtering C Motion tables and corresponding motion patterns C.1 Rotational motion C.2 Motion table with stepping motor for precise 1D motion patterns C.3 Motion table enabling relative target movement D Synchronisation of PET, motion and beam monitoring data E Sorting PET data by time or amplitude and calculating corresponding mean offsets Bibliography

Cette thèse s’inscrit dans une volonté d’améliorer la compréhension des mécanismes fondamentaux de radiolyse de biomolécules par des ions accélérés, à l’échelle moléculaire. Ainsi, les ions étudiés ont été de différentes nature (H+, He2+, C6+) et de différentes énergies, correspondant à une gamme de densité de dépôt d’énergie allant de 0,3 à 1000 eV/nm.Dans le vivant, l’eau ayant une place prépondérante, la compréhension de la radiolyse de l’eau est essentielle. L’espèce la plus réactive produite en milieu aéré, le radical hydroxyle (HO•) a été quantifiée en utilisant une sonde spécifique, l’acide 3-coumarine-carboxylique.Les dégâts indirects aux biomolécules, via les espèces issues de la radiolyse de l’eau, ont été étudiés en solution aqueuse diluée sur deux systèmes : un acide aminé, la phénylalanine et une protéine, la myoglobine. Les effets directs de radiolyse ont été étudiés sur la myoglobine en gels concentrés hydratés. Les phénomènes de radiolyse ont été caractérisés pour décrire les mécanismes en jeu et les produits issus de la radiolyse de la phénylalanine ont été systématiquement identifiées et quantifiées. / The goal of this thesis is to achieve a better understanding of fundamental mechanisms of the radiolysis of biomolecules by accelerated ions, at the molecular scale. To do so, different type of ions have been used (H+, He2+, C6+) at various energies, corresponding to densities of energy deposition from 0,3 to 1000 eV/nm.The main component in biological systems is water. Therefore, the comprehension of the water radiolysis under ions irradiation is essential. One of the most reactive species produced in aerated conditions, the hydroxyl radical (HO•), has been quantified using a specific probe, the 3- carboxylic acid coumarin.Indirect effects of radiolysis on biomolecules, involving water radiolysis species, have been studied in dilute aqueous solutions on two different systems: phenylalanine, an amino acid, and a protein, myoglobin. Direct radiolysis effect were studied on concentrated hydrogels of myoglobin ad other proteins. Elucidation of radiolysis mechanisms and quantification of phenylalanine radiolysis products were systematically performed.

An increasing number of cancer patients is treated with proton beams or other light ion beams which allow to deliver dose precisely to the tumor. However, the depth dose distribution of these particles, which enables this precision, is sensitive to deviations from the treatment plan, as e.g. anatomical changes. Thus, to assure the quality of the treatment, a non-invasive in-vivo dose verification is highly desired. This monitoring of particle therapy relies on the detection of secondary radiation which is produced by interactions between the beam particles and the nuclei of the patient’s tissue. Up to now, the only clinically applied method for in-vivo dosimetry is Positron Emission Tomography which makes use of the beta+-activity produced during the irradiation (PT-PET). Since from a PT-PET measurement the applied dose cannot be directly deduced, the simulated distribution of beta+-emitting nuclei is used as a basis for the analysis of the measured PT-PET data. Therefore, the reliable modeling of the production rates and the spatial distribution of the beta+-emitters is required. PT-PET applied during instead of after the treatment is referred to as in-beam PET. A challenge concerning in-beam PET is the design of the PET camera, because a standard full-ring scanner is not feasible. For instance, a double-head PET camera is applicable, but low count rates and the limited solid angle coverage can compromise the image quality. For this reason, a detector system which provides a time resolution allowing the incorporation of time-of-flight information (TOF) into the iterative reconstruction algorithm is desired to improve the quality of the reconstructed images. Secondly, Prompt Gamma Imaging (PGI), a technique based on the detection of prompt gamma-rays, is currently pursued. Concerning the emissions of prompt gamma-rays during particle irradiation, experimental data is not sufficiently available, making simulations necessary. Compton cameras are based on the detection of incoherently scattered photons and are investigated with respect to PGI. Monte Carlo simulations serve for the optimization of the camera design and the evaluation of criteria for the selection of measured events. Thus, for in-beam PET and PGI dedicated detection systems and, moreover, profound knowledge about the corresponding radiation fields are required. Using various simulation codes, this thesis contributes to the modelling of the beta+-emitters and photons produced during particle irradiation, as well as to the evaluation and optimization of hardware for both techniques. Concerning the modeling of the production of the relevant beta+-emitters, the abilities of the Monte Carlo simulation code PHITS and of the deterministic, one-dimensional code HIBRAC were assessed. The Monte Carlo tool GEANT4 was applied for an additional comparison. For irradiations with protons, helium, lithium, and carbon, the depth-dependent yields of the simulated beta+-emitters were compared to experimental data. In general, PHITS underestimated the yields of the considered beta+-emitters in contrast to GEANT4 which provided acceptable values. HIBRAC was substantially extended to enable the modeling of the depth-dependent yields of specific nuclides. For proton beams and carbon ion beams HIBRAC can compete with GEANT4 for this application. Since HIBRAC is fast, compact, and easy to modify, it could be a basis for the simulations of the beta+-emitters in clinical application. PHITS was also applied to the modeling of prompt gamma-rays during proton irradiation following an experimental setup. From this study, it can be concluded that PHITS could be an alternative to GEANT4 in this context. Another aim was the optimization of Compton camera prototypes. GEANT4 simulations were carried out with the focus on detection probabilities and the rate of valid events. Based on the results, the feasibility of a Compton camera setup consisting of a CZT detector and an LSO or BGO detector was confirmed. Several recommendations concerning the design and arrangement of the Compton camera prototype were derived. Furthermore, several promising event selection strategies were evaluated. The GEANT4 simulations were validated by comparing simulated to measured energy depositions in the detector layers. This comparison also led to the reconsideration of the efficiency of the prototype. A further study evaluated if electron-positron pairs resulting from pair productions could be detected with the existing prototype in addition to Compton events. Regarding the efficiency and the achievable angular resolution, the successful application of the considered prototype as pair production camera to the monitoring of particle therapy is questionable. Finally, the application of a PET camera consisting of Resistive Plate Chambers (RPCs) providing a good time resolution to in-beam PET was discussed. A scintillator-based PET camera based on a commercially available scanner was used as reference. This evaluation included simulations of the detector response, the image reconstructions using various procedures, and the analysis of image quality. Realistic activity distributions based on real treatment plans for carbon ion therapy were used. The low efficiency of the RPC-based PET camera led to images of poor quality. Neither visually nor with the semi-automatic tool YaPET a reliable detectability of range deviations was possible. The incorporation of TOF into the iterative reconstruction algorithm was especially advantageous for the considered RPC-based PET camera in terms of convergence and artifacts. The application of the real-time capable back projection method Direct TOF for the RPCbased PET camera resulted in an image quality comparable to the one achieved with the iterative algorihms. In total, this study does not indicate the further investigation of RPC-based PET cameras with similar efficiency for in-beam PET application. To sum up, simulation studies were performed aimed at the progress of in-vivo dosimetry. Regarding the modeling of the beta+-emitter production and prompt gamma-ray emissions, different simulation codes were evaluated. HIBRAC could be a basis for clinical PT-PET simulations, however, a detailed validation of the underlying cross section models is required. Several recommendations for the optimization of a Compton Camera prototype resulted from systematic variations of the setup. Nevertheless, the definite evaluation of the feasibility of a Compton camera for PGI can only be performed by further experiments. For PT-PET, the efficiency of the detector system is the crucial factor. Due to the obtained results for the considered RPC-based PET camera, the focus should be kept to scintillator-based PET cameras for this purpose.

L'interaction des rayonnements ionisants avec le vivant est marquée par des phénomènes stochastiques importants aussi bien en termes de dosimétrie physique que des effets biologiques induits. Cette thèse aborde trois problématiques de la thérapie et de l'estimation du risque des radiations pour la santé, à l'aide d'outils de modélisation et de simulations Monte Carlo. En effet, des calculs de l'énergie spécifique dans des volumes de différentes tailles ont montré que les amplitudes des fluctuations dépendent fortement de la taille de la cible. Elles sont particulièrement grandes dans le cas des cibles nanométriques. À partir de ces calculs, une étude sur la taille des dosimètres implantables pour le monitorage des traitements de radiothérapie a montré que des dimensions au moins micrométriques sont nécessaires pour assurer des mesures fiables. Les mêmes calculs ont permis l'analyse des effets de faibles doses d'irradiation, notamment la compatibilité de différentes tailles de cibles avec des données expérimentales d'aberrations chromosomiques. Les résultats suggèrent que l'activation du réseau mitochondrial peut être liée au déclenchement de mécanismes de radiorésistance dans les cellules CAL51. Finalement, un nouveau modèle (NanOx) de prédiction de l'efficacité de l'hadronthérapie (radiothérapie par faisceaux d'ions) est présenté et appliqué à la lignée cellulaire V79. Ce modèle est complètement stochastique et intègre les calculs de dosimétrie à plusieurs échelles pour modéliser des effets locaux et non locaux pouvant correspondre respectivement à des lésions de l'ADN et à un stress oxydatif / The interaction between ionizing radiation and living tissues is characterized by stochastic phenomena with non-negligible consequences both in terms of physical dosimetry and induced biological effects. The present work addresses three issues concerning radiotherapy and the estimation of radiation risks for health, by means of modeling tools and Monte Carlo simulations. Indeed, specific energy calculations in volumes of different sizes showed that the level of fluctuations strongly depends on the target size. Such fluctuations are especially high in the case of nanometric targets. Based on these calculations, a study about the size of implantable dosimeters employed in the monitoring of radiotherapy treatments demonstrated that these dosimeters should have at least micrometric dimensions in order to ensure reliable measurements. The same calculations have allowed the analysis of the effects of low doses of radiation, namely the compatibility between different target sizes and experimental data regarding chromosomal aberrations. The results suggest that the activation of the mitochondrial network may be linked to the triggering of radioresistance mechanisms for the CAL51 cell line. Finally, a new model (NanOx) to predict the effectiveness of particle therapy (radiotherapy with ion beams) is presented and applied to the V79 cell line. Such a model is completely stochastic and integrates the dosimetry calculations at multiple scales for modeling local and non-local effects, which can correspond respectively to DNA lesions and cellular oxidative stress

Hadron Therapy (HT) exploits properties of ion radiation to gain therapeutic advantages over existing photon-based forms of external radiation therapy. However, its relative superiority and cost-effectiveness have not been proven for all clinical situations. Establishing a robust evidence base for the development of best treatment practices is one of the major challenges for the field. This thesis investigates two research infrastructures for building this essential evidence. First, the thesis develops main components of a metadata-driven software architecture for the collection of clinical information and its analysis. This architecture acknowledges the diversity in the domain and supports data interoperability by sharing information models. Their compliance to common metamodels guarantees that primary data and analysis results can be interpreted outside of the immediate production context. This is a fundamental necessity for all aspects of the evidence creation process. A metamodel of data capture forms is developed with unique properties to support data collection and documentation in this architecture. The architecture's potential to support complex analysis processes is demonstrated with the help of a novel metamodel for Markov model based simulations, as used for the synthesis of evidence in health-economic assessments. The application of both metamodels is illustrated on the example of HT. Since the biological effect of particle radiation is a major source of uncertainty in HT, in its second part, this thesis undertakes first investigations towards a new research facility for bio-medical experiments with ion beams. It examines the feasibility of upgrading LEIR, an existing accelerator at the European Organisation for Nuclear Research (CERN), with a new slow extraction and investigates transport of the extracted beam to future experiments. Possible configurations for the slow-resonant extraction process are identified, and designs for horizontal and vertical beam transport lines developed. The results of these studies indicate future research directions towards a new ion beam facility for biomedical research.

An increasing number of cancer patients is treated with proton beams or other light ion beams which allow to deliver dose precisely to the tumor. However, the depth dose distribution of these particles, which enables this precision, is sensitive to deviations from the treatment plan, as e.g. anatomical changes. Thus, to assure the quality of the treatment, a non-invasive in-vivo dose verification is highly desired. This monitoring of particle therapy relies on the detection of secondary radiation which is produced by interactions between the beam particles and the nuclei of the patient’s tissue. Up to now, the only clinically applied method for in-vivo dosimetry is Positron Emission Tomography which makes use of the beta+-activity produced during the irradiation (PT-PET). Since from a PT-PET measurement the applied dose cannot be directly deduced, the simulated distribution of beta+-emitting nuclei is used as a basis for the analysis of the measured PT-PET data. Therefore, the reliable modeling of the production rates and the spatial distribution of the beta+-emitters is required. PT-PET applied during instead of after the treatment is referred to as in-beam PET. A challenge concerning in-beam PET is the design of the PET camera, because a standard full-ring scanner is not feasible. For instance, a double-head PET camera is applicable, but low count rates and the limited solid angle coverage can compromise the image quality. For this reason, a detector system which provides a time resolution allowing the incorporation of time-of-flight information (TOF) into the iterative reconstruction algorithm is desired to improve the quality of the reconstructed images. Secondly, Prompt Gamma Imaging (PGI), a technique based on the detection of prompt gamma-rays, is currently pursued. Concerning the emissions of prompt gamma-rays during particle irradiation, experimental data is not sufficiently available, making simulations necessary. Compton cameras are based on the detection of incoherently scattered photons and are investigated with respect to PGI. Monte Carlo simulations serve for the optimization of the camera design and the evaluation of criteria for the selection of measured events. Thus, for in-beam PET and PGI dedicated detection systems and, moreover, profound knowledge about the corresponding radiation fields are required. Using various simulation codes, this thesis contributes to the modelling of the beta+-emitters and photons produced during particle irradiation, as well as to the evaluation and optimization of hardware for both techniques. Concerning the modeling of the production of the relevant beta+-emitters, the abilities of the Monte Carlo simulation code PHITS and of the deterministic, one-dimensional code HIBRAC were assessed. The Monte Carlo tool GEANT4 was applied for an additional comparison. For irradiations with protons, helium, lithium, and carbon, the depth-dependent yields of the simulated beta+-emitters were compared to experimental data. In general, PHITS underestimated the yields of the considered beta+-emitters in contrast to GEANT4 which provided acceptable values. HIBRAC was substantially extended to enable the modeling of the depth-dependent yields of specific nuclides. For proton beams and carbon ion beams HIBRAC can compete with GEANT4 for this application. Since HIBRAC is fast, compact, and easy to modify, it could be a basis for the simulations of the beta+-emitters in clinical application. PHITS was also applied to the modeling of prompt gamma-rays during proton irradiation following an experimental setup. From this study, it can be concluded that PHITS could be an alternative to GEANT4 in this context. Another aim was the optimization of Compton camera prototypes. GEANT4 simulations were carried out with the focus on detection probabilities and the rate of valid events. Based on the results, the feasibility of a Compton camera setup consisting of a CZT detector and an LSO or BGO detector was confirmed. Several recommendations concerning the design and arrangement of the Compton camera prototype were derived. Furthermore, several promising event selection strategies were evaluated. The GEANT4 simulations were validated by comparing simulated to measured energy depositions in the detector layers. This comparison also led to the reconsideration of the efficiency of the prototype. A further study evaluated if electron-positron pairs resulting from pair productions could be detected with the existing prototype in addition to Compton events. Regarding the efficiency and the achievable angular resolution, the successful application of the considered prototype as pair production camera to the monitoring of particle therapy is questionable. Finally, the application of a PET camera consisting of Resistive Plate Chambers (RPCs) providing a good time resolution to in-beam PET was discussed. A scintillator-based PET camera based on a commercially available scanner was used as reference. This evaluation included simulations of the detector response, the image reconstructions using various procedures, and the analysis of image quality. Realistic activity distributions based on real treatment plans for carbon ion therapy were used. The low efficiency of the RPC-based PET camera led to images of poor quality. Neither visually nor with the semi-automatic tool YaPET a reliable detectability of range deviations was possible. The incorporation of TOF into the iterative reconstruction algorithm was especially advantageous for the considered RPC-based PET camera in terms of convergence and artifacts. The application of the real-time capable back projection method Direct TOF for the RPCbased PET camera resulted in an image quality comparable to the one achieved with the iterative algorihms. In total, this study does not indicate the further investigation of RPC-based PET cameras with similar efficiency for in-beam PET application. To sum up, simulation studies were performed aimed at the progress of in-vivo dosimetry. Regarding the modeling of the beta+-emitter production and prompt gamma-ray emissions, different simulation codes were evaluated. HIBRAC could be a basis for clinical PT-PET simulations, however, a detailed validation of the underlying cross section models is required. Several recommendations for the optimization of a Compton Camera prototype resulted from systematic variations of the setup. Nevertheless, the definite evaluation of the feasibility of a Compton camera for PGI can only be performed by further experiments. For PT-PET, the efficiency of the detector system is the crucial factor. Due to the obtained results for the considered RPC-based PET camera, the focus should be kept to scintillator-based PET cameras for this purpose.