Search

High-Energy Physics experiments are facing a multi-fold data increase with every new iteration. This is certainly the case for the upcoming High-Luminosity LHC upgrade. Such increased data processing requirements forces revisions to almost every step of the data processing pipeline. One such step in need of an overhaul is the task of particle track reconstruction, a.k.a., tracking. A Machine Learning-assisted solution is expected to provide significant improvements, since the most time-consuming step in tracking is the assignment of hits to particles or track candidates. This is the topic of this paper. We take inspiration from large language models. As such, we consider two approaches: the prediction of the next word in a sentence (next hit point in a track), as well as the one-shot prediction of all hits within an event. In an extensive design effort, we have experimented with three models based on the Transformer architecture and one model based on the U-Net architecture, performing track association predictions for collision event hit points. In our evaluation, we consider a spectrum of simple to complex representations of the problem, eliminating designs with lower metrics early on. We report extensive results, covering both prediction accuracy (score) and computational performance. We have made use of the REDVID simulation framework, as well as reductions applied to the TrackML data set, to compose five data sets from simple to complex, for our experiments. The results highlight distinct advantages among different designs in terms of prediction accuracy and computational performance, demonstrating the efficiency of our methodology. Most importantly, the results show the viability of a one-shot encoder-classifier based Transformer solution as a practical approach for the task of tracking.

In this work, we present a novel approach to transform supervised classifiers into effective unsupervised anomaly detectors. The method we have developed, termed Discriminatory Detection of Distortions (DDD), enhances anomaly detection by training a discriminator model on both original and artificially modified datasets. We conducted a comprehensive evaluation of our models on the Dark Machines Anomaly Score Challenge channels and a search for 4-top quark events, demonstrating the effectiveness of our approach across various final states and beyond the Standard Model scenarios. We compare the performance of the DDD method with the Deep Robust One-Class Classification method (DROCC), which incorporates signals in the training process, and the Deep Support Vector Data Description (DeepSVDD) method, a well established and well performing method for anomaly detection. Results show that the effectiveness of each model varies by signal and channel, with DDD proving to be a very effective anomaly detector. We recommend the combined use of DeepSVDD and DDD for purely unsupervised applications, with the addition of flow models for improved performance when resources allow. Findings suggest that network architectures that excel in supervised contexts, such as the particle transformer with standard model interactions, also perform well as unsupervised anomaly detectors. We also show that with these methods, it is likely possible to recognize 4-top quark production as an anomaly without prior knowledge of the process. We argue that the Large Hadron Collider community can transform supervised classifiers into anomaly detectors to uncover potential new physical phenomena in each search.

Track reconstruction is a vital aspect of High-Energy Physics (HEP) and plays a critical role in major experiments. In this study, we delve into unexplored avenues for particle track reconstruction and hit clustering. Firstly, we enhance the algorithmic design effort by utilising a simplified simulator (REDVID) to generate training data that is specifically composed for simplicity. We demonstrate the effectiveness of this data in guiding the development of optimal network architectures. Additionally, we investigate the application of image segmentation networks for this task, exploring their potential for accurate track reconstruction. Moreover, we approach the task from a different perspective by treating it as a hit sequence to track sequence translation problem. Specifically, we explore the utilisation of Transformer architectures for tracking purposes. Our preliminary findings are covered in detail. By considering this novel approach, we aim to uncover new insights and potential advancements in track reconstruction. This research sheds light on previously unexplored methods and provides valuable insights for the field of particle track reconstruction and hit clustering in HEP.

Subatomic particle track reconstruction (tracking) is a vital task in High-Energy Physics experiments. Tracking is exceptionally computationally challenging and fielded solutions, relying on traditional algorithms, do not scale linearly. Machine Learning (ML) assisted solutions are a promising answer. We argue that a complexity-reduced problem description and the data representing it, will facilitate the solution exploration workflow. We provide the REDuced VIrtual Detector (REDVID) as a complexity-reduced detector model and particle collision event simulator combo. REDVID is intended as a simulation-in-the-loop, to both generate synthetic data efficiently and to simplify the challenge of ML model design. The fully parametric nature of our tool, with regards to system-level configuration, while in contrast to physics-accurate simulations, allows for the generation of simplified data for research and education, at different levels. Resulting from the reduced complexity, we showcase the computational efficiency of REDVID by providing the computational cost figures for a multitude of simulation benchmarks. As a simulation and a generative tool for ML-assisted solution design, REDVID is highly flexible, reusable and open-source. Reference data sets generated with REDVID are publicly available. Data generated using REDVID has enabled rapid development of multiple novel ML model designs, which is currently ongoing.

