Source author record

Andrew Fowlie

Andrew Fowlie 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

27works
10topics
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

27 published item(s)

preprint2022arXiv

Comment on "Accumulating Evidence for the Associate Production of a Neutral Scalar with Mass around 151 GeV"

A recent paper [2109.02650] accumulates evidence for a new fundamental particle by combining several CMS and ATLAS searches for the Standard Model Higgs boson. The putative particle is a neutral scalar, $S$, with a mass of about 151 GeV. The reported significances are $5.1σ$ local and $4.8σ$ global. This nearly reaches the $5σ$ threshold for a discovery in high-energy physics. In this brief note we cast doubt on the strength of the evidence for a new particle. After taking into account the fact that signals were fitted to six different channels, we find that the significances are only $4.1σ$ local and $3.5σ$ global. The code and instructions for reproducing our calculations are available at https://github.com/andrewfowlie/accumulating_evidence.

preprint2022arXiv

Global fits of SUSY at future Higgs factories

In this work, we study the impact of electroweak and Higgs precision measurements at future electron-positron colliders on several typical supersymmetric models, including the Constrained Minimal Supersymmetric Standard Model (CMSSM), Non-Universal Higgs Mass generalisations (NUHM1, NUHM2), and the 7-dimensional Minimal Supersymmetric Standard Model (MSSM7). Using publicly-available data from the \textsf{GAMBIT} community, we post-process previous SUSY global fits with additional likelihoods to explore the discovery potential of Higgs factories, such as the Circular Electron Positron Collider (CEPC), the Future Circular Collider (FCC) and the International Linear Collider (ILC). We show that the currently allowed parameter space of these models will be further tested by future precision measurements. In particular, dark matter annihilation mechanisms may be distinguished by precise measurements of Higgs observables.

preprint2022arXiv

Nested sampling for frequentist computation: fast estimation of small $p$-values

We propose a novel method for computing $p$-values based on nested sampling (NS) applied to the sampling space rather than the parameter space of the problem, in contrast to its usage in Bayesian computation. The computational cost of NS scales as $\log^2{1/p}$, which compares favorably to the $1/p$ scaling for Monte Carlo (MC) simulations. For significances greater than about $4σ$ in both a toy problem and a simplified resonance search, we show that NS requires orders of magnitude fewer simulations than ordinary MC estimates. This is particularly relevant for high-energy physics, which adopts a $5σ$ gold standard for discovery. We conclude with remarks on new connections between Bayesian and frequentist computation and possibilities for tuning NS implementations for still better performance in this setting.

preprint2022arXiv

Nested sampling for physical scientists

We review Skilling's nested sampling (NS) algorithm for Bayesian inference and more broadly multi-dimensional integration. After recapitulating the principles of NS, we survey developments in implementing efficient NS algorithms in practice in high-dimensions, including methods for sampling from the so-called constrained prior. We outline the ways in which NS may be applied and describe the application of NS in three scientific fields in which the algorithm has proved to be useful: cosmology, gravitational-wave astronomy, and materials science. We close by making recommendations for best practice when using NS and by summarizing potential limitations and optimizations of NS.

preprint2022arXiv

Objective Bayesian approach to the Jeffreys-Lindley paradox

We consider the Jeffreys-Lindley paradox from an objective Bayesian perspective by attempting to find priors representing complete indifference to sample size in the problem. This means that we ensure that the prior for the unknown mean and the prior predictive for the $t$-statistic are independent of the sample size. If successful, this would lead to Bayesian model comparison that was independent of sample size and ameliorate the paradox. Unfortunately, it leads to an improper scale-invariant prior for the unknown mean. We show, however, that a truncated scale-invariant prior delays the dependence on sample size, which could be practically significant. Lastly, we shed light on the paradox by relating it to the fact that the scale-invariant prior is improper.

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

Nested sampling with plateaus

It was recently emphasised by Riley (2019); Schittenhelm & Wacker (2020) that that in the presence of plateaus in the likelihood function nested sampling (NS) produces faulty estimates of the evidence and posterior densities. After informally explaining the cause of the problem, we present a modified version of NS that handles plateaus and can be applied retrospectively to NS runs from popular NS software using anesthetic. In the modified NS, live points in a plateau are evicted one by one without replacement, with ordinary NS compression of the prior volume after each eviction but taking into account the dynamic number of live points. The live points are replenished once all points in the plateau are removed. We demonstrate it on a number of examples. Since the modification is simple, we propose that it becomes the canonical version of Skilling's NS algorithm.

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.

