[A survey of performance of public health risk assessment in emergencies of institutions for disease control and prevention at different levels in China].

Objective: To understand the performance of public health risk assessment in emergencies of institutions for disease control and prevention at different levels in China, and provide suggestions for the improvement of public health risk assessment. Methods: A self-administered survey was conducted in professionals involved in public health risk assessment in emergencies from national institution, provincial institutions and some prefectural institutions for disease control and prevention (1-2 prefectural institutions were selected using convenience sampling in each province) between March and April in 2021. Results: A total of 79 institutions for disease control and prevention were investigated, including 1 national institution, 32 provincial institutions and 46 prefectural institutions. By April 2021, all the 79 institutions surveyed had conducted risk assessment of public health emergencies, in which 61 (77.2%) had established departments responsible for the public health risk assessment, i.e. emergency management office or communicable disease prevention and control office (section), and regular risk assessment mechanisms. The main sources of information for public health risk assessment were public health surveillance systems, including the National Notifiable Diseases Reporting System (100.0%) and Public Health Emergencies Management Information System (97.5%). Compared with the provincial institutions, the prefectural institutions were more likely to use specific disease surveillance systems (84.8% vs. 62.5%; χ2=5.09, P=0.024). The risk management recommendations made by 43 institutions for disease control and prevention (54.4%) after the risk assessment were accepted by the superior health administrative departments and used in epidemic prevention and control. Conclusions: Public health risk assessment in emergencies has been widely carried out by national, provincial and prefectural institutions for disease control and prevention in China. Specialized departments and mechanisms have been established, but the information sources are still confined to public health surveillance systems and the application of the risk assessment results still needs to be further improved.

Assessment of geometrical accuracy of magnetic resonance images for radiation therapy of lung cancers.

The purpose of this research was to investigate the geometrical accuracy of magnetic resonance (MR) images used in the radiation therapy treatment planning for lung cancer. In this study, the capability of MR imaging to acquire dynamic two-dimensional images was explored to access the motion of lung tumors. Due to a number of factors, including the use of a large field-of-view for the thorax, MR images are particularly subject to geometrical distortions caused by the inhomogeneity and gradient nonlinearity of the magnetic field. To quantify such distortions, we constructed a phantom, which approximated the dimensions of the upper thorax and included two air cavities. Evenly spaced vials containing contrast agent could be held in three directions with their cross-sections in the coronal, sagittal, and axial planes, respectively, within the air cavities. MR images of the phantom were acquired using fast spin echo (FSE) and fast gradient echo (fGRE) sequences. The positions of the vials according to their centers of mass were measured from the MR images and registered to the corresponding computed tomography images for comparison. Results showed the fGRE sequence exhibited no errors >2.0 mm in the sagittal and coronal planes, whereas the FSE sequence produced images with errors between 2.0 and 4.0 mm along the phantom's perimeter in the axial plane. On the basis of these results, the fGRE sequence was considered to be clinically acceptable in acquiring images in all sagittal and coronal planes tested. However, the spatial accuracy in periphery of the axial FSE images exceeded the acceptable criteria for the acquisition parameters used in this study.

A method of simulating dynamic multileaf collimators using Monte Carlo techniques for intensity-modulated radiation therapy.

A method of modelling the dynamic motion of multileaf collimators (MLCs) for intensity-modulated radiation therapy (IMRT) was developed and implemented into the Monte Carlo simulation. The simulation of the dynamic MLCs (DMLCs) was based on randomizing leaf positions during a simulation so that the number of particle histories being simulated for each possible leaf position was proportional to the monitor units delivered to that position. This approach was incorporated into an EGS4 Monte Carlo program, and was evaluated in simulating the DMLCs for Varian accelerators (Varian Medical Systems, Palo Alto. CA, USA). The MU index of each segment, which was specified in the DMLC-control data, was used to compute the cumulative probability distribution function (CPDF) for the leaf positions. This CPDF was then used to sample the leaf positions during a real-time simulation, which allowed for either the step-shoot or sweeping-leaf motion in the beam delivery. Dose intensity maps for IMRT fields were computed using the above Monte Carlo method, with its accuracy verified by film measurements. The DMLC simulation improved the operational efficiency by eliminating the need to simulate multiple segments individually. More importantly, the dynamic motion of the leaves could be simulated more faithfully by using the above leaf-position sampling technique in the Monte Carlo simulation.

