Source author record

Jan Conrad

Jan Conrad appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

30 published item(s)

preprint2022arXiv

Simple and statistically sound recommendations for analysing physical theories

Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by multiple experiments, and random or grid sampling of model parameters. Whilst these methods are easy to apply, they exhibit pathologies even in low-dimensional parameter spaces, and quickly become problematic to use and interpret in higher dimensions. In this article we give clear guidance for going beyond these procedures, suggesting where possible simple methods for performing statistically sound inference, and recommendations of readily-available software tools and standards that can assist in doing so. Our aim is to provide any physicists lacking comprehensive statistical training with recommendations for reaching correct scientific conclusions, with only a modest increase in analysis burden. Our examples can be reproduced with the code publicly available at https://doi.org/10.5281/zenodo.4322283.

preprint2021arXiv

Publishing statistical models: Getting the most out of particle physics experiments

The statistical models used to derive the results of experimental analyses are of incredible scientific value and are essential information for analysis preservation and reuse. In this paper, we make the scientific case for systematically publishing the full statistical models and discuss the technical developments that make this practical. By means of a variety of physics cases -- including parton distribution functions, Higgs boson measurements, effective field theory interpretations, direct searches for new physics, heavy flavor physics, direct dark matter detection, world averages, and beyond the Standard Model global fits -- we illustrate how detailed information on the statistical modelling can enhance the short- and long-term impact of experimental results.

preprint2019arXiv

Searching for new physics with profile likelihoods: Wilks and beyond

Particle physics experiments use likelihood ratio tests extensively to compare hypotheses and to construct confidence intervals. Often, the null distribution of the likelihood ratio test statistic is approximated by a $χ^2$ distribution, following a theorem due to Wilks. However, many circumstances relevant to modern experiments can cause this theorem to fail. In this paper, we review how to identify these situations and construct valid inference.

preprint2018arXiv

A global fit of the MSSM with GAMBIT

We study the seven-dimensional Minimal Supersymmetric Standard Model (MSSM7) with the new GAMBIT software framework, with all parameters defined at the weak scale. Our analysis significantly extends previous weak-scale, phenomenological MSSM fits, by adding more and newer experimental analyses, improving the accuracy and detail of theoretical predictions, including dominant uncertainties from the Standard Model, the Galactic dark matter halo and the quark content of the nucleon, and employing novel and highly-efficient statistical sampling methods to scan the parameter space. We find regions of the MSSM7 that exhibit co-annihilation of neutralinos with charginos, stops and sbottoms, as well as models that undergo resonant annihilation via both light and heavy Higgs funnels. We find high-likelihood models with light charginos, stops and sbottoms that have the potential to be within the future reach of the LHC. Large parts of our preferred parameter regions will also be accessible to the next generation of direct and indirect dark matter searches, making prospects for discovery in the near future rather good.

preprint2018arXiv

DarkBit: A GAMBIT module for computing dark matter observables and likelihoods

We introduce DarkBit, an advanced software code for computing dark matter constraints on various extensions to the Standard Model of particle physics, comprising both new native code and interfaces to external packages. This release includes a dedicated signal yield calculator for gamma-ray observations, which significantly extends current tools by implementing a cascade decay Monte Carlo, as well as a dedicated likelihood calculator for current and future experiments (gamlike). This provides a general solution for studying complex particle physics models that predict dark matter annihilation to a multitude of final states. We also supply a direct detection package that models a large range of direct detection experiments (DDcalc), and provides the corresponding likelihoods for arbitrary combinations of spin-independent and spin-dependent scattering processes. Finally, we provide custom relic density routines along with interfaces to DarkSUSY, micrOMEGAs, and the neutrino telescope likelihood package nuLike. DarkBit is written in the framework of the Global And Modular Beyond the Standard Model Inference Tool (GAMBIT), providing seamless integration into a comprehensive statistical fitting framework that allows users to explore new models with both particle and astrophysics constraints, and a consistent treatment of systematic uncertainties. In this paper we describe its main functionality, provide a guide to getting started quickly, and show illustrative examples for results obtained with DarkBit (both as a standalone tool and as a GAMBIT module). This includes a quantitative comparison between two of the main dark matter codes (DarkSUSY and micrOMEGAs), and application of DarkBit's advanced direct and indirect detection routines to a simple effective dark matter model.

preprint2018arXiv

Global fits of GUT-scale SUSY models with GAMBIT

