Advancements in Very-High-Energy Gamma-Ray Astronomy with Applications to the Study of Cosmic Rays

This work aims to contribute to the study of the origins of cosmic rays, and broadly, to the advancement of both data analysis methods and instrumentation for very-high-energy γ-ray astronomy. First, reviewing the state of γ-ray astronomy, we show how gains in sensitivity can be achieved through sophisticated data analyses and improved instrumental designs. We then develop such an improved analysis method for the Very Energetic Radiation Imaging Telescope Array System (VERITAS) by combining Image Template Method (ITM) with Boosted Decision Trees (BDT), and study its performance, attaining a 30-50% improvement in integral sensitivity over the instrument’s standard analysis. Systematic issues in spectral reconstruction that the analysis displays are resolved satisfactorily by imposing a more stringent condition on the selection of its energy threshold. We employ the newly developed analysis to measure the γ-ray energy spectrum of the starburst galaxy M82, and combining our result with a measurement from the Fermi Large Area Telescope (Fermi-LAT), we find that a single power law fits the spectrum well between 100 GeV and 10 TeV, with no evidence for a spectral break or a cutoff. We conclude that this is in line with the current understanding that M82 is not a good proton calorimeter. Finally, we detail the design, implementation, and performance of the optical alignment system of the prototypeSchwarzschild-Couder Telescope (pSCT) for the Cherenkov Telescope Array (CTA), a novel two-mirror design that addresses many shortcomings of current instruments

Computational efficiency of staggered Wilson fermions: A first look

Results on the computational efficiency of 2-flavor staggered Wilson fermions compared to usual Wilson fermions in a quenched lattice QCD simulation on $16^3\times32$ lattice at $\beta=6$ are reported. We compare the cost of inverting the Dirac matrix on a source by the conjugate gradient (CG) method for both of these fermion formulations, at the same pion masses, and without preconditioning. We find that the number of CG iterations required for convergence, averaged over the ensemble, is less by a factor of almost 2 for staggered Wilson fermions, with only a mild dependence on the pion mass. We also compute the condition number of the fermion matrix and find that it is less by a factor of 4 for staggered Wilson fermions. The cost per CG iteration, dominated by the cost of matrix-vector multiplication for the Dirac matrix, is known from previous work to be less by a factor 2-3 for staggered Wilson compared to usual Wilson fermions. Thus we conclude that staggered Wilson fermions are 4-6 times cheaper for inverting the Dirac matrix on a source in the quenched backgrounds of our study.Comment: v2: Major correction and revisions: we had overlooked a factor 1/4 in the cost estimate for matrix-vector multiplication with the staggered Wilson Dirac matrix. This gives an increased speed-up by a factor 4 for the overall computation cost. 7 pages, 3 figures, presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, German