Issue 
A&A
Volume 614, June 2018



Article Number  A36  
Number of page(s)  6  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201730872  
Published online  11 June 2018 
Quantifying systematics from the shear inversion on weaklensing peak counts
^{1}
Institute for Astronomy, University of Edinburgh, Royal Observatory,
Blackford Hill,
Edinburgh
EH9 3HJ,
UK
email: calin@roe.ac.uk
^{2}
Service d’Astrophysique, CEA Saclay, Orme des Merisiers,
Bât. 709,
91191
GifsurYvette,
France
^{3}
Fenglin Veteran Hospital,
2 Zhongzheng Rd. Sec. 1, Fenglin Township,
Hualien
97544,
Taiwan
^{4}
Sorbonne Universités, UPMC Univ. Paris 6 et CNRS, UMR 7095, Institut d’Astrophysique de Paris,
98 bis boulevard Arago,
75014
Paris,
France
Received:
27
March
2017
Accepted:
13
February
2018
Weaklensing peak counts provide a straightforward way to constrain cosmology by linking local maxima of the lensing signal to the mass function. Recent applications to data have already been numerous and fruitful. However, the importance of understanding and dealing with systematics increases as data quality reaches an unprecedented level. One of the sources of systematics is the convergenceshear inversion. This effect, inevitable when carrying out a convergence field from observations, is usually neglected by theoretical peak models. Thus, it could have an impact on cosmological results. In this paper, we study the bias from neglecting (mismodeling) the inversion. Our tests show a small but nonnegligible bias. The cosmological dependence of this bias seems to be related to the parameter Σ_{8} ≡ (Ω_{m}∕(1 − α))^{1 − α}(σ_{8}∕α)^{α}, where α = 2∕3. When this bias propagates to the parameter estimation, we discovered that constraint contours involving the dark energy equation of state can differ by 2σ. Such an effect can be even larger for future highprecision surveys and we argue that the inversion should be properly modeled for theoretical peak models.
Key words: gravitational lensing: weak / largescale structure of Universe / methods: numerical
© ESO 2018
1 Introduction
The observation of background light sources results in distorted images due to light deflection. This effect is called gravitational lensing. The lensing signal tracks the matter distribution of the Universe and provides evidence of cosmic history. In the weaklensing (WL) regime, the information about structure formation is sensitive to both linear and nonlinear scales, providing a crucial tool to constrain cosmology (see for example Kilbinger 2015 for a review). With upcoming large surveys such as Euclid (Laureijs et al. 2011) and the Large Synoptic Survey Telescope (LSST; LSST Science Collaboration 2017), WL has been considered as an important probe to find out, notably, the nature of dark energy and the gravity law at cosmic scales.
Peak counts are powerful nonGaussian statistics that allow us to constrain cosmology from WL. Peaks of high signaltonoise ratio (S/N) are shown to be good tracers of massive halos (Yang et al. 2011; Lin & Kilbinger 2015a, hereafter Paper I). The more massive halos there are, the more likely it is to find high S/N peals in a lensing map. While older studies (Kruse & Schneider 1999, 2000; Reblinsky et al. 1999; Bartelmann et al. 2002; Hamana et al. 2004; Wang et al. 2004) seem to only focus on reducing false positive and false negative detections, recent studies (Jain & Van Waerbeke 2000; Dietrich & Hartlap 2010; Kratochvil et al. 2010; Fan et al. 2010; Paper I among others) tend to model true and false detections together in order to probe cosmology. Indeed, WL peak counts are sensitive to changes of shape of the mass function. In the literature, it has been shown that peak counts alone constrain cosmology more strictly than twopointcorrelation functions and the combination of both statistics improves the constraints (Dietrich & Hartlap 2010; Liu et al. 2015a). Applications to observational data are also multiple. The data of the CanadaFranceHawaiiTelescope Lensing Survey (CFHTLenS) have been analyzed by Liu et al. (2015a) and Liu et al. (2016), the CFHT Stripe82 Survey by Liu et al. (2015b), Dark Energy Survey Science Verification (DESSV) data by Kacprzak et al. (2016), and the KiloDegree Survey (KiDS) by Martinet et al. (2018) and Shan et al. (2018).
For future large and deep lensing surveys, the importance of dealing with systematics increases significantly. To cite a few, potential sources of systematics are shape measurement errors (Cardone et al. 2014), photometric redshift uncertainty (Cunha et al. 2014; Bonnett et al. 2016; Choi et al. 2016; Gruen & Brimioulle 2017), intrinsic alignment of galaxies (Chisari et al. 2014; Codis et al. 2015; Schaefer & Merkel 2015; Schrabback et al. 2015; Krause et al. 2016), baryon physics (Mohammed et al. 2014; HarnoisDéraps et al. 2015), and instrumental responses (Gurvich & Mandelbaum 2016; Okura et al. 2016; Kannawadi et al. 2016; Plazas et al. 2016). Concerning WL peak counts, studies of systematics have been done by Yang et al. (2013) and Osato et al. (2015) for baryon physics, Liu et al. (2014) for masking, and Liu et al. (2014) for magnification bias (see also Lin 2016 for a review of WLpeakrelated studies). These existing studies are not sufficient for modeling the peak statistics at high precision.
One of the systematics that has not been addressed is the convergenceshear inversion. In WL, the (reduced) shear is observable but the convergence field is not. In order to obtain the convergence, one common way is to invert its relation to the shear via the lensing potential. This is called the KaiserSquires inversion (Kaiser & Squires 1993). It assumes that the galaxy shape noise could be modeled as a Gaussian random field, which conserves the same properties before and after the inversion. However, in practice, this hypothesis is rarely fulfilled for two reasons. First, the field is probed only on irregular galaxy samples. Second, masks introduce missing data and border effects. Because of these facts, the noise level is not spatially uniform. It can even be far from uniform for the galaxy density and the filter size used in realistic survey scenarios. In addition, inversion methods does not account for the reduced shear without nonlinear effects. As a result, the convergenceshear inversion can lead to a systematic bias to WL observables. Alternatively, Seitz & Schneider (1995) proposed a nonlinear method to properly take the reduced shear into account. However, this approach still relies on the assumption of uniform noise, which needs to be dealt with carefully.
Up to now, there exist three modeling approaches for predicting WL peak counts. The first approach consists in generating Nbody runs to simulate structures that cause lensing (for example, Kratochvil et al. 2010 and its series including Liu et al. 2015a). Then, peaksare counted from the derived lensing maps. The second approach is a pure analytical calculation, which computes the probability of different lensing signal levels given a line of sight (for example, Fan et al. 2010 and its series including Liu et al. 2015b). After that, the probability of identifying a peak can be deduced with random field theory. The third approach adopts a semianalytical model (this work). Based on some simple assumptions, it generates simulations in a much faster way than Nbody runs and count peaks from the resulting lensing maps.
In the literature, Fan et al. (2010) and Shirasaki (2017) have proposed theoretical models for WL peaks. These models directly predict peak counts on a convergence field from theory^{1}. Therefore, such approaches do not account for the inversion that is required for establishing a convergence field from data. The purpose of this paper is to quantify this modeling bias from theoretical WL peak models when neglecting the inversion. Two algorithms, the KaiserSquire and SeitzSchneider inversions, are considered in this study. We use our stochastic model developed in Paper I, Lin & Kilbinger (2015b, Paper II), and Lin et al. (2016, Paper III) to examine this effect. The stochasticity of the model makes including different inversions and different mapmaking methods straightforward. This makes our model an ideal tool to achieve our objective.
After the introduction, the theoretical formalism of different inversion methods is introduced in Sect. 2. Then, we explain the methodology, including the construction of different comparison cases and the choice of some detailed settings in Sect. 3. The results are presented in Sect. 4. Then, we will summarize the paper and conclude with a discussion in Sect. 5.
2 Inversion formalism
In this paper, we aim to test the impact of the convergenceshear inversion on WL peak counts. The most commonly used technique is the KaiserSquires (KS) inversion (Kaiser & Squires 1993) in which the convergence κ is given by (1)
where γ is the shear, the Fourier transform operator, and ℓ the Fourier mode. The arguments of κ and γ are omitted. To recover κ, the inverse transform leaves an undetermined constant term in direct space, corresponding to ℓ_{1} = ℓ_{2} = 0. This global constant is usually set to zero as the expected mean of the convergence field is null. In the WL regime, the reduced shear g ≡ γ∕(1 − κ), which has an unbiased estimator using the observed galaxy ellipticities (Seitz & Schneider 1997), is often approximated by the shear, so that g ≈ γ and Eq. (1) is applied on g directly.
In order to account for the factor 1∕(1 − κ) properly, Seitz & Schneider (1995) proposed an alternative inversion method (SS inversion). It follows the iterative process below: (2)
where the index ^{(i)} stands for the ith iterate, and . In the case where is always positive (i.e., g < 1), the second line of Eq. (2) is equivalent to γ^{(i)} = (1 − κ^{(i−1)})g. The KS and SS inversions will be the two methods to study in this paper.
3 Methodology
We used the CAMELUS code (proposed in Paper I) to model peakcount predictions. This semianalytical model adopted a halo approach, deriving peak number counts from a mass function. It first drew halos with mass and redshift with respect to the input mass function. Then, it randomized the halos’ angular positions. After that it assigned source galaxies and computed their respective lensing signal. Finally, the model created a lensing map and counted peaks by their S/N. In this way, massive darkmatter halos were considered as the major source of WL peaks.
Concerning how CAMELUS carried out halo sampling and galaxy assignment in this study, readers are invited to read Paper III for details. In the following, we explain how lensing maps were made and how peaks were defined. In order to test the impact from the convergenceshear inversion, we compared peak counts in three cases:

