Temporal light field reconstruction for rendering distribution effects

Traditionally, effects that require evaluating multidimensional integrals for each pixel, such as motion blur, depth of field, and soft shadows, suffer from noise due to the variance of the high-dimensional integrand. In this paper, we describe a general reconstruction technique that exploits the anisotropy in the temporal light field and permits efficient reuse of samples between pixels, multiplying the effective sampling rate by a large factor. We show that our technique can be applied in situations that are challenging or impossible for previous anisotropic reconstruction methods, and that it can yield good results with very sparse inputs. We demonstrate our method for simultaneous motion blur, depth of field, and soft shadows

Decoupling algorithms from schedules for easy optimization of image processing pipelines

Using existing programming tools, writing high-performance image processing code requires sacrificing readability, portability, and modularity. We argue that this is a consequence of conflating what computations define the algorithm, with decisions about storage and the order of computation. We refer to these latter two concerns as the schedule, including choices of tiling, fusion, recomputation vs. storage, vectorization, and parallelism. We propose a representation for feed-forward imaging pipelines that separates the algorithm from its schedule, enabling high-performance without sacrificing code clarity. This decoupling simplifies the algorithm specification: images and intermediate buffers become functions over an infinite integer domain, with no explicit storage or boundary conditions. Imaging pipelines are compositions of functions. Programmers separately specify scheduling strategies for the various functions composing the algorithm, which allows them to efficiently explore different optimizations without changing the algorithmic code. We demonstrate the power of this representation by expressing a range of recent image processing applications in an embedded domain specific language called Halide, and compiling them for ARM, x86, and GPUs. Our compiler targets SIMD units, multiple cores, and complex memory hierarchies. We demonstrate that it can handle algorithms such as a camera raw pipeline, the bilateral grid, fast local Laplacian filtering, and image segmentation. The algorithms expressed in our language are both shorter and faster than state-of-the-art implementations.National Science Foundation (U.S.) (Grant 0964004)National Science Foundation (U.S.) (Grant 0964218)National Science Foundation (U.S.) (Grant 0832997)United States. Dept. of Energy (Award DE-SC0005288)Cognex CorporationAdobe System

Evidence for proton acceleration up to TeV energies based on VERITAS and Fermi-LAT observations of the Cas A SNR

We present a study of $\gamma$-ray emission from the core-collapse supernova remnant Cas~A in the energy range from 0.1GeV to 10TeV. We used 65 hours of VERITAS data to cover 200 GeV - 10 TeV, and 10.8 years of \textit{Fermi}-LAT data to cover 0.1-500 GeV. The spectral analysis of \textit{Fermi}-LAT data shows a significant spectral curvature around $1.3 \pm 0.4_{stat}$ GeV that is consistent with the expected spectrum from pion decay. Above this energy, the joint spectrum from \textit{Fermi}-LAT and VERITAS deviates significantly from a simple power-law, and is best described by a power-law with spectral index of $2.17\pm 0.02_{stat}$ with a cut-off energy of $2.3 \pm 0.5_{stat}$ TeV. These results, along with radio, X-ray and $\gamma$-ray data, are interpreted in the context of leptonic and hadronic models. Assuming a one-zone model, we exclude a purely leptonic scenario and conclude that proton acceleration up to at least 6 TeV is required to explain the observed $\gamma$-ray spectrum. From modeling of the entire multi-wavelength spectrum, a minimum magnetic field inside the remnant of $B_{\mathrm{min}}\approx150\,\mathrm{\mu G}$ is deduced.Comment: 33 pages, 9 Figures, 6 Table

Measurement of Cosmic-ray Electrons at TeV Energies by VERITAS

Cosmic-ray electrons and positrons (CREs) at GeV-TeV energies are a unique probe of our local Galactic neighborhood. CREs lose energy rapidly via synchrotron radiation and inverse-Compton scattering processes while propagating within the Galaxy and these losses limit their propagation distance. For electrons with TeV energies, the limit is on the order of a kiloparsec. Within that distance there are only a few known astrophysical objects capable of accelerating electrons to such high energies. It is also possible that the CREs are the products of the annihilation or decay of heavy dark matter (DM) particles. VERITAS, an array of imaging air Cherenkov telescopes in southern Arizona, USA, is primarily utilized for gamma-ray astronomy, but also simultaneously collects CREs during all observations. We describe our methods of identifying CREs in VERITAS data and present an energy spectrum, extending from 300 GeV to 5 TeV, obtained from approximately 300 hours of observations. A single power-law fit is ruled out in VERITAS data. We find that the spectrum of CREs is consistent with a broken power law, with a break energy at 710 $\pm$ 40$_{stat}$ $\pm$ 140$_{syst}$ GeV.Comment: 17 pages, 2 figures, accepted for publication in PR

Direct measurement of stellar angular diameters by the VERITAS Cherenkov Telescopes

