BEAM ENERGY MEASUREMENT SYSTEM

Abstract

A time-of-fight measurement system for measuring energy of a pulsed hadron beam, wherein each pulse of the beam is structured into a series of bunches of charged particles, said bunches being repeated according to a repetition rate of the order of magnitude of radiofrequency. The system comprises a first detector, a second detector and a third detector arranged along a beam path, each of the detectors being configured to detect the passage of a bunch of charged particles and provide an output signal dependent on phase of the detected bunch, wherein the second detector is spaced apart from the first detector by a first distance and wherein the third detector is spaced apart from the second detector by a second distance, wherein the first distance is set out in such a way as that time of flight of the bunch from the first detector to the second detector is approximately equal to, or lower than a repetition period of the bunches, and wherein the second distance is set out in such a way as that time of flight of the bunch from the second detector to the third detector is greater than a multiple of the repetition period of the bunches, and a processing unit configured to a) calculate phase shifts between the output signals of the detectors, and b) calculate energy of the pulse based on the calculated phase shifts.

Claims

1. A time-of-flight (TOF) measurement system for measuring energy of a pulsed hadron beam, wherein each pulse of the beam is structured into a series of bunches (B) of charged particles, said bunches being repeated according to a repetition rate of the order of magnitude of radiofrequency, said system comprising a first detector (1), a second detector (2) and a third detector (3) arranged along a beam path (10), each of said detectors being configured to detect the passage of a bunch (B) of charged particles and provide an output signal (.sub.νPP.1, νPP.2, νPP.3) dependent on phase of the detected bunch (B), wherein the second detector (2) is spaced apart from the first detector (1) by a first distance (L.sub.12) and wherein the third detector (3) is spaced apart from the second detector (2) by a second distance (L.sub.23), wherein said first distance is set out in such a way as that time of flight (t.sub.12) of the bunch (B) from the first detector (1) to the second detector (2) is approximately equal to, or lower than a repetition period (T.sub.RFQ) of the bunches (B), and wherein said second distance is set out in such a way as that time of flight (T.sub.23) of the bunch (B) from the second detector (2) to the third detector (3) is greater than a multiple of the repetition period (T.sub.RFQ) of the bunches (B), and processing means (7) configured to a) calculate phase shifts (Δϕ.sub.12, Δϕ.sub.13, Δϕ.sub.23) between the output signals (.sub.νPP.1, νPP.2, νPP.3) of the detectors (1, 2, 3), and b) calculate energy (E) of the pulse based on the calculated phase shifts.

2. A system according to claim 1, wherein said step a) comprises detecting a frequency (f.sub.g) of the output signals (.sub.νPP.1, νPP.2, νPP.3) of the detectors (1, 2, 3), and performing an I/Q method on the output signals (.sub.νPP.1, νPP.2, νPP.3), based on the detected frequency (f.sub.g), to calculate the amplitude (A.sub.PP.1, A.sub.PP.2, A.sub.PP.3) and phase (ϕ.sub.PP.1, ϕ.sub.PP.2, ϕ.sub.PP.3) of each output signal (.sub.νPP.1, νPP.2, νPP.3).

3. A system according to claim 1, wherein said detectors are detectors responsive to an electric field or magnetic field of the pulsed hadron beam passing thereby.

4. A system according to claim 1, wherein said detectors are detectors intercepting a fraction of the pulsed hadron beam.

5. A system according to claim 1, wherein the repetition rate of the bunches (B) is of an order of magnitude comprised between 100 MHz and 3 GHz, and preferably comprised between 100 MHz and 1 GHz.

6. Radiotherapy apparatus comprising at least one linear accelerator configured to produce and accelerate a hadron beam, and further comprising a beam energy measurement system according to claim 1.

7. Apparatus according to claim 6, wherein said linear accelerator is configured to produce and accelerate a proton beam.

