We present the results of the most comprehensive survey of the Galactic plane in very highenergy VHE rays including a public release of Galactic sky maps a catalog of VHE sources and the discovery of 16 new sources of VHE rays. The High Energy Spectroscopic System H.E.S.S. Galactic plane survey HGPS was a decadelong observation program carried out by the H.E.S.S. I array of Cherenkov telescopes in Namibia from 2004 to 2013. The observations amount to nearly 2700 h of qualityselected data covering the Galactic plane at longitudes from   250 to 65 and latitudes b 3. In addition to the unprecedented spatial coverage the HGPS also features a relatively high angular resolution 0.08  5 arcmin mean point spread function 68 containment radius sensitivity 1.5 Crab flux for pointlike sources and energy range 0.2100 TeV. We constructed a catalog of VHE ray sources from the HGPS data set with a systematic procedure for both source detection and characterization of morphology and spectrum. We present this likelihoodbased method in detail including the introduction of a model component to account for unresolved largescale emission along the Galactic plane. In total the resulting HGPS catalog contains 78 VHE sources of which 14 are not reanalyzed here for example due to their complex morphology namely shelllike sources and the Galactic center region. Where possible we provide a firm identification of the VHE source or plausible associations with sources in other astronomical catalogs. We also studied the characteristics of the VHE sources with source parameter distributions. 16 new sources were previously unknown or unpublished and we individually discuss their identifications or possible associations. We firmly identified 31 sources as pulsar wind nebulae PWNe supernova remnants SNRs composite SNRs or gammaray binaries. Among the 47 sources not yet identified most of them 36 have possible associations with cataloged objects notably PWNe and energetic pulsars that could power VHE PWNe. P The source catalog is available at the CDS via anonymous ftp to A hrefhttpcdsarc.ustrasbg.frhttpcdsarc.ustrasbg.frA A hrefhttpcdsarc.ustrasbg.frhttp130.79.128.5A or via A hrefhttpcdsarc.ustrasbg.frvizbinqcatJAA612A1httpcdsarc.ustrasbg.frvizbinqcatJAA612A1A
['gamma rays: general', 'surveys', 'Galaxy: general', 'Astrophysics - High Energy Astrophysical Phenomena']
The H.E.S.S. Galactic plane survey
{'The H.E.S.S. Galactic plane survey': "H.E.S.S. Collaboration, H. Abdalla 1 , A. Abramowski 2 , F. Aharonian 3 ; 4 ; 5 , F. Ait Benkhali 3 , E.O. Angüner 21 , M. Arakawa 43 , M. Arrieta 15 , P. Aubert 24 , M. Backes 8 , A. Balzer 9 , M. Barnard 1 , Y. Becherini 10 , J. Becker Tjus 11 , D. Berge 12 , S. Bernhard 13 , K. Bernlöhr 3 , R. Blackwell 14 , M. Böttcher 1 , C. Boisson 15 , J. Bolmont 16 , S. Bonnefoy 37 , P. Bordas 3 , J. Bregeon 17 , F. Brun GLYPH<3> 26 , P. Brun 18 , M. Bryan 9 , M. Büchele 36 , T. Bulik 19 , M. Capasso 29 , S. Carrigan 3 ; 48 , S. Caro GLYPH<11> 30 , A. Carosi 24 , S. Casanova 21 ; 3 , M. Cerruti 16 , N. Chakraborty 3 , R.C.G. Chaves GLYPH<3> 17 ; 22 , A. Chen 23 , J. Chevalier 24 , S. Colafrancesco 23 , B. Condon 26 , J. Conrad 27 ; 28 , I.D. Davids 8 , J. Decock 18 , C. Deil GLYPH<3> 3 , J. Devin 17 , P. deWilt 14 , L. Dirson 2 , A. Djannati-Ataï 31 , W. Domainko 3 , A. Donath GLYPH<3> 3 , L.O'C. Drury 4 , K. Dutson 33 , J. Dyks 34 , T. Edwards 3 , K. Egberts 35 , P. Eger 3 , G. Emery 16 , J.-P. Ernenwein 20 , S. Eschbach 36 , C. Farnier 27 ; 10 , S. Fegan 30 , M.V. Fernandes 2 , A. Fiasson 24 , G. Fontaine 30 , A. Förster 3 , S. Funk 36 , M. Füßling 37 , S. Gabici 31 , Y.A. Gallant 17 , T. Garrigoux 1 , H. Gast 3 ; 49 , F. Gaté 24 , G. Giavitto 37 , B. Giebels 30 , D. Glawion 25 , J.F. Glicenstein 18 , D. Gottschall 29 , M.-H. Grondin 26 , J. Hahn 3 , M. Haupt 37 , J. Hawkes 14 , G. Heinzelmann 2 , G. Henri 32 , G. Hermann 3 , J.A. Hinton 3 , W. Hofmann 3 , C. Hoischen 35 , T. L. Holch 7 , M. Holler 13 , D. Horns 2 , A. Ivascenko 1 , H. Iwasaki 43 , A. Jacholkowska 16 , M. Jamrozy 38 , D. Jankowsky 36 , F. Jankowsky 25 , M. Jingo 23 , L. Jouvin 31 , I. Jung-Richardt 36 , M.A. Kastendieck 2 , K. Katarzy'nski 39 , M. Katsuragawa 44 , U. Katz 36 , D. Kerszberg 16 , D. Khangulyan 43 , B. Khélifi 31 , J. King 3 , S. Klepser 37 , D. Klochkov 29 , W. Klu'zniak 34 , Nu. Komin 23 , K. Kosack 18 , S. Krakau 11 , M. Kraus 36 , P.P. Krüger 1 , H. La GLYPH<11> on 26 , G. Lamanna 24 , J. Lau 14 , J.-P. Lees 24 , J. Lefaucheur 15 , A. Lemière 31 , M. Lemoine-Goumard 26 , J.-P. Lenain 16 , E. Leser 35 , T. Lohse 7 , M. Lorentz 18 , R. Liu 3 , R. López-Coto 3 , I. Lypova 37 , V. Marandon GLYPH<3> 3 , D. Malyshev 29 , A. Marcowith 17 , C. Mariaud 30 , R. Marx 3 , G. Maurin 24 , N. Maxted 14 ; 45 , M. Mayer 7 , P.J. Meintjes 40 , M. Meyer 27 , A.M.W. Mitchell 3 , R. Moderski 34 , M. Mohamed 25 , L. Mohrmann 36 , K. Morå 27 , E. Moulin 18 , T. Murach 37 , S. Nakashima 44 , M. de Naurois 30 , H. Ndiyavala 1 , F. Niederwanger 13 , J. Niemiec 21 , L. Oakes 7 , P. O'Brien 33 , H. Odaka 44 , S. Ohm 37 , M. Ostrowski 38 , I. Oya 37 , M. Padovani 17 , M. Panter 3 , R.D. Parsons 3 , M. Paz Arribas 7 , N.W. Pekeur 1 , G. Pelletier 32 , C. Perennes 16 , P.-O. Petrucci 32 , B. Peyaud 18 , Q. Piel 24 , S. Pita 31 , V. Poireau 24 , H. Poon 3 , D. Prokhorov 10 , H. Prokoph 12 , G. Pühlhofer 29 , M. Punch 31 ; 10 , A. Quirrenbach 25 , S. Raab 36 , R. Rauth 13 , A. Reimer 13 , O. Reimer 13 , M. Renaud 17 , R. de los Reyes 3 , F. Rieger 3 ; 41 , L. Rinchiuso 18 , C. Romoli 4 , G. Rowell 14 , B. Rudak 34 , C.B. Rulten 15 , S. Safi-Harb 50 , V. Sahakian 6 ; 5 , S. Saito 43 , D.A. Sanchez 24 , A. Santangelo 29 , M. Sasaki 36 , M. Schandri 36 , R. Schlickeiser 11 , F. Schüssler 18 , A. Schulz 37 , U. Schwanke 7 , S. Schwemmer 25 , M. Seglar-Arroyo 18 , M. Settimo 16 , A.S. Sey GLYPH<11> ert 1 , N. Shafi 23 , I. Shilon 36 , K. Shiningayamwe 8 , R. Simoni 9 , H. Sol 15 , F. Spanier 1 , M. Spir-Jacob 31 , Ł. Stawarz 38 , R. Steenkamp 8 , C. Stegmann 35 ; 37 , C. Steppa 35 , I. Sushch 1 , T. Takahashi 44 , J.-P. Tavernet 16 , T. Tavernier 31 , A.M. Taylor 37 , R. Terrier 31 , L. Tibaldo 3 , D. Tiziani 36 , M. Tluczykont 2 , C. Trichard 20 , M. Tsirou 17 , N. Tsuji 43 , R. Tu GLYPH<11> s 3 , Y. Uchiyama 43 , D.J. van der Walt 1 , C. van Eldik 36 , C. van Rensburg 1 , B. van Soelen 40 , G. Vasileiadis 17 , J. Veh 36 , C. Venter 1 , A. Viana 3 ; 46 , P. Vincent 16 , J. Vink 9 , F. Voisin 14 , H.J. Völk 3 , T. Vuillaume 24 , Z. Wadiasingh 1 , S.J. Wagner 25 , P. Wagner 7 , R.M. Wagner 27 , R. White 3 , A. Wierzcholska 21 , P. Willmann 36 , A. Wörnlein 36 , D. Wouters 18 , R. Yang 3 , D. Zaborov 30 , M. Zacharias 1 , R. Zanin 3 , A.A. Zdziarski 34 , A. Zech 15 , F. Zefi 30 , A. Ziegler 36 , J. Zorn 3 , and N. ˙ Zywucka 38  \n(A GLYPH<14> liations can be found after the references)  \nApril 10, 2018", 'ABSTRACT': "We present the results of the most comprehensive survey of the Galactic plane in very high-energy (VHE) GLYPH<13> -rays, including a public release of Galactic sky maps, a catalog of VHE sources, and the discovery of 16 new sources of VHE GLYPH<13> -rays. The High Energy Spectroscopic System (H.E.S.S.) Galactic plane survey (HGPS) was a decade-long observation program carried out by the H.E.S.S. I array of Cherenkov telescopes in Namibia from 2004 to 2013. The observations amount to nearly 2700 h of quality-selected data, covering the Galactic plane at longitudes from ' = 250 GLYPH<14> to 65 GLYPH<14> and latitudes j b j 6 3 GLYPH<14> . In addition to the unprecedented spatial coverage, the HGPS also features a relatively high angular resolution (0 : 08 GLYPH<14> GLYPH<25> 5 arcmin mean point spread function 68% containment radius), sensitivity ( . 1 : 5%Crab flux for point-like sources), and energy range (0.2 to 100 TeV). We constructed a catalog of VHE GLYPH<13> -ray sources from the HGPS data set with a systematic procedure for both source detection and characterization of morphology and spectrum. We present this likelihood-based method in detail, including the introduction of a model component to account for unresolved, large-scale emission along the Galactic plane. In total, the resulting HGPS catalog contains 78 VHE sources, of which 14 are not reanalyzed here, for example, due to their complex morphology, namely shell-like sources and the Galactic center region. Where possible, we provide a firm identification of the VHE source or plausible associations with sources in other astronomical catalogs. We also studied the characteristics of the VHE sources with source parameter distributions. The 16 new sources were previously unknown or unpublished, and we individually discuss their identifications or possible associations. We firmly identified 31 sources as pulsar wind nebulae (PWNe), supernova remnants (SNRs), composite SNRs, or gamma-ray binaries. Among the 47 sources not yet identified, most of them (36) have possible associations with cataloged objects, notably PWNe and energetic pulsars that could power VHE PWNe.  \nKey words. H.E.S.S., Galactic plane survey, Cherenkov telescopes, gamma-ray sources", '1. Introduction': "In this paper, we present the results from the High Energy Spectroscopic System (H.E.S.S.) Galactic plane survey (HGPS), the deepest and most comprehensive survey of the inner Milky Way Galaxy undertaken so far in very high-energy (VHE; 0 : 1 < GLYPH<24> E < GLYPH<24> 100 TeV) GLYPH<13> -rays. Results include numerous sky images (maps) and a new source catalog that is the successor of two previous HGPS releases. The first release (Aharonian et al. 2005b) was based on GLYPH<24> 140 h of observations with the imaging atmospheric Cherenkov telescope (IACT) array H.E.S.S. and contained eight previously unknown sources of VHE GLYPH<13> -rays. In the second release (Aharonian et al. 2006d), we used 230 h of data, covering ' = 330 GLYPH<14> to 30 GLYPH<14> in Galactic longitude and j b j GLYPH<20> 3 GLYPH<14> in latitude. In total, we detected 22 sources of GLYPH<13> -rays in that data set. Since then, the HGPS data set enlarged by more than one order of magnitude in observation time, now comprising roughly 2700 h of high-quality data recorded in the years 2004 - 2013. The spatial coverage is also significantly larger, now encompassing the region from ' = 250 GLYPH<14> to 65 GLYPH<14> in longitude. H.E.S.S. provided periodic updates on this progress by publishing new unidentified sources (Aharonian et al. 2008a) and through conference proceedings (Hoppe 2008b; Chaves et al. 2008a; Chaves 2009; Gast et al. 2011; Deil 2012; Carrigan et al. 2013a,b).  \nCompared to the first HGPS releases over a decade ago, the deeper exposure over a much larger sky area of the Galaxy, combined with improved GLYPH<13> -ray reconstruction, analysis, and modeling techniques, now results in a new catalog containing 78 VHE GLYPH<13> -ray sources. Figure 1 illustrates the HGPS region and compares this region to the structure of the Galaxy, represented by an all-sky Planck CO(1-0) map, and the smaller regions of previous surveys performed by the IACT arrays HEGRA (High-EnergyGamma-Ray Astronomy, Aharonian et al. 2002) and VERITAS (Very Energetic Radiation Imaging Telescope Array System, Weinstein 2009). Even though the HGPS covers only a few percent of the entire sky, this region contains the vast majority of the known Galactic Fermi -LAT 2FHL GLYPH<13> -ray sources (Ackermann et al. 2016) 1 . The figure also shows the measured integral VHE GLYPH<13> -ray flux and the HGPS observation times. As can be seen from the map of observation times (Fig. 1, lower panel), the HGPS data set is not homogeneous. Nonetheless, the HGPS features on average a point-source sensitivity better than 1.5% Crab 2 in the core survey region within 60 GLYPH<14> in longitude of the Galactic center (see Fig. 4, lower panel).  \n2 Throughout this paper, and as is generally the case in VHE GLYPH<13> -ray astronomy, we use the Crab Nebula flux as a standard candle reference: 1 Crab unit is defined here as GLYPH<8> ( > 1TeV) = 2 : 26 GLYPH<1> 10 GLYPH<0> 11 cm GLYPH<0> 2 s GLYPH<0> 1 (Aharonian et al. 2006b).  \nIn this paper, we aim to present the entire data set of the HGPS in a way that is accessible and useful for the whole astronomical community. We have made the maps of VHE GLYPH<13> -ray significance, flux, upper limits, and sensitivity available online 3 for the first time in FITS format (Pence et al. 2010). We developed a semi-automatic analysis pipeline to construct a catalog by detecting and modeling discrete sources of VHE GLYPH<13> -ray emission present in these survey maps. We applied a standardized methodology to the characterization of the GLYPH<13> -ray sources to measure their morphological and spectral properties. The goal was to perform a robust analysis of sources in the survey region with as little manual intervention as possible. With such a generic approach, the catalog pipeline is not optimal for the few very bright and extended sources with complex (non-Gaussian) morphology. For these sources, dedicated analyses are more appropriate, and in all cases, they have already been performed and published elsewhere. We therefore exclude these sources, which are listed in Table 1 below, from the pipeline analysis but include the results from the dedicated analysis in the HGPS catalog for completeness.  \nWe have structured the present paper as follows: we describe the H.E.S.S. telescope array, the data set, and the analysis techniques in Sect. 2. We provide the maps of the VHE GLYPH<13> -ray sky in various representations and details of their production in Sect. 3. Section 4 explains how the HGPS catalog of GLYPH<13> -ray sources was constructed, then Sect. 5 presents and discusses the results, including source associations and identifications with other astronomical objects. Section 6 concludes the main paper with a summary of the HGPS and its results. In Sect. 7, we describe the supplementary online material (maps and catalog in electronic form), including caveats concerning measurements derived from the maps and catalog.", '2.1. The High Energy Stereoscopic System (H.E.S.S.)': "H.E.S.S. is an array of five IACTs located at an altitude of 1800 m above sea level in the Khomas highland of Namibia. It detects Cherenkov light emitted by charged particles in an electromagnetic extensive air shower (EAS) initiated when a primary photon ( GLYPH<13> -ray) of su GLYPH<14> cient energy enters Earth's atmosphere. This array consists of four smaller telescopes, built and operated in the first phase of the experiment (H.E.S.S. Phase I ) and a fifth much larger telescope, which was added to the center of the array in 2012 to launch the second phase (H.E.S.S. Phase II ) of the experiment.  \nH.E.S.S. accumulated the data presented here exclusively with the H.E.S.S. array during its first phase. These four H.E.S.S. Phase I telescopes have tessellated mirrors with a total area of 107 m 2 and cameras consisting of 960 photomultipli-  \nFig. 1. Illustration of HGPS region superimposed an all-sky image of Planck CO(1-0) data (Planck Collaboration et al. 2016) in Galactic coordinates and Hammer-Aito GLYPH<11> projection. For comparison, we overlay the HEGRA Galactic plane survey (Aharonian et al. 2002) and VERITAS Cygnus survey (Weinstein 2009) footprints. Triangles denote the Fermi -LAT 2FHL GLYPH<13> -ray sources (Ackermann et al. 2016) identified as Galactic, and stars indicate the 15 Galactic VHE GLYPH<13> -ray sources outside the HGPS region. H.E.S.S. has detected three of these, which are labeled SN 1006 (Acero et al. 2010a), the Crab Nebula (Aharonian et al. 2006a; H.E.S.S. Collaboration et al. 2014b), and HESS J0632 + 057 (Aharonian et al. 2007; Aliu et al. 2014a). The gray shaded regions denote the part of the sky that cannot be observed from the H.E.S.S. site at reasonable zenith angles (less than 60 GLYPH<14> ). The lower panels show the HGPS GLYPH<13> -ray flux above 1 TeV for regions where the sensitivity is better than 10% Crab (correlation radius R c = 0 : 4 GLYPH<14> ; see Sect. 3) and observation time, both also in Galactic coordinates. The white contours in the lower panels delineate the boundaries of the survey region; the HGPS has little or no exposure beyond Galactic latitudes of j b j GLYPH<20> 3 GLYPH<14> at most locations along the Galactic plane.  \n<!-- image -->  \ners. The energy threshold of the four-telescope array is roughly 200 GeV at zenith and increases with increasing zenith angle. We can reconstruct the arrival direction and energy of the primary photon with accuracies of GLYPH<24> 0 : 08 GLYPH<14> and GLYPH<24> 15%, respectively. Because of its comparatively large field of view (FoV), 5 GLYPH<14> in diameter, the H.E.S.S. Phase I array is well suited for survey operations. The relative acceptance for GLYPH<13> -rays is roughly uniform for the innermost 2 GLYPH<14> of the FoV and gradually drops toward the edges to 40% of the peak value at 4 GLYPH<14> diameter (Aharonian et al. 2006b).", '2.2. Observations, quality selection, and survey region': "The HGPS data set covers the period from January 2004 to January 2013. H.E.S.S. acquired this data set by pointing the IACT array to a given position in the sky for a nominal duration of 28 min (referred to as an observation run hereafter). We considered all runs with zenith angles up to 65 GLYPH<14> and observation positions centered in the Galactic coordinate range ' = 244 : 5 GLYPH<14> to 77 : 5 GLYPH<14> and j b j < 7 : 0 GLYPH<14> . To reduce systematic e GLYPH<11> ects arising from imperfect instrument or atmospheric conditions, we carefully selected good-quality runs as close as possible to the nominal description of the instrument used in the Monte Carlo (MC) sim-  \nulations (see Aharonian et al. 2006b). For example, the IACT cameras su GLYPH<11> er from occasional hardware problems a GLYPH<11> ecting individual or groups of camera pixels, so we did not use observation runs with significant pixel problems. In addition, we only used those runs with at least three operational telescopes.  \nFurthermore, despite the very good weather conditions at the H.E.S.S. site, both nightly and seasonal variations of the atmospheric transparency occur and require monitoring. Layers of dust or haze in the atmosphere e GLYPH<11> ectively act as a filter of the Cherenkov light created in an EAS, thereby raising the energy threshold for triggering the IACTs. Since we calculated the instrument response tables describing the performance of the instrument (e.g., the e GLYPH<11> ective areas) with MC simulations, deviations from the atmospheric conditions assumed in the simulations lead to systematic uncertainties in the determination of energy thresholds, reconstructed energies, and GLYPH<13> -ray fluxes. To account for this, we applied a further quality cut using only observations where the Cherenkov transparency coe GLYPH<14> cient T (Hahn et al. 2014), which characterizes the atmospheric conditions, falls within the range 0 : 8 < T < 1 : 2 (for clear skies, T = 1).  \nAfter applying the aforementioned data quality selection cuts, 6239 observation runs remain, GLYPH<24> 77% of which are runs with four telescopes operational. The total observation time is 2864 h, corresponding to a total livetime of 2673 h (6.7% average dead time). The third panel of Fig. 1 is a map of the observation time over the survey region, clearly showing a non-uniform exposure. This is a result of the HGPS observation strategy, summarized as follows:  \n- -Dedicated survey observations, taken with a typical spacing between pointings of 0 : 7 GLYPH<14> in longitude and in di GLYPH<11> erent latitude bands located between b = GLYPH<0> 1 : 8 GLYPH<14> and b = 1 GLYPH<14> . In addition, for the longitude bands ' = 355 GLYPH<14> to 5 GLYPH<14> and ' = 38 GLYPH<14> to 48 GLYPH<14> , we extended the survey observations in latitude, adding observation pointings from b = GLYPH<0> 3 : 5 GLYPH<14> to b = 3 : 5 GLYPH<14> to explore the possibility of high-latitude emission.\n- -Deeper follow-up observations of source candidates ('hot spots') seen in previous survey observations.\n- -Exploratory and follow-up observations of astrophysical objects located inside the survey region that were promising candidates for emitting VHE GLYPH<13> -rays.\n- -Observations to extend the HGPS spatial coverage and fillup observations to achieve a more uniform sensitivity across the Galactic plane.  \nCombining all of these observations, we achieved a more uniform, minimum 2% Crab flux sensitivity in the region between ' = 283 GLYPH<14> to 58 GLYPH<14> and b = GLYPH<0> 0 : 3 GLYPH<14> GLYPH<6> 0 : 7 GLYPH<14> (see the sensitivity map in Fig. 4).  \nArticle number, page 4 of 79", '2.3. Event reconstruction and selection': 'We first converted the camera pixel data to amplitudes measured in units of photoelectrons (p.e.), identifying the non-operational pixels for a given observation following the procedures described by Aharonian et al. (2004a). We then applied standard H.E.S.S. techniques for the analysis of the camera images: image cleaning, Hillas moment analysis, and the stereo reconstruction of the direction of the primary photon, described by Aharonian et al. (2006b). To suppress the hadronic background and select photon candidate events, we used a multivariate machine learning technique using boosted decision trees based on EAS and image shape parameters (Ohm et al. 2009). For the generation of the survey maps (Sect. 3), we applied the hard cuts configuration whereas for the extraction of source spectra (Sect. 5) we used the standard cuts. The most important distinguishing cut is a minimum of 160 p.e. for hard cuts and 60 p.e. for standard cuts, but there are other di GLYPH<11> erences. See Ohm et al. (2009) for further information; specifically, we used the GLYPH<16> analysis cuts listed in Table 2(a) for the HGPS.  \nWe cross-checked the results presented in this paper with an alternative calibration, reconstruction, and gamma-hadron separation method based on a semi-analytical description of the EAS development (de Naurois & Rolland 2009) with hard cuts of 120 p.e. for maps and standard cuts of 60 p.e. for spectra.  \nFor the energy reconstruction of the primary photons, we compared the image amplitudes in the cameras to the mean amplitudes found in MC simulations of the array (Bernlöhr 2008). Those simulations, which were analyzed with the same chain as the real data for the sake of consistency, include the detailed optical and electronic response of the instrument. The range of optical e GLYPH<14> ciencies encountered in the HGPS data set is large; e GLYPH<14> ciencies start at 100% of the nominal value and drop to almost 50% for some telescopes prior to the mirror refurbishments conducted in 2009-2011. Therefore, we produced several sets of MC simulations, each with optical e GLYPH<14> ciencies of the four telescopes corresponding to their states at suitably chosen times: at the start of H.E.S.S. operations; at the point when e GLYPH<14> ciencies had dropped to GLYPH<24> 70%, before the first mirror refurbishment campaign; and after the mirror refurbishment of each telescope. We then chose the set of simulations most closely matching the state of the system at a given time. Finally, we corrected the remaining di GLYPH<11> erence between simulated and actual optical e GLYPH<14> ciencies using a calibration technique based on the intensity of ring-shaped images from individual muons producing Cherenkov radiation above a telescope (Bolz 2004; Leroy 2004).', '3. HGPS sky maps': "In this section, we describe the methods used to produce the HGPS sky maps. We used the sky maps as the basis for subsequent construction of the HGPS source catalog; this catalog is  \nalso a data product that we release to the community along with this work.  \nWe first computed sky maps for each individual observation run. We then summed these maps over all observations. We chose to use a Cartesian projection in Galactic coordinates, covering the region from ' = 70 GLYPH<14> to 250 GLYPH<14> and b = GLYPH<6> 5 GLYPH<14> , and we set the pixel size to 0 : 02 GLYPH<14> / pixel.  \nIn Sect. 3.1, we describe the production of the map containing the detected events (events map). In Sect. 3.2, we describe the map of expected background events (acceptance map, Sect. 3.2.1), the estimation of a refined background map by introducing exclusion regions (Sect. 3.2.2), and the usage of the adaptive ring background method (Sect. 3.2.3). We then continue in Sect. 3.3 by describing the computation of the significance map, and, in Sect. 3.4, the exposure map (Sect. 3.4.1), which is used to derive quantities such as flux (Sect. 3.4.2), flux error and upper limits (Sect. 3.4.3), and sensitivities (Sect. 3.4.4).", '3.1. Events map': 'The events map consists of the reconstructed positions of the primary GLYPH<13> -ray photons from all events in the sky. To avoid systematic e GLYPH<11> ects near the edge of the FoV in each observation run, we only include events for which the direction of the primary photon is reconstructed within 2 GLYPH<14> of the center of the FoV. This choice results in an e GLYPH<11> ective analysis FoV of 4 GLYPH<14> diameter.  \nAt the lowest energies, the energy reconstruction is biased by EASs with upward fluctuations in the amount of detected Cherenkov light; downward fluctuations do not trigger the cameras. In order to derive reliable flux maps (see Sect. 3.4.2), we only kept events with an energy reconstructed above a defined safe energy threshold. We chose the level of this safe threshold such that, for each run, the energy bias as determined by MC simulations is below 10% across the entire FoV. This conservative approach (together with the use of hard analysis cuts defined in Sect. 2.3) leads to energy threshold values ranging from GLYPH<24> 400 GeV, where the array observed close to zenith, up to 2 TeV at 65 GLYPH<14> from zenith. Figure 2 plots the variation of the safe energy threshold with Galactic longitude, showing the energy threshold for each observation together with the minimum value for each longitude. The variations observed are mainly due to the zenith angle dependency, and regions of di GLYPH<11> erent Galactic longitude generally are observable at di GLYPH<11> erent zenith angles.', '3.2. Background estimation': 'Events passing the event reconstruction and selection procedure are considered GLYPH<13> -ray candidate events. Since these events are still dominantly from EASs induced by GLYPH<13> -ray-like cosmic rays and electrons or positrons, we estimated the amount of remaining background events on a statistical basis using a ring model (Berge et al. 2007) as detailed further below. For each test  \nposition, we counted the photon candidates found in a suitable ring-shaped region around that position in the same FoV. This yields an estimate of the background level after proper normalization and after excluding regions with actual GLYPH<13> -ray emission from the background estimate.', '3.2.1. Acceptance map': 'The acceptance map represents the number of expected events from cosmic-ray backgrounds estimated from runs of sky regions at similar zenith angles but without VHE GLYPH<13> -ray sources. As for the events map (see Sect. 3.1), we computed the acceptance map for energies above the safe energy threshold. To account for the di GLYPH<11> erences in optical e GLYPH<14> ciency and observation time between these runs and those under analysis, we normalized the acceptance map such that, outside the exclusion regions (see Sect. 3.2.2), the number of expected counts matches the number of measured counts. The acceptance maps are used to derive the normalization coe GLYPH<14> cient between the region of interest and the background region (see Sect. 3.3).', '3.2.2. Exclusion regions': 'The background estimation method described above only works if regions with VHE GLYPH<13> -ray emission are excluded from the background estimation region. We defined exclusion regions automatically using an iterative algorithm to avoid potential observer bias and to treat the entire data set in a uniform way. The procedure starts with the significance maps (see Sect. 3.3) produced for the two standard correlation radii R c = 0 : 1 GLYPH<14> and 0 : 2 GLYPH<14> . These radii define the circular region over which a quantity (e.g., GLYPH<13> -ray excess) is integrated. The procedure identifies regions above 5 GLYPH<27> and expands them by excluding an additional 0 : 3 GLYPH<14> beyond the 5 GLYPH<27> contour. This procedure is conservative; it minimizes the amount of surrounding signal that could potentially contaminate the background estimation. A first estimation of the exclusion regions is then included in the significance map production and a new set of exclusion regions is derived. We iterated this procedure until stable regions are obtained, which typically occurs after three iterations. The resulting regions are shown in Fig. 39 below.', '3.2.3. Adaptive ring method': 'In the HGPS, often exclusion regions cover a significant fraction of the FoV; therefore, we could not use the standard ring background method (Berge et al. 2007). For example, using a typical outer ring radius of GLYPH<24> 0.8 GLYPH<14> would lead to numerous holes in the sky maps at positions where the entire ring would be contained inside an exclusion region (i.e., where no background estimation was possible). A much larger outer radius (e.g., GLYPH<24> 1.5 GLYPH<14> ) would be necessary to prevent these holes but would lead to unnecessarily  \nFig. 2. HGPS minimum safe energy threshold as a function of Galactic longitude for a latitude of b = 0 GLYPH<14> . The blue curve shows the minimum threshold for hard cuts (used for maps), and the green curve indicates standard cuts (used for spectra). The black dots represent the safe threshold for each observation run obtained for the hard cuts configuration. The few black dots below the blue line correspond to runs at Galactic latitude j b j > 2 GLYPH<14> .  \n<!-- image -->  \nFig. 3. Illustration of the adaptive ring method for background estimation for a single observation (see Sect. 3.2.3). The HGPS significance image is shown in inverse grayscale and exclusion regions as blue contours. The analysis FoV for one observation is shown as a black circle with 2 GLYPH<14> radius and a black cross at the observation pointing position. The red rings illustrate the regions in which the background is estimated for two positions in the FoV (illustrated as red squares). Only regions in the ring inside the FoV and outside exclusion regions are used for background estimation. For the position in the lower right, the ring was adaptively enlarged to ensure an adequate background estimate (see text).  \n<!-- image -->  \nlarge uncertainties in the background estimation in regions without, or with small, exclusion regions where smaller ring radii are feasible.  \nArticle number, page 6 of 79  \nTo address the limitations of the standard method, we do not use a static ring geometry but rather adaptively change the inner and outer ring radii, as illustrated in Fig. 3, depending on the exclusion regions present in a given FoV. For a given test position within a FoV, we begin with a minimum inner ring radius of 0 : 7 GLYPH<14> and constant ring thickness 0 : 44 GLYPH<14> and enlarge the inner radius if a large portion of the ring area overlaps with exclusion regions. We do this until the acceptance integrated in the ring (but outside exclusion regions) is more than four times the acceptance integrated at the test position. A maximum outer radius of 1 : 7 GLYPH<14> avoids large uncertainties in the acceptance toward the edge of the FoV.', '3.3. Significance maps': "We produced significance maps to determine the exclusion regions (see Sect. 3.2.2). For each grid position ( '; b ) in a significance map, we counted the number of photon candidates N ON in the circular ON region, defined a priori by the correlation radius R c. We determined the background level by counting the number of photon candidates N OFF in the ring centered at ( '; b ). The background normalization factor is GLYPH<11> GLYPH<17> GLYPH<24> ON =GLYPH<24> OFF, where GLYPH<24> ON is the integral of the acceptance map within R c and GLYPH<24> OFF is the integral of the acceptance map within the ring. The number of excess events N GLYPH<13> within R c is then  \nN GLYPH<13> = N ON GLYPH<0> GLYPH<11> N OFF : (1)  \nWe computed the significance of this GLYPH<13> -ray excess according to Eq. 17 of Li & Ma (1983) without correcting further for trials.", '3.4. High-level maps': 'We can derive additional high-level maps based on the measurement of N GLYPH<13> within a given R c and the instrument response functions. In this work, we computed flux, flux error, sensitivity, and upper limit maps, starting from the formula  \nF = N GLYPH<13> N exp Z E 2 E 1 GLYPH<30> ref( E ) d E ; (2)  \nwhere F is the integral flux computed between the energies E 1 and E 2, N GLYPH<13> is the measured excess, and N exp is the total predicted number of excess events, also called exposure (see Sect. 3.4.1).', '3.4.1. Exposure maps': 'The exposure N exp in Eq. 2 is given by  \nN exp GLYPH<17> E = X R 2 runs T R Z 1 E min GLYPH<30> ref( E r) A e GLYPH<11> ( E r ; q R) d E r : (3)  \nHere, E r is the reconstructed energy, T R is the observation livetime, q R symbolizes the observation parameters for a specific run (zenith, o GLYPH<11> -axis, and azimuth angle; pattern of telescopes participating in the run; and optical e GLYPH<14> ciencies); A e GLYPH<11> is the e GLYPH<11> ective area obtained from MC simulations, which is assumed constant during a 28 min run; and E min is the safe threshold energy appropriate for the observation (as described in Sect. 3.1). We computed the quantity N exp for each position in the sky to create the expected GLYPH<13> -ray count map, also referred to as the exposure map E in the following. The function GLYPH<30> ref( E ) is the reference di GLYPH<11> erential photon number GLYPH<13> -ray source flux, assumed to be following a power law (PL) with a predefined spectral index, i.e.,  \nGLYPH<30> ref( E ) = GLYPH<30> 0 ( E = E 0) GLYPH<0> GLYPH<0> : (4)', '3.4.2. Flux maps': 'In Eq. 2, the flux value F is completely determined by the scaling factor N GLYPH<13> = N exp once the spectral shape is fixed. We chose to use E 1 = 1 TeV and E 2 = 1 . We stress that E 1 is not the threshold energy used in the analysis, but the energy above which the integral flux is given. In Eq. 4, one can choose the flux normalization GLYPH<30> 0 arbitrarily, since it cancels out in the computation of the flux. We also chose the spectral index GLYPH<0> = 2 : 3 in the released maps to be compatible with the average index of known Galactic VHE GLYPH<13> -ray sources. To test the impact of this latter assumption, we performed tests that show that, on average, flux variations are less than 5% if the assumed spectral index is varied by GLYPH<6> 0 : 2 (our systematic uncertainty of the spectral index).  \nThe released flux maps contain values of integral flux above 1 TeV, calculated according to Eq. 2, in units of cm GLYPH<0> 2 s GLYPH<0> 1 . This should be interpreted as the flux of a potential source, assuming  \na spectrum GLYPH<30> ref( E ), that is centered on a given pixel position in the map and fully enclosed within R c.  \nFigures 1 and 34 show two example flux maps computed with R c = 0 : 4 GLYPH<14> and 0 : 1 GLYPH<14> , respectively. The maps contain nonzero values only in regions in which the sensitivity is better than 2.5% Crab to prevent very large (positive and negative) values due to statistical fluctuations in low-exposure regions.', '3.4.3. Flux error and upper limit maps': 'Statistical uncertainties on the flux were computed by replacing N GLYPH<13> in Eq. 2 by N GLYPH<6> 1 GLYPH<27> GLYPH<13> , which are the upper and lower boundaries of the measured excess for a 68% confidence level. Those errors were computed with a Poisson likelihood method described in Rolke et al. (2005), using the same N ON and N OFF integrated within the circle of radius R c used when computing the excess maps. The values reported in the flux-error maps are the average of the upper and lower error bars.  \nSimilarly, an upper-limit map can be calculated by replacing N GLYPH<13> in Eq. 2 by N UL GLYPH<13> , that is, the upper limit on the excess found for a predefined confidence level of 95%; we used the same profile likelihood method as for the error bar.', '3.4.4. Sensitivity maps': "The sensitivity is defined as the minimal flux needed for a source with the assumed spectrum and fully contained within the correlation circle R c to be detected above the background at 5 GLYPH<27> statistical significance. Alternatively this can be thought of as a measure of ˆ N GLYPH<13> , the number of photons needed to reach such a significance level above the background determined by N OFF and GLYPH<11> . To compute the sensitivity map, N GLYPH<13> in Eq. 2 is replaced by ˆ N GLYPH<13> , which is determined by numerically solving Eq. 17 of Li & Ma (1983) for N ON (related to ˆ N GLYPH<13> by Eq. 1 above). We note that possible background systematics are not taken into account in this computation.  \nThe point-source sensitivity level reached by H.E.S.S. at all points in the HGPS data set is depicted in Fig. 4, where a projection of the sensitivity map along Galactic longitude at a Galactic latitude of b = 0 GLYPH<14> is also shown. It is typically at the level of 1% to 2% Crab. The deepest observations were obtained around interesting objects for which additional pointed observations were performed. Examples include the Galactic center region (around ' = 0 GLYPH<14> , where the best sensitivity of GLYPH<24> 0 : 3% Crab is reached), the Vela region ( ' = 266 GLYPH<14> ), the regions around HESS J1825 GLYPH<0> 137 and LS 5039 ( ' = 17 GLYPH<14> ), or around HESS J1303 GLYPH<0> 631 and PSR B1259 GLYPH<0> 63 ( ' = 304 GLYPH<14> ).  \nSimilarly, the sensitivity values along Galactic latitude for two values of longitude are shown in Fig. 11. For most of the surveyed region, the sensitivity decreases rapidly above j b j > 2 GLYPH<14> due to the finite FoV of the H.E.S.S. array and the observation pattern taken, except for a few regions, such as at ' = 0 GLYPH<14> where high  \nFig. 4. HGPS point-source sensitivity map and profile along the Galactic plane at a Galactic latitude b = 0 GLYPH<14> . The sensitivity is given in % Crab, for a correlation radius R c = 0 : 1 GLYPH<14> , assuming a spectral index GLYPH<0> = 2 : 3. This sensitivity is computed under the isolated point source assumption and is thus better than the actual sensitivity achieved for the HGPS source catalog (see Sect. 4.12).  \n<!-- image -->  \nlatitude observations were performed (see Sect. 2). The best sensitivity is obtained around b = GLYPH<0> 0 : 3 GLYPH<14> , reflecting the H.E.S.S. observation strategy; the latitude distribution of the sources peaks in this region.  \nWe note that the sensitivity shown in Fig. 4 does not correspond to the completeness of the HGPS source catalog. One major e GLYPH<11> ect is that the HGPS sensitivity is dependent on source size; it is less sensitive for larger sources, as shown in Fig. 13 and discussed at the end of Sect. 5.3. Other e GLYPH<11> ects that reduce the e GLYPH<11> ective sensitivity or completeness limit of HGPS are the detection threshold, which corresponds to GLYPH<24> 5 : 5 GLYPH<27> ; the large-scale emission model; and source confusion, as discussed in the following Sect. 4", '4.1. Introduction and overview': 'The HGPS source catalog construction procedure intends to improve upon previous H.E.S.S. survey publications both in sensitivity and homogeneity of the analysis performed. The previous iteration, the second H.E.S.S. survey paper of 2006 (Aharonian et al. 2006d), used a 230 h data set with inhomogeneous exposure that was limited to the innermost region of the Galaxy. This survey detected a total of 14 sources by locating peaks in significance maps on three di GLYPH<11> erent spatial scales: 0 : 1 GLYPH<14> , 0 : 22 GLYPH<14> , and 0 : 4 GLYPH<14> . It then modeled the sources by fitting two-dimensional symmetric Gaussian morphological models to determine the position, size and flux of each source, using a Poissonian maximumlikelihood method.  \nSince 2006, H.E.S.S. has increased its exposure tenfold and enlarged the survey region more than twofold, while also improving the homogeneity of the exposure. As illustrated in the  \nupper panel of Fig. 5, the data now show many regions of complex emission, for example, overlapping emission of varying sizes and multiple sources with clearly non-Gaussian morphologies. Apart from discrete emission, the Galactic plane also exhibits significant emission on large spatial scales (Abramowski et al. 2014a). For these reasons, we needed to develop a more complex analysis procedure to construct a more realistic model of the GLYPH<13> -ray emission in the entire survey region. Based on this model, we compiled the HGPS source catalog.  \nWe first introduce the maximum-likelihood method used for fitting the emission properties (Sect. 4.2). Next, we describe the H.E.S.S. point spread function (PSF; Sect. 4.3) and the TS maps (Sect. 4.4), which are two important elements in the analysis and catalog construction. The procedure is then as follows:  \n- 1. Cut out the Galactic center (GC) region and shell-type supernova remnants from the data set because of their complex morphologies (Sect. 4.5).\n- 2. Model the large-scale emission in the Galactic plane globally (Sect. 4.6).\n- 3. Split the HGPS region into manageable regions of interest (ROIs) (Sect. 4.7).\n- 4. Model the emission in each ROI as a superposition of components with Gaussian morphologies (Sect. 4.8).\n- 5. Merge Gaussian components into astrophysical VHE GLYPH<13> -ray sources (Sect. 4.9).\n- 6. Determine the total flux, position, and size of each GLYPH<13> -ray source (Sect. 4.10).\n- 7. Measure the spectrum of each source (Sect. 4.11).  \n- 8. Associate the HGPS sources with previously published H.E.S.S. sources and multiwavelength (MWL) catalogs of possible counterparts (Sect. 5.1).', '4.2. Poisson maximum-likelihood morphology fitting': "To detect and characterize sources and to model the large-scale emission in the Galactic plane, we used a spatially-binned likelihood analysis based on the following generic model:  \nN Pred = N Bkg + PSF GLYPH<3> ( E GLYPH<1> S ) ; (5)  \nwhere N Pred represents the predicted number of counts, N Bkg the background model created with the adaptive ring method (described in Sect. 3.2.3), E the exposure map (see Eq. 3 in Sect. 3.4.2), and S a two-dimensional parametric morphology model that we fit to the data. Additionally, we took into account the angular resolution of H.E.S.S. by convolving the flux model with a model of the PSF of the instrument.  \nAssuming Poisson statistics per bin, the maximumlikelihood fit then minimizes the Cash statistic (Cash 1979),  \nC = 2 X i GLYPH<0> Mi GLYPH<0> Di log Mi GLYPH<1> ; (6)  \nwhere the sum is taken over all bins i , and Mi (model) represents the expected number of counts according to Eq. 5 and Di (data) the actual measured counts per bin.  \nTo determine the statistical significance of a best-fit source model compared to the background-only model, we use a likelihood ratio test with test statistic TS. This is defined by the likelihood ratio or equivalently as the di GLYPH<11> erence in C between both hypotheses,  \nTS = C 0 GLYPH<0> C S ; (7)  \nwhere C 0 corresponds to the value of the Cash statistic of the background-only hypothesis and C S the best-fit model that includes the source emission.  \nFor a large number of counts, according to Wilks' theorem (Wilks 1938), TS is asymptotically distributed as GLYPH<31> 2 N , where N is the number of free parameters defining the flux model. In this limit, the statistical significance corresponds approximately to sign(Flux) GLYPH<1> p j TS j , where the sign of the best-fit flux is needed to allow for negative significance values in regions where the number of counts is smaller than the background estimate (e.g., due to a statistical downward fluctuation).  \nWe performed the modeling and fitting described above in Eqs. 5, 6, and 7 in pixel coordinates using the HGPS maps in Cartesian projection. Spatial distortion of flux models are negligible as a result of the projection from the celestial sphere because the HGPS observations only cover a latitude range of j b j 6 3 GLYPH<14> . We implemented the analysis in Python using As-  \ntropy version 1.3 (Astropy Collaboration et al. 2013), Sherpa version 4.8 (Freeman et al. 2001), and Gammapy version 0.6 (Donath et al. 2015; Deil et al. 2017).", '4.3. Point spread function': 'For HGPS, the PSF was computed for a given sky position assuming a power-law point source with a spectral index of 2.3 (average index of known VHE GLYPH<13> -ray sources) and assuming rotational symmetry of the PSF. Since the H.E.S.S. PSF varies with GLYPH<13> -ray energy and observing parameters such as the number of participating telescopes, zenith angle, and o GLYPH<11> set angle in the field of view, an e GLYPH<11> ective PSF corresponding to the HGPS survey counts maps was computed by applying the same cuts (especially safe energy threshold) and exposure weighting the PSF of contributing runs (i.e., within the FoV of 2 GLYPH<14> ). The per-run PSF was computed by interpolating PSFs with similar observation parameters, using precomputed lookups from MC EAS simulations. All computations were carried out using twodimensional histograms with axes GLYPH<18> 2 , where GLYPH<18> is the o GLYPH<11> set between the MC source position and the reconstructed event position, and log( Er ), where Er is the reconstructed event energy; at the very end, the integration over energy was performed, resulting in a one-dimensional histogram with axis GLYPH<18> 2 , which was fitted by a triple-exponential analytical function to obtain a smooth distribution,  \nd P d GLYPH<18> 2 ( GLYPH<18> 2 ) = 3 X i = 1 Ai exp 0 B B B B @ GLYPH<0> GLYPH<18> 2 2 GLYPH<27> 2 i 1 C C C C A ; (8)  \nwhere P is the event probability, and Ai and GLYPH<27> i are the weights and widths of the corresponding components, respectively. This ad hoc model corresponds to a triple-Gaussian, two-dimensional, PSF model when projected onto a sky map.  \nFor the HGPS catalog, the 68% containment radius of the PSF model adopted is typically GLYPH<18> GLYPH<24> 0 : 08 GLYPH<14> and varies by approximately GLYPH<6> 20% at the locations of the HGPS sources. For observations with large FoV o GLYPH<11> sets, the 68% containment increases by almost a factor of two to GLYPH<18> GLYPH<24> 0.15 GLYPH<14> , which is mostly relevant for high Galactic latitude sources at the edge of the HGPS survey region. The HGPS PSF has a 95% containment radius of GLYPH<18> GLYPH<24> 0.2 GLYPH<14> and approximately varies by GLYPH<6> 20% at the locations of the HGPS sources. The PSF at large FoV o GLYPH<11> sets (corresponding to high-GLAT regions in the survey map) is more tail heavy; there the 95% to 68% containment radius ratio increases from GLYPH<24> 2.5 up to 4. Section 4.10.2 discusses systematic uncertainties related to the PSF model in connection with upper limits on source sizes.', '4.4. Test statistics maps': 'In addition to the standard Li & Ma significance maps described in Sect. 3.3, we also used TS maps in the analysis. The TS de-  \nFig. 5. Illustration of the catalog model construction in the region of 350 GLYPH<14> to 328 GLYPH<14> in Galactic longitude. The upper panel shows the GLYPH<13> -ray excess counts smoothed by the PSF, the middle panel the PSF-convolved and smoothed excess model, and the lower panel the significance map of the residuals for a point-like source hypothesis (given in sign(Flux) p TS). The middle panel shows examples of the steps taken in the excess map modeling part of the source catalog procedure (see Sect. 4 for details). It starts by cutting out shell-type supernova remnants (SNRs; RXJ1713.7 GLYPH<0> 3946 and the SNR candidate HESS J1614 GLYPH<0> 518 in this region) and by assuming a fixed large-scale emission component. Then a multiGaussian model was fitted with the significant components shown in the middle panel as thin transparent circles. Some of these were discarded and are not part of the emission attributed to HGPS catalog sources. White circles show examples of single-component as well as multicomponent sources. For a complete overview of all analysis regions (ROIs) and excluded sources, see Fig. 39.  \n<!-- image -->  \nnotes the likelihood ratio of the assumed source hypothesis versus the null hypothesis (i.e., background only) for every position (pixel) in the map. We computed these maps assuming various spatial templates: a point-like source morphology (i.e., PSF only), and PSF-convolved Gaussian morphologies with widths 0 : 05 GLYPH<14> , 0 : 10 GLYPH<14> , and 0 : 20 GLYPH<14> . During the computation of each map, at the center of each map pixel, we performed a single-parameter likelihood fit of the amplitude of the template, according to Eq. 5. We then filled the map with the TS value defined in Eq. 7.  \nWeused the resulting TS maps primarily to compute residual maps and residual distributions. The main advantage over stan-  \ndard Li & Ma significance maps is that source morphology and PSF information can be taken into account. Additionally, this paper uses TS maps when presenting sky maps because they contain uniform statistical noise everywhere in the map. In contrast, flux or excess maps that are smoothed with the same spatial templates still show increased noise in regions of low exposure. We implemented the TS map algorithm available in Gammapy; see also Stewart (2009) for a more detailed description of TS maps.', '4.5. Sources not reanalyzed': "H.E.S.S. observations have revealed many sources with complex morphology, e.g., RX J0852.0 GLYPH<0> 4622 (also known as Vela Junior), which has a very pronounced shell-like structure (H.E.S.S. Collaboration et al. 2016a), or the Galactic center region, which has multiple point-sources embedded in a very elongated ridge-like emission (H.E.S.S. Collaboration et al. 2017g). Dedicated studies model such regions of emission using complex parametric models, for example, model templates based on molecular data, shell-like models, asymmetric Gaussian models, and combinations thereof. It is challenging to systematically model the emission across the entire Galactic plane using these more complex models, which tend to yield unstable or non-converging fit results because of the large number of free and often poorly constrained parameters. This can be especially problematic in ROIs with multiple, complex overlapping sources.  \nGiven the di GLYPH<14> culties with modeling complex source morphologies, we decided to restrict the HGPS analyses to a symmetrical Gaussian model assumption and exclude all firmly identified shell-like sources and the very complex GC region from reanalysis. A complete list of the ten excluded (or cut-out) sources in the HGPS region is given in Table 1. The table also contains four sources that were not significant in the current HGPS analysis but were found to be significant in other dedicated, published analyses; these cases are discussed in detail in Sect. 5.4.3. We refer to these 14 sources in total listed in Table 1 as 'EXTERN' HGPS sources and have included these sources in the HGPS source catalog because we wanted to give a complete list of sources in the HGPS region. We also have these sources included in the various distributions, histograms, and other plots exploring the global properties of the HGPS sources in Sect. 5.3. The morphological and spectral parameters for those sources were adapted from the most recent H.E.S.S. publication (listed in Table 1). 4", '4.6. Large-scale emission model': 'We previously demonstrated that there exists VHE GLYPH<13> -ray emission that is large scale and di GLYPH<11> use along the Galactic plane (Abramowski et al. 2014a). In that paper, we constructed a mask to exclude the regions of the plane where significant emission was detected. The latitude profile of excess GLYPH<13> -rays outside this mask clearly showed the presence of significant large-scale GLYPH<13> -ray emission. We do not extend the analysis of this di GLYPH<11> use emission  \nany further here. Whether the emission originates from interactions of di GLYPH<11> use cosmic rays in the interstellar medium or from faint, unresolved GLYPH<13> -ray sources (or a combination thereof) is not investigated. Instead, we take a pragmatic approach and model the large-scale emission present in the HGPS empirically as described in the following.  \nThe presence of a large-scale component of GLYPH<13> -ray emission along the Galactic plane complicates the extraction of the Gaussian GLYPH<13> -ray source components. This large-scale emission can mimic the presence of spurious degree-scale sources in some regions of the plane and it also tends to broaden the Gaussian components that describe otherwise well-defined sources. It is therefore necessary to model the large-scale GLYPH<13> -ray emission to measure the flux and morphology of the HGPS sources more accurately.  \nTo do so, we built an empirical surface brightness model of the large-scale emission (see Fig. 6), where the latitude profile is Gaussian and defined by three parameters: the peak position in latitude, the width, and amplitude of the Gaussian. We estimated the parameters using a maximum-likelihood fit in regions where no significant emission is measurable on small scales, i.e., outside the exclusion regions defined for the ring background model, taking exposure into account. Regardless of the physical origin of the large-scale emission, it is likely to be structured along the plane and not constant.  \nTo estimate the variable parameters of the model, we fit the Gaussian parameters in rectangular regions of width 20 GLYPH<14> in longitude and height 6 GLYPH<14> in latitude. We excluded all pixels inside the standard exclusion regions used to produce the background maps (see Sect. 3.2). The Gaussian parameters were dependent on the size of both the exclusion regions and rectangular regions. We found that the typical variations were GLYPH<24> 25%. To obtain a smooth sampling of the variations, we followed a sliding-window approach, distributing the centers of the rectangular regions every 2 : 5 GLYPH<14> in longitude and interpolating between these points.  \nThe maximum-likelihood fit compares the description of the data between the cosmic-ray (CR) background only and the CR background plus the model. We used the likelihood ratio test to estimate the significance of adding the large-scale component in each 20-deg-wide window, finding it to be larger than 3 GLYPH<27> (TS di GLYPH<11> erence of 9) over most of the HGPS region. Figure 6 shows the resulting best-fit Gaussian parameters together with the associated uncertainty intervals estimated from the likelihood error function. After this fit, we froze the parameters of the model for use in the GLYPH<13> -ray source detection and morphology fitting procedure.  \nWhile the approach presented here provides an estimate of the large-scale emission present in the HGPS maps, it does not comprise a measurement of the total Galactic di GLYPH<11> use emission (see discussion in Sect. 5.2).  \nTable 1. Fourteen EXTERN sources in the HGPS catalog, i.e., VHE sources in the HGPS region previously detected by H.E.S.S. that were not reanalyzed in this paper. For each source, we list the reason why it was not reanalyzed and give the reference that was used to fill the parameters in the HGPS source catalog. See Sect. 4.5 and for sources not significant in the HGPS analysis also Sect. 5.4.3.', '4.7. Regions of interest': 'To search for sources, we divided the whole HGPS region into smaller overlapping ROIs. This was necessary to limit both the number of simultaneously fit parameters and the number of pixels involved in the fit.  \nWe manually applied the following criteria to define the ROIs:  \n- (a) All significant emission (above 5 GLYPH<27> ) in the HGPS region should be contained in at least one ROI.\n- (b) No significant emission should be present close to the edges of an ROI.\n- (c) The width of each ROI should not exceed GLYPH<24> 10 GLYPH<14> in longitude to limit the number of sources involved in the fit.\n- (d) ROIs should cover the full HGPS latitude range from GLYPH<0> 5 GLYPH<14> to 5 GLYPH<14> .  \nIn cases in which criterion (b) could not be fulfilled, we excluded the corresponding emission from the ROI and assigned it to a di GLYPH<11> erent, overlapping ROI. Figure 39 illustrates the boundaries of the 18 ROIs defined with these criteria. Some of the ROIs show regions without any exposure; these regions were masked out and ignored in the subsequent likelihood fit.', '4.8. Multi-Gaussian source emission model': "After excluding shell-type supernova remnants (SNRs) and the GC region from reanalysis and adding a model for large-scale emission to the background, we modeled all remaining emission as a superposition of Gaussian components. We took the following model as a basis:  \nN Pred = N Bkg + PSF GLYPH<3> 0 B B B B B @ E GLYPH<1> X i S Gauss ; i 1 C C C C C A + E GLYPH<1> S LS ; (9)  \nArticle number, page 12 of 79  \nwhere N Pred corresponds to the predicted number of counts, N Bkg to the number of counts from the background model, S LS the contribution of the large-scale emission model, P i S Gauss ; i the sum of the Gaussian components, and E the exposure as defined in Eq. 3.  \nFor a given set of model parameters, we integrated the surface brightness distribution S over each spatial bin, multiplied it by the exposure E , and convolved it with the PSF to obtain the predicted number of counts per pixel. For every ROI, we took the PSF at the position of the brightest emission and assumed it to be constant within the ROI.  \nFor the Gaussian components, we chose the following parametrization:  \nS Gauss( r j GLYPH<30>; GLYPH<27> ) = GLYPH<30> 1 2 GLYPH<25>GLYPH<27> 2 exp GLYPH<0> r 2 2 GLYPH<27> 2 ! ; (10)  \nwhere S Gauss is the surface brightness, GLYPH<30> the total spatially integrated flux, and GLYPH<27> the width of the Gaussian component. The o GLYPH<11> set r = p ( ' GLYPH<0> ' 0) 2 + ( b GLYPH<0> b 0) 2 is defined with respect to the position ( ' 0 ; b 0) of the component measured in Galactic coordinates.  \nWe conducted the manual fitting process following a stepby-step procedure. Starting with one Gaussian component per ROI, we added Gaussian components successively and refit all of the parameters simultaneously until no significant residuals were left. In each step, we varied the starting parameters of the fit to avoid convergence toward a local minimum. The significance of the tested component was estimated from  \nTS = C (with component) GLYPH<0> C (best solution without component) :  \n(11)  \nWe considered the component to be statistically significant and kept it in the model when the TS value exceeded a threshold of  \nFig. 6. Distribution of the fit large-scale emission model parameters with Galactic longitude. The first panel gives the peak brightness of the large-scale emission model in units of 10 GLYPH<0> 9 cm GLYPH<0> 2 s GLYPH<0> 1 sr GLYPH<0> 1 ( GLYPH<25> 1 : 3% Crab deg GLYPH<0> 2 ). The second panel shows the peak position of the Gaussian along the Galactic latitude axis in degrees and the third panel shows the width ( GLYPH<27> ) of the Gaussian in degrees. The solid lines are the result of fitting each set of parameters every 2 : 5 GLYPH<14> in longitude and interpolating. The light blue bands show the 1 GLYPH<27> error region obtained from the covariance matrix of the likelihood function. The lower panel illustrates the 20 GLYPH<14> wide sliding-window method (red rectangle) that was used to determine the large-scale emission model in areas (shown in light blue) where the HGPS sensitivity is better than 2.5% Crab but outside exclusion regions (shown in dark blue); this is explained in further detail in the main text.  \n<!-- image -->  \nTS = 30. The probability of having one false detection in the HGPS survey from statistical background fluctuations is small ( p = 0 : 03). This number was determined by simulating 100 HGPS survey counts maps as Poisson-fluctuated background model maps, followed by a multi-Gaussian peak finding method, resulting in three peaks with TS GLYPH<21> 30. However, we note that this assessment of expected false detections lies on the assumption that the hadronic background as well as the large-scale and source gamma-ray emission model are perfect. In HGPS, as in any other Galactic plane survey with complex emission features, this is not the case. Several components with TS GLYPH<21> 30 are not confirmed by the cross-check analysis (see Sect. 4.9).  \nThe definition of TS above di GLYPH<11> ers slightly from the definition given in Eq. 7. For a single, isolated component, both values are identical. However, if a second, overlapping component exists, some of the emission of the first source is modeled by the second source, reducing the significance of the first. We therefore estimated the significance of a component from the TS di GLYPH<11> er-  \nence in the total model of the ROI and not from the TS di GLYPH<11> erence compared to the background-only model.  \nApplied to real data, we found a total of 98 significant Gaussian components using this procedure and TS threshold. Figure 7 depicts the residual p TS distributions over the entire HGPS region. These distributions demonstrate that there is approximate agreement with a normal Gaussian distribution; in particular, we find no features above the p TS = p 30 detection threshold. Inherent imperfections in the background, large-scale emission models and source emission models lead to a slight broadening of the distributions with respect to a normal distribution, as expected.  \nFor reference, the 98 Gaussian components have been assigned identifiers in the format HGPSC NNN , where NNN is a three-digit number (counting starts at 1), sorted by right ascension (which is right to left in the survey maps). The complete list of components is provided in the electronic catalog table (see Table 6).  \n√  \nFig. 7. Residual significance distribution after taking the HGPS emission model into account (see Fig. 5, middle panel). The significance was computed using a Gaussian source morphology of size GLYPH<27> = 0 : 05 GLYPH<14> , 0.10 GLYPH<14> and 0.20 GLYPH<14> . A vertical line at p TS = p 30 is shown, corresponding to the detection threshold for the HGPS multi-Gaussian modeling. The sky region corresponding to this distribution includes pixels inside exclusion regions, except for the Galactic center and shell-type SNRs, which were not modeled for the HGPS (see Table 1, lower panel of Fig. 5 and Fig. 39).  \n<!-- image -->", '4.9. Component selection, merging, and classification': "We repeated the entire modeling procedure described in the previous section with a second set of maps produced with an independent analysis framework (see Sect. 2.3). Five of the 98 HGPS components were not significant in the cross-check analysis and were therefore discarded (see Fig. 5 and Table 6). Those components we labeled with Discarded Small in the column Component\\_Class of the FITS table.  \nWe observed two other side e GLYPH<11> ects of the modeling procedure. Firstly, very bright VHE sources, even some with centerfilled morphologies such as Vela X, decomposed into several Gaussian components, modeling various morphological details of the source. Figure 5 illustrates this e GLYPH<11> ect: there are two multicomponent sources shown. Therefore in cases where overlapping components were not clearly resolved into separate emission peaks, we merged them into a single source in the HGPS catalog. In total, we found 15 such multicomponent sources: ten consisting of two Gaussian components and five consisting of three Gaussian components. It would be intriguing to analyze the complex morphology of these multicomponent sources in greater detail, but this kind of analysis is beyond the scope of this survey paper. We labeled components that are part of a multicomponent source as Source Multi . We used the label Source Single, respectively, if there is only one component modeling the source.  \nThe second side e GLYPH<11> ect was that some of the Gaussian components appeared to have very large sizes coupled with very low surface brightness. We interpret these components as artifacts of  \nthe modeling procedure, which picks up additional di GLYPH<11> use GLYPH<13> -ray emission that is not covered by our simple large-scale emission model (Sect. 4.6). For example, as shown in Fig. 5, the emission around ' GLYPH<24> 345 GLYPH<14> initially comprised three model components: two components that clearly converged on the two discrete emission peaks visible in the excess map and one very large and faint component that appeared to be modeling large-scale emission along the Galactic plane in between the two and not clearly related to either of the two peaks. In total, we found ten such largescale components (see Table 6), which we discarded and did not include in the final HGPS source catalog as they are likely lowbrightness di GLYPH<11> use emission. We labeled this class of components as Discarded Large in the component list.", '4.10.1. Position, size, and flux': 'For HGPS sources that consist of several components, we determined the final catalog parameters of the sources as follows:', 'Flux': 'Wefound the flux of 42 (out of 45) sources to be compatible with the previous published value, according to Eq. 29.  \nThe average di GLYPH<11> erence in flux for extended sources was 3 % with a width of 43 % for the distribution of values. While the average value is compatible with previous analyses, we still found a large scatter (albeit compatible to the systematic and statistical errors) of the distribution.  \nA fair comparison between flux values obtained with the current method and earlier analyses proved to be di GLYPH<14> cult again because of fundamental di GLYPH<11> erences between the methods used. In previous publications, aperture photometry was mostly used, while in this analysis the main flux measurement was based on a model fit, taking the PSF and morphology of the source and large-scale emission into account. Flux estimate di GLYPH<11> erences with these two methods are shown in Fig. 9 (both measures from the HGPS analysis, not with respect to previous publications). Many of the di GLYPH<11> erences in spectra and fluxes measured in the HGPS analysis and previous publications are the result of changes in the spectral extraction region (position and size).', 'Position': 'We found the position of 43 (out of 45) sources to be compatible with the previously published value, according to Eq. 29. For point-like sources we found an average shift of 0 : 02 GLYPH<6> 0 : 01 deg, while for extended sources the value was 0 : 06 GLYPH<6> 0 : 05 deg. Both values agree well with the expected scatter considering the statistical and systematic uncertainties on the measurements. As an additional check, we also verified that the positions of the identified GLYPH<13> -ray binaries (known point sources) HESS J1826 GLYPH<0> 148 and HESS J1302 GLYPH<0> 638 are in good agreement (within 40 00 ) with the reference positions of the corresponding objects LS 5039 and PSR B1259 GLYPH<0> 63 as listed in SIMBAD 10 .', 'Size': "Comparing the sizes of the extended sources we found 30 (out of 35) sources to be compatible with the previously published value. The average size di GLYPH<11> erence for the extended sources was on the order of GLYPH<24> 18 %, the distribution of values having a width of GLYPH<24> 40 %. This indicates that with the current analysis we measured slightly larger sizes of the sources on average, but the distribution is dominated by a large scatter. We expect the scatter to result mainly from di GLYPH<11> erences in the analysis procedure. Previous analyses mainly fitted single Gaussian morphologies, while in this analysis we allowed for multiple Gaussian components. Further di GLYPH<11> erences are the addition of the large-scale emission model and the systematic modeling of emission from neighboring sources.  \nPrevious publications found seven sources to be compatible with a point-like source. In the current analysis we found all these sources to be compatible with a point-like source again. Additionally, we identified the following three cases that are compatible with a point-like source according to Eq. 18, which were previously found to be extended:  \n- 1. For HESS J1427 GLYPH<0> 608 we measured a size of 0 : 048 GLYPH<6> 0 : 009 GLYPH<14> , compared to 0 : 063 GLYPH<6> 0 : 010 GLYPH<14> in Aharonian et al. (2008a). This  \nsource is a borderline case that just meets our criterion for a point-like source.  \n- 2. For HESS J1714 GLYPH<0> 385 we found a size of 0 : 034 GLYPH<6> 0 : 011 GLYPH<14> compared to 0 : 067 GLYPH<6> 0 : 017 GLYPH<14> in Aharonian et al. (2008c). With the current analysis, a smaller size was found because underlying emission was modeled by separate emission components (see Fig. 5).\n- 3. We now measure the size of HESS J1808 GLYPH<0> 204 to be 0 : 058 GLYPH<6> 0 : 014 GLYPH<14> (consistent with point-like, in the definition of Eq. 18), compared to the previously measured size 0 : 095 GLYPH<6> 0 : 015 GLYPH<14> (extended) (Abdalla et al. 2016). This discrepancy is due to the HGPS's inclusion of a large-scale emission component that now models GLYPH<13> -ray excess previously accounted for in the source component itself.", '4.11. Source spectra': 'In the morphology fit, we did not take into account uncertainties in the PSF model. However, studies using H.E.S.S. data (e.g., Stycz 2016) have revealed a systematic bias on the size of pointlike extragalactic sources on the order of GLYPH<27> syst = 0 : 03 GLYPH<14> , so we have adopted this number as the systematic uncertainty of the PSF.  \nGiven a measured source extension GLYPH<27> Source and corresponding uncertainty GLYPH<1> GLYPH<27> Source, we used the following criterion to claim a significant extension beyond the PSF:  \nGLYPH<27> Source GLYPH<0> 2 GLYPH<1> GLYPH<27> Source > GLYPH<27> syst ; (17)  \ni.e., if the extension of a source is 2 GLYPH<1> GLYPH<27> Source beyond the systematic minimum GLYPH<27> syst. If this criterion is not met, we consider the source to be compatible with being point-like and define an upper limit on the source size as follows:  \nGLYPH<27> UL = max( GLYPH<27> syst ; GLYPH<27> Source + 2 GLYPH<1> GLYPH<27> Source) : (18)', '4.10.3. Localization': "The HGPS source location error is characterized by error circles with radius R GLYPH<11> at confidence levels GLYPH<11> = 0 : 68 and GLYPH<11> = 0 : 95, computed as  \nR GLYPH<11> = f GLYPH<11> GLYPH<2> q GLYPH<1> ' 2 stat + GLYPH<1> ' 2 syst + GLYPH<1> b 2 stat + GLYPH<1> b 2 syst : (19)  \nThe values GLYPH<1> ' stat and GLYPH<1> b stat are the statistical errors on Galactic longitude ' and latitude b , respectively, from the morphology fit. For the H.E.S.S. systematic position error, a value of GLYPH<1> ' syst = GLYPH<1> b syst = 20 00 = 0 : 0056 GLYPH<14> per axis was assumed, following the method and value in (Acero et al. 2010b).  \nAssuming a Gaussian probability distribution, the factor f GLYPH<11> is chosen as f GLYPH<11> = p GLYPH<0> 2 log(1 GLYPH<0> GLYPH<11> ) for a given confidence level GLYPH<11> (see Eq. 1 in Abdo et al. 2009b).", '4.10.4. Source naming': 'The 78 HGPS catalog sources have been assigned source names in the format HESS JHHMM GLYPH<6> DDd , where HHMM and GLYPH<6> DDd are the source coordinates in right ascension and declination, respectively. For new sources, the source name is based on the source location reported in this paper. For sources that had been assigned names in previous H.E.S.S. publications or conference presentations, the existing name was kept for the HGPS catalog, even if the position in the HGPS analysis would have led to a di GLYPH<11> erent name. Similarly, the source candidates (or hotspots, see Sect. 5.6.17) have been assigned names in the format HOTS JHHMM GLYPH<6> DDd .  \nAfter detection and subsequent morphological analysis of the sources, we measured a spectrum for each of the sources using an aperture photometry method. In this method we sum the ON counts within an aperture defined as a circular region centered on the best-fit position of each source. We fit a spectral model within that aperture using an ON-OFF likelihood method (Piron et al. 2001), where the OFF background is estimated using reflected regions defined on a run-by-run basis (Fomin et al. 1994; Berge et al. 2007). Based on the morphology model, we then corrected the measured flux for containment and contamination from other nearby sources. For the spectral analysis, we applied standard cuts, resulting in energy thresholds in the range 0.2-0.5 TeV, lower than the thresholds achieved using hard cuts in the detection and morphology steps. Figure 2 shows the variation of the threshold with longitude. In the following sections, we describe the spectral analysis process in more detail.', '4.11.1. Circular apertures and reflected region background estimate': 'The optimal choice for the size for the spectral extraction region is a balance between including a large percentage of flux from the source and limiting the contamination of the measurement by hadronic background events, large-scale emission, and other nearby sources. Following these requirements, we chose the aperture radius R spec as follows:  \n- -R spec = R 70 for 34 medium-size sources, where R 70 is the 70% containment radius measured on the PSF-convolved excess model image ( R70 in the catalog),\n- -minimum R spec = 0 : 15 GLYPH<14> for 21 small ( R70 < 0 : 15 GLYPH<14> ) sources,\n- -maximum R spec = 0 : 5 GLYPH<14> for 9 very large ( R70 > 0 : 5 GLYPH<14> ) sources.  \nAminimal aperture radius of 0 : 15 GLYPH<14> was imposed to make the measurement of the source spectrum more robust against systematic uncertainties of the PSF and the source morphology assumption.  \nThe aperture radius was limited to a maximum radius of R spec = 0 : 50 GLYPH<14> to limit the fraction of observations that cannot be used for the spectrum measurement because no background estimate could be obtained.  \nAs illustrated in Fig. 8, the background is estimated using the reflected region method (Fomin et al. 1994; Berge et al. 2007). For every spectral extraction region (ON region), corresponding OFF regions with the same shape and o GLYPH<11> set to the pointing position are chosen outside exclusion regions.  \nThe method works well for small, isolated GLYPH<13> -ray sources such as active galactic nuclei (AGNs) or the Crab nebula, where typically GLYPH<24> 10 OFF regions are found in every observation. This results in a well-constrained background, and all the exposure can  \nFig. 8. Illustration of reflected region background estimation for spectra (Sect. 4.11.1). The HGPS significance image is shown in inverse grayscale and exclusion regions as blue contours. The analysis FoV for one observation is shown as a black circle with 2 GLYPH<14> radius and a black cross at the observation pointing position. The non-filled red circle illustrates the ON region for spectral analysis; the filled red circles indicate the OFF regions.  \n<!-- image -->  \nbe used for the spectral measurement. Because of the high density of sources in the Galactic plane, large areas of emission are excluded and only few reflected regions can be found. This effectively results in a loss of exposure for the spectrum measurement compared to the map measurement. For the HGPS analysis this is a large problem because of the very extensive exclusion regions used: 64% of the livetime is lost for spectral analysis compared to the total available livetime that is used in the mapbased analysis. For each source, see the Livetime\\_Spec and Livetime information in the source catalog. In cases where the loss of exposure is very high, the background cannot be well constrained, which consequently results in spectral parameters that are not well constrained. The following sources are a GLYPH<11> ected by this issue:  \n- -Sources located in or near large exclusion regions (see Fig. 39). An area of width GLYPH<24> 2 GLYPH<14> is often excluded along the Galactic plane, and this covers a significant portion of the analysis FoV, which has a diameter of 4 GLYPH<14> .\n- -Sources with large ON regions.\n- -Sources observed with too small or too large o GLYPH<11> sets because they are located close to other sources that were covered with dedicated observations.', '4.11.2. Flux containment and contamination correction': 'By construction and because of additional e GLYPH<11> ects such as PSF leakage or source morphologies featuring tails, the spectral extraction region does not contain the full flux of the source. Additionally, the large-scale emission model and other nearby Gaussian components bias the flux measurement within the spectral region. Based on this emission model, we separate the contributions from the di GLYPH<11> erent components and derive a correction factor for the spectral flux measurement.  \nThe total flux in the spectral measurement region is  \nF ON Total = F ON Source + F ON LS + F ON Other ; (20)  \nwhere F ON Source is the contribution from the source itself, F ON LS is the contribution from the large-scale emission model, and F ON Other is the contribution from nearby sources and other, discarded Gaussian emission components.  \nAssuming F Source is the flux measurement from the morphology fit, we define the correction factor as  \nC Correction = F Source = F ON Total : (21)  \nTo summarize the contributions from the large-scale emission model and other sources in close (angular) proximity, we define a quantity called contamination. This quantity measures the fraction of flux within the spectral region that does not originate from the source itself and is written as  \nC Contamination = F ON LS + F ON Other F ON Total : (22)  \nAdditionally, we define the containment of a source as the ratio between the flux of the source within the spectral measurement region F ON Source (taking the morphology model into account) and the total flux obtained from the morphology fit F Source as follows:  \nC Containment = F ON Source = F Source : (23)  \nThe HGPS catalog provides all the quantities mentioned in this section, and all aperture-photometry based flux measurements in the HGPS catalog (see Table 5) are corrected by the factor given in Eq. 21 (see Sect. 4.11.3 and 4.11.4).  \nWe note that this region-based spectral analysis method with a single integral flux correction factor assumes energyindependent source morphology. The spectra obtained for sources with energy-dependent morphology does not correspond to the correct total emission spectra of the sources. Currently energy-dependent morphology has been clearly established for two sources (HESS J1303 GLYPH<0> 631 and HESS J1825 GLYPH<0> 137), and there are hints of energy-dependent morphology for a few more. Furthermore, using an integral flux correction factor is not fully correct because the HGPS PSF is somewhat dependent on en-  \nergy (smaller PSF at higher energies). The resulting inaccuracy on HGPS spectral results is small, because we have chosen a minimal spectral aperture radius of 0.15 GLYPH<14> , which contains most of the emission for point sources at all energies. Generally, spectra for sources with large correction factors are likely to be less accurate, because the morphology models used to compute the correction are only approximations.', '4.11.3. Spectral model fit': 'We performed the spectral fits on the stacked 5 observations, using the ON-OFF Poisson likelihood function, referred to as the W statistic (WSTAT) in XSPEC 6 . For each observation, we applied a safe energy threshold (see Sect. 3.1) cut at low energies, and the maximum energy was chosen at the highest event energy in the stacked counts spectrum for the on region (resulting in a maximum energy of 30 TeV to 90 TeV). Energy dispersion was not taken into account via a matrix, but in an approximate way in which the e GLYPH<11> ective area is computed in such a way that it results in fully correct spectral results for power-law spectra with spectral index 2, and, given the good energy resolution of H.E.S.S., only small errors are made for other spectral shapes (Hoppe 2008a).  \nTo describe the spectral shape of the VHE GLYPH<13> -ray emission, we fit a PL model to the data, i.e.,  \nGLYPH<30> ( E ) = d N d E = GLYPH<30> 0 E E 0 ! GLYPH<0> GLYPH<0> ; (24)  \nwhere GLYPH<30> 0 is the di GLYPH<11> erential flux at a reference (pivot) energy E 0 and GLYPH<0> is the spectral index. In addition, we also fit an exponential cuto GLYPH<11> power-law (ECPL) model,  \nGLYPH<30> ( E ) = GLYPH<30> 0 E E 0 ! GLYPH<0> GLYPH<0> exp( GLYPH<0> GLYPH<21> E ) ; (25)  \nwhich additionally contains the inverse cuto GLYPH<11> energy GLYPH<21> = 1 = E cuto GLYPH<11> as a third, free parameter. The reference (pivot) energy E 0 is not a free parameter in either model; we compute this parameter on a source-by-source basis to minimize the correlation between the other fit parameters.  \nWe computed integral fluxes as  \nF ( E 1 ; E 2) = Z E 2 E 1 GLYPH<30> ( E ) d E ; (26)  \nusually for the energy band above 1 TeV, with integral flux errors computed using Gaussian error propagation. We computed  \nenergy fluxes for a given energy band as  \nG ( E 1 ; E 2) = Z E 2 E 1 E GLYPH<30> ( E ) d E : (27)  \nThe source catalog provides the PL fit results (see Table 5 for a description of columns) for every source and the ECPL parameters where the ECPL model is more likely (TS = W ECPL GLYPH<0> W PL > 9). All aperture-photometry based flux measurements are corrected by the factor given in Eq. 21.', '4.11.4. Flux points': 'Flux points are estimates of the di GLYPH<11> erential flux GLYPH<30> at a given set of reference energies Eref . To compute flux points for the HGPS catalog, we chose a method similar to that used for the Fermi -LAT catalogs (see, e.g., Sect. 5.3 in Acero et al. 2015a). For every source we selected a total number of six bins ( E 1 ; E 2) in reconstructed energy, logarithmically spaced between the safe energy threshold and a maximum energy of 50 TeV. The reference energy for the flux point estimation was set to the logarithmic bin center Eref = p E 1 E 2. The di GLYPH<11> erential flux GLYPH<30> was computed via a one-parameter likelihood fit (same method as described in Sect. 4.11.3), under the assumption of the global best-fit PL and using only the data within the bin of reconstructed energy ( E 1 ; E 2). An 1 GLYPH<27> asymmetric error on GLYPH<30> was computed from the likelihood profile, and for spectral points of small significance ( TS < 1), in addition an upper limit on GLYPH<30> was computed at 95% confidence level. All spectral point measurements in the HGPS catalog are corrected by the factor given in Eq. 21.', '4.12. Method discussion': 'The sensitivity profile and map shown in Fig. 4 were computed assuming a point-like source morphology and using the Li & Ma significance estimation. The likelihood fit method including the large-scale emission model component used for the catalog production fundamentally di GLYPH<11> ers from that. We qualitatively discuss below the most important di GLYPH<11> erences and their influence on the e GLYPH<11> ective sensitivity with which the catalog was produced.  \nIn Sect. 3.4.4, the sensitivity was defined as the minimum required flux for a source to be detected with a certain level of confidence. Assuming the source is extended, which applies to most of the Galactic sources found by H.E.S.S., the total flux of the source is distributed over a larger area on the sky. Given a fixed background level, the signal-to-noise ratio is decreased and the sensitivity scales with the size of the source as  \nF min( GLYPH<27> source) / q GLYPH<27> 2 source + GLYPH<27> 2 PSF ; (28)  \nwhere GLYPH<27> source is the size of the source and GLYPH<27> PSF the size of the PSF (Hinton & Hofmann 2009). It is constant for sources smaller  \nthan the PSF and increases linearly with source size for sources much larger than the PSF.  \nFor low surface brightness sources close to the Galactic plane, high levels of contamination (defined as in Eq. 22) from the large-scale emission model were observed. This e GLYPH<11> ectively reduces the sensitivity close to the Galactic plane and even caused a few previously detected H.E.S.S. sources to fall below the detection threshold (see also Sect. 5.4.3) chosen for the HGPS analysis. For sources far from the Galactic plane, however, the influence of the large-scale emission can be neglected.  \nSystematic and statistical background uncertainties, which are neglected in this analysis, bias the sensitivity for large, extended sources. Neglecting background fluctuations in the likelihood fit can lead to an overestimation of the significance of large sources, which can lead to unreliable detections of large emission components. In addition, the adaptive ring method (Sect. 3.2.3), which has a minimal inner ring radius of 0 : 7 GLYPH<14> , does not provide a reliable background estimate for those large emission components.  \nSystematic uncertainties of various origins a GLYPH<11> ect the spectral parameters of the sources. In addition to the transparency of the atmosphere, calibration, and event reconstruction (see Sect. 2), the analysis method itself can introduce uncertainties. In particular, the background and large-scale emission emission model, and the source extraction and measurement method (multi-Gaussian morphology and aperture photometry) influence the flux and spectral index measurement. We estimate the relative systematic uncertainties of the flux extracted from the maps (Sect. 3) and from the spectrum (Sect. 4.11) to be 30%; for the spectral index (Sect. 4.11) we estimate an absolute systematic uncertainty of 0.2. This estimate is based on the scatter seen in the cross-check analysis and other analyses (e.g., a source catalog extracted without a large-scale emission model component). For individual di GLYPH<14> cult sources (poor containment, large contamination, complex, and marginally significant morphology features), larger systematics may arise (see Sects. 5.4 and 5.5). We note that the systematic uncertainties quoted here are the same as in the previous HGPS publication (Aharonian et al. 2006d), and, as expected for a population of extended sources in the Galactic plane, these values are slightly larger than the systematic uncertainties previously estimated for isolated point-like sources such as the Crab nebula (Aharonian et al. 2006b).  \nAcomparison of the two methods presented in this paper for calculating HGPS source integral flux ( E > 1 TeV) was performed as a diagnostic test (see scatter plot in Fig. 9). The flux show on the x-axis is the total source flux estimated from the source morphology fit on the maps (given by Eqn. 12), assuming a power-law spectrum with index GLYPH<0> = 2 : 3. The flux estimate on the y-axis was obtained from a spectral analysis (see Eqn. 26 in Sect. 4.11), using a PL or ECPL spectral model assuption (best-fit model) and an aperture photometry method that  \nFig. 9. Comparison of integral source flux measurements above 1 TeV as calculated with two di GLYPH<11> erent methods. The flux estimate from maps is the total source flux according to the morphology model fit, assuming a spectral index of GLYPH<0> = 2 : 3 (the Flux\\_Map column in the catalog). The flux estimate from spectra is computed from the total source best-fit spectral model extracted using aperture photometry and aperture correction (the Flux\\_Spec\\_Int\\_1TeV column in the catalog). The gray band in the background illustrates a systematic uncertainty of 30% on the flux values.  \n<!-- image -->  \nincludes a containment and contamination correction according to the HGPS multi-Gaussian plus large-scale emission model. One can see that the two flux estimates agree very well for most sources within the statistical errors and the 30% flux systematic uncertainty that we quote above. There are exceptions however, which can to a large degree be attributed to di GLYPH<11> erences in the underlying morphology and spectral model assumptions of the two flux estimators. We note that when comparing either of these HGPS source flux estimates against the cross-check, the level of agreement is similar, but not quite as good (see Sect. 5.5 for a discussion of individual cases). When comparing against previous publications, the scatter is even larger (flux di GLYPH<11> erences up to a factor of 2 in a few cases), which can in many cases be understood to be the result of di GLYPH<11> erences in morphology model (of the source itself, of nearby overlapping sources, or the large-scale emission model) or the spectral extraction region; most previous publications did not apply containment or contamination corrections.', '5. Results and discussion': 'This section presents the results and a discussion of the HGPS based on the data set (Sect. 2), maps (Sect. 3), and catalog (Sect. 4).', '5.1. Source associations and firm identifications': 'Determining the physical nature of a VHE GLYPH<13> -ray source often requires detailed spectral and morphological characterization of the VHE emission and availability of complementary MWL information. Finding a likely counterpart of a point-like VHE source is generally easy thanks to the limited region of the sky to investigate. For an extended source, such as the vast majority of the HGPS sources, the procedure is often much more involved because of multiple spatial associations, unless the VHE morphology is similar to that observed at other wavelengths (e.g., for a large shell-type SNR).  \nWe therefore make a distinction between source associations and firm identifications of sources. The former is a list of astronomical objects, extracted from catalogs of plausible counterparts, which are are found to be spatially coincident with the HGPS source. When particularly solid evidence exists that connects one of these associated objects to the VHE emission, such as variability or shell-type morphology, we consider the HGPS source to be firmly identified.  \nIn Sect. 5.1.1 we first describe the systematic association procedure, followed by the discussion of the results of this search for plausible counterparts in Sect. 5.1.2. Finally we present the list of firmly identified HGPS sources in Sect. 5.1.3.', '5.1.1. Source association procedure': 'Our objective is to associate each HGPS source with plausible counterparts found among nearby objects in the most relevant counterpart catalogs (that is catalogs of objects already identified as VHE emitters such as SNRs and pulsar wind nebulae and high energy GLYPH<13> -ray sources, see Table 2). We search for these counterparts in a region larger than the nominal HGPS source size (its 68% containment radius); we associate all the objects whose cataloged positions are at an angular distance smaller than the source spectral extraction radius, R SPEC (Sect. 4.11.1; see also Fig. 35 GLYPH<11> .). This spatial criterion is motivated by the fact that often the origin of the relativistic particles is significantly o GLYPH<11> set from the VHE centroid, for example, when VHE pulsar wind nebulae (PWNe) are o GLYPH<11> set from energetic pulsars or extended well beyond their X-ray PWNe counterparts. We expect this procedure to be a GLYPH<11> ected by source confusion (multiple associations that are di GLYPH<14> cult to disentangle), especially for larger VHE sources.  \nThis criterion is a compromise between the number of spurious associations and the number of missed associations. A spurious association would be one with a counterpart that is physically unrelated to the HGPS source (e.g., a chance spatial coincidence in the same region of the sky). A missed association would be a real counterpart that is not selected by the procedure (e.g., a pulsar significantly o GLYPH<11> set from a VHE source could be missed even though it is known to generate a PWN). As a con-  \nsequence of this spatial criterion, larger sources naturally has a larger number of associated objects. The criterion is intended to be loose (inclusive) to minimize missed associations at the expense of including potentially spurious associations. Nonetheless, this procedure has certain limitations, for example, di GLYPH<14> culties in associating VHE emission with an SNR if the emission was produced in o GLYPH<11> set molecular clouds illuminated by cosmic rays that escaped from the SNR.  \nIn the following paragraphs, we briefly describe the catalogs used for the automatic association procedure applied to search for counterparts. We also describe a list of additional objects that have been associated with HGPS sources in previous publications but are not present in the counterpart search catalogs. We note that some of these catalogs contain a single, specific type of object (e.g., supernova remnants), whereas other catalogs contain multiple types of physical objects because they are the result of broad surveys at energies relevant to the HGPS (e.g., the Fermi -LAT catalogs).', 'High-energy GLYPH<13> -ray sources': 'We searched for associated high-energy (HE) GLYPH<13> -ray sources in the Fermi -LAT 2FHL source catalog (Ackermann et al. 2016) and the full 3FGL catalog (Acero et al. 2015a). The 2FHL catalog covers the 50 GeV to GLYPH<24> 2 TeV energy range, and the 3FGL catalog covers the 0.1-300 GeV range. They contain 44 and 352 sources in the HGPS region, respectively. We expect the Fermi -LAT catalogs to contain a significant number of HGPS sources. In the case of 2FHL, this is due to its energy range, which partially overlaps that of the HGPS, and its sensitivity, which reaches GLYPH<24> 3-4% Crab in the HGPS region (Ackermann et al. 2016). But even without such overlaps, we expect to find many Fermi -LAT associations, since many objects emit GLYPH<13> -rays following a spectrum that extends from the HE to the VHE range. Even for noncontinuous spectra we expect to find numerous associations, for example, when a pulsar emits GeV emission detected by Fermi -LAT and its wind nebula emits TeV emission detected by H.E.S.S.', 'Supernova remnants and pulsar wind nebulae': 'Supernova remnants and PWNe are among the most common particle accelerators in the Galaxy and are well-known VHE GLYPH<13> -ray emitters. Nonetheless, it is often challenging to establish associations between SNRs and VHE sources. For example, only specific regions of an SNR shell could be emitting or neighboring molecular clouds could be illuminated by multi-TeV particles that escaped the shock front of the SNR. Pulsar wind nebulae evolve as their pulsar ages and the available rotational energy (spin-down power) decreases. Since the X-ray synchrotron radiation from PWNe arises from higher energy electrons than the IC radiation in the VHE gamma-ray band, and the cooling time of the electrons decreases with their energy ( tc = E = ( dE = dt ), for ra-  \nTable 2. Results of the automatic association procedure for each catalog used (see main text for details and selections applied). The second column lists the numbers of objects in the HGPS survey region for each catalog. The third column gives the total number of associations found. The last column gives the number of HGPS sources having at least one associated object of a given category. The di GLYPH<11> erence between the two last columns is only large for 3FGL because 3FGL is the only counterpart catalog for which the source density is so high that many HGPS sources are associated with multiple 3FGL sources. Out of the 78 HGPS sources, only 11 are left without any association.  \ndiative losses tc / 1 = E ) we expect PWNe to shine longer in VHE gamma rays. Furthermore, a decreasing magnetic field with age can limit the emission time in radio and X-rays without a GLYPH<11> ecting the VHE emission. As a result, some old PWNe should be undetectable outside the VHE GLYPH<13> -ray domain (see, e.g., Aharonian et al. (1997); de Jager & Djannati-Ataï (2009); H.E.S.S. Collaboration et al. (2017f)). For such old PWNe only the detection of a middle-aged energetic pulsar in the vicinity of a VHE source can provide evidence toward the true nature of the VHE emission.  \nTo search for SNR and PWN associations, we take the most complete catalog of SNRs and PWNe to date into account, SNRcat 7 (Ferrand & Safi-Harb 2012). The SNRcat is a census of Galactic supernova remnants and their high-energy observations. It is based on the radio Catalogue of Galactic Supernova Remnants (Green 2014) but additionally includes other types of remnants in an e GLYPH<11> ort to be as complete and up-to-date as possible. In particular, it contains plerionic objects, PWNe with no observed shell. The possible presence of a PWN is usually assessed based on the presence of di GLYPH<11> use, nonthermal emission in radio, X-rays, or even GLYPH<13> -rays. Several of these cataloged objects have been classified by SNRcat as candidate PWNe solely because of the presence of VHE emission in the vicinity of an energetic pulsar. We removed those objects from the catalog used in our association procedure to avoid cases in which we might misleadingly selfassociate.  \nFor the association procedure, we split the SNRcat objects into three subsets based on their apparent type. The first subset consists of objects that have no evidence of nebular emission and mostly belong to the shell or filled-center types in SNRcat; this subset contains 211 objects within the HGPS region. The second subset consists of objects that are listed in SNRcat as PWNe (or PWNe candidates) showing no evidence for shell-like emission; this subset contains 29 objects within the HGPS region. The third subset consists of objects showing evidence of both shell and nebular emission, which we refer to as composite objects; this subset contains 42 objects within the HGPS region. For a fur-  \nther discussion of a potential PWN nature of these objects see the population study presented in H.E.S.S. Collaboration et al. (2017f).', 'Energetic pulsars': 'We selected energetic pulsars from version 1.54 of the ATNF catalog of radio pulsars (Manchester et al. 2005). We excluded millisecond pulsars because they are not expected to power VHE PWNe and applied a cut on the spin-down energy flux ˙ E = d 2 > 10 33 erg s GLYPH<0> 1 kpc GLYPH<0> 2 on the remaining pulsars. In addition, to take into account energetic pulsars of unknown distance, we included all objects with a spin-down luminosity ˙ E > 10 34 erg s GLYPH<0> 1 , resulting in a total of 222 pulsars used in the association procedure. We did not take into account pulsars that do not have a measured ˙ E . It is important to note that pulsars represent indirect associations: the associated pulsars are not directly emitting the unpulsed VHE GLYPH<13> -ray emission found in the HGPS, but rather indicate that they could be powering a PWN that directly emits such emission.', 'HE GLYPH<13> -ray sources': 'Of the 352 3FGL sources present in the HGPS region, we find 64 to be associated with an HGPS source. As expected, we also find a large portion of the 44 2FHL sources in the HGPS region to be associated with HGPS sources: only 13 of these have no HGPS counterpart. One of these sources is notably coincident with the VHE source candidate HOTS J1111 GLYPH<0> 611 (Sect. 5.6.17). Many of the other 2FHL sources lacking an HGPS association tend to be located in low-sensitivity parts of the HGPS region. Only four 2FHL sources in parts of the HGPS with good sensitivity show no significant VHE emission in the HGPS: Puppis A (H.E.S.S. Collaboration et al. 2015c), 2FHL J0826.1 GLYPH<0> 4500, GLYPH<17> Carinae, and the composite SNR G326.3 GLYPH<0> 1.8 (Temim et al. 2013).', 'Sources without physical associations': 'We find 24 of the 78 HGPS sources to be associated with shelllike SNRs. Given the large number of such objects in the HGPS region (211) and given their sizes, the number of chance coincidences is non-negligible. This is to be expected since we have not tried to specifically match SNR and HGPS source positions and sizes as in Acero et al. (2015b). Nonetheless, as discussed below, we find six known shells in the HGPS to be firmly identified and two more to be VHE shell candidates based on their specific morphologies (H.E.S.S. Collaboration et al. 2017a). We study the population of known SNRs in the HGPS further in a companion paper (H.E.S.S. Collaboration et al. 2017b).', 'Pulsar wind nebulae and composites': 'We find 37 of the SNRcat objects (in the HGPS region) containing a PWN or PWN candidate to be associated with an HGPS source. Conversely, we find more than 40% of HGPS sources to have at least one associated object in the PWN or composite classes. This supports the notion that systems containing PWNe are prolific VHE emitters. As discussed below, we are able to firmly identify about half of these associations using additional observational evidence such as similar MWL morphology or energy-dependent GLYPH<13> -ray morphology.', 'Pulsars': 'We find 47 of all the HGPS sources to be associated with an energetic pulsar. This suggests that the population of HGPS sources contains numerous PWNe. However, we selected a relatively low threshold ˙ E in our association criteria to minimize missed associations. We quantitatively study such selection e GLYPH<11> ects in a companion paper (H.E.S.S. Collaboration et al. 2017f) that provides a detailed look at the physical characteristics of firmly identified PWNe and a list of candidate PWN identifications based on various expected characteristics.', 'Extra associations': "For completeness, in addition to the associations obtained through the catalog-based, automatic procedure, we add a list of 20 extra associated objects that are plausible counterparts for some HGPS sources and are not covered by the limited set of catalogs we use. Previous publications had proposed most of these associations, often thanks to dedicated MWL observations of VHE source regions (e.g., the X-ray source XMMU J183245 GLYPH<0> 0921539 and HESS J1832 GLYPH<0> 093). We propose other associations in this work for some of the new sources (Sect. 5.6). We also include the original identifiers of VHE sources discovered first by other instruments (e.g., VER J1930 + 188, which corresponds to HESS J1930 + 188). Table 12 includes all of these extra associations, labeled 'EXTRA'.  \nEleven HGPS sources do not have any associations with known physical objects, although some are associated with HE GLYPH<13> -ray sources. We list and discuss these briefly here (the new VHE sources are discussed in Sect. 5.6):  \n- 1. HESS J1457 GLYPH<0> 593 is one of the new sources detected in the HGPS analysis. Although the automatic association procedure does not find any counterparts, the VHE GLYPH<13> -ray emission may originate in a molecular cloud illuminated by CRs that escaped from the nearby but o GLYPH<11> set SNR G318.2 + 0.1. This scenario is briefly described in Sect. 5.6.2.\n- 2. HESS J1503 GLYPH<0> 582 is also a new HGPS source and does not have any compelling association except for the HE GLYPH<13> -ray sources 3FGL J1503.5 GLYPH<0> 5801 and 2FHL J1505.1 GLYPH<0> 5808, neither of which is of a firmly identified nature. We describe this enigmatic source in Sect. 5.6.4.\n- 3. HESS J1626 GLYPH<0> 490 has only one association, with the HE GLYPH<13> -ray source 3FGL J1626.2 GLYPH<0> 4911. A dedicated XMM-Newton observation did not reveal any compelling X-ray counterpart either (Eger et al. 2011).\n- 4. HESS J1702 GLYPH<0> 420 is near the point-like source 2FHL J1703.4 GLYPH<0> 4145. The elongation of the VHE GLYPH<13> -ray emission prevented the automated procedure from making the association, but a connection between the objects seems plausible. The small size SNR G344.7 GLYPH<0> 0.1 (about 8 0 in diameter) is also in the vicinity and in good positional coincidence with the (point-like) 2FHL source.\n- 5. HESS J1708 GLYPH<0> 410 has no compelling association, even though this source was the target of dedicated X-ray observations to look for associated emission (Van Etten et al. 2009). Given the brightness and relatively steep spectrum of this VHE source ( GLYPH<0> = 2 : 57 GLYPH<6> 0 : 09), the absence of a counterpart at lower GLYPH<13> -ray energies in the Fermi -LAT catalogs is surprising and suggests the emission peaks in the hundreds of GeV range.\n- 6. HESS J1729 GLYPH<0> 345 is north of the nearby SNR HESS J1731 GLYPH<0> 347 (H.E.S.S. Collaboration et al. 2011a). An investigation into a potential connection between the two suggests the VHE emission from the former could be from a molecular cloud illuminated by hadronic particles that escaped from the SNR (Cui et al. 2016).\n- 7. HESS J1741 GLYPH<0> 302 is the subject of a dedicated companion paper (H.E.S.S. Collaboration et al. 2017c) discussing potential PWNe and SNR-related association scenarios, among others. These aspects are therefore not discussed here.\n- 8. HESS J1745 GLYPH<0> 303 is close to, but o GLYPH<11> set from, SNR G359.1 GLYPH<0> 0.5. Suzaku observations have revealed neutral iron  \nline emission in the region, suggesting the presence of molecular matter and making this object another possible case of a CR-illuminated cloud (Bamba et al. 2009). We find this object also to be associated with the HE GLYPH<13> -ray sources 2FHL J1745.1 GLYPH<0> 3035 and 3FGL J1745.1 GLYPH<0> 3011.  \n- 9. HESS J1828 GLYPH<0> 099 is a new HGPS source described in Sect. 5.6.8.\n- 10. HESS J1832 GLYPH<0> 085 is also a new HGPS source, described in Sect. 5.6.9.\n- 11. HESS J1858 + 020 has an association with the HE GLYPH<13> -ray source 3FGL J1857.9 + 0210 and is close to, but o GLYPH<11> set from, SNR G35.6 GLYPH<0> 0.4. A dedicated study (Paredes et al. 2014) did not find any compelling X-ray counterpart, although multiple possible scenarios were investigated, including CRilluminated molecular clouds.", '5.1.3. Firmly identified HGPS sources': 'In this section, we go one step further and treat those HGPS sources for which the physical origin of the VHE GLYPH<13> -ray emission has been firmly identified. Whereas the association criteria were principally based on positional evidence (angular o GLYPH<11> -set), we also perform a census of the additional evidence that is available to reinforce spatial associations and arrive at firm identifications. The supplementary observables we consider are correlated MWL variability, matching MWL morphology, and energy-dependent GLYPH<13> -ray morphology (Hinton & Hofmann 2009). Table 3 summarizes the results, along with the respective references for the additional evidence. Among the 78 sources in the HGPS region, we determine 31 to be firmly identified.  \nFirm identifications rely on di GLYPH<11> erent forms of evidence that vary depending on the source class. The VHE GLYPH<13> -ray emission from compact binary systems is always point-like and should exhibit variability that is also seen at lower energies. In contrast, the VHE emission from shell-type SNRs is extended (provided the SNR is su GLYPH<14> ciently large and close) and nonvariable, but can be identified based on the specific shell morphology and correlated morphology at lower energies.  \nComposite SNRs have both a shell and an interior PWN detected at lower energies and can be more complex to identify correctly. If the angular size of the shell emission is larger than the size of the VHE emission, we can identify the VHE emission as coming from the PWN filling the SNR. This is the case, for example, for HESS J1747 GLYPH<0> 281 (PWN in SNR G0.9 + 0.1) and HESS J1554 GLYPH<0> 550 (PWN in SNR G327.1 GLYPH<0> 1.1). In other cases, we are only able to identify the HGPS source with the composite SNR as a whole, i.e., we are confident that the VHE emission originates in the composite object but cannot disentangle whether it comes predominantly from the PWN or the shell (usually due to PSF limitations).  \nMore evolved stellar remnant systems are di GLYPH<14> cult to identify firmly. We can make a firm PWN identification when there is a PWN of comparable size and compatible position detected at lower energies. This is the case, for example, for HESS J1420 GLYPH<0> 607 (PWN G313.54 + 0.23) and HESS J1356 GLYPH<0> 645 (PWN G309.92 GLYPH<0> 2.51). In the absence of any clear PWN, or when its size at lower energies is much smaller than the VHE source, we have to rely on other evidence. The clearest such evidence is the detection of energy-dependent morphology, expected in PWNe because of the cooling of energetic electrons as they are transported away from the pulsar. At higher energies, the extent of the emission shrinks and its barycenter moves closer to the pulsar. This is the case for two sources thus far, HESS J1303 GLYPH<0> 631 (PWN G304.10 GLYPH<0> 0.24) and HESS J1825 GLYPH<0> 137 (PWN G18.00 GLYPH<0> 0.69). In the absence of such evidence, the identification of a VHE source as a PWN remains tentative when the only evidence is an energetic pulsar in the vicinity. Candidate PWN identifications are evaluated in detail in a companion paper (H.E.S.S. Collaboration et al. 2017f).  \nA large percentage (39%) of the 31 firmly identified sources are PWNe. The next largest source classes identified are SNR shells (26%) and composite SNRs (26%). Finally, GLYPH<13> -ray binary systems are also identified in the HGPS. It is not yet possible to identify firmly more than half of the total 78 HGPS sources with the conservative criteria we adopted, although the vast majority have one or more promising spatial associations that could prove to be real identifications following more in-depth studies beyond the scope of this work. We do not find any physical associations for 11 of the VHE sources in the HGPS, although for some of these, potentially related emission is seen in HE GLYPH<13> -rays, and for others, o GLYPH<11> set counterparts are present but simply not found by the automated association procedure adopted (see previous section). Figure 10 summarizes these identifications.  \nWe note that one source in HGPS, HESS J1943 + 213, is likely an extragalactic object. It has no measured extension and a radio counterpart that many recent studies tend to classify as a BL-Lac object (Peter et al. 2014; Straal et al. 2016; Akiyama et al. 2016). However, its VHE flux has not revealed any variability so far, which is unusual for such an object (Shahinyan & the VERITAS Collaboration 2016).', '5.2. Large-scale emission': "In Sect. 4.6, we introduced an empirical spatial model to account for the large-scale VHE GLYPH<13> -ray emission we observed along the Galactic plane to detect and characterize accurately the discrete VHE GLYPH<13> -ray sources. This model provides an estimate of the spatial distribution of the large-scale VHE emission discovered by Abramowski et al. (2014a). We find that the fit amplitude, latitudinal width, and position of this model, shown on Fig. 6, are consistent with the latitude profile of that previous work. The  \nTable 3. Table of 31 firmly-identified objects among the HGPS sources. The object classes are GLYPH<13> -ray binary, shell-type supernova remnant (SNR), pulsar wind nebula (PWN), and composite SNR (in cases where it is not possible to distinguish between the shell and interior nebula). The evidence used to identify the VHE GLYPH<13> -ray emission include position, morphology, variability, and energy-dependent morphology (ED Morph.).  \nFig. 10. Source identification summary pie chart. See Table 3 and Sect. 5.1.3.  \n<!-- image -->  \nwidth is also comparable to the HGPS source latitude distribution (Fig. 11, GLYPH<11> .) but smaller than that of molecular gas traced by CO emission (Dame et al. 2001).  \nOwing to the observational constraints and analysis used, the large-scale emission model cannot be considered a measurement of the total Galactic di GLYPH<11> use emission. The large-scale emission model provides an estimate of the di GLYPH<11> use emission present in the HGPS maps. Its parameter values depend on the map construction technique, in particular the exclusion region mask used in the analysis (Sect. 3.2.2), i.e., changes in the mask can alter the parameters of the model. For instance, the peak observed at ' GLYPH<24> 340 GLYPH<14> in Fig. 6 is due to the presence of low-level emission that is just below the threshold to be covered by the exclusion mask we use for the HGPS. While a significant percentage of the large-scale emission is expected to be truly interstellar diffuse emission, it is very likely that emission from discrete but unresolved sources contributes as well. Finally, some features in the HGPS large-scale emission model are likely artifacts of errors in the estimation of the background model of gamma-like cosmic ray EAS events (see Sect. 3.2); these events are the dominating model component in the HGPS counts maps, thus small  \nrelative errors in that background model can lead to significant changes in the excess model of the HGPS sources, but even more so the HGPS large-scale emission model.", '5.3. Source parameter distributions': "In the following section we study the global properties of the VHE GLYPH<13> -ray sources in the HGPS catalog. We compare certain key source parameters against each other and briefly discuss the implications in the context of the Galactic VHE source population, survey sensitivity, and firmly identified MWL source classes.  \nThe latitude distribution of the 78 HGPS sources is shown in Fig. 11. The distribution has a mean of b = GLYPH<0> 0 : 41 GLYPH<14> and a width of 0 : 87 GLYPH<14> . For visual comparison, the latitude distributions of the main classes of associated counterparts (Sect. 5.1) - SNRs, energetic pulsars, 3FGL sources, and 2FHL sources - are shown in this figure. Also shown for reference is an estimate of the matter density profile as traced by Planck measurements of CO(1-0) line emission (Planck Collaboration et al. 2016). It should be kept in mind throughout this section that the HGPS sensitivity is not uniform as a function of longitude or latitude (Sect. 3.4.4).  \nThe HGPS latitude distribution of sources correlates well with both potential counterparts and tracers of matter density. The distribution is somewhat skewed toward negative latitudes even though the HGPS sensitivity has a relatively wide and flat coverage in latitude. In Fig. 11, the sensitivity is illustrated by two curves showing regions of relatively good sensitivity (e.g., at ' = 0 GLYPH<14> ) and relatively poor sensitivity (e.g., at ' = 333 GLYPH<14> ). These curves demonstrate that the HGPS sensitivity coverage in latitude is, in general, much wider than the HGPS source distribution. Although there are local exceptions at some longitudes, the latitude coverage is generally flat in the range GLYPH<0> 2 : 0 GLYPH<14> < b < 1 : 5 GLYPH<14> , at various locations even in GLYPH<0> 2 : 5 GLYPH<14> < b < 2 : 5 GLYPH<14> . However, the counterpart catalogs are known to su GLYPH<11> er from various selection biases and the Galactic disk itself is known to not be perfectly symmetric as observed across the spectrum.  \nIn addition, one might still argue that, given the narrow range of latitudes observed with respect to surveys at other wavelengths, the HGPS sources may not be representative of the underlying distribution of VHE GLYPH<13> -ray sources. However, in light of the counterpart distributions, in particular the 2FHL sources, it can be reasonably assumed that the limited latitude coverage only has a weak e GLYPH<11> ect on the observed source population distribution.  \nThe longitude distribution of the 78 HGPS sources is shown in Fig. 12, together with the molecular interstellar matter column density profile as traced by CO(1-0) line emission (same as in the previous figure). The latter, measured by Planck (Planck Collaboration et al. 2016), has a uniform exposure (sensitivity) over the sky, unlike the HGPS, adding caveats to potential detailed cor-  \nrelations seen in this figure. We can nevertheless robustly conclude that there is a very general correlation in longitude between the number of HGPS sources and the molecular matter column density and that the HGPS sources are mostly found in the inner GLYPH<24> 60 GLYPH<14> of the Galaxy. Additionally, the spiral arm tangents as traced by CO (Vallée 2014) are shown in Fig. 12. An increased number of sources could be expected in the directions of the near spiral arm tangents (see Fig. 16). In the longitude distribution, a slight excess of sources in the direction of Scutum and between Norma and Crux-Centaurus can be observed. However, because of the limited sample size of 1-6 sources per bin, no significant increased source density in the direction of spiral arm tangents can be observed.  \nFor comparison, we also added distributions for the FermiLAT catalogs 3FGL and 2FHL to Fig. 12. While Fermi-LAT has a roughly uniform exposure, their sensitivity in the HGPS region is reduced in the inner Galaxy where di GLYPH<11> use emission is brighter, and also the source extraction is very di GLYPH<11> erent from the HGPS approach, so that a direct comparison is not possible. Finally we have chosen not to show the SNR and pulsar distributions in the Galactic longitude distribution at all because the coverage of those catalogs is not uniform.  \nWe compare the HGPS source integral fluxes ( E > 1 TeV) to source sizes in panel A of Fig. 13 and show the distributions of fluxes and sizes separately in panels B and C, respectively. In the flux-size figure, we plot the approximate flux sensitivity limit of the HGPS as a function of source size. One can see that the sensitivity worsens as the source size increases, as expressed by Eq. 28. The HGPS sources indeed generally follow this trend. From Fig. 13, we therefore conclude that the HGPS can be considered complete down to GLYPH<24> 10% Crab for sources < 0 : 7 GLYPH<14> . For smaller sources ( < 0 : 1 GLYPH<14> ), the HGPS achieves completeness at a few % Crab (see also Fig. 4).  \nWe show the distribution of HGPS source integral fluxes ( E > 1 TeV), which are calculated assuming a spectral index of GLYPH<0> = 2 : 3, in panel B of Fig. 13. At higher fluxes, we naturally expect the number of sources to decrease. At the lowest fluxes, we also expect the number to be small, because we reached the sensitivity limit of the HGPS.  \nAs can be seen in panel C of Fig. 13 and despite the modest H.E.S.S. PSF (0 : 08 GLYPH<14> ), the majority of sources are not compatible with being point-like but rather found to be significantly extended and as large as 1 GLYPH<14> . Owing to the methods used for background subtraction (see Sect. 3.2.3), the HGPS is not sensitive to sources with larger sizes.  \nThe firmly identified HGPS sources (Sect. 5.1) are highlighted in Fig. 13. It can be seen that all identified binary systems appear as point-like sources in the HGPS, as expected. The PWNe appear to have various angular sizes, in agreement with the diversity observed in the VHE PWN population (H.E.S.S. Collaboration et al. 2017f). Most identified SNRs are  \nFig. 11. Galactic latitude distribution of the HGPS sources (gray histogram). The bin size of this histogram is 0 : 3 GLYPH<14> . The HGPS point source sensitivity is shown (in units of % Crab) at two di GLYPH<11> erent longitudes of 0 GLYPH<14> and 333 GLYPH<14> . For comparison, the pulsar (PSR), supernova remnant (SNR), 3FGL and 2FHL source distributions in the HGPS longitude range are shown as overlaid curves, smoothed with Gaussians of width 0 : 15 GLYPH<14> . The dashed line shows Planck measurements of CO(1-0) line emission as an estimate for matter density in the Galaxy and similarly smoothed. All curves are normalized to the area of the histogram.  \n<!-- image -->  \nFig. 12. Galactic longitude distribution of the HGPS sources (gray histogram). The bin size of this histogram is 5 GLYPH<14> . For comparison, the 3FGL and 2FHL source distributions (smoothed with a Gaussian of width 5 GLYPH<14> ) and the Planck measurements of CO(1-0) line emission as an estimate for matter density in the Galaxy (smoothed with a Gaussian of width 2 : 5 GLYPH<14> ) are shown for the range in Galactic latitude b GLYPH<20> 5 GLYPH<14> and normalized to the area of the histogram. Spiral arm tangent locations shown are from Vallée (2014).  \n<!-- image -->  \nextended, likely owing to selection bias (smaller SNRs are difficult to identify, e.g., through shell-like morphology) and the H.E.S.S. PSF. The identified composite SNRs, on the other hand, are typically smaller, owing to the di GLYPH<14> culty in disentagling VHE emission from the SNR shell and interior PWN, similarly related to the H.E.S.S. PSF. In any case, it does not seem possible to identify the nature of the many unidentified sources solely on the basis of their sizes or a flux-size comparison.  \nFigure 14 shows the distribution of the HGPS source powerlaw (PL) spectral indices GLYPH<0> . For consistency, the PL model spec-  \ntral index is used for all sources, even those for which an exponential cuto GLYPH<11> power law (ECPL) fits better. The index distribution has a mean GLYPH<0> = 2 : 4 GLYPH<6> 0 : 3. This is compatible with the index ( GLYPH<0> = 2 : 3) adopted in the production of the HGPS flux maps (Sect. 3.4.2) and the HGPS PSF computation (Sect. 4.3). We note that individual source indices have typical statistical uncertainties of order GLYPH<6> 0 : 2 and a similar systematic uncertainty; HGPS data are often not su GLYPH<14> cient to precisely constrain the index because the energy range covered with good statistical precision is typically only about one decade (1 < GLYPH<24> E < GLYPH<24> 10 TeV). Finally, the  \nHESS J1825-137  \nFig. 13. Panel A: Integral source flux ( E > 1 TeV) vs. source size scatter plot with colors representing the di GLYPH<11> erent classes of firmly identified sources. For HGPS sources modeled as single Gaussians, the size is its width ( GLYPH<27> ). For sources modeled as multiple Gaussians (indicated with a circle around the marker), the size is the RMS of the two-dimensional intensity distribution (see Eq. 16). For sources with shell-like morphology (SNRs), the size is the outer shell radius. To improve the visibility of the plot, we do not show the SNR Vela Junior (HESS J0852 GLYPH<0> 463) at a size of 1 GLYPH<14> and a flux of 103% Crab. We illustrate the approximate sensitivity limit of the HGPS, as defined in Eq. 28, with an assumed point-source sensitivity of 1% Crab and an uncertainty band with a factor GLYPH<6> 2 to represent the sensitivity variations in the survey region (see caveats in main text). Panel B: Distribution of the integral fluxes ( E > 1 TeV)of the HGPS sources; colors are shown as in panel A. Panel C: Distribution of the HGPS source sizes; colors shown as in panel A. The first bin contains 30 sources, of which 17 are compatible with point-like sources according to Eq. 18. As in panel A, we omit Vela Junior, at a size of 1 GLYPH<14> .  \n<!-- image -->  \nfigure also shows how the firmly identified HGPS sources are distributed in index, showing no strong tendency with respect to source class.  \nWe show the cumulative log N ( > S ) - log S distribution of HGPS source integral fluxes ( E > 1 TeV, obtained from the maps) in Fig. 15. The 78 HGPS sources span a range in flux from 0.6% Crab to 103% Crab; 32 sources are above 10% Crab. We performed an unbinned likelihood fit of a PL model to the log N - log S distribution (also shown in Fig. 15), using only  \nthe range S > 10% Crab where we consider the HGPS survey mostly complete. The best-fit value of the PL slope is GLYPH<0> 1 : 3 GLYPH<6> 0 : 2 (for the cumulative distribution), and the amplitude corresponds to 32 GLYPH<6> 5 sources above 10% Crab. This slope is consistent with Galactic models in which equal-luminosity sources are homogeneously distributed in a thin disk, which predict a slope of GLYPH<0> 1 : 0. 8  \nFig. 14. Distribution of the HGPS source power-law (PL) spectral indices. For consistency, the PL model spectral index is used for all sources, even those for which an exponential cuto GLYPH<11> power law (ECPL) fits better. Taking statistical and systematic uncertainties into account, all indices are compatible within 2 GLYPH<27> with the mean GLYPH<0> = 2 : 4 GLYPH<6> 0 : 3 of the distribution.  \n<!-- image -->  \nFig. 15. Cumulative log N ( > S ) - log S distribution for the HGPS sources, showing the number of sources N above given flux thresholds S (integral flux above 1 TeV in % Crab). The line and error band show the result of an unbinned PL fit above a flux threshold of 10% Crab; the dashed line in the 1-10% Crab flux range illustrates the extension of the PL to fluxes below 10% Crab (for comparison, not fitted in that range).  \n<!-- image -->  \nThe only robust statement that can be inferred from the log N - log S distribution of HGPS sources is that it provides a lower limit on the true log N - log S distribution; that is, there are at least, for example, 70 sources above 1% Crab. If one assumes that log N - log S distributions are always concave (which most 'reasonable' spatial distributions and source luminosity functions encountered in the literature are), then the extrapolation of the PL fit shown in Fig. 15 sets an upper limit of GLYPH<24> 600 sources above 1% Crab, with a statistical error of a factor of 2.  \nMore detailed analyses of the log N - log S distribution or of the flux-size distribution are possible in principle but in practice do not yield robust results because of the limited number of sources and the large uncertainties concerning the e GLYPH<11> ective sensitivity achieved. We emphasize that the catalog creation procedure is complex (special treatment of known shell-type sources, large-scale emission model component, 15 discarded and several merged components; see Sect. 4.9), with the net e GLYPH<11> ect that the sensitivities shown in Fig. 4 and panel A of Fig. 13 are not reliably achieved, because those sensitivity estimates assume isolated sources, there is no underlying large-scale emission or source confusion, and there is a detection threshold of 5 GLYPH<27> , whereas the component detection threshold of TS = 30 corresponds to GLYPH<24> 5 : 5 GLYPH<27> .  \nA representation of the Galaxy seen face-on is depicted in Fig. 16 to visualize how much of the Galaxy the HGPS has been able to probe at di GLYPH<11> erent sensitivity levels. Two limits are shown, illustrating the sensitivity detection limit (horizon) of the HGPS for potential point-like sources with presumed luminosity of 10 33 and 10 34 erg / s. Given the achieved sensitivity in the Galactic plane, it is clear that H.E.S.S. has only probed a small fraction of the Galaxy - just up to a median distance of 7 : 3 kpc for bright (10 34 erg / s) point-like sources (and less for extended sources). Furthermore, this illustrative look at survey completeness strengthens the hypothesis that the large-scale emission described in Sect. 4.6 could be partly explained by a population of unresolved sources, presumed to be distant.", '5.4. Comparison with previous VHE publications': 'In total, we reanalyzed 48 VHE GLYPH<13> -ray sources that have been the subject of past H.E.S.S. publications. In this section we present a systematic comparison of the present HGPS results, with the latest published results, as summarized in gamma-cat 9 , the open TeV source catalog.  \nWe associated HGPS sources with previous analyses simply by the name of the source, which was unique except for three cases: HESS J1800 GLYPH<0> 240, HESS J1746 GLYPH<0> 308, and HESS J1930 + 188, which we discuss in detail in Sect. 5.4.2. We excluded these sources from the systematic comparison in the first place.  \nTo further identify the cases for which we obtained significantly di GLYPH<11> erent results from previously published analyses, we compared the position, size, spectral index, and flux of the remaining uniquely associated sources, taking statistical and systematic errors of the measurements into account. For each of these parameters, we estimated the total uncertainty GLYPH<27> tot as the 1 GLYPH<27> statistical and systematic uncertainties added in quadrature.  \nFig. 16. Illustration of the location of identified H.E.S.S. sources in the Galaxy with respect to HGPS completeness (sensitivity limits). This is a face-on view; the spiral arms (Vallée 2014) are schematically drawn as gray bars. The HGPS horizons for source luminosities of 10 33 and 10 34 erg / s (for a putative 5 GLYPH<27> detection of a point-like source, same as Fig. 4) are depicted by light blue and light brown lines (and shaded regions therein), respectively. The source distances are from SNRcat (Ferrand & Safi-Harb 2012) and ATNF pulsar catalog (Manchester et al. 2005). When no distance uncertainties were available, we applied a generic uncertainty of factor two on the distance. The three labeled sources are the Galactic GLYPH<13> -ray sources outside the HGPS region detected by H.E.S.S.  \n<!-- image -->  \nWe estimated this quantity for both the HGPS-derived source parameters and previously published H.E.S.S. values.  \nThe systematic uncertainties on position and size are given in Sect. 4.10.3 and Sect. 4.10.2, respectively. Additionally, we assumed a systematic uncertainty GLYPH<1>GLYPH<0> syst = 0 : 2 on the spectral index and 30% on the flux of the source, in agreement with previous estimates (Aharonian et al. 2006b). We then defined the  \ncriterion for significant outliers as  \nGLYPH<1> HGPS GLYPH<0> H : E : S : S : > 2 q GLYPH<27> 2 tot ; HGPS + GLYPH<27> 2 tot ; H : E : S : S : ; (29)  \nwhere GLYPH<1> HGPS GLYPH<0> H : E : S : S : is the di GLYPH<11> erence between the corresponding parameter values. When comparing the position we chose the angular separation as comparison parameter. We note that for many  \nsources, the data sample used here is significantly di GLYPH<11> erent from that used in the publication, hence the correlation of statistical errors is usually not too large.  \nWe first discuss the general level of agreement between the current and previous analyses (excluding the outliers) in Sect. 5.4.1 and later discuss the outliers of the comparison individually in Sect. 5.4.2.', '5.4.1. Agreement with previous publications': 'For the vast majority of sources, we find that there is good agreement between the HGPS-derived position, morphology, and spectrum within the statistical and systematic uncertainties.', 'Spectral index': 'For all sources we found the spectral power-law indices to be compatible with the previously published values. The mean difference in spectral index was 0 : 04 with a width of 0 : 23 for the distribution. This is well compatible with the expected scatter taking statistical and systematic uncertainties of the measured spectral indices into account.', '5.4.2. Differences with previous publications': 'In the following paragraphs, we list and discuss the outliers as identified by Eq.29.', 'HESS J1729 GLYPH<0> 345': 'This source (Vela X) exhibits complex morphology, and the HGPS analysis best models the VHE emission as a superposition of three Gaussian components with an average size 0 : 58 GLYPH<14> GLYPH<6> 0 : 052 GLYPH<14> . This value is somewhat larger than the value published first in Aharonian et al. (2006c), where it was modeled as a single asymmetric Gaussian of size 0 : 48 GLYPH<14> GLYPH<6> 0 : 03 GLYPH<14> GLYPH<2> 0 : 36 GLYPH<14> GLYPH<6> 0 : 03 GLYPH<14> . However, a more recent H.E.S.S. publication (Abramowski et al. 2012a) studied the complex emission more thoroughly. It fit profiles of the emission along two perpendicular axes, the main one aligned with the primary orientation of the emission. Along the major axis, the study measured a Gaussian size 0 : 52 GLYPH<14> GLYPH<6> 0 : 02 GLYPH<14> , and along the minor axis, two Gaussians (sizes 0 : 12 GLYPH<14> GLYPH<6> 0 : 02 GLYPH<14> and 0 : 60 GLYPH<14> GLYPH<6> 0 : 04 GLYPH<14> ) were required to best fit the emission. The HGPS model of the emission from HESS J0835 GLYPH<0> 455 is thus largely compatible with the most recent dedicated study of the VHE emission, and the apparent discrepancy is simply a result of comparing two di GLYPH<11> erent multi-component models with our general outlier criterion (Eq. 29).', 'HESS J1646 GLYPH<0> 458': 'HESS J1646 GLYPH<0> 458 is a complex emission region located in the vicinity of the stellar cluster Westerlund 1. Its morphology suggests it consists of multiple sources. Abramowski et al. (2012b) separated the emission into at least two distinct features (with radii 0 : 35 GLYPH<14> and 0 : 25 GLYPH<14> , respectively) as well as some structured extended emission, distributed over the signal region of 2.2 GLYPH<14> diameter, and even extending beyond. A flux above 1 TeV in the signal region of 7 : 6 GLYPH<6> 1 : 3 GLYPH<6> 1 : 5 GLYPH<2> 10 GLYPH<0> 12 cm GLYPH<0> 2 s GLYPH<0> 1 was derived, and a spectral index of 2 : 19 GLYPH<6> 0 : 08 GLYPH<6> 0 : 20. An ON-OFF background estimation technique was used to cope with the large source size.  \nIn the HGPS analysis, this complex emission is modeled by a single Gaussian component of 0.5 GLYPH<14> size shifted by 0.47 GLYPH<14> from the center of the region used in Abramowski et al. (2012b), with a lower flux above 1 TeV of 5 : 48 GLYPH<6> 0 : 46 GLYPH<2> 10 GLYPH<0> 12 cm GLYPH<0> 2 s GLYPH<0> 1 , and steeper index of 2 : 54 GLYPH<6> 0 : 13. Given the complex morpology and the large scale of the spectral extraction region used in Abramowski et al. (2012b), significant di GLYPH<11> erences in source parameters are to be expected; in the HGPS analysis part of the flux is absorbed in the large-scale di GLYPH<11> use background.', 'HESS J1708 GLYPH<0> 410': 'The flux above 1 TeV of HESS J1708 GLYPH<0> 410 is found to be smaller in the HGPS analysis than in Aharonian et al. (2008a). While the size of the source is similar in both cases, the di GLYPH<11> erent approaches used in the HGPS analysis lead to di GLYPH<11> erent integration radii used to derive the source spectrum. The HGPS analysis uses an integration radius about two times smaller than in the dedicated analysis, which explains the apparent discrepancy.  \nFor HESS J1729 GLYPH<0> 345, the HGPS analysis finds a flux above 1 TeV larger than in H.E.S.S. Collaboration et al. (2011a). Because of the HGPS morphology modeling of the source and its procedure to define the integration radius, the spectrum of this source is derived in a region with a radius about two times larger than in the dedicated publication, accounting for the observed di GLYPH<11> erence.', 'HESS J1745 GLYPH<0> 303': 'HESS J1745 GLYPH<0> 303 was studied in Aharonian et al. (2008b) with 80 h of data. Its morphology is complex and three subregions, called A, B, and C, were discussed. In the HGPS analysis, with more than 160 h on the region, two distinct sources are detected: HESS J1745 GLYPH<0> 303 and HESS J1746 GLYPH<0> 308. The former encloses the hotspots A and C and a fraction of region B. A second source is now detected at b = GLYPH<0> 1 : 11 GLYPH<14> latitude. This source contains part of hotspot B and emission at large latitudes that was not significant before, likely due to the additional livetime obtained since 2008. It is fainter and its spectrum is very steep but poorly constrained. There is also a third extended ( GLYPH<27> GLYPH<24> 0 : 5 GLYPH<14> ) Gaussian component in the region. It is currently considered to be a di GLYPH<11> use component. The association of the two sources and the extended component is unclear and the exact morphology of the VHE emission in the region will require dedicated studies.', 'HESS J1800 GLYPH<0> 240': 'In Aharonian et al. (2008d) the emission in the region of W 28 was found to be split into two components: HESS J1801 GLYPH<0> 233 (addressed below), which is not significant in the HGPS analysis and is coincident with the W 28 SNR itself, and a complex region HESS J1800 GLYPH<0> 240 o GLYPH<11> set by 0 : 5 GLYPH<14> to the south. The latter was previously found to be resolved into three hotspots dubbed HESS J1800 GLYPH<0> 240 A, B, and C (Aharonian et al. 2008d). Since sources HESS J1800 GLYPH<0> 240 A and B are spatially coincident with molecular clouds, Aharonian et al. (2008d) suggested that they were produced by CRs that had escaped the SNR and had illuminated ambient gas clouds, making this system an archetype of CR escape from evolved SNRs (see, e.g., Aharonian & Atoyan 1996; Slane et al. 2015; Gabici & Montmerle 2015).  \nIn the HGPS analysis, however, only one source is redetected, HESS J1800 GLYPH<0> 240, as one large Gaussian component centered on the hotspot B. The separation into several components does not result in a high enough TS to separate it into several significant sources in the analysis shown here.', 'HESS J1825 GLYPH<0> 137': 'HESS J1825 GLYPH<0> 137 is a large PWN with a bright core surrounded by extended, asymmetric emission. The HGPS analysis finds it has a size of 0 : 46 GLYPH<14> GLYPH<6> 0 : 03 GLYPH<14> , using three Gaussian components to  \nmodel the VHE entire GLYPH<13> -ray emission. This is significantly larger than the 0 : 24 GLYPH<14> GLYPH<6> 0 : 02 GLYPH<14> obtained with a single symmetric Gaussian model or the 0 : 23 GLYPH<14> GLYPH<6> 0 : 02 GLYPH<14> GLYPH<2> 0 : 26 GLYPH<14> GLYPH<6> 0 : 02 GLYPH<14> with a single asymmetric Gaussian in Aharonian et al. (2006g). These models were stated to have rather poor GLYPH<31> 2 goodness-of-fit values. The more complex approach taken for the morphology modeling in the HGPS improves the description of the GLYPH<13> -ray emission from this PWN and accounts for the di GLYPH<11> erences with respect to previous, simpler modeling.', 'HESS J1837 GLYPH<0> 069': 'The HGPS analysis finds HESS J1837 GLYPH<0> 069 to have a size of 0 : 36 GLYPH<14> GLYPH<6> 0 : 03 GLYPH<14> based on modeling the VHE GLYPH<13> -ray emission as three Gaussian components. This is larger than the size previously derived using a single asymmetric Gaussian (Aharonian et al. 2006d), i.e., 0 : 12 GLYPH<14> by 0 : 05 GLYPH<14> ; and using a single Gaussian (Marandon et al. 2008), i.e., 0 : 22 GLYPH<14> . The more complex modeling of the HGPS, which also takes into account more of the extended nebular emission from this identified PWN, explains the apparent discrepancy. Consequently, we used a larger region (twice the radius compared to Aharonian et al. (2006d)) to derive the spectrum, leading to an integral flux above 1 TeV that is larger by a factor of GLYPH<24> 3 than in the dedicated publication.', 'HESS J1857 + 026': 'The size of the source HESS J1857 + 026 is significantly larger in this analysis than previously published in Aharonian et al. (2008a). In the latter, the source is fit with an asymmetric Gaussian (0 : 11 GLYPH<14> GLYPH<6> 0 : 08 GLYPH<14> GLYPH<2> 0 : 08 GLYPH<14> GLYPH<6> 0 : 03 GLYPH<14> ), whereas the HGPS analysis best models the source with two Gaussian components for an approximate size of 0 : 26 GLYPH<14> GLYPH<6> 0 : 06 GLYPH<14> . The di GLYPH<11> erence in size is explained by the multicomponent approach of the HGPS that better takes into account the larger scale emission underneath the central bright core.', 'HESS J1908 + 063': "The position and size published in Aharonian et al. (2009) are significantly di GLYPH<11> erent from those obtained in the HGPS analysis. The position is o GLYPH<11> set by 0 : 17 GLYPH<14> and the size is found to be 0 : 48 GLYPH<14> GLYPH<6> 0 : 04 GLYPH<14> , which is 0 : 14 GLYPH<14> larger. We note that the size we find is consistent with that measured by the VERITAS Collaboration (Aliu et al. 2014b), even though the positions di GLYPH<11> er by 0 : 3 GLYPH<14> . A plausible cause for these discrepancies is that this is a large source likely composed of multiple components, where results are expected to be sensitive to the morphology assumptions and to details in background modeling techniques, in particular, if those tend to absorb large-scale features.  \nHESS J1923 + 141  \nThe VHE GLYPH<13> -ray source HESS J1923 + 141 (preliminary H.E.S.S.  \nresults published in Fiasson et al. 2009) is spatially coincident with the W 51 region, studied in detail with the MAGIC IACT (Aleksi'c et al. 2012). The HGPS results are generally compatible with those from MAGIC. However, the latter shows evidence for a GLYPH<13> -ray source composed of two components above 1 TeV, which cannot yet be confirmed by H.E.S.S. One component is coincident with the interaction region between W 51C and W 51B, while the other is coincident with the potential PWN CXOU J192318.5 + 140305 (Koo et al. 2005), suggesting that HESS J1923 + 141 may be a composite of VHE emission of different astrophysical origins.", 'HESS J1930 + 188': 'The VHE GLYPH<13> -ray source, discovered with VERITAS (with the identifier VER J1930 + 188, Acciari et al. 2010), is coincident with the composite SNR G54.1 + 0.3 and the pulsar PSR J1930 + 1852. We report on the H.E.S.S. observations of this source for the first time here. The HGPS source is found to have a slightly displaced position from the pulsar and the VERITAS best fit (by 0 : 04 GLYPH<14> ). Despite the agreement with the VERITAS spectral index, the integral flux above 1 TeV found in our analysis is GLYPH<24> 40% lower than their published flux. We note, however, that the apparent discrepancy with VERITAS is not confirmed by our cross-check analysis, which yields a flux for this source that is larger by more than the nominal 30% systematic flux uncertainty, and is in agreement with the VERITAS measurement.', '5.4.3. Sources not redetected': 'In total, there are four previously published VHE GLYPH<13> -ray sources that are not redetected with the current HGPS analysis. All of these are rather faint sources which, for the HGPS analysis, yield significances close to the HGPS detection threshold of TS = 30. We consider these as real sources of GLYPH<13> -rays; the nondetection in the HGPS is primarily a result of di GLYPH<11> erences between the HGPS analysis and specific analysis methods. We found that some of the most relevant di GLYPH<11> erences are  \n- 1. event reconstruction and GLYPH<13> -ray-hadron separation cuts that are less sensitive compared to more specialized methods that have been used in individual source analyses;\n- 2. higher energy threshold in the HGPS analysis, in conjunction with a soft spectrum of the tested source;\n- 3. use of the 2 GLYPH<14> FoV o GLYPH<11> set cut (see Sect. 3.1), which is tighter than the value used in many previous H.E.S.S. publications (2 : 5 GLYPH<14> or even 3 GLYPH<14> ).  \nIn addition, the use of a large-scale emission model and the modeling of nearby extended components and overlapping sources modifies the measured flux and hence the significance of a source compared to previous analyses, where larger scale  \nbackground features were accounted for in di GLYPH<11> erent ways (e.g., partly absorbed in the ring background). Given these di GLYPH<11> erences, it is not surprising that few faint sources fail the HGPS detection criteria.  \nIn the following paragraphs, we describe the individual cases in more detail. For completeness, we added all missing sources to the final HGPS catalog; the source parameters were taken from the corresponding publication (see also Table 1).', 'HESS J1718 GLYPH<0> 374 and HESS J1911 + 090': 'The VHE GLYPH<13> -ray sources HESS J1718 GLYPH<0> 374 and HESS J1911 + 090 (Figs. 36 and 35) were previously detected toward the SNR G349.7 + 0.2 and W 49B / SNR G43.3 GLYPH<0> 0.2, respectively. Both are thought to result from interactions with molecular clouds and exhibit correspondingly steep (soft) spectra, which have PL indices GLYPH<0> = 2 : 80 GLYPH<6> 0 : 27 (H.E.S.S. Collaboration et al. 2015b) and 3 : 14 GLYPH<6> 0 : 24 (H.E.S.S. Collaboration et al. 2016b), respectively. The energy threshold of the analyses is therefore key to detecting these sources. As described in Sect. 3, the maps that serve as a starting point for the source catalog have been produced using the hard cuts configuration and a conservative safe energy threshold, explaining the lack of detection of these sources in the HGPS analysis.', 'HESS J1741 GLYPH<0> 302': 'The unidentified source HESS J1741 GLYPH<0> 302 is located on the Galactic plane ( b = 0 : 05 GLYPH<14> ) and GLYPH<24> 1.7 GLYPH<14> away from the Galactic center. With an integral flux of GLYPH<24> 1% Crab above 1 TeV it is one of the faintest H.E.S.S. sources detected so far (H.E.S.S. Collaboration et al. 2017c). Because of the addition of the large-scale emission model in the HGPS analysis, HESS J1741 GLYPH<0> 302 does not reach the HGPS TS = 30 detection threshold.', 'HESS J1801 GLYPH<0> 233': 'HESS J1801 GLYPH<0> 233 is part of the HESS J1800 GLYPH<0> 240 and HESS J1801 GLYPH<0> 233 source complex discussed above, characterizing emission features of the SNR W 28 region (Aharonian et al. 2008d). The emission was found to be split into two components: HESS J1801 GLYPH<0> 233, which is coincident with the northeastern boundary of W 28 where the shockwave is interacting with a molecular cloud, and a complex region HESS J1800 GLYPH<0> 240 o GLYPH<11> -set by 0 : 5 GLYPH<14> to the south. HESS J1801 GLYPH<0> 233 does not reach the TS = 30 threshold and is therefore not found to be significant in the HGPS analysis. We note that the GLYPH<13> -ray emission from W 28 is bright in the GeV range and is clearly detected above 50 GeV (Ackermann et al. 2016). It has a steep spectral index of 2 : 7 GLYPH<6> 0 : 3 at VHE (Aharonian et al. 2008d). It is therefore not detected here because of our higher analysis energy threshold (about 400 GeV at a longitude of 7 GLYPH<14> , see Fig. 2) and because of the inclusion of the large-scale emission model in our analysis, which reduces the  \nsignificance of such a faint source. Furthermore, we reiterate that HESS J1800 GLYPH<0> 240 is detected in the HGPS as one large Gaussian source, see Sect. 5.4.2, rather than three individual hotspots as in Aharonian et al. (2008d). This potentially also contributes to a reduction of the significance of this previously established source HESS J1801 GLYPH<0> 233.', '5.5. Comparison with the cross-check analysis': 'For most sources, the spectral fit results reported in this catalog agree with those obtained from the independent cross-check analysis (see Sect. 4.12). For the following sources, however, larger di GLYPH<11> erences, exceeding the systematic errors, are observed. Several factors could explain these di GLYPH<11> erences, such as the lower energy threshold in the cross-check analysis, the di GLYPH<11> erences in the morphology models, or the fact that the cross-check spectrum analysis is run for the positions and sizes obtained with the main analysis.  \n- -HESS J1503 GLYPH<0> 582 (see Sect. 5.6.4): While the spectral indices are compatible, the derived integral flux above 1 TeV is about two times higher in the main analysis than in the cross-check analysis.\n- -HESS J1646 GLYPH<0> 458 (Westerlund 1): the cross-check analysis gives a spectrum that is about two times brighter around 1 TeV, with a curvature or cuto GLYPH<11> leading to similar fluxes as the main analysis at low and high energies. We would like to stress again that the HGPS analysis for this source is not very reliable, because the source size is similar to the H.E.S.S. field of view and a more careful individual study and background estimation is needed, as explained in Sect. 5.4.2 which points out the di GLYPH<11> erences with the previously published measurement.\n- -HESS J1718 GLYPH<0> 385: For both the main and cross-check analyses, the preferred spectral model is a power law with an exponential cuto GLYPH<11> . The cuto GLYPH<11> energies are compatible and the spectra are in agreement above GLYPH<24> 3 TeV. However, below this energy, some discrepancy is observed as the main analysis spectral fit yields a spectral index that is harder than in the cross-check analysis, resulting in an integral flux above 1 TeV about two times lower in the main analysis than in the cross-check analysis.\n- -HESS J1729 GLYPH<0> 345: While the derived spectral indices are compatible, the integral flux above 1 TeV is about two times higher in the cross-check analysis than in the main analysis.\n- -HESS J1746 GLYPH<0> 308: The large spectral index derived from the main analysis could not be confirmed by the cross-check analysis. The di GLYPH<11> erential flux values at 1 TeV are compatible, but the discrepancy in the obtained spectral indices leads  \nto an integral flux above 1 TeV about two times higher in the cross-check analysis than in the main analysis.  \n- -HESS J1852 GLYPH<0> 000 (see Sect. 5.6.16): The derived spectral indices are compatible, but the integral flux above 1 TeV is about two times higher in the cross-check analysis than in the main analysis.  \nSpectral model results for these six sources should therefore be treated with caution.', '5.6. New VHE sources': 'During the construction of the HGPS catalog, statistically significant VHE GLYPH<13> -ray emission was detected from 16 sources which were not previously known or for which only preliminary detections had been published (e.g. in conference proceedings). All of these new sources are confirmed by the cross-check analysis -we do not expect any of these new sources to be a false detection (see Sect. 4.8 and 4.9) The morphological and spectral properties of these new, confirmed VHE sources are provided in Tables 10 and 11, their spectra are shown in Fig. 17. Each new source is also briefly described in the following sections, in the context of its MWL environment and possible origin of the VHE GLYPH<13> -rays.', '5.6.1. HESS J1119 GLYPH<0> 614': "We confirm the discovery of VHE GLYPH<13> -ray emission from HESS J1119 GLYPH<0> 614 (Fig. 18) and identify it as the composite SNR G292.2 GLYPH<0> 0.5. We base the firm identification on the basis of spatial coincidence with the SNR and its associated PWN G292.15 GLYPH<0> 0.54 and highly magnetized pulsar PSR J1119 GLYPH<0> 6127. H.E.S.S. previously published (Djannati-Ataï et al. 2009) preliminary source properties that are compatible with the HGPS results.  \nA compact (size 6 00 GLYPH<2> 15 00 ), nonthermal PWN has been detected in X-rays (Gonzalez & Safi-Harb 2003; Safi-Harb & Kumar 2008) and is considered a candidate PWN in HE GLYPH<13> -rays (Acero et al. 2013). It is powered by the energetic pulsar PSR J1119 GLYPH<0> 6127, with spin-down luminosity ˙ E = 2 : 3 GLYPH<2> 10 36 erg s GLYPH<0> 1 and distance d = 8 : 4 GLYPH<6> 0 : 4 kpc (Caswell et al. 2004). The pulsar has been detected in radio (Camilo et al. 2000) and HE GLYPH<13> -rays (Parent et al. 2011; Acero et al. 2015a, as 3FGL J1119.1 GLYPH<0> 6127 in the latter) and is characterized by a relatively high surface Bfield (4 : 1 GLYPH<2> 10 13 G). Despite it being a rotation-powered pulsar, it has recently joined the other high-B pulsar PSR J1846 GLYPH<0> 0258 in revealing a magnetar-like behavior (Gö˘gü¸s et al. 2016; Younes et al. 2016; Antonopoulou et al. 2016). It is further notable for being among the handful of pulsars for which braking indices have been measured, in this case n = 2 : 684 GLYPH<6> 0 : 002 (Weltevrede et al. 2011), as opposed to simply assuming n = 3, giving a more precise characteristic age GLYPH<28> c = P ( n GLYPH<0> 1) ˙ P = 1 : 9 kyr, where P and ˙ P  \nare the currently measured period and period derivative, respectively.  \nConsidering the luminosity of HESS J1119 GLYPH<0> 614, L GLYPH<13> (1 GLYPH<0> 10 TeV) = 2 : 4 GLYPH<2> 10 34 ( d = 8 : 4kpc) 2 erg s GLYPH<0> 1 , the apparent e GLYPH<14> ciency of converting the pulsar's rotational energy to GLYPH<13> -rays, GLYPH<15> 1 GLYPH<0> 10 TeV GLYPH<17> L GLYPH<13> = ˙ E = 1 : 1%, is compatible with the e GLYPH<14> ciencies ( < GLYPH<24> 10%) of other VHE sources that have been identified as PWNe (Kargaltsev et al. 2013). The o GLYPH<11> set of the VHE emission from this young pulsar, where the X-ray PWN is located, is not statistically significant with respect to the uncertainty on the best-fit VHE centroid ( GLYPH<6> 0 : 02 GLYPH<14> ).  \nThe age of SNR G292.2 GLYPH<0> 0.5 is in the range 4.2 GLYPH<0> 7.1 kyr (Kumar et al. 2012). This can be reconciled with the characteristic age of the pulsar if the braking index n was much smaller than the current value until recently. This assumption is reasonable in light of recent evidence for erratic radio timing behavior from the pulsar (Weltevrede et al. 2011). The X-ray emission from the SNR is predominantly thermal and has an additional hard, nonthermal, X-ray component. This nonthermal emission is likely from the PWN, although an origin in the SNR reverse shock could not be ruled out (Kumar et al. 2012).  \nThe X-ray spectral measurements suggest the SNR is generally expanding in a low-density medium, appearing to disfavor a hadronic origin for the VHE GLYPH<13> -rays (Drury et al. 1994). However, there is also evidence for localized, high-density regions near the eastern SNR shell, including dark clouds and CO features (Kumar et al. 2012). We cannot confirm the claim by Kumar et al. (2012), based on preliminary H.E.S.S. results (DjannatiAtaï et al. 2009), that no VHE emission is detected from the eastern SNR shell, as it is well within the VHE emission region in the HGPS analysis.  \nIn conclusion, while the identification with the composite SNR and PWN system is firm, it is not yet clear whether the VHE emission originates in the SNR shock, either leptonically, from the shell itself, or hadronically, from interactions with ambient media; the PWN; or some combination thereof.", '5.6.2. HESS J1457 GLYPH<0> 593': "VHE GLYPH<13> -ray emission from the new source HESS J1457 GLYPH<0> 593 (Fig. 19) is associated with the SNR G318.2 + 0.1, on the basis of a spatial coincidence with a shell-type SNR and lack of other potential MWL counterparts. Preliminary H.E.S.S. morphological properties were initially published by Hofverberg et al. (2010). The HGPS source position is compatible with the preliminary position; however, the size of the source in the catalog is different because of a di GLYPH<11> erence in the assumed morphological model. Previously, the source was modeled as an asymmetric Gaussian (0 : 31 GLYPH<14> GLYPH<6> 0 : 07 GLYPH<14> by 0 : 17 GLYPH<14> GLYPH<6> 0 : 05 GLYPH<14> ) whereas the HGPS source is modeled, like all HGPS sources, as a symmetric Gaussian (0 : 33 GLYPH<14> GLYPH<6> 0 : 04 GLYPH<14> ). Nonetheless, the spatial overlap between  \nFig. 17. Fitted power-law spectral models with uncertainty bands and flux points for new sources.  \n<!-- image -->  \nHESS J1457 GLYPH<0> 593 and the southern part of the SNR shell still holds.  \nG318.2 + 0.1 is observed as a relatively large (40 0 GLYPH<2> 35 0 ) shell in radio (e.g., Whiteoak & Green 1996), which is characterized by two arc-like, nonthermal filaments in the northwest and southeast (SE) that together form the shell. The VHE emission is much larger than the SNR shell, and the VHE centroid is significantly o GLYPH<11> set ( GLYPH<24> 0.4 GLYPH<14> ) from the SNR center, although it is partially coincident with the SE rim of the shell. Furthermore, there is evidence in 12 CO (Dame et al. 2001) of a giant molecular  \ncloud (GMC) at ( '; b ) GLYPH<25> (318 : 4 GLYPH<14> ; GLYPH<0> 0 : 5 GLYPH<14> ) coincident with both the VHE emission and the SE rim; this GMC is 1 : 8 GLYPH<14> GLYPH<2> 1 : 1 GLYPH<14> (average physical size 80 pc) in size and has mass GLYPH<24> 3 GLYPH<2> 10 5 M GLYPH<12> and density GLYPH<24> 40 cm GLYPH<0> 3 , assuming the near solution of the kinematic distance 3 : 5 GLYPH<6> 0 : 2 kpc (Hofverberg et al. 2010). Little is known about G318.2 + 0.1 itself, but assuming it is at the same distance as the GMC and further assuming a Sedov-Taylor model for the SNR evolution, its physical diameter would be GLYPH<24> 40 pc and its age GLYPH<24> 8 kyr. These data suggest a plausible SNR and molecular cloud interaction scenario (e.g., Gabici et al. 2007), where  \n<!-- image -->  \np  \nFig. 18. Significance ( GLYPH<25> TS ) of the VHE GLYPH<13> -ray excess, centered on the new source HESS J1119 GLYPH<0> 614, with the H.E.S.S. The PSF for this data set shown inset. The black circle at the center indicates the 68% uncertainty in the best-fit centroid of the VHE emission. The white circle represents the 70% containment region of the emission ( R\\_SPEC , used also for spectral measurement). The approximate size of the radio shell of SNR G292.2 GLYPH<0> 0.5 is shown as a green circle and the PWN G292.15 GLYPH<0> 0.54 as a green marker. The position of the pulsar PSR J1119 GLYPH<0> 6127 is denoted by a cyan diamond. The FoV is 1 : 5 GLYPH<14> GLYPH<2> 1 : 5 GLYPH<14> .Fig. 19. VHE GLYPH<13> -ray image: HESS J1457 GLYPH<0> 593. See Fig. 18 for a general description. Additionally, the SNR G318.2 + 0.1 is shown by plotting its 843-MHz radio intensity (Whiteoak & Green 1996) with contours at 4, 8, and 12 mJy / beam. The FoV is 2 : 8 GLYPH<14> GLYPH<2> 2 : 8 GLYPH<14> .  \n<!-- image -->  \nFig. 20. VHE GLYPH<13> -ray image: HESS J1458 GLYPH<0> 608. See Fig. 18 for a general description. Additionally, the ellipse represents the 95% uncertainty in the position of the HE GLYPH<13> -ray point source 3FGL J1456.7 GLYPH<0> 6046, and the cyan diamond indicates the position of the pulsar. The FoV is 2 : 6 GLYPH<14> GLYPH<2> 2 : 6 GLYPH<14> .  \n<!-- image -->  \nparticles are accelerated in the shell, escape, and interact with a nearby but o GLYPH<11> set MC, producing GLYPH<13> -rays via hadronic p-p collisions.  \nAn X-ray study of the SNR with BeppoSAX and ROSAT did not find evidence for shell-like, nonthermal emission, nor thermal X-ray emission that should trace the interaction between the SNR and ISM (Bocchino et al. 2001). However, several hard Xray sources were found, suggestive of at least localized nonthermal electron acceleration. Additional MWL observations and spectral modeling are required to further investigate the scenario responsible for the production of VHE GLYPH<13> -rays.", '5.6.3. HESS J1458 GLYPH<0> 608': "VHE GLYPH<13> -ray emission from the new source HESS J1458 GLYPH<0> 608 (Fig. 20) is associated with the pulsar PSR J1459 GLYPH<0> 6053 and can likely be identified as a heretofore undetected PWN, on the basis of a spatial coincidence with an energetic pulsar and the absence of other plausible MWL counterparts. Preliminary VHE morphological and spectral properties were first announced by de los Reyes et al. (2012). The updated morphological properties from the HGPS catalog di GLYPH<11> er from those preliminary ones, which had underestimated the extent of the large, complex emission region (0 : 37 GLYPH<14> GLYPH<6> 0 : 03 GLYPH<14> vs. 0 : 17 GLYPH<14> GLYPH<6> 0 : 07 GLYPH<14> ; both morphological models 2D symmetric Gaussian), likely due to the irregular shape of the emission. Previously there was a hint for additional structure, possibly a second source hidden in the tail of a dominant source,  \nbut this remains statistically insignificant in the HGPS analysis with respect to a single-source Gaussian morphology. Also of note, the best-fit centroid of the VHE emission is now located closer to the GLYPH<13> -ray pulsar (0 : 11 GLYPH<14> versus 0 : 16 GLYPH<14> o GLYPH<11> set), bolstering the scenario in which the VHE emission is interpreted as a PWN powered by the pulsar. As expected for such changes in morphological properties, the HGPS spectral results also di GLYPH<11> er from the previously derived preliminary values.  \nThe pulsar PSR J1459 GLYPH<0> 6053 (also 3FGL J1459.4 GLYPH<0> 6053) is a relatively old ( GLYPH<28> c = 65 kyr) but still very energetic HE GLYPH<13> -ray pulsar with a spin-down luminosity 9 : 1 GLYPH<2> 10 35 erg s GLYPH<0> 1 and unknown distance ( d < 22 kpc) (Abdo et al. 2013). As noted above, it is o GLYPH<11> set 0 : 11 GLYPH<14> from the VHE centroid, which is consistent with o GLYPH<11> sets observed in other PSR and VHE PWN systems (e.g., Kargaltsev et al. 2013). The putative PWN has not been detected in X-rays potentially because of the age of the system (Ray et al. 2011) or HE GLYPH<13> -rays (Acero et al. 2013).  \nThe new VHE spectrum ( E > 0 : 46 TeV) is consistent with the 31-316 GeV Fermi -LAT upper limits. However, the conclusion, made by Acero et al. (2013), that the peak of the PWN's inverse Compton emission is located in this energy range has to be revised as the peak can now only be inferred to be at higher energies.  \nApart from the HE GLYPH<13> -ray pulsar, there is a second HE source (3FGL J1456.7 GLYPH<0> 6046) in the FoV. However, it is unclear if it is related to the PSR and PWN scenario, since it exhibits a highly curved, log-parabolic spectrum typical of blazars and a TS that fluctuates strongly with the choice of di GLYPH<11> use model or analysis method (Acero et al. 2015a).", '5.6.4. HESS J1503 GLYPH<0> 582': 'HESS J1503 GLYPH<0> 582 (Fig. 21) is a new source for which the origin of the VHE GLYPH<13> -ray emission is unidentified. H.E.S.S. earlier announced preliminary morphological and spectral properties for this source (Renaud et al. 2008), which are now superseded by those in this paper. The VHE emission appears to be one of the softest ( GLYPH<0> = 2 : 68 GLYPH<6> 0 : 08stat) in the HGPS, although both its morphological and spectral properties are a GLYPH<11> ected by systematic uncertainties larger than nominal (see, e.g., Sect. 4.12 and Sect. 5.5).  \nA point-like HE ( E > 50 GeV) GLYPH<13> -ray source, 2FHL J1505.1 GLYPH<0> 5808 (Ackermann et al. 2016), is spatially coincident with the VHE emission region. A comparison of the VHE and HE ( E > 50 GeV) spectra suggests that it may be a PWN (Ackermann et al. 2016), although no PWN or energetic pulsar has been detected so far. Another, di GLYPH<11> erent, point-like HE ( E > 100 MeV) GLYPH<13> -ray source, 3FGL J1503.5 GLYPH<0> 5801 (Acero et al. 2015a), is also within the VHE region. Its nature is unknown, but its log-parabolic spectrum suggests it may not be directly related to HESS J1503 GLYPH<0> 582.  \nFig. 21. VHE GLYPH<13> -ray image: HESS J1503 GLYPH<0> 582. See Fig. 18 for a general description. Additionally, the ellipse represents the 95% uncertainty in the position of the HE GLYPH<13> -ray point source 3FGL J1503.5 GLYPH<0> 5801; the circle represents the 68% uncertainty in the position of the HE ( E > 50 GeV) GLYPH<13> -ray point source 2FHL J1505.1 GLYPH<0> 5808; and the star represents the location of the X-ray point source. The FoV is 2 : 3 GLYPH<14> GLYPH<2> 2 : 3 GLYPH<14> .  \n<!-- image -->  \nFaint X-ray emission (AX J1504.6 GLYPH<0> 5824, Sugizaki et al. 2001) is present toward the edge of the VHE emission. Nominally cataloged as a cataclysmic variable, its X-ray properties are not well known owing to the low ASCA sensitivity. Analysis of more sensitive data from other X-ray telescopes is needed to investigate the possibility it may be a PWN; this is the case despite a lack of an energetic pulsar in the vicinity, but bearing in mind the unknown nature of the nearby 3FGL source.  \nA relatively comprehensive search of MWL archives (Renaud et al. 2008) led to the investigation of an atypical scenario where the VHE emission could be linked with a forbidden velocity wing (FVW): faint, characteristic 21-cm H i line emission structures seen as deviations from the canonical Galactic rotation curve (Kang & Koo 2007). The hypothesis is that this FVW, FVW 319.8 + 0.3, may be related to an older SNR in its radiative phase, as was the case for two other FVWs (Koo et al. 2006; Kang et al. 2012). Although the SNR would no longer have sufficient shock velocity to accelerate particles responsible for producing VHE GLYPH<13> -rays (Ptuskin & Zirakashvili 2005), it could nevertheless be indicative of increased or more recent activity in the region (stellar winds and / or supernova explosions). A large H i shell, the result of such activity, is nearby (see McClure-Gri GLYPH<14> ths et al. 2002, GSH 319 GLYPH<0> 01 + 13); however, its centroid is substantially o GLYPH<11> set by more than 1 GLYPH<14> from HESS J1503 GLYPH<0> 582 and its  \nFig. 22. VHE GLYPH<13> -ray image: HESS J1554 GLYPH<0> 550. See Fig. 18 for a general description. Additionally, the composite SNR is shown by plotting its 843-MHz radio intensity (Whiteoak & Green 1996) with contours at 1, 8, and 15 mJy / beam. The FoV is 0 : 9 GLYPH<14> GLYPH<2> 0 : 9 GLYPH<14> .  \n<!-- image -->  \nextent considerably larger than the VHE emission region, so it seems unlikely to be related.  \nOn the other hand, VERITAS also searched for VHE emission from an FVW, one which does show clear shell-type emission in H i (FVW 190.2 + 1.1). Despite observations that reached a sensitivity better than 1% Crab, VERITAS did not detect any significant VHE emission (Holder 2009). Furthermore, there is no definitive identification of VHE emission from young stellar clusters, with the possible exception of the superbubble 30 Dor C in the LMC (H.E.S.S. Collaboration et al. 2015d). Therefore, this FVW scenario remains speculative.', '5.6.5. HESS J1554 GLYPH<0> 550': 'VHE GLYPH<13> -ray emission from the new source HESS J1554 GLYPH<0> 550 (Fig. 22) is firmly identified with the PWN G327.15 GLYPH<0> 1.04 within the composite SNR G327.1 GLYPH<0> 1.1, on the basis of both a spatial coincidence with the PWN and the size of the VHE emission region, which can be constrained to less than 0 : 035 GLYPH<14> . Preliminary H.E.S.S. morphological and spectral properties of the VHE source were first published by Acero et al. (2012) and are compatible with the HGPS results. However, while previously the source size was given as 0 : 03 GLYPH<14> GLYPH<6> 0 : 01 GLYPH<14> , the more conservative HGPS analysis procedure used here (Sect. 4.10.2) finds the source to be compatible with being point-like, and there is a limit on the size that is nonetheless compatible.  \nFig. 23. VHE GLYPH<13> -ray image: HESS J1813 GLYPH<0> 126. See Fig. 18 for a general description. Additionally, the position of the pulsar is denoted by a cyan diamond. The FoV is 1 : 7 GLYPH<14> GLYPH<2> 1 : 7 GLYPH<14> .  \n<!-- image -->  \nThe VHE size limit rules out significant emission from the outer shell of the SNR and is compatible with the compact, nonthermal PWN, which is observed in both radio and X-rays (Temim et al. 2009, 2015) but not HE GLYPH<13> -rays (Acero et al. 2013; Ackermann et al. 2016). Furthermore, the VHE centroid is compatible with the peak of the radio emission from the PWN and the tail of X-ray PWN. Although pulsed emission from the putative pulsar at the heart of the composite SNR has not been detected in radio, X-ray, nor HE GLYPH<13> -ray bands, the X-ray data provide evidence for the existence of a powerful pulsar, that has an estimated ˙ E = 3 : 1 GLYPH<2> 10 36 erg s GLYPH<0> 1 (Temim et al. 2015). The distance to the SNR is not well determined, but has been estimated to be roughly 9 kpc (Sun et al. 1999). Assuming this distance, the VHE luminosity of HESS J1554 GLYPH<0> 550 is L GLYPH<13> (1 GLYPH<0> 10 TeV) = 1 : 0 GLYPH<2> 10 34 ( d = 9 kpc) 2 erg s GLYPH<0> 1 and the apparent e GLYPH<14> ciency GLYPH<15> 1 GLYPH<0> 10TeV GLYPH<17> L GLYPH<13> = ˙ E = 0 : 3%, which is compatible with the e GLYPH<14> ciencies ( < GLYPH<24> 10%) of other VHE sources that have been identified as PWNe (Kargaltsev et al. 2013).', '5.6.6. HESS J1813 GLYPH<0> 126': 'The HGPS catalog analysis has revealed an intriguing new source of VHE GLYPH<13> -rays (Fig. 23) not previously detected, one of the few o GLYPH<11> -plane VHE sources ( b = 2 : 5 GLYPH<14> ). The only plausible MWL counterpart associated with this emission is the energetic pulsar PSR J1813 GLYPH<0> 1246 (Abdo et al. 2009a), marginally coincident with the VHE best-fit centroid. This suggests the VHE emission originates in a PWN powered by the pulsar, which has  \na spin-down luminosity ˙ E = 6 : 3 GLYPH<2> 10 36 erg s GLYPH<0> 1 and characteristic age GLYPH<28> c = 43 kyr. The pulsar is one of the brightest GLYPH<13> -ray pulsars (3FGL J1813.4 GLYPH<0> 1246, Acero et al. 2015a) and the secondmost energetic, radio-quiet pulsar. This pulsar also been found to exhibit strong X-ray pulsations, and its distance has been recently constrained to d > 2 : 5 kpc (Marelli et al. 2014). This implies a lower limit on the VHE luminosity L GLYPH<13> (1 GLYPH<0> 10 TeV) > 2 : 9 GLYPH<2> 10 33 erg s GLYPH<0> 1 and a corresponding limit on the apparent efficiency GLYPH<15> 1 GLYPH<0> 10TeV > 0 : 05%.  \nIn other energy bands, no o GLYPH<11> -pulse emission (e.g. emission from the putative PWN) is detected in HE GLYPH<13> -rays (0.1-100 GeV) based on the analysis of five years of Fermi -LAT data (Marelli et al. 2014), dismissing earlier hints for a GeV PWN (Ackermann et al. 2011; Abdo et al. 2013), likely owing to the larger data set and improved models for di GLYPH<11> use emission used in the new analysis. In X-rays (0 : 3 GLYPH<0> 10 keV), despite relatively deep XMM-Newton (130 ks) and Chandra (50 ks) observations, no PWN is detected beyond 1 GLYPH<0> 1 : 5 00 from the pulsar (Marelli et al. 2014). This is very unusual for a pulsar this energetic; that is the derived upper limits in X-rays are only marginally compatible with known relations between PWN and pulsar luminosities (Kargaltsev & Pavlov 2008) and between PSR luminosity and distance to the PWN termination shock (Gaensler & Slane 2006).  \nTherefore, HESS J1813 GLYPH<0> 126 appears to be a rare case of a relic PWN (de Jager & Djannati-Ataï 2009) currently detected exclusively in the VHE domain. Observations in the hard Xray domain with NuSTAR would be useful to investigate the hint of a signal seen at 30-520 keV with INTEGRAL (Marelli et al. 2014) and to determine if there is an unpulsed, nebular component visible at those energies. Regardless, further work modeling the MWL spectral energy distribution is necessary to fully investigate this intriguing system.', '5.6.7. HESS J1826 GLYPH<0> 130': "The HGPS catalog analysis reveals a distinct new source of VHE GLYPH<13> -rays, HESS J1826 GLYPH<0> 130 (Fig. 24), in what was previously considered extended emission from the nearby PWN HESS J1825 GLYPH<0> 137 (Aharonian et al. 2006g). Because of the very close proximity to its bright neighbor, the spectral measurement is highly contaminated (41%). Angüner et al. (2017) reported preliminary findings for this new source.  \nHESS J1826 GLYPH<0> 130 is associated with the 'Eel' PWN 11 (PWN G18.5 GLYPH<0> 0.4), an elongated, nonthermal, X-ray source observed with Chandra (Roberts et al. 2007), and the energetic pulsar PSR J1826 GLYPH<0> 1256 (Abdo et al. 2009a), on the basis of a spatial coincidence. The best-fit VHE centroid is compatible with the Eel, while the pulsar is somewhat o GLYPH<11> set (0 : 09 GLYPH<14> ) from the centroid  \nFig. 24. VHE GLYPH<13> -ray image: HESS J1826 GLYPH<0> 130. See Fig. 18 for a general description. Additionally, the approximate centroid of the PWN is marked by a green cross, and the pulsar position is marked by a cyan diamond. The FoV is 2 : 0 GLYPH<14> GLYPH<2> 2 : 0 GLYPH<14> .  \n<!-- image -->  \nbut well within the VHE emission region (size 0 : 15 GLYPH<14> GLYPH<6> 0 : 02 GLYPH<14> ). The pulsar is now notable for being one of the brightest radioquiet GLYPH<13> -ray pulsars (3FGL J1826.1 GLYPH<0> 1256; Acero et al. 2015a). The distance of the pulsar is unfortunately not known, which precludes conclusions on the energetics, but its position, ˙ E = 3 : 6 GLYPH<2> 10 36 erg s GLYPH<0> 1 , and GLYPH<28> c = 14 kyr suggest it is probably powering the Eel. The PWN is not detected in HE GLYPH<13> -rays (Ackermann et al. 2011; Ackermann et al. 2016). Finally, we note that dense molecular gas was also found overlapping HESS J1826 GLYPH<0> 130 at a distance matching that of the dispersion measure of the pulsar (Voisin et al. 2016), suggesting a possible hadronic origin for this VHE source.  \nThe SNR G18.6 GLYPH<0> 0.2 (Brogan et al. 2006) is also coincident with the VHE emission region, although it is significantly smaller in size (0 : 1 GLYPH<14> diameter). Very little is known about this SNR, except that a partial shell-type morphology has been observed so far only in radio and IR and that its distance is estimated to be 4.0-5.2 kpc (Johanson & Kerton 2009).  \nA firm identification of the VHE source as a PWN is not possible at this time, in part resulting from the unknown distance to the Eel PWN and PSR system and the poorly studied SNR. We are currently preparing more advanced VHE spectral analysis methods that can account for contamination in crowded FoVs. These methods will enable more accurate modeling of the SED.  \nFig. 25. VHE GLYPH<13> -ray image: HESS J1828 GLYPH<0> 099. See Fig. 18 for a general description. The FoV is 1 : 0 GLYPH<14> GLYPH<2> 1 : 0 GLYPH<14> .  \n<!-- image -->", '5.6.8. HESS J1828 GLYPH<0> 099': 'HESS J1828 GLYPH<0> 099 is a new source of VHE GLYPH<13> -rays (Fig. 25), which is unique because it appears to be completely dark at lower energies with no apparent associations (see Table 12). It is also notable for being one of the 17 point-like sources in the HGPS catalog with a size (Gaussian std. dev.) less than 0 : 07 GLYPH<14> . The detection of a spatially coincident HE GLYPH<13> -ray source has been claimed (Neronov & Semikoz 2010) but not confirmed with the latest, significantly larger Fermi -LAT data sets, i.e. there is no respective source neither in the 3FGL catalog (Acero et al. 2015a) nor 2FHL catalog ( E > 50 GeV; Ackermann et al. 2016). Deeper follow-up observations in especially the radio and X-ray bands are strongly encouraged to probe nonthermal emission.', '5.6.9. HESS J1832 GLYPH<0> 085': 'HESS J1832 GLYPH<0> 085 (Fig. 26) is an unidentified source of VHE GLYPH<13> -rays. It is notable for its point-like morphology, which is measured to be less than 0 : 05 GLYPH<14> in extension, and its scarcity of promising MWL counterparts.  \nAn interesting object that is spatially coincident with HESS J1832 GLYPH<0> 085 is the pulsar PSR J1832 GLYPH<0> 0827 (Clifton & Lyne 1986), which has so far only been detected in radio wavelengths. The pulsar is likely at a distance of GLYPH<25> 4.9 kpc (Cordes & Lazio 2002), in agreement with other estimates in the range  \nFig. 26. VHE GLYPH<13> -ray image: HESS J1832 GLYPH<0> 085. See Fig. 18 for a general description. Additionally, the position of the two pulsars are denoted by cyan diamonds. The FoV is 0 : 7 GLYPH<14> GLYPH<2> 0 : 7 GLYPH<14> .  \n<!-- image -->  \n4.4-6.1 kpc (Frail et al. 1991) and has a spin-down luminosity 12 ˙ E = 9 : 3 GLYPH<2> 10 33 erg s GLYPH<0> 1 (Hobbs et al. 2004b). It is one of the few pulsars with a measured braking index, n = 2 : 5 GLYPH<6> 0 : 9 (Johnston & Galloway 1999), providing a characteristic age GLYPH<28> c GLYPH<25> 200 kyr. Another very intriguing object in the FoV is the energetic millisecond pulsar PSR J1832 GLYPH<0> 0836, which has a 2.7 ms period (Burgay et al. 2013). It has a spin-down luminosity ˙ E = 1 : 7 GLYPH<2> 10 34 erg s GLYPH<0> 1 , a very large characteristic age (typical of millisecond pulsars) GLYPH<28> c = 5 GLYPH<2> 10 9 kyr, and distance 1.1 kpc (Cordes & Lazio 2002).  \nThere are no known PWNe associated with these two pulsars nor close to HESS J1832 GLYPH<0> 085. If either or both of these pulsars are powering VHE PWNe, a relatively large conversion e GLYPH<14> ciency of GLYPH<15> 1 GLYPH<0> 10 TeV GLYPH<24> 23% would be required for PSR J1832 GLYPH<0> 0827, and a more reasonable GLYPH<15> 1 GLYPH<0> 10 TeV GLYPH<24> 0 : 6% for PSR J1832 GLYPH<0> 0836. The older ages are at odds with the inferred small sizes of the VHE PWNe, constrained to be less than GLYPH<25> 4 ( d = 4 : 9 kpc) pc and GLYPH<25> 1 ( d = 1 : 1 kpc) pc, respectively. These circumstances, plus the borderline low spin-down luminosity of PSR J1832 GLYPH<0> 0827, combine to disfavor a PSR and PWN scenario as the origin of the VHE emission in light of the known VHE PWN population (H.E.S.S. Collaboration et al. 2017f). The millisecond pulsar scenario is even more uncertain. That pulsar is slightly more energetic and much closer, but thus far millisecond pulsars, with ages of billions of years, are not known to produce PWNe that emit detectable levels of GLYPH<13> -rays at TeV energies.  \nFig. 27. VHE GLYPH<13> -ray image: HESS J1833 GLYPH<0> 105. See Fig. 18 for a general description. Additionally, the composite SNR is shown by plotting a green circle approximating the radio shell, and the pulsar position is denoted by a cyan diamond. The FoV is 0 : 8 GLYPH<14> GLYPH<2> 0 : 8 GLYPH<14> .  \n<!-- image -->  \nTherefore, the origin of the emission from this new, enigmatic, VHE GLYPH<13> -ray source is still very much unclear.', '5.6.10. HESS J1833 GLYPH<0> 105': 'VHE GLYPH<13> -ray emission from the new source HESS J1833 GLYPH<0> 105 (Fig. 27) can now be firmly identified with the composite SNR G21.5 GLYPH<0> 0.9 (Wilson & Weiler 1976), which contains a Crab-like PWN. Preliminary H.E.S.S. source properties were previously shown (Djannati-Ataï et al. 2008) and are compatible with the HGPS results, although at the time it was not yet possible to disentangle the possible contributions to the VHE emission from the PWN and SNR shell.  \nThe new identification is supported by a positional coincidence between the VHE emission centroid and the PWN center, but most importantly by the lack of extension of the VHE emission region; this region is constrained to be less than 0.03 GLYPH<14> , which is our systematic limit for source sizes. This implies that we cannot claim significant VHE emission from the forward shock of the spherical, faint SNR shell at a radius of 0.038 GLYPH<14> (Bocchino et al. 2005).  \nThe PWN has also been detected in X-rays (Safi-Harb et al. 2001; Bocchino et al. 2005) and IR (Zajczyk et al. 2012) although not in HE GLYPH<13> -rays (Acero et al. 2013), and its distance has been estimated to be d GLYPH<25> 4 : 8 kpc (Tian & Leahy 2008). It is powered by the very energetic PSR J1833 GLYPH<0> 1034, currently the fifth most energetic pulsar known in the Galaxy, and has a spin-  \nFig. 28. VHE GLYPH<13> -ray image: HESS J1843 GLYPH<0> 033. See Fig. 18 for a general description. Additionally, the SNR is shown by plotting a green circle approximating the radio shell. The 3FGL ellipses represents the 95% uncertainty in the position of the HE GLYPH<13> -ray point sources. The FoV is 1 : 4 GLYPH<14> GLYPH<2> 1 : 4 GLYPH<14> .  \n<!-- image -->  \ndown luminosity ˙ E GLYPH<25> 3 : 4 GLYPH<2> 10 37 erg s GLYPH<0> 1 . The pulsar has been detected in radio (Gupta et al. 2005; Camilo et al. 2006) and HE GLYPH<13> -rays (as 3FGL J1833.5 GLYPH<0> 1033; Acero et al. 2015a). The age of the system has been argued to be 870 + 200 GLYPH<0> 150 yr (Bietenholz & Bartel 2008), which is significantly less than the GLYPH<28> c = 4 : 9 kyr of the pulsar.  \nConsidering the luminosity of HESS J1833 GLYPH<0> 105, L GLYPH<13> (1 GLYPH<0> 10 TeV) = 2 : 6 GLYPH<2> 10 33 ( d / 4.9 kpc) 2 erg s GLYPH<0> 1 , the apparent e GLYPH<14> -ciency converting the rotational energy of the pulsar to GLYPH<13> -rays, GLYPH<15> 1 GLYPH<0> 10TeV GLYPH<17> L GLYPH<13> = ˙ E = 0 : 08%, is compatible with the e GLYPH<14> ciencies ( < GLYPH<24> 10%) of other VHE sources that have been identified as PWNe (Kargaltsev et al. 2013).  \nThe HGPS results confirm predictions that the PWN would emit VHE GLYPH<13> -rays at the level of a few percent of the Crab Nebula and exhibit a relatively hard spectrum (de Jager et al. 1995).', '5.6.11. HESS J1843 GLYPH<0> 033': "An extended region of VHE emission, called HESS J1843 GLYPH<0> 033, was first published by Hoppe (2008b). This emission is resolved by the HGPS catalog analysis into three components that were merged into two distinct sources: HESS J1843 GLYPH<0> 033 and HESS J1844 GLYPH<0> 030.  \nHESS J1843 GLYPH<0> 033 consists of two merged o GLYPH<11> set components (HGPSC 83 and HGPSC 84) and is therefore highly structured (see Fig. 28). The image of the source shows two peaks separated by GLYPH<24> 0.2 GLYPH<14> . The first Gaussian component is clearly as-  \nsociated with the upper peak. The second Gaussian component is larger and o GLYPH<11> set with respect to the lower peak. This is due to more di GLYPH<11> use, low-brightness emission around ( '; b ) = (28 : 6 GLYPH<14> ; GLYPH<0> 0 : 1 GLYPH<14> ), suggesting the presence of another currently unresolved source that shifts the position of the second component. HESS J1843 GLYPH<0> 033 is therefore most probably a complex region with overlapping sources that were merged in the HGPS analysis.  \nTwo GeV sources, 3FGL J1843.7 GLYPH<0> 0322 and 3FGL J1844.3 GLYPH<0> 0344, are found within the R80 extension of the source. The former is found in the main region of emission but does not seem to correlate well with any of the two main peaks. The latter Fermi -LAT source is located in the low-brightness region around ( '; b ) = (28 : 6 GLYPH<14> ; GLYPH<0> 0 : 1 GLYPH<14> ).  \nNo compelling radio counterpart was found in the VLA Galactic Plane Survey (Stil et al. 2006). Dedicated X-ray observations show the presence of a faint absorbed extended source with a nonthermal spectrum that is coincident with the HGPSC 83 component. No compelling counterpart for the second component has been found. We note however that the nearby radio SNR G28.6 GLYPH<0> 0.1 (Helfand et al. 1989) is filled with nonthermal X-rays (Ueno et al. 2003). If this emission is due to synchrotron X-rays produced by energetic electrons, IC emission at VHE is likely contributing to the low-brightness emission that is visible around the SNR position in Fig. 28.", '5.6.12. HESS J1844 GLYPH<0> 030': 'HESS J1844 GLYPH<0> 030 is a faint VHE GLYPH<13> -ray source that compatible with being point-like and located in the vicinity of the complex region of HESS J1843 GLYPH<0> 033. It is positionally coincident with a number of distinct objects, most notably the radio source PMN J1844 GLYPH<0> 0306 (cyan diamond in Fig. 29). The nature of the latter is ambiguous. Its elongated, jet-like morphology is very reminiscent of a radio galaxy, which is supported by 6 cm VLA observations revealing polarization along the structure (Helfand et al. 1989). This elongated radio feature is surrounded by a partial ring visible in the 21 cm VLA continuum image. The object is therefore classified as a SNR candidate in the MAGPIS catalog (Helfand et al. 2006), G29.37 + 0.10. It is also coincident with the X-ray source AX J1844.7 GLYPH<0> 0305 (Vasisht et al. 2000; Sugizaki et al. 2001).  \nThe association of the jet radio feature and the SNR candidate is unclear. Although rare, SNRs with jets are plausible, for example PWN structures such as MSH 15 GLYPH<0> 52 (Gaensler et al. 2002) or the SS433 / W 50 microquasar SNR system with its radio jets and lobes (Dubner et al. 1998; H.E.S.S. Collaboration et al. 2017d). The jet structure could also be a background radio galaxy aligned by chance with a faint radio shell. We note, however, thanks to H i absorption, Johanson & Kerton (2009) place the source at a distance between 5 and GLYPH<24> 15 kpc. Interestingly, a  \nFig. 29. VHE GLYPH<13> -ray image: HESS J1844 GLYPH<0> 030. See Fig. 18 for a general description. Additionally, the composite SNR is shown by plotting a green circle approximating the radio shell. The position of the X-ray source is denoted by a cross, while the position of PMN J1844 GLYPH<0> 0306 is indicated by a diamond. The FoV is 1 : 0 GLYPH<14> GLYPH<2> 1 : 0 GLYPH<14> .  \n<!-- image -->  \nheavily absorbed X-ray PWN, dubbed G29.4 + 0.1, is present in SNRcat and overlaps with a part of PMN J1844 GLYPH<0> 0306. Further MWL observations will be necessary to assess the nature of the system and the origin of the VHE emission.', '5.6.13. HESS J1846 GLYPH<0> 029': 'VHE GLYPH<13> -ray emission from the new source HESS J1846 GLYPH<0> 029 (see Fig. 30) is spatially coincident with G29.7 GLYPH<0> 0.3 (also known as Kes 75), one of the youngest composite SNRs in the Galaxy, which contains the nebula of PSR J1846 GLYPH<0> 0258. Preliminary results were presented in Djannati-Ataï et al. (2008) and are compatible with those obtained in the HGPS analysis.  \nPSR J1846 GLYPH<0> 0258 is a young, high magnetic-field pulsar. This source has a rotation period of 324 ms and a spin-down power of 8 : 3 GLYPH<2> 10 36 erg s GLYPH<0> 1 . It is among the youngest pulsars in the Galaxy with a characteristic age of only 723 years (Livingstone et al. 2006). It has experienced a strong increase in its pulsed flux in June 2006 associated with spectral (Kumar & SafiHarb 2008) and timing (Gavriil et al. 2008) changes in a similar manner to magnetars. The result of the search for variations in the VHE source flux at various timescales was negative (Terrier et al. 2008a).  \nA nebula of 20 00 in radius surrounds the pulsar in radio and X-ray wavelengths and Chandra high-resolution observations have revealed a jet and torus (Ng et al. 2008). A 3 0 diameter  \nFig. 30. VHE GLYPH<13> -ray image: HESS J1846 GLYPH<0> 029. See Fig. 18 for a general description. Additionally, the composite SNR is shown by plotting a green circle approximating the radio shell. The position of the pulsar is indicated by a cyan diamond. The FoV is 0 : 8 GLYPH<14> GLYPH<2> 0 : 8 GLYPH<14> .  \n<!-- image -->  \nasymmetric radio shell surrounds the PSR and PWN system. It consists mainly of two lobes to the south of the pulsar. These lobes are emitting X-rays from heated swept-up interstellar matter and ejecta (Morton et al. 2007). Infrared measurements suggest that the shock is in a region of typical density of 60 cm GLYPH<0> 3 (Temim et al. 2012). Su et al. (2009) found a bubble in the molecular matter in good coincidence with the SNR. They proposed that this structure is the wind blown bubble of the SNR progenitor.  \nThe extension of the VHE emission from HESS J1846 GLYPH<0> 029 is compatible with that of a point-like source. The upper limit on the size is 0 : 03 GLYPH<14> , that is, comparable with the SNR shell size. The position of this object is compatible with the position of PSR J1846 GLYPH<0> 0258, within localization uncertainties. Therefore, we are not able to distinguish between emission from the shell and emission from the PWN in this composite object.  \nAssuming a distance of 6 kpc (Leahy & Tian 2008), which yields a luminosity of L GLYPH<13> (1 GLYPH<0> 10 TeV) = 6 : 9 GLYPH<2> 10 33 ( d / 6 kpc) 2 erg s GLYPH<0> 1 , the apparent conversion e GLYPH<14> ciency of the rotational energy of the pulsar to GLYPH<13> -rays is GLYPH<15> 1 GLYPH<0> 10TeV GLYPH<17> L GLYPH<13> = ˙ E = 0 : 08%. The VHE emission is therefore completely consistent with an origin in the PWN (see also, e.g., Tanaka & Takahara 2011; Torres et al. 2014). Yet, given the uncertainties on extension it is not possible to exclude a contribution from GLYPH<13> -rays produced by particles accelerated at the SNR shock, in particular from collisions of hadrons with ambient and swept-up matter at the shock, or even  \nFig. 31. VHE GLYPH<13> -ray image: HESS J1848 GLYPH<0> 018. See Fig. 18 for a general description. Additionally, the 3FGL ellipse represents the 95% uncertainty in the position of the HE GLYPH<13> -ray point source. The central stellar cluster of the mini-starburst W 43 is denoted by a diamond. The FoV is 2 : 1 GLYPH<14> GLYPH<2> 2 : 1 GLYPH<14> .  \n<!-- image -->  \na contribution of escaping particles with the molecular shell revealed by Su et al. (2009).', '5.6.14. HESS J1848 GLYPH<0> 018': 'For the new source HESS J1848 GLYPH<0> 018 (Fig. 31) preliminary H.E.S.S. source properties were previously shown (Chaves et al. 2008b). These properties are compatible with the HGPS results except for the source size and flux; these were overestimated because the earlier analysis did not include a model for the di GLYPH<11> use emission (see Sect. 4.6), which is particularly bright in this region.  \nThe origin of the VHE GLYPH<13> -ray emission of HESS J1848 GLYPH<0> 018 is not yet firmly identified. No SNR or energetic pulsar is currently detected in the proximity, although we have associated the VHE source with 3FGL J1848.4 GLYPH<0> 0141 (Acero et al. 2015a). This unidentified HE GLYPH<13> -ray point source is significantly o GLYPH<11> set from the VHE GLYPH<13> -ray centroid (by GLYPH<24> 0.2 GLYPH<14> ) but well within the VHE emission region. Studies attempting to relate the HE with the (preliminary) VHE morphology and spectra remained inconclusive (Tam et al. 2010; Acero et al. 2013). A potential PSR and PWNscenario cannot be confirmed due to the lack of a detected pulsar (at any wavelength), although the HE spectrum does exhibit curvature typical of pulsars (Acero et al. 2015a). Furthermore, there is no known PWN nearby, although one study has shown marginal statistical evidence for an extension of the HE source (Lemoine-Goumard et al. 2011), which is expected if the  \nHE emission is from a PWN or the combination of a pulsar and PWN.  \nAn extensive search for other MWL counterparts found the VHE GLYPH<13> -ray emission to be in the direction of the massive starforming region W 43, a very active mini-starburst located at a distance of 6 : 2 GLYPH<6> 0 : 6 kpc (Russeil 2003). It is one of the closest and most luminous star-forming regions in the Galaxy (Motte et al. 2003), hosting a giant H ii region (G30.8 GLYPH<0> 0.2), a giant molecular cloud, and the Wolf-Rayet binary star system WR121a in the central stellar cluster together with O-type stars. The massive stars in the dense central cluster exhibit strong stellar winds with extreme mass loss rates, in particular the WN7subtype WR 121a (Blum et al. 1999).  \nThis unique MWL environment is of interest because the central cluster of W 43 could be the site of e GLYPH<14> cient particle acceleration in various plausible hadronic scenarios involving the high-velocity (up to 2000 km s GLYPH<0> 1 ) stellar winds (e.g., Reimer et al. 2006; Romero 2010). Furthermore, the very large amount of molecular gas present in W 43 ( GLYPH<24> 7 GLYPH<2> 10 6 M GLYPH<12> ; Nguyen Luong et al. 2011) provides a natural target for accelerated cosmic rays (regardless of their potential acceleration site), which would lead to GLYPH<13> -ray production via hadronic p-p collisions (e.g., Aharonian 1991).  \nIt is not yet possible to confirm the W 43 hadronic scenario for the origin of the VHE emission, in part because of the very complex morphologies present and the challenges in correlating features observed in radio and infrared observations at arcsecond scales with the GLYPH<24> 5 0 resolution in VHE. The VHE centroid, in particular, is significantly o GLYPH<11> set from the central cluster by GLYPH<24> 0.2 GLYPH<14> , although the extended VHE emission is generally coincident with the W 43 complex. This scenario remains under investigation, especially in light of the recent detection of the superbubble 30 Dor C in the LMC (H.E.S.S. Collaboration et al. 2015d), which suggests that particle acceleration occurring in the collective winds of massive stars can indeed produce VHE emission.', '5.6.15. HESS J1849 GLYPH<0> 000': 'The faint, slightly extended source HESS J1849 GLYPH<0> 000 (see Fig. 32) was first reported by Terrier et al. (2008b). It was found to be spatially coincident with the hard X-ray source IGR J18490 GLYPH<0> 0000 (Krivonos et al. 2012). XMM-Newton observations revealed a nonthermal, point-like, X-ray source surrounded by a nebula, making this object a solid PWN candidate. Follow-up observations of the hard X-ray source with RXTE have confirmed this hypothesis with the discovery of a 38.5 ms periodicity of the X-ray signal (Gotthelf et al. 2011). The associated pulsar, PSR J1849 GLYPH<0> 0001, was found to have a spindown luminosity ˙ E = 9 : 8 GLYPH<2> 10 36 erg s GLYPH<0> 1 and a characteristic age GLYPH<28> c = 42 : 9 kyr.  \nFig. 32. VHE GLYPH<13> -ray image: HESS J1849 GLYPH<0> 000. See Fig. 18 for a general description. Additionally, the position of the pulsar is indicated by a cyan diamond and the PWN by a green cross. The FoV is 1 : 5 GLYPH<14> GLYPH<2> 1 : 5 GLYPH<14> .  \n<!-- image -->  \nThe HGPS analysis confirms the existence of a source coincident with PSR J1849 GLYPH<0> 0001. The best-fit position of HESS J1849 GLYPH<0> 000 is located less than 0 : 03 GLYPH<14> from the X-ray pulsar position (cyan diamond on Fig. 32), well within statistical uncertainties in both source localizations. The best-fit size of the VHE emission is 0 : 09 GLYPH<14> , which is about a factor of two larger than that of the extended X-ray component (Gotthelf et al. 2011; Kuiper & Hermsen 2015).  \nThe source has an energy flux GLYPH<24> 2.1 GLYPH<2> 10 GLYPH<0> 12 erg cm GLYPH<0> 2 s GLYPH<0> 1 , in the range 1-10 TeV, a factor of GLYPH<24> 2 above the X-ray nebula energy flux in the range 2-10 keV (Kuiper & Hermsen 2015). This confirms the likely nature of HESS J1849 GLYPH<0> 000 as a PWN in transition between a young, synchrotron-dominated phase and an evolved, IC-dominated phase.', '5.6.16. HESS J1852 GLYPH<0> 000': "The new source of VHE GLYPH<13> -ray emission HESS J1852 GLYPH<0> 000 (Fig. 33) is currently unidentified due to multiple source counterpart confusion. It is spatially associated with the partial shell-type SNR G32.8 GLYPH<0> 0.1 (also known as Kes 78; Kesteven 1968; Velusamy & Kundu 1974), the incomplete shell-type SNR G33.2 GLYPH<0> 0.6 (Reich 1982), and two energetic pulsars, PSR J1853 GLYPH<0> 0004 and PSR J1853 + 0011 (Hobbs et al. 2004a). Preliminary H.E.S.S. source properties were previously shown (Kosack et al. 2011) and are compatible with the HGPS results. As mentioned in Sect. 5.5, the spectral properties of HESS J1852 GLYPH<0> 000 are a GLYPH<11> ected by systematic uncertainties larger than nominal.  \nFig. 33. VHE GLYPH<13> -ray image: HESS J1852 GLYPH<0> 000. See Fig. 18 for a general description. Additionally, the position of the pulsars are indicated by cyan diamonds, and the SNR is shown by plotting a green circle approximating the radio shell. The FoV is 2 : 0 GLYPH<14> GLYPH<2> 2 : 0 GLYPH<14> .  \n<!-- image -->  \nThe VHE emission is located along the eastern edge of SNR Kes 78 but extends well beyond the SNR. The SNR itself is characterized by an elongated and partial nonthermal shell seen in radio and X-rays (Zhou & Chen 2011; Bamba et al. 2016). It is interacting with adjacent molecular clouds, evidenced by the detection of a shock-excited OH(1720 MHz) maser on the shell (Koralesky et al. 1998) and studies of the CO molecular environment (Zhou & Chen 2011). The distance of the SNR is estimated to be GLYPH<24> 5 kpc (Koralesky et al. 1998; Zhou & Chen 2011), although GLYPH<24> 8.8 kpc has also been suggested (e.g., Xu & Zhang 2009). A hadronic origin of the VHE emission has been briefly discussed (Kosack et al. 2011), involving escaped cosmic rays from Kes 78 (e.g., Aharonian 1991; Gabici et al. 2007). However, the scenario remains unconfirmed in the absence of a more detailed study of the gas environment and its potential correlation with the complex VHE morphology.  \nThe presence of two radio pulsars, PSR J1853 GLYPH<0> 0004 and PSR J1853 + 0011, within the VHE emission region also suggests that the VHE GLYPH<13> -rays could originate in one of the PWN or could even be a result of superimposed emission from two PWNe. Although there are currently no known PWNe at other energies, the pulsars' spin-down luminosities ˙ E = 2 : 1 GLYPH<2> 10 35 erg s GLYPH<0> 1 and 2 : 1 GLYPH<2> 10 34 erg s GLYPH<0> 1 , respectively, and distances d = 6 : 6 kpc and 7.5 kpc, are reasonable in the context of other pulsars thought to be powering VHE PWNe (H.E.S.S. Collaboration et al. 2017f). The pulsars have so far only been detected in radio, although PSR J1853 GLYPH<0> 0004 has been associated with the HE GLYPH<13> -ray source  \n3FGL J1853.2 + 0006, which is itself a source whose existence and properties are currently uncertain (subject to analysis Flags 3 and 4 in Acero et al. 2015a).  \nIn conclusion, it is not yet clear whether the VHE emission originates from a hadronic SNR and molecular cloud interaction, previously undetected PWNe associated with one or both of the spatially coincident pulsars, or some other yet unknown source.", '5.6.17. Source candidates': 'Three VHE GLYPH<13> -ray source candidates (hotspots) were found above the TS = 30 detection threshold in one HGPS analysis (primary or cross-check), but these candidates had TS < 30 in the other analysis. These should be considered unconfirmed, or candidate, VHE sources to be confirmed by deeper VHE observations.', 'HOTS J1111 GLYPH<0> 611': "The VHE emission from the source candidate HOTS J1111 GLYPH<0> 611 has a significance of TS = 22 (cross-check TS = 41). It is located at ( ' , b ) = (291 : 18 GLYPH<14> GLYPH<6> 0 : 03 GLYPH<14> , GLYPH<0> 0 : 54 GLYPH<14> GLYPH<6> 0 : 03 GLYPH<14> ), has a measured integral flux F ( E > 1 TeV) = 3 : 8 GLYPH<2> 10 GLYPH<0> 13 cm GLYPH<0> 2 s GLYPH<0> 1 , and a size 0 : 09 GLYPH<14> GLYPH<6> 0 : 03 GLYPH<14> . It is located near ( GLYPH<24> 0.1 GLYPH<14> ) the very energetic GLYPH<28> c = 32 : 7 kyr pulsar PSR J1112 GLYPH<0> 6103 (Manchester et al. 2001) emitting in radio and HE GLYPH<13> -rays (3FGL J1111.9 GLYPH<0> 6058, Abdo et al. 2013). The pulsar has a high spin-down luminosity ˙ E = 4 : 5 GLYPH<2> 10 36 erg s GLYPH<0> 1 and a distance of 12.2 kpc (Cordes & Lazio 2002). Moreover, a significant HE GLYPH<13> -ray source (2FHL J1112.1 GLYPH<0> 6101e, Ackermann et al. 2016) above 50 GeV has been reported at 0 : 04 GLYPH<14> from the pulsar, which makes this HE source likely to be a PWN. The characteristics of this pulsar, the apparent e GLYPH<14> ciency GLYPH<15> E > 1 TeV GLYPH<24> 1%, and the presence of a HE component in its vicinity suggests that it could plausibly power a VHE PWN.", 'HESS J1831 GLYPH<0> 098': 'The source candidate HESS J1831 GLYPH<0> 098 is found to have TS = 59 in the main HGPS analysis but only TS = 17 in the crosscheck analysis and is therefore considered a source candidate. HESS J1831 GLYPH<0> 098 is located in a complex region with nearby di GLYPH<11> use components, which might explain the discrepancy observed for this source candidate. Preliminary VHE morphological and spectral properties on HESS J1831 GLYPH<0> 098 were announced by Sheidaei (2011). This source candidate is coincident with the energetic pulsar PSR J1831 GLYPH<0> 0952, which exhibits a spindown luminosity of 1 : 1 GLYPH<2> 10 36 erg s GLYPH<0> 1 and a characteristic age of 128 kyr. According to Sheidaei (2011), a GLYPH<15> 1 GLYPH<0> 20 TeV GLYPH<24> 1% conversion e GLYPH<14> ciency from rotational energy to GLYPH<13> -rays would be required to power a PWN; this is similar to values observed in other VHE PWN.', 'HOTS J1907 + 091': "The VHE emission from source candidate HOTS J1907 + 091 has a significance of only TS = 18 (cross-check TS = 43). It is located at ( ' , b ) = (42 : 88 GLYPH<14> GLYPH<6> 0 : 08 GLYPH<14> , 0 : 69 GLYPH<14> GLYPH<6> 0 : 08 GLYPH<14> ), has a measured integral flux F ( E > 1 TeV) = 4 : 3 GLYPH<2> 10 GLYPH<0> 13 cm GLYPH<0> 2 s GLYPH<0> 1 , and an extension of 0 : 17 GLYPH<14> GLYPH<6> 0 : 04 GLYPH<14> . Two potential counterparts are found to be spatially coincident with this source candidate: the magnetar SGR 1900 + 14 (Mazets et al. 1979) and the SNR G42.8 + 0.6 (Fuerst et al. 1987). The former has an age GLYPH<28> c = 0 : 90 kyr and a spin-down luminosity ˙ E = 2 : 6 GLYPH<2> 10 34 erg s GLYPH<0> 1 . It is assumed to be at a distance 12 : 5 GLYPH<6> 1 : 7 kpc (Davies et al. 2009) based on an association of the magnetar with a massive star cluster (Wachter et al. 2008). SGR 1900 + 14 has similar properties to those of another magnetar, SGR 1806 GLYPH<0> 20, that is associated with the VHE source HESS J1808 GLYPH<0> 204 (Abdalla et al. 2016). It underwent a major burst of soft GLYPH<13> -ray emission in 1998 (Hurley et al. 1999; Frail et al. 1999) and, similar to SGR 1806 GLYPH<0> 20, it might also be emitting VHE GLYPH<13> -rays. Little is known about SNR G42.8 + 0.6. The centroid of the VHE emission is marginally coincident with the magnetar, while the bulk of the emission overlaps the northeastern half of the SNR shell.", '6. Summary and conclusions': "The H.E.S.S. Collaboration has completed its Galactic plane survey, which is an observation and analysis program that spanned over a decade. This paper presents the final results of the survey. The four-telescope H.E.S.S. Phase I array was used for the observations, which features a 5 GLYPH<14> FoV that is well suited to scanning large regions of the sky like the Galactic plane. The Phase I array has a typical sensitivity to point-like GLYPH<13> -ray sources of 1% Crab Nebula integral flux ( E > 1 TeV) in less than 25 h.  \nThe H.E.S.S. Collaboration added a fifth, larger telescope to the array in 2012 (H.E.S.S. Phase II) to extend its sensitivity to lower energies as well as its ability to rapidly reposition to observe transient phenomena. However, it also features a smaller FoV than the four Phase I telescopes, making it much less suited for scanning large regions. In addition, the HGPS had improved the uniformity of its exposure and achieved a target sensitivity of 2% Crab flux in the inner Galaxy. Primarily for these reasons, as well as the diminishing gains stemming from source significance scaling approximately as the square root of livetime, the H.E.S.S. Collaboration in 2013 decided not to continue the HGPS observation program.  \nFirst early results from the HGPS were published in 2005 (Aharonian et al. 2005b). The observations at the time amounted to just 120 h yet led to the detection of ten VHE GLYPH<13> -ray sources in the inner Galaxy, eight of which were not previously known. Further results followed in 2006 (Aharonian et al. 2006d), using 230 h of data and discovering additional four GLYPH<13> -ray sources. Since then, we provided the community with periodic updates,  \nwhich had steadily increasing exposure and additional source discoveries, by releasing new unidentified sources (Aharonian et al. 2008a) and via published conference proceedings (Hoppe 2008b; Chaves et al. 2008a; Chaves 2009; Gast et al. 2011; Deil 2012; Carrigan et al. 2013a,b).  \nThe HGPS data set (Sect. 2) is now over a factor of ten larger than in 2006, comprising 2673 h of observations accumulated over the period 2004 to 2013. These data come from a variety of observations: observations from the initially published surveys, targeted observations of known sources, follow-up observations of newly discovered source candidates, observations to extend the HGPS spatial coverage, and fill-up observations to achieve a more uniform sensitivity across the Galactic plane (Fig. 4). The energy threshold of the HGPS varies with the longitude observed but is typically lower than 0.8 TeV for detections and maps (0.5 TeV for spectral analyses) and as low as 0.4 TeV (0.2 TeV) in many regions, especially the innermost Galaxy (Fig. 2).  \nCompared to the previous publication, the HGPS was also expanded to cover a much wider range of both longitude and latitude (Fig. 1). In the first Galactic quadrant, the HGPS now extends in longitude from the Galactic center to nearly ' = 65 GLYPH<14> , the northern limit of visibility from the southern-hemisphere H.E.S.S. site. In the fourth Galactic quadrant, the HGPS coverage is continuous and even extends beyond Vela to ' = 250 GLYPH<14> in the third quadrant. In latitude, the coverage varies but is generally b = GLYPH<6> 3 GLYPH<14> and as large as b = GLYPH<6> 5 GLYPH<14> in some regions to explore areas of particular interest o GLYPH<11> the plane. The point-source sensitivity is better than 2% Crab along the Galactic plane ( b = 0 GLYPH<14> ) over most of the longitudes covered by the HGPS (Figs. 4, 11). However, the flux sensitivity varies significantly owing to the mix of observations comprising the HGPS. It is better than 1% Crab in numerous regions but at a more modest level of 210% Crab o GLYPH<11> -plane. The HGPS achieves the best sensitivity at the Galactic center, reaching 0.3% of the Crab flux.  \nTo ensure robust results, the HGPS relies on results that agree between two independent software frameworks (chains; Sect. 2.3) used to calibrate raw Cherenkov data as well as reconstruct and analyze the GLYPH<13> -ray images and spectra. The primary software chain used the Hillas method for event reconstruction, and an event classification method using boosted decision trees. The secondary (cross-check) chain uses an alternative event reconstruction and classification based on EAS models, returning results that are in very good agreement globally although there are some variations on a source-by-source basis (discussed in Sect. 5.4). Monte Carlo simulations provide the instrument response functions that describe the performance of the instrument. The mean angular resolution of H.E.S.S. (68% containment radius of the PSF) is GLYPH<24> 0 : 08 GLYPH<14> and varies by approximately 10% across the survey region.  \nWe have generated a number of sky maps (images; Sect. 3), which are public data products 13 and also form the basis for the HGPS source catalog construction. To accumulate su GLYPH<14> cient signal and search for GLYPH<13> -ray emission of di GLYPH<11> erent sizes, we generated three di GLYPH<11> erent sets of maps with events spatially correlated over radii of 0 : 1 GLYPH<14> (point-like), 0 : 2 GLYPH<14> and 0 : 4 GLYPH<14> , respectively. To subtract background from hadronic CRs passing GLYPH<13> -ray selections in the FoV, we developed an adaptive version of the classic ring background method that is more flexible and can compensate for large exclusion regions that minimize signal contaminating background regions (Sect. 3.2.3). For each point in a sky map, we calculated the GLYPH<13> -ray excess, statistical significance of the GLYPH<13> -ray excess, and GLYPH<13> -ray flux (or flux upper limit). The main map products released (Sect. 7.1) are those of significance (Sect. 3.3), flux (Fig. 34) and upper limits. Auxiliary maps include flux errors and sensitivity (Fig. 4).  \nTo detect and characterize the VHE GLYPH<13> -ray sources, we developed a semi-automatic analysis pipeline to construct a source catalog (see Sect. 4). To disentangle individual GLYPH<13> -ray sources in complex regions of overlapping emission, we implemented morphological modeling based on two-dimensional maximum-likelihood estimation. We fit the GLYPH<13> -ray excess by twodimensional symmetric Gaussian components, keeping components with TS > 30. To arrive at the HGPS catalog, Gaussian components that that did not correspond to a clear emission peak in the main and cross-check analysis were rejected. Some components that strongly overlapped were merged into a single source for which position, extension, and flux were characterized by the moments of the multi-Gaussian emission model. In this process, it was necessary to model the underlying large-scale GLYPH<13> -ray emission along the Galactic plane to improve the modeling of the discrete sources (Sect. 4.6). We chose to use an empirical model derived with a sliding window method, 20 GLYPH<14> wide in longitude and Gaussian in latitude, whose Gaussian center, amplitude, and width were fit to the excess outside exclusion regions. We calculated source spectra using the reflected region background method when possible, fitting PL spectral models and determining the best-fit normalization and spectral index. Flux information is also available from the aforementioned maps albeit assuming a spectral index of GLYPH<0> = 2 : 3.  \nThe HGPS source catalog includes 78 sources of VHE GLYPH<13> -rays. Of these, 64 were detected with the HGPS pipeline analysis. For completeness, the catalog includes an additional 14 H.E.S.S. sources from regions excluded from the HGPS pipeline, for example, because of their complexity, such as the Galactic center region and sources with shell-like morphologies. H.E.S.S. has previously published the discovery of most of the HGPS sources, although in many cases the available observation time used for the HGPS analysis is considerably larger. Of the total 78 sources, 16 are new  \ndiscoveries published here for the first time. Five of these new sources are firmly identified objects: HESS J1554 GLYPH<0> 550 and HESS J1849 GLYPH<0> 000 are PWNe (Sect. 5.6.5, 5.6.15), and HESS J1119 GLYPH<0> 614, HESS J1833 GLYPH<0> 105, and HESS J1846 GLYPH<0> 029 are composite SNRs (Sect. 5.6.1, 5.6.10, 5.6.13). Three more of the new sources are spatially coincident with HE GLYPH<13> -ray pulsars, recently discovered in Fermi -LAT data, and are thus plausible PWN candidates.  \nThe HGPS sources have diverse characteristics (Sect. 5). Apart from the shell-like sources, most source morphologies are generally well-modeled as symmetric two-dimensional Gaussians, but their sizes range from point-like ( < GLYPH<24> 0 : 1 GLYPH<14> ) to 0 : 6 GLYPH<14> (Fig. 13). Their fluxes cover a wide range as well from 0.6% Crab to 103% Crab of which the majority are in the range 1-20% Crab (Fig. 13). The cumulative log N - log S distribution above 10% Crab (containing 32 sources) is well described by a power law of slope GLYPH<0> 1 : 3 GLYPH<6> 0 : 2 (Fig. 15), matching the expectation of a power law of slope GLYPH<0> 1 from a population of equalluminosity sources homogeneously distributed in the Galactic disk. Below 10% Crab, the HGPS source catalog is incomplete and can only provide a lower limit on the true number of fainter VHE sources (70 above 1% Crab). Spectral indices range from hard ( GLYPH<0> GLYPH<25> 2 : 0) to very soft ( GLYPH<0> GLYPH<25> 3 : 0) in an approximately normal distribution centered at 2 : 4 GLYPH<6> 0 : 3 (Fig. 14). The VHE sources cluster narrowly along the Galactic plane (median b = GLYPH<0> 0 : 20 GLYPH<14> , with a spread of 0 : 51 GLYPH<14> ), in good agreement with the distributions of SNRs, energetic pulsars, molecular gas, and HE GLYPH<13> -ray sources (Fig. 11). Their distribution in longitude (Fig. 12) shows a general correlation with molecular gas.  \nTo study the origin of the VHE GLYPH<13> -rays, we performed a systematic search to associate the HGPS sources with known or suspected VHE source classes, based largely on spatial compatibility with objects in the SNR and PWN catalog SNRcat, the ATNF pulsar catalog, and the Fermi -LAT 3FGL and 2FHL catalogs (Sect. 5.1). By comparing the HGPS catalog to plausible MWLcounterpart catalogs, we come to one of the main conclusions of the HGPS program: the majority (67 , or 86%) of the HGPS sources are associated with at least one astronomical object that could potentially account for the production of GLYPH<13> -rays at TeV energies. The unassociated sources (11 , 14%) are not necessarily dark, i.e. emitting exclusively in the VHE domain; it is also possible that their counterparts were missed by our association procedure. In short, most HGPS sources have either firm associations, plausible or potential counterparts in other wavelength regimes. Whether there remains a population of truly dark VHE sources in the HGPS can only be figured out with deeper MWLstudies.  \nWe then used additional, stricter criteria, such as shell-like morphology or variability, to establish firm identifications for 31 sources (Fig. 10). We found the largest identified VHE source class to be PWNe (12 sources, or 39% of identified sources),  \nfollowed by shell-type SNRs (8 , 26%); composite SNRs (8 , 26%), where both the interior PWN and SNR shell may contribute to the emission; and high-energy binary systems (3 , 10%). At present, only 40% of the HGPS sources can be firmly identified. This is typically due to di GLYPH<14> culties resolving ambiguity among competing scenarios involving multiple associated objects in large part because of the large intrinsic sizes of VHE GLYPH<13> -ray sources.  \nThe HGPS data set allows for population studies of sources. An early study of 15 globular clusters was published before the HGPS was completed (H.E.S.S. Collaboration et al. 2013). Two further such studies, on the primary Galactic VHE source classes of PWNe and SNRs, are published as companion articles to this paper (H.E.S.S. Collaboration et al. 2017f,b, respectively), together with more specific studies on a number of microquasars (H.E.S.S. Collaboration et al. 2016d) and bow shocks of runaway stars (H.E.S.S. Collaboration et al. 2017e). With the public release of the HGPS catalog along with sky maps, more comprehensive such population studies will become possible.  \nFurther insights into the Galactic VHE source population and di GLYPH<11> use emission in the coming years can be expected. H.E.S.S., Fermi -LAT and HAWC are surveying the Milky Way; the analysis methods for the individual gamma-ray data sets, and joint analysis methods combining multiple data sets are improving; and new surveys at lower wavelengths (especially those detecting nonthermal emission in the radio and X-ray bands) will be come available soon. The next major leap forward will be achieved by the Galactic plane survey of the Cherenkov Telescope Observatory (CTA), which will consist of two arrays in the northern and southern hemisphere (The Cherenkov Telescope Array Consortium et al. 2017). The Galactic plane survey is a key science project of CTA, and is planned to cover the whole Galactic plane, over a wider energy band and with better angular resolution and sensitivity compared to HGPS (Dubus et al. 2013; The Cherenkov Telescope Array Consortium et al. 2017).  \nIn conclusion, the additional exposure obtained since 2006, plus significant improvements in analysis and reconstruction methods, allowed us to probe much more of the Galaxy, whether it be more distant sources, fainter nearby sources, or regions never before observed at TeV energies. The HGPS program clearly demonstrates that sources of VHE GLYPH<13> -ray emission are common in the Galaxy and are linked to diverse sites of highenergy particle acceleration.", '7. Online material': 'In this section, we provide further information about the public data products released in electronic format. We also provide some guidance and caveats regarding the correct use of these products. The HGPS survey maps and source catalog presented in this paper are available for download at:', 'https://www.mpi-hd.mpg.de/hfm/HESS/hgps': 'In addition to the figures and tables present in Sect. 5, there are a series of HGPS maps and tables available online:  \n- -Figure 34: HGPS flux map\n- -Figures 35-38: Four-panel HGPS significance maps with all VHE sources and MWL associations labeled\n- -Table 10: HGPS catalog source morphology summary\n- -Table 11: HGPS catalog source spectrum summary\n- -Table 12: HGPS catalog source associations', 'Description': 'The HGPS source catalog (construction described in Sect. 4) and a number of other tables are available as BINTABLE FITS extensions in the hgps\\_catalog\\_v1.fits.gz file.  \nAn overview of the available tables (including links to the tables in this paper describing the columns in detail) is given in Table 4. Here is some further information on the content of the tables:  \n- -HGPS\\_Sources : The HGPS catalog, one source per row, identified via the Source\\_Name column, which is in the format HESS JHHMM GLYPH<6> DDd .\n- -HGPS\\_Gauss\\_Components : The HGPS Gaussian component list, one component per row. Reference back to HGPS\\_Sources catalog via the Source\\_Name column (if the component is part of a source).\n- -HGPS\\_Associations : The HGPS association list, one association per row. Reference back to HGPS\\_Sources catalog via the Source\\_Name column. A given HGPS catalog source can appear zero, one, or multiple times in this table.', 'Usage notes and caveats': 'For reasons of reproducibility, we decided to release the complete emission model on which the analysis is based. This includes the full list of Gaussian components and parameters of the large-scale emission model. When working with this data, beware of following usage notes and caveats:  \n- -Some of the components are unstable and are not confirmed by the cross-check analysis. Use the source catalog, not the component list, for studies based on the HGPS.\n- -For the HGPS catalog, we did not perform detailed persource systematic error estimates. In general, when using spectra from the HGPS catalog, we recommend assuming a systematic error of 30% on the absolute flux and 0.2 on the spectral index.\n- -In Fig. 9, there are a few sources where the integral flux estimate di GLYPH<11> ers by more than 30% when using the two methods discussed in this paper. As discussed in Sect. 4.12, the estimate of a source spectrum is a GLYPH<11> ected by the assumed source morphology, di GLYPH<11> use gamma and atmospheric hadronic background model and uncertainties in the instrument response functions. In particular, the integral flux estimate may be uncertain by more than 30% for sources with relatively low significance that are not spatially isolated from other sources. In those cases, one can assume the di GLYPH<11> erence between Flux\\_Map and Flux\\_Spec\\_Int\\_1TeV to be a lower limit on the systematic error.  \nTable 4. HGPS catalog FITS tables. See Sect. 7.2.  \nTable 5. HGPS FITS table columns for HGPS\\_Sources , the main source catalog. The column descriptions link back to sections and equations in the main text where needed.  \nTable 5. continued.  \nTable 5. continued.Table 6. HGPS FITS table columns for HGPS\\_Gauss\\_Components . See Sects. 4.8 and 4.10.  \nTable 7. HGPS FITS table columns for HGPS\\_Associations . See Sect. 5.1.  \nTable 8. HGPS FITS table columns for HGPS\\_Identifications . See Table 5.1.3.  \nTable 9. HGPS FITS table columns for HGPS\\_Large\\_Scale\\_Component See Sect. 4.6.', 'Acknowledgements': "This work made extensive use of gamma-cat 14 , SNRcat 15 (Ferrand & Safi-Harb 2012), ATNF 16 (Manchester et al. 2005), SIMBAD 17 (Wenger et al. 2000) and NASA's Astrophysics Data System Bibliographic Services. For data analysis, we made extensive use of the Python packages Gammapy 18 (Donath et al. 2015; Deil et al. 2017), Astropy 19 (Astropy Collaboration et al. 2013) and Sherpa 20 (Freeman et al. 2001), as well as Numpy (Dubois et al. 1996), Scipy (Jones et al. 2001) and Matplotlib (Hunter 2007).  \nThe support of the Namibian authorities and the University of Namibia in facilitating the construction and operation of H.E.S.S. is gratefully acknowledged, as is the support by the German Ministry for Education and Research (BMBF), the Max Planck Society, the German Research Foundation (DFG), the French Ministry for Research, the CNRS-IN2P3 and the Astroparticle Interdisciplinary Programme of the CNRS, the U.K. Science and Technology Facilities Council (STFC), the IPNP of the Charles University, the Czech Science Foundation, the Polish Ministry of Science and Higher Education, the South African Department of Science and Technology and National Research Foundation, the University of Namibia, the Innsbruck University, the Austrian Science Fund (FWF), and the Austrian Federal Ministry for Science, Research and Economy, and by the University of Adelaide and the Australian Research Council. We appreciate the excellent work of the technical support sta GLYPH<11> in Berlin, Durham, Hamburg, Heidelberg, Palaiseau, Paris, Saclay, and in Namibia in the construction and operation of the equipment. This work benefited from services provided by the H.E.S.S. Virtual Organisation, supported by the national resource providers of the EGI Federation.", 'List of Objects': "'HESS J1825 GLYPH<0> 137' on page 7  \n'LS 5039' on page 7  \n'HESS J1303 GLYPH<0> 631' on page 7  \n'PSR B1259 GLYPH<0> 63' on page 7  \n'RX J0852.0 GLYPH<0> 4622' on page 11  \n'Vela Junior' on page 11  \n'HESS J0852 GLYPH<0> 463' on page 12  \n'HESS J1442 GLYPH<0> 624' on page 12  \n'RCW 86' on page 12  \n'HESS J1534 GLYPH<0> 571' on page 12  \n'HESS J1614 GLYPH<0> 518' on page 12  \n'HESS J1713 GLYPH<0> 397' on page 12  \n'RX J1713.7 GLYPH<0> 3946' on page 12  \n'HESS J1731 GLYPH<0> 347' on page 12  \n'HESS J1912 + 101' on page 12  \n'HESS J1745 GLYPH<0> 290' on page 12  \n'HESS J1746 GLYPH<0> 285' on page 12  \n'HESS J1747 GLYPH<0> 281' on page 12  \n'HESS J1718 GLYPH<0> 374' on page 12  \n'HESS J1741 GLYPH<0> 302' on page 12  \n'HESS J1801 GLYPH<0> 233' on page 12  \n'W 28' on page 12  \n'HESS J1911 + 090' on page 12  \n'W 49B' on page 12  \n'HOTS J1111 GLYPH<0> 611' on page 20  \n'Puppis A' on page 20  \n'2FHL J0826.1 GLYPH<0> 4500' on page 20  \nArticle number, page 54 of 79", 'References': 'Abdalla, H., Abramowski, A., Aharonian, F., et al. 2016, ArXiv e-prints [ arXiv:1606.05404 ]  \nAbdo, A. A., Ackermann, M., Ajello, M., et al. 2009a, Science, 325, 840  \nAbdo, A. A., Ackermann, M., Ajello, M., et al. 2009b, ApJS, 183, 46  \nAbdo, A. A., Ajello, M., Allafort, A., et al. 2013, ApJS, 208, 17  \nAbeysekara, A. U., Albert, A., Alfaro, R., et al. 2017, ApJ, 843, 40  \nArticle number, page 56 of 79', 'H.E.S.S. Collaboration et al.: The H.E.S.S. Galactic plane survey': "de Naurois, M. & Rolland, L. 2009, Astroparticle Physics, 32, 231  \nDeil, C. 2012, in IAU Symposium, Vol. 284, The Spectral Energy Distribution of Galaxies - SED 2011, ed. R. J. Tu GLYPH<11> s & C. C. Popescu, 365-370  \nDeil, C., Zanin, R., Lefaucheur, J., et al. 2017, ArXiv e-prints [ arXiv:1709.01751 ]  \nDjannati-Ataï, A., de Jager, O. C., Terrier, R., Gallant, Y. A., & Hoppe, S. 2008, International Cosmic Ray Conference, 2, 823  \nDjannati-Ataï, A., Marandon, V., Chaves, R. C. G., et al. 2009, HESS discovery of VHE gamma-ray emission of a remarkable young composite SNR, Online presentations from the Workshop on Supernova Remnants and Pulsar Wind Nebulae in the Chandra Era, http: // cxc.harvard.edu / cdo / snr09 / program.html  \nDonath, A., Deil, C., Arribas, M. P., et al. 2015, ArXiv e-prints [ arXiv:1509.07408 ]  \nDrury, L. O., Aharonian, F. A., & Voelk, H. J. 1994, A&A, 287, 959  \nDubner, G. M., Holdaway, M., Goss, W. M., & Mirabel, I. F. 1998, AJ, 116, 1842  \nDubois, P. F., Hinsen, K., & Hugunin, J. 1996, Computers in Physics, 10  \nDubus, G., Contreras, J. L., Funk, S., et al. 2013, Astroparticle Physics, 43, 317  \nEger, P., Rowell, G., Kawamura, A., et al. 2011, A&A, 526, A82  \nFerrand, G. & Safi-Harb, S. 2012, Advances in Space Research, 49, 1313  \nFiasson, A., Marandon, V., Chaves, R. C. G., Tibolla, O., & for the H.E.S.S. Collaboration. 2009, International Cosmic Ray Conference  \nFomin, V. P., Stepanian, A. A., Lamb, R. C., et al. 1994, Astroparticle Physics, 2, 137  \nFrail, D. A., Cordes, J. M., Hankins, T. H., & Weisberg, J. M. 1991, ApJ, 382, 168  \nFrail, D. A., Kulkarni, S. R., & Bloom, J. S. 1999, Nature, 398, 127  \nFreeman, P., Doe, S., & Siemiginowska, A. 2001, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 4477, Society of PhotoOptical Instrumentation Engineers (SPIE) Conference Series, ed. J.-L. Starck & F. D. Murtagh, 76-87  \nFuerst, E., Reich, W., Reich, P., Handa, T., & Sofue, Y. 1987, A&AS, 69, 403  \nFunk, S., Hinton, J. A., Moriguchi, Y., et al. 2007, A&A, 470, 249  \nGabici, S., Aharonian, F. A., & Blasi, P. 2007, Ap&SS, 309, 365  \nGabici, S. & Montmerle, T. 2015, International Cosmic Ray Conference, 29  \nGaensler, B. M., Arons, J., Kaspi, V. M., et al. 2002, ApJ, 569, 878  \nGaensler, B. M. & Slane, P. O. 2006, ARA&A, 44, 17  \nGast, H., Brun, F., Carrigan, S., et al. 2011, International Cosmic Ray Conference, 7, 158  \nGavriil, F. P., Gonzalez, M. E., Gotthelf, E. V., et al. 2008, Science, 319, 1802  \nGonzalez, M. & Safi-Harb, S. 2003, ApJ, 591, L143  \nGotthelf, E. V. & Halpern, J. P. 2009, ApJ, 700, L158  \nGotthelf, E. V., Halpern, J. P., Terrier, R., & Mattana, F. 2011, ApJ, 729, L16  \nGotthelf, E. V., Tomsick, J. A., Halpern, J. P., et al. 2014, ApJ, 788, 155  \nGö˘gü¸s, E., Lin, L., Kaneko, Y., et al. 2016, ApJ, 829, L25  \nGreen, D. A. 2014, Bulletin of the Astronomical Society of India, 42, 47  \nGupta, Y., Mitra, D., Green, D. A., & Acharyya, A. 2005, Current Science, 89, 853  \nHahn, J., de los Reyes, R., Bernlöhr, K., et al. 2014, Astroparticle Physics, 54, 25  \nHelfand, D. J., Becker, R. H., White, R. L., Fallon, A., & Tuttle, S. 2006, AJ, 131, 2525  \nHelfand, D. J., Velusamy, T., Becker, R. H., & Lockman, F. J. 1989, ApJ, 341, 151  \nH.E.S.S. Collaboration, Abdalla, H., Abramowski, A., et al. 2017a, A search for new supernova remnant shells in the Galactic plane with H.E.S.S., A&A forthcoming  \nH.E.S.S. Collaboration, Abdalla, H., Abramowski, A., et al. 2017b, Galactic Supernova Remnants Population Study at Very High Gamma-Ray Energies with H.E.S.S., A&A forthcoming  \nH.E.S.S. Collaboration, Abdalla, H., Abramowski, A., et al. 2017c, HESS J1741-302: a hidden accelerator in the Galactic plane, A&A forthcoming  \nH.E.S.S. Collaboration, Abdalla, H., Abramowski, A., et al. 2017d, Observations of SS433 with the MAGIC and H.E.S.S. telescopes, A&A forthcoming  \nH.E.S.S. Collaboration, Abdalla, H., Abramowski, A., et al. 2016a, ArXiv e-prints [ arXiv:1611.01863 ]  \nH.E.S.S. Collaboration, Abdalla, H., Abramowski, A., et al. 2017e, ArXiv e-prints [  \narXiv:1705.02263  \n]  \nH.E.S.S. Collaboration, Abdalla, H., Abramowski, A., et al. 2016b, ArXiv e-prints [ arXiv:1609.00600 ]  \nH.E.S.S. Collaboration, Abdalla, H., Abramowski, A., et al. 2017f, ArXiv e-prints [ arXiv:1702.08280 ]  \nH.E.S.S. Collaboration, Abdalla, H., Abramowski, A., et al. 2016c, ArXiv e-prints [ arXiv:1609.08671 ]  \nH.E.S.S. Collaboration, Abdalla, H., Abramowski, A., et al. 2016d, ArXiv e-prints [ arXiv:1607.04613 ]  \nH.E.S.S. Collaboration, Abdalla, H., Abramowski, A., et al. 2017g, ArXiv e-prints [ arXiv:1706.04535 ]  \nH.E.S.S. Collaboration, Abramowski, A., Acero, F., et al. 2012, A&A, 548, A46  \nH.E.S.S. Collaboration, Abramowski, A., Acero, F., et al. 2013, A&A, 551, A26  \nH.E.S.S. Collaboration, Abramowski, A., Acero, F., et al. 2011a, A&A, 531, A81  \nH.E.S.S. Collaboration, Abramowski, A., Acero, F., et al. 2011b, A&A, 533, A103  \nH.E.S.S. Collaboration, Abramowski, A., Aharonian, F., et al. 2015a, A&A, 574, A27  \nH.E.S.S. Collaboration, Abramowski, A., Aharonian, F., et al. 2014a, A&A, 562, A40  \nH.E.S.S. Collaboration, Abramowski, A., Aharonian, F., et al. 2014b, A&A, 562, L4  \nH.E.S.S. Collaboration, Abramowski, A., Aharonian, F., et al. 2015b, A&A, 574, A100  \nH.E.S.S. Collaboration, Abramowski, A., Aharonian, F., et al. 2015c, A&A, 575, A81  \nH.E.S.S. Collaboration, Abramowski, A., Aharonian, F., et al. 2015d, Science, 347, 406  \nH.E.S.S. Collaboration, Abramowski, A., Aharonian, F., et al. 2015e, A&A, 577, A131  \nH.E.S.S. Collaboration, Abramowski, A., Aharonian, F., et al. 2016e, ArXiv e-prints [ arXiv:1601.04461 ]  \nH.E.S.S. Collaboration, Abramowski, A., Aharonian, F., et al. 2016f, Nature, 531, 476  \nHinton, J. A. & Hofmann, W. 2009, ARA&A, 47, 523  \nHobbs, G., Faulkner, A., Stairs, I. H., et al. 2004a, MNRAS, 352, 1439  \nHobbs, G., Lyne, A. G., Kramer, M., Martin, C. E., & Jordan, C. 2004b, MNRAS, 353, 1311  \nHofverberg, P., Chaves, R. C. G., Fiasson, A., et al. 2010, in 25th Texas Symposium on Relativistic Astrophysics, 196  \nHolder, J. 2009, ArXiv e-prints [ arXiv:0907.3918 ]  \nHoppe, S. 2008a, PhD thesis, Universität Heidelberg, Max-Planck-Institut für Kernphysik  \nHoppe, S. 2008b, International Cosmic Ray Conference, 2, 579  \nHunter, J. D. 2007, Computing In Science & Engineering, 9, 90  \nHurley, K., Cline, T., Mazets, E., et al. 1999, Nature, 397, 41  \nJohanson, A. K. & Kerton, C. R. 2009, AJ, 138, 1615  \nJohnston, S. & Galloway, D. 1999, MNRAS, 306, L50  \nJones, E., Oliphant, T., Peterson, P., et al. 2001, SciPy: Open source scientific tools for Python  \nKang, J.-H. & Koo, B.-C. 2007, ApJS, 173, 85  \nKang, J.-H., Koo, B.-C., & Salter, C. 2012, AJ, 143, 75  \nArticle number, page 60 of 79  \nReich, W. 1982, A&A, 106, 314  \nReimer, A., Pohl, M., & Reimer, O. 2006, ApJ, 644, 1118  \nRenaud, M., Goret, P., & Chaves, R. C. G. 2008, in American Institute of Physics Conference Series, Vol. 1085, American Institute of Physics Conference Series, ed. F. A. Aharonian, W. Hofmann, & F. Rieger, 281-284  \nRoberts, M. S. E., Gotthelf, E. V., Halpern, J. P., Brogan, C. L., & Ransom, S. M. 2007, in WE-Heraeus Seminar on Neutron Stars and Pulsars 40 years after the Discovery, ed. W. Becker & H. H. Huang, 24  \nRolke, W. A., López, A. M., & Conrad, J. 2005, Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 551, 493  \nRomero, G. E. 2010, Mem. Soc. Astron. Italiana, 81, 181  \nRusseil, D. 2003, A&A, 397, 133  \nSafi-Harb, S., Harrus, I. M., Petre, R., et al. 2001, ApJ, 561, 308  \nSafi-Harb, S. & Kumar, H. S. 2008, ApJ, 684, 532  \nShahinyan, K. & the VERITAS Collaboration. 2016, ArXiv e-prints [ arXiv:1610.05799 ]  \nSheidaei, F. 2011, International Cosmic Ray Conference, 7, 244  \nSlane, P., Bykov, A., Ellison, D. C., Dubner, G., & Castro, D. 2015, Space Sci. Rev., 188, 187  \nStewart, I. M. 2009, A&A, 495, 989  \nStil, J. M., Taylor, A. R., Dickey, J. M., et al. 2006, AJ, 132, 1158  \nStraal, S. M., Gabányi, K. É., van Leeuwen, J., et al. 2016, ApJ, 822, 117  \nStycz, K. 2016, PhD thesis, Humboldt-Universität Berlin  \nSu, Y., Chen, Y., Yang, J., et al. 2009, ApJ, 694, 376  \nSugizaki, M., Mitsuda, K., Kaneda, H., et al. 2001, ApJS, 134, 77  \nSun, M., Wang, Z.-r., & Chen, Y. 1999, ApJ, 511, 274  \nTam, P. H. T., Wagner, S. J., Tibolla, O., & Chaves, R. C. G. 2010, A&A, 518, A8  \nTanaka, S. J. & Takahara, F. 2011, ApJ, 741, 40  \nTemim, T., Slane, P., Arendt, R. G., & Dwek, E. 2012, ApJ, 745, 46  \nTemim, T., Slane, P., Castro, D., et al. 2013, ApJ, 768, 61  \nTemim, T., Slane, P., Gaensler, B. M., Hughes, J. P., & Van Der Swaluw, E. 2009, ApJ, 691, 895  \nTemim, T., Slane, P., Kolb, C., et al. 2015, ApJ, 808, 100  \nTerrier, R., Djannati-Atai, A., Hoppe, S., et al. 2008a, in American Institute of Physics Conference Series, Vol. 1085, American Institute of Physics Conference Series, ed. F. A. Aharonian, W. Hofmann, & F. Rieger, 316-319  \nTerrier, R., Mattana, F., Djannati-Atai, A., et al. 2008b, in American Institute of Physics Conference Series, Vol. 1085, American Institute of Physics Conference Series, ed. F. A. Aharonian, W. Hofmann, & F. Rieger, 312-315  \nThe Cherenkov Telescope Array Consortium, Acharya, B. S., Agudo, I., et al. 2017, ArXiv e-prints [ arXiv:1709.07997 ]  \nThe Fermi-LAT Collaboration. 2017, ArXiv e-prints [ arXiv:1702.00664 ]  \nTian, W. W. & Leahy, D. A. 2008, MNRAS, 391, L54  \nTorres, D. F., Cillis, A., Martín, J., & de Oña Wilhelmi, E. 2014, Journal of High Energy Astrophysics, 1, 31  \nUeno, M., Bamba, A., Koyama, K., & Ebisawa, K. 2003, ApJ, 588, 338  \nVallée, J. P. 2014, ApJS, 215, 1  \nVan Etten, A., Funk, S., & Hinton, J. 2009, ApJ, 707, 1717  \nArticle number, page 62 of 79  \nVasisht, G., Gotthelf, E. V., Torii, K., & Gaensler, B. M. 2000, ApJ, 542, L49  \nVelusamy, T. & Kundu, M. R. 1974, A&A, 32, 375  \nVoisin, F., Rowell, G., Burton, M. G., et al. 2016, MNRAS, 458, 2813  \nWachter, S., Ramirez-Ruiz, E., Dwarkadas, V. V., et al. 2008, Nature, 453, 626  \nWeinstein, A. 2009, ArXiv e-prints [ arXiv:0912.4492 ]  \nWeltevrede, P., Johnston, S., & Espinoza, C. M. 2011, MNRAS, 411, 1917  \nWenger, M., Ochsenbein, F., Egret, D., et al. 2000, A&AS, 143, 9  \nWhiteoak, J. B. Z. & Green, A. J. 1996, A&AS, 118, 329  \nWilks, S. S. 1938, The Annals of Mathematical Statistics, 9, 60  \nWilson, A. S. & Weiler, K. W. 1976, A&A, 53, 89  \nXu, J. & Zhang, H. 2009, ArXiv e-prints [ arXiv:0909.0394 ]  \nYounes, G., Archibald, R., Kouveliotou, C., et al. 2016, The Astronomer's Telegram, 9378  \nZajczyk, A., Gallant, Y. A., Slane, P., et al. 2012, A&A, 542, A12  \nZhou, P. & Chen, Y. 2011, ApJ, 743, 4  \n- 1 Centre for Space Research, North-West University, Potchefstroom 2520, South Africa\n- 2 Universität Hamburg, Institut für Experimentalphysik, Luruper Chaussee 149, D 22761 Hamburg, Germany\n- 3 Max-Planck-Institut für Kernphysik, P.O. Box 103980, D 69029 Heidelberg, Germany\n- 4 Dublin Institute for Advanced Studies, 31 Fitzwilliam Place, Dublin 2, Ireland\n- 5 National Academy of Sciences of the Republic of Armenia, Marshall Baghramian Avenue, 24, 0019 Yerevan, Republic of Armenia\n- 6 Yerevan Physics Institute, 2 Alikhanian Brothers St., 375036 Yerevan, Armenia\n- 7 Institut für Physik, Humboldt-Universität zu Berlin, Newtonstr. 15, D 12489 Berlin, Germany\n- 8 University of Namibia, Department of Physics, Private Bag 13301, Windhoek, Namibia\n- 9 GRAPPA, Anton Pannekoek Institute for Astronomy, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands\n- 10 Department of Physics and Electrical Engineering, Linnaeus University, 351 95 Växjö, Sweden\n- 11 Institut für Theoretische Physik, Lehrstuhl IV: Weltraum und Astrophysik, Ruhr-Universität Bochum, D 44780 Bochum, Germany\n- 12 GRAPPA, Anton Pannekoek Institute for Astronomy and Institute of High-Energy Physics, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands\n- 13 Institut für Astro- und Teilchenphysik, Leopold-Franzens-Universität Innsbruck, A-6020 Innsbruck, Austria\n- 14 School of Physical Sciences, University of Adelaide, Adelaide 5005, Australia\n- 15 LUTH, Observatoire de Paris, PSL Research University, CNRS, Université Paris Diderot, 5 Place Jules Janssen, 92190 Meudon, France\n- 16 Sorbonne Universités, UPMC Université Paris 06, Université Paris Diderot, Sorbonne Paris Cité, CNRS, Laboratoire de Physique Nucléaire et de Hautes Energies (LPNHE), 4 place Jussieu, F-75252, Paris Cedex 5, France\n- 17 Laboratoire Univers et Particules de Montpellier, Université Montpellier, CNRS / IN2P3, CC 72, Place Eugène Bataillon, F-34095 Montpellier Cedex 5, France\n- 18 IRFU, CEA, Université Paris-Saclay, F-91191 Gif-sur-Yvette, France\n- 19 Astronomical Observatory, The University of Warsaw, Al. Ujazdowskie 4, 00-478 Warsaw, Poland\n- 20 Aix Marseille Université, CNRS / IN2P3, CPPM, Marseille, France\n- 21 Instytut Fizyki Ja¸drowej PAN, ul. Radzikowskiego 152, 31-342 Kraków, Poland\n- 22 Funded by EU FP7 Marie Curie, grant agreement No. PIEF-GA-2012-332350,\n- 23 School of Physics, University of the Witwatersrand, 1 Jan Smuts Avenue, Braamfontein, Johannesburg, 2050 South Africa\n- 24 Laboratoire d'Annecy-le-Vieux de Physique des Particules, Université Savoie Mont-Blanc, CNRS / IN2P3, F-74941 Annecy-le-Vieux, France\n- 25 Landessternwarte, Universität Heidelberg, Königstuhl, D 69117 Heidelberg, Germany\n- 26 Université Bordeaux, CNRS / IN2P3, Centre d'Études Nucléaires de Bordeaux Gradignan, 33175 Gradignan, France\n- 27 Oskar Klein Centre, Department of Physics, Stockholm University, Albanova University Center, SE-10691 Stockholm, Sweden\n- 28 Wallenberg Academy Fellow,\n- 29 Institut für Astronomie und Astrophysik, Universität Tübingen, Sand 1, D 72076 Tübingen, Germany\n- 30 Laboratoire Leprince-Ringuet, Ecole Polytechnique, CNRS / IN2P3, F-91128 Palaiseau, France  \n- 31 APC, AstroParticule et Cosmologie, Université Paris Diderot, CNRS / IN2P3, CEA / Irfu, Observatoire de Paris, Sorbonne Paris Cité, 10, rue Alice Domon et Léonie Duquet, 75205 Paris Cedex 13, France\n- 32 Univ. Grenoble Alpes, CNRS, IPAG, F-38000 Grenoble, France\n- 33 Department of Physics and Astronomy, The University of Leicester, University Road, Leicester, LE1 7RH, United Kingdom\n- 34 Nicolaus Copernicus Astronomical Center, Polish Academy of Sciences, ul. Bartycka 18, 00-716 Warsaw, Poland\n- 35 Institut für Physik und Astronomie, Universität Potsdam, Karl-Liebknecht-Strasse 24 / 25, D 14476 Potsdam, Germany\n- 36 Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen Centre for Astroparticle Physics, Erwin-Rommel-Str. 1, D 91058 Erlangen, Germany\n- 37 DESY, D-15738 Zeuthen, Germany\n- 38 Obserwatorium Astronomiczne, Uniwersytet Jagiello'nski, ul. Orla 171, 30-244 Kraków, Poland\n- 39 Centre for Astronomy, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziadzka 5, 87-100 Torun, Poland\n- 40 Department of Physics, University of the Free State, PO Box 339, Bloemfontein 9300, South Africa\n- 41 Heisenberg Fellow (DFG), ITA Universität Heidelberg, Germany\n- 42 GRAPPA, Institute of High-Energy Physics, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands\n- 43 Department of Physics, Rikkyo University, 3-34-1 Nishi-Ikebukuro, Toshima-ku, Tokyo 171-8501, Japan\n- 44 Japan Aerpspace Exploration Agency (JAXA), Institute of Space and Astronautical Science (ISAS), 3-1-1 Yoshinodai, Chuo-ku, Sagamihara, Kanagawa 229-8510, Japan\n- 45 Now at The School of Physics, The University of New South Wales, Sydney, 2052, Australia\n- 46 Now at Instituto de Física de São Carlos, Universidade de São Paulo, Av. Trabalhador São-carlense, 400 - CEP 13566-590, São Carlos, SP, Brazil\n- 47 Now with the German Aerospace Center (DLR), Earth Observation Center (EOC), 82234 Wessling, Germany\n- 48 Now at Technische Universität Kaiserslautern, D-67663 Kaiserslautern, Germany\n- 49 Now at I. Physikalisches Institut B, RWTH Aachen University, 52056 Aachen, Germany\n- 50 Department of Physics and Astronomy, University of Manitoba, Winnipeg, MB R3T 2N2, Canada", 'A & A proofs: manuscript no. hgps\\_paper': "Akiyama, K., Stawarz, Ł., Tanaka, Y. T., et al. 2016, ApJ, 823, L26  \nAleksi'c, J., Alvarez, E. A., Antonelli, L. A., et al. 2012, A&A, 541, A13  \nAliu, E., Archambault, S., Aune, T., et al. 2014a, ApJ, 780, 168  \nAliu, E., Archambault, S., Aune, T., et al. 2014b, ApJ, 787, 166  \nAngüner, E. O., Aharonian, F., Bordas, P., et al. 2017, ArXiv e-prints [ arXiv:1701.07002 ]  \nAntonopoulou, D., Vasilopoulos, G., & Espinoza, C. M. 2016, The Astronomer's Telegram, 9282  \nAstropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33  \nBamba, A., Terada, Y., Hewitt, J., et al. 2016, ApJ, 818, 63  \nBamba, A., Yamazaki, R., Kohri, K., et al. 2009, ApJ, 691, 1854  \nBerge, D., Funk, S., & Hinton, J. 2007, A&A, 466, 1219  \nBernlöhr, K. 2008, Astroparticle Physics, 30, 149  \nBietenholz, M. F. & Bartel, N. 2008, MNRAS, 386, 1411  \nBlum, R. D., Damineli, A., & Conti, P. S. 1999, AJ, 117, 1392  \nBocchino, F., Parmar, A. N., Mereghetti, S., et al. 2001, A&A, 367, 629  \nBocchino, F., van der Swaluw, E., Chevalier, R., & Bandiera, R. 2005, A&A, 442, 539  \nBolz, O. 2004, PhD thesis, Ruprecht-Karls-Universität Heidelberg  \nBrogan, C. L., Gelfand, J. D., Gaensler, B. M., Kassim, N. E., & Lazio, T. J. W. 2006, ApJ, 639, L25  \nBurgay, M., Bailes, M., Bates, S. D., et al. 2013, MNRAS, 433, 259  \nCalabretta, M. R. & Greisen, E. W. 2002, A&A, 395, 1077  \nCamilo, F., Kaspi, V. M., Lyne, A. G., et al. 2000, ApJ, 541, 367  \nCamilo, F., Ransom, S. M., Gaensler, B. M., et al. 2006, ApJ, 637, 456  \nCarrigan, S., Brun, F., Chaves, R. C. G., et al. 2013a, ArXiv e-prints [ arXiv:1307.4690 ]  \nCarrigan, S., Brun, F., Chaves, R. C. G., et al. 2013b, ArXiv e-prints [ arXiv:1307.4868 ]  \nCash, W. 1979, ApJ, 228, 939  \nCaswell, J. L., McClure-Gri GLYPH<14> ths, N. M., & Cheung, M. C. M. 2004, MNRAS, 352, 1405  \nChaves, R. C. G. 2009, ArXiv e-prints [ arXiv:0907.0768 ]  \nChaves, R. C. G., de Oña Wilhemi, E., & Hoppe, S. 2008a, in American Institute of Physics Conference Series, Vol. 1085, American Institute of Physics Conference Series, ed. F. A. Aharonian, W. Hofmann, & F. Rieger, 219-222  \nChaves, R. C. G., Renaud, M., Lemoine-Goumard, M., & Goret, P. 2008b, in American Institute of Physics Conference Series, Vol. 1085, American Institute of Physics Conference Series, ed. F. A. Aharonian, W. Hofmann, & F. Rieger, 372-375  \nClifton, T. R. & Lyne, A. G. 1986, Nature, 320, 43  \nCordes, J. M. & Lazio, T. J. W. 2002, ArXiv Astrophysics e-prints [ astro-ph/0207156 ]  \nCui, Y., Pühlhofer, G., & Santangelo, A. 2016, A&A, 591, A68  \nDame, T. M., Hartmann, D., & Thaddeus, P. 2001, ApJ, 547, 792  \nDavies, B., Figer, D. F., Kudritzki, R.-P., et al. 2009, ApJ, 707, 844  \nde Jager, O. C. & Djannati-Ataï, A. 2009, in Astrophysics and Space Science Library, Vol. 357, Astrophysics and Space Science Library, ed. W. Becker, 451  \nde Jager, O. C., Harding, A. K., Baring, M. G., & Mastichiadis, A. 1995, International Cosmic Ray Conference, 2, 528  \nde los Reyes, R., Zajczyk, A., Chaves, R. C. G., & for the H.E.S.S. Collaboration. 2012, ArXiv e-prints [ arXiv:1205.0719 ]  \nArticle number, page 58 of 79", '270 275 Galactic Longitude (deg)': 'Galactic Latitude (deg) Galactic Latitude (deg) Galactic Latitude (deg) Galactic Latitude (deg)  \n-4 -3 -2 -1 +0 +1 +2 +3  \n-4 -3 -2 -1 +0 +1 +2 +3  \n-4 -3 -2 -1 +0 +1 +2 +3  \n-4 -3 -2 -1 +0 +1 +2 +3  \nFig. 34. Integral flux above 1 TeV using a correlation radius R c = 0 : 1 GLYPH<14> and assuming spectral index GLYPH<0> = 2 : 3, in units of % Crab. The map is only filled where the point-source sensitivity for a 5GLYPH<27> detection (c.f. Fig. 4) is better (lower) than 2.5% Crab. Article number, page 65 of 79  \n<!-- image -->  \nGalactic Longitude (deg)  \nFig. 35. HGPS sources and associations in context (1 of 4). The background image shows TS of the VHE GLYPH<13> -ray excess in the Galactic plane assuming a point source morphology. All HGPS catalog sources are shown on top with transparent circles that correspond to the measured size of the source. Source names are labeled above the Galactic plane. Associations for the sources are shown with markers in white and corresponding labels. Pulsar (PSR) associations from the ATNF catalog are shown with triangles; SNRs as white circles with the radius representing the size, PWN and Composite (COMP) associations from the SNRcat catalog are marked with squares; associations with Fermi -LAT 3FGL and 2FHL sources are shown with diamonds; extra associations are shown with hexagonal markers. 3FGL sources are not labeled when they are identical to the pulsar. 2FHL sources are not labeled, when they are identical to the 3FGL source. The various object categories and the association criteria applied are detailed in Sect. 5.1.  \n<!-- image -->  \np', 'Galactic Latitude (deg)': 'Fig. 36. HGPS in a MWL context (2 of 4). Fig. 35 continued.  \n<!-- image -->  \n<!-- image -->  \nFig. 37. HGPS in a MWL context (3 of 4). Fig. 35 continued.  \n<!-- image -->  \nArticle number, page 68 of 79  \n254  \nFig. 38. HGPS in a MWL context (4 of 4). Fig. 35 continued.  \n<!-- image -->  \n#14  \nArticle number, page 70 of 79  \nFig. 39. Illustration of cut-out regions (red), exclusion regions (blue) and modeling regions (ROIs; gray, numbered) in the HGPS map. The background image displays significance ( R c = 0 : 1 GLYPH<14> ) in inverse grayscale.  \n<!-- image -->  \np  \nTable 10. HGPS source catalog - summary of map-based measurements. This is a small excerpt of the information available in the FITS catalog. The data shown here correspond to the following catalog columns (described in Table 5): Source\\_Name , Spatial\\_Model , GLON , GLAT , Pos\\_Err\\_95 , Size , Size\\_Err , Size\\_UL , Flux\\_Map , Flux\\_Map\\_Err , and Sqrt\\_TS . The values for sources indicated with an asterisk are taken from external references; see Table 1 for details.  \np  \nTable 10. continued.  \nTable 11. HGPS source catalog - summary of spectral measurements. This is a small excerpt of the information available in the FITS catalog. The data shown here correspond to the following catalog columns (described in Table 5): Source\\_Name , RSpec , Energy\\_Range\\_Spec\\_Lo , Flux\\_Spec\\_Int\\_1TeV ( Flux\\_Map for extern sources), Index\\_Spec\\_PL , Index\\_Spec\\_PL\\_Err (for PL spectrum sources), Index\\_Spec\\_ECPL , Index\\_Spec\\_ECPL\\_Err , Lambda\\_Spec\\_ECPL , Lambda\\_Spec\\_ECPL\\_Err (for ECPL spectrum sources), Containment\\_RSpec , Contamination\\_RSpec and Flux\\_Correction\\_RSpec\\_To\\_Total . The values for sources indicated with an asterisk are taken from external references, see Table 1 for details.  \nTable 11. continued.  \nTable 12. HGPS source associations (see Sect. 5.1).  \nTable 12. continued.  \nTable 12. continued.  \nTable 12. continued.  \nTable 12. continued.'}