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.

Optimum Design for External Seismic Stability of Geosynthetic Reinforced Soil Walls: Reliability Based Approach

In this paper, an analytical study considering the effect of uncertainties in the seismic analysis of geosynthetic-reinforced soil (GRS) walls is presented. Using limit equilibrium method and assuming sliding wedge failure mechanism, analysis is conducted to evaluate the external stability of GRS walls when subjected to earthquake loads. Target reliability based approach is used to estimate the probability of failure in three modes of failure, viz., sliding, bearing, and eccentricity failure. The properties of reinforced backfill, retained backfill, foundation soil, and geosynthetic reinforcement are treated as random variables. In addition, the uncertainties associated with horizontal seismic acceleration and surcharge load acting on the wall are considered. The optimum length of reinforcement needed to maintain the stability against three modes of failure by targeting various component and system reliability indices is obtained. Studies have also been made to study the influence of various parameters on the seismic stability in three failure modes. The results are compared with those given by first-order second moment method and Monte Carlo simulation methods. In the illustrative example, external stability of the two walls, Gould and Valencia walls, subjected to Northridge earthquake is reexamined.

Nonlinear vibration analysis of composite laminated and sandwich plates with random material properties

Nonlinear vibration analysis is performed using a C-0 assumed strain interpolated finite element plate model based on Reddy's third order theory. An earlier model is modified to include the effect of transverse shear variation along the plate thickness and Von-Karman nonlinear strain terms. Monte Carlo Simulation with Latin Hypercube Sampling technique is used to obtain the variance of linear and nonlinear natural frequencies of the plate due to randomness in its material properties. Numerical results are obtained for composite plates with different aspect ratio, stacking sequence and oscillation amplitude ratio. The numerical results are validated with the available literature. It is found that the nonlinear frequencies show increasing non-Gaussian probability density function with increasing amplitude of vibration and show dual peaks at high amplitude ratios. This chaotic nature of the dispersion of nonlinear eigenvalues is also r

Optimum Design for External Seismic Stability of Geosynthetic Reinforced Soil Walls: Reliability Based Approach

In this paper, an analytical study considering the effect of uncertainties in the seismic analysis of geosynthetic-reinforced soil (GRS) walls is presented. Using limit equilibrium method and assuming sliding wedge failure mechanism, analysis is conducted to evaluate the external stability of GRS walls when subjected to earthquake loads. Target reliability based approach is used to estimate the probability of failure in three modes of failure, viz., sliding, bearing, and eccentricity failure. The properties of reinforced backfill, retained backfill, foundation soil, and geosynthetic reinforcement are treated as random variables. In addition, the uncertainties associated with horizontal seismic acceleration and surcharge load acting on the wall are considered. The optimum length of reinforcement needed to maintain the stability against three modes of failure by targeting various component and system reliability indices is obtained. Studies have also been made to study the influence of various parameters on the seismic stability in three failure modes. The results are compared with those given by first-order second moment method and Monte Carlo simulation methods. In the illustrative example, external stability of the two walls, Gould and Valencia walls, subjected to Northridge earthquake is reexamined.

Essays on Risk Modeling. Applications to Portfolio and Risk Management (summary section only)

In this thesis we deal with the concept of risk. The objective is to bring together and conclude on some normative information regarding quantitative portfolio management and risk assessment. The first essay concentrates on return dependency. We propose an algorithm for classifying markets into rising and falling. Given the algorithm, we derive a statistic: the Trend Switch Probability, for detection of long-term return dependency in the first moment. The empirical results suggest that the Trend Switch Probability is robust over various volatility specifications. The serial dependency in bear and bull markets behaves however differently. It is strongly positive in rising market whereas in bear markets it is closer to a random walk. Realized volatility, a technique for estimating volatility from high frequency data, is investigated in essays two and three. In the second essay we find, when measuring realized variance on a set of German stocks, that the second moment dependency structure is highly unstable and changes randomly. Results also suggest that volatility is non-stationary from time to time. In the third essay we examine the impact from market microstructure on the error between estimated realized volatility and the volatility of the underlying process. With simulation-based techniques we show that autocorrelation in returns leads to biased variance estimates and that lower sampling frequency and non-constant volatility increases the error variation between the estimated variance and the variance of the underlying process. From these essays we can conclude that volatility is not easily estimated, even from high frequency data. It is neither very well behaved in terms of stability nor dependency over time. Based on these observations, we would recommend the use of simple, transparent methods that are likely to be more robust over differing volatility regimes than models with a complex parameter universe. In analyzing long-term return dependency in the first moment we find that the Trend Switch Probability is a robust estimator. This is an interesting area for further research, with important implications for active asset allocation.