preprint2021arXiv

Thermal WIMPs and the Scale of New Physics: Global Fits of Dirac Dark Matter Effective Field Theories

We assess the status of a wide class of WIMP dark matter (DM) models in light of the latest experimental results using the global fitting framework $\textsf{GAMBIT}$. We perform a global analysis of effective field theory (EFT) operators describing the interactions between a gauge-singlet Dirac fermion and the Standard Model quarks, the gluons and the photon. In this bottom-up approach, we simultaneously vary the coefficients of 14 such operators up to dimension 7, along with the DM mass, the scale of new physics and several nuisance parameters. Our likelihood functions include the latest data from $\mathit{Planck}$, direct and indirect detection experiments, and the LHC. For DM masses below 100 GeV, we find that it is impossible to satisfy all constraints simultaneously while maintaining EFT validity at LHC energies. For new physics scales around 1 TeV, our results are influenced by several small excesses in the LHC data and depend on the prescription that we adopt to ensure EFT validity. Furthermore, we find large regions of viable parameter space where the EFT is valid and the relic density can be reproduced, implying that WIMPs can still account for the DM of the universe while being consistent with the latest data.

preprint2020arXiv

Global analyses of Higgs portal singlet dark matter models using GAMBIT

We present global analyses of effective Higgs portal dark matter models in the frequentist and Bayesian statistical frameworks. Complementing earlier studies of the scalar Higgs portal, we use GAMBIT to determine the preferred mass and coupling ranges for models with vector, Majorana and Dirac fermion dark matter. We also assess the relative plausibility of all four models using Bayesian model comparison. Our analysis includes up-to-date likelihood functions for the dark matter relic density, invisible Higgs decays, and direct and indirect searches for weakly-interacting dark matter including the latest XENON1T data. We also account for important uncertainties arising from the local density and velocity distribution of dark matter, nuclear matrix elements relevant to direct detection, and Standard Model masses and couplings. In all Higgs portal models, we find parameter regions that can explain all of dark matter and give a good fit to all data. The case of vector dark matter requires the most tuning and is therefore slightly disfavoured from a Bayesian point of view. In the case of fermionic dark matter, we find a strong preference for including a CP-violating phase that allows suppression of constraints from direct detection experiments, with odds in favour of CP violation of the order of 100:1. Finally, we present DDCalc 2.0.0, a tool for calculating direct detection observables and likelihoods for arbitrary non-relativistic effective operators.

preprint2020arXiv

Nested sampling cross-checks using order statistics

Nested sampling (NS) is an invaluable tool in data analysis in modern astrophysics, cosmology, gravitational wave astronomy and particle physics. We identify a previously unused property of NS related to order statistics: the insertion indexes of new live points into the existing live points should be uniformly distributed. This observation enabled us to create a novel cross-check of single NS runs. The tests can detect when an NS run failed to sample new live points from the constrained prior and plateaus in the likelihood function, which break an assumption of NS and thus leads to unreliable results. We applied our cross-check to NS runs on toy functions with known analytic results in 2 - 50 dimensions, showing that our approach can detect problematic runs on a variety of likelihoods, settings and dimensions. As an example of a realistic application, we cross-checked NS runs performed in the context of cosmological model selection. Since the cross-check is simple, we recommend that it become a mandatory test for every applicable NS run.

preprint2020arXiv

PhaseTracer: tracing cosmological phases and calculating transition properties

We present a C++ software package called PhaseTracer for mapping out cosmological phases, and potential transitions between them, for Standard Model extensions with any number of scalar fields. PhaseTracer traces the minima of effective potential as the temperature changes, and then calculates the critical temperatures, at which the minima are degenerate. PhaseTracer is constructed with modularity, flexibility and practicality in mind. It is fast and stable, and can receive potentials provided by other packages such as FlexibleSUSY. PhaseTracer can be useful analysing cosmological phase transitions which played an important role in the very early evolution of the Universe. If they were first order they could generate detectable gravitational waves and/or trigger electroweak baryogenesis to generate the observed matter anti-matter asymmetry of the Universe. The code can be obtained from https://github.com/PhaseTracer/PhaseTracer.

preprint2019arXiv

Combined collider constraints on neutralinos and charginos