8. A time-of-flight (TOF) measurement method for measuring energy of a pulsed hadron beam, wherein each pulse of the beam is structured into a series of bunches (B) of charged particles, said bunches being repeated according to a repetition rate of the order of magnitude of radiofrequency, wherein a first detector (1), a second detector (2) and a third detector (3) arranged along a beam path (10) are used, each of said detectors being configured to detect the passage of a bunch (B) of charged particles and provide an output signal (.sub.νPP.1, νPP.2, νPP.3) dependent on phase of the detected bunch (B), wherein the second detector (2) is spaced apart from the first detector (1) by a first distance (L.sub.12) and wherein the third detector (3) is spaced apart from the second detector (2) by a second distance (L.sub.23), wherein said first distance is set out in such a way as that time of flight (t.sub.12) of the bunch (B) from the first detector (1) to the second detector (2) is approximately equal to, or lower than a repetition period (T.sub.RFQ) of the bunches (B), and wherein said second distance is set out in such a way as that time of flight (T.sub.23) of the bunch (B) from the second detector (2) to the third detector (3) is greater than a multiple of the repetition period (T.sub.RFQ) of the bunches (B), and wherein said method comprises: a) calculating phase shifts (Δϕ.sub.12, Δϕ.sub.13, Δϕ.sub.23) between the output signals (.sub.νPP.1, νPP.2, νPP.3) of the detectors (1, 2, 3), and b) calculating energy (E) of the pulse based on the calculated phase shifts.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0030] Some preferred, but non-limiting, embodiments of the invention will now be described, with reference to the attached drawings, in which:

[0031] FIG. 1 shows a diagram of a detector arrangement of a beam energy measurement system according to the invention,

[0032] FIG. 2 shows a graph reporting limits on the distance and phase-shift errors for achieving 0.03% precision in the energy measurement at E=230 MeV,

[0033] FIG. 3 shows a hardware design of a prototype system according to the invention,

[0034] FIGS. 4 to 9 show flow diagrams of an exemplary beam energy measurement method, and

[0035] FIG. 10 shows a HEBT in which the invention may be situated.

DETAILED DESCRIPTION

[0036] With reference to FIG. 1, the energy measurement system according to the invention comprises a first detector 1, a second detector 2 and a third detector 3 arranged along a beam path 10, for example along a beam pipe of an accelerator HEBT. Each of the detectors 1, 2, 3 is configured to detect the passage of a bunch of charged particles and provide an output signal dependent on the phase of the detected bunch. Beam bunches are designated with B in FIG. 1. Each separate train of bunches constitutes a pulse. The second detector 2 is spaced apart from the first detector 1 by a first distance L.sub.12. The third detector 3 is spaced apart from the second detector 2 by a second distance L.sub.23. Distance between the first detector 1 and the third detector 3 is designated with L.sub.13.

[0037] In an embodiment the detectors 1-3 are capacitive pickups and in a specific embodiment these are phase probes. In place of phase probes a beam position monitor, a beam current transformer or a wall current monitor might be used. Alternatively a resonant cavity or an electro-optic crystal may be used. In general any device which measures the electric or magnetic fields of the particle beam is suitable. In an alternative embodiment a beam loss monitor or a device that intercepts part of the beam halo could be used.

[0038] Phase probes are capacitive sensors that can be used to detect in a non-interceptive way the passage of a bunch of charged particles. Their main component is a metallic ring, placed around the beam or beam pipe, on which a charge develops when a beam bunch passes inside it. This charge can be collected to get a current proportional to the variation of charge inside the ring.

[0039] t.sub.12 is the time taken by a particle bunch B to travel the distance L.sub.12, which can be used to compute the particles energy:

[00001]$\beta =\frac{L_{12}}{t_{12}}.Math.\frac{1}{c}$$\gamma =\frac{1}{\sqrt{1-{\beta }^2}}$$E=E_0.Math.\left(\gamma -1\right)$

[0040] where E is the kinetic energy of the particle and E.sub.0 is the rest energy of the particle (for protons: E.sub.0=938.272 MeV); c is the speed of light.

[0041] The system is designed to measure the phase shift Δϕ between the output signal of the probes 1-3. To be able to compute t.sub.12 from Δϕ.sub.12 (i.e. the phase shift between the output signal of the first detector 1 and the output signal of the second detector 2) the relation between the two has to be unambiguous. To this end, the first distance L.sub.12 is set out in such a way that the time of flight t.sub.12 of a bunch B from the first detector 1 to the second detector 2 is equal to, or lower than a repetition period T.sub.RFQ of the bunches B. This poses a limit on the maximum value of L.sub.12. For an energy range from 5 MeV to 230 MeV this limit is around 48 mm.

[0042] In a particular example the detectors of the invention are situated in a HEBT layout with distance L.sub.12=255 mm and distance L.sub.13=3595 mm. These distances provide a 0.03% E resolution for beams ranging from 70 up to 230 MeV, as shown in FIG. 2.