Neutral particles capable of travelling cosmic distances from a source to detectors on Earth are limited to photons and neutrinos. Examination of the Dark Matter annihilation/decay spectra for these particles reveals the presence of continuum spectra (e.g. due to fragmentation and W or Z decay) and peaks (due to direct annihilations/decays). However, when one explores extensions of the Standard Model (BSM), unexplored spectra emerge that differ significantly from those of the Standard Model (SM) for both neutrinos and photons. In this paper, we argue for the inclusion of important spectra that include peaks as well as previously largely unexplored entities such as boxes and combinations of box, peak and continuum decay spectra., Comment: 7 pages, 2 figures

Scientific understanding is a fundamental goal of science, allowing us to explain the world. There is currently no good way to measure the scientific understanding of agents, whether these be humans or Artificial Intelligence systems. Without a clear benchmark, it is challenging to evaluate and compare different levels of and approaches to scientific understanding. In this Roadmap, we propose a framework to create a benchmark for scientific understanding, utilizing tools from philosophy of science. We adopt a behavioral notion according to which genuine understanding should be recognized as an ability to perform certain tasks. We extend this notion by considering a set of questions that can gauge different levels of scientific understanding, covering information retrieval, the capability to arrange information to produce an explanation, and the ability to infer how things would be different under different circumstances. The Scientific Understanding Benchmark (SUB), which is formed by a set of these tests, allows for the evaluation and comparison of different approaches. Benchmarking plays a crucial role in establishing trust, ensuring quality control, and providing a basis for performance evaluation. By aligning machine and human scientific understanding we can improve their utility, ultimately advancing scientific understanding and helping to discover new insights within machines.

The Galactic Center Excess (GCE) in GeV gamma rays has been debated for over a decade, with the possibility that it might be due to dark matter annihilation or undetected point sources such as millisecond pulsars (MSPs). This study investigates how the gamma-ray emission model ($\gamma$EM) used in Galactic center analyses affects the interpretation of the GCE's nature. To address this issue, we construct an ultra-fast and powerful inference pipeline based on convolutional Deep Ensemble Networks. We explore the two main competing hypotheses for the GCE using a set of $\gamma$EMs with increasing parametric freedom. We calculate the fractional contribution ($f_{\mathrm{src}}$) of a dim population of MSPs to the total luminosity of the GCE and analyze its dependence on the complexity of the $\gamma$EM. For the simplest $\gamma$EM, we obtain $f_{\mathrm{src}} = 0.10 \pm 0.07$, while the most complex model yields $f_{\mathrm{src}} = 0.79 \pm 0.24.$ In conclusion, we find that the statement about the nature of the GCE (dark matter or not) strongly depends on the assumed $\gamma$EM. The quoted results for $f_{\mathrm{src}}$ do not account for the additional uncertainty arising from the fact that the observed gamma-ray sky is out-of-distribution concerning the investigated $\gamma$EM iterations. We quantify the reality gap between our $\gamma$EMs using deep-learning-based One-Class Deep Support Vector Data Description networks, revealing that all employed $\gamma$EMs have gaps to reality. Our study casts doubt on the validity of previous conclusions regarding the GCE and dark matter, and underscores the urgent need to account for the reality gap and consider previously overlooked ''out of domain'' uncertainties in future interpretations., Comment: 56 pages, 25 figures; comments welcome! Accepted for submission to JCAP; text coincides with the published version

A major task in particle physics is the measurement of rare signal processes. These measurements are highly dependent on the classification accuracy of these events in relation to the huge background of other Standard Model processes. Reducing the background by a few tens of percent with the same signal efficiency can already increase the sensitivity considerably. This work demonstrates the importance of incorporating physical information into deep learning-based event selection. The paper includes this information into different methods for classifying events, in particular Boosted Decision Trees, Transformer Architectures (Particle Transformer) and Graph Neural Networks (Particle Net). In addition to the physical information previously proposed for jet tagging, we add particle measures for energy-dependent particle-particle interaction strengths as predicted by the leading order interactions of the Standard Model (SM). We find that the integration of physical information into the attention matrix (transformers) or edges (graphs) notably improves background rejection by $10\%$ to $40\%$ over baseline models (a graph network), with about $10\%$ of this improvement directly attributable to what we call the SM interaction matrix. In a simplified statistical analysis, we find that such architectures can improve the significance of signals by a significant factor compared to a graph network (our base model)., Comment: Version 2 of this paper supersedes arXiv:2211.05143v1. 26 pages, 5 figures

