PoS(ICRC2017)407. Instrumentation for sflash

Similar documents
[CRI299] SFLASH: Absolute Measurement of Fluorescence Yield from Shower Particles

Energy Spectrum of Ultra-High-Energy Cosmic Rays Measured by The Telescope Array

Search for EeV Protons of Galactic Origin

Burst Shower Events Observed by the Telescope Array Surface Detector

PoS(ICRC2015)349. TA Spectrum Summary

PoS(ICRC2015)657. Telescope Array extension: TA 4

FLASH. FLuorescence in Air from SHowers (SLAC E-165) Pisin Chen SLAC. Report to DOE HEP Review SLAC, June 2-4, 2004

Imaging and non-imaging Cherenkov hybrid reconstruction

The Cosmic Ray Air Fluorescence Fresnel lens Telescope (CRAFFT) for the next generation UHECR observatory

PoS(PD07)031. General performance of the IceCube detector and the calibration results

ARTICLE IN PRESS. Nuclear Instruments and Methods in Physics Research A

PoS(ICRC2017)548. Telescope Array anisotropy summary

Measurements of the Fluorescence Light Yield in Electromagnetic Showers

Search for isotropic microwave radiation from electron beam in the atmosphere

PoS(ICRC2015)641. Cloud Monitoring using Nitrogen Laser for LHAASO Experiment. Z.D. Sun 1,Y. Zhang 2,F.R. Zhu 1 for the LHAASO Collaboration

PoS(ICRC2015)635. The NICHE Array: Status and Plans. Douglas R Bergman

A measurement of the air fluorescence yield

arxiv:physics/ v1 3 Aug 2006

Study of Performance Improvement for Air Shower Array with Surface Water Cherenkov Detectors

SLAC-PUB September 2006

Development of a prototype for Fluorescence detector Array of Single-pixel Telescopes (FAST)

Measurement of nuclear recoil responses of NaI(Tl) crystal for dark matter search

PoS(ICRC2017)462. Cherenkov differential detector at the Yakutsk extensive air shower array

PoS(ICRC2017)945. In-ice self-veto techniques for IceCube-Gen2. The IceCube-Gen2 Collaboration

PoS(ICRC2015)326. TA anisotropy summary

John Ellison University of California, Riverside. Quarknet 2008 at UCR

Direct WIMP Detection in Double-Phase Xenon TPCs

XMASS 1.5, the next step of the XMASS experiment

Muon track reconstruction and veto performance with D-Egg sensor for IceCube-Gen2

Optic Detectors Calibration for Measuring Ultra-High Energy Extensive Air Showers Cherenkov Radiation by 532 nm Laser

Characteristics of Forbush decreases measured by means of the new scintillation muon hodoscope ScMH

An introduction to AERA. An Introduction to the Auger Engineering Radio Array (AERA)

Primary cosmic ray mass composition above 1 PeV as measured by the PRISMA-YBJ array

Timing calibration of the LHAASO-KM2A electromagnetic particle detectors

Interpretation of the energy spectrum observed with the Telescope Array surface detectors

AIRFLY: Measurement of the Air Fluorescence induced by electrons

Determination of parameters of cascade showers in the water calorimeter using 3D-distribution of Cherenkov light

arxiv: v1 [physics.ins-det] 3 Feb 2011

Mass Composition Study at the Pierre Auger Observatory

Accurate Measurement of the Cosmic Ray Proton Spectrum from 100TeV to 10PeV with LHAASO

Muon track reconstruction and veto performance with D-Egg sensor for IceCube-Gen2

GRAINE Project: a balloon-borne emulsion gamma-ray telescope

The path towards measuring Cosmic Ray Air Showers from Space

Precise Measurement of the Absolute Yield of Fluorescence Photons in Atmospheric Gases

Full author list:

Zero degree neutron energy spectra measured by LHCf at s = 13 TeV proton-proton collision

Visualization of Xe and Sn Atoms Generated from Laser-Produced Plasma for EUV Light Source

Implications for the Radio Detection of Cosmic Rays

Development of Cs 2 Te photocathode RF gun system for compact THz SASE-FEL

Neutron flux measurement using fast-neutron activation of 12 B and 12 N isotopes in hydrocarbonate scintillators

Nuclear Physics Laboratory. Gamma spectroscopy with scintillation detectors. M. Makek Faculty of Science Department of Physics