[0043] Given the limit on L.sub.12, it is impossible to achieve 0.03% of relative error on the measurement of E with only two close probes, having only one bunch travelling through them, because this would require a precision in the phase shift measurement which is nowadays unachievable. This is the reason behind the use of a third probe. L.sub.23 is much greater than L.sub.12, so more than one bunch can be positioned along L.sub.23. In other words, time of flight T.sub.2 of the bunch from the second detector 2 to the third detector 3 is greater than a multiple of the repetition period T.sub.RFQ of the bunches. N.sub.13 and N.sub.23 represent the number of whole bunches present between, respectively, detectors 1 and 3, 2 and 3. In an embodiment, the repetition rate T.sub.RFQ is given by the RFQ of the linear accelerator. Note that using only two distant probes is not sufficient as this would not allow an unambiguous energy measurement over the range from 5 MeV to 230 MeV. This is because using the two distant probes only (for example 1 and 3 of FIG. 1) the train of bunches (typically around 1000 of them into a beam pulse) passing through the phase probes, induces a train of signals reaching the acquisition system in the same acquisition window (typically 1 microsecond). Excluding the tails of the bunched pulse (affected by noise), it is not possible to recognize which is the integer number of bunches to skip before computing the ToF (N.sub.13 of FIG. 1.). We could for example in this instance assume N.sub.13+1 or N.sub.13−1, resulting then in a wrong measurement of energy. If instead an extra detector between detectors 1 and 3 of FIG. 1 is included as previously explained, it is possible to determine the approximate energy and thus N.sub.13.

[0044] Expressed in another way, if only two detectors were used placed so closely together that only one bunch was present in the inter-detector beamline at any one instant, then the measurement made by both detectors could be interpreted unambiguously as the measurement of the same bunch. However the measurement made, while unambiguous, would be inaccurate because of the phase difference. In fact the real ToF between the two extremely closely spaced detectors would be extremely small and the relative error would be large (calculated from error in measurement=error of instrument/distance between detectors). If we increase the distance between the two detectors so that more than one bunch can now fit simultaneously between the two detectors then the measurement will now result in the detection of two trains of bunches, shifted from each other, and it will not be possible to predict which detected bunch in each train should become the basis of the energy measurement, or in other words how many bunches N should be skipped. To overcome this problem, i.e. to know the correct value of N, it is necessary to already know the energy of the bunches, but calculating energy is the point of the measurement so we face a conundrum.

[0045] As an example of this, if the energy of the beam is 100 MeV and:

[0046] L.sub.12=225 mm,

[0047] L.sub.13=3595 mm,

[0048] delta_L=0.1 mm, and

[0049] delta_phi=0.2 deg

[0050] implies:

[0051] energy_error_12=0.12%,

[0052] energy_error_13=0.01%

[0053] In an alternative embodiment, if:

[0054] L.sub.12=40 mm,

[0055] L.sub.13=1000 mm,

[0056] delta_L=0.1 mm, and

[0057] delta_phi=0.2 deg

[0058] implies:

[0059] energy_error_12=0.7%,

[0060] energy_error_13=0.025%

[0061] However, a worst case occurs if the energy is 230 MeV (a maximum energy in some systems), in which case the corresponding error values are:

[0062] in the first case.

[0063] energy_error_12 is 0.21% (instead of 0.12%),

[0064] energy_error_13 is still 0.01%;

[0065] in the second case,

[0066] energy_error_12 is 1.16% (instead of 0.7%),

[0067] energy_error_13 is 0.05% (instead of 0.025%).

[0068] Therefore we measure the approximate beam energy using the signals from detectors 1 and 2 (which allow unambiguous but inaccurate calculation of energy) and use this approximation to calculate the number N of bunches which must be skipped to allow an unambiguous measurement between detectors 1 and 3 (which suffer from ambiguity but provide for a more accurate calculation of energy). By doing this we simultaneously reduce inaccuracy while maintaining unambiguity of calculation. Therefore three detectors are needed to produce a measurement which is both unambiguous and accurate.

[0069] This layout greatly improves the energy measurement precision; given the precision in the distance measurement δL and in the phase shift measurement δΔϕ, the relative error on the energy measurement using only detectors 1 and 2 is

[00002]$\frac{{\beta }^2{\gamma }^3}{\gamma -1}.Math.\sqrt{{\left(\frac{\delta L}{L_{12}}\right)}^2+{\left(\frac{\delta \mathrm{\Delta \phi }}{2\pi }.Math.\frac{T_{\mathrm{RFQ}}}{t_{12}}\right)}^2}$

[0070] While when using also detector 3 it is

