1. Introduction
The ability to perform in vivo dosimetry for external beam radiotherapy (EBRT) is instinctively desirable, but not mandated in many countries. A review article by Ford and Evans on incident learning systems (ILS) for radiation oncology outlined the recent occurrences of rare but fatal radiation incidents [
1]. Ideal in vivo dosimetry systems would identify and quantify a therapeutic radiation incident, and, as a part of an ILS, would be valuable to track treatment outcomes in long term databases, even if errors in delivery were not fatal. A report into the implementation of in vivo dosimetry for EBRT in the UK [
2], as well as two recent review articles [
3,
4], all conclude that electronic portal imaging (EPID) dosimetry is the primary commercially available candidate to measure transit dose, as part of an in vivo dosimetry solution.
EPID dosimetry utilises the flat panel imager, designed to aid in patient positioning, that is available on most modern linacs. There are EPID dosimetry systems that have been used to verify the in vivo dose of patient treatment on a forward [
5] or back-projection method [
6,
7,
8]. The hope is that either through calculating the expected dose at the detector level or back projecting the dose through the planning CT, that errors in treatment will be detected. A direct comparison of the forward and back-projected dose distributions reconstructed by the EPID concluded that both can be useful in the identification of in vivo dosimetry errors [
9]. Through the implementation of an in vivo EPID dosimetry system in the Netherlands, relevant dosimetry errors, predominately due to changes in patient anatomy, were identified for 1 in every 300 patients treated with intensity-modulated radiation therapy (IMRT) or volumetric modulated arc therapy (VMAT) [
10]. Despite this preliminary success, EPID based in vivo dosimetry has yet to have a widespread implementation.
A major fundamental limitation of EPID dosimetry systems is the high Z materials the imager comprises of which induce a non-uniform dose–response as a function of the photon energy spectrum [
11]. To mitigate this effect, large amounts of commissioning data are required to form correction factor maps associated with a variable energy spectrum, caused by primary beam attenuation in the patient and beam scatter contributions from the patient and detector itself [
3]. The lack of vendor information on how to commission EPID dosimetry, along with the need for better software to perform corrections, has been suggested to be partly responsible for the absence of widespread clinical uptake of EPID dosimetry [
3]. Despite this, one of the main motives for recommending an EPID dosimetry system is that it is often already installed on modern linacs [
4], making it appear to be a convenient option. The complex commissioning and data correction process could be significantly reduced if vendors offered a water equivalent EPID, but no such system is commercially available. A major limitation to the production of a water equivalent EPID is that attempting to reduce the use of high Z material results in a reduced sensitivity of the imager [
12].
Two previous studies have proposed dual detector systems [
13,
14] where an additional detector, optimised for dosimetry, would be mechanically attached on top of the EPID. Both identified that accurate dosimetry could be performed in these conditions without the EPID negatively impacting the response of the detector, while the additional detector had minimal impact on the quality of EPID images. The first study adopted an array of ionization chambers [
13], while the second used a pixelated silicon detector [
14]. The use of the pixelated silicon was shown to be capable of performing accurate dosimetry, that could be directly compared with the treatment planning system (TPS) for a plan delivered to a heterogeneous phantom. It was also capable of correctly identifying a 5 mm shift in target positioning and a 5.2% overdose. While this system required no field size or scatter corrections, the detector used did have a high dose rate dependence, a limited sensitive area with a limited sensitivity from each diode and a poor spatial resolution. To make this dual detector system advantageous over the EPID, these issues needed to be addressed. The development of a large area, high sensitivity and spatial resolution, pixelated silicon detector array, with minimal dose rate dependence would allow advancement of transit in vivo dosimetry systems, that avoids the large commissioning work required for EPID dosimetry and allow for co-registration of patient imaging and dose distribution.
2. Materials and Methods
In this work, a new large area pixelated silicon detector, named “Magic Plate-987” (MP987) is presented and used to perform transit dosimetry. To verify the suitability of the MP987 detector and its readout system as a dosimetry tool in EBRT, a comprehensive set of basic characterisation measurements were undertaken to evaluate the performance of the device. Once completed, the MP987 was tested to evaluate its performance in transit dosimetry measurements.
The MP987 detector presented in
Figure 1a, comprises 987 diodes, ultrasonically bonded to a 0.5 mm thick fibreglass printed circuit board (PCB). Assembling was achieved by using the patented “drop in” technology adopted from the “Magic Plate” detector [
15]. The diodes used to populate the detector have a surface area of 1.5 × 1.5 mm
2 and are presented in
Figure 1b. They have a low resistivity p-type substrate with a 40 µm thick n-type epitaxial layer. The n-p junction extends as close as possible to the edges of the die, having a sensitive surface area of 1.27 × 1.27 mm
2. The large surface area of the sensitive volume and shallow epitaxial depth has three intended purposes: increase the signal to noise ratio for low dose rate measurements, increase radiation hardness and significantly reduce the effect of trapped charge at Si-SiO
2 interfaces after exposure to radiation [
16].
The MP987 utilises two different detector pitches to maximise the spatial resolution of the detector. The central area consists of a 13 × 13 diode formation at a 5 mm pitch, with the surrounding area having a pitch of 7.5 mm, giving a total detector area of 22.5 × 21 cm
2, as displayed in
Figure 1a. This allows a balance of high spatial resolution, for measuring small fields at the centre while extending the active area for larger and more conventional field sizes. The dual pitch avoids having to add multiple layers to the PCB, which would be required for routing the pixel connections to the electronic interface if the full surface was at the smaller pitch. Additional layers would increase the PCB thickness, requiring ground planes to avoid cross-talk, which could reduce the water equivalence of the detector and increase the attenuation for the EPID imager. The new PCB results in the MP987 having more than 8 times the number of active pixels and greater than 4.7 times the active area of the original “Magic Plate”.
To handle the large number of active pixels a custom data acquisition system (DAQ) was developed to be capable of reading out up to 1024 channels, on a pulse-by-pulse basis. The system developed was based on the 64-channel analogue front end chip, the AFE0064 (Texas Instrument, Dallas, TX, USA) and the field programable gate array (FPGA) Spartan-6 from Xilinx (USA) [
17]. A custom front-end board was developed for the MP987 which consisted of four AFE0064 chips, making each module capable of handling 256 channels. Four 256 channels modules are connected in a star-like network and readout simultaneously by a single FPGA using a common master bus for communication (
Figure 2). The FPGA was connected to a computer through a USB 2.0 interface. Custom software and firmware were developed to communicate with the FPGA, receive, decode, and display data in real-time.
2.1. Percentage Depth Dose
The percentage depth dose (PDD) measurement, along with all other measurements presented, used a 6 MV flattened beam delivered by a Clinac IX (Varian Medical Systems, Palo Alto, CA, USA) linear accelerator. The PDD measurements were performed in Solid Water (SW) (Promis Electro Optic, Wijchen, The Netherlands) with the source to surface distance (SSD) kept constant at 100 cm. The response of the MP987 central pixel irradiated with a field size of 10 × 10 cm2 at each measured depth (between 1.5 and 20 cm) was normalised to its response at dmax (1.5 cm depth, reference conditions). The same measurements were made with a Farmer (Scanditronix/Wellhofer FC23-C) ionisation chamber (IC), used as a reference detector.
2.2. Dose Per Pulse
A characterisation of the dose per pulse (DPP) dependence is essential when taking measurements at extended SSD. To characterise the MP987 DPP dependence, the response at varying SSD between 80 to 300 cm, with a constant 1.5 cm of SW buildup used. The response of three central pixels at each measurement distance was normalised to their response at 100 cm SSD, creating the vector S(SSD). Only pixels at the centre were considered, to avoid any introduction of energy dependence due to the varying energy spectrum of the flattened beam. The same set of data were acquired by the Farmer chamber and named
Sref (SSD). The dose per pulse dependence is then calculated as the ratio
S(SSD)/Sref (SSD). A more detailed description of the method adopted is outlined by Brace et al. [
14]. Five repetitions of the measurement were taken, and error bars were calculated as one standard deviation.
2.3. Radiation Damage
The epitaxial diode, adopted to manufacture the original MP121, suffered an uncharacteristic response variation with irradiation by Co-60 gamma rays [
15,
16]. The new topology of the diode has been designed to mitigate that large variation and obtain a more stable pixel response. To quantify the response stability of the new diode, five test structures have been irradiated in steps of 10 kGy at the ANSTO’s GATRI facility using a 1.2 kGy/h Cobalt-60 source. An initial measurement of the collected charge was taken at the linac’s reference conditions. The same measurement was taken after 10, 20 and 30 kGy of accumulated dose. The charge collected in each subsequent measurement was normalised to the initial measurement and the mean and standard deviation were calculated over the five samples.
2.4. Uniformity
A measure of the uniformity of diode response was performed by delivering a beam of 30 × 30 cm
2 with the detector placed at 10 cm depth in SW. Three irradiations of 200 MU each were delivered, the cumulative response of each diode after each delivery was averaged. Each pixel’s average response after three measurements was then normalised to the average value of all the pixels. The standard deviation of these normalised values was taken as a measure of the uniformity. The normalised values were used to equalise future measurements, to account for any variation in pixel sensitivity. This is the same procedure adopted by Aldosari et al. [
17]. A separate dark field measurement is not required as the pre-amplifier allows for an automatic subtraction of the baseline sampled before each measurement assuming that temperature is not changing.
2.5. Output Factor
Output factors (OF) for the MP987 were performed in SW and compared to the ionisation chamber CC04 (IBA Dosimetry, Schwarzenbruck, Germany). Field sizes between 3 × 3 cm2 and 20 × 20 cm2 were delivered three times each and averaged. The detectors were placed at 10 cm depth, 100 cm SSD with 10 cm of backscatter material. The response at all field sizes was normalised to the response of the 10 × 10 cm2 field.
2.6. Profiles
Central axis (CAX) profiles for 3 × 3, 5 × 5 and 10 × 10 cm2 fields were acquired and compared with the response of a PFD3G Photon diode (IBA Dosimetry, Schwarzenbruck, Germany). These measurements were taken in the same conditions adopted for the OF.
2.7. Linearity
With the diode at reference conditions a range of deliveries between 5 and 1000 MU were delivered. Each measurement is repeated three times and the average collected charge versus nominal MU was plotted. The linearity was quantified by a linear fit and its R2 value. The linearity can also provide the calibration factor assuming the measurements were taken at reference conditions (6 MV photon beam, field size of 10 × 10 cm2, SSD of 100 cm and 1.5 cm depth in a solid water phantom).
Calibration of the charge collected by the MP987 to the dose used in the proceeding measurements was performed with a different in-phantom measurement. The MP987 was placed at 10 cm depth inside SW at 100 cm SSD, and 100 MU was delivered three times with a 10 × 10 cm2 field. A reference measurement with a CC04 was used to find the ratio between the dose at 10 cm depth to the reference condition. The charge calibration for the MP987 was not done directly at the reference condition because at dmax the MP987 can saturate, unless using a very small integration time. Therefore, this methodology allows for the use of a larger integration that reduces noise in the system that is important for application with EPID.
2.8. Transit Dosimetry Characterisation
To validate the use of the MP987 as a transit dosimeter, a series of square fields from 1 × 1 to 20 × 20 cm
2 were delivered to a 15 cm thick SW phantom.
Figure 2 shows the MP987 placed and secured above the EPID with 7 mm of build-up (a standard 30 × 30 cm
2 solid water slab) above the detector and 5 mm between the detector and the EPID surface. The surface of the MP987 was set to 150 cm SSD and the centre of the SW was placed at isocentre.
The Philips Pinnacle
3 treatment planning system (Version 16.0) was used to calculate the predicted dose–response at the detector level using the same methodology adopted in [
14]. The phantom was CT scanned and the DICOM image was padded with additional voxels to extend the field of view. The full volume of the detector is represented in the extended CT dataset by adding a 12 mm thick water equivalent slab (7 mm of build up plus 5 mm of MP987 PCB and supporting material) at SSD of 150 cm. The dose was then scored at the plane located at 7 mm depth inside the water block and extracted from the treatment plan after a complete simulation of the treatment. As the dose gradient in the buildup region above the detector is large and there is a margin of uncertainty of the real position of the MP987 in respect to the extended CT dataset, an evaluation of the error produced by this uncertainty has been simulated using the TPS. The error bars are calculated as half the difference in dose between the slices at 6 and 8 mm deep.
2.9. IMRT Transit Dosimetry Study
A four-segment, IMRT lung treatment was delivered to a heterogeneous lung phantom with a 2 cm diameter spherical target made of SW. All deliveries were delivered from a gantry angle of 0 degrees (collapsed IMRT). The MP987 was used to measure the transit dose from each of the four segments, in the dual detector configuration. A screenshot of the lung phantom CT scan with the imported detector block is displayed in
Figure 3. A full description of the phantom used can be found in [
14]. A prescription of 36 Gy in 6 fractions was set to the PTV. The dose overlay extracted from the TPS was used for comparison with the MP987, using the same procedure adopted for the homogeneous SW phantom of
Section 2.8. The dose calculated by the TPS for each segment was output separately to allow dose verification to be performed.
The dose–response was compared between the MP987 and the TPS through gamma analysis [
18] using two separate criteria, 3%/3 mm and 5%/5 mm (dose/distance to agreement). To allow gamma analysis between the two data sets, the MP987 data were interpolated to a 1 mm grid size using a radial basis function [
19], to match the dose grid size of the TPS. The quoted distance to agreement used is in reference to the distance of the field, as it would be at the plane of the isocentre.
The four segments were first delivered as planned, with the target at isocentre, then again with a deliberate target misalignment introduced. The target was shifted laterally 5 mm off isocentre for the error treatment, to see if the dosimetry was sensitive enough to detect this target misalignment.
4. Discussion
The MP987 has shown to be suitable for performing a wide variety of phantom dosimetry measurements. The results not only demonstrate the dosimetric characteristics of the detector but the capability of the custom DAQ to perform real-time readout on a pulse-by-pulse basis for close to 1000 channels. Measurements carried out in phantom proved that most of the limitations from the previous generation of diodes adopted in the “Magic Plate” detector, associated with dose per pulse dependence and stability with radiation damage, are negligible in the new MP987 diodes. The reduction in these effects brings the technology closer towards a correction free in vivo dosimetry solution. Additionally, it has been demonstrated that MP987 is uniquely suited to the intended application of transit dosimetry. When a homogeneous phantom is used as the target, measured and predicted dose at the EPID level for field sizes between 3 × 3 and 10 × 10 cm2 agreed within error. While correction factors may be required for larger field sizes, the discrepancy remains within 3% for field sizes up to 20 × 20 cm2. In the first pre-clinical test pass rates for absolute dose gamma analysis of the four segments of a step and shoot IMRT plan showed between 94.2% and 97.4% pass fractions with a 5%/5 mm criteria. While this collapsed IMRT study is not a comprehensive assessment of multiple patient plans on a truly anthropomorphic phantom, it is still considered by the authors to be a significant result that will motive future in vivo dosimetry studies with the MP987.
The spatial resolution of MP987 may be too sparse for beamlets with their penumbra falling in the peripheral area of the detector. The high sensitivity of the diodes and the real-time multichannel electrometer performing pulse by pulse acquisition makes the detector uniquely suited to perform analysis at each control point of the delivery. Future developments of the MP987 include the mechanical setup to secure the system above the EPID without triggering the collision sensor and allowing the rotation of the linac’s gantry and allow delivery in clinical settings.
A potential additional application of the MP987 system is in high dose rate (HDR) brachytherapy source tracking. The compact mechanical setup, the large area and the high sensitivity make it uniquely suitable for this application. Therefore, both a fully integrated dual detector for in vivo dosimetry in external beam radiotherapy and HDR brachytherapy source tracking will be the subject of future work with the MP987.