COMPUTER IMPLEMENTED METHOD FOR REDUCING THE RISK OF INTERRUPTING AN IRRADIATION TREATMENT SESSION DUE TO A DEVIATION FROM A PLANNED VALUE OF AN OPERATING PARAMETER OF A PARTICLE ACCELERATING SYSTEM

Abstract

A computer implemented method for optimizing tolerance values of operating parameters of a particle accelerating system allowing a plurality of beamlets of particles accelerated along an irradiation axis, to deposit doses by pencil beam scanning into a structure of interest of a patient according to a treatment plan. The method calculates the dose (rate) volume histograms for a statistically representative number N of values randomly selected within a defined confidence level in preselected tentative statistical distributions of the operating parameters and compares the obtained calculated dose (rate) volume histogram with an acceptable band of variation of a target dose (rate) volume histogram. Once a tentative statistical distribution yields N calculated dose (rate) volume histograms which all fall within the acceptable band of variation, it is set as the final statistical distribution, and the particle accelerating system can be programmed with the final statistical distribution.

Claims

1. A method for optimizing tolerance values of operating parameters of a particle accelerating system providing a beam formed by a plurality of beamlets of particles accelerated along an irradiation axis, to deposit doses by pencil beam scanning into a structure of interest of a patient according to a treatment plan, the method comprising: (a) providing an input including: a definition of an array of beamlets in the treatment plan, wherein the definition comprises a plurality of planned parameters, including: a planned position of each beamlet on a plane normal to the irradiation axis; a planned monitor of each beamlet; and a planned beamlet scanning sequence over the planned positions; a planned starting time and end time at which each beamlet is to be delivered; a definition of the structure of interest, defining one or more tissues being traversed by a number of the beamlets; values of one or more target dose (rate) distribution histograms for the structure of interest obtained with a treatment with the planned parameters; and acceptable bands of variation within which the one or more target dose (rate) distribution histograms are allowed to vary; (b) providing a tentative statistical distribution of operating parameters of the particle accelerating system configured to center on corresponding average values representative of the performance of the particle accelerating system and define a confidence level of the tentative statistical distribution, wherein the operating parameters comprise, the monitor of each beamlet; the position of each beamlet; and the starting time and end time of delivery of each beamlet; (c) randomly selecting from the corresponding tentative statistical distributions within the predefined confidence levels a value of the monitor, a value of the position of the beamlet, and a value of each of the starting time and end time; (d) calculating one or more calculated dose (rate) distribution histograms with the values randomly selected; (e) repeating steps (c) and (d) a predetermined number of times to yield calculated distributions characterizing the one or more calculated dose (rate) distribution histograms for all the values randomly selected of the operating parameters; comparing the calculated distributions of the one or more calculated dose (rate) distribution histograms with the corresponding acceptable bands of variation; and (g) determining whether the calculated distributions of the one or more calculated dose (rate) distribution histograms are comprised in the corresponding acceptable bands of variation within the pre-defined confidence level.

2. The method of claim 1, further comprising setting a final statistical distribution with a final confidence level for the corresponding operating parameters, wherein: if the calculated distributions of the one or more calculated dose (rate) distribution histograms are all comprised within the corresponding acceptable bands of variation with the pre-defined confidence level, for the given treatment plan, setting the tentative statistical distribution as the final statistical distribution and the confidence level as the corresponding final confidence level to define the corresponding operating parameters; if any one of the one or more calculated dose (rate) distribution histograms calculated with one set of randomly selected values within the corresponding confidence levels of the statistical distribution of the operating parameters extends beyond the corresponding acceptable bands of variation, repeating steps (b) to (f), with new tentative statistical distributions of the operating parameters and/or selecting new confidence levels, until the calculated distributions of the one or more calculated dose (rate) distribution histograms are all comprised within the corresponding acceptable bands of variation, setting the new tentative statistical distribution as the final statistical distribution, and setting the new corresponding confidence level as the final confidence level to define the corresponding operating parameters.

3. The method of claim 2, wherein setting the final statistical distribution (Tf) to the corresponding operating parameters is carried out by a human operator or by a processor.

4. The method of claim 2, wherein the new tentative statistical distributions have lower standard deviations than the original tentative statistical distributions.