[00003]$\frac{{\beta }^2{\gamma }^3}{\gamma -1}.Math.\sqrt{{\left(\frac{\delta L}{L_{13}}\right)}^2+{\left(\frac{\delta \mathrm{\Delta \phi }}{2\pi }.Math.\frac{T_{\mathrm{RFQ}}}{t_{13}}\right)}^2}$

[0071] which can be reduced by increasing L.sub.13 (T.sub.13 also increases consequently). The same reasoning can be applied to the opposite situation, i.e. in the case in which they are arranged in ‘reverse order’, and in this case the distance L.sub.23 is set out in such a way as that time of flight t.sub.23 of the bunch B from the second detector 2 to the third detector 3 is equal to, or lower than a repetition period T.sub.RFQ of the bunches B, and then having L.sub.12 much greater than L.sub.23 such that the time of flight T.sub.12 of the bunch from the first detector 1 to the second detector 2 is greater than a multiple of the repetition period T.sub.RFQ of the bunches. In such a case N.sub.23 will always be 0 and N.sub.12 has to be used in its place.

[0072] FIG. 2 shows another important aspect of the previous formulas; in the design phase of the accelerator layout the distance L.sub.13 should be chosen in such way that the measurement errors on both the phase and the distance still allow to achieve the beam energy resolution needed by a medical accelerator. δΔϕ is given by the electronics and from eventual interferences from high power sources in the accelerator, while δL is the result of the accuracy of survey instruments commonly used in the accelerator alignment, like a laser tracker (typical accuracy better than ±50 μm) combined with the uncertainty on the ring electrode alignment with respect to external references (typical accuracy ±50 μm). FIG. 2 depicts the upper limits on δL and δΔφ for achieving 0.03% precision in the energy measurement of a beam at 230 MeV at different values for L.sub.13. The allowed combinations of δΔϕ and δL are those sitting below the curves, and the dot corresponding to the actual system performance is indicated.

[0073] An example of hardware design which allows to measure the phase shifts between the output signals of the phase probes is shown in FIG. 3. The signals formed on the detectors 1, 2, 3 at beam passage enter limiters, devices used only to protect downstream electronics from unwanted spikes on the signals and successively they pass through RF coaxial relays, where they are routed to preamplification stages (at 5), where it is eventually possible to remotely change the gain. The amplified signals are then mixed down from 750 MHz to about 50 MHz (at 6), before being transmitted along coaxial cables that exit the machine room and enters directly a fast ADC plus FPGA (at 7) for acquisition and processing. The down-mix is not necessary, it is an implementation choice to relax requirements on the ADC and processing. The coaxial relays above mentioned are used to inject at the beginning of the three electronic chains a calibration signal at the same frequency of the bunches repetition rate (1/T.sub.RFQ), so to acquire and subtract in the data processing the fixed contribution in phase offset between the three channels due to mismatches in the cable lengths and in the installed electronics blocks (amplifiers, mixers, etc). The calibration of the system should run only once for a few milliseconds and it can be repeated any time during ‘no beam’ operation (for example before a new treatment starts) so to eventually update the phase offset; in this way it is possible to compensate any kind of long-term effect on the system like temperature drifts in the accelerator room.