Bootstrap and Fast Double Bootstrap Tests of Cointegration Rank with Financial Time Series

The likelihood ratio test of cointegration rank is the most widely used test for cointegration. Many studies have shown that its finite sample distribution is not well approximated by the limiting distribution. The article introduces and evaluates by Monte Carlo simulation experiments bootstrap and fast double bootstrap (FDB) algorithms for the likelihood ratio test. It finds that the performance of the bootstrap test is very good. The more sophisticated FDB produces a further improvement in cases where the performance of the asymptotic test is very unsatisfactory and the ordinary bootstrap does not work as well as it might. Furthermore, the Monte Carlo simulations provide a number of guidelines on when the bootstrap and FDB tests can be expected to work well. Finally, the tests are applied to US interest rates and international stock prices series. It is found that the asymptotic test tends to overestimate the cointegration rank, while the bootstrap and FDB tests choose the correct cointegration rank.

Markov random fields and spatial smoothing in improving of forest inventory estimates

Markov random fields (MRF) are popular in image processing applications to describe spatial dependencies between image units. Here, we take a look at the theory and the models of MRFs with an application to improve forest inventory estimates. Typically, autocorrelation between study units is a nuisance in statistical inference, but we take an advantage of the dependencies to smooth noisy measurements by borrowing information from the neighbouring units. We build a stochastic spatial model, which we estimate with a Markov chain Monte Carlo simulation method. The smooth values are validated against another data set increasing our confidence that the estimates are more accurate than the originals.

Spectral element approach to wave propagation in uncertain beam structures

This paper presents a study on the uncertainty in material parameters of wave propagation responses in metallic beam structures. Special effort is made to quantify the effect of uncertainty in the wave propagation responses at high frequencies. Both the modulus of elasticity and the density are considered uncertain. The analysis is performed using a Monte Carlo simulation (MCS) under the spectral finite element method (SEM). The randomness in the material properties is characterized by three different distributions, the normal, Weibull and extreme value distributions. Their effect on wave propagation in beams is investigated. The numerical study shows that the CPU time taken for MCS under SEM is about 48 times less than for MCS under a conventional one-dimensional finite element environment for 50 kHz loading. The numerical results presented investigate effects of material uncertainties on high frequency modes. A study is performed on the usage of different beam theories and their uncertain responses due to dynamic impulse load. These studies show that even for a small coefficient of variation, significant changes in the above parameters are noticed. A number of interesting results are presented, showing the true effects of uncertainty response due to dynamic impulse load.

Hysteresis in model spin system

The aim of this paper is to construct a nonequilibrium statistical‐mechanics theory to study hysteresis in ferromagnetic systems. We study the hysteretic response of model spin systems to periodic magnetic fields H(t) as a function of the amplitude H0 and frequency Ω. At fixed H0, we find conventional, squarelike hysteresis loops at low Ω, and rounded, roughly elliptical loops at high Ω, in agreement with experiments. For the O(N→∞), d=3, (Φ2)2 model with Langevin dynamics, we find a novel scaling behavior for the area A of the hysteresis loop, of the form A∝H0.660Ω0.33. We carry out a Monte Carlo simulation of the hysteretic response of the two‐dimensional, nearest‐neighbor, ferromagnetic Ising model. These results agree qualitatively with the results obtained for the O(N) model.

Effect of uncertainty on helicopter performance predictions