5. The method of claim 1, wherein the tentative statistical distributions of each of the operating parameters are Gaussian distributions; and the value of the confidence level is comprised between 68% and 99.7% of the tentative statistical distribution, wherein a confidence level of 68% corresponds to μj±σj, a confidence level of 95% corresponds to μj±2σj, and a confidence level of 99.7% corresponds to μj±3σj, wherein μj is an average value and σj is the standard deviation of the corresponding tentative statistical distributions.

6. The method of claim 5, wherein the average values and standard deviations of each operational parameter are different for each beamlet.

7. The method of claim 5, wherein the value of the confidence level is comprised between 95.5% and 99% of the tentative statistical distribution.

8. The method of claim 1, wherein the particle accelerating system is equipped with a cyclic checker configured to measure at different intervals, or continuously, actual values of the operating parameters including the monitor, the position, and of the starting time and end time of the beamlets emitted by the particle accelerating system.

9. The method of claim 8, wherein the particle accelerating system includes a processor configured to compare the actual values of the operating parameters with the corresponding confidence level and to stop a treatment session when one actual value of an operating parameter falls outside of the corresponding final confidence level.

10. The method of claim 1, wherein the planned parameters also comprise a planned beamlet size, whose values used for calculating the one or more calculated dose (rate) distribution histograms are randomly selected within a corresponding tentative statistical distribution of the beamlet size.

11. The method of claim 1, wherein the calculated distributions of the one or more calculated dose (rate) distribution histograms are defined by a corresponding area comprised between a minimum calculated dose (rate) distribution histogram and a maximum calculated dose (rate) distribution histogram, wherein the minimum calculated dose (rate) distribution histogram is defined by the lowest values of the calculated dose (rate) distribution histograms calculated with the predefined confidence level from the N randomly selected values of the monitor, position of the beamlets, and starting time and end time; and the maximum calculated dose (rate) distribution histogram is defined by the highest values of the calculated dose (rate) distribution histograms calculated with the predefined confidence level from the N randomly selected values of the monitor, position of the beamlets, and starting time and end time.

12. The method of claim 1, wherein the treatment plan includes depositing doses into at least a portion of the structure of interest at ultra-high deposition rates (UHDR) defined as a deposition rate greater than or equal to 1 Gy/s.

13. The method of claim 1, wherein the planned parameters also comprise a planned beam current, whose values used for calculating the calculated dose (rate) distribution histogram are randomly selected within a corresponding tentative statistical distribution of the beam current.

14. The method of claim 1, wherein the dose distribution histogram is a dose volume histogram, and wherein the dose rate distribution histogram is a dose rate volume histogram or a differential dose rate histogram.

15. An error predictor for implementing the method of claim 1, the error predictor comprising: a memory storing a plurality of tentative statistical distributions for each operating parameter centered on a plurality of corresponding average values; a user interface configured to receive user input to: enter the treatment plan including one or more target dose (rate) distribution histograms and the corresponding acceptable bands of variation; select from the memory or enter a planned starting time and end time at which each beamlet is to be delivered; select from the memory or enter for each beamlet, a first tentative statistical distribution for each operating parameter, including the monitor, the position of the beamlet, and starting time and end time; and enter a confidence level on the operational parameters; a processor configured to: (i) randomly select a value of the monitor, a value of the position of the beamlet, values of the starting time and end time comprised within the predefined confidence levels of the corresponding tentative statistical distributions; (ii) calculate a calculated one or more of a dose distribution histogram and a calculated dose rate distribution histogram with the values randomly selected; and (iii) repeate steps (i) and (ii) a predetermined number of times to yield the calculated distributions of the calculated dose distribution histogram and calculated dose rate distribution histogram for each beamlet.

16. The error predictor of claim 15, wherein the processor is further configured, if any one of the one or more of the calculated distributions of a corresponding calculated dose rate distribution histogram are not comprised within the corresponding acceptable bands of variation with the pre-defined confidence level, to repeat (i) to (iii) with new tentative statistical distributions of the operating parameters, until the calculated distributions of the one or more of calculated dose rate distribution histogram are all comprised within the corresponding acceptable bands of variation with the pre-defined confidence level.

Description

BRIEF DESCRIPTION OF THE FIGURES

[0063] For a fuller understanding of the nature of the present disclosure, reference is made to the following detailed description taken in conjunction with the accompanying drawings in which:

[0064] FIGS. 1(a) to 1(c) show examples of tDVH, tDRVH, and tDDRH curves with their corresponding acceptable bands of variation (BV).

[0065] FIG. 2 shows a Gaussian distribution of an operating parameter of a given particle accelerating system.

[0066] FIGS. 3(a) to 3(c) show various steps of one embodiment of the present disclosure, with a positive result.

[0067] FIGS. 4(a) to 4(c) show various steps of one embodiment of the present disclosure, with a negative result.

[0068] FIGS. 5(a) to 5(d) show N runs of calculations yielding calculated distributions (CDj) of the calculated cD(R)DHj shown in FIGS. 3(b) and 3(c).

[0069] FIG. 6 shows a flowchart with various steps of a method of the present disclosure.

DETAILED DESCRIPTION

[0070] The present disclosure provides a computer implemented method and an error predicting module which considerably reduce the risk of carrying out a treatment session which does not respect a corresponding treatment plan, for reasons of equipment.

[0071] The present disclosure also provides a computer implemented method for optimizing tolerance values of operating parameters of a particle accelerating system allowing a beam formed by a plurality of beamlets of accelerated particles to deposit doses by pencil beam scanning (=PBS) into a structure of interest of a patient according to a treatment plan (=TP). The method allows determining the confidence level (CLj) for a given particle accelerating system to deliver beamlets which will fulfil the TP as a function of a set of operating parameters. If the confidence level obtained is too low, an alternative set of operating parameters needs be evaluated. These and other advantages of the present disclosure are presented below.

Method for Optimizing Tolerance Values of Operating Parameters

[0072] The present disclosure provides a method implemented by a computer for optimizing tolerance values of operating parameters of a particle accelerating system allowing a beam formed by a plurality of beamlets of particles accelerated along an irradiation axis (Z), to deposit doses by pencil beam scanning (=PBS) to a patient according to a treatment plan (=TP). The particles may be protons, but they can be electrons, heavy ions beams, but also waves formed by accelerated particles interacting with a converting material, such as x-rays (or y-rays). As illustrated in FIG. 6, the computer implemented method receives the input of a number of values of planned parameters. These can be provided by a treatment plan (TP) including the following planned parameters, [0073] a planned position (Xpi) of each beamlet (bi) on a plane (X, Y) normal to the irradiation axis (Z), [0074] a planned monitor unit (MUpi) of each beamlet, [0075] a planned beamlet scanning sequence over the planned positions (Xpi),

[0076] The method also receives a planned starting time (t0pi) and end time (t1pi) at which each beamlet is to be delivered. This is particularly useful in case of FLASH-RT. The starting and end times may or may not be part of the TP.

[0077] A structure of interest is defined, characterizing one or more tissues, of the patient, being traversed by or interacting with the beamlets (bi). The structure of interest comprises the target comprising the tumoral cells to be killed, but also healthy tissues traversed or somehow touched by one or more beamlets. These include, for example, the tissues located upstream from the target along the irradiation axis (Z), i.e., between the nozzle of the particle accelerating system and the target, or also the tissues adjacent to and surrounding the target.

[0078] An objective of the treatment is to yield at the end of a session, given values of a target dose (rate) distribution histogram (tD(R)DH) for the structure of interest, within an acceptable band of variation (BV) in which a target dose (rate) volume histogram tD(R)VH is allowed to vary. As shown in FIGS. 1(a) to 1(c), tD(R)DH can include for example a target dose volume histogram (=tDVH) (show in FIG. 1(a)), a target dose rate volume histogram (=tDRVH) (shown in FIG. 1(b)), or a target differential dose rate histogram (tDDRH) (shown in FIG. 1(c)). Note that other histogram representations are possible and included in the term tD(R)DH, to display the information on the dose and dose rate distributions. FIGS. 1(a) to 1(c) show examples of tDVH, tDRVH, and tDDRH represented by the solid lines, as well as of the corresponding acceptable bands of variation (BV) represented by the area comprised between the dashed lines. Similar acceptable bands of variation can be defined for any type of alternative dose (rate) distribution histograms (D(R)DH) representations. If the beamlets delivered by the particle accelerating system achieve a dose distribution histogram (DDH) and a dose rate distribution histogram (DDRH) comprised within the acceptable band of variation as shown in FIG. 1(c), the session is successful. In case the DDH or DRDH falls outside of the acceptable band of variation (BV), the treatment session must be interrupted, to the great discomfort of the patient who must wait for a free slot in the particle accelerator system's schedule to resume the treatment session. The present disclosure allows anticipating within a given confidence level the occurrence of such an event, as shown in FIGS. 4(b) and 4(c) and modifying the tentative statistical distribution of operating parameters of the particle accelerating system to reduce the probability of having to interrupt a session.

