CN113536651A - Radiation source intensity reconstruction method based on reverse particle transport - Google Patents

Radiation source intensity reconstruction method based on reverse particle transport Download PDF

Info

Publication number
CN113536651A
CN113536651A CN202110672890.0A CN202110672890A CN113536651A CN 113536651 A CN113536651 A CN 113536651A CN 202110672890 A CN202110672890 A CN 202110672890A CN 113536651 A CN113536651 A CN 113536651A
Authority
CN
China
Prior art keywords
particle
particles
energy
transport
reverse
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.)
Granted
Application number
CN202110672890.0A
Other languages
Chinese (zh)
Other versions
CN113536651B (en
Inventor
不公告发明人
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhongke Chaojing Nanjing Technology Co ltd
Original Assignee
Zhongke Chaojing Nanjing Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhongke Chaojing Nanjing Technology Co ltd filed Critical Zhongke Chaojing Nanjing Technology Co ltd
Priority to CN202110672890.0A priority Critical patent/CN113536651B/en
Publication of CN113536651A publication Critical patent/CN113536651A/en
Application granted granted Critical
Publication of CN113536651B publication Critical patent/CN113536651B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/25Design optimisation, verification or simulation using particle-based methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/90Details of database functions independent of the retrieved data types
    • G06F16/901Indexing; Data structures therefor; Storage structures
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/08Probabilistic or stochastic CAD

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Measurement Of Radiation (AREA)

Abstract

The invention discloses a radiation source intensity reconstruction method based on reverse particle transport, which comprises the following steps: establishing a particle reverse transport nuclear database; acquiring particle energy and intensity distribution information at measuring equipment, and establishing a reverse Monte Carlo radiation transport calculation model by combining initial source spatial distribution information, medium image information and measuring equipment information; sampling to obtain source particles according to the particle energy and intensity distribution information at the measuring equipment, and simulating the particles by using a Monte Carlo reverse transport method until the particles are transported to a radiation source; and performing reverse transport simulation on all the particles, and performing normalization processing on the particle intensity distribution information of the radiation source obtained by statistics to obtain the particle source intensity distribution of the radiation source. The invention can effectively obtain the strong distribution of the particle source and can carry out accurate radiation source inversion.

Description

Radiation source intensity reconstruction method based on reverse particle transport
Technical Field
The invention relates to the application fields of nuclear energy and nuclear technology such as reactor, radiotherapy, radiation protection and the like, in particular to a radiation source intensity reconstruction method based on reverse particle transport.
Background
The intensity distribution of the radiation source directly determines the distribution condition of a radiation dose field, and accurately obtaining the intensity distribution of the radiation source is a key problem in the fields of radiation protection, reactor detection, radiotherapy dose monitoring and the like. In the practical application process, the problems that the radiation source intensity is high, the space is limited, the detection equipment is difficult to place and the like exist, the radiation source intensity is difficult to obtain through direct measurement, the transmission distribution can be formed after rays pass through the medium, the transmission distribution reflects the radiation source intensity distribution to a certain extent, therefore, the method for obtaining the radiation source intensity distribution through reverse reconstruction of the transmission distribution is researched, and the method has important significance on radiation protection, reactor detection, radiotherapy dosage monitoring and the like.
At present, a scattering kernel deconvolution iteration method is mostly adopted for radiation source intensity reconstruction, the method obtains radiation source intensity distribution through repeated iteration optimization based on analysis scattering kernel, although the calculation speed is high, the problem of low reconstruction precision exists, and therefore a set of accurate radiation source intensity reconstruction method needs to be developed.
Disclosure of Invention
The purpose of the invention is as follows: the invention aims to provide a radiation source intensity reconstruction method based on reverse particle transport, which aims to solve the problem that a radiation source is difficult to reconstruct accurately in the fields of radiation protection, reactor detection, radiotherapy dosage monitoring and the like.
The technical scheme is as follows: the invention discloses a radiation source intensity reconstruction method based on reverse particle transport, which comprises the following steps:
(1) establishing a particle reverse transport nuclear database;
(2) acquiring particle energy and intensity distribution information at measuring equipment, and establishing a reverse Monte Carlo radiation transport calculation model by combining initial source spatial distribution information, medium image information and measuring equipment information;
(3) sampling to obtain source particles according to the particle energy and intensity distribution information at the measuring equipment, and simulating the particles by using a Monte Carlo reverse transport method until the particles are transported to a radiation source;
(4) and performing reverse transport simulation on all the particles, and performing normalization processing on the particle intensity distribution information of the radiation source obtained by statistics to obtain the particle source intensity distribution of the radiation source.
In the step (1), the particle reverse transport nuclear database comprises incident particle energy and angle distribution information corresponding to particles after reaction with different elements.
In the step (1), the establishing of the particle reverse transport nuclear database specifically comprises the following steps:
(1.1) dividing the particle energy of 0 MeV-20 MeV into N groups, wherein N is an integer greater than 0;
(1.2) building a unit ball, filling different elements in the unit ball, setting a point source at the center of the unit ball, taking each energy group in the step (1.1) as source energy distribution, simulating by using a Monte Carlo method, and counting particle energy/angle distribution after one collision;
(1.3) after all energy groups are simulated, the energy/angle distribution data of the collided particles are processed according to different elements, the energy/angle distribution data of the incident particles corresponding to the energy of the emergent particles are obtained, and a particle reverse transport nuclear database is established.
In the step (2), the establishing of the reverse Monte Carlo radiation transport calculation model specifically comprises the following steps:
(2.1) calibrating the detection data of the measuring equipment by using the calibration information of the measuring data to obtain the energy and intensity distribution of particles at the measuring equipment;
and (2.2) establishing a reverse Monte Carlo radiation transport calculation model by taking the measuring equipment as a reverse particle transport source and the initial source space as a counting space.
In the step (3), the source particles are obtained by sampling according to the particle energy and intensity distribution information at the measuring equipment, and the particles are simulated by using a Monte Carlo reverse transport method until the particle transport reaches the radiation source, which specifically comprises the following steps:
(3.1) sampling to obtain source photon states including positions, energies and angles according to a reverse Monte Carlo radiation transport calculation model;
(3.2) according to the position, the direction and the energy of the particles, combining the information of the medium material, and sampling based on a probability density function to obtain the position and the post-collision state of the next collision point of the particles;
(3.3) simulating the current particles until the particles are transported to the space where the radiation source is located, and adding the weight of the particles transported to the space where the radiation source is located into the intensity distribution of the source particles;
(3.4) if the current particle is not transported to the radiation source space, the current particle weight is evenly distributed to the particle contributions generated by the same count space.
In step (3.2), the expression of the probability density function P (E, Ω, r → E ', Ω ', r ') is:
Figure BDA0003119426990000021
in the formula, E is the current energy of the particle, omega is the current angle of the particle, r is the current position of the particle, E ' is the energy before the particle collision, omega ' is the angle before the collision, r ' is the position of the last collision point of the particle, and sigmat(r, E') is the total cross-section of the reaction at r for particles of energy Es(r ', E ') is the scattering cross section of the particle with energy E ' at r ', and f (r '; E, omega → E ', omega ') is the inverse transport nuclear data of the particle at r
Has the advantages that: compared with the prior art, the invention has the beneficial effects that: (1) the radiation source intensity is reconstructed by adopting a reverse Monte Carlo method, a particle reverse transport nuclear database is established through Monte Carlo simulation, accurate nuclear data support is provided for particle reverse transport, and the problem that accurate particle reverse transport nuclear data are difficult to obtain is solved; (2) the source intensity reconstruction is carried out based on a reverse Monte Carlo method, and compared with the traditional analytic iteration method, the high-precision radiation source intensity distribution reconstruction can be realized; (3) the invention can effectively improve the reconstruction precision of the intensity distribution of the radiation source in the fields of radiation protection, reactor detection, radiotherapy dosage monitoring and the like.
Drawings
FIG. 1 is a schematic flow diagram of the present invention;
fig. 2 is a comparison graph of input picture and reconstruction source intensity obtained by the method of the present invention.
Detailed Description
The technical solution of the present invention is further described in detail below with reference to the detailed description and the accompanying drawings.
As shown in fig. 1, the present invention includes a radiation source intensity reconstruction method based on inverse particle transport, including the following steps:
(1) establishing a particle reverse transport nuclear database; the particle reverse transport nuclear database comprises incident particle energy and angle distribution information corresponding to particles after reaction with different elements; the method specifically comprises the following steps:
(1.1) dividing the particle energy of 0 MeV-20 MeV into N groups, wherein N is an integer greater than 0;
(1.2) establishing a unit sphere, filling different elements in the unit sphere, setting a point source at the center of the unit sphere, taking each energy group in the step (1.1) as source energy distribution, simulating by using a Monte Carlo method, and counting particle energy/angle distribution f (E → E ', mu) after one collision, wherein E is incident particle energy, E' is emergent particle energy, and mu is the cosine of an emergent angle;
(1.3) after all energy groups are simulated, the energy/angle distribution data of the collided particles are processed aiming at different elements, the energy/angle distribution data of the incident particles corresponding to the energy of the emergent particles are obtained, and a particle reverse transport nuclear database f (E' → E, mu) is established.
(2) Acquiring particle energy and intensity distribution information at measuring equipment, and establishing a reverse Monte Carlo radiation transport calculation model by combining initial source spatial distribution information, medium image information and measuring equipment information; the method for establishing the reverse Monte Carlo radiation transport calculation model specifically comprises the following steps:
(2.1) calibrating the detection data of the measuring equipment by using the calibration information of the measuring data to obtain the energy and intensity distribution of particles at the measuring equipment;
and (2.2) establishing a reverse Monte Carlo radiation transport calculation model by taking the measuring equipment as a reverse particle transport source and the initial source space as a counting space.
(3) Sampling to obtain source particles according to the particle energy and intensity distribution information at the measuring equipment, and simulating the particles by using a Monte Carlo reverse transport method until the particles are transported to a radiation source; the method specifically comprises the following steps:
(3.1) sampling based on a probability density function p (E, omega, r → E ', omega ', r ') according to the position, the direction and the energy of the particle and by combining medium material information to obtain the position and the post-collision state of the next collision point of the particle;
Figure BDA0003119426990000041
in the formula, E is the current energy of the particle, omega is the current angle of the particle, r is the current position of the particle, E ' is the energy before the particle collision, omega ' is the angle before the collision, r ' is the position of the last collision point of the particle, and sigmat(r, E') is the total cross-section of the reaction at r for particles of energy Es(r ', E') is the scattering cross section of the particle with energy E 'at r', and f (r '; E, omega → E', omega ') is the inverse transport nuclear data of the particle at r';
(3.2) simulating the current particles until the particles are transported to the space where the radiation source is located, and adding the weight of the particles transported to the space where the radiation source is located into the intensity distribution of the source particles;
(3.3) if the current particle is not transported to the radiation source space, the current particle weight is evenly distributed to the particle contributions generated by the same count space.
(4) And performing reverse transport simulation on all the particles, and performing normalization processing on the particle intensity distribution information of the radiation source obtained by statistics to obtain the particle source intensity distribution of the radiation source.
Examples
A typical lung cancer patient is subjected to radiotherapy radiation source intensity reconstruction. The lung cancer patient receives intensity modulated radiotherapy, and a total of 7 beams are irradiated, wherein the beam 1 is subjected to radiation source intensity reconstruction, the beam 1 is vertically incident to the surface of a human body, the distance between a source point and an isocenter of the human body is 100cm, the distance between the source point and a portal imaging device is 140cm, and a measurement image acquired by the portal imaging device is shown in fig. 2 (a). The specific reconstruction steps are as follows:
1. photon reverse transport nuclear database establishment
(1.1) dividing photon energy of 0 MeV-20 MeV into 200 groups, wherein the energy interval of each group is 0.1 MeV; (1.2) establishing a unit sphere, filling different elements (main elements related to human bodies) in the unit sphere, setting a point source at the center of the unit sphere, taking each energy group in the step (1.1) as source energy distribution, simulating by using a Monte Carlo method, and counting particle energy/angle distribution f (E → E ', mu) after one collision, wherein E is incident particle energy, E' is emergent particle energy, and mu is the cosine of an emergent angle;
(1.3) after simulating all energy group photons, processing photon energy/angle distribution data after collision aiming at different elements to obtain incident photon energy/angle distribution data f (E' → E, mu) corresponding to emergent photon energy, and establishing a photon reverse transport nuclear database.
2. Establishment of reverse Monte Carlo radiation transport calculation model
(2.1) calibrating the detection data shown in the figure 2(a) by using the calibration information of the measurement data to obtain photon energy and intensity distribution at the measuring equipment;
and (2.2) establishing a reverse photon transport model by taking the measuring equipment as a reverse photon transport source and the initial source space as a counting space.
3. Monte Carlo reverse photon transport simulation
(3.1) sampling to obtain source photon states including positions, energies and angles according to a reverse Monte Carlo radiation transport calculation model;
(3.2) according to the position, the direction and the energy of the particles, combining the information of the medium material, and sampling based on a probability density function to obtain the position and the post-collision state of the next collision point of the particles; wherein, the expression of the probability density function P (E, Ω, r → E ', Ω ', r ') is:
Figure BDA0003119426990000051
in the formula, E is the current energy of the particle, omega is the current angle of the particle, r is the current position of the particle, E ' is the energy before the particle collision, omega ' is the angle before the collision, r ' is the position of the last collision point of the particle, and sigmat(r, E ') is the total cross-section of the reaction at r for particles with energy E',∑s(r ', E') is the scattering cross section at r 'of the particle with energy E', and f (r '; E, Ω → E', Ω ') is the inverse transport nuclear data at r' of the particle.
(3.3) simulating the current particles until the particles are transported to the space where the radiation source is located, and adding the weight of the particles transported to the space where the radiation source is located into the intensity distribution of the source particles;
(3.4) if the current particle is not transported to the radiation source space, the current particle weight is evenly distributed to the particle contributions generated by the same count space.
4. Radiation source intensity distribution reconstruction
The photon intensity distribution information at the radiation source obtained by statistics is normalized to obtain the photon source intensity distribution at the radiation source, as shown in fig. 2 (b).

Claims (6)

1. A radiation source intensity reconstruction method based on reverse particle transport is characterized in that: the method comprises the following steps:
(1) establishing a particle reverse transport nuclear database;
(2) acquiring particle energy and intensity distribution information at measuring equipment, and establishing a reverse Monte Carlo radiation transport calculation model by combining initial source spatial distribution information, medium image information and measuring equipment information;
(3) sampling to obtain source particles according to the particle energy and intensity distribution information at the measuring equipment, and simulating the particles by using a Monte Carlo reverse transport method until the particles are transported to a radiation source;
(4) and performing reverse transport simulation on all the particles, and performing normalization processing on the particle intensity distribution information of the radiation source obtained by statistics to obtain the particle source intensity distribution of the radiation source.
2. The radiation source intensity reconstruction method based on reverse particle transport as claimed in claim 1, wherein in step (1), the particle reverse transport kernel database includes incident particle energy and angle distribution information corresponding to particles after reaction with different elements.
3. The radiation source intensity reconstruction method based on reverse particle transport according to claim 1, wherein in the step (1), the establishing of the particle reverse transport nuclear database specifically includes the following steps:
(1.1) dividing the particle energy of 0 MeV-20 MeV into N groups, wherein N is an integer greater than 0;
(1.2) building a unit ball, filling different elements in the unit ball, setting a point source at the center of the unit ball, taking each energy group in the step (1.1) as source energy distribution, simulating by using a Monte Carlo method, and counting particle energy/angle distribution after one collision;
(1.3) after all energy groups are simulated, the energy/angle distribution data of the collided particles are processed according to different elements, the energy/angle distribution data of the incident particles corresponding to the energy of the emergent particles are obtained, and a particle reverse transport nuclear database is established.
4. The radiation source intensity reconstruction method based on inverse particle transport according to claim 1, wherein in the step (2), the establishing of the inverse monte carlo radiation transport calculation model specifically includes the following steps:
(2.1) calibrating the detection data of the measuring equipment by using the calibration information of the measuring data to obtain the energy and intensity distribution of particles at the measuring equipment;
and (2.2) establishing a reverse Monte Carlo radiation transport calculation model by taking the measuring equipment as a reverse particle transport source and the initial source space as a counting space.
5. The radiation source intensity reconstruction method based on reverse particle transport according to claim 1, wherein in the step (3), the source particles are obtained by sampling according to the particle energy and intensity distribution information at the measuring device, and the monte carlo reverse transport method is used to simulate the particles until the particle transport reaches the radiation source, which specifically includes the following steps:
(3.1) sampling to obtain source photon states including positions, energies and angles according to a reverse Monte Carlo radiation transport calculation model;
(3.2) according to the position, the direction and the energy of the particles, combining the information of the medium material, and sampling based on a probability density function to obtain the position and the post-collision state of the next collision point of the particles;
(3.3) simulating the current particles until the particles are transported to the space where the radiation source is located, and adding the weight of the particles transported to the space where the radiation source is located into the intensity distribution of the source particles;
(3.4) if the current particle is not transported to the radiation source space, the current particle weight is evenly distributed to the particle contributions generated by the same count space.
6. The inverse particle transport-based radiation source intensity reconstruction method according to claim 5, wherein in step (3.2), the probability density function P (E, Ω, r → E ', Ω ', r ') is expressed as:
Figure FDA0003119426980000021
in the formula, E is the current energy of the particle, omega is the current angle of the particle, r is the current position of the particle, E ' is the energy before the particle collision, omega ' is the angle before the collision, r ' is the position of the last collision point of the particle, and sigmat(r, E') is the total cross-section of the reaction at r for particles of energy Es(r ', E') is the scattering cross section at r 'of the particle with energy E', and f (r '; E, Ω → E', Ω ') is the inverse transport nuclear data at r' of the particle.
CN202110672890.0A 2021-06-17 2021-06-17 Radiation source intensity reconstruction method based on reverse particle transport Active CN113536651B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110672890.0A CN113536651B (en) 2021-06-17 2021-06-17 Radiation source intensity reconstruction method based on reverse particle transport

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110672890.0A CN113536651B (en) 2021-06-17 2021-06-17 Radiation source intensity reconstruction method based on reverse particle transport

Publications (2)

Publication Number Publication Date
CN113536651A true CN113536651A (en) 2021-10-22
CN113536651B CN113536651B (en) 2024-04-02

Family

ID=78125098

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110672890.0A Active CN113536651B (en) 2021-06-17 2021-06-17 Radiation source intensity reconstruction method based on reverse particle transport

Country Status (1)

Country Link
CN (1) CN113536651B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114861510A (en) * 2022-03-31 2022-08-05 西安交通大学 Particle omni-directional energy spectrum incidence method based on Geant4

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10253448A (en) * 1997-03-06 1998-09-25 Fuji Electric Co Ltd Method for calculating and detecting light source position and device for detecting light source position
US5870697A (en) * 1996-03-05 1999-02-09 The Regents Of The University Of California Calculation of radiation therapy dose using all particle Monte Carlo transport
US6285969B1 (en) * 1998-05-22 2001-09-04 The Regents Of The University Of California Use of single scatter electron monte carlo transport for medical radiation sciences
CN103065018A (en) * 2013-01-13 2013-04-24 中国科学院合肥物质科学研究院 Reverse Monte Carlo particle transporting and simulating system
CN104376217A (en) * 2014-11-20 2015-02-25 中国科学院合肥物质科学研究院 Radiation shielding calculation method based on monte carlo self-adaptive variance reduction
CN106354896A (en) * 2015-07-17 2017-01-25 上海联影医疗科技有限公司 Method and device for determining intensity distribution of particles passing through beam adjusting device
CN112949156A (en) * 2021-03-25 2021-06-11 中科超精(南京)科技有限公司 Monte Carlo de-scattering correction method

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5870697A (en) * 1996-03-05 1999-02-09 The Regents Of The University Of California Calculation of radiation therapy dose using all particle Monte Carlo transport
JPH10253448A (en) * 1997-03-06 1998-09-25 Fuji Electric Co Ltd Method for calculating and detecting light source position and device for detecting light source position
US6285969B1 (en) * 1998-05-22 2001-09-04 The Regents Of The University Of California Use of single scatter electron monte carlo transport for medical radiation sciences
CN103065018A (en) * 2013-01-13 2013-04-24 中国科学院合肥物质科学研究院 Reverse Monte Carlo particle transporting and simulating system
CN104376217A (en) * 2014-11-20 2015-02-25 中国科学院合肥物质科学研究院 Radiation shielding calculation method based on monte carlo self-adaptive variance reduction
CN106354896A (en) * 2015-07-17 2017-01-25 上海联影医疗科技有限公司 Method and device for determining intensity distribution of particles passing through beam adjusting device
CN112949156A (en) * 2021-03-25 2021-06-11 中科超精(南京)科技有限公司 Monte Carlo de-scattering correction method

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114861510A (en) * 2022-03-31 2022-08-05 西安交通大学 Particle omni-directional energy spectrum incidence method based on Geant4
CN114861510B (en) * 2022-03-31 2024-04-09 西安交通大学 Geant 4-based particle omnidirectional energy spectrum incidence method

Also Published As

Publication number Publication date
CN113536651B (en) 2024-04-02

Similar Documents

Publication Publication Date Title
Draeger et al. 3D prompt gamma imaging for proton beam range verification
CN105204056B (en) Online dose monitoring and the method for checking in radiotherapy
Hurley et al. Water‐equivalent path length calibration of a prototype proton CT scanner
Yepes et al. A GPU implementation of a track-repeating algorithm for proton radiotherapy dose calculations
Pennazio et al. Carbon ions beam therapy monitoring with the INSIDE in-beam PET
EP3428690A1 (en) Radiation dose measuring method
WO2022078175A1 (en) Boron neutron capture therapy system and treatment plan generation method therefor
Sommerer et al. In-beam PET monitoring of mono-energetic 16O and 12C beams: experiments and FLUKA simulations for homogeneous targets
CN103091699A (en) Device and method for measuring strong gamma ray energy spectrum using scattering method
Piersimoni et al. The effect of beam purity and scanner complexity on proton CT accuracy
Robert et al. PET-based dose delivery verification in proton therapy: a GATE based simulation study of five PET system designs in clinical conditions
Schultze et al. Particle-tracking proton computed tomography—data acquisition, preprocessing, and preconditioning
Jones et al. Characterization of Compton-scatter imaging with an analytical simulation method
Yao et al. Rapid compton camera imaging for source terms investigation in the nuclear decommissioning with a subset-driven origin ensemble algorithm
Vanstalle et al. Benchmarking Geant4 hadronic models for prompt‐γ monitoring in carbon ion therapy
Masuda et al. ML-EM algorithm for dose estimation using PET in proton therapy
CN113536651B (en) Radiation source intensity reconstruction method based on reverse particle transport
Sohrabi et al. Photoneutron spectrometry by novel multi-directional spherical neutron spectrometry system
Zhang et al. Convolution neural network toward Monte Carlo photon dose calculation in radiation therapy
Brost et al. A mathematical deconvolution formulation for superficial dose distribution measurement by Cerenkov light dosimetry
CN112949156B (en) Monte Carlo scattering correction method
Beebe-Wang et al. Simulation of proton therapy treatment verification via PET imaging of induced positron-emitters
Manabe et al. Proof-of-concept study on a water phantom-based neutron spectrometer: Experimental test with 252Cf and 241Am-Be sources
Nasirzadeh et al. Modeling GE advance PET-scanner using FLUKA simulation code
Li et al. Photon energy spectrum reconstruction based on Monte Carlo and measured percentage depth dose in accurate radiotherapy

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
GR01 Patent grant
GR01 Patent grant