Geant4 Hadronic Cascade Models and CMS Data Analysis : Computational Challenges in the LHC era

This work belongs to the field of computational high-energy physics (HEP). The key methods used in this thesis work to meet the challenges raised by the Large Hadron Collider (LHC) era experiments are object-orientation with software engineering, Monte Carlo simulation, the computer technology of clusters, and artificial neural networks. The first aspect discussed is the development of hadronic cascade models, used for the accurate simulation of medium-energy hadron-nucleus reactions, up to 10 GeV. These models are typically needed in hadronic calorimeter studies and in the estimation of radiation backgrounds. Various applications outside HEP include the medical field (such as hadron treatment simulations), space science (satellite shielding), and nuclear physics (spallation studies). Validation results are presented for several significant improvements released in Geant4 simulation tool, and the significance of the new models for computing in the Large Hadron Collider era is estimated. In particular, we estimate the ability of the Bertini cascade to simulate Compact Muon Solenoid (CMS) hadron calorimeter HCAL. LHC test beam activity has a tightly coupled cycle of simulation-to-data analysis. Typically, a Geant4 computer experiment is used to understand test beam measurements. Thus an another aspect of this thesis is a description of studies related to developing new CMS H2 test beam data analysis tools and performing data analysis on the basis of CMS Monte Carlo events. These events have been simulated in detail using Geant4 physics models, full CMS detector description, and event reconstruction. Using the ROOT data analysis framework we have developed an offline ANN-based approach to tag b-jets associated with heavy neutral Higgs particles, and we show that this kind of NN methodology can be successfully used to separate the Higgs signal from the background in the CMS experiment.

Root system traits of Norway spruce, Scots pine, and silver birch in mixed boreal forests: an analysis of root architecture, morphology, and anatomy

The aim of this thesis was to unravel the functional-structural characteristics of root systems of Betula pendula Roth., Picea abies (L.) Karst., and Pinus sylvestris L. in mixed boreal forest stands differing in their developmental stage and site fertility. The root systems of these species had similar structural regularities: horizontally-oriented shallow roots defined the horizontal area of influence, and within this area, each species placed fine roots in the uppermost soil layers, while sinker roots defined the maximum rooting depth. Large radial spread and high ramification of coarse roots, and the high specific root length (SRL) and root length density (RLD) of fine roots indicated the high belowground competitiveness and root plasticity of B. pendula. Smaller radial root spread and sparser branching of coarse roots, and low SRL and RLD of fine roots of the conifers could indicate their more conservative resource use and high association with and dependence on ectomycorrhiza-forming fungi. The vertical fine root distributions of the species were mostly overlapping, implying the possibility for intense belowground competition for nutrients. In each species, conduits tapered and their frequency increased from distal roots to the stem, from the stem to the branches, and to leaf petioles in B. pendula. Conduit tapering was organ-specific in each species violating the assumptions of the general vascular scaling model (WBE). This reflects the hierarchical organization of a tree and differences between organs in the relative importance of transport, safety, and mechanical demands. The applied root model was capable of depicting the mass, length and spread of coarse roots of B. pendula and P. abies, and to the lesser extent in P. sylvestris. The roots did not follow self-similar fractal branching, because the parameter values varied within the root systems. Model parameters indicate differences in rooting behavior, and therefore different ecophysiological adaptations between species.

Automatic Lameness Detection in a Milking Robot : Instrumentation, measurement software, algorithms for data analysis and a neural network model