MEASUREMENT OF TEMPORAL RESOLUTION AND DETECTION EFFICIENCY OF X-RAY STREAK CAMERA BY SINGLE PHOTON IMAGES

Dependence of the GRAPES-3 EAS particle density and trigger rate on atmospheric pressure and temperature

GEANT4 simulation about the laboratory astrophysics in Taiwan

PoS(TIPP2014)033. Upgrade of MEG Liquid Xenon Calorimeter. Ryu SAWADA. ICEPP, the University of Tokyo

Depth of maximum of air-shower profiles at the Pierre Auger Observatory: Measurements above ev and Composition Implications

Nuclear Instruments and Methods in Physics Research A

PoS(ICRC2017)330. Influence of weather conditions in the Yakutsk array region on Cherenkov light flux measurements. Igor Petrov

Some Thoughts on Laboratory Astrophysics for UHE Cosmic Rays. Pierre Sokolsky University of Utah SABRE Workshop SLAC, March, 2006

Photoemission Properties of and Under Various Temperature and Incident Photon Energy Conditions

PoS(PhotoDet 2012)010

First results from the full-scale prototype for the Fluorescence detector Array of Single-pixel Telescopes

Scintillation Detector

arxiv: v1 [astro-ph] 30 Jul 2008

PoS(ICRC2017)522. Testing the agreement between the X max distributions measured by the Pierre Auger and Telescope Array Observatories

hν' Φ e - Gamma spectroscopy - Prelab questions 1. What characteristics distinguish x-rays from gamma rays? Is either more intrinsically dangerous?

A method of observing Cherenkov light at the Yakutsk EAS array

arxiv: v1 [astro-ph.he] 25 Mar 2015

Antifreeze design for Muon Detector of LHAASO

Measurement of air shower maxima and p-air cross section with the Telescope Array

Calibration of photo sensors for the space-based cosmic ray telescope JEM-EUSO

The air fluorescence yield. Fernando Arqueros Universidad Complutense de Madrid Spain

PoS(ICRC2017)1052. Overview and Performance of the Wavelength-shifting Optical Module (WOM) The IceCube-Gen2 collaboration

Muon detector and muon flux measurement at Yangyang underground laboratory for the COSINE-100 experiment

Studies of Hadron Calorimeter

Results of the Search for Ultra High-Energy Neutrinos with ANITA-II

The Fluorescence Detector of the Pierre Auger Observatory

Emittance and Quantum Efficiency Measurements from a 1.6 cell S- Band Photocathode RF Gun with Mg Cathode *

EAS spectrum in thermal neutrons measured with PRISMA-32

Air Fluorescence Photon Yield In Cosmic Ray Showers

Performance of the MCP-PMT for the Belle II TOP counter

The NUCLEON Space Experiment Preliminary Results. Skobeltsyn Institute of Nuclear Physics, Moscow State University, Moscow, , Russia

The Physics of Ultrahigh Energy Cosmic Rays. Example Poster Presentation Physics 5110 Spring 2009 Reminder: Posters are due Wed April 29 in class.

Cherenkov Detector. Cosmic Rays Cherenkov Detector. Lodovico Lappetito. CherenkovDetector_ENG - 28/04/2016 Pag. 1

IceCube: Ultra-high Energy Neutrinos

Distillation purification and radon assay of liquid xenon

High quantum efficiency S-20 photocathodes for photon counting applications

Analytical Chemistry II

Feasibility of antenna array experiment for Earth skimming tau-neutrino detection in Antarctica

PoS(ICRC2017)326. The influence of weather effects on the reconstruction of extensive air showers at the Pierre Auger Observatory

The ExaVolt Antenna (EVA)

Development of Analysis Method using GEANT4 for Cosmic Ray Radiography

Confocal Microscope Imaging of Single emitter fluorescence and Observing Photon Antibunching Using Hanbury Brown and Twiss setup. Lab.

Position-sensitive SiPM detector for separation of Cherenkov and fluorescent light of EAS

Chapter 4 Scintillation Detectors

UHE Cosmic Rays in the Auger Era

Ultra-High-Energy Cosmic Rays: A Tale of Two Observatories

Because light behaves like a wave, we can describe it in one of two ways by its wavelength or by its frequency.

Highlights from the Telescope Array Experiment

Results from the Telescope Array Experiment

Transcription:

Instrumentation for sflash B.K. Shin a *, S. Atwood b, K. Belov c, J. Belz b, P. Chen d, C. Field e, M. Fukushima f,g, C. Hast e, J. Huang d, H. Huey d, T. Liu d, D. Ivanov b, K. Jobe e, C. Jui b, J. Nam d, C. Naudet c, J. Matthews b, M. Potts b, K. Reil e, D. Saltzberg h, P. Sokolsky b, S. Thomas b, G. B. Thomson b, S. Wang d, and for the sflash a Faculty of Science, Osaka City University, Osaka, 558-8585, Japan. b Department of Physics and Astronomy, University of Utah, Salt Lake City, UT 84108, USA. c Jet Propulsion Laboratory, Pasadena, CA 91109, USA. d Dept. of Physics, Grad. Inst. of Astrophys.,& Leung Center for Cosmology and Particle Astrophysics, National Taiwan University, Taipei, Taiwan. e SLAC National Accelerator Laboratory, Menlo Park, CA, 94025, USA. f Institute for Cosmic Ray Research, University of Tokyo, Kashiwa, Chiba, Japan. g Kavli Institute for the Physics and Mathematics of the Universe, University of Tokyo, Kashiwa, Chiba, Japan. h Dept. of Physics and Astronomy, Univ. of California, Los Angeles, CA 90095, USA. E-mail: bkshin@sci.osaka-cu.ac.jp We will report on the setup and calibration of the instrumentation for sflash. The sflash experiment is a measurement of the air fluorescence from 10 21 ev artificial air showers developed in an alumina target by an electron beam provided from End Station A of the Stanford Linear Accelerator Center (SLAC). sflash employs Photomultiplier tubes (PMTs) to detect the fluorescence photons and an Integrating Charge Transformer (ICT) to measure the beam current. The PMTs are positioned ~10m perpendicular to shower axis are used to measure the air fluorescence photons. The absolute gain of the PMTs was measured using CRAYS at the Institute of Cosmic-Ray Research, Japan. CRAYS is a PMT photo-gain calibration system consisting of a vessel containing nitrogen gas and a 14µJ nitrogen (337.1 nm) laser. The gain of each PMT was monitored during the course of the experiment using a YAP pulser (YAlO: Ce + 241 Am) potted into an ultraviolet optical filter attached to the surface of the PMT window. We expect to achieve an overall experimental uncertainty of ~6%. The ICT (Bergoz) consists of 20 turns of coil to measure the beam current. In sflash experiment we set it ~1 m away from beam exit to measure the beam current just prior to entering the target material. The ICT was calibrated with 2 ns pulse generator at Jet Propulsion Laboratory. The calibration uncertainty of the ICT is ~1%. 35th International Cosmic Ray Conference ICRC2017 10-20 July, 2017 Bexco, Busan, Korea *Speaker ã Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License (CC BY-NC-ND 4.0). http://pos.sissa.it/

Instruments of sflash 1.Introduction Figure 1. A picture of sflash experiment. The shower was generated by beam interaction with alumina target. The PMT for measurement of AFY photons was set 10 m away from shower axis. We covered the path of AFY by block clothes to prevent light from building. 2. Overview of sflash Configuration The sflash measured the AFY produced by an air shower of particles developed in alumina from an incident electron beam provided by ESA in SLAC. The definition of the AFY is a number of air fluorescence photons per unit energy deposited in the air. We measured a number of photons using photo-multiplier tubes (PMTs) and estimated the energy deposition by Geant4 simulation with the beam energy obtained from SLAC and beam charge as monitored by Integral Charge Transformer (ICT). Figure 1 and 2 show a picture and a schematic of final setup and configuration of sflash during September 2016 in ESA. We installed the PMTs ~10 meters perpendicular from beam line and ICT between the movable table with alumina brick targets and the beam dump. To reject unexpected photons from extraneous sources and gamma rays from targets and beam dump, we covered the area including the PMTs and air shower volume with black cloth to create a dark box and set concrete shields between the PMTs and potential noise sources. A window of width 1000 ± 1 mm in beam direction was set ~94 cm away from shower axis to precisely limit the volume 2 The air fluorescence yield (AFY) is the primary parameter required for the measurement of the energy of Ultra High Energy Cosmic Rays (UHECR) observed using the air fluorescence technique. The Extensive Air Shower (EAS) is a phenomenon in which the interaction between a high energy cosmic rays (>1015 ev) and the atmosphere creates a shower of daughter particles. We can observe EAS via fluorescence photons emitted from nitrogen molecules excited by the passage of charged daughter particles. The AFY is typically defined as the number of air fluorescence photons produced per unit of energy deposited in the atmosphere by an EAS. The spectrum of air fluorescence has roughly 20 emission lines in the near ultraviolet. Approximately 26% of the light is in 337 nm line (2P 0 0) of nitrogen molecules. Previous attempts to measure the AFY have measured the air fluorescence produced by particle beams in air. To better understand the real AFY from atmospheric EAS, we started the sflash experiment to measure air fluorescence produced by actual air showers initiated in alumina by ~1021 ev equivalent electron beam provided from End Station A (ESA) of the Stanford Linear Accelerator Center (SLAC). We succeeded with taking data in September in 2016. We described the configuration and instruments of last sflash experiments in this paper.