We present the most comprehensive global fits to date of three supersymmetric models motivated by grand unification: the Constrained Minimal Supersymmetric Standard Model (CMSSM), and its Non-Universal Higgs Mass generalisations NUHM1 and NUHM2. We include likelihoods from a number of direct and indirect dark matter searches, a large collection of electroweak precision and flavour observables, direct searches for supersymmetry at LEP and Runs I and II of the LHC, and constraints from Higgs observables. Our analysis improves on existing results not only in terms of the number of included observables, but also in the level of detail with which we treat them, our sampling techniques for scanning the parameter space, and our treatment of nuisance parameters. We show that stau co-annihilation is now ruled out in the CMSSM at more than 95\% confidence. Stop co-annihilation turns out to be one of the most promising mechanisms for achieving an appropriate relic density of dark matter in all three models, whilst avoiding all other constraints. We find high-likelihood regions of parameter space featuring light stops and charginos, making them potentially detectable in the near future at the LHC. We also show that tonne-scale direct detection will play a largely complementary role, probing large parts of the remaining viable parameter space, including essentially all models with multi-TeV neutralinos.

preprint2018arXiv

Status of the scalar singlet dark matter model

One of the simplest viable models for dark matter is an additional neutral scalar, stabilised by a $\mathbb{Z}_2$ symmetry. Using the GAMBIT package and combining results from four independent samplers, we present Bayesian and frequentist global fits of this model. We vary the singlet mass and coupling along with 13 nuisance parameters, including nuclear uncertainties relevant for direct detection, the local dark matter density, and selected quark masses and couplings. We include the dark matter relic density measured by Planck, direct searches with LUX, PandaX, SuperCDMS and XENON100, limits on invisible Higgs decays from the Large Hadron Collider, searches for high-energy neutrinos from dark matter annihilation in the Sun with IceCube, and searches for gamma rays from annihilation in dwarf galaxies with the Fermi-LAT. Viable solutions remain at couplings of order unity, for singlet masses between the Higgs mass and about 300 GeV, and at masses above $\sim$1 TeV. Only in the latter case can the scalar singlet constitute all of dark matter. Frequentist analysis shows that the low-mass resonance region, where the singlet is about half the mass of the Higgs, can also account for all of dark matter, and remains viable. However, Bayesian considerations show this region to be rather fine-tuned.

preprint2017arXiv

Comparison of statistical sampling methods with ScannerBit, the GAMBIT scanning module

We introduce ScannerBit, the statistics and sampling module of the public, open-source global fitting framework GAMBIT. ScannerBit provides a standardised interface to different sampling algorithms, enabling the use and comparison of multiple computational methods for inferring profile likelihoods, Bayesian posteriors, and other statistical quantities. The current version offers random, grid, raster, nested sampling, differential evolution, Markov Chain Monte Carlo (MCMC) and ensemble Monte Carlo samplers. We also announce the release of a new standalone differential evolution sampler, Diver, and describe its design, usage and interface to ScannerBit. We subject Diver and three other samplers (the nested sampler MultiNest, the MCMC GreAT, and the native ScannerBit implementation of the ensemble Monte Carlo algorithm T-Walk) to a battery of statistical tests. For this we use a realistic physical likelihood function, based on the scalar singlet model of dark matter. We examine the performance of each sampler as a function of its adjustable settings, and the dimensionality of the sampling problem. We evaluate performance on four metrics: optimality of the best fit found, completeness in exploring the best-fit region, number of likelihood evaluations, and total runtime. For Bayesian posterior estimation at high resolution, T-Walk provides the most accurate and timely mapping of the full parameter space. For profile likelihood analysis in less than about ten dimensions, we find that Diver and MultiNest score similarly in terms of best fit and speed, outperforming GreAT and T-Walk; in ten or more dimensions, Diver substantially outperforms the other three samplers on all metrics.

preprint2017arXiv

GAMBIT: The Global and Modular Beyond-the-Standard-Model Inference Tool

We describe the open-source global fitting package GAMBIT: the Global And Modular Beyond-the-Standard-Model Inference Tool. GAMBIT combines extensive calculations of observables and likelihoods in particle and astroparticle physics with a hierarchical model database, advanced tools for automatically building analyses of essentially any model, a flexible and powerful system for interfacing to external codes, a suite of different statistical methods and parameter scanning algorithms, and a host of other utilities designed to make scans faster, safer and more easily-extendible than in the past. Here we give a detailed description of the framework, its design and motivation, and the current models and other specific components presently implemented in GAMBIT. Accompanying papers deal with individual modules and present first GAMBIT results. GAMBIT can be downloaded from gambit.hepforge.org.