Case 1: κ is computed directly from raytracing simulations.

Case 2: g is computed from raytracing simulations and κ is given by the KS inversion.

Case 3: g is computed from raytracing simulations and κ is given by the SS inversion.
The κ and g could be given by a halo’s projected mass. Here, κ was computed following Eqs. (26) and (27) of Takada & Jain (2003a) and g by combining these previous equations with Eqs. (16) and (17) of Takada & Jain (2003b). The convergence fields of all three cases were generated from the same series of CAMELUS simulations, so that there was no statistical fluctuation caused by different halo or Galaxy samples. In this setup, Case 1 closely follows the logic of theoretical models from Fan et al. (2010) and Shirasaki (2017) where the inversion is omitted, whereas Cases 2 and 3 represented the realworld scenario. If the inversion caused a misestimation of peak counts, then a comparison of Case 1 to Case 2 or 3 should be able to separate this contribution from other systematic sources.
For this paper, we performed two sets of fast simulations. The first was to quantify the bias from mismodeling the inversion using the following diagnostic: , where was the number of peaks obtained directly from the convergence or from the shear with inversion. The denominator was the peaks from the inverted shear since this corresponded to the observation. Therefore, this indicator showed the peakcount deviation due to mismodeling. The choice of the denominator was due to the fact that the bias that this paper studies does not refer to “the bias caused by the inversion” but “the bias of mismodeling which does not account for the inversion”. Hereafter, we will use “mismodeling bias” as a shortened term to refer to the bias from mismodeling of the inversion. For this test, 2000 independent realizations were carried out for nine cosmologies varying Ω _{m} and σ_{8}. The dark energy equation of state was fixed at . A realization was a 450 × 450 map with a pixel width of 0.8 arcmin, such that the field area was 36 deg^{2}. Three Gaussian filters of widths 1.2, 2.4, and 4.8 arcmin were applied. A characteristic mask taken from the W1 field of CFHTLenS data was also applied.
The second set of fast simulations was to measure how cosmological constraints could be affected. Here, we chose to study cosmological constraints in a threedimensional space composed of Ω _{m}, σ_{8}, and . We mimicked the observational data vector by peaks from Case 2 under a reference cosmology , while the likelihood was computed in both Cases 1 and 2. By doing so, the difference between likelihoods could visualize how mismodeling would propagate to parameter constraints. This second set was the same fast simulations that were used in Paper III, which consisted of 37 536 cosmologies with 500 independent realizations each. The pixel size, field area, and smoothing filters were configured in the same way as the first simulation set.
Weaklensing peaks from Cases 2 and 3 were defined in the same way as in Paper III: shape noise had been added for all galaxies, which were binned later; after the inversion and the filtering, the noise level for S/N was estimated locally based on the effective number of neighboring galaxies covered under the filter, in the same way as Eq. (19) of Paper III. However, for Case 1, the simulations have been done differently. For Case 1, we first binned galaxies to generate a noisefree convergence map, then added a constant pixel noise before filtering the map. The noise variance is defined as (3)
where was the sum of the variance of the intrinsic ellipticity distribution, n_{gal} the source galaxy number density, and A_{pix} the area of a pixel. The motivation for this setting, though different from other cases, is that theoretical models are restricted to such a configuration. Only if this configuration is respected do we start to follow the same modeling logic as them. As a result, the definition of S/N should also change. Since the noise was constant for Case 1, the expected noise level in S/N was simply defined as (4)
where θ_{G} was the smoothing size.
From our cosmological parameter space, we projected constraint contours on all of the twodimensional plans to study numerically the constraint diagnostic. To take into account degeneracies, the constraint contours were fitted with three indicators,
which were respectively degeneracy lines for the Ω_{m} – σ_{8}, σ_{8} – , and Ω _{m} – planes. This definition of Σ _{8} allowed a good measurement of contour width independent from α. Concerning , Eqs. (6) and (7) assumed that was connected to two other considered parameters by an affine relation. An analysis of these indicators, already used in Paper III, would allow us to examine how degeneracy lines vary. Here, we set α = 2∕3, β = 0, a_{1} = 0.108, and a_{2} = 0.128 and computed the fitted Σ_{8}, I_{1}, and I_{2} respectively in Cases 1 and 2. This choice of α and β made Σ _{8} a functional form similar to ; and the choices of a_{1} and a_{2} were the bestfit results taken from Paper III.
4 Results
Figure 1 shows the peak function from the reference cosmology with the filter of 2.4 arcmin, using the first simulation set mentioned in Sect. 3. The error bars are rescaled to correspond to a survey of 1080 deg^{2}. We see that the biases of the convergence case compared to both shear cases are relatively small. They are globally contained within ±30% of Cases 2 and 3. The difference between the KS and SS inversions is even smaller, revealing that the inversion method has little importance for peaks. This figure also shows that the inversion decreases the peak counts under this configuration. However, we find that this cannot be generalized for all cosmologies and all filter sizes.
Interpreting the mechanism of this bias is challenging. One may think that the bias is caused by using g = γ∕(1 − κ) in the calculation instead of the true shear γ. If we follow this reasoning, for high peaks, the reduced shear g should be larger than γ since κ is positive. After a linear KS inversion, the deduced convergence should also be larger than the true convergence, so the peak function from Case 2 should be situated at the righthand side (which also means on top) of the one from Case 1. However, this is exactly the opposite of what Fig. 1 indicates, meaning that the reduced shear is not the origin of this bias. Actually, if we compare Case 3 (where the factor 1∕(1 − κ) has been taken into account) to Case 1, we see that the difference between two peak functions is still present. As a result, we suggest the irregularity of galaxy distribution to be the major origin of this bias.
We observe that high and low peaks do not seem to have the same variation. In Fig. 1, the convergence modeling underestimates the number of low peaks and overestimates high ones. However, their amplitude depends on the applied filter. This is shown more clearly in Table 1 where the peakcount bias in a specific bin ν ∈ [4.0, 4.5] with all three filter sizes and both inversion methods are presented. The values in the parentheses are the uncertainties of the biases estimated using the jackknife technique. We omitted one of the 2000 realizations at a time and computed the variance over these new derived subsamples. The jackknife errors confirm a bias of the order of 20% for the two first filters. Neglecting the inversion effect does not affect the bias in the same way for different filters, and we find the same ambiguity for other S/N bins which are not shown.
Table 1 also shows an ambiguous dependence on cosmology. It is not clear that high Ω _{m} or high σ_{8} would yield a larger or smaller bias. To visualize better the cosmological dependence of the bias, we computed the bias for each cosmology from the second fast simulation set (Sect. 3) and obtained a “bias map” for the KS inversion (Case 1 compared to Case 2). For example, the bias maps for the bin ν ∈ [4.0, 4.5] and the filter sizes of 1.2 and 2.4 arcmin are presented in Fig. 2. Here, we focus on the Ω _{m} – σ_{8} plane by fixing at 0.96. Although noisy, the contours of bias levels are visible, and within the studied range of cosmologies, the bias can vary from –4% to 38%. The level of noise in Fig. 2 is about two times larger than that of Table 1, since the former possesses four times less realizations than the latter. The figure shows the complex dependence of the peakcount bias on cosmology. Moreover, when we draw the bias maps for other bins and filter sizes, both contour shape and value range vary a lot. This variation is difficult to explain. It seems as if these contours still have a “banana” shape but are not centered on the same region at all. Therefore, the reference cosmology can be sometimes in the very biased region, sometimes very far from the very biased region. We argue that cosmology affects the mismodeling bias in a very complex way that depends on the filter size and the bin range.
Following the geometry of contours displayed in Fig. 2, we explored the link between the mismodeling bias and the reduced parameter Σ_{8}, defined by Eq. (5). Figure 3 shows the variation of the bias for different Σ _{8} when three filters are applied. It has the same S/N bin and the same inversion as the previous figure such that the blue dots in Fig. 3 correspond to the left panel in Fig. 2 and red pluses to the right panel. Despite a large scatter, we observe a clear pattern of variations. The filter size and the bin choice should have affected the amplitude of the bias and the value of Σ _{8} where the bias reaches its maximum. The scatter would have been reduced if α or β parameters had been chosen differently. However, the optimal values for one case are not necessarily the optimal ones for another. Therefore, we did not look for a minimization of the scatter in Fig. 3. As a remark, the artifactlike cut at Σ _{8} = 1.3 for all filters is due to the lack of sample points in the parameter space (gray zone in Fig. 2).
How does the mismodeling bias affect cosmological constraints? Figure 4 shows a comparison where solid and dashed contours are obtained from the inverted shear of Case 2 and shaded areas from Case 1. More precisely, contours correspond to a case where the inversion has been properly considered, while shaded areas are interpreted as coming from a simplified modeling that ignores the inversion effect. Although the bias in individual bins is relatively small as shown in Fig. 1, the impact on constraints is not negligible. While the Ω _{m} – σ_{8} constraint shifts approximately along the degeneracy line, both joint constraints with are more affected, resulting in a biased estimation of parameters. For example, the reference parameter has been excluded at 2σ on the leftbottom panel in Fig. 4.
The shift of contours observed above are quantified by the diagnostics mentioned in Sect. 3. The values of indicators are shown in Table 2. On the plane Ω _{m} – σ_{8}, Table 2 shows that both cases exclude each other mutually by 1.5–2σ, whereas the top panel in Fig. 4 does not suggest such a result. This is due to a bad fitting of Σ _{8} when the dragging parameter β is not relaxed. For this family of degeneracy lines, which have a hyperbolic form, there are at least two degrees of freedom. We could have relaxed both α and β to make Eq. (5) loyally describe the degeneracy lines. However, the optimal values of (α, β) are not necessarily the optimal ones for other cases, and these parameters have to be set to the same values, otherwise the comparison between Σ_{8} from different cases or surveys would be meaningless. Moreover, there is no physicallymotivated choice. Therefore, we were forced to fix α and β arbitrarily. Eventually, α = 2∕3 and β = 0 are coherent to most studies in the literature. We also tried the commonly used form for characterizing Σ_{8} and the same difficulty occurred. This explains why twodimensional contours between two cases (Fig. 4) appear more inconsistent with onedimensional intervals (Table 2).
From the two lowers panels in Fig. 4, we find that not considering the mismodeling bias does yield a 2σ tension compared to the realistic case, which justifies our finding that is more affected. Knowing that the bias depends on the S/N bin and the filter size in a nontrivial way, readers should keep in mind that another choice of the data vector might yield a contour shift different from what we have observed. In the end, the mismodeling bias is not completely negligible in cosmological constraints.
Fig. 1
Peak function from different studied cases. Blue circles, red squares, and green diamonds respectively show the mean peak density over 2000 realizations in Cases 1–3. The error bars are rescaled to correspond to a survey of 1080 deg^{2}. The filter size is 2.4 arcmin. The input cosmology is . We can see that the impact of not modeling the inversion is small. Between two inversion techniques, there is only a slight difference. 

