Calculating Dilepton Rates from Monte Carlo Simulations of Parton Production

To calculate dilepton rates in a Monte Carlo simulation of ultrarelativistic heavy ion collisions, one usually scales the number of similar QCD processes by a ratio of the corresponding differential probabilities. We derive the formula for such a ratio especially for dilepton bremsstrahlung processes. We also discuss the non-triviality of including higher order corrections to direct Drell-Yan process. The resultant mass spectra from our Monte Carlo simulation are consistent with the semi-analytical calculation using dilepton fragmentation functions.Comment: 14 pages in RevTex, 3 figures in uuencoded files, LBL-3466

Compton Scattering of Fe K\alpha Lines from Accreting White Dwarfs

Compton scattering in the bulk accretion flow of the accretion column in magnetic cataclysmic variables (mCVs) can significantly shift photon energies in the X-ray emission lines resulting from accretion shocks. In particular, Compton recoil can potentially broaden the 6.7 and 6.97 keV Fe K$\alpha$ emission lines produced in the post-shock region, and contaminate the fluorescent 6.4 keV neutral Fe K$\alpha$ line reflected off the white dwarf surface. We present nonlinear Monte Carlo simulations demonstrating these effects, and we discuss the interpretation of observed Fe K$\alpha$ linewidths in mCVs in light of these new results. The implications for other accreting compact objects are also discussed.Comment: 7 pages, 3 figs, LaTeX, accepted by PAS

Dosimetric evaluation of Acuros XB Advanced Dose Calculation algorithm in heterogeneous media

<p>Abstract</p> <p>Background</p> <p>A study was realised to evaluate and determine relative figures of merit of a new algorithm for photon dose calculation when applied to inhomogeneous media.</p> <p>Methods</p> <p>The new Acuros XB algorithm implemented in the Varian Eclipse treatment planning system was compared against a Monte Carlo method (VMC++), and the Analytical Anisotropic Algorithm (AAA). The study was carried out in virtual phantoms characterized by simple geometrical structures. An insert of different material and density was included in a phantom built of skeletal-muscle and HU = 0 (setting "A"): Normal Lung (lung, 0.198 g/cm³); Light Lung (lung, 0.035 g/cm³); Bone (bone, 1.798 g/cm³); another phantom (setting "B") was built of adipose material and including thin layers of bone (1.85 g/cm³), adipose (0.92 g/cm³), cartilage (1.4745 g/cm³), air (0.0012 g/cm³). Investigations were performed for 6 and 15 MV photon beams, and for a large (13 × 13 cm²) and a small (2.8 × 13 cm²) field.</p> <p>Results</p> <p>Results are provided in terms of depth dose curves, transverse profiles and Gamma analysis (3 mm/3% and 2 mm/2% distance to agreement/dose difference criteria) in planes parallel to the beam central axis; Monte Carlo simulations were assumed as reference. Acuros XB gave an average gamma agreement, with a 3 mm/3% criteria, of 100%, 86% and 100% for Normal Lung, Light Lung and Bone settings, respectively, and dose to medium calculations. The same figures were 86%, 11% and 100% for AAA, where only dose rescaled to water calculations are possible.</p> <p>Conclusions</p> <p>In conclusion, Acuros XB algorithm provides a valid and accurate alternative to Monte Carlo calculations for heterogeneity management.</p

Countering beam divergence effects with focused segmented scintillators for high DQE megavoltage active matrix imagers

The imaging performance of active matrix flat-panel imagers designed for megavoltage imaging (MV AMFPIs) is severely constrained by relatively low x-ray detection efficiency, which leads to a detective quantum efficiency (DQE) of only ∼1%. Previous theoretical and empirical studies by our group have demonstrated the potential for addressing this constraint through the utilization of thick, two-dimensional, segmented scintillators with optically isolated crystals. However, this strategy is constrained by the degradation of high-frequency DQE resulting from spatial resolution loss at locations away from the central beam axis due to oblique incidence of radiation. To address this challenge, segmented scintillators constructed so that the crystals are individually focused toward the radiation source are proposed and theoretically investigated. The study was performed using Monte Carlo simulations of radiation transport to examine the modulation transfer function and DQE of focused segmented scintillators with thicknesses ranging from 5 to 60 mm. The results demonstrate that, independent of scintillator thickness, the introduction of focusing largely restores spatial resolution and DQE performance otherwise lost in thick, unfocused segmented scintillators. For the case of a 60 mm thick BGO scintillator and at a location 20 cm off the central beam axis, use of focusing improves DQE by up to a factor of ∼130 at non-zero spatial frequencies. The results also indicate relatively robust tolerance of such scintillators to positional displacements, of up to 10 cm in the source-to-detector direction and 2 cm in the lateral direction, from their optimal focusing position, which could potentially enhance practical clinical use of focused segmented scintillators in MV AMFPIs.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/98594/1/0031-9155_57_16_5343.pd