The massive use of artificial neural networks (ANNs), increasingly popular in many areas of scientific computing, rapidly increases the energy consumption of modern high-performance computing systems. An appealing and possibly more sustainable alternative is provided by novel neuromorphic paradigms, which directly implement ANNs in hardware. However, little is known about the actual benefits of running ANNs on neuromorphic hardware for use cases in scientific computing. Here we present a methodology for measuring the energy cost and compute time for inference tasks with ANNs on conventional hardware. In addition, we have designed an architecture for these tasks and estimate the same metrics based on a state-of-the-art analog in-memory computing (AIMC) platform, one of the key paradigms in neuromorphic computing. Both methodologies are compared for a use case in quantum many-body physics in two dimensional condensed matter systems and for anomaly detection at 40 MHz rates at the Large Hadron Collider in particle physics. We find that AIMC can achieve up to one order of magnitude shorter computation times than conventional hardware, at an energy cost that is up to three orders of magnitude smaller. This suggests great potential for faster and more sustainable scientific computing with neuromorphic hardware., Comment: 7 pages, 4 figures, submitted to and accepted by APL Machine Learning

Model-independent search strategies have been increasingly proposed in recent years because on the one hand there has been no clear signal for new physics and on the other hand there is a lack of a highly probable and parameter-free extension of the standard model. For these reasons, there is no simple search target so far. In this work, we try to take a new direction and ask the question: bearing in mind that we have a large number of new physics theories that go beyond the Standard Model and may contain a grain of truth, can we improve our search strategy for unknown signals by using them "in combination"? In particular, we show that a signal hypothesis based on a large, intermingled set of many different theoretical signal models can be a superior approach to find an unknown BSM signal. Applied to a recent data challenge, we show that "mixture-of-theories training" outperforms strategies that optimize signal regions with a single BSM model as well as most unsupervised strategies. Applications of this work include anomaly detection and the definition of signal regions in the search for signals of new physics., Comment: 40 pages, 10 figures

Neutrino telescope experiments are rapidly becoming more competitive in indirect detection searches for dark matter. Neutrino signals arising from dark matter annihilations are typically assumed to originate from the hadronisation and decay of Standard Model particles. Here we showcase a supersymmetric model, the BLSSMIS, that can simultaneously obey current experimental limits while still providing a potentially observable non-standard neutrino spectrum from dark matter annihilation., Comment: 10 pages, 6 figures

First-principle simulations are at the heart of the high-energy physics research program. They link the vast data output of multi-purpose detectors with fundamental theory predictions and interpretation. This review illustrates a wide range of applications of modern machine learning to event generation and simulation-based inference, including conceptional developments driven by the specific requirements of particle physics. New ideas and tools developed at the interface of particle physics and machine learning will improve the speed and precision of forward simulations, handle the complexity of collision data, and enhance inference as an inverse simulation problem., Comment: Review article based on a Snowmass 2021 contribution

$\textbf{Aims}$. With the ever-increasing survey speed of optical wide-field telescopes and the importance of discovering transients when they are still young, rapid and reliable source localization is paramount. We present AutoSourceID-Light (ASID-L), an innovative framework that uses computer vision techniques that can naturally deal with large amounts of data and rapidly localize sources in optical images. $\textbf{Methods}$. We show that the AutoSourceID-Light algorithm based on U-shaped networks and enhanced with a Laplacian of Gaussian filter (Chen et al. 1987) enables outstanding performances in the localization of sources. A U-Net (Ronneberger et al. 2015) network discerns the sources in the images from many different artifacts and passes the result to a Laplacian of Gaussian filter that then estimates the exact location. $\textbf{Results}$. Application on optical images of the MeerLICHT telescope demonstrates the great speed and localization power of the method. We compare the results with the widely used SExtractor (Bertin & Arnouts 1996) and show the out-performances of our method. AutoSourceID-Light rapidly detects more sources not only in low and mid crowded fields, but particularly in areas with more than 150 sources per square arcminute.

We propose a new method to define anomaly scores and apply this to particle physics collider events. Anomalies can be either rare, meaning that these events are a minority in the normal dataset, or different, meaning they have values that are not inside the dataset. We quantify these two properties using an ensemble of One-Class Deep Support Vector Data Description models, which quantifies differentness, and an autoregressive flow model, which quantifies rareness. These two parameters are then combined into a single anomaly score using different combination algorithms. We train the models using a dataset containing only simulated collisions from the Standard Model of particle physics and test it using various hypothetical signals in four different channels and a secret dataset where the signals are unknown to us. The anomaly detection method described here has been evaluated in a summary paper [1] where it performed very well compared to a large number of other methods. The method is simple to implement and is applicable to other datasets in other fields as well., Comment: 7 pages, 5 figures; improved text

