key: cord-0653925-c86qzisj authors: Winston, Elaine; Hora, Joseph; Tolls, Volker title: A Census of Star Formation in the Outer Galaxy II: The GLIMPSE360 Field date: 2020-06-04 journal: nan DOI: nan sha: 15681d7c393bb9de6af81b335fc7582398539563 doc_id: 653925 cord_uid: c86qzisj We have conducted a study of star formation in the outer Galaxy from 65degr$$180degr in Galactic longitude. Using the $DBSCAN$ method on the combined catalog, we identify 618 clusters or aggregations of YSOs having 5 or more members. We identify 10,476 Class I, 29,604 Class II, and 7,325 anemic Class II/Class III YSOs. The ratio of YSOs identified as members of clusters was 25,528/47,338, or 54%. We found 100 of the clusters identified have previously measured distances in the {it WISE} ion{H}{2} survey. We used these distances in our spectral energy distribution (SED) fitting of the YSOs in these clusters, of which 96 had YSOs with $<3sigma$ fits. We used the derived masses from the SED model fits to estimate the initial mass function (IMF) in the inner and outer Galaxy clusters: dividing the clusters by Galactocentric distances, the slopes were $Gamma = 1.87 pm 0.31$ above 3~M$_{odot}$ for $R_{Gal}<11.5$~kpc and $Gamma = 1.15 pm 0.24$ above 3~M$_{odot}$ for $R_{Gal}>11.5$~kpc. The slope of the combined IMF was found to be $Gamma = 1.92 pm 0.42$ above 3~M$_{odot}$. These values are consistent with each other within the uncertainties, and with literature values in the inner Galaxy high-mass star formation regions. The slopes are likely also consistent with a universal Salpeter IMF. The study of star formation has been revolutionized with the launch of the Spitzer Space Telescope (Spitzer; Werner et al. 2004 ). Much of the published work to date has focused on regions at Galactocentric radii less than the sun, and in nearby clouds. The Spitzer Legacy programs 1 , executed early in the mission, provided a large dataset for studies of star formation. The original GLIMPSE survey (Benjamin et al. 2003 ) and subsequent follow-on programs (Churchwell et al. 2009 ) mapped the inner Galactic plane (295 • < l < 65 • ) with IRAC at 3.6, 4.5, 5.8, and 8 µm, and the MIPSGAL programs (Carey et al. 2009 ) mapped the regions with MIPS at 24 and 70 µm. The c2d program (Evans et al. 2003 ) scanned large areas in nearby molecular clouds for low luminosity sources to obtain a sample of nearby solar-type stars for debris disk studies. The FEPS program (Meyer et al. 2004 ) studied a large sample of young nearby solartype stars to trace the evolution of circumstellar gas and dust from primordial planet-building stages in young circumstellar disks through to older collisionally generated debris disks. Later Legacy and other large programs further contributed to our knowledge of star formation in nearby molecular clouds. The Gould's Belt program (Allen et al. 2006) completed the observations of all prominent star-forming regions within 500 pc. Megeath et al. (2012 Megeath et al. ( , 2016 conducted studies of the Orion A and B clouds, identifying thousands of YSO candidates in this nearby massive star-forming region. Other programs mapped the massive star formation complexes such as Cygnus-X (Hora et al. 2009 ) and Vela-Carina (Majewski et al. 2007 ) at distances of ∼1 -2 kpc. These and many other individual programs have produced a wealth of data on nearby star-forming regions which have been utilized in thousands of papers and will continue to be mined for years to come. The outer Galaxy is a distinctly different environment to that of the inner Galaxy, with conditions seemingly less likely to efficiently form stars.(In this paper, we will use the term 'outer Galaxy' to refer to clusters with a Galactocentric radius greater than ∼8 kpc and a Galactic longitude between roughly 65 • < l < 265 • .) The efficiency with which a molecular cloud forms stars is thought to be dependent on its density, temperature, and chemical abundances (e.g. Evans 1999) . The metallicity of the Milky Way is believed to decline as a function of Galacto-centric radius (Rudolph et al. 1997) . Average temperatures in molecular clouds are found to be lower (Mead & Kutner 1988) , as is the cosmic ray flux (Bloemen et al 1984) . Further, the volume density of molecular clouds in the outer Galaxy is lower, and so interaction rates and incidence of spiral arm crossings will be lower compared to inner Galaxy regions over the star forming lifetime of an individual cloud. The Spitzer SMOG survey (Carey et al 2008) was designed to help fill in our knowledge of star formation in the outer Galaxy by providing deep coverage of a field in the outer Galaxy in the IRAC and MIPS bands. Utilizing data from this survey, we presented in our previous paper (Winston et al. 2019 , Paper I) an initial study of the outer regions of our Galaxy, where environmental factors may impact on the star formation occurring there. In this paper, we present a census of star formation across the entire plane of our outer Milky Way galaxy using Spitzer's GLIMPSE360 survey (Whitney et al. 2008; Whitney & GLIMPSE360 Team 2009 ). This survey allows us to identify the young stellar populations of these clusters to better answer the question of whether the colder, less dense, and lower metallicity environment of the outer Galaxy affects the formation and evolution of young stars. With this work we will expand on our previous study by applying the techniques outlined in Paper I to the outer Galaxy covered by the GLIMPSE360 Galactic plane survey from 65 • < l < 265 • , in a 3 • wide strip that follows the warp in the outer Galactic disk. Here we identify young stellar objects (YSOs) across this ∼600 deg 2 strip from their excess infrared emission, locate clusters of YSOs indicating new regions of star formation, determine the evolutionary class of the YSOs to analyze the protostellar ratio of the clusters, and make a preliminary assessment of the initial mass function (IMF) across clusters with known distances in the outer Galaxy. This paper is organized as follows. In Section 2 we discuss the origins of the infrared catalogs. In Section 3 we discuss contamination removal, and the identification of the YSOs and their evolutionary classification. We then discuss the spatial distribution of the young stars and the identification of stellar clusters in the outer Galactic fields in Section 4. We discuss the fits to the SEDs of YSOs in clusters with known distances in Section 5.1, and in Section 6 compare our results to previous YSO catalogs constructed for the outer Galaxy. Finally, a brief summary is presented in Section 7. The GLIMPSE360 survey completes the coverage of the Galactic plane that began with the earlier GLIMPSE surveys. The observations were taken as part of a Spitzer Warm Mission Exploration Science program 2 , and were performed using the two short-wavelength IRAC bands at 3.6 and 4.5 µm (Fazio et al. 2004 ). The GLIMPSE360 point source archive covering Galactic longitudes from 65 • < l < 265 • was downloaded from the Infrared Science Archive (IRSA). This catalog did not contain the SMOG field data which were previously reported in Paper I (102 • < l < 109 • ), and which were observed for the cryogenic Spitzer Mapping of the Outer Galaxy (SMOG) program. It also excludes the central field surrounding the Cygnus star forming complex (76 • < l < 82 • ), which forms part of the Cygnus-X survey, though the flanking fields in the Galactic plane are included in the archive. Two types of catalog are available: a highly reliable point source Catalog and a highly complete point source Archive. The Archive includes sources with spatial positions as close as 0. 5, while the Catalog excludes sources closer than 2 in position. As in Paper I, the more complete Archive was used in this paper, to aid in the detection of fainter, more embedded YSOs. The IRAC observations at 3.6 and 4.5 µm were obtained in High Dynamic Range mode, which was comprised of three visits per mosaic position with 0.6 and 12 second integrations, similar to the SMOG survey. These data products were produced using the GLIMPSE team's pipeline. The GLIMPSE360 catalog contained a total of 49,378,042 sources. The mid-IR photometry was supplemented by J, H and K-band photometry from the 2MASS point source catalog (Skrutskie et al. 2006) . The photometric catalogues were merged using a 1. 6 match-ing radius. Documentation describing the GLIMPSE360 survey and the reduction process in detail is available on the IRSA website 3 . The Wide-field Infrared Survey Explorer (WISE; Wright et al. 2010 ) provides mid-IR photometry at 3.5, 4.6, 12, and 22 µm. The AllWISE catalog is an all-sky survey combining the cryogenic WISE All Sky survey and the NEOWISE post-cryogenic survey (Mainzer et al. 2011 ). The catalog is available via the NASA/IPAC Infrared Science Archive (IRSA) archive. 4 Selecting all the sources in the GLIMPSE360 field was made, resulting in a regional catalog of 14,483,596 point sources. The WISE satellite has considerably lower spatial resolution when compared to Spitzer , with a highest resolution of ∼6. 1 at 3.5 µm. However, the astrometric accuracy and matching to the 2MASS catalog are to within 1. 5, and so we performed catalog matching to the GLIMPSE360 catalog at 1. 5. Young stellar objects are most frequently identified by their excess emission at IR wavelengths. This emission arises from reprocessed stellar radiation in the dusty material of their natal envelopes or circumstellar disks. The infrared identification of YSOs is carried out by identifying sources that possess colors indicative of IR excess and distinguishing them from reddened and/or cool stars (Winston et al. 2007; Allen et al. 2004; Gutermuth et al. 2004) . A full description of the criteria for identification and selection of YSO and non-YSO sources as applied to the datasets here is given in the Appendices. The following subsections outline the removal of background extragalactic objects, YSO selection methods for each dataset, the complete YSO catalog, and the evolutionary classification of the YSOs based on their excess IR emission. The GLIMPSE360 field covers the majority of plane of the outer Galaxy, where the background sources suffer negligible extinction from the Galactic bar, and thus it is expected that many of the point sources detected will be active galactic nuclei (AGN) or star-forming galaxies (polycyclic aromatic hydrocarbon (PAH) galaxies). Knots of emission in the structure of molecular clouds may also be mistaken for YSOs. A further source of confusion in the YSO sample comes from sources with photometric contamination of the apertures by PAH emission. Such sources, which we will refer to as contaminants, were identified in the SMOG field as outlined in Paper I. The same methods could not be applied to the GLIMPSE360 data, due to the lack of 5.8 and 8 µm photometry. In order to constrain the photometric and color cuts to be applied to the GLIMPSE360 data to remove contaminants, the SMOG field data were utilized and cuts were determined from the characteristics of the contaminants in that field. A description of the applied cuts is given in Appendix A. Sources in the color/magnitude spaces where the contaminants were located were removed before the selection of YSOs was carried out. Figure 1 shows the criteria for the removal of contaminants in the GLIMPSE360 data for sources in a section of the field in the Galactic longitude range 180 • < l < 185 • . It was necessary to show only a subsection in order to enhance the clarity of the plot. Young stellar objects were selected using a combination of color-color diagrams (CCDs) with IRAC and 2MASS+IRAC colors, as described in Appendix 3.4. Photometric uncertainties of < 0.2 mags and magni-tudes fainter than the saturation limit were required in all the bands used for a particular color-color diagram to select YSOs. The full GLIMPSE360 catalog contains 49,378,042 sources, of which 28,837 were identified as YSO candidates using the combined 2MASS and IRAC photometry. Hereafter we will refer to the YSO candidates as "YSOs", however a definitive classification would require a more detailed analysis of the spectra and other characteristics of each object. Figure 2 shows the three source selection color-color diagrams used for the identification of YSOs for sources in the section of the field from 180 • < l < 185 • . Remaining contaminants in the YSO sample that are not accounted for here include: galactic asymptotic giant branch (AGB) stars, and highly/unusually reddened field stars that may be confused for anemic Class II (Class IIa) or Class III objects. Such objects are expected to be scattered randomly over the field. This issue is discussed further in Section 4.3.3 on AGB contamination. The WISE catalog covering the GLIMPSE360 field comprised ∼29% the number of detections as the IRAC catalog. The long edges of the GLIMPSE360 field exhibit a sawtooth pattern due to the IRAC mapping procedure. In our extraction of the WISE data from the all-sky catalog, we used a simple Galactic latitude cutoff that matches the largest extent of the IRAC data, so the edges of the WISE source distribution are smooth compared to the sawtooth edge of the GLIMPSE360 data. For this reason, a number of WISE sources fall into regions that are not included in the IRAC catalog and are identified solely by the WISE photometry on the edges of the field. Contaminant removal was undertaken using a procedure similar to the one used in Paper I and as outlined in Appendix B. The resulting catalog contained 543,457 point sources after the removal of spurious sources and extragalactic contaminants. A search for YSOs was performed using the four band WISE and 2MASS photometry; a detailed discussion is given in Appendix B. Of the 543,457 point sources, a total of 20,892 were identified as YSOs. Paper I compared the photometry in the IRAC bands 1 & 2 and WISE bands 1 & 2 for sources matched to within 1. 5. It was found that the photometric magnitudes matched well to within the uncertainties, which include measurement errors and possible variability of the sources themselves. In the GLIMPSE360 field, we used photometry from the two available IRAC bands to replace the appropriate WISE band fluxes where the sources matched to within 1. 5, and the WISE selection criteria were applied again to locate YSOs that may have reliable longer wavelength detections in WISE, but due to the shallower survey coverage or source confusion may not have been reliably detected at the shorter wavelengths. In these cases, the 2MASS selection criteria were not applied, since they replicate the selection made using the GLIMPSE360 data directly. The combined IRAC+WISE source list contained 318,588 objects. A total of 11,196 YSOs were identified following our methodology, as outlined in Appendix C. Figure 3 shows the two source selection color-color diagrams used in the IRAC+WISE identification of YSOs for sources in the section of the field from 180 • < l < 185 • . The three sets of YSOs selections, 2MASS+IRAC (28,837 sources), IRAC+WISE (11, 196 sources) , and WISE (20,892 sources), were merged and common objects combined based on the source spatial position. Sources in different sets that were within 1. 0 of each other were assumed to be the same source and combined. The 2MASS+IRAC and IRAC+WISE catalogs also required a detection in the GLIMPSE360 catalog, which uniquely identifies each source. The unique list of YSOs in the GLIMPSE360 field contained 42,757 objects. The GLIMPSE360 YSOs were then combined with the SMOG field YSOs, to create the full catalog, which contained 47,405 candidate YSOs. We call this the Star Formation in the Outer Galaxy, or SFOG catalog, and we refer to the total field covered by all of the component surveys as the SFOG field. The SFOG catalog was then used for the following analysis. Table 1 lists a selection of the column identifiers of the photometry table for the full list of identified YSOs. This table is a stacked table of subsets from the standard IPAC tables available on the IPAC website, for the GLIMPSE360, SMOG, and WISE datasets. It is available in its entirety in the online version of the paper. The online SIMBAD catalog was searched for matches to the SFOG catalog within 2 , with 7,343 YSOs (∼17%) found to have a previous identification. Table 2 lists the object identifiers, positions, and by which selection method the YSO was identified. Table 3 lists the YSOs that were found to have matches in the SIM-BAD catalog along with their alternate identifications and object types. Figure 4 shows the spatial distribu- tion of the identified YSOs over the outer plane of the Milky Way. Young stars evolve through a number of broad stages from the embedded core phase, through the protostellar phase where stellar accretion is still dominant, to the circumstellar disk phase where the envelope has dissipated and processing of disk material is ongoing, to the weak disk regime where the disk has dissipated and planets will have formed. A general evolutionary classification of YSOs in the SFOG field was carried out by measuring the slope, α, of the spectral energy distribution (SED) across the mid-IR bandpasses (Lada & Wilking 1984) . In our SFOG catalog, many sources are detected with only the two shorter IRAC bands, while a smaller number also have the longer wavelength WISE data. The SED slope was calculated based on the available photometric bands longward of the 2MASS K-band inclusive, for each source by performing a least squares polynomial fit to the data. Protostellar objects (Class 0 and I) have a rising slope, α > 0; Class II sources are characterized by decreasing slopes between −1.6 < α < 0, while Class IIa sources lack optically thick emission from a disk and possess decreasing slopes α < −1.6, consistent with a weak emission above a stellar photosphere. Truly diskless Class III objects show slopes between −2.7 < α < −2.0 (Lada 1987; Lada et al. 2006) . Many previous publications using Spitzer data combine these two categories as Class III objects (e.g. Winston et al. 2019; Saral et al. 2017; Hora et al. 2009 ). As we do not have ancillary data, such as X-ray observations, to separate young completely diskless Class III YSOs from field stars, we cannot reliably differentiate between the weak disk bearing Class IIa YSOs identified here and truly mid-IR diskless Class III YSOs (Winston et al. 2011 . For this reason, we list the class of all objects with slope α < −1.6 as being Class IIa/III in this paper. The identified YSOs were assigned classes as follows: 10,476 Class 0/I, 29,604 Class II sources, and 7,325 weak emission Class IIa/III stars. Further, some of the WISE identified sources were classified from the color-color diagram as candidate transition disk sources. This clas-sification was not used in the subsequent analysis, with only the classification based on SED slope reported. Figure 5 shows the spatial distribution of the identified YSOs by evolutionary classification over the plane. The classification of each YSO is listed in Table 2 . The SFOG catalog is susceptible to the same difficulties in assessing incompleteness in the sample as the SMOG field, discussed in Paper I. In brief, the large spatial distribution and lack of data on the distances in the majority of the clusters means that estimates of the minimum mass YSOs detected are not possible. Incompleteness by evolutionary class is also difficult to quan- 19h50m06.311s +28d26m02.9688s 19h42m53.2805s +29d23m11.598s 1 0 0 0 -1 2 a A "1" in these YSO columns means the object was identified as a YSO based on this dataset (see Section 3). b The cluster number of the YSO is indicated in this column (see Table 4 ). A "-1" indicates that there was no cluster affiliation identified. c The numerical YSO class, as described in Section 3.4.1. Note- Table 2 is published in its entirety in the machine readable format. A subset of the column identifiers are shown here for guidance regarding its form and primary content. Table 3 is published in its entirety in the machine readable format. A portion is shown here for guidance regarding its form and content. tify: because of the wavelength range of the 2MASS and IRAC data from 1.2 to 4.5 µm, the survey is most sensitive to Class II pre-main sequence objects. The WISE long wavelength bands at 12 and 22 µm makes it more sensitive than GLIMPSE360 for the detection of embedded protostars, but it is hampered by lower resolution and sensitivity at the shorter wavelengths. A cursory visual examination of the spatial distribution of the YSOs over the ∼600 deg 2 SFOG field shows evidence of clustering/clumping, as can be seen in Figure 4. Following the method described in Paper I, we used the DBSCAN (Ester et al. 1996) density based algorithm to identify over-densities in the spatial distri-bution of the YSOs. The values of the two free parameters, (the scaling size for clustering) and M inP ts (the minimum number of points required to define a dense region) were determined following the approach used by Joncour et al. (2018) in the Taurus region. Figure 6 (a) shows the one point correlation function, which gives the ratio of the cumulative distributions of the identified YSOs and a random distribution over the same field. From this, the value of = 0. • 1 was selected. Figure 6 (b) examines the cumulative distribution of three nearest neighbour distributions at 8, 9, 10 th nearest neighbors for the random distribution showing that a probability of 0.001 occurs at 0. • 1 for a minimum cluster size with 9 members. With these values, 621 clusters were identified in the SFOG field. Figure 7 shows the spatial distribution of the identified clusters of YSOs over the plane. The black YSOs are unclustered, with the clustered YSOs colorcoded by identified cluster. The clustered YSOs represented 54% of the whole catalog, with 21,810 YSOs not clustered. The minimum cluster size was 5 members, the largest cluster identified contained 1,177 members. The median cluster size is 17 members. Of the 621 clusters, 133 have ten members or less, 25 clusters have 100-200 members, and 22 cluster have more than 200 identified members. There is a possibility that some of the smaller clusters identified are not genuine stellar groupings, but chance over-densities in the field. To determine the statistical likelihood that the identified groupings are true clusters, we repeated the DBSCAN analysis for 1000 iterations of randomly distributed points covering a similar area as the SFOG field. Using the same values of = 0. • 1 and M inP ts = 9, we varied the number of fake 'stars' as follows: 21810, 29604, 36000, 43000, 47405. These are the numbers of unclustered YSOs, Class II YSOs, two intermediate values, and the total number of YSOs in the SFOG catalog. This range was used to take into account the known clustering of objects, which reduces the overall density in the field, and thus reduces the likelihood of a 'fake' cluster being detected. The average number of fake clusters detected for each run was: 0.04, 0.43, 1.9, 6.9, 13.8. We can therefore assume that between 0-14 of the 621 clusters found by DBSCAN may be misidentified random alignments. The 'fake' clusters in our random trials all have ≤10 members, and so if any clusters in our SFOG list in Table 4 are chance overdensities rather than true clusters, they would most likely be in the subset of 133 clusters with ≤10 YSOs. There is the further possibility that the unclustered YSOs are contaminants; foreground or background Galactic field dwarfs or AGB stars with either a high extinction or a small amount of dusty material surrounding them that leads to an excess of flux in the mid-IR or extragalactic sources. Spectroscopic analysis would be necessary to secure the identification of all sources. The approximate area of each cluster was quantified by measuring the convex hull of the associated cluster members. The convex hull is the set of points whose vertices include all of the points in the set. Figures 9 and 10 show eight examples of clusters identified across the SFOG field, showcasing the range of sizes and environments over which clusters were identified. In many cases it is clear that the clusters form part of a larger association or star forming complex. They are often surrounded by non-clustered YSOs that are likely to be associated with the cluster but did not satisfy the DB-SCAN criteria. We present these clusters as starting points for future studies. The cluster identification of each YSO is listed in Table 2. Table 4 lists the properties of the largest clusters (number of YSOs > 100) found in the SFOG catalog, sorted by decreasing number of YSOs in the cluster. The table gives the cluster number and name, the number of YSOs, coordinates of the cluster central point, Note-The clusters are shown in order of decreasing number of YSO members. Only the clusters with more than 100 YSOs are shown here, Table 4 is published in the electronic version in its entirety in the machine readable format. the circular radius based on the separation of the most distant YSOs, and association with previously identified star forming regions. The complete version of the table is available in the electronic version of this paper. The electronic version also includes a list of WISE H II counterparts, and a full listing of the SIMBAD objects located within the convex hull of each cluster. These sources are not assumed to be physically associated with the cluster, and no attempt has been made to filter the lists or to match them to the YSOs. They are provided as a reference for more in-depth studies. Images of each of the individual clusters were constructed using mosaics with a radius of 5 times the estimated convex hull radius of each cluster. The IRAC mosaics were downloaded from the IPAC servers, and the WISE Coadder 5 was used to generate mosaics of the fields in the WISE bands. For cases where a cluster was near the edge of the IPAC mosaics of the IRAC data, adjoining mosaics were combined to form images that could be used to display the full cluster. These mosaics were then used to create a selection of 3-color images and a WISE 12 µm band grayscale image overlaid with the cluster convex hull and the locations of the YSOs by evolutionary class. The full set of mosaics for every cluster are available on the Harvard Dataverse SFOG page (Winston et al. 2020) . Figure 8 shows the spatial distribution of the identified clusters overlaid on a schematic view of the Milky Way Galaxy (Hurt 2008) . As can be seen from Table 4 , all of the larger clusters are associated with previously known star forming regions, identified for example through surveys of H II regions, CO emission, radio continuum emission, detection of dark clouds, or the infrared emission from compact sources or the surrounding nebula. However, we also find 6 of the smaller clusters do not have any SIM-BAD sources within their convex hulls and may be newly identified. The true level of contaminants remaining in the catalog used to select YSOs is difficult to estimate precisely 5 https://irsa.ipac.caltech.edu/applications/ICORE given the broad range in spatial coverage and distances covered by the SFOG field. Spectral typing of the YSOs would provide confirmation of their nature, and a program to obtain spectra and X-ray observations of a selection of the clusters is currently underway. We discuss in the following sections possible sources of contamination in the YSO catalog and their effects on the clustering analysis. To examine the effectiveness of the contamination cuts used in the GLIMPSE360 analysis, a comparison was made to the SMOG field of Paper I. By applying the same contamination cuts and selection criteria used in this paper, 1,512 YSOs were selected in the SMOG field. Of these, 1,102 objects were matched at 1 between the two catalogs. A majority of the non-matching sources were included in the SFOG catalog due to the differences between the selection criteria required because of the lack of IRAC 5.8 and 8 µm and MIPS 24 µm coverage in the GLIMPSE360 field. Plotting these sources on an IRAC four-color CCD, 22 objects exhibited colors similar to those of galaxies. This would imply that 22/1,512 or <2% of the candidate YSOs are likely to be extragalactic contaminants. By applying this percentage to the 28,837 objects in the GLIMPSE360 component of our catalog, we would estimate that 420 of these may possibly be contaminants. These would likely be found predominantly in the non-clustered population since both foreground and background contaminants tend to be randomly distributed across the field. In order to assess the validity of the clustering algorithm, each of the regions was visually checked to look for interesting or spurious clusters. From this assessment, a total of three of the identified clusters, #515, #516, and #564, were found to be invalid due to spurious WISE sources. The WISE photometry contaminant removal process removes confirmed diffraction spike ob-jects, but those objects with tentative identification as photometric diffraction spikes are not removed. In each of these cases, a visual inspection of the images containing the 'cluster' found them to contain a number of diffraction spike objects surrounding a saturated bright foreground star. In removing these three spurious clusters, the number of clusters is reduced to 618 in total. The number of YSOs is reduced from 47,405 to 47,338 YSOs. The breakdown by evolutionary class becomes: 10,461 Class I, 29,552 Class II, and 7,325 Class IIa/III YSOs. The Besançon Galactic population synthesis models were used to estimate the AGB population at four points along the Galactic plane (Robin et al. 2003 (Robin et al. , 2014 . The points chosen were at l = 70 • , 105 • , 180 • , 250 • and b = 1 • , 1 • , 0 • , −1 • , respectively, with a 1 deg 2 field at each location. The models predicted that these fields contained either one or two AGB stars each. By scaling these numbers up to the approximate 600 sq. deg size of the SFOG field, we estimate that there may be between 600-1,200 AGB stars in the SFOG field. This would imply that ∼1,200/47,338 or ∼2.5% of the candidate YSOs may be AGB stars. The AGB stars would tend to be randomly distributed over the field, and thus less likely to be included in a cluster, so these contaminants are not likely to affect our clustering analysis significantly. The Python SEDFitter 6 package of Robitaille et al. (2007) was used to provide an estimate of the mass, age, disk and accretion properties of the YSOs in the SFOG field. The code uses a sample grid of YSO model SEDs with varying age, mass, inclination etc., to compare to the input photometry, with the scale factor S (dependent on source distance and luminosity) and the extinction A V as free parameters. The code returns a sample of best fit models and their associated parameters. From our previous SMOG study, it was determined that within each cluster the range in distances for the best fit models covered the entire provided distance range, and that therefore the results could not be used in any meaningful way when the distance to the YSO is unknown. Therefore only those clusters with reliable distance estimates from WISE H II regions (Anderson et al. 2014) located within the cluster's convex hull were used in the SEDFitter modeling. Reliable distances were found for 100 of the 618 clusters identified in the SFOG field. The H II regions associated with each cluster and their distances are listed in the full online version of Table 4 . The range of Galactocentric distances for these 100 clusters ranged from 7.8 kpc to 18.1 kpc. The SEDFitter routine was run using a fixed distance range based on the distance to the cluster, allowing the A V to vary from 0-40 A V , for each cluster separately. YSOs lacking photometry across a sufficient number of bands were not fit successfully. Further, of those YSOs fit, not all had fits with low χ 2 values. Of the 100 clusters run, 96 contained some YSOs with 'good' fits (χ 2 < 3.0). In general, four photometric points were required for a 'good' fit to the SED model. In total, 6,234 YSOs had reliable model fits across the 96 clusters. Figure 11 shows three examples of SEDFitter model fits to YSOs in the SFOG catalog. Table 5 lists a sample of the results of the SEDFitter routine for the YSOs in the 96 clusters for which reliable fits were obtained. A number of parameters are presented for each model fit, including the best fit to the object mass, disk mass, age, A V , central temperature, disk accretion rate, etc., and the χ 2 value. The weighted average values of all parameters for fits with χ 2 < 3 were calculated and are presented for each source. The upper and lower limit model fit parameters are also supplied in each case. The full data table, including all columns, is available online in electronic format. Given the uncertainties in the distances and the sparse photometry for each source, we do not consider the individual ages and masses derived from the model fits to be entirely reliable and will not discuss them further here. However, cumulatively for all clusters, they can provide an insight into the relative ages and masses of the YSOs in different regions of the outer Galaxy. Fig-315°00 Figure 9 . Four examples of clusters identified in the SFOG field: #147, #217, #228, and #247 from top to bottom. In the left column are 3-color images in WISE 12 µm (red), IRAC 4.5 µm (green), and IRAC 3.6 µm (blue); the center column contains 3-color images with WISE 22 µm (red), WISE 12 µm (green), and IRAC 4.5 µm (blue); and the right column shows the WISE 12 µm in reverse grayscale with the identified YSOs and the calculated convex hulls for each cluster overlaid. The symbols show the positions of the Class I (red circles), Class II (green squares), and Class IIa/III (cyan pentagons) YSOs. All of the cluster images and associated FITS files are available from Winston et al. (2020) . Best fit 2 = 1.861 AV = 6.8 Scale = 0.18 Figure 11 . Examples of three SEDFitter model fits to YSOs in the SFOG catalog with detections in 2MASS, IRAC, and WISE: SRC28966, SRC34028, and SRC37530 in clusters 15, 85, and 106, respectively. Note- Table 5 is published in its entirety in the machine readable format. A portion is shown here for guidance regarding its form and content. ure 12 shows the age and relative mass of the YSOs in three clusters by their spatial distribution. The size of the circles indicates the mass of the YSO relative to the most massive object that we identified in that cluster. The color indicates the age, the range in ages in Myr is shown by the colorbar for each plot. There were no strong trends in the distribution of age or mass across the 618 clusters. The initial mass function (IMF) of the clusters in the outer Galaxy is of great importance to determine whether the environment of the outer Galaxy has had an affect on the star formation efficiency and rate. Because for each cluster we have the SED-derived masses for only a small sample of YSOs, we examined the IMF of all 96 clusters combined, as shown in Figure 13 (left). We then split the clusters into two groups based on the mass of the most massive identified member, with the cut at 10M , as shown in Figure 13 (right). There were 58 low mass clusters (2,929 YSOs), and 38 high mass clusters (3,305 YSOs). The clusters were then split into two groups based on their Galactocentric radius, with the division between inner and outer Galaxy placed at 11.5 kpc, based on Huang et al. (2015) . There were 46 clusters with 4,533 YSOs with R Gal < 11.5 kpc, and 50 clusters with 1,701 YSOs with R Gal > 11.5 kpc. In the figures, the red line plots the power law slope with Γ = 1.35, roughly the value of the Salpeter slope (Salpeter 1955 ). The blue line shows the Miller & Scalo (1979) broken power law fit with slopes of Γ 1 = 1.7 and Γ 2 = 2.3. The Kroupa (2001) IMF with three distinct values of Γ for the low, solar, and high mass regions is shown in yellow, as presented by Weisz et al. (2015) who used a slope of Γ = 1.45 above two solar masses. Linear regression fits were made to the high-mass end of the resultant histograms in the range 3 M < M * < 10 M . The slopes of the fits are listed in Table 6 . The slopes for the complete sample of YSOs, and the nearby and low mass clusters, are broadly consistent with the values of Γ presented in the literature by Parravano et al. (2018) with Γ = 1.7 − 2.1, as discussed in detail in Paper I. The high mass and more distant clusters have shallower slopes than the others, though still consistent with the Salpeter value. The 1-σ uncertainties overlap for the high/low mass cluster slopes, though they do not quite overlap for the near/far cluster slopes. The near/far distance appears to be due to the Figure 12 . Examples of the SEDFitter results for spatial distributions of ages and masses within each cluster. The relative sizes of each symbol indicate how massive the YSO is with respect to the most massive object in that cluster. The color bar shows how the color scaling relates to the calculated age of the YSO (in Myr), with the bluer objects being younger, and the redder objects being older. Figure 13 . Examination of the Initial Mass Function for the outer Galactic regions identified in the SFOG field. Left: The IMF as determined by combining the SEDFitter calculated masses for the 96 clusters with distance estimates and with members with χ 2 < 3 fits. The green histogram shows the best fit model masses for all the YSO. The data show reasonable correlation to both the Salpeter slope of the IMF (m∼-1.35) and more recent estimates (m∼-2.7, m∼-2.3) for a broken power law fit. Right: The IMF split into the 58 clusters with a highest mass member lower than 10 M and the 38 clusters with a highest mass member greater than 10 M . 1.275 0.345 R Gal < 11.5 kpc 1.873 0.308 R Gal > 11.5 kpc 1.146 0.241 difference in the modeled populations of the high and low mass clusters, and not necessarily to a difference in the environment of the outer Galaxy. To further investigate the comparison of the SEDfitter generated IMFs of the inner and outer galaxies, we undertook a comparison of the Cygnus-X Legacy Survey catalog of IRAC and MIPS identified YSOs (Winston et al., in preparation) . This catalog contains a total of 30,646 YSOs -2,029 Class I, 27,672 Class II, and 945 Class IIa/III young stars -across the Cygnus X North and South fields. These YSOs were run through the SEDFitter routine in the same way as the SFOG Figure 14 . Examination of the Initial Mass Function for the outer Galactic regions identified in the SFOG field. Left: The IMF for the 46 clusters within a Galactocentric radius of 11.5 kpc. Right: The IMF for the 50 clusters at Galactocentric radius greater than 11.5 kpc. data, with a fixed distance of 1.4 kpc, and the IMF constructed. The linear regression fit to the data between 3 M < M * < 10 M was Γ = 2.315 ± 0.298. This value is consistent with the literature values and with the slope measured for all the SFOG clusters and those with M axM < 10 or R Gal < 11.5 kpc. The IMF is generally presumed to be universal, however the precise origins of the IMF, its relation to the Core Mass Function (CMF), and the effects of galactic and local star-forming environment on its evolution are not fully understood. Paladini et al. (2019) provide a useful overview of our current understanding in their recent white paper. An accurate and in-depth study of the IMF in the outer Galaxy is beyond the scope of this paper. The results presented here should be taken as an indication that the IMF in the Outer Galaxy is similar to that of the Inner Galaxy, and that this may indicate that metallicity does not greatly impact the IMF. However, there are a number of caveats, both astrophysical and analytical, which affect the SEDFitter derived IMF that we will now outline. Variations in the IMF of an astrophysical origin could be due to e.g. metallicity (e.g. Kroupa 2020), feedback from massive stars in the cluster either triggering or impeding star formation (e.g. Walch et al. 2013; Krumholz et al. 2016) , jets and outflows (e.g. Cunningham et al. 2018) , disk evolution (e.g. Povich et al. 2016) , and clustered vs. distributed star formation (e.g. Bonnell et al. 2011) . Recent theory suggests that massive stars may form a few Myrs after initial low mass star formation begins, which would also lead to a variation in IMF with cluster age (Vázquez-Semadeni et al. 2019) . Faster disk evolution in intermediate and massive stars has also been reported by Povich et al. (2016) in the M17 SWex star forming region that would also bias the slope of the IMF in this mass regime. The limitations of the data can impact the accuracy of the SEDFitter results. The distance estimates to the clusters are derived from kinematic distances, which can have uncertainties of at least 10-20% and are subject to systematic uncertainties (Anderson et al. 2014 ). Cantat-Gaudin & Anders (2020) find in their Gaia DR2 study of open clusters that not all previously identified members are physically associated with the regions, this would also affect the SEDFitter results, which assume a single distance to each cluster. Further, each cluster also has a spatial 'depth' and so a YSO at the near edge of the region may lie a few hundred parsecs nearer than one at the far edge. The clusters are all of different ages, and within each cluster the subclusters may also have a range in age. Thus the disk fraction, and hence percentage of the population detected with Spitzer, will vary -this should be accounted for when stacking the clusters in the IMF, but cannot be done here. There are also the issues of varying levels of diffuse mid-IR emission and source crowding in both Spitzer and to a greater extent the WISE data. Diffuse emission reduces sensitivity in the IRAC bands, meaning fainter sources will not be detected. Source crowding is primarily a spatial resolution issue and is most prevalent in the central cores of clusters. This leads to an underestimation of the disk bearing population, particularly of the lower mass YSOs. The mid-IR data also does not provide a measure of the diskless population of the clusters. These populations are generally identified using spectroscopy or Xray observations. The coverage of X-ray studies is highly constrained in comparison to IR surveys and is generally focused on the centers of known star-forming regions (Winston et al. 2009 (Winston et al. , 2011 Broos et al. 2013) . These studies show large populations of YSOs without disks, and trends in disk fractions with proximity to OB stars. Our selection method is not very sensitive to massive stars, and may not detect all the of known massive stars in a region. This would lead to a steeper IMF slope at higher masses. A comparison with Skiff (2014) shows that of the 3302 objects in their catalog with spectral types 'O' or 'B' that were located within the convex hulls of the clusters, 591 were matched to a YSO in the SFOG catalog within 2. 0, or 18% of the known massive population. The accuracy of the SEDFitter analysis is limited by the small number of near-and mid-IR bands available for model fitting and the reddening across those bands. Povich et al. (2019) find a degeneracy between stellar effective temperature and extinction that can lead to inaccurately modeled masses. The presence of disk material adds to this issue, especially as disk/envelope masses are poorly constrained without far-IR and longer photometric data. The SFOG catalog of YSOs and identified clusters was compared to a number of other published surveys and databases covering the outer Milky Way Galaxy. The individual YSO candidates were matched to the online SIMBAD database to within a 2 radius for all 47,338 YSOs, with 7,343 matches. Of these, 4,428 had been previously identified as a type of young stellar object, leaving 42,910 possibly new YSO candidates (some previous studies may not be listed on the SIMBAD database). Table 3 provides a sample listing of the YSOs with matches within 2 of a SIMBAD source, giving the source identification and the object type. The full table is available online in electronic format, and includes further selected information pertaining to the SIMBAD objects. In their recent paper, Armentrout et al. (2020) list 166 H II regions previously thought to be radio quiet, that they confirm to be weak radio sources. Of the 166, 37 were found to lie within the convex hull of one of the 618 clusters in the SFOG field. Twenty of these matches were found between 90 • and 115 • Galactic longitude. There was no significant difference in the median number of YSOs or the median effective hull radius between the clusters in the radio quiet Armentrout sample or the WISE H II sample. Further, when matching the H II regions to the catalog of YSOs, it was found that 58 had a YSO located within 10 and that 107 had a YSO located within 1 . The positions of the clusters were also compared to those of Avedisova (2002) , who reported 66,887 clusters over the whole Galactic plane. Of these, 24,101 were within the boundaries of the SFOG field. Of the 618 clusters we found in the SFOG field, 260 were matched to one or more of the Avedisova (2002) clusters, while 358 were new to this survey. The Gaia DR2 data release was also searched for matches to the SFOG catalog. Of the 47,338 YSOs identified by SFOG, 25,919 had a Gaia counterpart within 0. 5. These matched sources were then compared to the Bailer-Jones et al. (2018) A comparison was then made to the 100 identified clusters with known distances to determine if the distances of the YSOs in those clusters corresponded to the WISE H II catalog distances. The clustered YSOs showed both a large range in individual uncertainty of a few kiloparsecs and a spread within the cluster YSOs of a few kiloparsecs. While some individual YSOs had distances corresponding to that from WISE, the median and average values for the cluster did not. It can be noted that the Gaia distance uncertainties increase greatly beyond ∼0.5 kpc, and that this may account for the lack of consistency in cluster distance estimates. Given the lack of consistency in the cluster distances, it was decided not to use these in the SEDFitter routine or to report them in the catalog. Cantat-Gaudin & Anders (2020) calculated Gaia DR2 distances to previously identified open clusters, using the mode of distance likelihoods of known members. Of these, 27 matched to an SFOG cluster and 12 had distances in the WISE HII catalog. The differences in the measured distances varied widely (from 250-7300 pc). Given that young clusters are embedded and Gaia is more likely to detect foreground stars or the front cluster population, we decided to use the kinematic distances in this paper. The first comparison was made to the SMOG field of Paper I, to assess the effect of the loss of the 5.8 and 8 µm IRAC bands on the detection rate of YSOs. The SMOG field 2MASS+IRAC selection technique identified 3835 YSOs. The same GLIMPSE360 criteria used in this paper identified 1,512 YSOs. Of these, 1,102 objects matched between the two catalogs. This suggests that we are detecting 29% of the original catalog with the GLIMPSE360 criteria. Much of this difference in the number of YSOs detected is attributable to the more stringent 3.6 µm cut around 14 mag in the SFOG selection criteria; many of the undetected SMOG YSOs are fainter and redder than those selected by the SFOG criteria. This would imply that the total YSO population may be at least as high as 73,247 across the field to the depth of the SMOG observations. Of the remaining 410 unmatched YSOs that were identified with the new method and not in the original catalog, about seventy fall on the edges of the field where there was no overlap between the two IRAC FOVs, and thus were not in the original 'cleaned' source catalog after contaminant removal. A further ∼30 SFOG objects fall in the galactic contaminants color-space of the IRAC four-band CCD. The rest show weak [3.6-4.5] color excess that was not deemed sufficient for selection in the original SMOG criteria. The spatial distribution of the unmatched sources shows that the majority (>75%) are associated with matched YSOs and/or identified clusters. The catalog of Tóth et al. (2014) AKARI YSOs was compared to the SFOG YSOs. The AKARI catalog was based on the Far-Infrared Surveyor All-Sky Survey catalog composed of photometry at four IR wavelength bands centered at 65, 90, 140 and 160 µm. Of the 44,001 AKARI YSOs, 14,986 were located in the SFOG field. Of these, 49 matched to the SFOG catalog within 2 and 261 to within 5 . However, given the AKARI spatial resolution of 1 to 1. 5 at these wavelengths, the AKARI sources contain many individual YSOs and likely significant emission from the clouds surrounding these clusters, as well as deeply embedded YSOs, and are therefore not likely true or unique matches to the IRAC-identified YSOs in the SFOG survey. We Marton et al. (2016) sample. Including the full SFOG catalog, we match 6,425 of 47,338, or 14%, of our YSOs to the Marton sample. A comparison of the CCDs of the two samples of YSOs, shows that ours has a more conservative color cut, with the Marton sample selecting YSOs with WISE1−WISE2 < 0, which accounts for most of the difference in the catalogs. We compared our detection of YSOs in the W5 starforming region to that reported in Koenig et al. (2008) , who published a Spitzer IRAC and MIPS survey of the region conducted during the cryogenic mission phase. They identify 3 deeply embedded objects, 171 Class I, 1,809 Class II, 79 transition disks, (for a total of 2,062 IR-excess YSOs) and 15,709 Class IIa/III stars. Of these, we match 1,091 objects within a 1 radius; none of the deeply embedded objects, 88 Class I (52%), 947 Class II (52%), 18 transition disks (23%), and 38 Class IIa/III (0.2%). The Koenig et al. (2008) Class IIa/III objects have weak excess and are positionally associated with the region, hence the very low detection rate in the SFOG sample. The 50% detection rate for the IR-excess YSOs roughly corresponds to the brighter YSOs in the two shorter IRAC bands. The spatial distribution of the YSOs and the identified clusters match closely to those found in this paper. The SFOG catalog was also compared to the Rivera-Ingraham et al. (2011) (1.1%) of these classes, respectively. The SFOG catalog does not identify the deeply embedded YSOs; some of these were selected using Spitzer MIPS photometry and so may not show excess emission in the shorter wavelengths, the remainder were below 14th magnitude in the IRAC 3.6 µm channel and so were not selected in SFOG. However, the spatial distribution of the W3 YSOs in the SFOG catalog traces a similar cluster distribution to that of Rivera-Ingraham et al. (2011) . We made a further comparison of the SFOG catalog to the Rebull et al. (2011) Spitzer IRAC and MIPS survey of the North American & Pelican nebula star-forming regions. They report a total of 2,196 YSOs in the field, with 262 that lie in the overlap region with the SFOG catalog field. Of these, 132 are matched to within a 1 radius of an SFOG YSO, being again the brighter 3.6 and 4.5 µm sources. We thus draw two conclusions from the comparisons to these well-studied regions and other IR-based YSO catalogs. First, we are finding roughly half of the previously identified YSOs; those that are less embedded and brighter at 3.6 µm. We also used more conservative color cuts and brighter magnitude limits to minimize the number of spurious YSO identifications. Secondly, from the YSOs we do detect, we find that we are reliably identifying the overall spatial structure and main clusters in these regions that were found in the other surveys that used deeper integrations or had coverage in the IRAC 5.8 and 8 µm and MIPS 24 µm bands. We have undertaken a study of the 600 deg 2 SFOG field comprising the GLIMPSE360 and SMOG survey regions. We combined the Spitzer data with 2MASS near-IR photometry, and used the WISE catalog of the field to identify more embedded YSOs. • We identify 42,757 YSOs with IR-excess emission in the GLIMPSE360 and WISE data. When combined with the SMOG field and after removing the spurious WISE sources, we find a total of 47,338 YSOs. • The evolutionary class of the YSOs was determined from the SED slope: 10,461 Class I, 29,552 Class II, and 7,325 Class IIa/III. • We identify 618 reliable clusters in the SFOG field. The ratio of YSOs identified as members of clusters was 25,528/47,338, or 54%. The smallest cluster has 5 members, and the largest has 1,177 members, with a median size of 17 YSOs. Of the 618 clusters, 47 have more than 100 members, and 22 have more than 200 members. • One hundred clusters had a distance estimate from H II regions within their convex hulls. The SEDs of the YSOs in these clusters were fitted using the SEDFitter routine, of these 96 had reliable fits. • From the modeled masses, the IMF was constructed for the clusters across the SFOG field. The slope of the combined IMF was found to be Γ = 2.38 ± 0.20 above 3 M . Dividing the clusters by Galactocentric distances, the slopes were Γ = 1.87 ± 0.31 above 3 M for R Gal < 11.5 kpc and Γ = 1.15 ± 0.24 above 3 M for R Gal > 11.5 kpc. These values are consistent with each other within the uncertainties, and with those obtained in the inner Galaxy high-mass SFRs. The slopes are likely also consistent with a universal Salpeter IMF. Background galaxy contaminants, including candidate AGNs and PAH galaxies, and saturated sources were removed using a cut in color-magnitude space. [3.6] < 6.0 and The full GLIMPSE360 Archive contained 49,378,049 sources. The catalog contained 7,527,352 objects after contaminants were removed. The SFOG fields lie in the direction of the outer Galaxy where the level of shielding from the Galactic center from both the stellar population and dust component is reduced and thus the extragalactic background is expected to be higher. Three combinations of 2MASS and IRAC bands were used to select for less extincted objects. The photometry was first dereddened, and only sources with good values of extinction were included for selection. For these three selection criteria: 20,339, 24,560, and 16,941 candidate YSOs were selected respectively, for a combined total of 28,837 YSOs identified using the GLIMPSE360 photometry. Following the process laid out by Fischer et al. (2016) , spurious detections were cleaned from the catalog. The first step was to remove those sources with uppercase flags in bands W1, W2, and W3. Upper limits in bands W1, W2, W3 are then removed, and a saturation cut-off of W 1 > 5 applied to the data. The initial catalog contained 14,483,596 sources. The remaining catalog contained 543,457 sources, or ≈ 4% of the original catalog. The contaminating background galaxies and source selection criteria for the WISE data were taken from Koenig et al. (2014); Fischer et al. (2016) and adapted to the requirements of the GLIMPSE360 field. We adjusted the criteria for removal of AGN and star-forming galaxy contaminants slightly from those of Koenig et al. (2014) Of those 14,483,596 sources, 13,940,139 were identified as contaminants, leaving a cleaned catalog of 543,457 sources. The clean catalog was then searched for YSOs according to criteria taken from the Fischer et al. (2016) and Koenig et al. (2014) papers. From these there were 2,245 transition disk candidates and 7,094 YSOs identified from the four-band WISE diagram, and 7,244 Class I and 10,447 Class II sources identified from the WISE 3-band diagram. There was a total of 20,892 candidate YSOs identified with WISE. American Astronomical Society Meeting Abstracts #213 213 Astronomical Data Analysis Software and Systems XII The Build-up of Galaxies Through Multiple Tracers and Facilities IAU Symp. 115, Star Forming Regions VizieR Online Data Catalog & GLIMPSE360 Team This work is based on observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology under NASA contract 1407. We gratefully acknowledge funding support for this work from NASA ADAP grant NNX16AF37G. This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. Support for the IRAC instrument was provided by NASA through contract 960541 issued by JPL. This publication makes use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of The removal of contaminating sources and the selection of YSOs in the GLIMPSE360 field were undertaken in a different manner to those implemented for the SMOG field, with IRAC four band coverage which was based on the methods of Gutermuth et al. (2008b) and Gutermuth et al. (2009) . The location of the SMOG field contaminants in color-space were used as the basis for identifying the locations of the contaminants in the two band IRAC data of GLIMPSE360. The IRAC+WISE YSOs were selected following the same selection cut-offs as were applied to the WISE sources. The cleaned WISE catalog was matched to the contaminant-removed IRAC catalog, and sources with a 1. 5 or closer spatial coincidence were considered to be the same object, giving 318,588 objects in the joined catalog. The YSOs were then selected by replacing the WISE bands 1 & 2 with IRAC bands 1 & 2, and replacing this photometry in the WISE color-color diagram YSO selections. The 2MASS selection criteria were not applied here as they replicate the IRAC+2MASS selection from the GLIMPSE360 data. A total of 11,196 candidate YSOs were selected using this method.