[0074] The diagrams in FIGS. 4-9 represent an example of the computation that takes place inside the FPGA 7, which extracts the energy information from the detectors signals. The algorithm can be resumed in the following steps (FIGS. 4 and 5): [0075] Acquire the signals ν.sub.PP.1, ν.sub.PP.2, ν.sub.PP.3 coming from the detectors 1-3 (step 20). The exact time at which the acquisition should start (step 10 is computed based on a system-level trigger signal (which tells when the next pulse will occur). In a further embodiment it may also be based on the requirements that the detector signals are acquired for a duration which is slightly shorter than the pulse duration and starting slightly after the pulse so that the central part of the actual pulse is sampled. [0076] Perform a time-of-flight analysis of the signals (step 30). As will be explained in the following, this analysis comprises: [0077] Detect the exact frequency f.sub.g of the signals ν.sub.PP.1, ν.sub.PP.2, ν.sub.PP.3. This is necessary because the down-mixing introduces some uncertainty in the signal frequency. [0078] Use the detected frequency f.sub.g to perform an I/Q method on the signals, which gives the amplitude A.sub.PP.1, A.sub.PP.2, A.sub.PP.3 and the phase φ.sub.PP.1, φ.sub.PP.2, φ.sub.PP.3 of each signal. [0079] Check the signals amplitudes A.sub.PP.1, A.sub.PP.2, A.sub.PP.3 to see if the pulse has been correctly detected: If not, stop. [0080] Use the signals phases φ.sub.PP.1, φ.sub.PP.2, φ.sub.PP.3 to compute the energy E, and broadcast this information across the control system (step 40).

[0081] In FIGS. 4-9, the following blocks are used: [0082] Start/Stop: Where the (sub)algorithm flow starts or stops. [0083] Input/Read/Acquire: n the main diagram, this means that a signal is acquired by the hardware. In sub-diagrams, this means that a variable that has previously been set by the caller is required for this sub-algorithm to work. [0084] Output/Write/Send: In the main diagram, this means that some kind of value is transmitted by the algorithm to some other processing unit that may be interested in that value. In sub-diagrams, this means that a variable that has been set by the sub-algorithm is made available to the caller. [0085] Computation: Perform some kind of computation. [0086] Decision/Branching: A point in the diagram where the algorithm flow can take different paths based on a predicate. [0087] Sub-diagram: Execute the specified sub-diagram.

[0088] Every variable should be set before it is used or alternatively it is set as an a priori known value; the latter is true for constants known from physics and for the following variables: A.sub.min: Minimum amplitude required for the beam pulse to be considered correctly detected. f.sub.sampling: Sampling frequency.

[0089] A further explanation is required for the correct interpretation of the flow diagrams: Different flows enclosed by black horizontal lines represents operations that can be performed in parallel.

[0090] A label put at the top of a parallel branch has to be considered as additional subscript to every variable that appears in that branch, including input and output variables in sub-diagrams. If the label is at the bottom, only output variables gains the subscript.

[0091] With reference to FIG. 5, the frequencies f.sub.PP.1, f.sub.PP.2, f.sub.PP.3 of the output signals of the detectors 1-3 are detected (at 100).

[0092] FIG. 6 show this detection procedure in detail. The output signal ν.sub.PP of each detector is sampled with the sampling frequency f.sub.sampling (at 101). Then, the following variables are

[0093] calculated at 102:

[0094] N=number of samples in ν.sub.PP

[0095] G=Fast Fourier Transform of ν.sub.PP.

[0096] b=arg max.sub.i|G(i)| between f.sub.min and f.sub.max, wherein f.sub.min and f.sub.max are minimum and maximum values which can be set to constrain the search for the maximum in the Fast Fourier Transform of ν.sub.PP. They can be used when unknown frequencies are present in the Transform, although this should not be the normal situation.

[00004]$\beta =b+\frac{N}{\pi }.Math.\mathrm{atan}\left[\frac{\sin \left(\frac{\pi }{N}\right)}{\cos \left(\frac{\pi }{N}\right)+\frac{.Math.G_b.Math.}{.Math.G_{b+1}.Math.}}\right].$

[0097] Then, the frequency f.sub.PP of each signal is calculated as

[00005]$f_{\mathrm{PP}}=\frac{\beta }{N}.Math.f_{\mathrm{sampling}}.$

[0098] The frequency f.sub.g of the signal is then calculated (at 200 in FIG. 5) as f.sub.g=mean(f.sub.PP.1, f.sub.PP.2, f.sub.PP.3).

[0099] Then, an I/Q method is performed on each of the output signals ν.sub.PP (at 300). The I/Q method is shown in detail in FIG. 7. The I and Q signals are calculated at 301 as

I=Σ.sub.n=0.sup.N-1ν.sub.PP(n).Math.sin(2πf.sub.gnT.sub.s),

Q=Σ.sub.n=0.sup.N-1ν.sub.PP(n).Math.cos(2πf.sub.gnT.sub.s).

[0100] where T.sub.s=f.sup.1.sub.sampling is the sampling period.

[0101] Then, phase φ.sub.PP and amplitude A.sub.PP of each signal are calculated at 302 as

[00006]${\phi }_{\mathrm{PP}}=\mathrm{atan}\left(\frac{Q}{I}\right),A_{\mathrm{PP}}=\frac{2}{N}.Math.\sqrt{Q^2+I^2}.$

[0102] The signal amplitudes A.sub.PP are then compared with A.sub.min (at 400 in FIG. 5). If at least one of these amplitudes is lower than A.sub.min, the process is stopped.

[0103] Otherwise, the phase shifts Δφ of the output signals are calculated at 500, and subjected to wrapping at 600 (see also FIG. 8).

[0104] Then, energy values E.sub.13 and E.sub.23 of the pulse is calculated (at 700) based on time-of-flight measurements between detectors 1 and 3, and between detectors 2 and 3, respectively. This calculation is shown in detail in FIG. 9. The following variables are calculated at 701 (here, the subscript C takes place of subscript 12, while subscript F takes place of subscript 13 or 23, depending on whether E.sub.13 or E.sub.23 is calculated; T.sub.RF corresponds to T.sub.RFQ mentioned above):

[00007]$N=\left[\left(\frac{L_F}{L_C}.Math.{\mathrm{\Delta \phi }}_C-{\mathrm{\Delta \phi }}_F\right).Math.\frac{1}{2\pi }\right],T_F=T_{\mathrm{RF}}.Math.\left(N+\frac{{\mathrm{\Delta \phi }}_F}{2\pi }\right),\beta =\frac{L_F}{T_F}.Math.\frac{1}{c},\gamma =\frac{1}{\sqrt{1-{\beta }^2}},E_{\left(13\mathrm{or}23\right)}=A.Math.E_0.Math.\left(\gamma -1\right),\frac{\delta E}{.Math.E.Math.}=\frac{{\beta }^2{\gamma }^3}{\gamma -1}.Math.\sqrt{{\left(\frac{\delta L}{L_F}\right)}^2+{\left(\frac{\mathrm{\delta \Delta }\phi }{{\mathrm{\Delta \phi }}_C}.Math.\frac{L_C}{L_F}\right)}^2}$

[0105] The energy E of the pulse is calculated as a mean value between E.sub.13 and E.sub.23 (at 800).

[0106] According to an alternative embodiment, it would be sufficient to use E.sub.13 or E.sub.23 to provide the beam energy. The mean value between E.sub.13 and E.sub.23 is used to improve the measurement accuracy. The accuracy might be even further improved using a fourth detector/phase probe or more, but this would add complexity to the system.

[0107] The ToF beam energy measurement system allows a high accuracy beam energy measurement to be made at a very high measurement rate (up to 200 Hz) and provides the result typically within 1 ms from the passage of the beam pulse, making it suitable to be used not only as Beam Diagnostics device but also in the Beam Delivery System to monitor each beam pulse average energy, which has been delivered to the patient.

[0108] Such a highly responsive system is fast enough to allow the system to take actions to disable generation of the next beam pulse.

[0109] The system according to the invention makes no assumptions about the speed at which the beam energy can be changed, thus it poses no restrictions on the energy change rate. This is an improvement over the state-of-the-art because current beam energy measurement systems are either destructive or they do not allow the measurement of fast beam energy changes.

[0110] A particular embodiment of the invention is shown in FIG. 10, which shows a HEBT in which the invention may be situated. The last accelerating unit of the linear accelerator (901) issues the proton beam which comprises pulses, each of which comprises bunches, which bunches then pass the first detector (1), the second detector (2) and finally the third detector (3), all of which measure a bunch of the proton beam.

[0111] In a particular embodiment the last accelerating unit may be a CCL.

[0112] The detectors (1, 2, 3) share space in the HEBT with other components (902, 903) for example quadrupoles, ACCTs, BPMs, vacuum pumps, etc. The actual components present will depend on the particular HEBT layout, which will depend on the particular geometry of the installation.

[0113] The distances between the detectors are:

[0114] L.sub.12=255 mm and

[0115] L.sub.13=3595 mm.

[0116] These distances allow achievement of 0.03% E resolution for beams ranging from 70 up to 230 MeV, as shown in FIG. 2.

[0117] After passing through the third detector (3) the proton pulse will continue along the remainder of the HEBT (1000) which leads the beam towards a nozzle and then, finally, the patient placed in a treatment room.

REFERENCES

[0118] 1 first detector [0119] 2 second detector [0120] 3 third detector [0121] 5 amplifier [0122] 6 mixer and filter [0123] 7 FPGA [0124] 10 procedural step [0125] 20 procedural step [0126] 30 procedural step [0127] 40 procedural step [0128] 100 detection step [0129] 101 sampling step [0130] 102 calculation step [0131] 200 frequency calculation step [0132] 300 I/Q method step [0133] 301 I, Q signal calculation step [0134] 302 phase and amplitude calculation step [0135] 400 comparison step [0136] 500 phase shift calculation step [0137] 600 wrapping step [0138] 700 Energy value calculation step [0139] 701 calculation step [0140] 800 mean value calculation step