The aim of this thesis is to develop a fully automatic lameness detection system that operates in a milking robot. The instrumentation, measurement software, algorithms for data analysis and a neural network model for lameness detection were developed. Automatic milking has become a common practice in dairy husbandry, and in the year 2006 about 4000 farms worldwide used over 6000 milking robots. There is a worldwide movement with the objective of fully automating every process from feeding to milking. Increase in automation is a consequence of increasing farm sizes, the demand for more efficient production and the growth of labour costs. As the level of automation increases, the time that the cattle keeper uses for monitoring animals often decreases. This has created a need for systems for automatically monitoring the health of farm animals. The popularity of milking robots also offers a new and unique possibility to monitor animals in a single confined space up to four times daily. Lameness is a crucial welfare issue in the modern dairy industry. Limb disorders cause serious welfare, health and economic problems especially in loose housing of cattle. Lameness causes losses in milk production and leads to early culling of animals. These costs could be reduced with early identification and treatment. At present, only a few methods for automatically detecting lameness have been developed, and the most common methods used for lameness detection and assessment are various visual locomotion scoring systems. The problem with locomotion scoring is that it needs experience to be conducted properly, it is labour intensive as an on-farm method and the results are subjective. A four balance system for measuring the leg load distribution of dairy cows during milking in order to detect lameness was developed and set up in the University of Helsinki Research farm Suitia. The leg weights of 73 cows were successfully recorded during almost 10,000 robotic milkings over a period of 5 months. The cows were locomotion scored weekly, and the lame cows were inspected clinically for hoof lesions. Unsuccessful measurements, caused by cows standing outside the balances, were removed from the data with a special algorithm, and the mean leg loads and the number of kicks during milking was calculated. In order to develop an expert system to automatically detect lameness cases, a model was needed. A probabilistic neural network (PNN) classifier model was chosen for the task. The data was divided in two parts and 5,074 measurements from 37 cows were used to train the model. The operation of the model was evaluated for its ability to detect lameness in the validating dataset, which had 4,868 measurements from 36 cows. The model was able to classify 96% of the measurements correctly as sound or lame cows, and 100% of the lameness cases in the validation data were identified. The number of measurements causing false alarms was 1.1%. The developed model has the potential to be used for on-farm decision support and can be used in a real-time lameness monitoring system.

Analysis of Heterogeneous IP Traffic in Wireless Environment

Wireless access is expected to play a crucial role in the future of the Internet. The demands of the wireless environment are not always compatible with the assumptions that were made on the era of the wired links. At the same time, new services that take advantage of the advances in many areas of technology are invented. These services include delivery of mass media like television and radio, Internet phone calls, and video conferencing. The network must be able to deliver these services with acceptable performance and quality to the end user. This thesis presents an experimental study to measure the performance of bulk data TCP transfers, streaming audio flows, and HTTP transfers which compete the limited bandwidth of the GPRS/UMTS-like wireless link. The wireless link characteristics are modeled with a wireless network emulator. We analyze how different competing workload types behave with regular TPC and how the active queue management, the Differentiated services (DiffServ), and a combination of TCP enhancements affect the performance and the quality of service. We test on four link types including an error-free link and the links with different Automatic Repeat reQuest (ARQ) persistency. The analysis consists of comparing the resulting performance in different configurations based on defined metrics. We observed that DiffServ and Random Early Detection (RED) with Explicit Congestion Notification (ECN) are useful, and in some conditions necessary, for quality of service and fairness because a long queuing delay and congestion related packet losses cause problems without DiffServ and RED. However, we observed situations, where there is still room for significant improvements if the link-level is aware of the quality of service. Only very error-prone link diminishes the benefits to nil. The combination of TCP enhancements improves performance. These include initial window of four, Control Block Interdependence (CBI) and Forward RTO recovery (F-RTO). The initial window of four helps a later starting TCP flow to start faster but generates congestion under some conditions. CBI prevents slow-start overshoot and balances slow start in the presence of error drops, and F-RTO reduces unnecessary retransmissions successfully.

Ecological consequences of genetic modifications - an invasion analysis approach