Incorporating dynamic collimator motion in Monte Carlo simulations: an application in modelling a dynamic wedge.

In radiation therapy, new treatment modalities employing dynamic collimation and intensity modulation increase the complexity of dose calculation because a new dimension, time, has to be incorporated into the traditional three-dimensional problem. In this work, we investigated two classes of sampling technique to incorporate dynamic collimator motion in Monte Carlo simulation. The methods were initially evaluated for modelling enhanced dynamic wedges (EDWs) from Varian accelerators (Varian Medical Systems, Palo Alto, USA). In the position-probability-sampling or PPS method, a cumulative probability distribution function (CPDF) was computed for the collimator position, which could then be sampled during simulations. In the static-component-simulation or SCS method, a dynamic field is approximated by multiple static fields in a step-shoot fashion. The weights of the particles or the number of particles simulated for each component field are computed from the probability distribution function (PDF) of the collimator position. The CPDF and PDF were computed from the segmented treatment tables (STTs) for the EDWs. An output correction factor had to be applied in this calculation to account for the backscattered radiation affecting monitor chamber readings. Comparison of the phase-space data from the PPS method (with the step-shoot motion) with those from the SCS method showed excellent agreement. The accuracy of the PPS method was further verified from the agreement between the measured and calculated dose distributions. Compared to the SCS method, the PPS method is more automated and efficient from an operational point of view. The principle of the PPS method can be extended to simulate other dynamic motions, and in particular, intensity-modulated beams using multileaf collimators.

Prevalence and predictors of highly active antiretroviral therapy use in patients with HIV infection in the united states. HCSUS Consortium. HIV Cost and Services Utilization.

BACKGROUND: Highly active antiretroviral therapy (HAART) became standard for HIV in 1996. Studies at that time showed that most people infected with HIV had initiated HAART, but that members of minority groups and poor people had lower HAART use. It is not known whether high levels of HAART use have been sustained or whether socioeconomic and racial disparities have diminished over time. OBJECTIVES: To determine the proportion of patients who had received and were receiving HAART by January 1998, and to evaluate predictors of HAART receipt. DESIGN AND PARTICIPANTS: Prospective cohort study of a national probability sample of 2267 adults receiving HIV care who completed baseline, first follow-up, and second follow-up interviews from January 1996 to January 1998. MAIN OUTCOME VARIABLES: Proportion currently using HAART at second follow-up (August 1997 to January 1998), contrasted with the cumulative proportions using HAART at any time before January 1998 and before December 1996. ANALYSES: Bivariate and multiple logistic regression analysis of population characteristics predicting current use of HAART at the time of the second follow-up interview. RESULTS: The proportion of patients ever having received HAART increased from 37% in December 1996 to 71% by January 1998, but only 53% of people were receiving HAART at the time of the second follow-up interview. Differences between sociodemographic groups in ever using HAART narrowed after 1996. In bivariate analysis, several groups remained significantly less likely to be using HAART at the time of the second follow-up interview: blacks, male and female drug users, female heterosexuals, people with less education, those uninsured and insured by Medicaid, those in the Northeast, and those with CD4 counts of >/=500 cells/microl (all p <.05). Using multiple logistic regression analysis, low CD4 count (for CD4 <50 cells/microl: odds ratio [OR], 3.20; p <.001) remained a significant predictor of current HAART use at the time of the second follow-up interview, but lack of insurance (OR, 0.71; p <.05) predicted not receiving HAART. CONCLUSIONS: The proportion of persons under HIV care in the United States who had ever received HAART increased to over 70% of the affected population by January 1998 and the disparities in use between groups narrowed but did not disappear. However, nearly half of those eligible for HAART according to the U.S. Department of Health and Human Services guidelines were not actually receiving it nearly 2 years after these medications were first introduced. Strategies to promote the initiation and continuation of HAART are needed for those without contraindications and those who can tolerate it.

Backscatter towards the monitor ion chamber in high-energy photon and electron beams: charge integration versus Monte Carlo simulation.

