Statistics of gammaray point sources below the Fermi detection limit
Abstract
An analytic relation between the statistics of photons in pixels and the number counts of multiphoton point sources is used to constrain the distribution of gammaray point sources below the Fermi detection limit at energies above 1 GeV and at latitudes below and above . The derived sourcecount distribution is consistent with the distribution found by the Fermi collaboration based on the first Fermi point source catalogue. In particular, we find that the contribution of resolved and unresolved active galactic nuclei (AGN) to the total gammaray flux is below 20 to 25%. In the best fit model, the AGNlike point source fraction is %. Using the fact that the Galactic emission varies across the sky while the extragalactic diffuse emission is isotropic, we put a lower limit of 51% on Galactic diffuse emission and an upper limit of 32% on the contribution from extragalactic weak sources, such as starforming galaxies. Possible systematic uncertainties are discussed.
Keywords: galaxies: active; gamma rays: diffuse background; gamma rays: general; methods: statistical; quasars: general
1 Introduction
Faint gammaray point sources cannot be detected individually, however their presence affects the statistics of photons across the sky. One can use the observed statistics of photons to infer some general properties about the population of point sources below the Fermi detection limit. The idea of using the statistics of photon counts or intensity maps in the study of faint point sources is familiar in radio and Xray observations (e.g., Scheuer, 1957; Condon, 1974; Scheuer, 1974; Hasinger et al., 1993; Miyaji & Griffiths, 2002; Sołtan, 2011). A closely related subject is the use of extreme statistics to constrain the nonGaussianity of cosmic microwave background data, e.g., Colombi et al. (2011) and references therein. In this paper, we extend the prior analysis (Slatyer & Finkbeiner, 2010) and use the points in cells statistics to understand the population of gammaray point sources.
In a larger context, the problem is to separate different sources of gammaray emission. The standard strategy is either to use templates (Dobler et al., 2010) that trace the sources or else to simulate the cosmic ray propagation and gammaray production using, e.g., Galprop (Strong et al., 2007, 2009). In this paper we would like to adopt a different methodology: we assume some general properties of the sources in order to obtain modelindependent constraints on the contribution from these sources.
We separate the sources based on their statistical properties. In the paper we consider three sources of gammarays at high latitudes. The first source is a diffuse source that varies over large angles. The main contribution to this source comes from the Galactic diffuse emission ( production, ICS photons, and bremsstrahlung) and, possibly, a population of faint Galactic point sources, such as millisecond pulsars (FaucherGiguère & Loeb, 2010; Malyshev et al., 2010; ?). We will call this source nonisotropic Galactic diffuse emission. It puts a lower limit on the actual Galactic emission.
The second source corresponds to an isotropic distribution of gammarays, i.e., the statistics of photons across the sky is consistent with the Poisson distribution for this source. The isotropic flux has contributions from the homogeneous part of the Galactic diffuse emission, from diffuse intergalactic emission (references can be found in Stecker & Venters, 2010), and from a population of very weak extragalactic sources that on average emit much less than one photon during the time of observation, such as starforming galaxies (e.g., Pavlidou & Fields, 2002; Fields et al., 2010), starburst galaxies (Thompson et al., 2007), galaxy clusters (Berrington & Dermer, 2003), etc. The value of the isotropic emission puts an upper limit on the gammaray flux from very weak extragalactic sources, as well as an upper limit on the contamination from cosmic rays (Abdo et al., 2010b).
The third source is a population of point sources modeled by a broken powerlaw sourcecount distribution. We assume that the point sources are distributed homogeneously over the sky. The statistics of photons coming from these point sources has a nontrivial form (derived in Appendix A) different from the Poisson statistics. These sources model a population of AGNlike point sources (Stecker et al., 1993; Padovani et al., 1993; Chiang et al., 1995; Stecker & Salamon, 1996; Mücke & Pohl, 2000; Abazajian et al., 2010; Neronov & Semikoz, 2011).
We note that in our approach it is impossible to separate the isotropic part of Galactic emission from the diffuse extragalactic background (EGB). Consequently, all results will be quoted either with respect to the total gammaray flux or in absolute values. We find that at high latitudes (below and above ) for energies above 1 GeV the contribution of AGNlike point sources (both resolved and unresolved) to the total gammaray flux is %, the contribution from Galactic diffuse emission is above 51% and the contribution from extragalactic weak sources is below 32%. Using these values, we estimate that the contribution of unresolved AGNlike point sources to the EGB model in Abdo et al. (2010b) is below , which is consistent with Abdo et al. (2010c).
The paper is organized as follows. In Section 2, we describe a general model and a fitting algorithm. In Section 3, we present an analysis of the Fermi data. Section 4 has discussion. In Appendix A, we derive the statistics of photons coming from a population of point sources. In Appendix B, we determine a model for nonisotropic Galactic diffuse emission. In Appendix C, we repeat the data analysis of Section 3 for different pixel sizes.
2 Model
In this section we describe the model that we later use to fit the Fermi gammaray data. At first, we present an analytic relation between the counts of photons in pixels and the statistics of point sources. Then we take into account the detector point spread function (PSF) and describe a model of Galactic diffuse gammaray emission. The model that we present in this section is rather general. An application of this analysis for the Fermi data is presented in Section 3.
2.1 Statistics of Photon Counts and Point Sources
We will assume a pixelation of the sphere with pixels of equal area. Usually a model of gammaray emission is represented in terms of the expected number of gammarays in pixel . The problem is that the number and the positions of weak point sources are not known. The position of a point source can be found only by a “detection” of the source. However, if one is interested only in the total number of sources of a given flux, then a “detection” of sources is not necessary. In this paper, we use the statistics of photon counts in pixels, in order to find the sourcecount distribution of point sources without actual identification of the sources.
In order to find the statistics of photon counts in pixels we calculate the total number of pixels that contain photons. In this calculation, the information about the position of the pixels on the sky is not preserved, but the remaining information may be sufficient to constrain general properties of the population of sources. There are two competing conditions for this method to work. On the one hand, there should be sufficiently many pixels to determine the statistics of photon counts, while, on the other hand, the size of the pixels cannot be too small compared to the PSF, otherwise the statistics corresponding to the presence of point sources cannot be distinguished from nonisotropic diffuse emission.
If the total number of pixels is and is the observed number of pixels with photons, then we can estimate the probability to find photons in a pixel as
(1) 
For large , the statistical uncertainty of is approximately (there is a small correction due to the fact that the total number of pixels is fixed, i.e., the process is multinomial). If there are no point sources and all the emission comes from an isotropic diffuse source, then the probability distribution is the Poisson probability of getting photons given the mean rate. In the presence of point sources the probability distribution is more complicated than the Poisson distribution. In the context of Xray point sources, the probability distribution is called a distribution (or a diagram) and is usually computed with the help of Monte Carlo simulations (e.g., Miyaji & Griffiths, 2002).
In this paper we find an analytic relation between the probability distribution of photon counts in pixels and the sourcecount distribution. We will use this relation to find a sourcecount distribution that has the best fit to the observed probability distribution determined in Equation (1).
In calculations we use the method of probability generating functions (an introduction can be found in, e.g., Hoel et al., 1971, Section 3.6). For a given discrete probability distribution , the corresponding generating function is defined as a power series in an auxiliary variable
(2) 
If the generating function is known, then the probabilities can be found by picking the coefficient in front of , or, equivalently, by differentiating with respect to
(3) 
An important property of probability generating functions is that a sum of two independent random variables corresponds to a product of the corresponding probability generating functions (Hoel et al., 1971). For example, if there are two independent sources of gammarays on the sky, e.g., Galactic and extragalactic, then the probability generating function for the photon counts with the two sources overlaid is given by the product of probability generating functions corresponding to the two sources (see Appendix A for more details).
Let denote the average number of sources inside a pixel that emit exactly photons during the time of observation. Using the sumproduct property of probability generating functions (Hoel et al., 1971, Section 3.6), we derive in Appendix A the following relation between the probability generating function for the photon counts and the expected number of photon sources
(4) 
This formula provides an analytic relation between the expected numbers of photon sources and the statistics of photons in pixels. The probabilities to observe photons in a pixel are determined by expanding the right hand side of this equation and picking the coefficient in front of .
If we substitute , then Equation (2) becomes a discrete Fourier transform analog of the probability characteristic functions (e.g., Hoel et al., 1971). The probability characteristic functions were extensively used in the study of radio point sources below detection limit (e.g., Scheuer, 1957). Similarly to probability generating functions, the characteristic function for a sum of two independent random variables corresponds to a product of the two characteristic functions.
The statistics of point sources can be described in terms of the differential source counts as a function of flux , denoted by . One of the simplest forms for the sourcecount distribution is a broken powerlaw
(5) 
where we have to require and
in order to have a finite total flux
Let be an average flux from a source (normalized to give the number of photons during the time of observation). Then the probability to observe photons from this source is
(6) 
If the sources are distributed isotropically, then the average number of photon sources inside a pixel is given by the Poisson probabilities for sources with flux to emit photons
(7) 
where is the angular area of the pixel.
2.2 Psf
In the case of nonzero PSF a source in some pixel may contribute gammarays to nearby pixels as well. In order to correctly estimate the number of point sources of a certain strength, we need to know how often a gammaray from a source in one pixel is detected in a different pixel.
We can represent the effect of the PSF as a smearing of a point source over some area, so that the flux from the point source is split between several pixels. We determine the average properties of the flux splitting from Monte Carlo simulations. In the simulation, we assume that the profile of the PSF is Gaussian. In order to calculate the average splitting, we put Gaussians at random positions on the sphere and integrate every Gaussian over the pixels to get a collection of fractions , such that (in practice, one can truncate at some minimal value of ). Denote by the number of fractions for Gaussians that fall within , then the average distribution of fractions is defined as
(8) 
This distribution is normalized as
(9) 
The case of zero PSF corresponds to . Sums over fractions are substituted by the integral . In particular, the contribution of point sources with intrinsic flux between and to the total number of photon sources is
(10) 
The expected number of photon sources inside a pixel including the effect of PSF is
(11) 
An example of the function relevant for the data analysis in Section 3 is presented in Figure 1. If we rescale the integration variable in Equation (11) by , then the effect of the PSF can be represented as a change of the sourcecount function
(12) 
The apparent number of sources with small fluxes in is larger than the corresponding number of physical sources in due to contribution from point sources in nearby pixels, whereas the apparent number of sources with large fluxes is smaller due to loss of the flux to nearby pixels. One can see the effect of PSF in Figure 2 on the right. On the same figure, we also plot the expected number of photon sources for discrete , derived from .
2.3 The Model
In Equation (4), the probabilities and the average expected number of photon sources may depend on the position of the pixel. The quantity that we use is the averaged probability , where is the number of pixels with photons. The generating function for averaged probabilities is given by
(13) 
This form of the generating function can be used when the distribution of sources is not isotropic or the exposure is nonuniform. If the exposure is sufficiently uniform, then we may assume that the distribution of extragalactic sources is independent of the pixel position and can be taken out of the sum over pixels. The diffuse emission corresponds to “onephoton” source density that depends on the pixel, . Thus, the expression on the righthand side of Equation (13) can be split into a product
(14) 
The second term in the product on the right hand side corresponds to the isotropic distribution of AGNlike point sources
(15) 
while the first term corresponds to diffuse emission or, indistinguishably, onephoton sources
(16) 
It is convenient to separate the diffuse emission into a nonisotropic part that puts a lower limit on Galactic diffuse emission and an isotropic part that consists of isotropic Galactic component and a possible population of weak extragalactic sources (additional to AGNlike sources)
(17) 
The main problem with the nonisotropic diffuse emission is a variability on large scales (on small scales the variability in Galactic emission is subdominant to Poisson noise and variability due to the presence of undetected point sources). In Appendix B we construct a model for the nonisotropic component that varies on large scales. We mask the known point sources and use low multipole spherical harmonics in order to filter out the smallscale variations while preserving the largescale ones. We add a constant to so that . In the following, the nonisotropic component of diffuse emission is fixed. The corresponding generating function is
(18) 
We denote the probability generating function for the isotropic emission as
(19) 
The total generating function for the probability distribution of photons in pixel is the product of the three components
(20) 
The parameters of the pointsource distribution in Equation (5) and the level of the isotropic flux are found from fitting the probability distribution given in Equation (20) to the observed probability distribution. The details of the fitting algorithm are presented in the next subsection.
This simple factorized form of the generating function is valid in every pixel. In general, this factorization is not possible if we take an average over pixels, unless at most one sources is nonisotropic and the other sources can be taken out of the average over pixels. The factorization leads to a significant simplification of calculations at the expense that, due to a variation in exposure, a part of isotropic EGB may be misinterpreted as a nonisotropic Galactic foreground. This effect is expected to be small, provided that the Fermi LAT exposure is sufficiently uniform (Atwood et al., 2009). Another possible source of systematic uncertainty is that the distribution of galaxies in the local neighborhood of the Milky Way is nonuniform and a part of emission from these galaxies can be misinterpreted as Galactic emission.
2.4 Fitting Algorithm
The algorithm has two main parts:

Determine the nonisotropic part of Galactic diffuse emission. We mask the bright point sources and use low multipoles of spherical harmonics to find the variable part of the remaining flux (Appendix B). Then we add a constant such that . In the following is fixed, i.e. we do not vary this component in fitting the pixel counts.

Use pixel counts to find and . The sourcecount distribution has four parameters: normalization, break, and two indices. Together with the level of isotropic contribution, , this gives five fitting parameters which we denote by , where denotes the normalization in the sourcecount distribution in Equation (5).
The model prediction for the pixel counts is determined by multiplying the righthand side of Equation (20) by
(21) 
The expected pixel counts are found by expanding the lefthand side of this equation in powers of and picking the coefficient in front of . Given the observed number of pixels with photons, the likelihood of is estimated by the Poisson probability
(22) 
The overall likelihood of the model is estimated as
(23) 
Here the product is over the number of photons , while the data are represented by the number of pixels that contain photons. This is different from a usual coordinate space fit of maps, where the product is over pixels and the data are the number of photons in a pixel .
We use this likelihood function to find the bestfit parameters for a given set of observed pixel counts . The significance of deviation of from can be estimated as
(24) 
In the case of large counts the likelihood is well approximated by the Gaussian distribution and is the deviation in sigma values.
3 Data analysis
We consider 11 months of Fermi data (August 4, 2008  July 4, 2009) that were used for the first Fermi pointsource catalog (Abdo et al., 2010a) and for the Fermi analysis of the sourcecount distribution (Abdo et al., 2010c).
In order to reduce the PSF, we use frontconverted gammarays, i.e.,
the gammarays converted in the front part of the FermiLAT tracker
(Atwood et al., 2009),
with energies between 1 GeV and 300 GeV.
The corresponding (quadratic) average of the PSF is .
The distribution of pixel counts is presented in Figure 2 on the left. We model this distribution by a combination of three components: nonisotropic Galactic diffuse emission derived in Appendix B, isotropic distribution of photons, and the distribution of AGNlike point sources described in Section 2.3. The bestfit models for the isotropic flux and the AGNlike point sources are presented in Figure 2 on the right.
In fitting we consider pixel counts with , i.e., the number of data points is . We use for the Galactic diffuse emission model. The sampling is performed with MetropolisHastings Markov Chain Monte Carlo method with 1500 steps. The average model is presented in Figure 2 and compared with Fermi models in Table 1. We find that the position of the break is almost a flat direction. The results of fitting with a set of fixed break positions are presented in Figure 3.
Analysis  

