CN109661860A - Particle accelerator and its determination method, apparatus and medium for drawing particle energy - Google Patents
Particle accelerator and its determination method, apparatus and medium for drawing particle energy Download PDFInfo
- Publication number
- CN109661860A CN109661860A CN201880002494.3A CN201880002494A CN109661860A CN 109661860 A CN109661860 A CN 109661860A CN 201880002494 A CN201880002494 A CN 201880002494A CN 109661860 A CN109661860 A CN 109661860A
- Authority
- CN
- China
- Prior art keywords
- particle
- pulse
- phase
- accelerator
- energy
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- H—ELECTRICITY
- H05—ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
- H05H—PLASMA TECHNIQUE; PRODUCTION OF ACCELERATED ELECTRICALLY-CHARGED PARTICLES OR OF NEUTRONS; PRODUCTION OR ACCELERATION OF NEUTRAL MOLECULAR OR ATOMIC BEAMS
- H05H7/00—Details of devices of the types covered by groups H05H9/00, H05H11/00, H05H13/00
- H05H7/001—Arrangements for beam delivery or irradiation
-
- H—ELECTRICITY
- H05—ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
- H05H—PLASMA TECHNIQUE; PRODUCTION OF ACCELERATED ELECTRICALLY-CHARGED PARTICLES OR OF NEUTRONS; PRODUCTION OR ACCELERATION OF NEUTRAL MOLECULAR OR ATOMIC BEAMS
- H05H7/00—Details of devices of the types covered by groups H05H9/00, H05H11/00, H05H13/00
- H05H7/001—Arrangements for beam delivery or irradiation
- H05H2007/008—Arrangements for beam delivery or irradiation for measuring beam parameters
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Plasma & Fusion (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Particle Accelerators (AREA)
Abstract
It discloses particle accelerator and its draws the determination method, apparatus and medium of particle energy.The particle accelerator by FCT have particle by when the pulse signal that generates, sampling pulse signal is formed after ADC is sampled, and when reaching reference pulse number by the step-by-step counting number of pulse counter sampling pulse signal, it obtains corresponding with reference to duration, and is sent to controller by reference pulse number and with reference to duration；Controller by reference to umber of pulse, with reference to duration and pulse movement length, determine movement velocity of the particle in accelerating cavity, and combine grain proton rest mass, determine kinetic energy of the particle in accelerating cavity as particle energy.
Description
Technical field
The invention relates to Medical Devices technologies, such as are related to a kind of particle accelerator and its draw particle energy
Determine method, apparatus and medium.
Background technique
In biomedicine, when carrying out radiotherapy to tumour, need to carry out charged particle using particle accelerator
Accelerate so that particle obtains energy.Since the energy of particle determines the depth that radiotherapy energy is deposited in human body, passing through
It before synchrotron draws particle progress tumor radiotherapy treatment, needs accurately to measure particle energy, so that accurately positioning is put
Depth is treated, farthest to reduce damage of the radiotherapy side effect to human health cell.
The method measured indirectly is generallyd use to the measurement of particle energy in the related technology: by measurement synchrotron
The magnetic field strength of main second level magnet is counter to push away particle energy, or the resonance frequency by the accelerating cavity in measurement synchrotron
It is counter to push away particle energy.However, in indirect measurement process, not due to the tuned frequency of main secondary magnetic field intensity or accelerating cavity
Stability, so that the accuracy of the anti-particle energy pushed away is lower, the energy for the particle for causing particle accelerator to be drawn and patient are real
The energy match of particle needed for border is poor, increases the risk of patient's treatment.
Summary of the invention
The application provides a kind of particle accelerator and its draws the determination method, apparatus and medium of particle energy, to improve
Determine that accelerator draws the accuracy of particle energy.
The embodiment of the present application provides a kind of determination method of accelerator extraction particle energy, comprising:
Obtain in pulse counter the pulse number of received sampling pulse signal reach corresponding when reference pulse number M
With reference to duration；Wherein, described to sum to obtain according to the pulse period of M pulse with reference to duration；
According to described with reference to duration, the reference pulse number and pulse movement length, determine particle in accelerating cavity
Movement velocity；Wherein, the pulse movement length is movement length of the particle under a pulse period；
According to the movement velocity and rest mass of the particle, determine that kinetic energy of the particle in the accelerating cavity is made
For particle energy.
The embodiment of the present application also provides the determining devices that a kind of accelerator draws particle energy, comprising:
With reference to duration obtain module, be set as obtain pulse counter in received sampling pulse signal pulse number
Reach corresponding with reference to duration when reference pulse number；Wherein, described all according to the pulse of the several pulses of reference pulse with reference to duration
Expect and obtains；
Movement velocity determining module is set as according to described with reference to duration, the reference pulse number and pulse movement
Length determines movement velocity of the particle in accelerating cavity；Wherein, the pulse movement length is the particle in a pulse
Movement length under period；
Particle energy determining module is set as movement velocity and rest mass according to the particle, determines the grain
Kinetic energy of the son in the accelerating cavity is as particle energy.
The embodiment of the present application also provides a kind of particle accelerator, including injector, beam device and the synchronous ring of accelerator are cut,
The injector is connected with described one end for cutting beam device, and the other end for cutting beam device ring synchronous with the accelerator is connected, also
Including fast current mutual inductor, analog-digital converter, pulse counter and controller；
The fast current mutual inductor is installed on the synchronous ring of the accelerator, is set as when there is particle to pass through, generates
Pulse signal, and the pulse signal is sent to the analog-digital converter；
The analog-digital converter is connected between the fast current mutual inductor and the pulse counter, be set as by
The pulse signal is sampled using frequency according to setting, obtains sampling pulse signal, and the sampling pulse signal is sent out
It send to the pulse counter；
The pulse counter is connected between the fast current mutual inductor and the controller, is set as to default
In reference time the pulse number of the received sampling pulse signal counted, using count results as reference pulse
Number, and the reference pulse number is sent to the controller；
The controller is connect with the pulse counter, is set as according to the preset reference time, received described
Reference pulse number and pulse movement length determine movement velocity of the particle in accelerating cavity；According to the movement of the particle
Speed and rest mass determine kinetic energy of the particle in the accelerating cavity as particle energy；Wherein, the pulse
Movement length is movement length of the particle under a pulse period.
The embodiment of the present application also provides a kind of computer readable storage mediums, are stored thereon with computer program, the journey
A kind of determination method of accelerator extraction particle energy as provided by first aspect embodiment is realized when sequence is executed by processor.
Detailed description of the invention
Figure 1A is the structure chart of one of the embodiment of the present application one particle accelerator；
Figure 1B is one of the embodiment of the present application one pulse signal schematic diagram；
Fig. 2 is the flow chart for the determination method that one of the embodiment of the present application two accelerator draws particle energy；
Fig. 3 is the flow chart for the determination method that one of the embodiment of the present application three accelerator draws particle energy；
Fig. 4 is the structure chart for the determining device that one of the embodiment of the present application four accelerator draws particle energy.
Specific embodiment
The application is described in further detail with reference to the accompanying drawings and examples.It is understood that this place is retouched
The specific embodiment stated is used only for explaining the application, rather than the restriction to the application.It also should be noted that in order to just
Part relevant to the application is illustrated only in description, attached drawing rather than entire infrastructure.
Embodiment one
Figure 1A is the structure chart of one of the embodiment of the present application one particle accelerator.Particle accelerator as shown in Figure 1A
Include: injector 110, cut the synchronous ring 130 of beam device 120 and accelerator, the injector 110 cuts the one of beam device 120 with described
End is connected, and the other end for cutting beam device 120 ring 130 synchronous with the accelerator is connected；The particle accelerator further includes quick
Current transformer (Fast Current Transformer, FCT) 140, analog-digital converter (Analog to Digital
Converter, ADC) 150, pulse counter 160 and controller 170.
The fast current mutual inductor 140 is installed on the synchronous ring 130 of the accelerator, is set as having particle to pass through
When, pulse signal is generated, and the pulse signal is sent to the analog-digital converter 150.
The analog-digital converter 150 is connected between the fast current mutual inductor 140 and the pulse counter 160,
It is set as sampling the pulse signal according to setting using frequency, obtains sampling pulse signal, and by the sampling arteries and veins
It rushes signal and is sent to the pulse counter 160.
The pulse counter 160 is connected between the fast current mutual inductor 140 and the controller 170, setting
When reaching reference pulse number M for the pulse number in sampling pulse signal, when by the pulse period summation of M pulse as reference
It is long, and the reference pulse number and the duration that refers to are sent to the controller 170.
The controller 170 is connect with the pulse counter 160, is set as referring to duration, the reference according to described
Umber of pulse and pulse movement length determine movement velocity of the particle in accelerating cavity；According to the movement velocity of the particle
And rest mass, determine kinetic energy of the particle in the accelerating cavity as particle energy；Wherein, the pulse movement
Length is movement length of the particle under a pulse period.
In one embodiment, injector 110, the element vapour particles for being set as needs to inject are ionized into ion, and will
Ion is formed by particle beam and is injected into the synchronous ring 130 of accelerator, wherein injector 110 determines the particle beams to be injected
Type and beam intensity etc..Beam device 120 is cut, is set as limiting the particle that injector 110 is injected into the synchronous ring 130 of accelerator
The line length of line；Accelerator synchronizes ring 130, is set as to the particle sync plus white injected.Wherein, particle can be matter
Son or heavy ion.
In one embodiment, the synchronous ring 130 of accelerator is hollow shell structure, wherein hollow space forms accelerating cavity.
In the shell structure of the synchronous ring 130 of accelerator, FCT140 is installed.In one embodiment, particle is synchronous in accelerating cavity adds
Speed, and when there is particle to pass through FCT140, FCT140, which can be corresponded to, generates a pulse signal, and generated pulse signal is sent out
It send to ADC150.Figure 1B representatively illustrates a kind of pulse signal figure, wherein a pulse period corresponds to particle in accelerating cavity
The time span of movement one week.
In one embodiment, when ADC150 receives pulse signal, it will according to setting sample frequency to pulse signal
It is sampled, obtains sampling pulse signal, and sampling pulse is sent to the pulse counter 160.Wherein, setting sampling frequency
Rate by particle accelerator operator or designer according to demand or empirical value sets itself.In one embodiment, setting is adopted
Sample frequency can be 100MHz.
In one embodiment, after pulse counter 160 receives sampling pulse signal, sampling pulse signal is counted
Number, and when number of pulses reaches reference pulse number, the pulse period of multiple pulses in reference pulse number is summed, is referred to
Duration.Wherein, reference pulse number is the pulse number for having the complete pulse period.Figure 1B representatively illustrates the side of step-by-step counting
Formula is the case where rising edge of a pulse counts.Wherein, reference pulse number is M, and the rising edge of a pulse detected is that (M+1) is a, phase
When the reference answered a length of (T2-T1), the sum of corresponding M pulse period, or the sum of (M+1) a pulse period.
Pulse counter 160 is sent to controller 170 by reference pulse number and with reference to duration.In one embodiment, pulse
Counter can be field programmable gate array (Field Programmable Gate Array, FPGA).Wherein, ginseng is preset
The time range that the time is preset step-by-step counting is examined, by particle accelerator operator or designer according to precision need
It wants or empirical value sets itself.In one embodiment, reference pulse number can be 500.
In one embodiment, controller 170 is according to preset reference time, reference pulse number and pulse movement length,
Determine movement velocity of the particle in accelerating cavity；And according to the movement velocity and rest mass of particle, determine that particle is accelerating
Kinetic energy in chamber is as particle energy.Wherein, preset reference time and pulse movement length preset and are stored in control
In the storage device of device.Wherein, pulse movement length is movement length of the particle under a pulse period.In an embodiment
In, pulse movement length is the perimeter of accelerating cavity in the synchronous ring 130 of accelerator.
The embodiment of the present application by FCT have particle by when the pulse signal that generates, sampling arteries and veins is formed after ADC is sampled
When rushing signal, and reaching reference pulse number by the step-by-step counting number of pulse counter sampling pulse signal, obtain corresponding
Controller is sent to reference to duration, and by reference pulse number and with reference to duration；Controller by reference to umber of pulse, with reference to duration with
And pulse movement length, it determines movement velocity of the particle in accelerating cavity, and combine grain proton rest mass, determines that particle exists
Kinetic energy in accelerating cavity is as particle energy.Particle energy is directly determined by the movement velocity of particle by adopting the above technical scheme
Amount is improved determining accelerator and draws the accuracy of particle energy, and then avoided and put due to accelerator particle energy with patient
Actual bodily harm caused by particle energy needed for treating mismatches, reduces the Operative risk of patient.
In one embodiment, particle accelerator further include: detection probe 180 and synchronous ring controller 190.
The controller 170 is connected with the fast current mutual inductor 140, and the synchronous ring controller 190 is connected to institute
It states between the ring 130 synchronous with the accelerator of controller 170.
The detection probe 180 is installed on the accelerator synchronization ring 130, is set as detecting the accelerator synchronization
The voltage-phase of ring 130 is sent to the controller 170 as accelerator phase, and by the accelerator phase.
The fast current mutual inductor 140, is also configured to when there is particle to pass through, and obtains the beam phase of particle, and will
The beam phase is sent to the controller 170.
The controller 170 is set as the particle energy being compared with default particle energy, and in the particle
When energy is less than the default particle energy, first phase adjustment letter is generated according to the beam phase and the accelerator phase
Number it is sent to the Accelerator control device 190；When the particle energy is greater than the default particle energy, according to the line
Phase and the accelerator phase generate second phase adjustment signal and are sent to the Accelerator control device 190.
The Accelerator control device 190 is set as adjusting the accelerator phase when receiving first phase adjustment signal
Position, so that the beam phase and the phase difference value of the accelerator phase are accelerating phase threshold value；Receiving second phase tune
When entire signal, the accelerator phase is adjusted, so that the beam phase and the phase difference value of the accelerator phase are to slow down
Phase threshold.
Wherein, the decelerating phase threshold value differs 180 ° with the accelerating phase threshold value.
In one embodiment, detection probe 180 is in real time or timing detects the voltage-phase that accelerator synchronizes ring 130, and will
The voltage-phase detected is sent in controller 170 as accelerator phase.FCT140 can also be detected when there is particle to pass through
The beam phase of particle, and beam phase detected is sent in controller 170.It has determined and has accelerated in controller 170
It, can be by determining particle energy and the default grain determined according to patient's lesion demand in chamber after the energy of the particle of sync plus white
Sub- energy compares.
When particle energy is less than default particle energy, show that the energy of the particle in the synchronous ring 130 of accelerator is too small, by
Required energy value when particle energy is unable to reach Patients Treated by Radiotherapy, when carrying out radiotherapy to patient using the particle of present energy,
It will release energy when particle will be made to be unable to reach required emplacement depth, therefore controller 170 can correspond to and generate the first phase
Position adjustment signal is sent in synchronous ring controller 190.Synchronous ring controller 190 adds according to first phase adjustment signal, adjustment
Fast device phase makes the phase difference value accelerating phase threshold value of beam phase and accelerator phase.Wherein, accelerating phase threshold value is by grain
The designer of sub- accelerator is according to test result sets itself.
When particle energy is greater than default particle energy, show that the energy of the particle in the synchronous ring 130 of accelerator is excessive, by
In particle energy beyond energy value needed for Patients Treated by Radiotherapy, when carrying out radiotherapy to patient using the particle of energy before current year, it will
Make when particle releases energy beyond required emplacement depth, therefore controller 170 can correspond to and generate second phase adjustment signal
It is sent in synchronous ring controller 190.Synchronous ring 190 adjusts accelerator phase, makes line phase according to second phase adjustment signal
The phase difference value of position and accelerator phase is decelerating phase threshold value.Wherein, decelerating phase threshold value is differed with accelerating phase threshold value
180°。
When particle energy is equal with default particle energy, show the energy of the particle in the synchronous ring 130 of accelerator just
For energy value needed for Patients Treated by Radiotherapy, therefore controller 170 will not generate synchronous ring controller 190 is made to adjust accelerator phase
Phase adjustment signal.
The embodiment of the present application detects the voltage-phase of the synchronous ring of accelerator as accelerator phase by setting detection probe
It is sent to controller；Beam phase when particle accelerates is obtained by FCT is sent to controller；By controller by determining grain
Sub- energy is compared with default particle energy, generates the first phase adjustment letter for making synchronous ring controller adjustment accelerator phase
Number or second phase adjustment signal；Synchronize ring controller first phase adjustment signal or second phase adjustment signal based on the received
The voltage-phase for adjusting the synchronous ring of accelerator, makes beam phase and accelerator phase difference value accelerating phase threshold value or decelerating phase
Threshold value.Accelerator phase is adjusted in particle energy and inconsistent default particle energy by adopting the above technical scheme, with
Realize that the energy to the particle of sync plus white in the synchronous ring of accelerator control effectively.
Embodiment two
Fig. 2 is the flow chart for the determination method that one of the embodiment of the present application two accelerator draws particle energy.This Shen
Please embodiment when being suitable for accelerating the synchronized ring of particle of the synchronous ring of particle booster, the case where determining particle energy, the party
Method is suited to speed up device and draws the determining device of particle energy to execute, and the device is real by least one of software and hardware
It is existing, and in the controller of concrete configuration particle accelerator provided by the technical solution such as above-mentioned multiple embodiments.
Accelerator as shown in Figure 2 draws the determination method of particle energy, comprising: step S210, step S220 and step
Rapid S210.
In step S210, obtain in pulse counter the pulse number of received sampling pulse signal reach with reference to arteries and veins
It rushes corresponding with reference to duration when number M.
Wherein, described to sum to obtain according to the pulse period of M pulse with reference to duration.
It, can be according to the step-by-step counting mode of setting, to sampling pulse after pulse counter receives sampling pulse signal
Pulse in signal is counted.When number of pulses reaches reference pulse number M, the sum of the pulse period of M pulse is calculated, is obtained
To referring to duration.Pulse counter will be sent in controller with reference to duration and reference pulse number.Certainly, reference pulse number may be used also
To be stored in advance in the controller, pulse counter will be only sent in controller with reference to duration, be deposited in advance in controller at this time
The reference pulse number of storage is consistent with the reference pulse number in pulse counter.In one embodiment, pulse counter can be
FPGA.In one embodiment, reference pulse number can be 500.
In one embodiment, reference pulse number is the pulse number for having the complete pulse period.In one embodiment, when adopting
When being counted with rising edge, when reference pulse number is M, the rising edge of a pulse detected is that (M+1) is a.
In step S220, according to described with reference to duration, the reference pulse number and pulse movement length, grain is determined
Movement velocity of the son in accelerating cavity.
Wherein, the pulse movement length is movement length of the particle under a pulse period.
According to pulse movement length and reference pulse number, determine with reference to movement overall length of the particle in accelerator in duration
Degree；According to movement total length and with reference to the ratio between duration, movement velocity of the particle in accelerating cavity is determined.
In one embodiment, movement velocity of the particle in accelerating cavity is determined according to the following formula:
Wherein, v is the movement velocity of the particle, l0For the pulse movement length, n is the reference
Umber of pulse, t are described with reference to duration.
In one embodiment, reference pulse number 500 are reached to the step-by-step counting of sampling pulse signal in pulse counter,
A length of 50.01 μ s when obtained corresponding reference, a pulse period corresponding pulse movement length are the synchronous ring of accelerator
Perimeter 18m, then, corresponding movement velocity are as follows: (18 × 500)/(50.01 × 10-6)=1.7996 × 108m/s。
In step S230, according to the movement velocity and rest mass of the particle, determine that the particle adds described
Kinetic energy in fast chamber is as particle energy.
According to the movement velocity and rest mass of particle, the moving-mass of particle is determined；According to the moving-mass and
The rest mass determines kinetic energy of the particle in the accelerating cavity as particle energy.
In one embodiment, according to formulaDetermine the moving-mass of particle；
According to formula E=mc2-m0c2, determine the particle energy；
Wherein, v is the movement velocity, m0For the rest mass, c is the light velocity, and E is the particle energy.
In one embodiment, when particle is proton, matter proton rest mass is 938MeV/c2, and the movement speed of proton
Degree is 1.7996 × 108When m/s, corresponding particle energy are as follows: { 1/sqrt [1- (1.7996 × 108)2/(3×108)2]–1}×
938=249.84MeV.
The embodiment of the present application by obtain in pulse counter the pulse number of received sampling pulse signal reach ginseng
It examines corresponding with reference to duration when umber of pulse；According to reference duration, reference pulse number and pulse movement length, determine that particle exists
Movement velocity in accelerating cavity；According to the movement velocity and rest mass of particle, determine that kinetic energy of the particle in accelerating cavity is made
For particle energy.Particle energy is directly determined by the movement velocity of particle by adopting the above technical scheme, improves determining acceleration
Device draws the accuracy of particle energy, and then does not avoid the particle energy as needed for accelerator particle energy and Patients Treated by Radiotherapy not
With caused actual bodily harm, the Operative risk of patient is reduced.
Embodiment three
Fig. 3 is the flow chart for the determination method that one of the embodiment of the present application three accelerator draws particle energy.
" according to the preset reference time, the reference pulse number and pulse movement length, particle is determined in operation
It after movement velocity in accelerating cavity ", further include " according to the sample frequency of the sampling pulse signal, determining and referring to duration
Accuracy value；According to described with reference to duration accuracy value and in conjunction with the reference pulse number and pulse movement length, determine described in
The movement velocity section of particle；According to the movement velocity section and the grain proton rest mass, the particle is determined
Particle energy section ", in allowable range of error, to determine section where particle energy.
In one embodiment, " according to the movement velocity and rest mass of the particle, determine that the particle exists in operation
Kinetic energy in the accelerating cavity is as particle energy " after, addition " is compared the particle energy and default particle energy
Compared with；According to the beam phase of Particles Moving in comparison result and accelerating cavity, accelerator phase is adjusted；Wherein, the accelerator
Phase is the voltage-phase of the synchronous ring of accelerator ", to realize the energy to the particle in accelerating cavity by determining particle energy
It adjusts.
Accelerator as shown in Figure 3 draws the determination method of particle energy, comprising: step S310 to step S370.
In step s310, obtain in pulse counter the pulse number of received sampling pulse signal reach with reference to arteries and veins
It rushes corresponding with reference to duration when number M.
Wherein, described to sum to obtain according to the pulse period of M pulse with reference to duration.
In step s 320, grain is determined with reference to duration, the reference pulse number and pulse movement length according to described
Movement velocity of the son in accelerating cavity.
Wherein, the pulse movement length is movement length of the particle under a pulse period.
In step S330, according to the sample frequency of the sampling pulse signal, determines and refer to duration accuracy value；According to institute
It states with reference to duration accuracy value and in conjunction with the reference pulse number and pulse movement length, determines the movement velocity of the particle
Section.
In one embodiment, by the integral multiple reciprocal of the sample frequency of sampling pulse signal, as reference duration precision
Value；The reference duration that will acquire, which is added, obtains maximum with reference to duration with reference to duration accuracy value；The reference duration that will acquire subtracts ginseng
It examines duration accuracy value and obtains minimum with reference to duration；Duration, reference pulse number and pulse movement length are referred to according to minimum, really
Determine the maximum movement speed of particle；Duration, reference pulse number and pulse movement length are referred to according to maximum, determines particle
Minimum movement speed；Using minimum movement speed and maximum movement speed as the endpoint value in section, continuous movement velocity is obtained
Section.Wherein, movement velocity section can be open interval, or closed interval.
In one embodiment, according to formulaWherein, vminFor minimum movement speed, vmaxFor maximum
Movement velocity, l0For the pulse movement length, n is the reference pulse number, and t is described with reference to duration, and f is sampling frequency
Rate, a are the error parameter allowed, and wherein a is positive integer.Preferably, 1 a.
In one embodiment, corresponding to refer to duration precision when error parameter is 1 when sample frequency is 100MHz
Value is 0.01 μ s.When reference pulse number is 500, when corresponding reference a length of 50.01 μ s, when pulse movement length is 18m,
Minimum refers to duration are as follows: 50.01-0.01=50 μ s；Maximum refers to duration are as follows: 50.01+0.01=50.02 μ s.Correspondingly, most
Big movement velocity are as follows: 18 × 500/ (50 × 10-6)=1.8 × 108m/s；Minimum movement speed are as follows: 18 × 500/ (50.02 ×
10-6)=1.79928 × 108m/s；Correspondingly, speed movement section is (1.79928 × 108M/s, 1.8 × 108M/s), or
It is [1.79928 × 108M/s, 1.8 × 108m/s]。
In step S340, according to the movement velocity section and the grain proton rest mass, the particle is determined
Particle energy section.
According to the minimum movement speed and grain proton rest mass in movement velocity section, the minimum of the particle is determined
Moving-mass；Determine kinetic energy of the particle in the accelerating cavity as most according to minimum movement quality and the rest mass
Small particles energy.
In one embodiment, according to formulaDetermine the minimum movement quality of particle；According to public affairs
Formula Emin=mminc2-m0c2, determine smallest particles energy；Wherein, vminFor minimum movement speed, m0For the rest mass, c is
The light velocity, EminFor smallest particles energy.
According to the maximum movement speed and grain proton rest mass in movement velocity section, the maximum of the particle is determined
Moving-mass；Determine kinetic energy of the particle in the accelerating cavity as most according to largest motion quality and the rest mass
Big particle energy.
In one embodiment, according to formulaDetermine the largest motion quality of particle；According to public affairs
Formula Emax=mmaxc2-m0c2, determine maximum particle energy；Wherein, vmaxFor maximum movement speed, m0For the rest mass, c is
The light velocity, EmaxFor maximum particle energy.
Using smallest particles energy and maximum particle energy as the endpoint value in section, continuous particle energy section is obtained.
With corresponding movement velocity section it is found that particle energy section can be open interval (Emin, Emax), or closed interval
[Emin, Emax]。
In one embodiment, when speed movement section is (1.79928 × 108M/s, 1.8 × 108M/s), Huo Zhewei
[1.79928×108M/s, 1.8 × 108M/s], maximum movement speed is 1.8 × 108M/s, minimum movement speed are 1.79928
×108m/s；The maximum particle energy of corresponding particle is 234.5Mev, and smallest particles energy is 234.24Mev；Particle energy
Section is (234.24Mev, 234.5Mev), or is [234.24Mev, 234.5Mev].
In step S350, according to the movement velocity and rest mass of the particle, determine that the particle adds described
Kinetic energy in fast chamber is as particle energy.
In step S360, the particle energy is compared with default particle energy.
Wherein, it presets particle energy to be calculated according to the lesion situation of different patients correspondence, corresponding Patients Treated by Radiotherapy when institute
The theoretical value of the particle energy needed.
In step S370, according to the beam phase of Particles Moving in comparison result and accelerating cavity, accelerator phase is adjusted
Position；Wherein, the accelerator phase is the voltage-phase of the synchronous ring of accelerator.
In one embodiment, if the particle energy is less than the default particle energy, first phase adjustment letter is generated
Number, so that the beam phase and the phase difference value of the accelerator phase are accelerating phase threshold value；If the particle energy is big
In the default particle energy, then second phase adjustment signal is generated, so that the beam phase and the accelerator phase
Phase difference value is decelerating phase threshold value；Wherein, the decelerating phase threshold value differs 180 ° with the accelerating phase threshold value.
Wherein, first phase adjustment signal is sent to synchronous ring controller by controller, and synchronous ring controller receives this
After first phase adjustment signal, corresponding adjustment accelerator phase accelerates beam phase and the phase difference value of accelerator phase
Phase difference value；Second phase adjustment signal is sent to synchronous ring controller by controller, synchronous ring controller receive this second
After phase adjustment signal, corresponding adjustment accelerator phase makes the phase difference value decelerating phase of beam phase Yu accelerator phase
Difference.Wherein, accelerating phase threshold value is by the designer of particle accelerator according to test result sets itself.
The embodiment of the present application makes by the way that after determining the movement velocity of particle, additional particle energy section determines step
It obtains and is determined according to the sample frequency of sampling pulse signal with reference to duration accuracy value；And it is combined according to reference duration accuracy value and refers to arteries and veins
Number and pulse movement length are rushed, determines the movement velocity section of particle；According to the static of movement velocity section and particle
Quality determines the particle energy section of particle.So that determining section where particle energy in allowable range of error.The application
Embodiment is by after determining particle energy, operation of the additional particle energy compared with default particle energy, and by comparing
As a result the voltage-phase of the synchronous ring of adjustment accelerator, and then realize the energy adjustment to the particle accelerated in accelerating cavity.
Example IV
Fig. 4 is the structure chart for the determining device that one of the embodiment of the present application four accelerator draws particle energy.This Shen
Please embodiment when being suitable for accelerating the synchronized ring of particle of the synchronous ring of particle booster, the case where determining particle energy, the dress
It sets and is realized by least one of software and hardware, and be configured at the particle as provided by the technical solution of above-mentioned multiple embodiments
In the controller of accelerator.Accelerator as shown in Figure 4 draws the determining device of particle energy, comprising: obtains mould with reference to duration
Block 410, movement velocity determining module 420 and particle energy determining module 430.
Wherein, obtain module 410 with reference to duration, be set as obtaining in pulse counter received sampling pulse signal
Pulse number reaches corresponding with reference to duration when reference pulse number M；Wherein, described all according to the pulse of M pulse with reference to duration
Expect and obtains.
Movement velocity determining module 420 is set as according to described with reference to duration, the reference pulse number and pulse fortune
Dynamic length, determines movement velocity of the particle in accelerating cavity；Wherein, the pulse movement length is the particle in an arteries and veins
Rush the movement length under the period.
Particle energy determining module 430, is set as movement velocity and rest mass according to the particle, determine described in
Kinetic energy of the particle in the accelerating cavity is as particle energy.
The embodiment of the present application obtains module by reference to duration and obtains the received sampling pulse signal of institute in pulse counter
Pulse number reach corresponding with reference to duration when reference pulse number；By movement velocity determining module according to reference duration, ginseng
Umber of pulse and pulse movement length are examined, determines movement velocity of the particle in accelerating cavity；Pass through particle energy determining module
According to the movement velocity and rest mass of particle, determine kinetic energy of the particle in accelerating cavity as particle energy.Using above-mentioned
Technical solution directly determines particle energy by the movement velocity of particle, improves determining accelerator and draws the accurate of particle energy
Degree, and then avoid actual bodily harm caused by the particle energy as needed for accelerator particle energy and Patients Treated by Radiotherapy mismatches, drop
The low Operative risk of patient.
In one embodiment, the movement velocity determining module 420, comprising: movement velocity determination unit.
Movement velocity determination unit is set as movement velocity of the determining particle in accelerating cavity according to the following formula:
Wherein, v is the movement velocity of the particle, l0For the pulse movement length, n is the reference
Umber of pulse, t are described with reference to duration.
In one embodiment, described device further includes particle energy section determining module, comprising: movement velocity section is true
Order member and particle energy interval determination unit.
Movement velocity interval determination unit, be set as it is described according to it is described with reference to duration, the reference pulse number and
Pulse movement length, after determining movement velocity of the particle in accelerating cavity, according to the sampling of sampling pulse signal frequency
Rate determines and refers to duration accuracy value；Pulse fortune is recalled according to the reference duration accuracy value and in conjunction with the reference pulse number scale
Dynamic length, determines the movement velocity section of the particle；
Particle energy interval determination unit is set as the static matter according to the movement velocity section and the particle
Amount, determines the particle energy section of the particle.
In one embodiment, the particle energy determining module 430, comprising: moving-mass determination unit and particle energy
Determination unit.
Moving-mass determination unit is set as movement velocity and rest mass according to the particle, determines the grain
The moving-mass of son.
Particle energy determination unit is set as determining the particle in institute according to the moving-mass and the rest mass
The kinetic energy in accelerating cavity is stated as particle energy.
In one embodiment, described device further includes accelerator phase adjusting module, comprising: energy comparison unit and phase
Position adjustment unit.
Energy comparison unit, is set as in the movement velocity and rest mass according to the particle, determine described in
Particle compares the particle energy and default particle energy after the kinetic energy in the accelerating cavity is as particle energy
Compared with.
Phasing unit, is set as the beam phase according to Particles Moving in comparison result and accelerating cavity, and adjustment adds
Fast device phase；Wherein, the accelerator phase is the voltage-phase of the synchronous ring of accelerator.
In one embodiment, the phasing unit, setting are as follows: if the particle energy is less than the default particle energy
Amount, then generate first phase adjustment signal, so that the beam phase and the phase difference value of the accelerator phase are to accelerate phase
Position threshold value；If the particle energy is greater than the default particle energy, second phase adjustment signal is generated, so that the line
The phase difference value of phase and the accelerator phase is decelerating phase threshold value；Wherein, the decelerating phase threshold value and the acceleration
Phase threshold differs 180 °.
Above-mentioned accelerator draws particle energy, and accelerator provided by the application any embodiment can be performed in fixed output quota product really
The determination method for drawing particle energy has the corresponding functional module of determination method for executing accelerator extraction particle energy and has
Beneficial effect.
Claims (10)
1. a kind of particle accelerator, including injector, cut beam device and the synchronous ring of accelerator, the injector and described cuts beam device
One end be connected, it further includes fast current mutual inductor, modulus that the other end for cutting beam device ring synchronous with the accelerator is connected
Converter, pulse counter and controller；
The fast current mutual inductor is installed on the synchronous ring of the accelerator, is set as when there is particle to pass through, generates pulse
Signal, and the pulse signal is sent to the analog-digital converter；
The analog-digital converter is connected between the fast current mutual inductor and the pulse counter, is set as according to setting
Surely the pulse signal is sampled using frequency, obtains sampling pulse signal, and the sampling pulse signal is sent to
The pulse counter；
The pulse counter is connected between the fast current mutual inductor and the controller, is set as in sampling pulse
When the pulse number of signal reaches reference pulse number M, the pulse period summation of the M pulse is used as and refers to duration, and will
The reference pulse number and the duration that refers to are sent to the controller；
The controller is connect with the pulse counter, be set as according to it is described with reference to duration, the reference pulse number and
Pulse movement length determines movement velocity of the particle in accelerating cavity；According to the movement velocity of the particle and static matter
Amount, determines kinetic energy of the particle in the accelerating cavity as particle energy；Wherein, the pulse movement length is described
Movement length of the particle under a pulse period.
2. particle accelerator according to claim 1 further includes detection probe and synchronous ring controller；The controller with
The fast current mutual inductor be connected, the synchronous ring controller be connected to controller ring synchronous with the accelerator it
Between；
The detection probe is installed on the synchronous ring of the accelerator, is set as detecting the voltage phase of the synchronous ring of the accelerator
Position is used as accelerator phase, and the accelerator phase is sent to the controller；
The fast current mutual inductor, is also configured to when there is particle to pass through, and obtains the beam phase of particle, and by the line
Phase is sent to the controller；
The controller is set as the particle energy being compared with default particle energy, and small in the particle energy
When the default particle energy, first phase adjustment signal is generated according to the beam phase and the accelerator phase and is sent
To the Accelerator control device；When the particle energy is greater than the default particle energy, according to the beam phase and institute
It states accelerator phase generation second phase adjustment signal and is sent to the Accelerator control device；
The Accelerator control device is set as the accelerator phase being adjusted, so that institute when receiving first phase adjustment signal
The phase difference value for stating beam phase and the accelerator phase is accelerating phase threshold value；When receiving second phase adjustment signal,
The accelerator phase is adjusted, so that the beam phase and the phase difference value of the accelerator phase are decelerating phase threshold value；
Wherein, the decelerating phase threshold value differs 180 ° with the accelerating phase threshold value.
3. a kind of accelerator draws the determination method of particle energy, it is applied to particle accelerator as claimed in claim 1 or 2,
Include:
Obtain in pulse counter the pulse number of received sampling pulse signal reach corresponding reference when reference pulse number M
Duration；Wherein, described to sum to obtain according to the pulse period of M pulse with reference to duration；
According to described with reference to duration, the reference pulse number and pulse movement length, fortune of the particle in accelerating cavity is determined
Dynamic speed；Wherein, the pulse movement length is movement length of the particle under a pulse period；
According to the movement velocity and rest mass of the particle, determine kinetic energy of the particle in the accelerating cavity as grain
Sub- energy.
4. according to the method described in claim 3, wherein, it is described according to described with reference to duration, the reference pulse number and list
Impulse movement length determines movement velocity of the particle in accelerating cavity, comprising:
Movement velocity of the particle in accelerating cavity is determined according to the following formula:
Wherein, v is the movement velocity of the particle, l0For the pulse movement length, n is the reference pulse
Number, t are described with reference to duration.
5. according to the method described in claim 3, described according to the reference duration, the reference pulse number and pulse
Movement length, after determining movement velocity of the particle in accelerating cavity, further includes:
According to the sample frequency of the sampling pulse signal, determines and refer to duration accuracy value；Duration accuracy value is referred to according to described
And pulse movement length is recalled in conjunction with the reference pulse number scale, determine the movement velocity section of the particle；
According to the movement velocity section and the grain proton rest mass, the particle energy section of the particle is determined.
6. according to the method described in claim 3, wherein, the movement velocity and rest mass according to the particle, really
Fixed kinetic energy of the particle in the accelerating cavity is as particle energy, comprising:
According to the movement velocity and rest mass of the particle, the moving-mass of the particle is determined；
Determine kinetic energy of the particle in the accelerating cavity as particle energy according to the moving-mass and the rest mass
Amount.
7. according to the described in any item methods of claim 3-6, in the movement velocity according to the particle and static matter
Amount, determines the particle after the kinetic energy in the accelerating cavity is as particle energy, further includes:
The particle energy is compared with default particle energy；
According to the beam phase of Particles Moving in comparison result and accelerating cavity, accelerator phase is adjusted；Wherein, the accelerator
Phase is the voltage-phase of the synchronous ring of accelerator.
It is described according to comparison result and the beam phase 8. according to the method described in claim 7, wherein, described in adjustment
Accelerator phase, comprising:
If the particle energy is less than the default particle energy, first phase adjustment signal is generated, so that the line phase
The phase difference value of position and the accelerator phase is accelerating phase threshold value；
If the particle energy is greater than the default particle energy, second phase adjustment signal is generated, so that the line phase
The phase difference value of position and the accelerator phase is decelerating phase threshold value；
Wherein, the decelerating phase threshold value differs 180 ° with the accelerating phase threshold value.
9. the determining device that a kind of accelerator draws particle energy, comprising:
Obtain module with reference to duration, be set as obtaining in pulse counter the pulse number of received sampling pulse signal reach
It is corresponding when reference pulse number M to refer to duration；Wherein, described to sum to obtain according to the pulse period of M pulse of ginseng with reference to duration；
Movement velocity determining module is set as according to the reference duration, the reference pulse number and pulse movement length,
Determine movement velocity of the particle in accelerating cavity；Wherein, the pulse movement length is the particle in a pulse period
Under movement length；
Particle energy determining module is set as movement velocity and rest mass according to the particle, determines that the particle exists
Kinetic energy in the accelerating cavity is as particle energy.
10. a kind of computer readable storage medium, is stored thereon with computer program, realized such as when which is executed by processor
A kind of determination method of accelerator extraction particle energy described in any one of claim 3-8.
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/CN2018/117407 WO2020107146A1 (en) | 2018-11-26 | 2018-11-26 | Particle accelerator, method and apparatus for determining energy of particles led out therefrom, and medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109661860A true CN109661860A (en) | 2019-04-19 |
Family
ID=66118506
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201880002494.3A Pending CN109661860A (en) | 2018-11-26 | 2018-11-26 | Particle accelerator and its determination method, apparatus and medium for drawing particle energy |
| CN201911102646.XA Active CN110708856B (en) | 2018-11-26 | 2019-11-12 | Particle accelerator, determination method and device for energy of extracted particles of particle accelerator and medium |
Family Applications After (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911102646.XA Active CN110708856B (en) | 2018-11-26 | 2019-11-12 | Particle accelerator, determination method and device for energy of extracted particles of particle accelerator and medium |
Country Status (2)
| Country | Link |
|---|---|
| CN (2) | CN109661860A (en) |
| WO (1) | WO2020107146A1 (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109661860A (en) * | 2018-11-26 | 2019-04-19 | 新瑞阳光粒子医疗装备(无锡)有限公司 | Particle accelerator and its determination method, apparatus and medium for drawing particle energy |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102263003A (en) * | 2011-06-03 | 2011-11-30 | 中国科学院西安光学精密机械研究所 | Method and apparatus for mapping flight time and momentum energy of refraction type charged particle |
| US20120286702A1 (en) * | 2011-05-09 | 2012-11-15 | Bazaz Gaurav | Apparatus and method for energy storage with relativistic particle acceleration |
| CN102971802A (en) * | 2010-06-23 | 2013-03-13 | Gsi亥姆霍兹重离子研究中心有限责任公司 | Particle beam generating device |
| CN103412319A (en) * | 2013-08-19 | 2013-11-27 | 中国科学院高能物理研究所 | Detector signal processing method and device based on pulse width spectrum |
| CN103458967A (en) * | 2012-04-11 | 2013-12-18 | 株式会社东芝 | Radiation therapy system and therapy planning device |
| CN103605147A (en) * | 2013-11-22 | 2014-02-26 | 中国航空工业集团公司北京航空制造工程研究所 | Method and system for measuring energy density of multi-dimensional electronic beams based on marginal integration |
| WO2014123591A2 (en) * | 2012-10-17 | 2014-08-14 | Cornell University | Generation and acceleration of charged particles using compact devices and systems |
| WO2015028602A1 (en) * | 2013-08-30 | 2015-03-05 | Uniwersytet Jagiellonski | A method for calibration of tof-pet detectors using cosmic radiation |
| CN104663003A (en) * | 2012-07-27 | 2015-05-27 | 麻省理工学院 | Synchrocyclotron beam orbit and RF drive synchrocyclotron |
| CN106596112A (en) * | 2016-12-16 | 2017-04-26 | 哈尔滨工业大学 | Hall thruster ion energy distribution measuring method and system |
| CN108883304A (en) * | 2018-06-22 | 2018-11-23 | 新瑞阳光粒子医疗装备(无锡)有限公司 | Synchrotron control method, device, equipment and storage medium |
| CN110708856A (en) * | 2018-11-26 | 2020-01-17 | 新瑞阳光粒子医疗装备(无锡)有限公司 | Particle accelerator, determination method and device for energy of extracted particles of particle accelerator and medium |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2017160758A1 (en) * | 2016-03-17 | 2017-09-21 | Alexey Radovinsky | Particle acceleration in a variable-energy synchrocyclotron by a single-tuned variable-frequency rf drive |
| CN106855522B (en) * | 2016-12-30 | 2019-07-16 | 中国科学院高能物理研究所 | White light neutron imaging method and the material composition lossless detection method for using it |
-
2018
- 2018-11-26 CN CN201880002494.3A patent/CN109661860A/en active Pending
- 2018-11-26 WO PCT/CN2018/117407 patent/WO2020107146A1/en active Application Filing
-
2019
- 2019-11-12 CN CN201911102646.XA patent/CN110708856B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102971802A (en) * | 2010-06-23 | 2013-03-13 | Gsi亥姆霍兹重离子研究中心有限责任公司 | Particle beam generating device |
| US20120286702A1 (en) * | 2011-05-09 | 2012-11-15 | Bazaz Gaurav | Apparatus and method for energy storage with relativistic particle acceleration |
| CN102263003A (en) * | 2011-06-03 | 2011-11-30 | 中国科学院西安光学精密机械研究所 | Method and apparatus for mapping flight time and momentum energy of refraction type charged particle |
| CN103458967A (en) * | 2012-04-11 | 2013-12-18 | 株式会社东芝 | Radiation therapy system and therapy planning device |
| CN104663003A (en) * | 2012-07-27 | 2015-05-27 | 麻省理工学院 | Synchrocyclotron beam orbit and RF drive synchrocyclotron |
| WO2014123591A2 (en) * | 2012-10-17 | 2014-08-14 | Cornell University | Generation and acceleration of charged particles using compact devices and systems |
| CN103412319A (en) * | 2013-08-19 | 2013-11-27 | 中国科学院高能物理研究所 | Detector signal processing method and device based on pulse width spectrum |
| WO2015028602A1 (en) * | 2013-08-30 | 2015-03-05 | Uniwersytet Jagiellonski | A method for calibration of tof-pet detectors using cosmic radiation |
| CN103605147A (en) * | 2013-11-22 | 2014-02-26 | 中国航空工业集团公司北京航空制造工程研究所 | Method and system for measuring energy density of multi-dimensional electronic beams based on marginal integration |
| CN106596112A (en) * | 2016-12-16 | 2017-04-26 | 哈尔滨工业大学 | Hall thruster ion energy distribution measuring method and system |
| CN108883304A (en) * | 2018-06-22 | 2018-11-23 | 新瑞阳光粒子医疗装备(无锡)有限公司 | Synchrotron control method, device, equipment and storage medium |
| CN110708856A (en) * | 2018-11-26 | 2020-01-17 | 新瑞阳光粒子医疗装备(无锡)有限公司 | Particle accelerator, determination method and device for energy of extracted particles of particle accelerator and medium |
Non-Patent Citations (2)
| Title |
|---|
| 姚则悟等: "测量迥旋加速器束流能量的一种简易方法", 《原子能科学技术》 * |
| 岳军会等: "ADS注入器I束流相位及能量测量系统设计", 《核技术》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110708856A (en) | 2020-01-17 |
| WO2020107146A1 (en) | 2020-06-04 |
| CN110708856B (en) | 2021-07-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102762023B (en) | Circular accelerator and operation method thereof | |
| Popov et al. | Noise-induced temporary threshold shift and recovery in Yangtze finless porpoises Neophocaena phocaenoides asiaeorientalis | |
| EP1530410A2 (en) | Particle therapy system | |
| JP5224454B2 (en) | Abnormal tissue detection device | |
| US8772733B2 (en) | Charged particle accelerator and particle beam therapy system | |
| CN107920412B (en) | Method and system for controlling ion beam pulse extraction | |
| CN109661860A (en) | Particle accelerator and its determination method, apparatus and medium for drawing particle energy | |
| CN104053475A (en) | System for irradiating charged particles and method for irradiating charged particles | |
| CA3166860A1 (en) | Time-domain analysis of signals for charge detection mass spectrometry | |
| CN211061706U (en) | Flight time measuring system for measuring energy of pulse hadron beam and radiotherapy equipment | |
| CN110320485B (en) | Device and method for measuring gradient delay and first-order field unevenness and storage medium | |
| Kühteubl et al. | Investigating alternative extraction methods at MedAustron | |
| US10357182B2 (en) | Method and device for controlling MRI noise | |
| CN104939867A (en) | Measurement apparatus, measurement method, and ultrasonic measurement apparatus | |
| Fischer et al. | Tune Modulation from Beam‐Beam Interaction and Unequal Radio Frequencies in RHIC | |
| Reiter et al. | Improved Signal Treatment for Capacitive Linac Pickups | |
| Tobiyama et al. | Bunch by bunch feedback systems for J-PARC MR | |
| Chapman-Hatchett et al. | Tune measurement for the CERN Proton Synchrotron Booster rings using DSP in VME | |
| Ohshima et al. | Timing synchronization system for beam injection from the SACLA linac to the SPring-8 storage ring | |
| Nakamura et al. | Intra-bunch feedback system for the J-PARC main ring | |
| Cavazzana et al. | Real time control of electron density on rfx-mod tokamak discharges | |
| Gilpatrick et al. | Proposed beam position and phase measurements for the LANSCE linac | |
| CN116318131A (en) | Circuit and method for solving analog scanning frequency continuity | |
| Zinov et al. | Time structure of particle beams obtained from the synchrocyclotron in the United Institute of Nuclear Research (OIYaI) | |
| Scheidt | First results with the prototypes of new BPM electronics for the booster of the ESRF |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190419 |
|
| WD01 | Invention patent application deemed withdrawn after publication |