of energy deposition viewed by the PMTs. The window was optionally blocked by a remotely controlled shutter to permit the measurement of the background contribution to the air fluorescence signal. The total energy in this experiment was close to maximum energy of real UHECR (~10 21 ev), the so-called GZK limit [1][2]. We develop the air shower by the electron beam interaction with various thicknesses of alumina (AlO3) targets to create 10 3 or 10 4 times more particles or lower energy than the original beam. The target blocks were arranged as shown in Figure 2 (green box) to prevent incident beam electrons from potentially passing through gaps between the target alumina blocks. The alumina blocks are on a remotely controlled moveable table to change the depth of target in discrete steps from 0 to 10 inches. We evaluated the beam background with PMT during June 2016 and succeeded with taking air fluorescence data in September in 2016 with this configuration. The details the instruments what used for the precise measurement of AFY in September 2016 are described in the following sections. Figure 2. Right) Map of sflash configuration, Left top) a picture of PMT setup, sflash used 7PMTs. Left bottom) a picture of beam exit, the target is ~2 m away from beam window and the ICT measured beam charge between beam window and target. 2.1 Photomultiplier Tubes Figure 2 left top is a picture of PMT setup of sflash. We used seven PMTs which are the same models used in the air fluorescence telescopes of the TA experiment [3]. The left four PMTs (R9508, HAMATSU) are used at the Black Rock (BR) station positioned at the south-east of TA array. The other three PMTs (XP30620-FL, Photonics) are used at the Middle Drum (MD) station position at the north of TA site. A 337.1 nm narrow bandpass filter was installed on one of BR PMTs to examine the fraction of air fluorescence observed at this wavelength as compared to the total AFY. Another PMT has 3.14 cm 2 mask to evaluate the systematic uncertainty by comparison of a fully exposed PMTs with one with only the central most sensitive area of the photocathode exposed. All of the MD PMTs used a 2.54 cm diameter of mask. Independent PMT calibrations to get a relation between a number of detected photons and output charge from PMT anode were performed for the two types of PMTs. The absolute gain of the BR PMTs were determined using the CRAYS system at the Institutes of Cosmic Ray Research, Japan. The MD PMTs were independently calibrated using a Roving Xenon Flasher (RXF) at The University of Utah, USA. The RXF calibration was performed primarily as a basic cross check. 3

The most recent absolute calibration of the RXF was performed in November 2003. The 2003 RXF calibration had an uncertainty of ~7%. The primary absolute calibration for sflash is the CRAYS calibration. The calibration constant of GPMT is defined number of detected photons (Np) per output charge from PMT anode (Q PMT ) which shown in below, G PMT = Np/Q PMT. (1) The CRAYS is system of PMT calibration which consists of laser and scattering chamber. A PMT was set on the chamber which is orthogonal to laser shooting direction. And 337.1 nm wavelength, 5 ns laser pulse is injected into a scattering chamber filled with 99.9% purity nitrogen gas. 90 degree scattered photons illuminate the PMT photocathode through a 16 cm 2 window which located at a distance of 4 cm perpendicular to the beam axis [4]. The PMTs have a small light source YAP (YAlO 3 ) potted on the UV filter to monitor any relative change in gain between the CRAYS calibration and the measurements at SLAC. The MD PMTs were calibrated using the Roving Xenon Flasher (RXF). The RXF is used for the absolute calibration at the MD FD and for relative calibration between TA FD telescopes. A PMT and RXF were set ~2.1 meters away each other in dark room. The known output of the RXF was used to obtain a calibration factor for the PMTs. The RXF produces ~11 photons per mm 2 at this distance and near room temperature. Figure 3. On site calibration using UVLED, The calibration was operated when beam off. The UVLED was position at lower than a path of AFY, so it was not interrupt with AFY measurement. For on-site monitoring of the PMT calibration, we used a temperature stabilized UV LED light source. The UV LED was developed to príovide a new type standard candle for cross calibration between telescopes in TA. The ~367 nm light output of the UV LED can be adjusted in frequency, pulse width, and amplitude. The temperature of the UV LED and associated electronics is stabilized to ±0.01 accuracy to provide a high stability output independent of the temperature of environment. Figure 3 shows setup in the sflash experiment. The UVLED was set ~4 meters away from PMTs and ~1.15 m from the ground for avoiding obscuration of air fluorescence photons from the window. The UV LED settings were temperature = 45.00, amplitude = 160 ma (~263.2 pj), and width = 1.00 µsec. The table 1 is a data set measured at sflash and comparison between BR PMT1 and MD PMT5. PMT Name HV Energy [pj] Width [us] Q [nvs] Gain [Np/nVs] λ eff AREA [cm2] Np/Cm2 Ratio BR PMT1 679 263.2 1.0 9.439 2540.40 1.007 31.06 7.77E+02 100.00% MD PMT5 886.6 263.2 1.0 13.670 258.57 5.06 6.85E+02 88.19% Table 1, the preliminary result of UVLED run in sflash. The discrepancy of photo density between two PMT was 12%. 4