This work  
Fermi MeV  
Fermi 1  10 GeV 
In units of .
Rescaled from the break in Abdo et al. (2010c)
according to energy spectrum .
We find that the contribution of point sources to the gammaray flux is . This fraction is obtained by integrating the bestfit sourcecount model presented in Figure 2 and in Table 1 from zero to infinity. This fraction is not model independent, i.e., it may change for different shapes of the sourcecount function used in fitting, but, provided that the broken powerlaw sourcecounts distribution gives a good fit to the data, we do not expect a large systematic uncertainty due to the change in the shape of the fitting function. In part this can be justified by looking at the models with fixed positions of the break: most of the models in Figure 3 have the pointsource fraction .
At smaller position of the break the contribution of point sources becomes larger ( for the left point in Figure 3) but the fit has smaller likelihood. It would be interesting to estimate constraints on models involving high redshift AGNs to explain the unresolved part of the EGB (e.g., Abazajian et al., 2010; Neronov & Semikoz, 2011) using the statistics of photon counts.
We can also put an approximate lower bound of 51% on Galactic diffuse emission, and an approximate upper bound of 32% on isotropic emission consistent with Poisson statistics. In terms of the absolute flux values, the total gammaray flux above 1 GeV for is , the diffuse Galactic flux is , the isotropic component of the flux is , the flux from AGNlike point sources with parameters given in Table 1 and in Figure 2 is .
Using this value of the total flux from AGNlike point sources, we can estimate the gammaray flux from undetected sources. The total flux in the energy range 1  100 GeV from point sources at detected by Fermi (Abdo et al., 2010a) is . Consequently, the flux from undetected AGNlike point sources above can be estimated as . This value is reasonable, since from Figure 4 it follows that the detected point sources give a good approximation to the bestfit distribution of sources down to fluxes below the break in , while the contribution of point sources to the gammaray background is saturated around the break. The extragalactic diffuse emission can be estimated from the model presented in Abdo et al. (2010b). In this model, the EGB flux above 1 GeV is . The contribution of unresolved AGN’s to this EGB flux is about . This value is consistent with the estimations in Abdo et al. (2010c). We would like to stress that our analysis is different from Abdo et al. (2010c). We do not count any sources, but rather do a “blind” statistical analysis. The consistency of results may serve as a nontrivial check of both methods.
4 Discussion
One of the most interesting problems in gammaray astrophysics is to understand the origin of EGB flux, which is defined as the total gammaray flux minus the contribution from resolved point sources, minus the Galactic foreground (Abdo et al., 2010b). In general, we can separate the sources of diffuse EGB into two big classes: nonPoisson sources and Poissonlike sources. In the first case, the contribution comes from sources below detection limit that may emit several photons during the observation time, these are the AGNlike point sources (e.g., Padovani et al., 1993; Chiang et al., 1995; Stecker & Salamon, 1996; Abazajian et al., 2010; Neronov & Semikoz, 2011). In this case the statistics of photon counts in pixels across the sky will be different from Poisson statistics due to correlation among photons coming from the same point source. In the second case, there is a large number of sources emitting on average much fewer than one photon during the observation time, e.g., starforming and starburst galaxies (Pavlidou & Fields, 2002; Thompson et al., 2007; Fields et al., 2010). The statistics of photons in this case will be very close to Poisson statistics.
These two possibilities are relatively hard to separate based only on energy spectrum arguments (see, however, Abazajian et al., 2010; Stecker & Venters, 2010; Neronov & Semikoz, 2011). Alternatively, one can use a correlation between gammaray emission and multiwavelength emission from AGNs (Urry & Padovani, 1995; Mücke & Pohl, 2000; Abdo et al., 2010d; Angelakis et al., 2010; Li & Cao, 2011) together with AGN population studies (e.g., Padovani & Giommi, 1995; Hopkins et al., 2007; Padovani et al., 2007; Rigby et al., 2008; Hodge et al., 2009; Smolčić et al., 2009) to infer the AGN contribution to gammaray background.
In this paper we show that the statistics of photons in pixels can also be used to constrain the distribution of point sources below the detection limit. We find that there is a slight preference for two populations of point sources: the AGNlike point sources with a break at relatively high flux and a population of faint sources such as starforming galaxies rather than a single population of AGNlike point sources with a break at a smaller flux. This statement is insensitive to the energy spectrum since we use the integrated flux above 1 GeV but it may depend on the form of the function used to fit the point source number counts. Possible sources of systematic uncertainty include the source smearing due to PSF, nonisotropic Galactic diffuse emission, nonisotropic emission (on large scales) from nearby galaxies, and nonhomogeneous exposure. The interpretation of the pointsource number counts may also be affected by clustering of sources on scales smaller than the PSF.
A stronger constrain on the population of gammaray point sources may come from an unpixelized analysis of the data. In this more general analysis, the likelihood of a model depends on the full gammaray data (positions on the sky, energy, arrival time) rather than the counts of gammarays in pixels. This likelihood would have the form of an integral over flux times the source number counts, times an integral over all possible positions of the sources on the sky. The photon counts approximation is simpler computationally and already gives rough constraints on the source population, while the full unbinned analysis may provide stronger constraints at the expense of computational complexity. Examples of unbinned analysis using some part of gammaray data include the gammaray twopoint correlation function (e.g., Ave et al., 2009; GeringerSameth & Koushiappas, 2010), angular power spectrum (Ando & Komatsu, 2006; Hensley et al., 2010), nearest neighbor statistics (Sołtan, 2011), etc.
We also believe that the techniques of generating functions in the study of photon counts statistics, known for some time in radio and Xray observations, may have further important applications in data analysis of current and future radio, IR, Xray, and gammaray observations.
Acknowledgments.
The authors are thankful to Kevork Abazajian, Marko Ajello, Jo Bovy, Christopher Kochanek, Jennifer SiegalGaskins, Tracy Slatyer, Tonia Venters, and Itay Yavin for stimulating discussions and useful comments. We are also thankful to Ilias Cholis for the help in selecting the Fermi gammaray data. This work is supported in part by the Russian Foundation of Basic Research under Grant No. RFBR 100201315 (D.M.), by the NSF Grant No. PHY0758032 (D.M.), by Packard Fellowship No. 19991462 (D.M.), by NASA grant NNX08AJ48G (D.W.H.), and by the NSF grant AST0908357 (D.W.H.). Some results in this paper have been obtained using the HEALPix package (Górski et al., 2005).
Appendix A Derivation of statistics
In this Appendix we derive the photon counts statistics in the presence of point sources. The derivation is based on the following well known property of generating functions for probabilities: the generating function of a sum of two independent sources is the product of the corresponding generating functions (e.g., Hoel et al., 1971, Section 3.6). Indeed, consider a box where we can put red balls with probabilities , , and blue balls with corresponding probabilities , , then the probability to find balls of any color is
(A1) 
The last equation is the rule for multiplication of polynomials or power series. Let
(A2) 
then
(A3) 
As before, denote the probability to observe photons in a pixel by . The corresponding generating function is
(A4) 
The knowledge of the generating function is equivalent to the knowledge of all , provided that every can be found from by picking the term in front of .
We will now derive the generating function for probabilities to observe photons from a collection of point sources. Denote by the average number of point sources inside a pixel that emit exactly photons during the time of observation. We assume that the probability to find photon sources in a pixel is given by the Poisson statistics
(A5) 
then the probability to find photons from photon sources is
(A6) 
The corresponding probability generating function is
(A7) 
The generating function of probabilities to observe photons from any sources is the product of the generating functions for every
(A8)  
If we denote and
(A9)  
(A10) 
then Equation (A8) can be represented as
(A11) 
which is the discrete Fourier transform analog of characteristic functions considered in, e.g., Scheuer (1957).
Appendix B Nonisotropic Galactic diffuse emission
In this appendix, we describe a model for the nonisotropic Galactic diffuse emission. This signal is “filtered” by low multipoles of the gammaray data. In this approach, the homogeneous part of Galactic emission is indistinguishable from extragalactic flux. Also some part of extragalactic emission from galaxies and galaxy clusters close to Milky Way may lead to a signal that varies over the sky and can be misinterpreted as nonisotropic Galactic emission. However, the majority of extragalactic emission comes from higher redshifts where the distribution of extragalactic sources is sufficiently uniform on large scales.
In Figure 5, we plot the angular power spectrum for the data at high latitudes () before and after masking the point sources in the first Fermi catalog (Abdo et al., 2010a). The angular power spectrum is
(B1) 
where ’s are the spherical harmonics coefficients of , where is the number of photon counts inside pixel and is the mask function. The mask function is equal to one (zero) for (). In the case of masked point sources, we also have when a Fermi point source is inside the pixel or within two PSF from the boundary of the pixel. We also subtract the average of the data within the unmasked region in order to avoid a nontrivial contribution from a constant source inside the window. We choose the normalization such that for the Poisson noise .
The algorithm for estimating the component that varies on large scales is as follows:

Calculate the ’s inside a window larger than needed (in order to avoid edge effects). We use for the data above . We also fill in the pixels with point sources by the average of nearest neighbors.

The largescale distribution of gammarays is defined as
(B2) where we find the coefficients and from the best fit of the diffuse model to the photon counts in the pixels without point sources (see Figure 7).

We represent so that
(B3) is nonnegative for . In fitting to full data, and are fixed, while is allowed to vary together with parameters describing the AGNlike point sources.
An example of the diffuse emission model for is presented in the middle plot of Figure 6. The top plot represents the counts of photons in pixels inside the window with masked Fermi gammaray point sources. The bottom plot represents the deviation of the model from the data in “sigma” values. In Figure 4 we study the effect of changing . The difference is rather small, i.e., already captures the largescale distribution of gammarays reasonably well.
Appendix C Variation of pixel size
In the data analysis in Section 3, we use the pixel size of about corresponding to HEALPix parameter nside = 32. In this appendix, we repeat the analysis of Section 3 for different sizes of pixels, in particular we consider nside = 64 and nside = 16 corresponding to pixel sizes and respectively.
The results for nside = 64 are presented in Figure 8. The total number of pixels inside the window below and above in latitude is about 25,000, i.e., the statistics of pixel counts is relatively good, but the effect of PSF becomes more significant than in the case of nside = 32. For PSF = , less than about 60% of the flux from a source can be inside one pixel. In particular, one can note that the number of sources with small fluxes is about two times larger than in the input model. Nonetheless, in the case of nside = 64 the bestfit model is similar to nside = 32 case. The position of the break, , is somewhat lower than the values in Table 1 but the overall source count distribution (Figure 9) is consistent with the source counts found by Fermi collaboration (Abdo et al., 2010a, c).
The results of fitting for nside = 16 are presented in Figure 10. In this case, the PSF can be neglected. However, the number of pixels inside the window is only about 1,500. As one can see on the left side of Figure 10, the error bars are rather large and the photon counts in most of the pixels are dominated by the contribution from diffuse emission. As a result, our method is not sensitive to the contribution from point sources for nside = 16. In particular, the MCMC is dominated by models with large . The best fit value is an order of magnitude larger than the values in Table 1. The reason is that models with large have a comparable likelihood to the models with small but there are many more ways to choose a large .
We conclude that our method is relatively robust to the change of pixel size: the cases of nside = 32 and nside = 64 give very similar results, unless the number of pixels is not sufficient to draw any conclusion about the statistics of sources, as in the case of nside = 16.
Footnotes
 affiliation: dm137 at nyu.edu
