The Lateral Trigger Probability function for the Ultra-High Energy Cosmic Ray Showers detected by the Pierre Auger Observatory
Pierre Auger Collaboration, P. Abreu, M. Aglietta, E. J. Ahn, I. F. M. Albuquerque, D. Allard, I. Allekotte, J. Allen, P. Allison, J. Alvarez Castillo, J. Alvarez-Muñiz, M. Ambrosio, A. Aminaei, L. Anchordoqui, S. Andringa, T. Antičić, A. Anzalone, C. Aramo, E. Arganda, F. Arqueros, H. Asorey, P. Assis, J. Aublin, M. Ave, M. Avenier, G. Avila, T. Bäcker, M. Balzer, K. B. Barber, A. F. Barbosa, R. Bardenet, S. L. C. Barroso, B. Baughman, J. Bäuml, J. J. Beatty, B. R. Becker, K. H. Becker, A. Bellétoile, J. A. Bellido, S. BenZvi, C. Berat, X. Bertou, P. L. Biermann, P. Billoir, F. Blanco, M. Blanco, C. Bleve, H. Blümer, M. Boháčová, D. Boncioli, C. Bonifazi, R. Bonino, N. Borodai, J. Brack, P. Brogueira, W. C. Brown, R. Bruijn, P. Buchholz, A. Bueno, R. E. Burton, K. S. Caballero-Mora, L. Caramete, R. Caruso, A. Castellina, O. Catalano, G. Cataldi, L. Cazon, R. Cester, J. Chauvin, S. H. Cheng, A. Chiavassa, J. A. Chinellato, A. Chou, J. Chudoba, R. W. Clay, M. R. Coluccia, R. Conceição, F. Contreras, H. Cook, M. J. Cooper, J. Coppens, A. Cordier, S. Coutu, C. E. Covault, A. Creusot, A. Criss, J. Cronin, A. Curutiu, S. Dagoret-Campagne, R. Dallier, S. Dasso, K. Daumiller, B. R. Dawson, R. M. de Almeida, M. De Domenico, C. De Donato, S. J. de Jong, G. De La Vega, W. J. M. de Mello Junior, J. R. T. de Mello Neto, et al. (398 additional authors not shown)
aa r X i v : . [ a s t r o - ph . H E ] N ov The Lateral Trigger Probability function for theUltra-High Energy Cosmic Ray Showersdetected by the Pierre Auger Observatory
The Pierre Auger Collaboration
P. Abreu , M. Aglietta , E.J. Ahn , I.F.M. Albuquerque , D. Allard ,I. Allekotte , J. Allen , P. Allison , J. Alvarez Castillo , J. Alvarez-Mu˜niz ,M. Ambrosio , A. Aminaei , L. Anchordoqui , S. Andringa , T. Antiˇci´c ,A. Anzalone , C. Aramo , E. Arganda , F. Arqueros , H. Asorey , P. Assis ,J. Aublin , M. Ave , M. Avenier , G. Avila , T. B¨acker , M. Balzer ,K.B. Barber , A.F. Barbosa , R. Bardenet , S.L.C. Barroso , B. Baughman ,J. B¨auml , , J.J. Beatty , B.R. Becker , K.H. Becker , A. Bell´etoile ,J.A. Bellido , S. BenZvi , C. Berat , X. Bertou , P.L. Biermann , P. Billoir ,F. Blanco , M. Blanco , C. Bleve , H. Bl¨umer , , M. Boh´aˇcov´a , ,D. Boncioli , C. Bonifazi , , R. Bonino , N. Borodai , J. Brack ,P. Brogueira , W.C. Brown , R. Bruijn , P. Buchholz , A. Bueno ,R.E. Burton , K.S. Caballero-Mora , L. Caramete , R. Caruso , A. Castellina ,O. Catalano , G. Cataldi , L. Cazon , R. Cester , J. Chauvin , S.H. Cheng ,A. Chiavassa , J.A. Chinellato , A. Chou , , J. Chudoba , R.W. Clay ,M.R. Coluccia , R. Concei¸c˜ao , F. Contreras , H. Cook , M.J. Cooper ,J. Coppens , , A. Cordier , U. Cotti , S. Coutu , C.E. Covault ,A. Creusot , , A. Criss , J. Cronin , A. Curutiu , S. Dagoret-Campagne ,R. Dallier , S. Dasso , , K. Daumiller , B.R. Dawson , R.M. de Almeida , ,M. De Domenico , C. De Donato , , S.J. de Jong , G. De La Vega , W.J.M. deMello Junior , J.R.T. de Mello Neto , I. De Mitri , V. de Souza , K.D. deVries , G. Decerprit , L. del Peral , O. Deligny , H. Dembinski , , N. Dhital ,C. Di Giulio , , J.C. Diaz , M.L. D´ıaz Castro , P.N. Diep , C. Dobrigkeit ,W. Docters , J.C. D’Olivo , P.N. Dong , , A. Dorofeev , J.C. dos Anjos ,M.T. Dova , D. D’Urso , I. Dutan , J. Ebr , R. Engel , M. Erdmann ,C.O. Escobar , A. Etchegoyen , P. Facal San Luis , I. Fajardo Tapia ,H. Falcke , , G. Farrar , A.C. Fauth , N. Fazzini , A.P. Ferguson , A. Ferrero ,B. Fick , A. Filevich , A. Filipˇciˇc , , S. Fliescher , C.E. Fracchiolla ,E.D. Fraenkel , U. Fr¨ohlich , B. Fuchs , R. Gaior , R.F. Gamarra ,S. Gambetta , B. Garc´ıa , D. Garc´ıa G´amez , D. Garcia-Pinto , A. Gascon ,H. Gemmeke , K. Gesterling , P.L. Ghia , , U. Giaccari , M. Giller ,H. Glass , M.S. Gold , G. Golup , F. Gomez Albarracin , M. G´omez Berisso ,P. Gon¸calves , D. Gonzalez , J.G. Gonzalez , B. Gookin , D. G´ora , ,A. Gorgi , P. Gouffon , S.R. Gozzini , E. Grashorn , S. Grebe , N. Griffith ,M. Grigat , A.F. Grillo , Y. Guardincerri , F. Guarino , G.P. Guedes ,A. Guzman , J.D. Hague , P. Hansen , D. Harari , S. Harmsma , ,J.L. Harton , A. Haungs , T. Hebbeker , D. Heck , A.E. Herve , C. Hojvat ,N. Hollon , V.C. Holmes , P. Homola , J.R. H¨orandel , A. Horneffer ,M. Hrabovsk´y , , T. Huege , A. Insolia , F. Ionita , A. Italiano , C. Jarne ,S. Jiraskova , K. Kadija , K.H. Kampert , P. Karhan , P. Kasper , B. K´egl ,B. Keilhauer , A. Keivani , J.L. Kelley , E. Kemp , R.M. Kieckhafer ,H.O. Klages , M. Kleifges , J. Kleinfeller , J. Knapp , D.-H. Koang ,K. Kotera , N. Krohm , O. Kr¨omer , D. Kruppke-Hansen , F. Kuehn ,D. Kuempel , J.K. Kulbartz , N. Kunka , G. La Rosa , C. Lachaud , Preprint submitted to Elsevier November 8, 2018 . Lautridou , M.S.A.B. Le˜ao , D. Lebrun , P. Lebrun , M.A. Leigui deOliveira , A. Lemiere , A. Letessier-Selvon , I. Lhenry-Yvon , K. Link ,R. L´opez , A. Lopez Ag¨uera , K. Louedec , J. Lozano Bahilo , A. Lucero , ,M. Ludwig , H. Lyberis , M.C. Maccarone , C. Macolino , S. Maldera ,D. Mandat , P. Mantsch , A.G. Mariazzi , J. Marin , , V. Marin , I.C. Maris ,H.R. Marquez Falcon , G. Marsella , D. Martello , L. Martin , H. Martinez ,O. Mart´ınez Bravo , H.J. Mathes , J. Matthews , , J.A.J. Matthews ,G. Matthiae , D. Maurizio , P.O. Mazur , G. Medina-Tanco , M. Melissas ,D. Melo , , E. Menichetti , A. Menshikov , P. Mertsch , C. Meurer ,S. Mi´canovi´c , M.I. Micheletti , W. Miller , L. Miramonti , S. Mollerach ,M. Monasor , D. Monnier Ragaigne , F. Montanet , B. Morales , C. Morello ,E. Moreno , J.C. Moreno , C. Morris , M. Mostaf´a , C.A. Moura , ,S. Mueller , M.A. Muller , G. M¨uller , M. M¨unchmeyer , R. Mussa ,G. Navarra † , J.L. Navarro , S. Navas , P. Necesal , L. Nellen , A. Nelles ,P.T. Nhung , L. Niemietz , N. Nierstenhoefer , D. Nitz , D. Nosek ,L. Noˇzka , M. Nyklicek , J. Oehlschl¨ager , A. Olinto , P. Oliva ,V.M. Olmos-Gilbaja , M. Ortiz , N. Pacheco , D. Pakk Selmi-Dei ,M. Palatka , J. Pallotta , N. Palmieri , G. Parente , E. Parizot , A. Parra ,R.D. Parsons , S. Pastor , T. Paul , M. Pech , J. P¸ekala , R. Pelayo ,I.M. Pepe , L. Perrone , R. Pesce , E. Petermann , S. Petrera , P. Petrinca ,A. Petrolini , Y. Petrov , J. Petrovic , C. Pfendner , N. Phan , R. Piegaia ,T. Pierog , P. Pieroni , M. Pimenta , V. Pirronello , M. Platino , V.H. Ponce ,M. Pontz , P. Privitera , M. Prouza , E.J. Quel , S. Querchfeld ,J. Rautenberg , O. Ravel , D. Ravignani , B. Revenu , J. Ridky , S. Riggi , ,M. Risse , P. Ristori , H. Rivera , V. Rizi , J. Roberts , C. Robledo ,W. Rodrigues de Carvalho , , G. Rodriguez , J. Rodriguez Martino , ,J. Rodriguez Rojo , I. Rodriguez-Cabo , M.D. Rodr´ıguez-Fr´ıas , G. Ros ,J. Rosado , T. Rossler , M. Roth , B. Rouill´e-d’Orfeuil , E. Roulet ,A.C. Rovero , C. R¨uhle , F. Salamida , , H. Salazar , G. Salina , F. S´anchez ,M. Santander , C.E. Santo , E. Santos , E.M. Santos , F. Sarazin , B. Sarkar ,S. Sarkar , R. Sato , N. Scharf , V. Scherini , H. Schieler , P. Schiffer ,A. Schmidt , F. Schmidt , T. Schmidt , O. Scholten , H. Schoorlemmer ,J. Schovancova , P. Schov´anek , F. Schr¨oder , S. Schulte , D. Schuster ,S.J. Sciutto , M. Scuderi , A. Segreto , M. Settimo , A. Shadkam ,R.C. Shellard , , I. Sidelnik , G. Sigl , H.H. Silva Lopez , A. ´Smia lkowski ,R. ˇSm´ıda , , G.R. Snow , P. Sommers , J. Sorokin , H. Spinka , ,R. Squartini , J. Stapleton , J. Stasielak , M. Stephan , E. Strazzeri ,A. Stutz , F. Suarez , T. Suomij¨arvi , A.D. Supanitsky , , T. ˇSuˇsa ,M.S. Sutherland , , J. Swain , Z. Szadkowski , , M. Szuba , A. Tamashiro ,A. Tapia , M. Tartare , O. Ta¸sc˘au , C.G. Tavera Ruiz , R. Tcaciuc ,D. Tegolo , , N.T. Thao , D. Thomas , J. Tiffenberg , C. Timmermans , ,D.K. Tiwari , W. Tkaczyk , C.J. Todero Peixoto , , B. Tom´e , A. Tonachini ,P. Travnicek , D.B. Tridapalli , G. Tristram , E. Trovato , M. Tueros , ,R. Ulrich , , M. Unger , M. Urban , J.F. Vald´es Galicia , I. Vali˜no , ,L. Valore , A.M. van den Berg , E. Varela , B. Vargas C´ardenas ,J.R. V´azquez , R.A. V´azquez , D. Veberiˇc , , V. Verzi , J. Vicha ,M. Videla , L. Villase˜nor , H. Wahlberg , P. Wahrlich , O. Wainberg ,D. Warner , A.A. Watson , M. Weber , K. Weidenhaupt , A. Weindl ,S. Westerhoff , B.J. Whelan , G. Wieczorek , L. Wiencke , B. Wilczy´nska ,H. Wilczy´nski , M. Will , C. Williams , T. Winchen , L. Winders , .G. Winnick , M. Wommer , B. Wundheiler , T. Yamamoto a , T. Yapici ,P. Younk , G. Yuan , A. Yushkov , , B. Zamorano , E. Zas ,D. Zavrtanik , , M. Zavrtanik , , I. Zaw , A. Zepeda , M. Ziolkowski Centro At´omico Bariloche and Instituto Balseiro (CNEA- UNCuyo-CONICET),San Carlos de Bariloche, Argentina Centro At´omico Constituyentes (Comisi´on Nacional de Energ´ıaAt´omica/CONICET/UTN-FRBA), Buenos Aires, Argentina Centro de Investigaciones en L´aseres y Aplicaciones, CITEFA and CONICET,Argentina Departamento de F´ısica, FCEyN, Universidad de Buenos Aires y CONICET,Argentina IFLP, Universidad Nacional de La Plata and CONICET, La Plata, Argentina Instituto de Astronom´ıa y F´ısica del Espacio (CONICET- UBA), Buenos Aires,Argentina Instituto de F´ısica de Rosario (IFIR) - CONICET/U.N.R. and Facultad deCiencias Bioqu´ımicas y Farmac´euticas U.N.R., Rosario, Argentina National Technological University, Faculty Mendoza (CONICET/CNEA),Mendoza, Argentina Pierre Auger Southern Observatory, Malarg¨ue, Argentina Pierre Auger Southern Observatory and Comisi´on Nacional de Energ´ıa At´omica,Malarg¨ue, Argentina University of Adelaide, Adelaide, S.A., Australia Centro Brasileiro de Pesquisas Fisicas, Rio de Janeiro, RJ, Brazil Pontif´ıcia Universidade Cat´olica, Rio de Janeiro, RJ, Brazil Universidade de S˜ao Paulo, Instituto de F´ısica, S˜ao Carlos, SP, Brazil Universidade de S˜ao Paulo, Instituto de F´ısica, S˜ao Paulo, SP, Brazil Universidade Estadual de Campinas, IFGW, Campinas, SP, Brazil Universidade Estadual de Feira de Santana, Brazil Universidade Estadual do Sudoeste da Bahia, Vitoria da Conquista, BA, Brazil Universidade Federal da Bahia, Salvador, BA, Brazil Universidade Federal do ABC, Santo Andr´e, SP, Brazil Universidade Federal do Rio de Janeiro, Instituto de F´ısica, Rio de Janeiro, RJ,Brazil Universidade Federal Fluminense, EEIMVR, Volta Redonda, RJ, Brazil Rudjer Boˇskovi´c Institute, 10000 Zagreb, Croatia Charles University, Faculty of Mathematics and Physics, Institute of Particle andNuclear Physics, Prague, Czech Republic Institute of Physics of the Academy of Sciences of the Czech Republic, Prague,Czech Republic Palacky University, RCATM, Olomouc, Czech Republic Institut de Physique Nucl´eaire d’Orsay (IPNO), Universit´e Paris 11,CNRS-IN2P3, Orsay, France Laboratoire AstroParticule et Cosmologie (APC), Universit´e Paris 7,CNRS-IN2P3, Paris, France Laboratoire de l’Acc´el´erateur Lin´eaire (LAL), Universit´e Paris 11, CNRS-IN2P3,Orsay, France Laboratoire de Physique Nucl´eaire et de Hautes Energies (LPNHE), Universit´esParis 6 et Paris 7, CNRS-IN2P3, Paris, France Laboratoire de Physique Subatomique et de Cosmologie (LPSC), Universit´eJoseph Fourier, INPG, CNRS-IN2P3, Grenoble, France SUBATECH, CNRS-IN2P3, Nantes, France Bergische Universit¨at Wuppertal, Wuppertal, Germany Karlsruhe Institute of Technology - Campus North - Institut f¨ur Kernphysik,Karlsruhe, Germany Karlsruhe Institute of Technology - Campus North - Institut f¨urProzessdatenverarbeitung und Elektronik, Karlsruhe, Germany Karlsruhe Institute of Technology - Campus South - Institut f¨ur ExperimentelleKernphysik (IEKP), Karlsruhe, Germany Max-Planck-Institut f¨ur Radioastronomie, Bonn, Germany RWTH Aachen University, III. Physikalisches Institut A, Aachen, Germany Universit¨at Hamburg, Hamburg, Germany Universit¨at Siegen, Siegen, Germany Dipartimento di Fisica dell’Universit`a and INFN, Genova, Italy Universit`a dell’Aquila and INFN, L’Aquila, Italy Universit`a di Milano and Sezione INFN, Milan, Italy Dipartimento di Fisica dell’Universit`a del Salento and Sezione INFN, Lecce, Italy Universit`a di Napoli ”Federico II” and Sezione INFN, Napoli, Italy Universit`a di Roma II ”Tor Vergata” and Sezione INFN, Roma, Italy Universit`a di Catania and Sezione INFN, Catania, Italy Universit`a di Torino and Sezione INFN, Torino, Italy Dipartimento di Ingegneria dell’Innovazione dell’Universit`a del Salento andSezione INFN, Lecce, Italy Istituto di Astrofisica Spaziale e Fisica Cosmica di Palermo (INAF), Palermo, Italy Istituto di Fisica dello Spazio Interplanetario (INAF), Universit`a di Torino andSezione INFN, Torino, Italy INFN, Laboratori Nazionali del Gran Sasso, Assergi (L’Aquila), Italy Universit`a di Palermo and Sezione INFN, Catania, Italy Benem´erita Universidad Aut´onoma de Puebla, Puebla, Mexico Centro de Investigaci´on y de Estudios Avanzados del IPN (CINVESTAV), M´exico,D.F., Mexico Universidad Michoacana de San Nicolas de Hidalgo, Morelia, Michoacan, Mexico Universidad Nacional Autonoma de Mexico, Mexico, D.F., Mexico IMAPP, Radboud University, Nijmegen, Netherlands Kernfysisch Versneller Instituut, University of Groningen, Groningen, Netherlands NIKHEF, Amsterdam, Netherlands ASTRON, Dwingeloo, Netherlands Institute of Nuclear Physics PAN, Krakow, Poland University of L´od´z, L´od´z, Poland LIP and Instituto Superior T´ecnico, Lisboa, Portugal J. Stefan Institute, Ljubljana, Slovenia Laboratory for Astroparticle Physics, University of Nova Gorica, Slovenia Instituto de F´ısica Corpuscular, CSIC-Universitat de Val`encia, Valencia, Spain Universidad Complutense de Madrid, Madrid, Spain Universidad de Alcal´a, Alcal´a de Henares (Madrid), Spain Universidad de Granada & C.A.F.P.E., Granada, Spain Universidad de Santiago de Compostela, Spain Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford,United Kingdom School of Physics and Astronomy, University of Leeds, United Kingdom Argonne National Laboratory, Argonne, IL, USA Case Western Reserve University, Cleveland, OH, USA Colorado School of Mines, Golden, CO, USA Colorado State University, Fort Collins, CO, USA Colorado State University, Pueblo, CO, USA Fermilab, Batavia, IL, USA Louisiana State University, Baton Rouge, LA, USA Michigan Technological University, Houghton, MI, USA New York University, New York, NY, USA Northeastern University, Boston, MA, USA Ohio State University, Columbus, OH, USA Pennsylvania State University, University Park, PA, USA
Southern University, Baton Rouge, LA, USA
University of Chicago, Enrico Fermi Institute, Chicago, IL, USA
University of Nebraska, Lincoln, NE, USA
University of New Mexico, Albuquerque, NM, USA
University of Wisconsin, Madison, WI, USA
University of Wisconsin, Milwaukee, WI, USA
Institute for Nuclear Science and Technology (INST), Hanoi, Vietnam( † ) Deceased(a) at Konan University, Kobe, Japan Abstract
In this paper we introduce the concept of Lateral Trigger Probability (LTP)function, i.e., the probability for an extensive air shower (EAS) to trigger an in-dividual detector of a ground based array as a function of distance to the showeraxis, taking into account energy, mass and direction of the primary cosmic ray.We apply this concept to the surface array of the Pierre Auger Observatoryconsisting of a 1.5 km spaced grid of about 1600 water Cherenkov stations. Us-ing Monte Carlo simulations of ultra-high energy showers the LTP functions arederived for energies in the range between 10 and 10 eV and zenith anglesup to 65 ◦ . A parametrization combining a step function with an exponentialis found to reproduce them very well in the considered range of energies andzenith angles. The LTP functions can also be obtained from data using eventssimultaneously observed by the fluorescence and the surface detector of thePierre Auger Observatory (hybrid events). We validate the Monte-Carlo resultsshowing how LTP functions from data are in good agreement with simulations. Keywords:
Ultra-High Energy Cosmic Rays, Pierre Auger Observatory,Extensive Air Showers, Trigger performance, Surface detector, Hybriddetector. 5 . Introduction
The Pierre Auger Observatory has been conceived to study the origin andthe nature of ultra high-energy cosmic rays. Because of the scarcity of the fluxat the highest energies, their direct measurement from space is technically un-feasible and the use of very large detectors is required at the ground. What canbe observed is the extensive air shower of secondary particles produced in thepropagation through the atmosphere. The Pierre Auger Observatory is locatednear Malarg¨ue, Argentina, at 1400 m a.s.l. and it employes two independent andcomplementary measurement techniques [1]. The surface array (SD), consistingof about 1600 water Cherenkov detectors on a triangular grid of 1.5 km spacingcovering an area of approximately 3000 km , records the secondary particles atthe ground and thus samples their lateral density distribution. The fluorescencedetector (FD), consisting of 24 telescopes at four sites, overlooks the surfacearray and observes the longitudinal profile of air showers by collecting the flu-orescence light emitted along the path through the atmosphere [2]. Unlike thesurface detector array with its nearly 100% duty cycle, the FD can only operateon clear and moonless nights giving an overall duty cycle of about 13% [3]. Asa consequence, only a fraction of showers are observed by both detectors. Forthese events, called hereafter “hybrid”, the combination of information fromthe surface array and the fluorescence telescopes enhances the reconstructioncapability. Energy and direction reconstruction accuracy of hybrid events is infact better than the one the SD and FD could achieve independently.One of the main goals of the Pierre Auger Observatory is to measure the fluxof cosmic rays at the highest energies. This task relies on an accurate determi-nation of the detector exposure for SD-only [4] and hybrid [3] operation modes.The hybrid exposure is calculated using the simultaneous simulation of FD andSD response. Besides the dependence on energy and distance to an FD-site,the hybrid exposure is influenced by several factors including the atmosphericconditions, the trigger status of all active detectors and their instantaneous datataking configuration. The calculation of the SD response is based on the deepknowledge of the array capability to trigger once a shower with a given energyand zenith angle hits the ground. Since the trigger in an EAS array is alwaysa combination of trigger states of neighboring detectors, the acceptance of anyEAS array is directly connected to the probability that an individual detectortriggers when a shower lands at a certain distance from it. This defines the con-cept of “Lateral Trigger Probability” function. This function has been used as apowerful tool for simulations in the analysis for the measurement of the hybridenergy spectrum [5] and of the atmospheric depth at shower maximum [6].The trigger design of the Auger surface detector is described in detail in [4].Each water Cherenkov detector of the surface array has a 10 m water surfacearea and 1.2 m water depth, with three 9 in. photomultiplier tubes (PMTs)looking through optical coupling material into the water volume, which is con-tained in a Tyvek R (cid:13) reflective liner. The signals provided by each PMT aredigitised by 40 MHz 10-bit Flash Analog to Digital Converters (FADCs) [1].The achieved dynamical range is sufficient to cover with good precision both the6ignals produced in the detectors near to the shower axis ( ∼ µs )and those produced far from the shower axis ( ∼ µs ). We recall herethe basic structure of the used trigger algorithms. The two first levels (T1 andT2) are formed at each surface detector. Each trigger level can be divided intwo modes, a threshold trigger (TH) and a time-over-threshold trigger (ToT).The first level threshold trigger (TH-T1) requires the coincidence of the signalsfrom the three PMTs equipping each station, each PMT signal being above1.75 “Vertical Equivalent Muon” (VEM) . The TH-T1 trigger is used to reducethe rate due to atmospheric muons to ∼
100 Hz and can reach the second level,TH-T2, when the peak signal reaches at least 3.2 VEM in coincidence between3 PMTs signals, further reducing the rate to ∼
20 Hz. The second mode, theToT, requires at least 13 time bins (i.e. more than 325 ns) in 120 FADC binsof a sliding window of 3 µ s to be above a threshold of 0.2 VEM in coincidencein 2 out of 3 PMTs. Time-over-threshold trigger stations are automaticallypromoted to the second level. The threshold trigger is especially efficient atdetecting strong narrow signals, mostly encountered in horizontal showers orclose to the axis of vertical showers. On the other hand, the ToT is intended toselect sequences of small signals spread in time. This is typical of low energyvertical showers dominated by an electromagnetic component or of high energyshowers triggering stations at large distance from the shower axis because ofmuons produced high in the atmosphere.Higher level triggers are obtained by requiring the spatial and temporalcoincidence of at least three stations satisfying the T2 conditions. In particular,for zenith angles below 60 ◦ , the full efficiency for SD is reached at 10 . eV [4].In addition, if at least one FD telescope triggers in coincidence with one secondlevel trigger station, a hybrid trigger is formed. Since FD has a lower energythreshold, hybrid events are also detected below the minimum energy for anindependent SD trigger. For zenith angles below 60 ◦ , the hybrid detector reachesnearly full efficiency at 10 eV [3].In section 2, the concept of a Lateral Trigger Probability (LTP) functionis formalized and applied to the particular case of the surface detector of thePierre Auger Observatory. In section 3, the LTP functions for a single time-over-threshold trigger station are derived and parametrized for different primaryparticles (proton, iron, photon) and their dependence on energy and zenith angleis explored for zenith angles up to 65 ◦ and for energies between 10 and 10 eV.This energy range is relevant as it covers the interval in which the SD-only andthe hybrid detection mode become fully efficient. The dependence on the choiceof the hadronic interaction models is also discussed in section 3. In section4, hybrid data are finally used to validate the simulation and to estimate theimpact of weather conditions on the observed efficiency. The LTP functions are The distribution of measured light due to atmospheric muons produces a peak in thePMT charge distribution, Q peakVEM (or VEM in short), as well as a peak in that of the pulseheight, I peakVEM , both of them being proportional to those produced by a vertical through-goingmuon [1].
2. Concept of Lateral Trigger Probability
The trigger probability of a single water Cherenkov detector depends onseveral independent physical parameters : i) the characteristics of the primarycosmic ray that initiates an air shower, e.g., its energy and mass, ii) the typeand geometry of the detector used to observe air showers (in the following wewill only study water Cherenkov detectors used for the surface detector array ofthe Pierre Auger Observatory), iii) the trigger condition used to detect a signalfrom air showers iv) the geometry of the incoming shower, e.g. its incidencezenith angle and position with respect to the detector. To formalize these de-pendencies we define the Lateral Trigger Probability function Λ
A,E,θ,T r ( r, φ ) asthe probability to trigger on an air shower induced by a primary particle of en-ergy E , mass A and zenith angle θ [7]. Here, r and φ are the radial coordinatesof the single detector in the plane normal to the shower axis (shower frame).Using a trigger condition T r , this probability is simply given by:Λ
A,E,θ,T r ( r, φ ) = N N + N . (1)where N and N are respectively the number of triggered and un-triggereddetectors with coordinates r and φ in the shower frame.
3. Simulations
The LTP functions have been derived using detailed simulations of the EASdevelopment and of the detector response. The simulation sample consists ofabout 15000 CORSIKA [8] showers (proton, iron and photon primaries) withzenith angle distributed as sin θ cos θ ( θ < ◦ ) and energies ranging between10 and 10 eV in steps of 0.25 in the logarithmic scale. A “thin sam-pling” mechanism at the level of 10 − (optimal thinning) is applied followingthe standard method used for CORSIKA simulation with energies larger than10 eV [9]. The showers have been generated with the models QGSJETII [10]and FLUKA [11] for high and low energy hadronic interactions.In the simulation, the position of the shower core (i.e. the intersection of theshower axis with the ground) is uniformly distributed over the surface array andeach shower is used 5 times, each time with a different core position, in orderto cover different areas of the array and explore all the detector configurations.The surface detector response is simulated using GEANT4 [12] and adopting thesampling procedure to regenerate particles in a ground detector from thinnedair shower simulations as described in [13]. The entire detector simulation iscarried out within the framework provided by the Auger Offline software [14].The trigger status of SD stations is inspected within a radius of 3 km fromthe shower axis and the Lateral Trigger Probability is then derived according to8q. 1. At distances larger than 3 km, the trigger efficiency is negligibly small forthe class of events studied in this paper. All trigger modes of the surface detec-tor are simulated in detail at all levels. However, for events with zenith anglesbelow 65 ◦ , the majority of the stations forming a second level trigger satisfythe ToT condition. In particular, for the considered zenith angles, the fractionof TH-T2 trigger stations not being also ToT is about 1%, approximately inde-pendent of the energy. Thus, we focus the analysis on the ToT stations.The lateral trigger probability for a ToT station is shown in Fig. 1 at a givenenergy and for different ranges of the cosine of the zenith angle θ . The max-imum effective distance for detection increases with energy and, for a givenenergy, with the cosine of the zenith angle, i.e. events with larger zenith angletend to trigger less due to the attenuation of their electromagnetic component.For moderately inclined showers, an asymmetry is expected in the signal de-tected in the stations placed at the same distance to the shower axis but withdifferent azimuth in the shower frame [15]. Indeed, secondary particles arrivingearlier traverse less atmosphere and are less attenuated than the late ones. As aconsequence, early stations may exhibit larger trigger probabilities and producelarger signals. Actually, for zenith angles below 65 ◦ , this effect has been foundto have a quite low influence on the trigger probability, only noticeable above30 ◦ (in simulations as well as in the data). In the following we consider LTPfunctions averaged over all azimuths in the showers frame. A more detailedtreatment including the azimuthal dependence does not introduce measurabledifferences for acceptance calculations.A fit combining a step function (close to the axis) with an exponential (fur-ther away) reproduces reasonably well the full simulated data set. The form ofthe fit function used is:LTP(r) = e − r − R R r ≤ R e C · ( r − R ) r > R (2)where R , ∆ R and C are free fit parameters, with R being the distance whereLTP is equal to 0.5. A fit performed according to eq. 2 is superimposed on eachplot shown in Fig. 1. As an example, the ToT trigger probability at energy E =10 eV and for two angular bins (vertical showers on the left and showerswith larger zenith angle on the right) is shown in Fig. 2: the exponential canreproduce very well the tail of the probability distribution at large distancesfrom axis, in particular for inclined events. The dependences of fit parameters R , ∆ R and C on energy and zenith angle can be parametrized by quadraticpolynomials in the variables cos θ and log ( E/eV ). The corresponding coeffi-cients are tabulated in the Appendix for proton, iron and photon primaries. InFig. 3, the ToT trigger probability from parametrization has been superimposedon the simulation (proton primary, all zenith angles up to 65 ◦ are merged). Thecomparison is performed as in the following. For each simulated event, i.e. fora certain primary, energy and arrival direction, the LTP is calculated using theparametrization (lines) and shown together with the full simulation (points).9 tation distance to shower axis (km) LT P f o r a ToT s t a t i on (E/eV)=17.5 log <0.6 θ θ θ θ θ station distance to shower axis (km) LT P f o r a ToT s t a t i on (E/eV)=18 log <0.6 θ θ θ θ θ station distance to shower axis (km) LT P f o r a ToT s t a t i on (E/eV)=18.5 log <0.6 θ θ θ θ θ station distance to shower axis (km) LT P f o r a ToT s t a t i on (E/eV)=19 log <0.6 θ θ θ θ θ Figure 1: Lateral Trigger Probability from simulations (proton primary) for aToT station at a given energy, from 10 eV up to 10 eV in steps of 0.5 inthe logarithmic scale. Different bins of cos θ are also shown together with a fitperformed according to eq. 2, superimposed as a continuous line. station distance to shower axis (km) LT P f o r a ToT s t a t i on -2 -1 (E/eV) = 19 log < 1 θ station distance to shower axis (km) LT P f o r a ToT s t a t i on -2 -1 (E/eV) = 19 log < 0.6 θ Figure 2: Fit made with a step function in proximity of the shower axis (con-tinuous line) and by an exponential at larger distances (dashed line). The ToTprobability is shown for vertical (left) and inclined (right) showers at energy of10 eV. 10he agreement is remarkably good in the entire energy range for proton (shownin the figure) and for iron and photon primaries. The detector response to showers induced by different primary particles isshown in Fig. 4, for two classes of events, vertical (0 ◦ < θ < ◦ ) on the topand moderately inclined (38 ◦ < θ < ◦ ) on the bottom. Because of theirlarger number of muons, showers induced by iron nuclei provide a higher triggercapability at larger distances than those induced by protons, for all zenith angles.However, the difference between proton and iron is too small to give any hint formass composition analysis. On the other hand, the LTP functions for photonprimaries differ sensibly from those of hadrons (they vanish at shorter distances,about 500 m less at an energy of 10 eV). This is a consequence of the structureof the lateral distribution of photon showers, i.e. at a given energy, their effective footprint at the ground is smaller than the one of hadrons. Moreover, in photonshowers there is a much smaller number of muons.It is worth noting that the energy threshold corresponding to full efficiencyfor SD, derived from data and simulation in ref. [4], has been found to becompatible with the expectation for hadronic primaries. Different choices of high energy interaction models influence the simulation ofshower development and could affect the expected trigger efficiency. The depen-dence of the Lateral Trigger Probability on the assumptions for the hadronicinteraction model has been investigated using a sample of simulated showers(proton) produced with SIBYLL [16]. As shown in Fig. 5, the LTP functionsderived with the two hadronic interaction models differ only at large distancefrom the shower axis, in a range where the efficiency degrades rapidly. In thisregion, SIBYLL gives a lower LTP since this model predicts on average a smallernumber of muons. Those differences are however too small to imply an observ-able impact on the detector acceptance.
4. LTP functions from data and comparison with simulation
The LTP functions can be derived from data by calculating the ratio of trig-gered to active stations within a given distance from the reconstructed showeraxis. While doing this, the actual surface detector configuration must be accu-rately taken into account as a function of time. In addition, only high qualitydata are selected to avoid biases due to mis-reconstructed energies and/or ar-rival directions. The use of hybrid events allows to derive LTP functions alsofor energies below the threshold of an independent SD trigger. This is a benefitof the hybrid design that aims to fully exploit the distinctive potential offeredby the Pierre Auger Observatory. Two years of hybrid data collected betweenJune 2006 and May 2008 were used for this study. The events are selected asdescribed in [3] and this ensures an angular resolution of about 0.6 ◦ and a core11 tation distance to shower axis (km)0.5 1 1.5 2 2.5 3 LT P f o r a ToT s t a t i on (E/eV)=17 log (E/eV)=17.25 log (E/eV)=17.5 log (E/eV)=17.75 log (E/eV)=18 log (E/eV)=18.25 log (E/eV)=18.5 log (E/eV)=18.75 log (E/eV)=19 log ° < zenith < 65 ° Figure 3: Lateral Trigger Probability for a ToT station as a function of stationdistance to shower axis and for different energies (proton primary). The outcomeof the parametrization is superimposed as a line. All zenith angles up to 65 ◦ are merged. 12 tation distance to shower axis (km)0.5 1 1.5 2 2.5 3 LT P f o r a ToT s t a t i on iron proton photon ° < zenith < 38 ° (E/eV)=19 log (E/eV)=19 log (E/eV)=19 log (E/eV)=19 log (E/eV)=19 log station distance to shower axis (km)0.5 1 1.5 2 2.5 3 LT P f o r a ToT s t a t i on (E/eV)=19 log iron proton photon ° < zenith < 65 ° (E/eV)=19 log Figure 4: Lateral Trigger Probability for a ToT station. Proton, iron andphoton primaries of energy 10 eV for two zenith angle ranges, 0 ◦ - 38 ◦ (top)and 38 ◦ - 65 ◦ (bottom). The outcome of the parametrization is superimposedas a continuous line. station distance to shower axis (km)0.5 1 1.5 2 2.5 3 LT P f o r a ToT s t a t i on proton QGSJETII proton SIBYLL station distance to shower axis (km) QG S J E T II/ S I BY LL (E/eV)=19 log Figure 5: Lateral Trigger Probability for a ToT station (zenith angle between 0 ◦ and 65 ◦ ). Proton primary at energy of 10 eV with QGSJETII and SIBYLL.The ratio QGSJETII/SIBYLL is shown in the inset.13osition determination better than 70 m. Further requirements on the goodnessof the reconstructed longitudinal profile provide an energy resolution of about10% above 10 eV and less than 15% at lower energies [3].The LTP measured from data is shown in Fig. 6 for different energy intervals.To verify the performance of the parametrization described in section 3, foreach selected event, the LTP of any active station within 3 km from the showeraxis is calculated using the reconstructed energy and direction. The predictedprobability (dashed line) is then superimposed on data (points), see Fig. 6. Inthis way, data are compared to simulation taking into account the actual statusof the detector. The shaded area gives the interval of expected values assumingthat data are pure proton (lower edge) or pure iron (upper edge). A 50% protonand 50% iron mixed composition has been assumed for the parametrization(dashed line). The agreement is good over the entire energy range. This featureactually starts at very low energies, even below the range of full efficiency forthe hybrid detection [3]. In this case, whereas in data only events with at leastone SD ToT station are selected, in simulation also the events that did nottrigger at all are taken into account in the calculation of the probability. Asa consequence, the comparison between data and simulation could be biased.However, the good level of agreement actually reached reflects the fact that thehybrid detection is very close to fully efficient and the energy reconstructionremains reasonably good within the scope of this analysis down to energy ofabout 10 . eV.For each energy interval considered, the agreement between data and simu-lation has also proven to hold in two zenith angle bands (0 ◦ - 38 ◦ and 38 ◦ - 65 ◦ ),see Fig. 7. The effect of atmospheric variations (in pressure, temperature and air den-sity) on extensive air showers development has been extensively studied withthe surface detector data [17]. A significant modulation of the rate of eventswith the atmospheric variables, both on a seasonal scale ( ∼ ∼
2% on average during a day) has been observed. Thismodulation is mainly explained as due to the change with the air density of theMoli`ere radius near ground thus influencing the trigger probability and the rateof events above a fixed energy. Hybrid data in the energy range around 10 eVhave been used to investigate this effect on LTP. Data have been separated byseason and are shown, together with the parametrization, for austral winter andaustral summer, see Fig. 8, top panel. The ratio of summer and winter rela-tive to the parametrization is shown in the bottom panel. Results qualitativelymatch the expectation. Higher temperature at the ground, as for the australsummer, induces a reduction of the air density weakly enhancing the triggerprobability at a given distance relative to all other seasons. Nevertheless theeffect is almost negligible on the scale of the measurable trigger efficiency.14 tation distance to shower axis (km) LT P f o r a ToT s t a t i on (E/eV)<17.7 station distance to shower axis (km) LT P f o r a ToT s t a t i on (E/eV)<18.2 station distance to shower axis (km) LT P f o r a ToT s t a t i on (E/eV)<18.7 station distance to shower axis (km) LT P f o r a ToT s t a t i on (E/eV)<19.2 Figure 6: Comparison of simulation with hybrid data collected in two years. Allzenith angles up to 65 ◦ merged. The energy intervals are 10 . < E < . eV,10 . < E < . eV, 10 . < E < . eV, 10 . < E < . eV.
5. Summary and Conclusions
In the previous sections we have introduced the concept of Lateral TriggerProbability function as a tool to characterize the single detector trigger effi-ciency. We have derived LTP functions for the particular case of the surfacedetector of the Pierre Auger Observatory using simulations. We discussed theirevolution with different physical parameters of air showers such as the energy,zenith angle and nature of the primary particle. We also investigated the impactof choosing different hadronic interaction models in the simulations. Further-more, we estimated the LTP functions at different energies and zenith anglesusing hybrid data and showed that seasonal effects are visible in the triggerprobabilities retrieved from data as expected from previous studies [17].The good agreement between simulations and data over a wide energy range(between 10 . eV and 10 eV) demonstrates the accuracy of the differentaspects of the simulation procedure (i.e. air shower, detectors and trigger sim-ulation) as well as the quality of the reconstruction obtained for hybrid data.These comparisons support and validate the use of simulated LTP functionsin the estimate of the hybrid aperture described in [3]. Monitoring the LTPfunctions over a longer period of time can be used to study the long-term per-15 tation distance to shower axis (km) LT P f o r a ToT s t a t i on Hybrid data (2006-2008)Parametrization (E/eV) < 17.7 ° < zenith < 38 ° station distance to shower axis (km) LT P f o r a ToT s t a t i on ° < zenith < 65 ° station distance to shower axis (km) LT P f o r a ToT s t a t i on Hybrid data (2006-2008)Parametrization (E/eV) < 18.2 ° < zenith < 38 ° station distance to shower axis (km) LT P f o r a ToT s t a t i on ° < zenith < 65 ° station distance to shower axis (km) LT P f o r a ToT s t a t i on Hybrid data (2006-2008)Parametrization (E/eV) < 18.7 ° < zenith < 38 ° station distance to shower axis (km) LT P f o r a ToT s t a t i on ° < zenith < 65 ° station distance to shower axis (km) LT P f o r a ToT s t a t i on Hybrid data (2006-2008)Parametrization (E/eV) < 19.2 ° < zenith < 38 ° station distance to shower axis (km) LT P f o r a ToT s t a t i on ° < zenith < 65 ° Figure 7: Comparison of simulation with hybrid data collected in two years.Zenith angles are split in two ranges 0 ◦ - 38 ◦ (left) and 38 ◦ - 65 ◦ (right).From top to bottom the energy intervals are 10 . < E < . eV,10 . < E < . eV, 10 . < E < . eV, 10 . < E < . eV.16 tation distance to shower axis (km)0.5 1 1.5 2 2.5 3 LT P f o r a ToT s t a t i on austral summeraustral winterparam (E/eV)<18.2 station distance to shower axis (km)0.5 1 1.5 2 2.5 3 s u mm e r ( w i n t e r ) / p a r a m austral summeraustral winterparam Figure 8: LTP functions from hybrid data at energy of about 10 eV for australwinter and austral summer compared to the parametrization derived in section 3(top) and ratio relative to the parametrization (bottom).17ormance of the SD trigger for individual stations both above and below theacceptance saturation energy.As a final consideration, LTP functions can be derived at higher energiesusing SD-only data because, at energy above ∼ . eV, despite the statis-tics of hybrids becoming small, the surface detector is fully efficient and thegeometrical reconstruction is accurate. As mentioned in the Introduction, theprobability of a high level trigger for the surface detector is a combination ofsingle detector probabilities. Hence LTP functions provide a robust and simplemethod to estimate the energy or zenith angle dependence of SD acceptance forany arbitrary configuration. This makes this technique a valuable tool to designother experiments and future enhancements of the Pierre Auger Observatory.
6. Acknowledgements
The successful installation and commissioning of the Pierre Auger Observa-tory would not have been possible without the strong commitment and effortfrom the technical and administrative staff in Malarg¨ue.We are very grateful to the following agencies and organizations for finan-cial support: Comisi´on Nacional de Energ´ıa At´omica, Fundaci´on Antorchas,Gobierno De La Provincia de Mendoza, Municipalidad de Malarg¨ue, NDMHoldings and Valle Las Le˜nas, in gratitude for their continuing cooperationover land access, Argentina; the Australian Research Council; Conselho Na-cional de Desenvolvimento Cient´ıfico e Tecnol´ogico (CNPq), Financiadora deEstudos e Projetos (FINEP), Funda¸c˜ao de Amparo `a Pesquisa do Estado deRio de Janeiro (FAPERJ), Funda¸c˜ao de Amparo `a Pesquisa do Estado de S˜aoPaulo (FAPESP), Minist´erio de Ciˆencia e Tecnologia (MCT), Brazil; AVCR,AV0Z10100502 and AV0Z10100522, GAAV KJB300100801 and KJB100100904,MSMT-CR LA08016, LC527, 1M06002, and MSM0021620859, Czech Republic;Centre de Calcul IN2P3/CNRS, Centre National de la Recherche Scientifique(CNRS), Conseil R´egional Ile-de-France, D´epartement Physique Nucl´eaire etCorpusculaire (PNC-IN2P3/CNRS), D´epartement Sciences de l’Univers (SDU-INSU/CNRS), France; Bundesministerium f¨ur Bildung und Forschung (BMBF),Deutsche Forschungsgemeinschaft (DFG), Finanzministerium Baden-W¨urttemberg,Helmholtz-Gemeinschaft Deutscher Forschungszentren (HGF), Ministerium f¨urInnovation, Wissenschaft und Forschung, Nordrhein-Westfalen, Ministerium f¨urWissenschaft, Forschung und Kunst, Baden-W¨urttemberg, Germany; IstitutoNazionale di Fisica Nucleare (INFN), Istituto Nazionale di Astrofisica (INAF),Ministero dell’Istruzione, dell’Universit`a e della Ricerca (MIUR), Gran SassoCenter for Astroparticle Physics (CFA), Italy; Consejo Nacional de Ciencia yTecnolog´ıa (CONACYT), Mexico; Ministerie van Onderwijs, Cultuur en Weten-schap, Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO), Sticht-ing voor Fundamenteel Onderzoek der Materie (FOM), Netherlands; Ministryof Science and Higher Education, Grant Nos. 1 P03 D 014 30 and N N202207238, Poland; Funda¸c˜ao para a Ciˆencia e a Tecnologia, Portugal; Ministryfor Higher Education, Science, and Technology, Slovenian Research Agency,Slovenia; Comunidad de Madrid, Consejer´ıa de Educaci´on de la Comunidad18e Castilla La Mancha, FEDER funds, Ministerio de Ciencia e Innovaci´onand Consolider-Ingenio 2010 (CPAN), Generalitat Valenciana, Junta de An-daluc´ıa, Xunta de Galicia, Spain; Science and Technology Facilities Council,United Kingdom; Department of Energy, Contract Nos. DE-AC02-07CH11359,DE-FR02-04ER41300, National Science Foundation, Grant No. 0969400, TheGrainger Foundation USA; NAFOSTED, Vietnam; ALFA-EC / HELEN, Eu-ropean Union 6th Framework Program, Grant No. MEIF-CT-2005-025057, Eu-ropean Union 7th Framework Program, Grant No. PIEF-GA-2008-220240, andUNESCO.
References [1] J. Abraham et al. [Pierre Auger Collaboration], Nuclear Instruments andMethods in Physics Research A523 (2004) 50.[2] J. Abraham et al. [Pierre Auger Collaboration], Nuclear Instruments andMethods in Physics Research A620 (2010) 227.[3] P.Abreu et al. [Pierre Auger Collaboration], Astroparticle Physics 34 (2011)368.[4] J. Abraham et al. [Pierre Auger Collaboration], Nuclear Instruments andMethods in Physics Research A613 (2010) 29.[5] J. Abraham et al. [Pierre Auger Collaboration], Physics Letters B685 (2010)239.[6] J. Abraham et al. [Pierre Auger Collaboration], Physical Review Letters104 (2010) 091101.[7] D. Allard, 29th Int. Cosmic Ray Conf. (2005) arXiv:astro-ph/0511104v1.[8] D. Heck et al. , “CORSIKA: A Monte Carlo Code to Simulate ExtensiveAir Showers” FLUKA: a multi-particle transport code , CERN-2005-10(2005) INFN/TC 05/11, SLAC-R-773; A. Fass`o et al.
The physics mod-els of FLUKA: status and recent developments , Computing in High En-ergy and Nuclear Physics 2003 Conference (CHEP2003), La Jolla, CA,USA, March 24-28, 2003, (paper MOMT005) eConf C0303241 (2003)arXiv:hep-ph/0306267. 1912] S. Agostinelli et al., Nuclear Instruments and Methods in Physics ResearchA506 (2003) 250; IEEE Transactions on Nuclear Science 53 No. 1 (2006)270.[13] P. Billoir, Astroparticle Physics 30 (2008) 270.[14] S. Argir`o et al., Nuclear Instruments and Methods in Physics ResearchA580 (2007) 1485.[15] M.T. Dova, M.E. Manceido, A.G. Mariazzi, H. Wahlberg, F. Arqueros D.Garca-Pinto Astroparticle Physics 31 (2009) 312.[16] Eun-Joo Ahn et al., Physical Review D80, (2009) 094003; R.S. Fletcher etal., Physical Review D50 (1994) 5710.[17] J. Abraham et al. [Pierre Auger Collaboration], Astroparticle Physics 32(2009) 89. 20 ppendix: LTP parametrization
The LTP is fitted, as discussed in section 3, to the following function:LTP(r) = e − r − R R r ≤ R e C · ( r − R ) r > R (A.1)with R being the distance where the LTP is equal to 0.5.The dependences of fit parameters R , ∆ R and C on energy and zenith an-gle can be parametrized by quadratic polynomials in the variables cos θ and log ( E/eV ). The corresponding coefficients are given for proton, iron andphoton primaries (0 ◦ < θ < ◦ ), separately. Concerning the accuracy of theparameters, a change at the level of (1 ÷ Proton showers
The overall parametrization for proton primaries (0 ◦ < θ < ◦ ) is summa-rized in the following matrix equation: R km = θ cos θ T · . · − . · . · − − . · . · − . · − − . · − . · − . · − · ( E/eV )log ( E/eV ) ∆ R km = θ cos θ T · − . · . · − − . · − . · − . · . · − − . · . · − − . · − · ( E/eV )log ( E/eV ) C km − = θ cos θ T · − . · . · − . · − − . · . · − . · − · ( E/eV )log ( E/eV ) Iron showers
The overall parametrization for iron primaries (0 ◦ < θ < ◦ ) is summarizedin the following matrix equation: R km = θ cos θ T · . · − . · . · − − . · . · − . · − − . · . · − . · − · ( E/eV )log ( E/eV ) ∆ R km = θ cos θ T · − . · − . · − . · − . · − − . · − . · − · ( E/eV )log ( E/eV ) C km − = θ cos θ T · − . · . · − . · . · − . · . · · ( E/eV )log ( E/eV ) hoton showers The overall parametrization for photon primaries (0 ◦ < θ < ◦ ) is summa-rized in the following matrix equation: R km = θ cos θ T · . · − . · . · − − . · . · − . · − . · − . · . · − · ( E/eV )log ( E/eV ) ∆ R km = θ cos θ T · . · − . · . · − − . · . · − . · − . · − . · . · − · ( E/eV )log ( E/eV ) C km − = θ cos θ T · − . · . · − . · . · − . · . · − . · . · − . · · ( E/eV )log ( E/eV )