Biological invasions are considered as one of the greatest threats to biodiversity, as they may lead to disruption and homogenization of natural communities, and in the worst case, to native species extinctions. The introduction of gene modified organisms (GMOs) to agricultural, fisheries and forestry practices brings them into contact with natural populations. GMOs may appear as new invasive species if they are able to (1) invade into natural habitats or (2) hybridize with their wild relatives. The benefits of GMOs, such as increased yield or decreased use of insecticides or herbicides in cultivation, may thus be reduced due the potential risks they may cause. A careful ecological risk analysis therefore has to precede any responsible GMO introduction. In this thesis I study ecological invasion in relation to GMOs, and what kind of consequences invasion may have in natural populations. A set of theoretical models that combine life-history evolution, population dynamics, and population genetics were developed for the hazard identification part of ecological risks assessment of GMOs. In addition, the potential benefits of GMOs in management of an invasive pest were analyzed. In the first study I showed that a population that is fluctuating due to scramble-type density dependence (due to, e.g., nutrient competition in plants) may be invaded by a population that is relatively more limited by a resource (e.g., light in plants) that is a cause of contest-type density dependence. This result emphasises the higher risk of invasion in unstable environments. The next two studies focused on escape of a growth hormone (GH) transgenic fish into a natural population. The results showed that previous models may have given too pessimistic a view of the so called Trojan gene -effect, where the invading genotype is harmful for the population as a whole. The previously suggested population extinctions did not occur in my studies, since the changes in mating preferences caused by the GH-fish were be ameliorated by decreased level of competition. The GH-invaders may also have to exceed a threshold density before invasion can be successful. I also showed that the prevalence of mature parr (aka. sneaker) strategy among GH-fish may have clear effect on invasion outcome. The fourth study assessed the risks and developed methods against the invasion of the Colorado Potato Beetle (CPB, Leptinotarsa decemlineata). I showed that the eradication of CPB is most important for the prevention of their establishment, but the cultivation of transgenic Bt-potato could also be effective. In general, my results emphasise that invasion of transgenic species or genotypes to be possible under certain realistic conditions and resulting in competitive exclusion, population decline through outbreeding depression and genotypic displacement of native species. Ecological risk assessment should regard the decline and displacement of the wild genotype by an introduced one as a consequence that is as serious as the population extinction. It will also be crucial to take into account different kinds of behavioural differences among species when assessing the possible hazards that GMOs may cause if escaped. The benefits found of GMO crops effectiveness in pest management may also be too optimistic since CPB may evolve resistance to Bt-toxin. The models in this thesis could be further applied in case specific risk assessment of GMOs by supplementing them with detailed data of the species biology, the effect of the transgene introduced to the species, and also the characteristics of the populations or the environments in the risk of being invaded.

Effects of unilateral vocal fold paralysis treatment with fascia augmentation on airflow mechanics, tracheal sounds and voice acoustics

Data on the influence of unilateral vocal fold paralysis on breathing, especially other than information obtained by spirometry, are relatively scarce. Even less is known about the effect of its treatment by vocal fold medialization. Consequently, there was a need to study the issue by combining multiple instruments capable of assessing airflow dynamics and voice. This need was emphasized by a recently developed medialization technique, autologous fascia injection; its effects on breathing have not previously been investigated. A cohort of ten patients with unilateral vocal fold paralysis was studied before and after autologous fascia injection by using flow-volume spirometry, body plethysmography and acoustic analysis of breathing and voice. Preoperative results were compared with those of ten healthy controls. A second cohort of 11 subjects with unilateral vocal fold paralysis was studied pre- and postoperatively by using flow-volume spirometry, impulse oscillometry, acoustic analysis of voice, voice handicap index and subjective assessment of dyspnoea. Preoperative peak inspiratory flow and specific airway conductance were significantly lower and airway resistance was significantly higher in the patients than in the healthy controls (78% vs. 107%, 73% vs. 116% and 182% vs. 125% of predicted; p = 0.004, p = 0.004 and p = 0.026, respectively). Patients had a higher root mean square of spectral power of tracheal sounds than controls, and three of them had wheezes as opposed to no wheezing in healthy subjects. Autologous fascia injection significantly improved acoustic parameters of the voice in both cohorts and voice handicap index in the latter cohort, indicating that this procedure successfully improved voice in unilateral vocal fold paralysis. Peak inspiratory flow decreased significantly as a consequence of this procedure (from 4.54 ± 1.68 l to 4.21 ± 1.26 l, p = 0.03, in pooled data of both cohorts), but no change occurred in the other variables of flow-volume spirometry, body-plethysmography and impulse oscillometry. Eight of the ten patients studied by acoustic analysis of breathing had wheezes after vocal fold medialization compared with only three patients before the procedure, and the numbers of wheezes per recorded inspirium and expirium increased significantly (from 0.02 to 0.42 and from 0.03 to 0.36; p = 0.028 and p = 0.043, respectively). In conclusion, unilateral vocal fold paralysis was observed to disturb forced breathing and also to cause some signs of disturbed tidal breathing. Findings of flow volume spirometry were consistent with variable extra-thoracic obstruction. Vocal fold medialization by autologous fascia injection improved the quality of the voice in patients with unilateral vocal fold paralysis, but also decreased peak inspiratory flow and induced wheezing during tidal breathing. However, these airflow changes did not appear to cause significant symptoms in patients.

