THERMUS—A thermal model package for ROOT Wheaton, S.; Cleymans, J.; Hauer, M.
Computer physics communications,
2009, 2009-1-00, Letnik:
180, Številka:
1
Journal Article
Recenzirano
Odprti dostop
THERMUS is a package of C++ classes and functions allowing statistical-thermal model analyses of particle production in relativistic heavy-ion collisions to be performed within the ROOT framework of ...analysis. Calculations are possible within three statistical ensembles; a grand-canonical treatment of the conserved charges
B,
S and
Q, a fully canonical treatment of the conserved charges, and a mixed-canonical ensemble combining a canonical treatment of strangeness with a grand-canonical treatment of baryon number and electric charge. THERMUS allows for the assignment of decay chains and detector efficiencies specific to each particle yield, which enables sensible fitting of model parameters to experimental data.
Program title: THERMUS, version 2.1
Catalogue identifier: AEBW_v1_0
Program summary URL:
http://cpc.cs.qub.ac.uk/summaries/AEBW_v1_0.html
Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland
Licensing provisions: Standard CPC licence,
http://cpc.cs.qub.ac.uk/licence/licence.html
No. of lines in distributed program, including test data, etc.: 17 152
No. of bytes in distributed program, including test data, etc.: 93 581
Distribution format: tar.gz
Programming language: C++
Computer: PC, Pentium 4, 1 GB RAM (not hardware dependent)
Operating system: Linux: FEDORA, RedHat, etc.
Classification: 17.7
External routines: Numerical Recipes in C 1, ROOT 2
Nature of problem: Statistical-thermal model analyses of heavy-ion collision data require the calculation of both primordial particle densities and contributions from resonance decay. A set of thermal parameters (the number depending on the particular model imposed) and a set of thermalized particles, with their decays specified, is required as input to these models. The output is then a complete set of primordial thermal quantities for each particle, together with the contributions to the final particle yields from resonance decay. In many applications of statistical-thermal models it is required to fit experimental particle multiplicities or particle ratios. In such analyses, the input is a set of experimental yields and ratios, a set of particles comprising the assumed hadron resonance gas formed in the collision and the constraints to be placed on the system. The thermal model parameters consistent with the specified constraints leading to the best-fit to the experimental data are then output.
Solution method: THERMUS is a package designed for incorporation into the ROOT 2 framework, used extensively by the heavy-ion community. As such, it utilizes a great deal of ROOT's functionality in its operation. ROOT features used in THERMUS include its containers, the wrapper TMinuit implementing the MINUIT fitting package, and the TMath class of mathematical functions and routines. Arguably the most useful feature is the utilization of CINT as the control language, which allows interactive access to the THERMUS objects. Three distinct statistical ensembles are included in THERMUS, while additional options to include quantum statistics, resonance width and excluded volume corrections are also available. THERMUS provides a default particle list including all mesons (up to the
K
4
∗
(2045)) and baryons (up to the
Ω
−
) listed in the July 2002 Particle Physics Booklet 3. For each typically unstable particle in this list, THERMUS includes a text-file listing its decays. With thermal parameters specified, THERMUS calculates primordial thermal densities either by performing numerical integrations or else, in the case of the Boltzmann approximation without resonance width in the grand-canonical ensemble, by evaluating Bessel functions. Particle decay chains are then used to evaluate experimental observables (i.e. particle yields following resonance decay). Additional detector efficiency factors allow fine-tuning of the model predictions to a specific detector arrangement. When parameters are required to be constrained, use is made of the ‘Numerical Recipes in C’ 1 function which applies the Broyden globally convergent secant method of solving nonlinear systems of equations. Since the NRC software is not freely-available, it has to be purchased by the user. THERMUS provides the means of imposing a large number of constraints on the chosen model (amongst others, THERMUS can fix the baryon-to-charge ratio of the system, the strangeness density of the system and the primordial energy per hadron). Fits to experimental data are accomplished in THERMUS by using the ROOT TMinuit class. In its default operation, the standard
χ
2
function is minimized, yielding the set of best-fit thermal parameters. THERMUS allows the assignment of separate decay chains to each experimental input. In this way, the model is able to match the specific feed-down corrections of a particular data set.
Running time: Depending on the analysis required, run-times vary from seconds (for the evaluation of particle multiplicities given a set of parameters) to several minutes (for fits to experimental data subject to constraints).
References:
1
W.H. Press, S.A. Teukolsky, W.T. Vetterling, B.P. Flannery, Numerical Recipes in C: The Art of Scientific Computing, Cambridge University Press, Cambridge, 2002.
2
R. Brun, F. Rademakers, Nucl. Inst. Meth. Phys. Res. A 389 (1997) 81. See also
http://root.cern.ch/.
3
K. Hagiwara et al., Phys. Rev. D 66 (2002) 010001.
To date, projections of human migration induced by sea-level change (SLC) largely suggest large-scale displacement away from vulnerable coastlines. However, results from our model of Bangladesh ...suggest counterintuitively that people will continue to migrate toward the vulnerable coastline irrespective of the flooding amplified by future SLC under all emissions scenarios until the end of this century. We developed an empirically calibrated agent-based model of household migration decision-making that captures the multi-faceted push, pull and mooring influences on migration at a household scale. We then exposed ∼4800 000 simulated migrants to 871 scenarios of projected 21st-century coastal flooding under future emissions pathways. Our model does not predict flooding impacts great enough to drive populations away from coastlines in any of the scenarios. One reason is that while flooding does accelerate a transition from agricultural to non-agricultural income opportunities, livelihood alternatives are most abundant in coastal cities. At the same time, some coastal populations are unable to migrate, as flood losses accumulate and reduce the set of livelihood alternatives (so-called 'trapped' populations). However, even when we increased access to credit, a commonly-proposed policy lever for incentivizing migration in the face of climate risk, we found that the number of immobile agents actually rose. These findings imply that instead of a straightforward relationship between displacement and migration, projections need to consider the multiple constraints on, and preferences for, mobility. Our model demonstrates that decision-makers seeking to affect migration outcomes around SLC would do well to consider individual-level adaptive behaviors and motivations that evolve through time, as well as the potential for unintended behavioral responses.
This report provides a review of early and late effects of radiation in normal tissues and organs with respect to radiation protection. It was instigated following a recommendation in Publication 103 ...(ICRP, 2007), and it provides updated estimates of ‘practical’ threshold doses for tissue injury defined at the level of 1% incidence. Estimates are given for morbidity and mortality endpoints in all organ systems following acute, fractionated, or chronic exposure. The organ systems comprise the haematopoietic, immune, reproductive, circulatory, respiratory, musculoskeletal, endocrine, and nervous systems; the digestive and urinary tracts; the skin; and the eye.
Particular attention is paid to circulatory disease and cataracts because of recent evidence of higher incidences of injury than expected after lower doses; hence, threshold doses appear to be lower than previously considered. This is largely because of the increasing incidences with increasing times after exposure. In the context of protection, it is the threshold doses for very long follow-up times that are the most relevant for workers and the public; for example, the atomic bomb survivors with 40–50years of follow-up. Radiotherapy data generally apply for shorter follow-up times because of competing causes of death in cancer patients, and hence the risks of radiation-induced circulatory disease at those earlier times are lower.
A variety of biological response modifiers have been used to help reduce late reactions in many tissues. These include antioxidants, radical scavengers, inhibitors of apoptosis, anti-inflammatory drugs, angiotensin-converting enzyme inhibitors, growth factors, and cytokines. In many cases, these give dose modification factors of 1.1–1.2, and in a few cases 1.5–2, indicating the potential for increasing threshold doses in known exposure cases. In contrast, there are agents that enhance radiation responses, notably other cytotoxic agents such as antimetabolites, alkylating agents, anti-angiogenic drugs, and antibiotics, as well as genetic and comorbidity factors.
Most tissues show a sparing effect of dose fractionation, so that total doses for a given endpoint are higher if the dose is fractionated rather than when given as a single dose. However, for reactions manifesting very late after low total doses, particularly for cataracts and circulatory disease, it appears that the rate of dose delivery does not modify the low incidence. This implies that the injury in these cases and at these low dose levels is caused by single-hit irreparable-type events. For these two tissues, a threshold dose of 0.5Gy is proposed herein for practical purposes, irrespective of the rate of dose delivery, and future studies may elucidate this judgement further.
Gabapentin was used for the treatment of term and preterm infants with suspected visceral hyperalgesia caused by a variety of neurologic and gastrointestinal morbidities. Improved feeding tolerance ...and decreased irritability were seen, as well as decreased usage of opioids and benzodiazepines. Adverse events occurred with abrupt discontinuation of this medication.
Potassium (K) fertilization is important to maintain adequate concentrations of plant available K in agricultural soils to achieve best yields and improve crop stress tolerance and water‐use ...efficiency (WUE). Water‐use efficiency (WUE) can be expressed on various spatiotemporal scales, and it is known that responses of WUE to external stress are not uniform across scales. Multiscale evaluations of the impact of varying K fertilization on the WUE of C3 crops under field conditions are missing so far. In the present field study, we evaluated effects of K fertilization on WUE of sugar beet (Beta vulgaris L.) on short‐termed leaf‐ (WUELeaf) and canopy‐scales (WUECanopy) and as the agronomic ratio of white sugar yield (WSY) to in‐season water use (i.e. WUEWSY). In K‐fertilized plots, WUEWSY was enhanced by 15.9%. This effect is attributed to increased beet yield and WSY, as no differences in total in‐season water use between fertilized and unfertilized plots were observed. Potassium (K) fertilization significantly enhanced the leaf area index, resulting in a more efficient depletion of soil moisture by roots in K‐fertilized plots. As a consequence, WUELeaf was increased due to stomatal adjustment. Potassium (K) improved WUECanopy only by tendency. It is concluded that K fertilization improves the WUE of field‐grown sugar beet across scales, but processes that regulate WUE are highly scale dependent.
Children with severe impairment of the central nervous system (CNS) have a high incidence of distressing symptoms, with many experiencing frequent recurrent pain episodes.
The study objective was to ...describe presenting pain behaviors, daily dose, and response to gabapentin for the management of frequent recurrent pain in this population.
A retrospective analysis was performed with data from 22 children with severe impairment of the CNS residing at a long-term care facility, treated with gabapentin for recurrent pain behaviors. Response was considered significant if the frequency and severity of symptoms decreased by more than 50% as assessed by nursing staff.
Pain behaviors commonly reported included facial grimacing, crying, or moaning. Intermittent increase in muscle tone was identified in 86% (n=19). Gastrointestinal (GI) symptoms occurred in 64% (n=14), including pain localized to the GI tract and vomiting. All were assessed for nociceptive pain sources, many with repeated testing. Most were on medications for spasticity (n=20, 91%) and gastroesophageal reflux disease (GERD) (n=22, 100%) prior to gabapentin use. Of the 22 treated with gabapentin, 21 (91%) had a significant decrease in symptoms. No serious adverse events occurred. The mean gabapentin dose for children five years of age or less (n=11) was 50 mg/kg/day (95% CI 45-56) compared to children older than 11 years (n=11) with a mean dose of 36 mg/kg/day (95% CI 34-38).
Gabapentin appears to be an effective treatment for children with severe impairment of the CNS and recurrent pain behaviors, including intermittent changes in muscle tone. Dosing information can guide treatment trials and future prospective studies.
Radiotherapy of thoracic and chest wall tumors, if all or part of the heart was included in the radiation field, can lead to radiation-induced heart disease (RIHD), a late and potentially severe side ...effect. RIHD presents clinically several years after irradiation and manifestations include accelerated atherosclerosis, pericardial and myocardial fibrosis, conduction abnormalities, and injury to cardiac valves. The pathogenesis of RIHD is largely unknown, and a treatment is not available. Hence, ongoing pre-clinical studies aim to elucidate molecular and cellular mechanisms of RIHD. Here, an overview of recent pre-clinical studies is given, and based on the results of these studies, potential targets for intervention in RIHD are discussed.
We report on recurrent pancytopenia over five years in two children with severe impairment of the central nervous system. Assessment by hematology did not identify an etiology, including bone marrow ...biopsy in one. Both patients had sustained normalized blood cell counts following interventions to maintain or return to a temperature above 33 °C. Acute cytopenias following medically induced and environmental hypothermia have been reported. Recurrent pancytopenia due to centrally mediated hypothermia in patients with severe neurological impairment is often not recognized, putting such children at risk for unnecessary testing and transfusions. We provide a practical approach to management that is feasible for caregivers in the home setting with suggestions for monitoring.
Radiation-induced heart disease (RIHD) is a potentially severe side effect of radiotherapy of thoracic and chest wall tumors if all or part of the heart was included in the radiation field. RIHD ...presents clinically several years after irradiation and manifestations include accelerated atherosclerosis, pericardial and myocardial fibrosis, conduction abnormalities, and injury to cardiac valves. There is no method to prevent or reverse these injuries when the heart is exposed to ionizing radiation. This paper presents an overview of recent studies that address the role of microvascular injury, endothelial dysfunction, mast cells, and the renin angiotensin system in animal models of cardiac radiation injury. These insights into the basic mechanisms of RIHD may lead to the identification of targets for intervention in this late radiotherapy side effect.