In some linear accelerators, the charge collected by the monitor ion chamber is partly caused by backscattered particles from accelerator components downstream from the chamber. This influences the output of the accelerator and also has to be taken into account when output factors are derived from Monte Carlo simulations. In this work, the contribution of backscattered particles to the monitor ion chamber response of a Varian 2100C linac was determined for photon beams (6, 10 MV) and for electron beams (6, 12, 20 MeV). The experimental procedure consisted of charge integration from the target in a photon beam or from the monitor ion chamber in electron beams. The Monte Carlo code EGS4/BEAM was used to study the contribution of backscattered particles to the dose deposited in the monitor ion chamber. Both measurements and simulations showed a linear increase in backscatter fraction with decreasing field size for photon and electron beams. For 6 MV and 10 MV photon beams, a 2-3% increase in backscatter was obtained for a 0.5 x 0.5 cm2 field compared to a 40 x 40 cm2 field. The results for the 6 MV beam were slightly higher than for the 10 MV beam. For electron beams (6, 12, 20 MeV), an increase of similar magnitude was obtained from measurements and simulations for 6 MeV electrons. For higher energy electron beams a smaller increase in backscatter fraction was found. The problem is of less importance for electron beams since large variations of field size for a single electron energy usually do not occur.

Modeling photon output caused by backscattered radiation into the monitor chamber from collimator jaws using a Monte Carlo technique.

Dose per monitor unit in photon fields generated by clinical linear accelerators can be affected by the backscattered radiation into the monitor chamber from collimator jaws. Thus, it is necessary to account for the backscattered radiation in computing monitor unit setting for a treatment field. In this work, we investigated effects of the backscatter from collimator jaws based on Monte Carlo simulations of a clinical linear accelerator. The backscattered radiation scored within the monitor chamber was identified as originating either from the upper jaws (Y jaws), or from the lower jaws (X jaws). From the results of Monte Carlo simulations, ratios of the monitor-chamber-scored dose caused by the backscatter to the dose caused by the forward radiation, R(x,y), were modeled as functions of the individual X and Y jaw positions. The amount of the backscattered radiation for any field setting was then computed as a compound contribution from both the X and Y jaws. The dose ratios of R(x,y) were then used to calculate the change in photon output caused by the backscatter, Scb(x,y). Results of these calculations were compared with available measured data based on counting the electron pulses or charge from the electron target of an accelerator. Data from this study showed that the backscattered radiation contributes approximately 3% to the monitor-chamber-scored dose. A majority of the backscattered radiation comes from the upper jaws, which are located closer to the monitor chamber. The amount of the backscatter decreases approximately in a linear fashion with the jaw opening. This results in about a 2% increase of photon output from a 10 cm x 10 cm field to a 40 cm x 40 cm field. The off-axis location of the jaw opening does not have a significant effect on the magnitude of the backscatter. The backscatter effect is significant for monitor chambers using kapton windows, particularly for treatment fields using moving jaws. Applying the backscatter correction improves the accuracy of monitor-unit calculation using a model-based dose calculation algorithm such as the convolution method.

Calculating dose and output factors for wedged photon radiotherapy fields using a convolution/superposition method.

We have developed a convolution/superposition method to calculate dose distributions in photon treatment fields with beam modifiers such as physical wedges. The dose component due to wedge generated radiation was accounted for by using an extended phantom model, which integrated a wedge, an air gap, and a patient phantom as the calculation phantom. The inhomogeneities in the extended phantom and the effect of beam hardening by the wedge were both corrected for in the convolution dose calculation. The calculated dose was verified by Monte Carlo simulation of the same extended phantom. A new dual photon source model was also used in the convolution method to account for both primary photons from the target and extra-focal photons from the primary collimator and flattening filter. Thus, realistic photon energy fluence distributions in the extended phantom were used for the dose calculation. The calculated dose distributions and the wedge factors agreed with the measured data within 2% for a variety of treatment fields including asymmetric fields. Our results showed that the wedge-generated radiation could contribute a significant fraction of the total dose in patients. This dose component depends on a specific field configuration, thus wedge factor changes with photon energy, wedge angle, field size, depth, and patient phantom SSD. The variation of the wedge factor can be predicted accurately by our convolution approach with the extended phantom model, which allows for more accurate dose or monitor unit computation for photon fields with beam modifiers.