preprint2016arXiv

A method for comparing non-nested models with application to astrophysical searches for new physics

Searches for unknown physics and decisions between competing astrophysical models to explain data both rely on statistical hypothesis testing. The usual approach in searches for new physical phenomena is based on the statistical Likelihood Ratio Test (LRT) and its asymptotic properties. In the common situation, when neither of the two models under comparison is a special case of the other i.e., when the hypotheses are non-nested, this test is not applicable. In astrophysics, this problem occurs when two models that reside in different parameter spaces are to be compared. An important example is the recently reported excess emission in astrophysical $γ$-rays and the question whether its origin is known astrophysics or dark matter. We develop and study a new, simple, generally applicable, frequentist method and validate its statistical properties using a suite of simulations studies. We exemplify it on realistic simulated data of the Fermi-LAT $γ$-ray satellite, where non-nested hypotheses testing appears in the search for particle dark matter.

preprint2016arXiv

On methods for correcting for the look-elsewhere effect in searches for new physics

The search for new significant peaks over a energy spectrum often involves a statistical multiple hypothesis testing problem. Separate tests of hypothesis are conducted at different locations producing an ensemble of local p-values, the smallest of which is reported as evidence for the new resonance. Unfortunately, controlling the false detection rate (type I error rate) of such procedures may lead to excessively stringent acceptance criteria. In the recent physics literature, two promising statistical tools have been proposed to overcome these limitations. In 2005, a method to "find needles in haystacks" was introduced by Pilla et al. [1], and a second method was later proposed by Gross and Vitells [2] in the context of the "look elsewhere effect" and trial factors. We show that, for relatively small sample sizes, the former leads to an artificial inflation of statistical power that stems from an increase in the false detection rate, whereas the two methods exhibit similar performance for large sample sizes. We apply the methods to realistic simulations of the Fermi Large Area Telescope data, in particular the search for dark matter annihilation lines. Further, we discuss the counter-intutive scenario where the look-elsewhere corrections are more conservative than much more computationally efficient corrections for multiple hypothesis testing. Finally, we provide general guidelines for navigating the tradeoffs between statistical and computational efficiency when selecting a statistical procedure for signal detection.

preprint2016arXiv

Sensitivity of the Cherenkov Telescope Array to the Detection of Intergalactic Magnetic Fields

Very high energy (VHE; energy $E \gtrsim 100\,$GeV) $γ$-rays originating from extragalactic sources undergo pair production with low-energy photons of background radiation fields. These pairs can inverse-Compton-scatter background photons, initiating an electromagnetic cascade. The spatial and temporal structure of this secondary $γ$-ray signal is altered as the $e^+e^-$ pairs are deflected in an intergalactic magnetic field (IGMF). We investigate how VHE observations with the future Cherenkov Telescope Array with its high angular resolution and broad energy range, can potentially probe the IGMF. We identify promising sources and simulate $γ$-ray spectra over a wide range of values of the IGMF strength and coherence length using the publicly available ELMAG Monte Carlo code. Combining simulated observations in a joint likelihood approach, we find that current limits on the IGMF can be significantly improved. The projected sensitivity depends strongly on the time a source has been $γ$-ray active and on the emitted maximum $γ$-ray energy.

preprint2015arXiv

The 1st Fermi Lat Supernova Remnant Catalog

To uniformly determine the properties of supernova remnants (SNRs) at high energies, we have developed the first systematic survey at energies from 1 to 100 GeV using data from the Fermi Large Area Telescope. Based on the spatial overlap of sources detected at GeV energies with SNRs known from radio surveys, we classify 30 sources as likely GeV SNRs. We also report 14 marginal associations and 245 flux upper limits. A mock catalog in which the positions of known remnants are scrambled in Galactic longitude, allows us to determine an upper limit of 22% on the number of GeV candidates falsely identified as SNRs. We have also developed a method to estimate spectral and spatial systematic errors arising from the diffuse interstellar emission model, a key component of all Galactic Fermi LAT analyses. By studying remnants uniformly in aggregate, we measure the GeV properties common to these objects and provide a crucial context for the detailed modeling of individual SNRs. Combining our GeV results with multiwavelength (MW) data, including radio, X-ray, and TeV, demonstrates the need for improvements to previously sufficient, simple models describing the GeV and radio emission from these objects. We model the GeV and MW emission from SNRs in aggregate to constrain their maximal contribution to observed Galactic cosmic rays.