At GeV energies, the sky is dominated by the interstellar emission from the Galaxy. With limited statistics and spatial resolution, accurately separating point sources is therefore challenging. Here we present the first application of deep learning based algorithms to automatically detect and classify point sources from gamma-ray data. For concreteness we refer to this approach as AutoSourceID. To detect point sources, we utilized U-shaped convolutional networks for image segmentation and {\it k}-means for source clustering and localization. We also explored the Centroid-Net algorithm, which is designed to find and count objects. The training data are based on 9.5 years of exposure from The Fermi Large Area Telescope (Fermi-LAT) and we used source properties of active galactic nuclei (AGNs) and pulsars (PSRs) from the fourth Fermi-LAT source catalog (4FGL) in addition to several models of background interstellar emission. The results of the localization algorithm are fed into a classification neural network that is trained to separate the three general source classes (AGNs, PSRs, and FAKE sources). We compared our localization algorithms qualitatively with traditional methods and find them to have similar detection thresholds. We also demonstrate the robustness of our source localization algorithms to modifications in the interstellar emission models, which presents a clear advantage over traditional methods. The classification network is able to discriminate between the three classes with typical accuracy of $\sim$ 70%, as long as balanced data sets are used in classification training. In https://github.com/bapanes/AutoSourceID, we publish our training data sets and analysis scripts and invite the community to join the data challenge aimed to improve the localization and classification of gamma-ray point sources., Comment: 19 pages (two-column format), 14 figures, accepted for publication in Astronomy & Astrophysics

Optimisation problems are ubiquitous in particle and astrophysics, and involve locating the optimum of a complicated function of many parameters that may be computationally expensive to evaluate. We describe a number of global optimisation algorithms that are not yet widely used in particle astrophysics, benchmark them against random sampling and existing techniques, and perform a detailed comparison of their performance on a range of test functions. These include four analytic test functions of varying dimensionality, and a realistic example derived from a recent global fit of weak-scale supersymmetry. Although the best algorithm to use depends on the function being investigated, we are able to present general conclusions about the relative merits of random sampling, Differential Evolution, Particle Swarm Optimisation, the Covariance Matrix Adaptation Evolution Strategy, Bayesian Optimisation, Grey Wolf Optimisation, and the PyGMO Artificial Bee Colony, Gaussian Particle Filter and Adaptive Memory Programming for Global Optimisation algorithms., Comment: Experimental framework publicly available at http://www.github.com/darkmachines/high-dimensional-sampling