Augustine and the functions of concupiscence

This study analyses Augustine s concept of concupiscentia, or evil desire (together with two cognate terms, libido and cupiditas) in the context of his entire oeuvre. By the aid of systematic analysis, the concept and its development is explored in four distinct ways. It is claimed that Augustine used the concept of concupiscentia for several theological purposes, and the task of the study is to represent these distinct functions, and their connections to Augustine s general theological and philosophical convictions. The study opens with a survey on terminology. A general overview of the occurrences of the negatively connoted words for desire in Latin literature precedes a corresponding examination of Augustine s own works. In this introductory chapter it is shown that, despite certain preferences in the uses of the words, a sufficient degree of synonymy reigns so as to allow an analysis of the concept without tightly discriminating between the terms. The theological functions of concupiscentia with its distinct contexts are analysed in separate chapters. The function of concupiscentia as a divine punishment is explored first (Ch 3). It is seen how Augustine links together concupiscentia and ideas about divine justice, and finally suggests that in the inordinate, psychologically experienced sexual desire, the original theological disobedience of Adam and Eve can be perceived. Augustine was criticized for this solution already in his own times, and the analysis of the function of concupiscentia as a divine punishment ends in a discussion on the critical response of punitive concupiscentia by Julian of Aeclanum. Augustine also attached to concupiscentia another central theological function by viewing evil desire as an inward originating cause for all external evil actions. In the study, this function is analysed by surveying two formally distinct images of evil desire, i.e. as the root (radix) of all evil, and as a threefold (triplex) matrix of evil actions (Ch 4). Both of these images were based on a single verse of the Bible (1 Jn 2, 16 and 1 Tim 6, 10). This function of concupiscentia was formed both parallel to, and in answer to, Manichaean insights into concupiscentia. Being familiar with the traditional philosophical discussions on the nature and therapy of emotions, Augustine situated concupiscentia also into this context. It is acknowledged that these philosophical traditions had an obvious impact into his way of explaining psychological processes in connection with concupiscentia. Not only did Augustine implicitly receive and exploit these traditions, but he also explicitly moulded and criticized them in connection with concupiscentia. Eventually, Augustine conceives the philosophical traditions of emotions as partly useful but also partly inadequate to deal with concupiscentia (Ch 5). The role of concupiscentia in connection to divine grace and Christian renewal is analysed in the final chapter of the study. Augustine s gradual development in internalizing the effects of concupiscentia also into the life of a baptized Christian are elucidated, as are the strong limitations and mitigations Augustine makes to the concept when attaching it into the life under grace (sub gratia). A crucial part in the development of this function is played by Augustine s changing interpretation of Rom 7, and the way concupiscentia appears in Augustine s readings of this text is therefore also analysed. As a result of the analysis of these four distinct functions and contexts of concupiscentia, it is concluded that Augustine s concept of concupiscentia is fairly tightly and coherently connected to his views of central theological importance. Especially the functions of concupiscentia as a punishment and the function of concupiscentia in Christian renewal were both tightly interwoven into Augustine s view of God s being and God s grace. The study shows the importance of reading Augustine s discussions on evil desire with a constant awareness of their role in their larger context, that is, of their function in each situation. The study warns against too simplistic and unifying readings of Augustine s concupiscentia, emphasizing the need to acknowledge both the necessitating, sinful aspects of concupiscentia, and the domesticated features of concupiscentia during Christian renewal.

Transcriptional analysis of persistent Chlamydia pneumoniae infection in vitro