preprint2015arXiv

WIMP searches with gamma rays in the Fermi era: challenges, methods and results

The launch of the gamma-ray telescope Fermi Large Area Telescope (Fermi-LAT) started a pivotal period in indirect detection of dark matter. By outperforming expectations, for the first time a robust and stringent test of the paradigm of weakly interacting massive particles (WIMPs) is within reach. In this paper, we discuss astrophysical targets for WIMP detection and the challenges they present, review the analysis tools which have been employed to tackle these challenges, and summarize the status of constraints on and the claimed detections in the WIMP parameter space. Methods and results will be discussed in comparison to Imaging Air Cherenkov Telescopes. We also provide an outlook on short term and longer term developments.

preprint2014arXiv

Indirect Detection of WIMP Dark Matter: a compact review

Indirect detection of dark matter particles, i.e. the detection of annihilation or decay products of Weakly Interacting Massive Particles, has entered a pivotal phase as experiments reach sensitivities that probe the most interesting parameter space. This period is naturally accompanied by claims of detection. In this contribution I discuss and compare different probes (gamma-rays, neutrinos and charged cosmic rays) and review the status and prospects of constraints and recent detection claims.

preprint2014arXiv

On detecting oscillations of gamma rays into axion-like particles in turbulent and coherent magnetic fields

Background radiation fields pervade the Universe, and above a certain energy any $γ$-ray flux emitted by an extragalactic source should be attenuated due to $e^+e^-$ pair production. The opacity could be alleviated if photons oscillated into hypothetical axion-like particles (ALPs) in ambient magnetic fields, leading to a $γ$-ray excess especially at high optical depths that could be detected with imaging air Cherenkov telescopes (IACTs). Here, we introduce a method to search for such a signal in $γ$-ray data and to estimate sensitivities for future observations. Different magnetic fields close to the $γ$-ray source are taken into account in which photons can convert into ALPs that then propagate unimpeded over cosmological distances until they re-convert in the magnetic field of the Milky Way. Specifically, we consider the coherent field at parsec scales in a blazar jet as well as the turbulent field inside a galaxy cluster. For the latter, we explicitly derive the transversal components of a magnetic field with gaussian turbulence which are responsible for the photon-ALP mixing. To illustrate the method, we apply it to a mock IACT array with characteristics similar to the Cherenkov Telescope Array and investigate the dependence of the sensitivity to detect a $γ$-ray excess on the magnetic-field parameters.

preprint2014arXiv

Statistical Issues in Astrophysical Searches for Particle Dark Matter

In this review statistical issues appearing in astrophysical searches for particle dark matter, i.e. indirect detection (dark matter annihilating into standard model particles) or direct detection (dark matter particles scattering in deep underground detectors) are discussed. One particular aspect of these searches is the presence of very large uncertainties in nuisance parameters (astrophysical factors) that are degenerate with parameters of interest (mass and annihilation/decay cross sections for the particles). The likelihood approach has become the most powerful tool, offering at least one well motivated method for incorporation of nuisance parameters and increasing the sensitivity of experiments by allowing a combination of targets superior to the more traditional data stacking. Other statistical challenges appearing in astrophysical searches are to large extent similar to any new physics search, for example at colliders, a prime example being the calculation of trial factors. Frequentist methods prevail for hypothesis testing and interval estimation, Bayesian methods are used for assessment of nuisance parameters and parameter estimation in complex parameter spaces. The basic statistical concepts will be exposed, illustrated with concrete examples from experimental searches and caveats will be pointed out.

preprint2014arXiv

The sensitivity of Cherenkov telescopes to dark matter and astrophysical anisotropies in the diffuse gamma-ray background