[0079] FIG. 1(c) includes the actual values of the differential dose rate distribution histogram (aDDRH) obtained upon measuring the beamlets parameters during a treatment session (shown as long-dashed line). It can be seen that the curve aDDRH in FIG. 1(c) is fully enclosed within the acceptable band of variation (BV). The treatment session was therefore successful. Had the curve aDDRH extended beyond the acceptable band of variation (BV), however, the treatment session would have had to be interrupted. In order to decrease the risk of having to interrupt a treatment session because of the actual dose (rate) distribution histograms (D(R)DH) falling outside of the acceptable band of variation, the present disclosure implements the following actions with a computer.

[0080] First, as shown in FIGS. 3(a) and 4(a), a tentative statistical distribution (Tj) is selected of operating parameters (MUj, Xj, t0j, t1j) of the particle accelerating system centered on corresponding average values (μj) representative of the performance of the particle accelerating system. A confidence level (CLj) of the tentative statistical distribution (Tj) of each of the operating parameters is defined. The operating parameters include: [0081] monitor unit (MUj) of each beamlet, [0082] position of each beamlet (Xj), and [0083] starting time (t0j) and end time (t1j) of delivery of each beamlet.

[0084] The selection of tentative statistical distributions (Tj) is based on the performance of the particle accelerating system and is expected to generate a beam delivering the beamlets scanning sequence and generating calculated dose volume histograms (=cDVHj) and calculated dose rate volume histograms (=cDRVHj) in the structure of interest that are comprised within the acceptable band of variation (BV). For example, a treatment plan system (TPS) can determine average values (μj) achievable by the particle accelerating system which would yield the desired cDVH and cDRVH. The actual operating parameters of the treatment machine on a specific day, however, are not restricted to the corresponding average values (μj) but are distributed generally over a Gaussian curve, varying from day to day or during the course of a day (shown in FIG. 2). The resulting calculated dose (rate) distribution histograms (cD(R)DH) will therefore vary depending on the actual values of the operating parameters at the specific time of the treatment and must therefore be calculated taking account of the distribution. It would not be practical to consider all the possible values of each operating parameter defined by the corresponding distribution curves. A confidence level (CLj) can be defined, restricting the distribution to boundaries which are considered as acceptable. For example, the confidence level can be defined with respect to the standard deviation (σj), such as for example, μj±nσj, with n=1 to 3 (shown in FIG. 2).

[0085] A value (MUij) of the monitor unit, a value (Xij) of the position of the beamlet, and values (t0ij, t1ij) of starting time (t0ij) and end time (t1ij) are randomly selected from the corresponding tentative statistical distributions (Tj) within the corresponding confidence levels (CLj) as shown with the black circles in the Gaussian distribution curves (Tj) of FIGS. 3(a) and 4(a). The values are selected within the confidence level (CLj) previously defined.

[0086] As shown in FIG. 5(a), in a first Run (=Run 1), the calculated dose (rate) distribution histogram (=cD(R)DHj) is calculated with the values randomly selected. If the calculated cD(R)DHj is enclosed within the corresponding acceptable band of variation, new values of the operating parameters are randomly selected from the tentative statistical distributions (Tj) and the corresponding cD(R)DHj is calculated in a second Run (=Run 2) as illustrated in FIG. 5(b). These operations are repeated a statistical number of times (N) as illustrated in FIG. 5(c), “Run N”, as long as the calculated cD(R)DH-Run x are enclosed within the acceptable band of variation (BV). After N runs, the N curves “cD(R)DHj-Run x” with x=1 to N define in combination a calculated distribution (CDj), represented in FIGS. 3(b), 3(c), 4(b), 4(c), and 5(d) by the shaded areas bounded by the dotted lines of cD(R)DHj0 and cD(R)DHj 1.