Searches for supersymmetric electroweakinos have entered a crucial phase, as the integrated luminosity of the Large Hadron Collider is now high enough to compensate for their weak production cross-sections. Working in a framework where the neutralinos and charginos are the only light sparticles in the Minimal Supersymmetric Standard Model, we use gambit to perform a detailed likelihood analysis of the electroweakino sector. We focus on the impacts of recent ATLAS and CMS searches with 36 fb$^{-1}$ of 13 TeV proton-proton collision data. We also include constraints from LEP and invisible decays of the $Z$ and Higgs bosons. Under the background-only hypothesis, we show that current LHC searches do not robustly exclude any range of neutralino or chargino masses. However, a pattern of excesses in several LHC analyses points towards a possible signal, with neutralino masses of $(m_{\tildeχ_1^0}, m_{\tildeχ_2^0}, m_{\tildeχ_3^0}, m_{\tildeχ_4^0})$ = (8-155, 103-260, 130-473, 219-502) GeV and chargino masses of $(m_{\tildeχ_1^{\pm}}, m_{\tildeχ_2^{\pm}})$ = (104-259, 224-507) GeV at the 95% confidence level. The lightest neutralino is mostly bino, with a possible modest Higgsino or wino component. We find that this excess has a combined local significance of $3.3σ$, subject to a number of cautions. If one includes LHC searches for charginos and neutralinos conducted with 8 TeV proton-proton collision data, the local significance is lowered to 2.9$σ$. We briefly consider the implications for dark matter, finding that the correct relic density can be obtained through the Higgs-funnel and $Z$-funnel mechanisms, even assuming that all other sparticles are decoupled. All samples, gambit input files and best-fit models from this study are available on Zenodo.

preprint2019arXiv

Strong first-order phase transitions in the NMSSM --- a comprehensive survey

Motivated by the fact that the Next-to-Minimal Supersymmetric Standard Model is one of the most plausible models that can accommodate electroweak baryogenesis, we analyze its phase structure by tracing the temperature dependence of the minima of the effective potential. Our results reveal rich patterns of phase structure that end in the observed electroweak symmetry breaking vacuum. We classify these patterns according to the first transition in their history and show the strong first-order phase transitions that may be possible in each type of pattern. These could allow for the generation of the matter-antimatter asymmetry or potentially observable gravitational waves. For a selection of benchmark points, we checked that the phase transitions completed and calculated the nucleation temperatures. We furthermore present samples that feature strong first-order phase transitions from an extensive scan of the whole parameter space. We highlight common features of our samples, including the fact that the Standard Model like Higgs is often not the lightest Higgs in the model.

preprint2016arXiv

LHCO_reader: A new code for reading and analyzing detector-level events stored in LHCO format

We present a new Python module - LHCO_reader - for reading detector-level events in LHCO format generated from detector-simulators such as PGS and Delphes. Emphasis is placed upon ease of use. The module is installed via pip. Once installed, one can read and scrutinize events, with helpful functions and classes for cutting, plotting and inspecting events, and manipulating four-momenta by e.g. boosting. Implementing complicated analyses should be straight-forward by coding any cuts into Python. Furthermore, common kinematic variables, including $α_T$ and razor are included out of the box, and $m_{T2}$ and $m_{T2}^W$ are provided by interfaces with external libraries.

preprint2016arXiv

Naturalness of the relaxion mechanism

The relaxion mechanism is a novel solution to the hierarchy problem. In this first statistical analysis of the relaxion mechanism, we quantify the relative plausibility of a QCD and a non-QCD relaxion model versus the Standard Model with Bayesian statistics, which includes an automatic penalty for fine-tuning. We find that in light of the hierarchy between the weak and Planck scales, relaxion models are favoured by colossal Bayes-factors. Constraints upon \eg the vacuum energy during relaxation, however, shrink the Bayes-factors such that relaxion models are only slightly favoured. Including the bounds on $\left|θ_\text{QCD}\right|$ shatters the plausibility of the QCD relaxion model as it typically yields $\left|θ_\text{QCD}\right| \gg 0$. Finally, we augment our models with scalar-field inflation and consider measurements of inflationary observables from BICEP/Planck. We find that, all told, the Standard Model is favoured by huge Bayes-factors as the relaxion models require fine-tuning such that the Hubble parameter is less than the height of the periodic barriers. Thus, whilst we confirm that relaxion models could solve the hierarchy problem, we find that their unconventional cosmology demolishes their plausibility.

preprint2016arXiv

Superplot: a graphical interface for plotting and analysing MultiNest output