In this article, the capability of present (H.E.S.S., MAGIC, VERITAS) and planned (CTA) ground-based Cherenkov telescope systems for detecting angular anisotropies in the diffuse gamma-ray background is investigated. Following up on a study of the impact of instrumental characteristics (effective area, field of view, angular resolution, and background rejection efficiency), the first part examines the influence of different observational strategies, i.e. whether a single deep observation or a splitting over multiple shallow fields is preferred. In the second part, the sensitivity to anisotropies generated by self-annihilating dark matter is studied for different common dark matter models. We find that a relative contribution of ~10% from dark matter annihilation to the extra-galactic diffuse gamma-ray background can be detected with planned configurations of CTA. In terms of the thermally-averaged self-annihilation cross section, the sensitivity of CTA corresponds to values below the thermal freeze-out expectation <sigma v> = 3 x 10-26 cm3s-1 for dark matter particles lighter than ~200 GeV. We stress the importance of constraining anisotropies from unresolved astrophysical sources with currently operating instruments already, as a novel and complementary method for investigating the properties of TeV sources.

preprint2012arXiv

Fermi-LAT measurement of the diffuse gamma-ray emission and constraints on the Galactic Dark Matter signal

We study diffuse gamma-ray emission at intermediate Galactic latitudes measured by the Fermi Large Area Telescope with the aim of searching for a signal from dark matter annihilation or decay. In the absence of a robust dark matter signal, we set conservative dark matter limits requiring that the dark matter signal does not exceed the observed diffuse gamma-ray emission. A second set of more stringent limits is derived based on modeling the foreground astrophysical diffuse emission. Uncertainties in the height of the diffusive cosmic-ray halo, the distribution of the cosmic-ray sources in the Galaxy, the cosmic-ray electron index of the injection spectrum and the column density of the interstellar gas are taken into account using a profile likelihood formalism, while the parameters governing the cosmic-ray propagation have been derived from fits to local cosmic-ray data. The resulting limits impact the range of particle masses over which dark matter thermal production in the early Universe is possible, and challenge the interpretation of the PAMELA/Fermi-LAT cosmic ray anomalies as annihilation of dark matter.

preprint2012arXiv

Handling Systematic Uncertainties and Combined Source Analyses for Atmospheric Cherenkov Telescopes

In response to an increasing availability of statistically rich observational data sets, the performance and applicability of traditional Atmospheric Cherenkov Telescope analyses in the regime of systematically dominated measurement uncertainties is examined. In particular, the effect of systematic uncertainties affecting the relative normalisation of fiducial ON and OFF-source sampling regions - often denoted as α - is investigated using combined source analysis as a representative example case. The traditional summation of accumulated ON and OFF-source event counts is found to perform sub-optimally in the studied contexts and requires careful calibration to correct for unexpected and potentially misleading statistical behaviour. More specifically, failure to recognise and correct for erroneous estimates of α is found to produce substantial overestimates of the combined population significance which worsen with increasing target multiplicity. An alternative joint likelihood technique is introduced, which is designed to treat systematic uncertainties in a uniform and statistically robust manner. This alternate method is shown to yield dramatically enhanced performance and reliability with respect to the more traditional approach.

preprint2012arXiv

Investigating Gamma-Ray Lines from Dark Matter with Future Observatories

We study the prospects for studying line features in gamma-ray spectra with upcoming gamma-ray experiments, such as HESS-II, the Cherenkov Telescope Array (CTA), and the GAMMA-400 satellite. As an example we use the narrow feature at 130 GeV seen in public data from the Fermi-LAT satellite. We found that all three experiments should be able to confidently confirm or rule out the presence of this 130 GeV feature. If it is real, it should be confirmed with a confidence level higher than 5 sigma. Assuming it to be a spectral signature of dark matter origin, GAMMA-400, thanks to a projected energy resolution of about 1.5% at 100 GeV, should also be able to resolve both the γγ-line and a corresponding Zγ- or Hγ-feature, if the corresponding branching ratio is comparable to that into two photons. It will also allow to distinguish between a gamma-ray line and the similar feature resulting from internal bremsstrahlung photons.

preprint2012arXiv

Searches for Particle Dark Matter with gamma-rays

In this contribution I review the present status and discuss some prospects for indirect detection of dark matter with gamma-rays. Thanks to the Fermi Large Area Telescope, searches in gamma-rays have reached sensitivities that allow to probe the most interesting parameter space of the weakly interacting massive particles (WIMP) paradigm. This gain in sensitivity is naturally accompanied by a number of detection claims or indications, the most recent being the claim of a line feature at a dark matter particle mass of $\sim$ 130 GeV at the Galactic Centre, a claim which requires confirmation from the Fermi-LAT collaboration and other experiments, for example HESS II or the planned Gamma-400 satellite. Predictions for the next generation air Cherenkov telescope, Cherenkov Telescope Array (CTA), together with forecasts on future Fermi-LAT constraints arrive at the exciting possibility that the cosmological benchmark cross-section could be probed from masses of a few GeV to a few TeV. Consequently, non-detection would pose a challenge to the WIMP paradigm, but the reached sensitivities also imply that --optimistically-- a detection is in the cards.