On leave of absence from ITEP, Moscow, Russia, B. Cheremushkinskaya 25  A benchmark model is provided by a static universe filled homogeneously with sources that have intrinsic luminosity . The number of sources with flux larger than is . The differential number count is . Index provides a benchmark value for differential source count at large . A break or a cutoff at small flux (large distance to the source) is necessary given the finite size of the visible universe.
 http://wwwglast.slac.stanford.edu/software/IS/glast_lat_performance.htm
 We consider the cases of nside = 16 and nside = 64 in Appendix C.
References
 Abazajian, K. N., Blanchet, S., & Harding, J. P. 2010, arXiv:1012.1247
 Abdo, A. A., et al. 2010a, ApJS, 188, 405
 —. 2010b, Physical Review Letters, 104, 101101
 —. 2010c, ApJ, 720, 435
 —. 2010d, ApJ, 716, 30
 Ando, S., & Komatsu, E. 2006, Phys. Rev. D, 73, 023521
 Angelakis, E., Fuhrmann, L., Nestoras, I., Zensus, J. A., Marchili, N., Pavlidou, V., & Krichbaum, T. P. 2010, arXiv:1006.5610
 Atwood, W. B., et al. 2009, ApJ, 697, 1071
 Ave, M., et al. 2009, J. Cosmology Astropart. Phys, 7, 23
 Berrington, R. C., & Dermer, C. D. 2003, ApJ, 594, 709
 Chiang, J., Fichtel, C. E., von Montigny, C., Nolan, P. L., & Petrosian, V. 1995, ApJ, 452, 156
 Colombi, S., Davis, O., Devriendt, J., Prunet, S., & Silk, J. 2011, arXiv:1102.5707
 Condon, J. J. 1974, ApJ, 188, 279
 Dobler, G., Finkbeiner, D. P., Cholis, I., Slatyer, T., & Weiner, N. 2010, ApJ, 717, 825
 FaucherGiguère, C., & Loeb, A. 2010, J. Cosmology Astropart. Phys, 1, 5
 Fields, B. D., Pavlidou, V., & Prodanović, T. 2010, ApJ, 722, L199
 GeringerSameth, A., & Koushiappas, S. M. 2010, arXiv:1012.1873
 Górski, K. M., Hivon, E., Banday, A. J., Wandelt, B. D., Hansen, F. K., Reinecke, M., & Bartelmann, M. 2005, ApJ, 622, 759
 Hasinger, G., Burg, R., Giacconi, R., Hartner, G., Schmidt, M., Trumper, J., & Zamorani, G. 1993, A&A, 275, 1
 Hensley, B. S., SiegalGaskins, J. M., & Pavlidou, V. 2010, ApJ, 723, 277
 Hodge, J. A., Zeimann, G. R., Becker, R. H., & White, R. L. 2009, AJ, 138, 900
 Hoel, P. G., Port, S. C., & Stone, C. J. 1971, Introduction to Probability Theory (Boston, MA: Houghton Mifflin)
 Hopkins, P. F., Richards, G. T., & Hernquist, L. 2007, ApJ, 654, 731
 Li, F., & Cao, X. 2011, arXiv:1103.4545
 Malyshev, D., Cholis, I., & Gelfand, J. D. 2010, ApJ, 722, 1939
 Miyaji, T., & Griffiths, R. E. 2002, ApJ, 564, L5
 Mücke, A., & Pohl, M. 2000, MNRAS, 312, 177
 Neronov, A., & Semikoz, D. V. 2011, arXiv:1103.3484
 Padovani, P., Ghisellini, G., Fabian, A. C., & Celotti, A. 1993, MNRAS, 260, L21
 Padovani, P., & Giommi, P. 1995, ApJ, 444, 567
 Padovani, P., Giommi, P., Landt, H., & Perlman, E. S. 2007, ApJ, 662, 182
 Pavlidou, V., & Fields, B. D. 2002, ApJ, 575, L5
 Rigby, E. E., Best, P. N., & Snellen, I. A. G. 2008, MNRAS, 385, 310
 Scheuer, P. A. G. 1957, in Proceedings of the Cambridge Philosophical Society, Vol. 53, Proceedings of the Cambridge Philosophical Society, 764–773
 Scheuer, P. A. G. 1974, MNRAS, 166, 329
 Slatyer, T. R., & Finkbeiner, D. P. 2010, MNRAS, 405, 1777
 Smolčić, V., et al. 2009, ApJ, 696, 24
 Sołtan, A. M. 2011, arXiv:1101.0256
 Stecker, F. W., & Salamon, M. H. 1996, ApJ, 464, 600
 Stecker, F. W., Salamon, M. H., & Malkan, M. A. 1993, ApJ, 410, L71
 Stecker, F. W., & Venters, T. M. 2010, arXiv:1012.3678
 Strong, A. W., Moskalenko, I. V., Porter, T. A., Jóhannesson, G., Orlando, E., & Digel, S. W. 2009, arXiv:0907.0559
 Strong, A. W., Moskalenko, I. V., & Ptuskin, V. S. 2007, Annual Review of Nuclear and Particle Science, 57, 285
 Thompson, T. A., Quataert, E., & Waxman, E. 2007, ApJ, 654, 219
 Urry, C. M., & Padovani, P. 1995, PASP, 107, 803