We present an application, Superplot, for calculating and plotting statistical quantities relevant to parameter inference from a "chain" of samples drawn from a parameter space, produced by e.g. MultiNest. A simple graphical interface allows one to browse a chain of many variables quickly, and make publication quality plots of, inter alia, one- and two-dimensional profile likelihood, posterior pdf (with kernel density estimation), confidence intervals and credible regions. In this short manual, we document installation and basic usage, and define all statistical quantities and conventions. The code is fully compatible with Linux, Windows and Mac OSX. Furthermore, if preferred, all functionality is available through the command line rather than a graphical interface.

preprint2015arXiv

Examining a right-handed quark mixing matrix with $b$-tags at the LHC

Encouraged by a hint in a search for right-handed $W$ bosons at the LHC, we investigate whether the unitarity of a right-handed quark mixing matrix and the equality of the left- and right-handed quark mixing matrices could be tested at the LHC. We propose a particular test, involving counting the numbers of $b$-tags in the final state, and simulate the test at the event level with Monte-Carlo tools for the forthcoming $\sqrt{s}=13$ TeV LHC run. We find that testing unitarity with 20/fb will be challenging; our test successfully rejects unitarity if the right-handed quark mixing matrix is non-unitary, but only in particular cases. On the other hand, our test may provide the first opportunity to test the unitarity of a right-handed quark mixing matrix and with 3000/fb severely constrains possible departures from unitarity in the latter. We refine our previous work, testing the equality of quark mixing matrices, with full collider simulation. With 20/fb, we are sensitive to mixing angles as small as $30^\circ$, and with 3000/fb, angles as small as $7.5^\circ$, confirming our preliminary analysis. We briefly investigate testing the unitarity of the SM CKM matrix with a similar method by studying semileptonic $t\bar t$ production, concluding that systematics make it particularly difficult.

preprint2015arXiv

The little-hierarchy problem is a little problem: understanding the difference between the big- and little-hierarchy problems with Bayesian probability

Experiments are once again under way at the LHC. This time around, however, the mood in the high-energy physics community is pessimistic. There is a growing suspicion that naturalness arguments that predict new physics near the weak scale are faulty and that prospects for a new discovery are limited. We argue that such doubts originate from a misunderstanding of the foundations of naturalness arguments. In spite of the first run at the LHC, which aggravated the little-hierarchy problem, there is no cause for doubting naturalness or natural theories. Naturalness is grounded in Bayesian probability logic - it is not a scientific theory and it makes no sense to claim that it could be falsified or that it is under pressure from experimental data. We should remain optimistic about discovery prospects; natural theories, such as supersymmetry, generally predict new physics close to the weak scale. Furthermore, from a Bayesian perspective, we briefly discuss 't Hooft's technical naturalness and a contentious claim that the little-hierarchy problem hints that the Standard Model is a fundamental theory.

preprint2014arXiv

CMSSM, naturalness and the "fine-tuning price" of the Very Large Hadron Collider

The absence of supersymmetry or other new physics at the Large Hadron Collider (LHC) has lead many to question naturalness arguments. With Bayesian statistics, we argue that natural models are most probable and that naturalness is not merely an aesthetic principle. We calculate a probabilistic measure of naturalness, the Bayesian evidence, for the Standard Model (SM) with and without quadratic divergences, confirming that the SM with quadratic divergences is improbable. We calculate the Bayesian evidence for the Constrained Minimal Supersymmetric Standard Model (CMSSM) with naturalness priors in three cases: with only the $M_Z$ measurement; with the $M_Z$ measurement and LHC measurements; and with the $M_Z$ measurement, $m_h$ measurement and a hypothetical null result from a $\sqrt{s}=100\,\text{TeV}$ Very Large Hadron Collider (VLHC) with $3000/\text{fb}$. The "fine-tuning price" of the VLHC given LHC results would be $\sim400$, which is slightly less than that of the LHC results given the electroweak scale ($\sim500$).

preprint2014arXiv

Is the CNMSSM more credible than the CMSSM?