The effect of uncertainties on performance predictions of a helicopter is studied in this article. The aeroelastic parameters such as the air density, blade profile drag coefficient, main rotor angular velocity, main rotor radius, and blade chord are considered as uncertain variables. The propagation of these uncertainties in the performance parameters such as thrust coefficient, figure of merit, induced velocity, and power required are studied using Monte Carlo simulation and the first-order reliability method. The Rankine-Froude momentum theory is used for performance prediction in hover, axial climb, and forward flight. The propagation of uncertainty causes large deviations from the baseline deterministic predictions, which undoubtedly affect both the achievable performance and the safety of the helicopter. The numerical results in this article provide useful bounds on helicopter power requirements.

Free vibration analysis Of stochastic strings

The free vibration of strings with randomly varying mass and stiffness is considered. The joint probability density functions of the eigenvalues and eigenfunctions are characterized in terms of the solution of a pair of stochastic non-linear initial value problems. Analytical solutions of these equations based on the method of stochastic averaging are obtained. The effects of the mean and autocorrelation of the mass process are included in the analysis. Numerical results for the marginal probability density functions of eigenvalues and eigenfunctions are obtained and are found to compare well with Monte Carlo simulation results. The random eigenvalues, when normalized with respect to their corresponding deterministic values, are observed to tend to become first order stochastically stationary with respect to the mode count.

Composite material and piezoelectric coefficient uncertainty effects on structural health monitoring using feedback control gains as damage indicators

This article addresses uncertainty effect on the health monitoring of a smart structure using control gain shifts as damage indicators. A finite element model of the smart composite plate with surface-bonded piezoelectric sensors and actuators is formulated using first-order shear deformation theory and a matrix crack model is integrated into the finite element model. A constant gain velocity/position feedback control algorithm is used to provide active damping to the structure. Numerical results show that the response of the structure is changed due to matrix cracks and this change can be compensated by actively tuning the feedback controller. This change in control gain can be used as a damage indicator for structural health monitoring. Monte Carlo simulation is conducted to study the effect of material uncertainty on the damage indicator by considering composite material properties and piezoelectric coefficients as independent random variables. It is found that the change in position feedback control gain is a robust damage indicator.

Nonequilibrium phase transition in the kinetic Ising model: Critical slowing down and the specific-heat singularity

The nonequilibrium dynamic phase transition, in the kinetic Ising model in the presence of an oscillating magnetic field has been studied both by Monte Carlo simulation and by solving numerically the mean-field dynamic equation of motion for the average magnetization. In both cases, the Debye ''relaxation'' behavior of the dynamic order parameter has been observed and the ''relaxation time'' is found to diverge near the dynamic transition point. The Debye relaxation of the dynamic order parameter and the power law divergence of the relaxation time have been obtained from a very approximate solution of the mean-field dynamic equation. The temperature variation of appropriately defined ''specific heat'' is studied by the Monte Carlo simulation near the transition point. The specific heat has been observed to diverge near the dynamic transition point.

Nonequilibrium phase transition in the kinetic Ising model: Divergences of fluctuations and responses near the transition point

The nonequilibrium dynamic phase transition in the kinetic Ising model in the presence of an oscillating magnetic field is studied by Monte Carlo simulation. The fluctuation of the dynamic older parameter is studied as a function of temperature near the dynamic transition point. The temperature variation of appropriately defined ''susceptibility'' is also studied near the dynamic transition point. Similarly, the fluctuation of energy and appropriately defined ''specific heat'' is studied as a function of temperature near the dynamic transition point. In both cases, the fluctuations (of dynamic order parameter and energy) and the corresponding responses diverge (in power law fashion) near the dynamic transition point with similar critical behavior (with identical exponent values).

Analysis of LDPC Codes for Compression of Nonuniform Sources with Side Information Using Density Evolution

In this paper, we explore the use of LDPC codes for nonuniform sources under distributed source coding paradigm. Our analysis reveals that several capacity approaching LDPC codes indeed do approach the Slepian-Wolf bound for nonuniform sources as well. The Monte Carlo simulation results show that highly biased sources can be compressed to 0.049 bits/sample away from Slepian-Wolf bound for moderate block lengths.