preprint2011arXiv

How well will ton-scale dark matter direct detection experiments constrain minimal supersymmetry?

Weakly interacting massive particles (WIMPs) are amongst the most interesting dark matter (DM) candidates. Many DM candidates naturally arise in theories beyond the standard model (SM) of particle physics, like weak-scale supersymmetry (SUSY). Experiments aim to detect WIMPs by scattering, annihilation or direct production, and thereby determine the underlying theory to which they belong, along with its parameters. Here we examine the prospects for further constraining the Constrained Minimal Supersymmetric Standard Model (CMSSM) with future ton-scale direct detection experiments. We consider ton-scale extrapolations of three current experiments: CDMS, XENON and COUPP, with 1000 kg-years of raw exposure each. We assume energy resolutions, energy ranges and efficiencies similar to the current versions of the experiments, and include backgrounds at target levels. Our analysis is based on full likelihood constructions for the experiments. We also take into account present uncertainties on hadronic matrix elements for neutralino-quark couplings, and on halo model parameters. We generate synthetic data based on four benchmark points and scan over the CMSSM parameter space using nested sampling. We construct both Bayesian posterior PDFs and frequentist profile likelihoods for the model parameters, as well as the mass and various cross-sections of the lightest neutralino. Future ton-scale experiments will help substantially in constraining supersymmetry, especially when results of experiments primarily targeting spin-dependent nuclear scattering are combined with those directed more toward spin-independent interactions.

preprint2011arXiv

Implications for constrained supersymmetry of combined H.E.S.S. observations of dwarf galaxies, the Galactic halo and the Galactic Centre

In order to place limits on dark matter (DM) properties using $γ$-ray observations, previous analyses have often assumed a very simple parametrisation of the $γ$-ray annihilation yield; typically, it has been assumed that annihilation proceeds through a single channel only. In realistic DM models, annihilation may occur into many different final states, making this quite a rough ansatz. With additional processes like virtual internal bremsstrahlung and final state radiation, this ansatz becomes even more incorrect, and the need for scans of explicit model parameter spaces becomes clear. Here we present scans of the parameter space of the Constrained Minimal Supersymmetric Standard Model (CMSSM), considering $γ$-ray spectra from three dwarf galaxies, the Galactic Centre region and the broader Galactic halo, as obtained with the High-Energy Stereoscopic System (H.E.S.S.). We present a series of likelihood scans combining the H.E.S.S. data with other experimental results. We show that observations of the Sagittarius, Carina and Sculptor dwarf galaxies disfavour the coannihilation region of the CMSSM and models with large annihilation cross-sections. This is true even under reasonable assumptions about the DM density profiles, and constitutes the strongest constraint to date on coannihilation models within the CMSSM. The Galactic halo has a similar, but weaker, effect. The Galactic Centre search is complicated by a strong (unknown) $γ$-ray source, and we see no effect on the CMSSM parameter space when assuming a realistic Galactic Centre DM density profile.

preprint2011arXiv

Statistical coverage for supersymmetric parameter estimation: a case study with direct detection of dark matter

Models of weak-scale supersymmetry offer viable dark matter (DM) candidates. Their parameter spaces are however rather large and complex, such that pinning down the actual parameter values from experimental data can depend strongly on the employed statistical framework and scanning algorithm. In frequentist parameter estimation, a central requirement for properly constructed confidence intervals is that they cover true parameter values, preferably at exactly the stated confidence level when experiments are repeated infinitely many times. Since most widely-used scanning techniques are optimised for Bayesian statistics, one needs to assess their abilities in providing correct confidence intervals in terms of the statistical coverage. Here we investigate this for the Constrained Minimal Supersymmetric Standard Model (CMSSM) when only constrained by data from direct searches for dark matter. We construct confidence intervals from one-dimensional profile likelihoods and study the coverage by generating several pseudo-experiments for a few benchmark sets of pseudo-true parameters. We use nested sampling to scan the parameter space and evaluate the coverage for the benchmarks when either flat or logarithmic priors are imposed on gaugino and scalar mass parameters. The sampling algorithm has been used in the configuration usually adopted for exploration of the Bayesian posterior. We observe both under- and over-coverage, which in some cases vary quite dramatically when benchmarks or priors are modified. We show how most of the variation can be explained as the impact of explicit priors as well as sampling effects, where the latter are indirectly imposed by physicality conditions. For comparison, we also evaluate the coverage for Bayesian credible intervals, and observe significant under-coverage in those cases.