The Jacobian as a measure of planar dose congruence

We propose a new starting point for comparing dose distributions in therapeutic radiation physics using a Jacobian-based measure. The measure is normalization independent, free of tunable parameters, bounded and converges to a unique value when comparing unrelated dose distributions. We present a preliminary demonstration of the sensitivity and general characteristics of this measure.Comment: 9 pages, 2 figure

Automated algorithm for CBCT-based dose calculations of prostate radiotherapy with bilateral hip prostheses

ABSTRACTOBJECTIVE:Cone beam CT (CBCT) images contain more scatter than a conventional CT image and therefore provide inaccurate Hounsfield units (HUs). Consequently, CBCT images cannot be used directly for radiotherapy dose calculation. The aim of this study is to enable dose calculations to be performed with the use of CBCT images taken during radiotherapy and evaluate the necessity of replanning.METHODS:A patient with prostate cancer with bilateral metallic prosthetic hip replacements was imaged using both CT and CBCT. The multilevel threshold (MLT) algorithm was used to categorize pixel values in the CBCT images into segments of homogeneous HU. The variation in HU with position in the CBCT images was taken into consideration. This segmentation method relies on the operator dividing the CBCT data into a set of volumes where the variation in the relationship between pixel values and HUs is small. An automated MLT algorithm was developed to reduce the operator time associated with the process. An intensity-modulated radiation therapy plan was generated from CT images of the patient. The plan was then copied to the segmented CBCT (sCBCT) data sets with identical settings, and the doses were recalculated and compared.RESULTS:Gamma evaluation showed that the percentage of points in the rectum with γ < 1 (3%/3 mm) were 98.7% and 97.7% in the sCBCT using MLT and the automated MLT algorithms, respectively. Compared with the planning CT (pCT) plan, the MLT algorithm showed −0.46% dose difference with 8 h operator time while the automated MLT algorithm showed −1.3%, which are both considered to be clinically acceptable, when using collapsed cone algorithm.CONCLUSION:The segmentation of CBCT images using the method in this study can be used for dose calculation. For a patient with prostate cancer with bilateral hip prostheses and the associated issues with CT imaging, the MLT algorithms achieved a sufficient dose calculation accuracy that is clinically acceptable. The automated MLT algorithm reduced the operator time associated with implementing the MLT algorithm to achieve clinically acceptable accuracy. This saved time makes the automated MLT algorithm superior and easier to implement in the clinical setting.ADVANCES IN KNOWLEDGE:The MLT algorithm has been extended to the complex example of a patient with bilateral hip prostheses, which with the introduction of automation is feasible for use in adaptive radiotherapy, as an alternative to obtaining a new pCT and reoutlining the structures

Space-time Structure of Initial Parton Production in Ultrarelativistic Heavy Ion Collisions

The space and time evolution of initial parton production in ultrarelativistic heavy ion collisions is investigated within the framework of perturbative QCD which includes both initial and final state radiations. Uncertainty principle is used to relate the life time of a radiating parton to its virtuality and momentum. The interaction time of each hard or semihard parton scattering is also taken into account. For central $Au+Au$ collisions at $\sqrt{s}=200$ GeV, most of the partons are found to be produced within 0.5 fm/c after the total overlap of the two colliding nuclei. The local momentum distribution is approximately isotropical at that time. The implication on how to treat correctly the the secondary scattering in an ultimate parton cascading model is also discussed.Comment: 19 pages in REVTEX with 12 figures in separate uuencoded postscript files, LBL-3415

Photon beam relative dose validation of the DPM Monte Carlo code in lung‐equivalent media

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/134925/1/mp5671.pd

R&D on co-working transport schemes in Geant4

A research and development (R&D) project related to the extension of the Geant4 toolkit has been recently launched to address fundamental methods in radiation transport simulation. The project focuses on simulation at different scales in the same experimental environment; this problem requires new methods across the current boundaries of condensed-random-walk and discrete transport schemes. The new developments have been motivated by experimental requirements in various domains, including nanodosimetry, astronomy and detector developments for high energy physics applications.Comment: To be published in the Proceedings of the CHEP (Computing in High Energy Physics) 2009 conferenc