Chlamydia pneumoniae can cause acute respiratory infections including pneumonia. Repeated and persistent Chlamydia infections occur and persistent C. pneumoniae infection may have a role in the pathogenesis of atherosclerosis and coronary heart disease and may also contribute to the development of chronic inflammatory lung diseases like chronic obstructive pulmonary disease (COPD) and asthma. In this thesis in vitro models for persistent C. pneumonia infection were established in epithelial and monocyte/macrophage cell lines. Expression of host cell genes in the persistent C. pneumoniae infection model of epithelial cells was studied by microarray and RT-PCR. In the monocyte/macrophage infection model expression of selected C. pneumoniae genes were studied by RT-PCR and immunofluorescence microscopy. Chlamydia is able to modulate host cell gene expression and apoptosis of host cells, which may assist Chlamydia to evade the host cells' immune responses. This, in turn, may lead to extended survival of the organism inside epithelial cells and promote the development of persistent infection. To simulate persistent C. pneumoniae infection in vivo, we set up a persistent infection model exposing the HL cell cultures to IFN-gamma. When HL cell cultures were treated with moderate concentration of IFN-gamma, the replication of C. pneumoniae DNA was unaffected while differentiation into infectious elementary bodies (EB) was strongly inhibited. By transmission electron microscopy small atypical inclusions were identified in IFN-gamma treated cultures. No second cycle of infection was observed in cells exposed to IFN-gamma , whereas C. pneumoniae was able to undergo a second cycle of infection in unexposed HL cells. Although monocytic cells can naturally restrict chlamydial growth, IFN-gamma further reduced production of infectious C. pneumoniae in Mono Mac 6 cells. Under both studied conditions no second cycle of infection could be detected in monocytic cell line suggesting persistent infection in these cells. As a step toward understanding the role of host genes in the development and pathogenesis of persistent C. pneumoniae infection, modulation of host cell gene expression during IFN-gamma induced persistent infection was examined and compared to that seen during active C. pneumoniae infection or IFN-gamma treatment. Total RNA was collected at 6 to 150 h after infection of an epithelial cell line (HL) and analyzed by a cDNA array (available at that time) representing approximately 4000 human transcripts. In initial analysis 250 of the 4000 genes were identified as differentially expressed upon active and persistent chlamydial infection and IFN-gamma treatment. In persistent infection more potent up-regulation of many genes was observed in IFN-gamma induced persistent infection than in active infection or in IFN-gamma treated cell cultures. Also sustained up-regulation was observed for some genes. In addition, we could identify nine host cell genes whose transcription was specifically altered during the IFN-gamma induced persistent C. pneumoniae infection. Strongest up-regulation in persistent infection in relation to controls was identified for insulin like growth factor binding protein 6, interferon-stimulated protein 15 kDa, cyclin D1 and interleukin 7 receptor. These results suggest that during persistent infection, C. pneumoniae reprograms the host transcriptional machinery regulating a variety of cellular processes including adhesion, cell cycle regulation, growth and inflammatory response, all of which may play important roles in the pathogenesis of persistent C. pneumoniae infection. C. pneumoniae DNA can be detected in peripheral blood mononuclear cells indicating that the bacterium can also infect monocytic cells in vivo and thereby monocytes can assist the spread of infection from the lungs to other anatomical sites. Persistent infection established at these sites could promote inflammation and enhance pathology. Thus, the mononuclear cells are in a strategic position in the development of persistent infection. To investigate the intracellular replication and fate of C. pneumoniae in mononuclear cells we analyzed the transcription of 11 C. pneumoniae genes in Mono Mac 6 cells during infection by real time RT-PCR. Our results suggest that the transcriptional profile of the studied genes in monocytes is different from that seen in epithelial cells and that IFN-gamma has a less significant effect on C. pneumoniae transcription in monocytes. Furthermore, our study shows that type III secretion system (T3SS) related genes are transcribed and that Chlamydia possesses a functional T3SS during infection in monocytes. Since C. pneumoniae infection in monocytes has been implicated to have reduced antibiotic susceptibility, this creates opportunities for novel therapeutics targeting T3SS in the management of chlamydial infection in monocytes.

Analysis of forest soil chemistry and hydrology with a dynamic model ACIDIC.