preprint2010arXiv

A Profile Likelihood Analysis of the Constrained MSSM with Genetic Algorithms

The Constrained Minimal Supersymmetric Standard Model (CMSSM) is one of the simplest and most widely-studied supersymmetric extensions to the standard model of particle physics. Nevertheless, current data do not sufficiently constrain the model parameters in a way completely independent of priors, statistical measures and scanning techniques. We present a new technique for scanning supersymmetric parameter spaces, optimised for frequentist profile likelihood analyses and based on Genetic Algorithms. We apply this technique to the CMSSM, taking into account existing collider and cosmological data in our global fit. We compare our method to the MultiNest algorithm, an efficient Bayesian technique, paying particular attention to the best-fit points and implications for particle masses at the LHC and dark matter searches. Our global best-fit point lies in the focus point region. We find many high-likelihood points in both the stau co-annihilation and focus point regions, including a previously neglected section of the co-annihilation region at large m_0. We show that there are many high-likelihood points in the CMSSM parameter space commonly missed by existing scanning techniques, especially at high masses. This has a significant influence on the derived confidence regions for parameters and observables, and can dramatically change the entire statistical inference of such scans.

preprint2009arXiv

Direct Constraints on Minimal Supersymmetry from Fermi-LAT Observations of the Dwarf Galaxy Segue 1

The dwarf galaxy Segue 1 is one of the most promising targets for the indirect detection of dark matter. Here we examine what constraints 9 months of Fermi-LAT gamma-ray observations of Segue 1 place upon the Constrained Minimal Supersymmetric Standard Model (CMSSM), with the lightest neutralino as the dark matter particle. We use nested sampling to explore the CMSSM parameter space, simultaneously fitting other relevant constraints from accelerator bounds, the relic density, electroweak precision observables, the anomalous magnetic moment of the muon and B-physics. We include spectral and spatial fits to the Fermi observations, a full treatment of the instrumental response and its related uncertainty, and detailed background models. We also perform an extrapolation to 5 years of observations, assuming no signal is observed from Segue 1 in that time. Results marginally disfavour models with low neutralino masses and high annihilation cross-sections. Virtually all of these models are however already disfavoured by existing experimental or relic density constraints.

preprint2006arXiv

Fluctuation studies and energy reconstruction in a segmented calorimeter

In order to better understand energy estimation of electromagnetic showers in segmented calorimeters, detailed Geant4 simulation studies of electromagnetic showers in the energy range 1--100 GeV in CsI have been performed. When sampled in layers of 1.99 cm thickness, corresponding to 1.08 radiation lengths, the energy fluctuations in the samples show distributions that vary significantly with depth. The energy distributions are similar for incident electrons and photons and were found to change systematically along the shower, varying little with the initial energy. Three probability distributions have been fitted to the data: negative binomial, log-normal and Gaussian distributions, none of which gives a good fit over the full shower length. The obtained parameterizations might be useful in the use of the maximum likelihood method for estimating the energy. Two methods for estimating the energy have also been studied. One method is based on fitting individual longitudinal shower profiles with a $Γ$-function, the other one corrects the measured energy for leakage by utilizing the energy deposited in the last layer. Our simulations indicate that the last-layer correction method applied to photons and electrons of 1 and 10 GeV gives about a factor of 2 improvement in the energy resolution.

preprint2005arXiv

Beam Test Performance and Simulation of Prototypes for the ALICE Silicon Pixel Detector

The silicon pixel detector (SPD) of the ALICE experiment in preparation at the Large Hadron Collider (LHC) at CERN is designed to provide the precise vertex reconstruction needed for measuring heavy flavor production in heavy ion collisions at very high energies and high multiplicity. The SPD forms the innermost part of the Inner Tracking System (ITS) which also includes silicon drift and silicon strip detectors. Single assembly prototypes of the ALICE SPD have been tested at the CERN SPS using high energy proton/pion beams in 2002 and 2003. We report on the experimental determination of the spatial precision. We also report on the first combined beam test with prototypes of the other ITS silicon detector technologies at the CERN SPS in November 2004. The issue of SPD simulation is briefly discussed.