Correcting kernel tilting and hardening in convolution/superposition dose calculations for clinical divergent and polychromatic photon beams.

To account for clinical divergent and polychromatic photon beams, we have developed kernel tilting and kernel hardening correction methods for convolution dose calculation algorithms. The new correction methods were validated by Monte Carlo simulation. The accuracy and computation time of the our kernel tilting and kernel hardening correction methods were also compared to the existing approaches including terma divergence correction, dose divergence correction methods, and the effective mean kernel method with no kernel hardening correction. Treatment fields of 10 x 10-40 x 40 cm2 (field size at source to axis distance (SAD)) with source to source distances (SSDs) of 60, 80, and 100 cm, and photon energies of 6, 10, and 18 MV have been studied. Our results showed that based on the relative dose errors at a depth of 15 cm along the central axis, the terma divergence correction may be used for fields smaller than 10 x 10 cm2 with a SSD larger than 80 cm; the dose divergence correction with an additional kernel hardening correction can reduce dose error and may be more applicable than the terma divergence correction. For both these methods, the dose error increased linearly with the depth in the phantom; the 90% isodose lines at the depth of 15 cm were shifted by about 2%-5% of the field width due to significant underestimation of the penumbra dose. The kernel hardening effect was less prominent than the kernel tilting effect for clinical photon beams. The dose error by using nonhardening corrected kernel is less than 2.0% at a depth of 15 cm along the central axis, yet it increased with a smaller field size and lower photon energy. The kernel hardening correction could be more important to compute dose in the fields with beam modifiers such as wedges when beam hardening is more significant. The kernel tilting correction and kernel hardening correction increased computation time by about 3 times, and 0.5-1 times, respectively. This can be justified by more accurate dose calculations for the majority of clinical treatments.

A dual source photon beam model used in convolution/superposition dose calculations for clinical megavoltage x-ray beams.

A realistic model of photon beams generated by clinical linear accelerators has been incorporated in a convolution/superposition method to compute dose distributions in photon treatment fields. In this beam model, a primary photon source represents photons directly from the target, and an extra-focal photon source represents scattered photons from the primary collimator and the flattening filter. Monte Carlo simulation was used to study clinical linear accelerators producing photon beams. From the output of the Monte Carlo simulation, the fluence and spectral distributions of each photon component, as well as the geometrical characteristics of each photon source with respect to its distance to the isocenter and its source distribution, were analyzed. These quantities were used to reproduce realistic photon distributions in treatment fields, and thus to compute dose distributions using the convolution method. Our results showed that compared to the primary photon fluence, the extra-focal photon fluence from the primary collimator and the flattening filter was 11%-16% at the isocenter, among which 70% was contributed by the flattening filter. The variation of extra-focal photons in different treatment fields was predicted accurately by accounting for the finite size of the extra-focal source. Compared to measurements, dose distributions in photon treatment fields, including those of asymmetric jaw settings and at different SSDs were calculated accurately, particularly in the penumbral region, by using the convolution method with the new dual source photon beam model.

Calculating output factors for photon beam radiotherapy using a convolution/superposition method based on a dual source photon beam model.

A realistic photon beam model based on Monte Carlo simulation of clinical linear accelerators was implemented in a convolution/superposition dose calculation algorithm. A primary and an extra-focal sources were used in this beam model to represent the direct photons from the target and the scattered photons from other head structures, respectively. The effect of the finite size of the extra-focal source was modeled by a convolution of the source fluence distribution with the collimator aperture function. Relative photon output in air (Sc) and in phantom (Scp) were computed using the convolution method with this new photon beam model. Our results showed that in a 10 MV photon beam, the Sc, Sp (phantom scatter factor), and Scp factors increased by 11%, 10%, and 22%, respectively, as the field size changed from 3 x 3 cm2 to 40 x 40 cm2. The variation of the Sc factor was contributed mostly by an increase of the extra-focal radiation with field size. The radiation backscattered into the monitor chamber inside the accelerator head affected the Sc by about 2% in the same field range. The output factors in elongated fields, asymmetric fields, and blocked fields were also investigated in this study. Our results showed that if the effect of the backscattered radiation was taken into account, output factors in these treatment fields can be predicted accurately by our convolution algorithm using the dual source photon beam model.