The angular size of a star is a critical factor in determining its basic properties. Direct measurement of stellar angular diameters is difficult: at interstellar distances stars are generally too small to resolve by any individual imaging telescope. This fundamental limitation can be overcome by studying the diffraction pattern in the shadow cast when an asteroid occults a star, but only when the photometric uncertainty is smaller than the noise added by atmospheric scintillation. Atmospheric Cherenkov telescopes used for particle astrophysics observations have not generally been exploited for optical astronomy due to the modest optical quality of the mirror surface. However, their large mirror area makes them well suited for such high-time-resolution precision photometry measurements. Here we report two occultations of stars observed by the VERITAS Cherenkov telescopes with millisecond sampling, from which we are able to provide a direct measurement of the occulted stars' angular diameter at the $\leq0.1$ milliarcsecond scale. This is a resolution never achieved before with optical measurements and represents an order of magnitude improvement over the equivalent lunar occultation method. We compare the resulting stellar radius with empirically derived estimates from temperature and brightness measurements, confirming the latter can be biased for stars with ambiguous stellar classifications.Comment: Accepted for publication in Nature Astronom

Discovery of very-high-energy emission from RGB J2243+203 and derivation of its redshift upper limit

Very-high-energy (VHE; $>$ 100 GeV) gamma-ray emission from the blazar RGB J2243+203 was discovered with the VERITAS Cherenkov telescope array, during the period between 21 and 24 December 2014. The VERITAS energy spectrum from this source can be fit by a power law with a photon index of $4.6 \pm 0.5$, and a flux normalization at 0.15 TeV of $(6.3 \pm 1.1) \times 10^{-10} ~ \textrm{cm}^{-2} \textrm{s}^{-1} \textrm{TeV}^{-1}$. The integrated \textit{Fermi}-LAT flux from 1 GeV to 100 GeV during the VERITAS detection is $(4.1 \pm 0.8) \times 10^{\textrm{-8}} ~\textrm{cm}^{\textrm{-2}}\textrm{s}^{\textrm{-1}}$, which is an order of magnitude larger than the four-year-averaged flux in the same energy range reported in the 3FGL catalog, ($4.0 \pm 0.1 \times 10^{\textrm{-9}} ~ \textrm{cm}^{\textrm{-2}}\textrm{s}^{\textrm{-1}}$). The detection with VERITAS triggered observations in the X-ray band with the \textit{Swift}-XRT. However, due to scheduling constraints \textit{Swift}-XRT observations were performed 67 hours after the VERITAS detection, not simultaneous with the VERITAS observations. The observed X-ray energy spectrum between 2 keV and 10 keV can be fitted with a power-law with a spectral index of $2.7 \pm 0.2$, and the integrated photon flux in the same energy band is $(3.6 \pm 0.6) \times 10^{-13} ~\textrm{cm}^{-2} \textrm{s}^{-1}$. EBL model-dependent upper limits of the blazar redshift have been derived. Depending on the EBL model used, the upper limit varies in the range from z $<~0.9$ to z $<~1.1$

Dark Matter Constraints from a Joint Analysis of Dwarf Spheroidal Galaxy Observations with VERITAS

We present constraints on the annihilation cross section of WIMP dark matter based on the joint statistical analysis of four dwarf galaxies with VERITAS. These results are derived from an optimized photon weighting statistical technique that improves on standard imaging atmospheric Cherenkov telescope (IACT) analyses by utilizing the spectral and spatial properties of individual photon events. We report on the results of $\sim$230 hours of observations of five dwarf galaxies and the joint statistical analysis of four of the dwarf galaxies. We find no evidence of gamma-ray emission from any individual dwarf nor in the joint analysis. The derived upper limit on the dark matter annihilation cross section from the joint analysis is $1.35\times 10^{-23} {\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the bottom quark ($b\bar{b}$) final state, $2.85\times 10^{-24}{\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the tau lepton ($\tau^{+}\tau^{-}$) final state and $1.32\times 10^{-25}{\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the gauge boson ($\gamma\gamma$) final state.Comment: 14 pages, 9 figures, published in PRD, Ascii tables containing annihilation cross sections limits are available for download as ancillary files with readme.txt file description of limit

Very-high-energy observations of the binaries V 404 Cyg and 4U 0115+634 during giant X-ray outbursts

Transient X-ray binaries produce major outbursts in which the X-ray flux can increase over the quiescent level by factors as large as $10^7$. The low-mass X-ray binary V 404 Cyg and the high-mass system 4U 0115+634 underwent such major outbursts in June and October 2015, respectively. We present here observations at energies above hundreds of GeV with the VERITAS observatory taken during some of the brightest X-ray activity ever observed from these systems. No gamma-ray emission has been detected by VERITAS in 2.5 hours of observations of the microquasar V 404 Cyg from 2015, June 20-21. The upper flux limits derived from these observations on the gamma-ray flux above 200 GeV of F $< 4.4\times 10^{-12}$ cm$^{-2}$ s$^{-1}$ correspond to a tiny fraction (about $10^{-6}$) of the Eddington luminosity of the system, in stark contrast to that seen in the X-ray band. No gamma rays have been detected during observations of 4U 0115+634 in the period of major X-ray activity in October 2015. The flux upper limit derived from our observations is F $< 2.1\times 10^{-12}$ cm$^{-2}$ s$^{-1}$ for gamma rays above 300 GeV, setting an upper limit on the ratio of gamma-ray to X-ray luminosity of less than 4%.Comment: Accepted for publication in the Astrophysical Journa