In this talk, we discuss the physics modelling of particle spectra arising from dark matter (DM) annihilation or decay. In the context of the indirect searches of DM, the final state products will, in general, undergo a set of complicated processes such as resonance decays, QED/QCD radiation, hadronisation and hadron decays. This set of processes lead to stable particles (photons, positrons, anti-protons, and neutrinos among others) which travel for very long distances before reaching the detectors. The modelling of their spectra contains some uncertainties which are often neglected in the relevant analyses. We discuss the sources of these uncertainties and estimate their impact on photon energy spectra for benchmark DM scenarios with $m_\chi \in [10, 1000]\,$GeV. Instructions for how to retrieve complete tables from Zenodo are also provided., Comment: 14 pages, 5 figures, Contribution for the talk presented in Tools for High Energy Physics and Cosmology workshop (TOOLS2020), 2-6 November 2020, IP2I Lyon, France. Based on arXiv:1812.07424. Relevant data can be found on Zenodo (http://doi.org/10.5281/zenodo.3764809)

In this manuscript, we propose to expand the use of scientific repositories such as Zenodo and HEP Data, in particular in order to better examine multi-parametric solutions of physical models. The implementation of interactive web-based visualizations enables fast and comfortable re-analysis and comparisons of phenomenological data. In order to illustrate our point of view, we present some examples and demos for dark matter models, supersymmetry exclusions and LHC simulations.

The lack of evidence for new physics at the Large Hadron Collider so far has prompted the development of model-independent search techniques. In this study, we compare the anomaly scores of a variety of anomaly detection techniques: an isolation forest, a Gaussian mixture model, a static autoencoder, and a $\beta$-variational autoencoder (VAE), where we define the reconstruction loss of the latter as a weighted combination of regression and classification terms. We apply these algorithms to the 4-vectors of simulated LHC data, but also investigate the performance when the non-VAE algorithms are applied to the latent space variables created by the VAE. In addition, we assess the performance when the anomaly scores of these algorithms are combined in various ways. Using supersymmetric benchmark points, we find that the logical AND combination of the anomaly scores yielded from algorithms trained in the latent space of the VAE is the most effective discriminator of all methods tested.

Traditional simulations on High-Performance Computing (HPC) systems typically involve modeling very large domains and/or very complex equations. HPC systems allow running large models, but limits in performance increase that have become more prominent in the last 5-10 years will likely be experienced. Therefore new approaches are needed to increase application performance. Deep learning appears to be a promising way to achieve this. Recently deep learning has been employed to enhance solving problems that traditionally are solved with large-scale numerical simulations using HPC. This type of application, deep learning for high-performance computing, is the theme of this whitepaper. Our goal is to provide concrete guidelines to scientists and others that would like to explore opportunities for applying deep learning approaches in their own large-scale numerical simulations. These guidelines have been extracted from a number of experiments that have been undertaken in various scientific domains over the last two years, and which are described in more detail in the Appendix. Additionally, we share the most important lessons that we have learned., Comment: White paper consists of 36 pages and 15 figures

In this paper, we minimize and compare two different fine-tuning measures in four high-scale supersymmetric models that are embedded in the MSSM. In addition, we determine the impact of current and future dark matter direct detection and collider experiments on the fine-tuning. We then compare the low-scale electroweak measure with the high-scale Barbieri-Giudice measure, which generally do not agree. However, we find that they do reduce to the same value when the higgsino parameter drives the degree of fine-tuning. Depending on the high-scale model and fine-tuning definition, we find a minimal fine-tuning of $3-38$ (corresponding to $\mathcal{O}(10-1)\%$) for the low-scale measure, and $63-571$ (corresponding to $\mathcal{O}(1-0.1)\%$) for the high-scale measure. In addition, minimally fine-tuned spectra give rise to a dark matter relic density that is between $10^{-3} < \Omega h^2 < 1$, when $\mu$ determines the minimum of the fine-tuning. We stress that it is too early to conclude on the fate of supersymmetry, based only on the fine-tuning paradigm.

Constraining the parameters of physical models with $>5-10$ parameters is a widespread problem in fields like particle physics and astronomy. The generation of data to explore this parameter space often requires large amounts of computational resources. The commonly used solution of reducing the number of relevant physical parameters hampers the generality of the results. In this paper we show that this problem can be alleviated by the use of active learning. We illustrate this with examples from high energy physics, a field where simulations are often expensive and parameter spaces are high-dimensional. We show that the active learning techniques query-by-committee and query-by-dropout-committee allow for the identification of model points in interesting regions of high-dimensional parameter spaces (e.g. around decision boundaries). This makes it possible to constrain model parameters more efficiently than is currently done with the most common sampling algorithms and to train better performing machine learning models on the same amount of data. Code implementing the experiments in this paper can be found on GitHub.

We present a study for the generation of events from a physical process with deep generative models. The simulation of physical processes requires not only the production of physical events, but also to ensure these events occur with the correct frequencies. We investigate the feasibility of learning the event generation and the frequency of occurrence with Generative Adversarial Networks (GANs) and Variational Autoencoders (VAEs) to produce events like Monte Carlo generators. We study three processes: a simple two-body decay, the processes $e^+e^-\to Z \to l^+l^-$ and $p p \to t\bar{t} $ including the decay of the top quarks and a simulation of the detector response. We find that the tested GAN architectures and the standard VAE are not able to learn the distributions precisely. By buffering density information of encoded Monte Carlo events given the encoder of a VAE we are able to construct a prior for the sampling of new events from the decoder that yields distributions that are in very good agreement with real Monte Carlo events and are generated several orders of magnitude faster. Applications of this work include generic density estimation and sampling, targeted event generation via a principal component analysis of encoded ground truth data, anomaly detection and more efficient importance sampling, e.g. for the phase space integration of matrix elements in quantum field theories., Comment: 22 pages, 9 figures

Motivated by the recent galactic center gamma-ray excess identified in the Fermi-LAT data, we perform a detailed study of QCD fragmentation uncertainties in the modeling of the energy spectra of gamma-rays from Dark-Matter (DM) annihilation. When Dark-Matter particles annihilate to coloured final states, either directly or via decays such as $W^{(*)}\to q\bar{q}'$, photons are produced from a complex sequence of shower, hadronisation and hadron decays. In phenomenological studies, their energy spectra are typically computed using Monte Carlo event generators. These results have however intrinsic uncertainties due to the specific model used and the choice of model parameters, which are difficult to asses and which are typically neglected. We derive a new set of hadronisation parameters (tunes) for the \textsc{Pythia~8.2} Monte Carlo generator from a fit to LEP and SLD data at the $Z$ peak. For the first time, we also derive a conservative set of uncertainties on the shower and hadronisation model parameters. Their impact on the gamma-ray energy spectra is evaluated and discussed for a range of DM masses and annihilation channels. The spectra and their uncertainties are also provided in tabulated form for future use. The fragmentation-parameter uncertainties may be useful for collider studies as well., Comment: Tables for the particle spectra with their uncertainties can be found in https://zenodo.org/record/3764809#.XqvK0y-cZQJ; Discussion slightly improved, and some references are added. Comments are welcome!

We present a deep learning solution to the prediction of particle production cross sections over a complicated, high-dimensional parameter space. We demonstrate the applicability by providing state-of-the-art predictions for the production of charginos and neutralinos at the Large Hadron Collider (LHC) at the next-to-leading order in the phenomenological MSSM-19 and explicitly demonstrate the performance for $pp\to\tilde{\chi}^+_1\tilde{\chi}^-_1,$ $\tilde{\chi}^0_2\tilde{\chi}^0_2$ and $\tilde{\chi}^0_2\tilde{\chi}^\pm_1$ as a proof of concept which will be extended to all SUSY electroweak pairs. We obtain errors that are lower than the uncertainty from scale and parton distribution functions with mean absolute percentage errors of well below $0.5\,\%$ allowing a safe inference at the next-to-leading order with inference times that improve the Monte Carlo integration procedures that have been available so far by a factor of $\mathcal{O}(10^7)$ from $\mathcal{O}(\rm{min})$ to $\mathcal{O}(\mu\rm{s})$ per evaluation., Comment: 7 pages, 3 figures

SPOT is an open source and free visual data analytics tool for multi-dimensional data-sets. Its web-based interface allows a quick analysis of complex data interactively. The operations on data such as aggregation and filtering are implemented. The generated charts are responsive and OpenGL supported. It follows FAIR principles to allow reuse and comparison of the published data-sets. The software also support PostgreSQL database for scalability.

The Fermi Collaboration has recently updated their analysis of gamma rays from the center of the Galaxy. They reconfirm the presence of an unexplained emission feature which is most prominent in the region of $1-10$ GeV, known as the Galactic Center GeV excess (GCE). Although the GCE is now firmly detected, an interpretation of this emission as a signal of self-annihilating dark matter (DM) particles is not unambiguously possible due to systematic effects in the gamma-ray modeling estimated in the Galactic Plane. In this paper we build a covariance matrix, collecting different systematic uncertainties investigated in the Fermi Collaboration's paper that affect the GCE spectrum. We show that models where part of the GCE is due to annihilating DM can still be consistent with the new data. We also re-evaluate the parameter space regions of the minimal supersymmetric Standard Model (MSSM) that can contribute dominantly to the GCE via neutralino DM annihilation. All recent constraints from DM direct detection experiments such as PICO, LUX, PandaX and Xenon1T, limits on the annihilation cross section from dwarf spheroidal galaxies and Large Hadron Collider limits are considered in this analysis. Due to a slight shift in the energy spectrum of the GC excess with respect to the previous Fermi analysis, and the recent limits from direct detection experiments, we find a slightly shifted parameter region of the MSSM compared to our previous analysis that is consistent with the GCE. Neutralinos with a mass between $85-220$ GeV can describe the excess via annihilation into a pair of $W$-bosons or top quarks. Remarkably, there are low fine-tuning models among the regions that we have found. The complete set of solutions will be probed by upcoming direct detection experiments and with dedicated searches in the upcoming data of the Large Hadron Collider., Comment: Replaced with the published version in JCAP. Auxiliary material on the construction of the covariance matrix and the GCE spectrum is included

We present a new method to interpret the $\gamma$-ray data of our inner Galaxy as measured by the Fermi Large Area Telescope (Fermi LAT). We train and test convolutional neural networks with simulated Fermi-LAT images based on models tuned to real data. We use this method to investigate the origin of an excess emission of GeV $\gamma$-rays seen in previous studies. Interpretations of this excess include $\gamma$ rays created by the annihilation of dark matter particles and $\gamma$ rays originating from a collection of unresolved point sources, such as millisecond pulsars. Our new method allows precise measurements of the contribution and properties of an unresolved population of $\gamma$-ray point sources in the interstellar diffuse emission model., Comment: 24 pages, 11 figures

We identify the parameter regions of the phenomenological minimal supersymmetric standard model (pMSSM) with the minimal possible fine-tuning. We show that the fine-tuning of the pMSSM is not large, nor under pressure by LHC searches. Low sbottom, stop and gluino masses turn out to be less relevant for low fine-tuning than commonly assumed. We show a link between low fine-tuning and the dark matter relic density. Fine-tuning arguments point to models with a dark matter candidate yielding the correct dark matter relic density: a bino-higgsino particle with a mass of $35-155$ GeV. Some of these candidates are compatible with recent hints seen in astrophysics experiments such as Fermi-LAT and AMS-02. We argue that upcoming direct search experiments, such as XENON1T, will test all of the most natural solutions in the next few years due to the sensitivity of these experiments on the spin-dependent WIMP-nucleon cross section., Comment: 8 pages, 5 figures, published version

A key research question at the Large Hadron Collider (LHC) is the test of models of new physics. Testing if a particular parameter set of such a model is excluded by LHC data is a challenge: It requires the time consuming generation of scattering events, the simulation of the detector response, the event reconstruction, cross section calculations and analysis code to test against several hundred signal regions defined by the ATLAS and CMS experiment. In the BSM-AI project we attack this challenge with a new approach. Machine learning tools are thought to predict within a fraction of a millisecond if a model is excluded or not directly from the model parameters. A first example is SUSY-AI, trained on the phenomenological supersymmetric standard model (pMSSM). About 300,000 pMSSM model sets - each tested with 200 signal regions by ATLAS - have been used to train and validate SUSY-AI. The code is currently able to reproduce the ATLAS exclusion regions in 19 dimensions with an accuracy of at least 93 percent. It has been validated further within the constrained MSSM and a minimal natural supersymmetric model, again showing high accuracy. SUSY-AI and its future BSM derivatives will help to solve the problem of recasting LHC results for any model of new physics. SUSY-AI can be downloaded at http://susyai.hepforge.org/. An on-line interface to the program for quick testing purposes can be found at http://www.susy-ai.org/.

Global fit studies performed in the pMSSM and the photon excess signal originating from the Galactic Center seem to suggest compressed electroweak supersymmetric spectra with a $\sim$100 GeV bino-like dark matter particle. We find that these scenarios are not probed by traditional electroweak supersymmetry searches at the LHC. We propose to extend the ATLAS and CMS electroweak supersymmetry searches with an improved strategy for bino-like dark matter, focusing on chargino plus next-to-lightest neutralino production, with a subsequent decay into a tri-lepton final state. We explore the sensitivity for pMSSM scenarios with $\Delta m = m_{\rm NLSP} - m_{\rm LSP} \sim (5 - 50)$ GeV in the $\sqrt{s} = 14$ TeV run of the LHC. Counterintuitively, we find that the requirement of low missing transverse energy increases the sensitivity compared to the current ATLAS and CMS searches. With 300 fb$^{-1}$ of data we expect the LHC experiments to be able to discover these supersymmetric spectra with mass gaps down to $\Delta m \sim 9$ GeV for DM masses between 40 and 140 GeV. We stress the importance of a dedicated search strategy that targets precisely these favored pMSSM spectra., Comment: Published in JHEP

We present a new global fit of the 19-dimensional phenomenological Minimal Supersymmetric Standard Model (pMSSM-19) that comply with all the latest experimental results from dark matter indirect, direct and accelerator dark matter searches. We show that the model provides a satisfactory explanation of the excess of gamma-rays from the Galactic centre observed by the Fermi~Large Area Telescope, assuming that it is produced by the annihilation of neutralinos in the Milky Way halo. We identify two regions that pass all the constraints: the first corresponds to neutralinos with a mass ~80-100 GeV annihilating into WW with a branching ratio of 95% ; the second to heavier neutralinos, with mass ~180-200 GeV annihilating into t tbar with a branching ratio of 87%. We show that neutralinos compatible with the Galactic centre GeV excess will soon be within the reach of LHC run-II -- notably through searches for charginos and neutralinos, squarks and light smuons -- and of Xenon1T, thanks to its unprecedented sensitivity to spin-dependent cross-section off neutrons., Comment: Minor changes following referee reports. Main conclusions unchanged. Matches version published in JCAP

Observations with the Fermi Large Area Telescope (LAT) indicate a possible small photon signal originating from the dwarf galaxy Reticulum II that exceeds the expected background between 2 GeV and 10 GeV. We have investigated two specific scenarios for annihilating WIMP dark matter within the phenomenological Minimal Supersymmetric Standard Model (pMSSM) framework as a possible source for these photons. We find that the same parameter ranges in pMSSM as reported by an earlier paper to be consistent with the Galactic center excess, is also consistent with the excess observed in Reticulum II, resulting in a J-factor of $\log_{10}(J(\alpha_{int}=0.5 deg)) \simeq (20.3-20.5)^{+0.2}_{-0.3}$. This J-factor is consistent with $\log_{10}(J(\alpha_{int}=0.5 deg)) = 19.5^{+1.0}_{-0.6}$ GeV$^2$cm$^{-5}$, which is derived using an optimized spherical Jeans analysis of kinematic data obtained from the Michigan/Magellan Fiber System (M2FS)., Comment: 4 pages, 2 figures, accepted in JCAP

Observations with the Fermi Large Area Telescope (LAT) indicate an excess in gamma rays originating from the center of our Galaxy. A possible explanation for this excess is the annihilation of Dark Matter particles. We have investigated the annihilation of neutralinos as Dark Matter candidates within the phenomenological Minimal Supersymmetric Standard Model (pMSSM). An iterative particle filter approach was used to search for solutions within the pMSSM. We found solutions that are consistent with astroparticle physics and collider experiments, and provide a fit to the energy spectrum of the excess. The neutralino is a Bino/Higgsino or Bino/Wino/Higgsino mixture with a mass in the range $84-92$~GeV or $87-97$~GeV annihilating into W bosons. A third solutions is found for a neutralino of mass $174-187$~GeV annihilating into top quarks. The best solutions yield a Dark Matter relic density $0.06 < \Omega h^2 <0.13$. These pMSSM solutions make clear forecasts for LHC, direct and indirect DM detection experiments. If the MSSM explanation of the excess seen by Fermi-LAT is correct, a DM signal might be discovered soon., Comment: Large extension of previous paper: 2 more solutions found in the MSSM (Bino-Higgsino, Bino-Wino-Higgsino into WW and Bino into ttbar), added description on extra fit uncertainties, added description on flavor observables, added discussion on dwarf limits

We have systematically investigated the production of a Higgs boson with a mass of about $125$ GeV in the decays of supersymmetric particles within the phenomenological MSSM (pMSSM). We find regions of parameter space that are consistent with all world data and that predict a sizeable rate of anomalous Higgs, di-Higgs and even tri-Higgs events at the 14 TeV LHC. All relevant SUSY production processes are investigated. We find that Higgs bosons can be produced in a large variety of SUSY processes, resulting in a large range of different detector signatures containing missing transverse momentum. Such Higgs events are outstanding signatures for new physics already for the early 14 TeV LHC data. SUSY processes are also important to interprete deviations found in upcoming Standard Model Higgs and di-Higgs production measurements., Comment: Version submitted to JHEP

We study the impact of the jets and missing transverse momentum SUSY analyses of the ATLAS experiment on the phenomenological MSSM (pMSSM). We investigate sets of SUSY models with a flat and logarithmic prior in the SUSY mass scale and a mass range up to 1 and 3 TeV, respectively. These models were found previously in the study 'Supersymmetry without Prejudice'. Removing models with long-lived SUSY particles, we show that 99% of 20000 randomly generated pMSSM model points with a flat prior and 87% for a logarithmic prior are excluded by the ATLAS results. For models with squarks and gluinos below 600 GeV all models of the pMSSM grid are excluded. We identify SUSY spectra where the current ATLAS search strategy is less sensitive and propose extensions to the inclusive jets search channel.

So far, no supersymmetric particles have been detected at the Large Hadron Collider (LHC). However, the recent Higgs results have interesting implications for the SUSY parameter space. In this paper, we study the consequences of an LHC Higgs signal for a model with non-universal gaugino masses in the context of SU(5) unification. The gaugino mass ratios associated with the higher representations produce viable spectra that are largely inaccessible to the current LHC and direct dark matter detection experiments. Thus, in light of the Higgs results, the non-observation of SUSY is no surprise., Comment: supplementary file containing plots with log priors in ancillary files. v2: added some comments on more general settings and references, accepted for publication in JHEP

First ATLAS searches for signals of Supersymmetry in proton-proton collisions at the LHC are presented. These searches are performed in various channels containing different lepton and jet multiplicities in the final states; the full data sample recorded in the 2010 LHC run, corresponding to an integrated luminosity of 35 $\rm{pb}^{-1}$, has been analysed. The limits on squarks and gluinos are the most stringent to date., Comment: Prepared for the 46th Rencontres de Moriond on Electroweak Interactions and Unified Theories, La Thuile, Italy, 13 - 20 Mar 2011

We present the recent prospects and strategies of the ATLAS and CMS experiments in the search for signatures of Supersymmetry. The emphasis is placed on standard signatures with missing transverse energy and on the initial data taking periods., Comment: Parallel talk at ICHEP08, Philadelphia, USA, July 2008. 4 pages

Distributions from high-pT HERA event data analyzed in a general search for new physics at H1 have been incorporated into Quaero, an algorithm designed to automate tests of specific hypotheses with high energy collider data. The use of Quaero@H1 to search for leptoquarks, R-parity violating supersymmetry, and excited quarks provides examples to develop intuition for the algorithm's performance., Comment: Submitted to Eur. Phys. J. C

The H1 and ZEUS collaborations have used the HERA I data to search for physics beyond the SM and to test electroweak physics in electron-proton collisions. The new period of data taking (HERA II) has started and first HERA II analyses become available. An overview is given of recent highlights, including isolated lepton and multi-lepton events, tests of models of new physics and the first systematic search for deviations in all final states with high transverse momentum., Comment: Talk given at Recontres du VietNam 2004

This article reviews recent jet physics results from HERA.