[0087] If, on the other hand, any one of the cD(R)DHj-Run x calculated extends beyond the acceptable band of variation (BV), it could be concluded that there would be a risk higher than the predefined confidence level (CLj) of having to interrupt a treatment run with the operating parameters of the particle accelerating system according to the tentative statistical distributions (Tj). A new tentative statistical distribution (T(j+1)) of the operating parameters is then selected, and the cD(R)DH(j+1) is calculated in a new series of N Runs with randomly selected values of the new tentative statistical distributions (T(j+1)) at each successive Run. This operation is repeated with new tentative statistical distributions (T(j+k)) until the corresponding calculated distributions (CD(j+k)) are entirely comprised within the acceptable band of variation (BV).

[0088] Besides the planned position (Xpi), the planned monitor unit (MUpi), the planned beamlet scanning sequence over the planned positions (Xpi), and the planned starting time (t0pi) and end time (t1pi), the planned parameters can also comprise a planned beamlet size (dj) and a beam current (Ij), whose values (dij, Iij)) used for calculating the calculated dose (rate) distribution histogram (=cD(R)DHj) are randomly selected within a corresponding tentative statistical distribution (Tj) of the beamlet size (dj) and beam current (Ij).

[0089] The calculated distributions (CDj) of the cD(R)DHj can be defined by a corresponding area comprised between, [0090] a minimum calculated dose (rate) distribution histogram (=cD(R)DHj0) and, [0091] a maximum calculated dose (rate) distribution histogram (=cD(R)DHj1),
wherein [0092] cD(R)DHj0 is defined by the lowest values of cD(R)DHj calculated with the predefined confidence level (CLj) from the N randomly selected values of the monitor unit (Muij), position (Xij) of the beamlets, and starting time and end time (t0ij, t1ij) and wherein [0093] cD(R)DHj1 is defined by the highest values of cD(R)DHj calculated with the predefined confidence level (CLj) from the N randomly selected values of the monitor unit (Muij), position (X0i) of the beamlets, and starting time and end time (t0ij, t1ij).

[0094] By comparing the calculated distributions (CDj) of cD(R)DHj with the acceptable band of variation (BV), it can be determined whether or not the calculated distribution (CDj) of cD(R)DHjj is comprised in the corresponding acceptable band of variation for the pre-defined confidence level (CLj). FIGS. 3(b) and 3(c) show an example of tentative statistical distributions (Tj) with associated confidence level (CLj) which yields calculated distributions (CDj) of cD(R)DHj comprised within the corresponding acceptable bands of variation (BV). By contrast, FIGS. 4(b) and 4(c) show an example of tentative statistical distribution (Tj) with associated confidence level (CLj) which yield calculated distributions (CDj) of cD(R)DHj extending beyond the corresponding acceptable bands of variation (BV) as indicated by the shaded areas and black arrows. For example, it can be seen in FIG. 4(b), that the cD(R)DHj-Run 1 already extends beyond the acceptable band of variation (BV). Calculation with the corresponding tentative statistical distribution can be stopped and a new tentative statistical distribution (T(j+1)) can be selected to repeat the N runs with values of the operating parameters randomly selected from the new tentative statistical distribution (T(j+1)).

[0095] The method of the present disclosure may conclude by setting a final statistical distribution (Tf) to the corresponding operating parameters with a final confidence level (CLf). If the N runs of calculated distributions (CDj) of cD(R)DHj are comprised within the corresponding acceptable band of variation (BV) with the pre-defined confidence level (CLj), for the given treatment plan (TP), the tentative statistical distribution (Tj) can be set as the final statistical distribution (i.e., Tf=Tj) with the final confidence level (CLf=CLj) to define the corresponding operating parameters. It can be concluded that, by implementing the final statistical distributions (Tf) of operating parameters, the particle accelerating system has a probability equal to the final confidence level (CU) of delivering the beamlets fulfilling the treatment plan (TP) with corresponding actual dose (rate) distribution histogram(aD(R)DH) comprised within the acceptable band of variation (BV), as shown e.g., in FIG. 1(c).