With Bayesian statistics, we investigate the full parameter space of the constrained "next-to-minimal" supersymmetric Standard Model (CNMSSM) with naturalness priors, which were derived in a previous work. In the past, most Bayesian analyses of the CNMSSM ignored naturalness of the electroweak (EW) scale by making prejudicial assumptions for parameters defined at the EW scale. We test the CNMSSM against the CMSSM with Bayesian evidence, which, with naturalness priors, incorporates a penalty for fine-tuning of the EW scale. With the evidence, we measure credibility with respect to all measurements, including the EW scale and LHC direct searches. We find that the evidence in favor of the CNMSSM versus the CMSSM is "positive" to "strong" but that if one ignores the $μ$-problem, the evidence is "barely worth mentioning" to "positive." The $μ$-problem significantly influences our findings. Unless one considers the $μ$-problem, the evidence in favor of the CNMSSM versus the CMSSM is at best "positive," which is two grades below "very strong." We, furthermore, identify the most probable regions of the CMSSM and CNMSSM parameter spaces and examine prospects for future discovery at hadron colliders.

preprint2014arXiv

Prospects for constrained supersymmetry at $\sqrt{s}=33$ TeV and $\sqrt{s}=100$ TeV proton-proton super-colliders

Discussions are underway for a high-energy proton-proton collider. Two preliminary ideas are the $\sqrt{s}=33$ TeV HE-LHC and the $\sqrt{s}=100$ TeV VLHC. With Bayesian statistics, we calculate the probabilities that the LHC, HE-LHC and VLHC discover SUSY in the future, assuming that nature is described by the CMSSM and given the experimental data from the LHC, LUX and Planck. We find that the LHC with $300$/fb at $\sqrt{s}=14$ TeV has a $15$-$75%$ probability of discovering SUSY. Should that run fail to discover SUSY, the probability of discovering SUSY with $3000$/fb is merely $1$-$10%$. Were SUSY to remain undetected at the LHC, the HE-LHC would have a $35$-$85%$ probability of discovering SUSY with $3000$/fb. The VLHC, on the other hand, ought to be definitive; the probability of it discovering SUSY, assuming that the CMSSM is the correct model, is $100%$.

preprint2014arXiv

Reconstructing CMSSM parameters at the LHC with $\sqrt{s}=14$ TeV via the golden decay channel

We identify a benchmark point in the CMSSM's heavy stau-coannihilation region, which is favored by experiments, and demonstrate that it could be accessible to the LHC at $\sqrt{s}=14$ TeV with 300/fb of integrated luminosity via a golden decay measurement. With Monte-Carlo, we simulate sparticle production and subsequent golden decay at the event level and perform pseudo-measurements of sparticle masses from kinematic endpoints in invariant mass distributions. We find that two lightest neutralino masses and the first and second generation left-handed slepton and squark masses could be rather precisely measured with correlated uncertainties. We investigate whether from such measurements one could determine the CMSSM's Lagrangian parameters by including a likelihood from our pseudo-measurements of sparticle masses in a Bayesian analysis of the CMSSM's parameter space. We find that the CMSSM's parameters can be accurately determined, with the exception of the common trilinear parameter. Experimental measurements of the relic density by Planck and the Higgs boson's mass slightly improve this determination, especially for the common trilinear parameter. Finally, within our benchmark scenario, we show that the neutralino dark matter will be accessible to direct searches in future one tonne detectors.

preprint2014arXiv

Testing quark mixing in minimal left-right symmetric models with $b$-tags at the LHC

Motivated by a hint in a CMS search for right-handed $W$-bosons in $eejj$ final states, we propose an experimental test of quark-mixing matrices in a general left-right symmetric model, based on counting the numbers of $b$-tags from right-handed $W$-boson hadronic decays. We find that, with our test, differences between left- and right-handed quark-mixing matrices could be detected at the LHC with $\sqrt{s}=14$ TeV. With an integrated luminosity of about $20/\textrm{fb}$, our test is sensitive to right-handed quark-mixing angles as small as about $30^\circ$ and with $3000/\textrm{fb}$, our test's sensitivity improves to right-handed mixing angles as small as about $7.5^\circ$. Our test's sensitivity might be further enhanced by tuning $b$-tagging efficiency against purity.

preprint2013arXiv

Dark matter and collider signatures of the MSSM

