Confidence Level Computation for Combining Searches with Small Statistics

This article describes an efficient procedure for computing approximate confidence levels for searches for new particles where the expected signal and background levels are small enough to require the use of Poisson statistics. The results of many independent searches for the same particle may be combined easily, regardless of the discriminating variables which may be measured for the candidate events. The effects of systematic uncertainty in the signal and background models are incorporated in the confidence levels. The procedure described allows efficient computation of expected confidence levels.Comment: 17 pages, 2 figures, LaTeX, program available at http://home.cern.ch/~thomasj/searchlimits/ecl.html one replacement to fix a typo in a UR

Recent Heavy-Flavor Measurements from OPAL

A selection of recent heavy-flavor results from OPAL using the LEP1 data sample are presented. The average polarization of b baryons in hadronic Z^0 decay has been measured to be -0.56^{+0.20}_{-0.13} (stat.) +- 0.09(syst.) using semileptonic decays of Lambda_b baryons. A search has been conducted for the radially excited D*' and has produced a 95% CL upper limit on its production of f(Z^0 -> D*'+-(2629))xBr(D*'+- -> D*+- pi+ pi-) < 2.1x10^{-3}. Finally, the measurement of the product branching ratio $f(b -> Lambda_b)xBr(Lambda_b -> Lambda X)= (2.67 +- 0.38 (stat) ^{+0.67}_{-0.60}(syst.))% has been made. This measurement, along with an earlier measurement of the product branching ratio f(b -> Lambda_b)xBr(Lambda_b -> Lambda l X), has been used to compute an updated R_{Lambda l} = Br(Lambda_b -> Lambda l X)/Br(Lambda_b -> \Lambda X)= (8.0 +- 1.2 (stat.) +- 0.9 (syst.))%, consistent with the expected low semileptonic branching fraction of the Lambda_b inferred from its short lifetime compared to the other b hadrons.Comment: 9 Pages, LaTeX, 4 figures. Proceedings contributed to DPF'99, UCLA, Parallel session talk 3.0

Optimal low-dispersion low-dissipation LBM schemes for computational aeroacoustics

Lattice Boltmzmann Methods (LBM) have been proved to be very effective methods for computational aeroacoustics (CAA), which have been used to capture the dynamics of weak acoustic fluctuations. In this paper, we propose a strategy to reduce the dispersive and disspative errors of the two-dimensional (2D) multi-relaxation-time lattice Boltzmann method (MRT-LBM). By presenting an effective algorithm, we obtain a uniform form of the linearized Navier-Stokes equations corresponding to the MRT-LBM in wave-number space. Using the matrix perturbation theory and the equivalent modified equation approach for finite difference methods, we propose a class of minimization problems to optimize the free-parameters in the MRT-LBM. We obtain this way a dispersion-relation-preserving LBM (DRP-LBM) to circumvent the minimized dispersion error of the MRT-LBM. The dissipation relation precision is also improved.And the stability of the MRT-LBM with the small bulk viscosity is guaranteed. Von Neuman analysis of the linearized MRT-LBM is performed to validate the optimized dispersion/dissipation relations considering monochromatic wave solutions. Meanwhile, dispersion and dissipation errors of the optimized MRT-LBM are quantitatively compared with the original MRT-LBM . Finally, some numerical simulations are carried out to assess the new optimized MRT-LBM schemes.Comment: 33 page

Modelling Shear Flows with SPH and Grid Based Methods

Given the importance of shear flows for astrophysical gas dynamics, we study the evolution of the Kelvin-Helmholtz instability (KHI) analytically and numerically. We derive the dispersion relation for the two-dimensional KHI including viscous dissipation. The resulting expression for the growth rate is then used to estimate the intrinsic viscosity of four numerical schemes depending on code-specific as well as on physical parameters. Our set of numerical schemes includes the Tree-SPH code VINE, an alternative SPH formulation developed by Price (2008), and the finite-volume grid codes FLASH and PLUTO. In the first part, we explicitly demonstrate the effect of dissipation-inhibiting mechanisms such as the Balsara viscosity on the evolution of the KHI. With VINE, increasing density contrasts lead to a continuously increasing suppression of the KHI (with complete suppression from a contrast of 6:1 or higher). The alternative SPH formulation including an artificial thermal conductivity reproduces the analytically expected growth rates up to a density contrast of 10:1. The second part addresses the shear flow evolution with FLASH and PLUTO. Both codes result in a consistent non-viscous evolution (in the equal as well as in the different density case) in agreement with the analytical prediction. The viscous evolution studied with FLASH shows minor deviations from the analytical prediction.Comment: 16 pages, 17 figure

Reproducibility and Replication of Experimental Particle Physics Results

Recently, much attention has been focused on the replicability of scientific results, causing scientists, statisticians, and journal editors to examine closely their methodologies and publishing criteria. Experimental particle physicists have been aware of the precursors of non-replicable research for many decades and have many safeguards to ensure that the published results are as reliable as possible. The experiments require large investments of time and effort to design, construct, and operate. Large collaborations produce and check the results, and many papers are signed by more than three thousand authors. This paper gives an introduction to what experimental particle physics is and to some of the tools that are used to analyze the data. It describes the procedures used to ensure that results can be computationally reproduced, both by collaborators and by non-collaborators. It describes the status of publicly available data sets and analysis tools that aid in reproduction and recasting of experimental results. It also describes methods particle physicists use to maximize the reliability of the results, which increases the probability that they can be replicated by other collaborations or even the same collaborations with more data and new personnel. Examples of results that were later found to be false are given, both with failed replication attempts and one with alarmingly successful replications. While some of the characteristics of particle physics experiments are unique, many of the procedures and techniques can be and are used in other fields.Comment: 50 pages, 6 figure

On the interaction of a thin, supersonic shell with a molecular cloud

Molecular clouds (MCs) are stellar nurseries, however, formation of stars within MCs depends on the ambient physical conditions. MCs, over a free-fall time are exposed to numerous dynamical phenomena, of which, the interaction with a thin, dense shell of gas is but one. Below we present results from self-gravitating, 3-D smoothed particle hydrodynamics ({\small SPH}) simulations of the problem; seven realisations of the problem have been performed by varying the precollision density within the cloud, the nature of the post-collision shock, and the spatial resolution in the computational domain. Irrespective of the type of shock, a complex network of dense filaments, seeded by numerical noise, readily appears in the shocked cloud. Segregation of the dense and rarefied gas phases also manifests itself in a bimodal distribution of gas density. We demonstrate that the power-spectrum for rarefied gas is Kolomogorov like, while that for the denser gas is considerably steeper. As a corollary to the main problem, we also look into the possibly degenerative effect of the {\small SPH} artificial viscosity on the impact of the incident shell. It is observed that stronger viscosity leads to greater post-shock dissipation, that strongly decelerates the incident shock-front and promotes formation of contiguous structure, albeit on a much longer timescale. We conclude that too much viscosity is likely to enhance the proclivity towards gravitational boundedness of structure, leading to unphysical fragmentation.On the other hand, insufficient resolution appears to suppress fragmentation. Convergence of results is tested at both extremes, first by repeating the test case with more than a million particles and then with only half the number of particles in the original test case.Comment: 15 pages, 15 figures, and 1 Table; To appear in Monthly Notices to the RA