[0096] If, on the other hand, any one of the N runs of calculated dose (rate) distribution histogram (cD(R)DHj) extends beyond the boundaries of the corresponding acceptable bands of variation (BV) with the pre-defined confidence level (CLj) as shown e.g., in FIG. 4(b), a new tentative statistical distribution (T(j+1)) and/or a new confidence level (CL(j+1)) can be defined for one or more of the operating parameters and the cD(R)DHj can be calculated N times with N randomly selected values of each of the operating parameters. These operations can be repeated as often as necessary with new tentative statistical distributions (T(j+k)) of the operating parameters and/or with alternative confidence levels (CL(j+k)), until the calculated distributions (CD(j+k)) of cD(R)DHj(j+k) is comprised within the corresponding acceptable band of variation (BV) for the pre-defined confidence level (CLj) or (CL(j+k)). The final statistical distribution (Tf) can be set as being the tentative statistical distribution (T(j+k)) (i.e., Tf=T(j+k)) and the corresponding confidence level (CL(j+k)) can be set as the final confidence level (CLf=CL(j+k)), to define the corresponding operating parameters.

[0097] The new tentative statistical distributions (T(j+k)) defined in case the j+(k−1) preceding tentative statistical distributions did not yield calculated distributions included within the bands of variation (BV) surrounding target dose (rate) distribution histogram (tD(R)DH) for the pre-defined confidence level (CLj), can be selected as distributions having lower standard deviations (σj) (or variances (σj).sup.2) than the corresponding preceding tentative statistical distributions (Tj+(k−1)).

[0098] Setting the final statistical distribution (Tf) for the corresponding operating parameters can be carried out by a human operator or automatically, by a processor.

[0099] The tentative statistical distributions (Tj) of each of the operating parameters may be Gaussian distributions. The values of the confidence levels (CL) may be comprised between 68% and 99.7%, or may be between 95.5 and 99% of the tentative statistical distribution. As shown in FIG. 2, a confidence level (CLj) of 68% corresponds to μj±σj. A confidence level (CLj) of 95% corresponds to μj±2σj and a confidence level of 99.7% of the tentative statistical distribution, corresponds to μj±3σj, wherein σj is the standard deviation of the corresponding tentative statistical distributions. For example, in terms of absolute deviations, the position (Xj) of a spot may, for example, vary by ±1 mm from the average value Rj of the Xj. The monitor unit (MUj) can for example vary by about 0.5% around the average value (Rj) of MUj. Note that the average values (Rj) and standard deviations (σj) of each operational parameter can be different for each beamlet (bi).

[0100] The particle accelerating system can be equipped with a cyclic check module configured for measuring at different intervals or continuously actual values of the monitor unit (MUai), the position (Xai), and of the starting time and end time (t0ai, t1ai) of the beamlets emitted by the particle accelerating system. A processor can be configured to determine whether any actual operating parameter falls outside of the corresponding final confidence level (CU) (shown in FIGS. 3(a) and 4(a), showing the confidence levels (CLj) spanning over a portion only of the statistical distributions (Tj)). In such case, an alarm can be triggered informing an operator of this event. As a safety measure, the processor can also be configured to stop the treatment as soon as an operating parameter falling outside of the confidence level (CLj) is detected.

[0101] The fact that one or more values of the actual operating parameters fall outside of the corresponding confidence levels (CLj) does not mean that the resulting actual dose (rate) distribution histogram (aD(R)DH) necessarily falls outside of the acceptable band of variation (BV). It could therefore be over-shooting to interrupt the treatment session simply because one actual value of any operating parameter falls outside of the corresponding confidence level (CLj), as it could perfectly yield a dose (rate) distribution histogram (D(R)DH) comprised within the acceptable band of variation (BV). The processor may be configured to calculate the cD(R)DH as soon as a measured actual value of an operating parameter falls outside of the corresponding confidence level (CLj) to determine whether or not the calculated cD(R)DH are comprised within the corresponding acceptable bands of variation (BV). The calculated dose (rate) distribution histogram (cD(R)DH) can be calculated based on the actual values of the operating parameters measured on the beamlets already delivered, including the parameter falling outside of the confidence level (CLj), and on the average values (Rj) of the operating parameters according to the final statistical distribution (Tf) for the beamlets which remain to be delivered to end the treatment session. If the calculated cD(R)DH falls outside of the acceptable band of variation (BV), the treatment session is stopped. If, on the other hand, the calculated cD(R)DH is within the acceptable band of variation (BV), the treatment session can proceed further. With this functionality of the processor, the risk of stopping a treatment session is further reduced.