We explore the MSSM with 9 free parameters (p9MSSM) that have been selected as a minimum set that allows an investigation of neutralino dark matter and collider signatures while maintaining consistency with several constraints. These include measurement of the dark matter relic density from PLANCK, main properties of the discovered Higgs boson, LHC direct SUSY searches, recent evidence for a Standard Model-like BR(Bs -> μ+ μ-), and the measurement of δ(g-2), plus a number of other electroweak and flavor physics constraints. We perform a simulation of two LHC direct SUSY searches at sqrt(s)=8 TeV: the CMS inclusive α_T search for squarks and gluinos and the CMS electroweak production search with 3l+E_T^miss in the final state. We use the latter to identify the regions of the parameter space, consistent at 2σ with δ(g-2), that are not excluded by the direct limits from the electroweak production. We find that they correspond to a neutralino mass in the window 200 GeV<m_χ<500 GeV. We also implement the likelihood for the XENON100 exclusion bound, in which we consider for the first time the impact of a recent determination of the Σ_{πN} term from CHAOS data, Σ_{πN}=43\pm12 MeV. We show that in light of this measurement, the present statistical impact of the XENON100 bound is greatly reduced, although future sensitivities of the LUX and XENON1T experiments will have decisive impact on the mixed bino/higgsino composition of the neutralino. We point out some tension between the constraints from δ(g-2) and XENON100. Finally, we present prospects for various indirect searches of dark matter, namely γ-ray fluxes from dSphs and the Galactic Center at Fermi-LAT, and the positron flux at AMS02. We also show the 5-year sensitivity on the spin-dependent neutralino-proton cross section due to neutrino fluxes from the Sun at IceCube.

preprint2012arXiv

Bayesian Implications of Current LHC and XENON100 Search Limits for the Constrained MSSM

The CMS Collaboration has released the results of its search for supersymmetry, by applying an alphaT method to 1.1/fb of data at 7 TeV. The null result excludes (at 95% CL) a low-mass region of the Constrained MSSM's parameter space that was previously favored by other experiments. Additionally, the negative result of the XENON100 dark matter search has excluded (at 90% CL) values of the spin-independent scattering cross sections sigma^SI_p as low as 10^-8 pb. We incorporate these improved experimental constraints into a global Bayesian fit of the Constrained MSSM by constructing approximate likelihood functions. In the case of the alphaT limit, we simulate detector efficiency for the CMS alphaT 1.1/fb and validate our method against the official 95% CL contour. We identify the 68% and 95% credible posterior regions of the CMSSM parameters, and also find the best-fit point. We find that the credible regions change considerably once a likelihood from alphaT is included, in particular the narrow light Higgs resonance region becomes excluded, but the focus point/horizontal branch region remains allowed at the 1sigma level. Adding the limit from XENON100 has a weaker additional effect, in part due to large uncertainties in evaluating sigma^SI_p, which we include in a conservative way, although we find that it reduces the posterior probability of the focus point region to the 2sigma level. The new regions of high posterior favor squarks lighter than the gluino and all but one Higgs bosons heavy. The dark matter neutralino mass is found in the range 250 GeV <~ m_Chi1 <~ 343 GeV (at 1sigma) while, as the result of improved limits from the LHC, the favored range of sigma^SI_p is pushed down to values below 10^{-9} pb. We highlight tension between (g-2)_mu and BR(b->sg), which is exacerbated by including the alphaT limit; each constraint favors a different region of the CMSSM's mass parameters.

preprint2012arXiv

Constrained MSSM favoring new territories: The impact of new LHC limits and a 125 GeV Higgs boson

We present an updated and extended global analysis of the Constrained MSSM (CMSSM) taking into account new limits on supersymmetry from ~5/fb data sets at the LHC. In particular, in the case of the razor limit obtained by the CMS Collaboration we simulate detector efficiency for the experimental analysis and derive an approximate but accurate likelihood function. We discuss the impact on the global fit of a possible Higgs boson with mass near 125 GeV, as implied by recent data, and of a new improved limit on BR(B_s->μμ). We identify high posterior probability regions of the CMSSM parameters as the stau-coannihilation and the A-funnel region, with the importance of the latter now being much larger due to the combined effect of the above three LHC results and of dark matter relic density. We also find that the focus point region is now disfavored. Ensuing implications for superpartner masses favor even larger values than before, and even lower ranges for dark matter spin-independent cross section, σ^{SI}_p<10^{-9} pb. We also find that relatively minor variations in applying experimental constraints can induce a large shift in the location of the best-fit point. This puts into question the robustness of applying the usual chisquare approach to the CMSSM. We discuss the goodness-of-fit and find that, while it is difficult to calculate a p-value, the g-2 constraint makes, nevertheless, the overall fit of the CMSSM poor. We consider a scan without this constraint, and we allow μ to be either positive or negative. We find that the global fit improves enormously for both signs of μ, with a slight preference for μ<0 caused by a better fit to BR(b->sγ) and BR(B_s->μμ).