Open with DEXTER 
Inversion mismodeling bias obtained under different cosmologies and different inversion methods.
Fig. 2
Bias map under a wide range of cosmological parameters. For these two panels, the KS inversion is applied. The dark energy equation of state is . The S/N bin is ν = [4.0, 4.5]. The filter size is 1.2 arcmin for the left panel and 2.4 arcmin for the right panel. Contours of parameters sharing the same bias level can be clearly observed. However, they are not necessarily centered on the same parameters when the filter changes. The bias levels are not similar either. These contours vary even more when we examine other choices of bins. 

Open with DEXTER 
Fig. 3
Bias as a function of Σ_{8} for different filters. Blue dots, red pluses, and green crosses are respectively the bias of a specific S/N bin [4.0, 4.5] under different cosmologies for a Gaussian filter of 1.2, 2.4, and 4.8 arcmin. Instead of visualizing on the Ω _{m} – σ_{8} plane, here we chose a onedimensional representation by the reduced parameter Σ _{8}. As a result, blues dots contain exactly the same information as the left panel in Fig. 2 and red pluses the same as the right panel. 

Open with DEXTER 
Fig. 4
Impact of the inversion mismodeling bias on Ω_{m} – σ_{8} – constraints. Solid and dashed lines represent contours derived from the inverted shear, whereas shaded areas represent contours derived from the convergence. Although the bias in individual bins is relatively small, the shift of contours is not negligible. 

