Improvement in Fast Particle Track Reconstruction with Robust Statistics

ArXiv 1308.5501 (2013)

Authors:

MG Aartsen, R Abbasi, Y Abdou, M Ackermann, J Adams, JA Aguilar, M Ahlers, D Altmann, J Auffenberg, X Bai, M Baker, SW Barwick, V Baum, R Bay, JJ Beatty, S Bechet, J Becker Tjus, K-H Becker, ML Benabderrahmane, S BenZvi, P Berghaus, D Berley, E Bernardini, A Bernhard, DZ Besson, G Binder, D Bindig, M Bissok, E Blaufuss, J Blumenthal, DJ Boersma, S Bohaichuk, C Bohm, D Bose, S Böser, O Botner, L Brayeur, H-P Bretz, AM Brown, R Bruijn, J Brunner, M Carson, J Casey, M Casier, D Chirkin, A Christov, B Christy, K Clark, F Clevermann, S Coenders, S Cohen, DF Cowen, AH Cruz Silva, M Danninger, J Daughhetee, JC Davis, M Day, C De Clercq, S De Ridder, P Desiati, KD de Vries, M de With, T DeYoung, JC Díaz-Vélez, M Dunkman, R Eagan, B Eberhardt, J Eisch, S Euler, PA Evenson, O Fadiran, AR Fazely, A Fedynitch, J Feintzeig, T Feusels, K Filimonov, C Finley, T Fischer-Wasels, S Flis, A Franckowiak, K Frantzen, T Fuchs, TK Gaisser, J Gallagher, L Gerhardt, L Gladstone, T Glüsenkamp, A Goldschmidt, G Golup, JG Gonzalez, JA Goodman, D Góra, DT Grandmont, D Grant, A Groß, C Ha, A Haj Ismail, P Hallen, A Hallgren, F Halzen, K Hanson, D Heereman, D Heinen, K Helbing, R Hellauer, S Hickford, GC Hill, KD Hoffman, R Hoffmann, A Homeier, K Hoshina, W Huelsnitz, PO Hulth, K Hultqvist, S Hussain, A Ishihara, E Jacobi, J Jacobsen, K Jagielski, GS Japaridze, K Jero, O Jlelati, B Kaminsky, A Kappes, T Karg, A Karle, JL Kelley, J Kiryluk, J Kläs, SR Klein, J-H Köhne, G Kohnen, H Kolanoski, L Köpke, C Kopper, S Kopper, DJ Koskinen, M Kowalski, M Krasberg, K Krings, G Kroll, J Kunnen, N Kurahashi, T Kuwabara, M Labare, H Landsman, MJ Larson, M Lesiak-Bzdak, M Leuermann, J Leute, J Lünemann, O Macías, J Madsen, G Maggi, R Maruyama, K Mase, HS Matis, F McNally, K Meagher, M Merck, T Meures, S Miarecki, E Middell, N Milke, J Miller, L Mohrmann, T Montaruli, R Morse, R Nahnhauer, U Naumann, H Niederhausen, SC Nowicki, DR Nygren, A Obertacke, S Odrowski, A Olivas, A Omairat, A O'Murchadha, L Paul, JA Pepper, C Pérez de los Heros, C Pfendner, D Pieloth, E Pinat, J Posselt, PB Price, GT Przybylski, L Rädel, M Rameez, K Rawlins, C Ré, B Recht, P Redl, R Reimann, E Resconi, W Rhode, M Ribordy, M Richman, B Riedel, JP Rodrigues, C Rott, T Ruhe, B Ruzybayev, D Ryckbosch, SM Saba, T Salameh, H-G Sander, M Santander, S Sarkar, K Schatto, F Scheriau, T Schmidt, M Schmitz, S Schoenen, S Schöneberg, A Schönwald, A Schukraft, L Schulte, O Schulz, D Seckel, Y Sestayo, S Seunarine, R Shanidze, C Sheremata, MWE Smith, D Soldin, GM Spiczak, C Spiering, M Stamatikos, T Stanev, NA Stanisha, A Stasik, T Stezelberger, RG Stokstad, A Stößl, EA Strahler, R Ström, GW Sullivan, H Taavola, I Taboada, A Tamburro, A Tepe, S Ter-Antonyan, G Tešić, S Tilav, PA Toale, S Toscano, E Unger, M Usner, S Vallecorsa, N van Eijndhoven, A Van Overloop, J van Santen, M Vehring, M Voge, M Vraeghe, C Walck, T Waldenmaier, M Wallraff, Ch Weaver, M Wellons, C Wendt, S Westerhoff, N Whitehorn, K Wiebe, CH Wiebusch, DR Williams, H Wissing, M Wolf, TR Wood, K Woschnagg, DL Xu, XW Xu, JP Yanez, G Yodh, S Yoshida, P Zarzhitsky, J Ziemann, S Zierke, M Zoll