In this study we analyze how the ion concentrations in forest soil solution are determined by hydrological and biogeochemical processes. A dynamic model ACIDIC was developed, including processes common to dynamic soil acidification models. The model treats up to eight interacting layers and simulates soil hydrology, transpiration, root water and nutrient uptake, cation exchange, dissolution and reactions of Al hydroxides in solution, and the formation of carbonic acid and its dissociation products. It includes also a possibility to a simultaneous use of preferential and matrix flow paths, enabling the throughfall water to enter the deeper soil layers in macropores without first reacting with the upper layers. Three different combinations of routing the throughfall water via macro- and micropores through the soil profile is presented. The large vertical gradient in the observed total charge was simulated succesfully. According to the simulations, gradient is mostly caused by differences in the intensity of water uptake, sulfate adsorption and organic anion retention at the various depths. The temporal variations in Ca and Mg concentrations were simulated fairly well in all soil layers. For H+, Al and K there were much more variation in the observed than in the simulated concentrations. Flow in macropores is a possible explanation for the apparent disequilibrium of the cation exchange for H+ and K, as the solution H+ and K concentrations have great vertical gradients in soil. The amount of exchangeable H+ increased in the O and E horizons and decreased in the Bs1 and Bs2 horizons, the net change in whole soil profile being a decrease. A large part of the decrease of the exchangeable H+ in the illuvial B horizon was caused by sulfate adsorption. The model produces soil water amounts and solution ion concentrations which are comparable to the measured values, and it can be used in both hydrological and chemical studies of soils.

Immunochemical analysis of prolamins in gluten-free foods

People with coeliac disease have to maintain a gluten-free diet, which means excluding wheat, barley and rye prolamin proteins from their diet. Immunochemical methods are used to analyse the harmful proteins and to control the purity of gluten-free foods. In this thesis, the behaviour of prolamins in immunological gluten assays and with different prolamin-specific antibodies was examined. The immunoassays were also used to detect residual rye prolamins in sourdough systems after enzymatic hydrolysis and wheat prolamins after deamidation. The aim was to characterize the ability of the gluten analysis assays to quantify different prolamins in varying matrices in order to improve the accuracy of the assays. Prolamin groups of cereals consist of a complex mixture of proteins that vary in their size and amino acid sequences. Two common characteristics distinguish prolamins from other cereal proteins. Firstly, they are soluble in aqueous alcohols, and secondly, most of the prolamins are mainly formed from repetitive amino acid sequences containing high amounts of proline and glutamine. The diversity among prolamin proteins sets high requirements for their quantification. In the present study, prolamin contents were evaluated using enzyme-linked immunosorbent assays based on ω- and R5 antibodies. In addition, assays based on A1 and G12 antibodies were used to examine the effect of deamidation on prolamin proteins. The prolamin compositions and the cross-reactivity of antibodies with prolamin groups were evaluated with electrophoretic separation and Western blotting. The results of this thesis research demonstrate that the currently used gluten analysis methods are not able to accurately quantify barley prolamins, especially when hydrolysed or mixed in oats. However, more precise results can be obtained when the standard more closely matches the sample proteins, as demonstrated with barley prolamin standards. The study also revealed that all of the harmful prolamins, i.e. wheat, barley and rye prolamins, are most efficiently extracted with 40% 1-propanol containing 1% dithiothreitol at 50 °C. The extractability of barley and rye prolamins was considerably higher with 40% 1-propanol than with 60% ethanol, which is typically used for prolamin extraction. The prolamin levels of rye were lowered by 99.5% from the original levels when an enzyme-active rye-malt sourdough system was used for prolamin degradation. Such extensive degradation of rye prolamins suggest the use of sourdough as a part of gluten-free baking. Deamidation increases the diversity of prolamins and improves their solubility and ability to form structures such as emulsions and foams. Deamidation changes the protein structure, which has consequences for antibody recognition in gluten analysis. According to the resuts of the present work, the analysis methods were not able to quantify wheat gluten after deamidation except at very high concentrations. Consequently, deamidated gluten peptides can exist in food products and remain undetected, and thus cause a risk for people with gluten intolerance. The results of this thesis demonstrate that current gluten analysis methods cannot accurately quantify prolamins in all food matrices. New information on the prolamins of rye and barley in addition to wheat prolamins is also provided in this thesis, which is essential for improving gluten analysis methods so that they can more accurately quantify prolamins from harmful cereals.