Open with DEXTER 
Diagnostics of cosmological constraints.
5 Summary and discussion
In this paper, we performed fast simulations with CAMELUS to quantify the effect of neglecting the convergenceshear inversion in WL peakcount modeling. We called it “inversion mismodeling bias” since it is introduced when theoretical models fail to include a bias caused by the inversion while it is present in real data. To quantify this effect, on the one hand, we simulated the shear signal as from observations and applied two inversion methods; on the other hand, we simulated directly the convergence signal to count the number of peaks as analytical models do. The comparison of the two yielded an estimation of how theoretical peak models would deviate from the truth.
We have found that not accounting for the inversion has a relatively weak effect on WL peak counts, and that the difference between the KS and SS inversion methods is also small. For the reference cosmology , the bias from neglecting the inversion is contained within about 30%. Modeling this bias could be challenging. Its dependency on the S/N value and the filter size is not trivial. Its dependency on cosmology seems to have a link to the reduced parameter Σ _{8}, but this link cannot be summarized as a large scatter caused by the choice of the slope α of Σ _{8}.
By comparing cosmological constraints obtained from two different cases, we examined the propagation of the mismodeling bias during the process of estimating cosmological parameters. We have found that the Ω_{m} – σ_{8} contours are less affected. The change is mainly along the degeneracy line. However, for the dark energy equation of state , we have seen that both on the Ω_{m} – and σ_{8} – planes, contours can exclude each other by 2σ. The dark energy parameter seems to be more affected.
For the future WL surveys aimed at constraining cosmological parameters with high precision, modeling any possible systematic sources, including the inversion, will be indispensable. For cosmic shear, the inversion can be bypassed since linking the twopointcorrelation functions of the shear to theory is straightforward. For peaks, unless we only focus on the very high S/N regime, the link is less trivial. A possible option to bypass the inversion would be the aperture mass (Kaiser et al. 1994; Schneider 1996), which allows us to identify tangential shear peaks. However, we should remember that the aperture mass does not account for the nonlinearity of the reduced shear. Either usingthe aperture mass or the inversion, modeling the correction should be a common goal for future studies.
Acknowledgements
For their financial support, the authors would like to thank the DIMACAV thesis fellowship of the Région d’ÎledeFrance and the French national program for cosmology and galaxies (PNCG), and the computing support from the in2p3 Computing Centre. We are also thankful to the anonymous referee for highlighting ambiguities of writing and for her/his timeinvestment. We also thank Zuhui Fan and Xiangkun Liu for discussions about this paper. Linc is grateful to François Lanusse and Sandrine Pires for discussions on the buildup of tests used in this paper, to Marie Gay for her endless help on technical issues, and to LiWen Huang, Corentin Champseix, PinTing Lin, and YuJung Sun for their aid in Linc’s relocation to Edinburgh.
References
 Bartelmann, M., Perrotta, F., & Baccigalupi, C. 2002, A&A, 396, 21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bonnett, C., Troxel, M. A., Hartley, W., et al. 2016, Phys. Rev. D, 94, 042005 [NASA ADS] [CrossRef] [Google Scholar]
 Cardone, V. F., Martinelli, M., Calabrese, E., et al. 2014, MNRAS, 439, 202 [NASA ADS] [CrossRef] [Google Scholar]
 Chisari, N. E., Mandelbaum, R., Strauss, M. A., Huff, E. M., & Bahcall, N. A. 2014, MNRAS, 445, 726 [NASA ADS] [CrossRef] [Google Scholar]
 Choi, A., Heymans, C., Blake, C., et al. 2016, MNRAS, 463, 3737 [NASA ADS] [CrossRef] [Google Scholar]
 Codis, S., Gavazzi, R., Dubois, Y., et al. 2015, MNRAS, 448, 3391 [NASA ADS] [CrossRef] [Google Scholar]
 Cunha, C. E., Huterer, D., Lin, H., Busha, M. T., & Wechsler, R. H. 2014, MNRAS, 444, 129 [NASA ADS] [CrossRef] [Google Scholar]
 Dietrich, J. P., & Hartlap, J. 2010, MNRAS, 402, 1049 [NASA ADS] [CrossRef] [Google Scholar]
 Fan, Z., Shan, H., & Liu, J. 2010, ApJ, 719, 1408 [NASA ADS] [CrossRef] [Google Scholar]
 Gruen, D., & Brimioulle, F. 2017, MNRAS, 468, 769 [NASA ADS] [CrossRef] [Google Scholar]
 Gurvich, A., & Mandelbaum, R. 2016, MNRAS, 457, 3522 [NASA ADS] [CrossRef] [Google Scholar]
 Hamana, T., Takada, M., & Yoshida, N. 2004, MNRAS, 350, 893 [NASA ADS] [CrossRef] [Google Scholar]
 HarnoisDéraps, J., van Waerbeke, L., Viola, M., & Heymans, C. 2015, MNRAS, 450, 1212 [NASA ADS] [CrossRef] [Google Scholar]
 Jain, B., & Van Waerbeke L. 2000, ApJ, 530, L1 [NASA ADS] [CrossRef] [Google Scholar]
 Kacprzak, T., Kirk, D., Friedrich, O., et al. 2016, MNRAS, 463, 3653 [NASA ADS] [CrossRef] [Google Scholar]
 Kaiser, N., & Squires, G. 1993, ApJ, 404, 441 [NASA ADS] [CrossRef] [Google Scholar]
 Kaiser, N., Squires, G., Fahlman, G., & Woods, D. 1994, in Clusters of Galaxies, eds. F. Durret, A. Mazure, & J. Tran Thanh Van, 269 [Google Scholar]
 Kannawadi, A., Shapiro, C. A., Mandelbaum, R., et al. 2016, PASP, 128, 095001 [NASA ADS] [CrossRef] [Google Scholar]
 Kilbinger, M. 2015, Rep. Prog. Phys., 78, 086901 [NASA ADS] [CrossRef] [Google Scholar]
 Kratochvil, J. M., Haiman, Z., & May, M. 2010, Phys. Rev. D, 81, 043519 [NASA ADS] [CrossRef] [Google Scholar]
 Krause, E., Eifler, T., & Blazek, J. 2016, MNRAS, 456, 207 [NASA ADS] [CrossRef] [Google Scholar]
 Kruse, G., & Schneider, P. 1999, MNRAS, 302, 821 [NASA ADS] [CrossRef] [Google Scholar]
 Kruse, G., & Schneider, P. 2000, MNRAS, 318, 321 [NASA ADS] [CrossRef] [Google Scholar]
 Laureijs, R., Amiaux, J., Arduini, S., et al. 2011, ArXiv eprints [arXiv:1110.3193] [Google Scholar]
 Lin, C.A. 2016, PhD thesis, CEA Saclay, France, ArXiv eprints [arXiv:1612.04041] [Google Scholar]
 Lin, C.A., & Kilbinger, M. 2015a, A&A, 576, A24 (Paper I) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lin, C.A., & Kilbinger, M. 2015b, A&A, 583, A70 (Paper II) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lin, C.A., Kilbinger, M., & Pires, S. 2016, A&A, 593, A88 (Paper III) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Liu, J., Haiman, Z., Hui, L., Kratochvil, J. M., & May, M. 2014, Phys. Rev. D, 89, 023515 [NASA ADS] [CrossRef] [Google Scholar]
 Liu, J., Petri, A., Haiman, Z., et al. 2015a, Phys. Rev. D, 91, 063507 [NASA ADS] [CrossRef] [Google Scholar]
 Liu, X., Wang, Q., Pan, C., & Fan, Z. 2014, ApJ, 784, 31 [NASA ADS] [CrossRef] [Google Scholar]
 Liu, X., Pan, C., Li, R., et al. 2015b, MNRAS, 450, 2888 [NASA ADS] [CrossRef] [Google Scholar]
 Liu, X., Li, B., Zhao, G.B., et al. 2016, Phys. Rev. Lett., 117, 051101 [NASA ADS] [CrossRef] [Google Scholar]
 LSST Science Collaboration (Marshall, P., et al.) 2017, ArXiv eprints [arXiv:1708.04058] [Google Scholar]
 Martinet, N., Schneider, P., Hildebrandt, H., et al. 2018, MNRAS, 474, 712 [NASA ADS] [CrossRef] [Google Scholar]
 Mohammed, I., Martizzi, D., Teyssier, R., & Amara, A. 2014, MNRAS, submitted [arXiv:1410.6826] [Google Scholar]
 Okura, Y., Petri, A., May, M., Plazas, A. A., & Tamagawa, T. 2016, ApJ, 825, 61 [NASA ADS] [CrossRef] [Google Scholar]
 Osato, K., Shirasaki, M., & Yoshida, N. 2015, ApJ, 806, 186 [NASA ADS] [CrossRef] [Google Scholar]
 Plazas, A. A., Shapiro, C., Kannawadi, A., et al. 2016, PASP, 128, 104001 [NASA ADS] [CrossRef] [Google Scholar]
 Reblinsky, K., Kruse, G., Jain, B., & Schneider, P. 1999, A&A, 351, 815 [NASA ADS] [Google Scholar]
 Schaefer, B. M., & Merkel, P. M. 2015, MNRAS, submitted [arXiv:1506.07366] [Google Scholar]
 Schneider, P. 1996, MNRAS, 283, 837 [NASA ADS] [CrossRef] [Google Scholar]
 Schrabback, T., Hilbert, S., Hoekstra, H., et al. 2015, MNRAS, 454, 1432 [NASA ADS] [CrossRef] [Google Scholar]
 Seitz, C., & Schneider, P. 1995, A&A, 297, 287 [NASA ADS] [Google Scholar]
 Seitz, C., & Schneider, P. 1997, A&A, 318, 687 [NASA ADS] [Google Scholar]
 Shan, H., Liu, X., Hildebrandt, H., et al. 2018, MNRAS, 474, 1116 [NASA ADS] [CrossRef] [Google Scholar]
 Shirasaki, M. 2017, MNRAS, 465, 1974 [NASA ADS] [CrossRef] [Google Scholar]
 Takada, M., & Jain, B. 2003a, MNRAS, 340, 580 [NASA ADS] [CrossRef] [Google Scholar]
 Takada, M., & Jain, B. 2003b, MNRAS, 344, 857 [NASA ADS] [CrossRef] [Google Scholar]
 Wang, S., Khoury, J., Haiman, Z., & May, M. 2004, Phys. Rev. D, 70, 123008 [NASA ADS] [CrossRef] [Google Scholar]
 Yang, X., Kratochvil, J. M., Wang, S., et al. 2011, Phys. Rev. D, 84, 043529 [NASA ADS] [CrossRef] [Google Scholar]
 Yang, X., Kratochvil, J. M., Huffenberger, K., Haiman, Z., & May, M. 2013, Phys. Rev. D, 87, 023511 [NASA ADS] [CrossRef] [Google Scholar]