[0102] This function does not take excessive calculating power, as it would concern only 100%−CLj % of the actual values of the operating parameters. For a confidence level (CLj) of 95% i.e., μj±2σj, there would be a probability of only 5% that a value of an operating parameter should fall outside of the confidence level (CLj) whch the D(R)DH would have to be calculated for. For a confidence level of 99.7% of the tentative statistical distribution, i.e., μj±3σj, it would concern a probability of merely 0.3% where such calculation would be required.

[0103] In an embodiment, the treatment plan includes FLASH-RT, in that doses are to be deposited into at least a portion of the structure of interest at ultra-high deposition rates (UHDR) defined as a deposition rate greater than or equal to 1 Gy/s. In this embodiment, the planned parameters also comprise a beam current (I), whose values (Ij) used for calculating the calculated dose distribution histogram (=cDDHj) and, in particular, the calculated dose rate distribution histogram (=cDRDHj) are randomly selected within a corresponding tentative statistical distribution (Tj) of the beam current (Ij).

Error Predicting Module

[0104] The present disclosure also provides an error predicting module configured to implement the method described above. The error predicting module comprises a memory comprising a plurality of tentative statistical distributions (Tj) for each operating parameter centered on a plurality of corresponding average values (Rj). It also comprises a user interface configured to enter, [0105] planned operators including a treatment plan (TP) including a planned position (Xpi), a planned monitor unit (MUpi), a planned beamlet scanning sequence, and planned starting time (t0pi) and end time (t1pi), [0106] the target dose (rate) distribution histogram (=tD(R)DH) and corresponding acceptable band of variation (BV), and [0107] for each beamlet, a first tentative statistical distribution (Tj) selected among or defined by entering values of the average values (μj) and standard deviations (σj) from one or more of the plurality of tentative statistical distributions for each operating parameter, including at least a monitor unit (MUD, a position (Xj) of the beamlet, and starting time and end time (t0i, t1i),

[0108] The error predicting module comprises a processor configured to, [0109] randomly select a value of the monitor unit (MUij), a value of the position (Xij) of the beamlet, values of the starting time and end time (t0ij, t1ij) comprised within the predefined confidence levels (CLj) of the corresponding tentative statistical distributions (Tj), [0110] calculate a calculated dose (rate) distribution histogram (=cD(R)Dj) with the values randomly selected, [0111] repeat the last two steps a statistically representative number of times (N) to yield the calculated distributions (CDj) of the calculated cD(R)Dj for each beamlet.

[0112] The processor can be further configured, in case the calculated distribution (CDj) of cD(R)Dj is not comprised within the corresponding acceptable band of variation with the pre-defined confidence level (CLj), to repeat the foregoing three steps with new tentative statistical distributions (T(j+k)) of the operating parameters, until the calculated distributions (CDj) of cD(R)Dj is comprised within the corresponding acceptable bands of variation with the pre-defined confidence level (CLj).

Flowchart (FIG. 6)

[0113] FIG. 6 shows a flowchart illustrating various steps of an embodiment of the present disclosure. The requirements of a treatment plan (TP) (shown in step (1) of FIG. 6) are translated into planned parameters for all the beamlets (shown in step (2) of FIG. 6). One or more target dose (rate) distribution histograms (tD(R)DH) are defined (shown in step (3) of FIG. 6) with their corresponding acceptable bands of variation (BV) (shown in step (4) of FIG. 6). As defined above, the dose (rate) distribution histograms (D(R)DH) can include non-exhaustively, the dose volume histogram (DVDH), the dose rate volume histogram (DRVH), or the differential dose rate distribution histogram (DDRH), and the like.

[0114] The particle accelerating system is simulated (shown in step (5) of FIG. 6), by defining a first tentative statistical distribution (Tj) and a corresponding confidence level (CLj) (show in steps (6) and (7) of FIG. 6). A value of each operational parameter is randomly selected within the confidence level (CLj) of the tentative statistical distributions (Tj) (shown in step (8) of FIG. 6) and the resulting cD(R)DHj-Run (x=1) is calculated (shown in step (9) of FIG. 6). If the calculated CD(R)DH-Run (x=1) is not within the acceptable band of variation (BV) a new tentative statistical distribution (T(j+1)) is selected with a corresponding confidence level (CL(j+1)) (show in steps (10) and (11) of FIG. 6), and steps (7) to (10) are repeated. If, on the other hand, the calculated CD(R)DH-Run (x=1) is within the acceptable band of variation (BV), a new value of each operational parameter is randomly selected a statistically representative number of times (N) within the confidence level (CLj) of the tentative statistical distributions (Tj) (shown in steps (11), (12) and (8) of FIG. 6) and steps (8) to (10) are repeated. When all N Runs x yielding calculated CD(R)DH-Run (x=1 to N) are within the acceptable band of variation (BV), the corresponding tentative statistical distribution (Tj) can be set as the final statistical distribution (Tf) (show in step (13) of FIG. 6). The particle accelerating system can be programmed with the operational parameters according to the corresponding final statistical distributions (Tf). With such programming, the particle accelerating system has a probability of CLj % of delivering beamlets satisfying the treatment plan without interruption of the treatment session.

[0115] The method of the present disclosure can be implemented with a single calculated dose (rate) histogram (cD(R)DH) or with several histograms which must all be within the corresponding acceptable bands of variation (BV) to set the final statistical distribution (Tf). If a first tentative statistical distribution (Tj) yields one histogram (e.g., cDVH) enclosed within the corresponding acceptable band of variation (BV) but another histogram (e.g., cDRVH) extending beyond the acceptable band of variation (BV), a new tentative statistical distribution (T(j+1)) is selected and the method is carried out again, until a tentative statistical distribution (T(j+k)) is found that fits all the required histograms into the corresponding acceptable bands of variation (BV). Dose rate related histograms (e.g., DRVH, DDRH) are particularly useful for treatment plans comprising beamlets to be emitted in FLASH-mode for depositing doses at ultra-high deposition rates into selected spots of the structure of interest.

TABLE-US-00001 Ref Description BV Band of variation CDj Calculated distribution of cDVHj and cDVRHj D(R)DH Dose (rate) distribution histogram (D(R)DH) includes both “dose distribution histograms (DDH)” and “dose rate distribution histograms (DRDH).” cD(R)Dj Calculated dose (rate) distribution histogram of beamlet Si with tentative statistical distribution Tj cDVH Calculated dose volume histogram of beamlet Si with tentative statistical distribution Tj cDRVH Calculated dose rate volume histogram of beamlet Si with tentative statistical distribution Tj cD(R)Dj0 Minimum calculated dose (rate) distribution histogram cD(R)Dj1 Maximum calculated dose (rate) distribution histogram CLj Confidence level CLf Final confidence level Dj Beamlet size DRj Dose deposition rate x x = 1 − N Runs j, k Identifies one tentative statistical distribution of an operational parameter MUai Actual monitor unit of a beamlet Si as delivered by the particle accelerator MUj Monitor unit distribution of a beamlet Si according to the tentative statistical distribution MUij Random value from the monitor unit distribution of a beamlet Si MUpi Planned monitor unit of a beamlet Si PBS Pencil Beam Scanning Si Beamlet Tj, T(j + k) Tentative statistical distribution of an operational parameter Tf Final statistical distribution of an operational parameter TP Treatment plan t0aj Actual starting time of a bemlet Si t0j Starting time distribution of a beamlet Si according to the tentative statistical distribution t0ij Random value from the starting time distribution of a beamlet Si t0pi Planned starting time of a beamlet Si t1ai Actual end time of a bemlet Si t1i End time distribution of a beamlet Si according to the tentative statistical distribution t1ij Random value from the end time distribution of a beamlet Si t1pi Planned end time of a beamlet Si tD(R)DH Target dose (rate) distribution histogram of beamlet Si with tentative statistical distribution Tj tDVH Target dose volume histogram of beamlet Si with tentative statistical distribution Tj tDRVH Target dose rate volume histogram of beamlet Si with tentative statistical distribution Tj Xai Actual position of a beamlet Si Xi Position of a beamlet Si Xij Random value from the position of a beamlet Si Xpi Planned position of a beamlet Si (X, Y) Plane normal to the irradiation axis (Z) Z Irradiation axis μj Average value of an operational parameter distribution of beamlet Si σj Standard deviation of an operational parameter distribution of beamlet Si