Abstract:

The IceCube project has transformed one cubic kilometer of deep natural Antarctic ice into a Cherenkov detector. Muon neutrinos are detected and their direction inferred by mapping the light produced by the secondary muon track inside the volume instrumented with photomultipliers. Reconstructing the muon track from the observed light is challenging due to noise, light scattering in the ice medium, and the possibility of simultaneously having multiple muons inside the detector, resulting from the large flux of cosmic ray muons. This manuscript describes work on two problems: (1) the track reconstruction problem, in which, given a set of observations, the goal is to recover the track of a muon; and (2) the coincident event problem, which is to determine how many muons are active in the detector during a time window. Rather than solving these problems by developing more complex physical models that are applied at later stages of the analysis, our approach is to augment the detectors early reconstruction with data filters and robust statistical techniques. These can be implemented at the level of on-line reconstruction and, therefore, improve all subsequent reconstructions. Using the metric of median angular resolution, a standard metric for track reconstruction, we improve the accuracy in the initial reconstruction direction by 13%. We also present improvements in measuring the number of muons in coincident events: we can accurately determine the number of muons 98% of the time.

Colliding clusters and dark matter self-interactions

ArXiv 1308.3419 (2013)

Authors:

Felix Kahlhoefer, Kai Schmidt-Hoberg, Mads T Frandsen, Subir Sarkar

Abstract:

When a dark matter halo moves through a background of dark matter particles, self-interactions can lead to both deceleration and evaporation of the halo and thus shift its centroid relative to the collisionless stars and galaxies. We study the magnitude and time evolution of this shift for two classes of dark matter self-interactions, viz. frequent self-interactions with small momentum transfer (e.g. due to long-range interactions) and rare self-interactions with large momentum transfer (e.g. contact interactions), and find important differences between the two cases. We find that neither effect can be strong enough to completely separate the dark matter halo from the galaxies, if we impose conservative bounds on the self-interaction cross-section. The majority of both populations remain bound to the same gravitational potential and the peaks of their distributions are therefore always coincident. Consequently any apparent separation is mainly due to particles which are leaving the gravitational potential, so will be largest shortly after the collision but not observable in evolved systems. Nevertheless the fraction of collisions with large momentum transfer is an important characteristic of self-interactions, which can potentially be extracted from observational data and provide an important clue as to the nature of dark matter.

Colliding clusters and dark matter self-interactions

(2013)

Authors:

Felix Kahlhoefer, Kai Schmidt-Hoberg, Mads T Frandsen, Subir Sarkar

Reconstruction of the primordial power spectrum of curvature perturbations using multiple data sets

ArXiv 1308.2317 (2013)

Authors:

Paul Hunt, Subir Sarkar

Abstract:

Detailed knowledge of the primordial power spectrum of curvature perturbations is essential both in order to elucidate the physical mechanism (`inflation') which generated it, and for estimating the cosmological parameters from observations of the cosmic microwave background and large-scale structure. Hence it ought to be extracted from such data in a model-independent manner, however this is difficult because relevant cosmological observables are given by a convolution of the primordial perturbations with some smoothing kernel which depends on both the assumed world model and the matter content of the universe. Moreover the deconvolution problem is ill-conditioned so a regularisation scheme must be employed to control error propagation. We demonstrate that `Tikhonov regularisation' can robustly reconstruct the primordial spectrum from multiple cosmological data sets, a significant advantage being that both its uncertainty and resolution are then quantified. Using Monte Carlo simulations we investigate several regularisation parameter selection methods and find that generalised cross-validation and Mallow's $C_p$ method give optimal results. We apply our inversion procedure to data from the Wilkinson Microwave Anisotropy Probe, other ground-based small angular scale CMB experiments, and the Sloan Digital Sky Survey. The reconstructed spectrum (assuming the standard $\Lambda$CDM cosmology) is \emph{not} scale-free but has an infrared cutoff at $k \lesssim 5 \times 10^{-4}\; \mathrm{Mpc}^{-1}$ (due to the anomalously low CMB quadrupole) and several features with $\sim 2 \sigma$ significance at $k/\mathrm{Mpc}^{-1} \sim$ 0.0013--0.0025, 0.0362--0.0402 and 0.051--0.056, reflecting the `WMAP glitches'. To test whether these are indeed real will require more accurate data, such as from the Planck satellite and new ground-based experiments.

Reconstruction of the primordial power spectrum of curvature perturbations using multiple data sets

(2013)

Authors:

Paul Hunt, Subir Sarkar