The model of Fan et al. (2010) can estimate peaks directly from a shear field, but it is computationally expensive.
All Tables
Inversion mismodeling bias obtained under different cosmologies and different inversion methods.
All Figures
Fig. 1
Peak function from different studied cases. Blue circles, red squares, and green diamonds respectively show the mean peak density over 2000 realizations in Cases 1–3. The error bars are rescaled to correspond to a survey of 1080 deg^{2}. The filter size is 2.4 arcmin. The input cosmology is . We can see that the impact of not modeling the inversion is small. Between two inversion techniques, there is only a slight difference. 

Open with DEXTER  
In the text 
Fig. 2
Bias map under a wide range of cosmological parameters. For these two panels, the KS inversion is applied. The dark energy equation of state is . The S/N bin is ν = [4.0, 4.5]. The filter size is 1.2 arcmin for the left panel and 2.4 arcmin for the right panel. Contours of parameters sharing the same bias level can be clearly observed. However, they are not necessarily centered on the same parameters when the filter changes. The bias levels are not similar either. These contours vary even more when we examine other choices of bins. 

Open with DEXTER  
In the text 
Fig. 3
Bias as a function of Σ_{8} for different filters. Blue dots, red pluses, and green crosses are respectively the bias of a specific S/N bin [4.0, 4.5] under different cosmologies for a Gaussian filter of 1.2, 2.4, and 4.8 arcmin. Instead of visualizing on the Ω _{m} – σ_{8} plane, here we chose a onedimensional representation by the reduced parameter Σ _{8}. As a result, blues dots contain exactly the same information as the left panel in Fig. 2 and red pluses the same as the right panel. 

Open with DEXTER  
In the text 
Fig. 4
Impact of the inversion mismodeling bias on Ω_{m} – σ_{8} – constraints. Solid and dashed lines represent contours derived from the inverted shear, whereas shaded areas represent contours derived from the convergence. Although the bias in individual bins is relatively small, the shift of contours is not negligible. 

Open with DEXTER  
In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.