The comparison of measured photon-density is about ~ 12% discrepancy. This is consistent with the relative uncertainty of each calibration. We are presently evaluating the UV LED light yield and cross calibration between two types of PMTs in a laboratory of the University of Utah. The Number of detected photon (Np) for AFY measurement by PMT was calculated shown in below formula, Np = (Np open Np close )/ω (2), where NP open and NP close are the numbers of photons calculated from the product of Q PMT and G PMT in the condition of the shutter opened and closed, respectively, to reject background photons. The factor ω is a correction for the wavelength distribution of AFY. As a result of past FLASH experiment, AFY has ~20 wavelengths band spectrum in a range of 290~420 nm. So the ω is calibration factor of GPMT considering of AFY spectrum with the transmittance of UV filter (T Filter ) on the photocathode of PMT and quantum efficiency (QE) of each type of PMT shown in below. ω =Σ R λ T filter_λ /T filter_std QE λ /QE std (3), where R λ is a relative yield fraction for wavelength λ and std means the wavelength of light source for PMT calibration. The resulting set of AF Np measurement of two PMTs is shown in Figure 4. The x and y axis in Figure 4 indicate the target depth and detected photon density on the PMT cathode. The red and black dots are the results for BR and MD PMTs, respectively. We obtain a similar relative photon density ratio between MD and BR PMTs as compared to the UV LED. We will further examine the ~12% uncertainty between two types of PMTs with careful calibration of the UV LED. Figure 4. AFY photon density by two PMTs.The x-axis shows the depth of target and the y-axis shows photon density. The photon density is calculated Np over area of photo chathode. The result shows that a BR PMT is about 10% higher than a MD PMT. We will investgate with UVLED calibration. 2.3 Beam Monitoring We monitored the charge of the beam pulse using an Integrating Charge Transformer (ICT) and a horn. The ICT (In-air ICT, Bergoz) is designed to measure the beam charge with an accuracy 5

of ±0.4 pc in the air. The ICT was positioned ~50 cm from beam exit window between the window and the alumina target. The output of the ICT was read out using an oscilloscope (TDC694C, Tektronix) with a preamplifier. The number of electrons (Ne) was estimated from the measured pulse area from ICT (Q ICT ) multiplied by calibration factor (F ICT ). The FICT was calibrated using pulse generator in a laboratory of JPL. The ICT was calibrated with a pulse generator and an oscilloscope. For calibration we used 2 ns and 160 pc pulse which is the similar to conditions from the electron beam used in sflash. Then we measured output charge of ICT with preamp which was used in sflash and input pulse with an oscilloscope. The result of F ICT is determined 0.08016032 [pc/nvs]. Figure 5 shows number of electrons of a data set. The beam charge was 1.057 10 9 (or ~160 pc) and shoot to shoot fluctuation was 2%. The beam charge was stable while we measure the AFY. Figure 5. the distribution of number of electrons of beam. The beam charge was about 109 electrons and shot to shot fluctuation was 2% which is stable during beam shooting. The Beam spot size was monitored using two screens with web cameras. One screen was set on the beam exit window and the other screen was set on the beam dump shown in Figure 6. Beam spot size was ~ 1 cm in diameter at beam exit and dump which means the beam was well focused and does not coherence with ICT. Figure 6. The beam spots at beam window (left) and beam dump (right). The size of beam spots were ~1 cm at beam window and dump. It means beam was well focused and not interuppt with ICT. 6

2.4 Atmospheric condition The atmospheric conditions temperature, humidity, and pressure are important parameters of AFY [5]. And PMT gain depends on temperature. The temperature and humidity were monitored by thermometer (1922T, KN Lab) with 0.5 accuracy in temperature and 5% accuracy in relative humidity. We set the two thermometers in the AFY dark box and near the PMTs. Figure 7 shows a history temperature and humidity during sflash. The temperature and relative humidity were 23.0 ± 1 and ~40 ± 5%. The local atmospheric pressure was obtained from a SLAC weather station. This weather station is about 100 m away from ESA where sflash is operated. Figure 7. Atomosphereic condition in the ESA during sflash. The air termperature and humidity in building were very stable. The temperature was 23.0 ± 1 and relative humidity was ~40 ± 5 %. 4. Summary The main target of sflash experiment is to measure AFY precisely for the study of UHECRs. We used ~10 21 ev energy artificially generated air showers which is an energy of primary interest for the study of UHECRs. The shower was generated by electron beam provided from SLAC with various depths of alumina targets. The definition of AFY is a number of AF photons per unit energy deposition. The number of photons is determined by two types of PMTs and energy deposition is estimated by Geant4 simulation with the beam energy and charge measured by SLAC and ICT. To detect AF photons from energy deposition of the shower and to measure the charge of the beam, we set the PMTs 10 m way from shower center and ICT between a bunch of targets and beam window. As an uncertainty of preliminary result is higher than 12% which the main contribution from the discrepancy of photo density between two types PMT. However, we expect to solve this with UV LED calibration. After solving this we can achieve more accuracy result. We also report the preliminary result of sflash in 2016 and simulation of energy deposition in this conference. Acknowledgements The Telescope Array experiment and the sflash experiment are supported by the Japan Society for the Promotion of Science through Grants-in-Aids for Scientific Research on Specially Promoted Research (15H05693) and for Scientific Research (S) (15H05741), and the Inter-University Research Program of the Institute for Cosmic Ray Research; by the U.S. National Science Foundation awards PHY-0307098, PHY-0601915, PHY-0649681, PHY- 7

0703893, PHY-0758342, PHY-0848320, PHY-1069280, PHY-1069286, PHY-1404495 and PHY-1404502; by the National Research Foundation of Korea (2015R1A2A1A01006870, 2015R1A2A1A15055344, 2016R1A5A1013277, 2007-0093860, 2016R1A2B4014967); by the Russian Academy of Sciences, RFBR grant 16-02-00962a (INR), IISN project No. 4.4502.13, and Belgian Science Policy under IUAP VII/37 (ULB). The foundations of Dr. Ezekiel R. and Edna Wattis Dumke, Willard L. Eccles, and George S. and Dolores Dor\'e Eccles all helped with generous donations. The State of Utah supported the project through its Economic Development Board, and the University of Utah through the Office of the Vice President for Research. The experimental site became available through the cooperation of the Utah School and Institutional Trust Lands Administration (SITLA), U.S. Bureau of Land Management (BLM), and the U.S. Air Force. We appreciate the assistance of the State of Utah and Fillmore offices of the BLM in crafting the Plan of Development for the site. We also wish to thank the people and the officials of Millard County, Utah for their steadfast and warm support. We gratefully acknowledge the contributions from the technical staffs of our home institutions. An allocation of computer time from the Center for High Performance Computing at the University of Utah is gratefully acknowledged. The authors thank SLAC National Accelerator Center for providing facilities and support. Part of this research was funded through the JPL Internal Research and Technology Development program. References [1] Greisen, Physical Review Letters 16 (17): 748 750. (1966) [2] Zatsepin and Kuzmin, Journal of Experimental and Theoretical Physics Letters 4: 78 80. (1966) [3] H. Kawai et al., Nuclear Physics B (Proc. Suppl.) 175 176 (2008) 221 226 [4] S. ~Kawana et al., Nuclear Instruments and Methods in Physics Research A 681 (2012) 68-77 [5] B. Keilhauer et al., EPJ Web of Conferences 53, 01010 (2013) 8