this is for holding javascript data
Elizabeth edited wlenc.tex
about 10 years ago
Commit id: dc2a89a214f5653f70e41716416fe3abb243ec20
deletions | additions
diff --git a/wlenc.tex b/wlenc.tex
index d9624fa..c7c85fb 100644
--- a/wlenc.tex
+++ b/wlenc.tex
...
% mn2esample.tex
%
% v2.1 released 22nd May 2002 (G. Hutton)
%
% The mnsample.tex file has been amended to highlight
% the proper use of LaTeX2e code with the class file
% and using natbib cross-referencing. These changes
% do not reflect the original paper by A. V. Raveendran.
%
% Previous versions of this sample document were
% compatible with the LaTeX 2.09 style file mn.sty
% v1.2 released 5th September 1994 (M. Reed)
% v1.1 released 18th July 1994
% v1.0 released 28th January 1994
\documentclass[useAMS,usegraphicx,usenatbib]{mn2e}
\usepackage[letterpaper,top=1.5cm,bottom=2cm,left=2cm,right=2cm]{geometry}
%\usepackage[a4paper,top=1.5cm,bottom=3cm,left=1.5cm,right=1.5cm]{geometry}
\usepackage{graphicx}
\usepackage{psfig}
\usepackage{amsfonts,amssymb,amsmath}
% If your system does not have the AMS fonts version 2.0 installed, then
% remove the useAMS option.
%
% useAMS allows you to obtain upright Greek characters.
% e.g. \umu, \upi etc. See the section on "Upright Greek characters" in
% this guide for further information.
%
% If you are using AMS 2.0 fonts, bold math letters/symbols are available
% at a larger range of sizes for NFSS release 1 and 2 (using \boldmath or
% preferably \bmath).
%
% The usenatbib command allows the use of Patrick Daly's natbib.sty for
% cross-referencing.
%
% If you wish to typeset the paper in Times font (if you do not have the
% PostScript Type 1 Computer Modern fonts you will need to do this to get
% smoother fonts in a PDF file) then uncomment the next line
% \usepackage{Times}
%%%%% AUTHORS - PLACE YOUR OWN MACROS HERE %%%%%
\newcommand{\hmsun}{h^{-1}{\rm M}_\odot}
\newcommand{\hmpc}{h^{-1}{\rm Mpc}}
\newcommand{\hkpc}{h^{-1}{\rm kpc}}
\newcommand{\msun}{{\rm M}_\odot}
\newcommand{\mpc}{{\rm Mpc}}
\newcommand{\kpc}{{\rm kpc}}
\newcommand{\kms}{{\rm ~km~s^{-1}}}
\hyphenation{lu-mi-no-si-ty}
\hyphenation{co-va-rian-ce}
\DeclareMathOperator\arctanh{arctanh}
\DeclareMathOperator\arcosh{arcosh}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\title[Low X-ray Luminosity Galaxy Clusters. N: Weak Lensing Mass Reconstruction at 0.18 $\lt$ z $\lt$ 0.70]{Low X-ray Luminosity Galaxy Clusters. N: Weak Lensing Mass Reconstruction at 0.18 $\lt$ z $\lt$ 0.70}
\author[Elizabeth Johana Gonzalez, et al]{\parbox[t]{\textwidth}{\vspace{-1cm}
Elizabeth Johana Gonzalez$^{1,2,3}$\thanks{E-mail :[email protected]},
Mariano J. Dom\'{\i}nguez Romero$^{1,2,3,4}$,
Jos\'e Luis Nilo Castel\'on$^{5}$,
Diego Garc\'{\i}a Lambas$^{1,2,3}$,
}\vspace{0.2cm}\\
$^{1}$ Instituto de Astronom\'{\i}a Te\'orica y Experimental (CCT C\'ordoba, CONICET, UNC),
Laprida 854, X5000BGR, C\'ordoba, Argentina.\\
$^{2}$ Observatorio Astron\'omico de C\'ordoba, Universidad Nacional de C\'ordoba, Laprida 854, X5000BGR, C\'ordoba, Argentina.\\
$^{3}$ Consejo Nacional de Investigaciones Cient\'{\i}ficas y T\'ecnicas, Rivadavia 1917, C1033AAJ Buenos Aires, Argentina.\\
$^{4}$ Center for Gravitational Wave Astronmy, University of Texas at Brownsville, Fort Brown 80, Calvary 206, Brownsville, TX 78520, USA.\\
$^{5}$ Departamento de F´ısica, Facultad de Ciencias, Universidad de La Serena, Cisternas 1200, La Serena, Chile.\\
}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{document}
\date{Accepted XXX. Received March 2014; in original form XXX}
\pagerange{\pageref{firstpage}--\pageref{lastpage}} \pubyear{2002}
\maketitle
\label{firstpage}
\begin{abstract}
This is the Nth of a series of papers on low X-ray luminosity
galaxy clusters in which we present Weak Lensing Mass reconstructions
based on the photometric images ...
%
We develop a weak gravitational lensing analysis pipeline based on several
public domain packages ..
%
The Weak Lensing analysis pipeline is tested on a simulated images of
a galaxy cluster field developed for the DES survey, providing a recovery
of the (induced or injected) mass distribution with high ...
%
The high quality images allow us to reconstruct the projected mass density
profiles showing ...
%
Our results indicate ...
\end{abstract}
\begin{keywords}
dark matter -- clusters: galaxies -- X-rays.
\end{keywords}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Introduction:}
Weak and strong gravitational lensing are excellent tools to probe
the projected mass distribution of clusters, with strong lensing confined
to the central regions of clusters, whereas weak lensing can yield mass
measurements for larger radii.
%
Low luminosity X-ray clusters of galaxies have not been extensively studied
compared to more massive systems, nevertheless ..
% XXX citas XXX
Here we present a weak lensing analysis pipeline based on well stablished
techniques such as: ...
% XXX lista de soft usado XXX
Gravitational lensing analysis of low X-ray luminosity clusters of galaxies were
done by
% XXX check similar studies XXX
Low X-ray galaxy clusters could shed light on the assembly processes and
environmental effects on their galaxy population as those analyzed on paper I
% XXX check first paper Nilo conclusions XXX
In this work we adopt when necessary a standard cosmological model: $H_{0}$\,=\,70\,h$_{70}$\,km\,s$^{-1}$\,Mpc$
^{-1}$, \\$ \Omega_{m} $\,=\,0.3,
and $ \Omega_{\Lambda} $\,=\,0.7.
%
The plan of this paper is as follows. In Sec. \ref{sec:data}, we introduce
the details of the systems of galaxies observed ...
% XXX ... .
In Sec. \ref{sec:method}, we describe the WL pipeline methodology used ...
%XXX ..
and the methodology to estimate the corresponding errors.
%
In Sec. \ref{sec:results} we present and discuss the
results .
%
Finally, in Sec. \ref{sec:conclusions} we summarise the main results
of this work.
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Observational Dataset:}
\label{sec:data}
%
Our sample consist of selected seven galaxy clusters from the ROSAT PSPC
Pointed Observations (XXX Vikhlinin et al. 1998 and Mullis et al. 2003 XXX) for this project.
These had X-ray luminosities in the range $L_{X} \sim 0.5–45 × 10^{43} erg s^{-1}$,
angular core-radius smaller than 60 arcsec and redshifts between 0.18 to 0.70 (for
further details see Paper I, \citep{}).
%
The galaxy photometric properties of our clusters sample were obtained with
Gemini North and South telescopes using the Multi-Object Spectrograph (Hook et al.
2004, hereafter GMOS) in the Image mode were presented in paper II \citep{}.
%
% Discutir si ponemos la tabla 1 del paper del Nilo. En caso afirmativo
% copiar y referenciar ACA, XXX Tabla 1 incluir XXX.
% podemos incluir el límite en magnitud usado también si les parece
%
This detector is an array of three 2048 × 4608 pixels EEV
CCDs, and using a 2×2 binning, the pixel scale is 0.1454
arcsec per pixel which corresponds to a field of view (FOV)
of 5.5 arcmin2 of the sky. This strategy allowed us to study
about 1.2 Mpc of the cluster centre, regardless of the redshift.
%
The clusters were imaged through the \textit{r'} and \textit{g'}
Sloan filters (Fukugita et al. 1996) in optimal observing conditions
(seeing disk diameter $\thicksim$ 0.7''). Images were reduced using
GEMINI/IRAF package, following the standard procedure.\\
%
The galaxy population for our clusters of galaxies sample was
extensively analized in paper II \citep{}, At redshift > 0.4 a less
important Red Cluster Sequence with a smaller extension was found compared
against the low redshift sample, but with a clear presence of the blue population.
The total galaxy radial density profiles were well fitted by a single power law.
%
For the three galaxy clusters with spectroscopic measurements,
% XXX poner cúales XXX
no significant detection of substructures through
the Dressler Shectman Test \citep{} analysis was found.
%
The fractions of elliptical, lenticular and blue galaxies
were calculated as a function of cluster redshift, confirming an
increasing fraction of blue galaxiesfrom 0.1 to 0.5 with redshift,
a fact related to the Butcher-Oemler effect. As redshift increases,
the fraction of lenticulars decreases from about 0.4 to 0.1, whereas
the fraction of early-types remaines almost constant at about 0.35.
\citep{Nilo2}
% esta textual del paper de Nilo 2, si les parece poner esto,
% podemos reescribirlo.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Weak Lensing Analysis:}
\label{sec:method}
%
%
We perform a pipeline based on python language to make the
lensing analysis. Its detects and classifies the sources, determines the \textit{Point Spread Function} (henceforth PSF) at the position of every object, measures the shape of the galaxies accounting for the PSF deformation, selects the background galaxies and computes the shear profile. \\
In order to check the pipeline funcionality, we apply it to the \textit{DES Cluster Simulation} images, publically available (Gill et al. 2009).\\
The steps in the weak lensing analysis and the results of the aplication of it to the simulated data, are described in the next subsections.
%
\subsection{Object detection and classification}
%
We used SExtractor (Bertin \& Arnouts, 1996) for the detection and fotometry of the
sources, in a two-pass mode. A first run is made to detect bright objects, with a detection level of 5$ \sigma $ above the background, in order to estimate the seeing and the saturation level of each image. The seeing is estimated using the average of the FWHM parameter of the point-like objects, selected from the FWHM/MAG\_BEST diagram since for these objects the FWHM is independent of the magnitude (\textbf{se entiende, o les parece mejor agregar el grafico?)}). Determining the seeing is important for SExtractor to perform the star-galaxy classification. The saturation level is estimated as 0.8 times the maximum value of the FLUX\_MAX parameter. These parameters, \textit{seeing} and \textit{saturation level}, are taking into account in the SExtractor configuration file for the second run, with a lower threshold detection limit of 1.5$ \sigma $. Second run is made in dual mode, detecting objects on the \textit{r'} image, while astrometric and photometric parameters are measured on all individual images.\\
%
\begin{figure}
\centering
\includegraphics[width=.5\textwidth]{./../plots/new-ones/mu_mag.eps}~\\
\includegraphics[width=.5\textwidth]{./../plots/new-ones/mag_fwhm.eps}
\caption{Objects detected in the galaxy cluster [VMF98]102 \textit{r'-}band image, stars (\textit{blue points}), galaxies (\textit{green points}), and artifacts (\textit{black points}). In the MU\_MAX/MAG\_BEST plane (\textit{left}), stars are in the region marked by the red solid line $\pm$ 0.4 magnitudes. To check the selection, we plot the sources in the MAG\_BEST/FWHM plane (\textit{left}), where point like object selected as stars have almost the same FWHM despite the magnitude.}
\label{sources}
\end{figure}
%
For the object classification in stars, galaxies and false detections,
we considered a similar analysis as Bardeau\,et\,al.\,2005, taking into account the position of the source in the magnitude/central flux diagram, the FWHM respect to the seeing and the stellarity index, according to the CLASS\_STAR parameter. Objects that are more sharply peaked than the PSF (FWHM $<$ \textit{seeing} - 0.5 pixel) and with FLAG parameter $>$ 4, are considered as false detections. As the light distribution of a point source scales with magnitude, objects on the line magnitude/central flux, $\pm$ 0.4 magnitudes (Figure\,\ref{sources}), FWHM $<$ \textit{seeing} + 1 pixel and CLASS\_STAR $>$ 0.8 are considered as stars. The rest of the objects are considered as galaxies. \\
%
\subsection{Shape measurements}
%
For the shape measurements we used IM2SHAPE (Bridle et al. 2002). Its computes
the shapes parameters modeling the object as sums of Gaussians, convolved with a
PSF which is also a sum of Gaussians. For simplicity, both the PSF and the object
are modeled with a single elliptical Gaussian profile. \\
The PSF was determined measuring the shapes of the stars as they are intrinsically pointlike objects. Once we obtained the shape of the stars, we clean the catalogue removing those objects with ellipticity greater than 0.2, which mainly appear to be false detections or faint galaxies, and by looking at the 5 nearest stars at each star position and removing those that differ by more than 2$\sigma$ from the local average shape. \\
After cleaning the catalogue, we linearly interpolate the local PSF at each object position by averaging the shapes of the five closest stars. We checked that
our Im2shape implementation can recover point-like objects by applying this PSF correction to each star. Figure\,\ref{PSF} shows the results of the shape parameters measurements for these stars, with and wihtout taking into account the PSF in the shape measurement: the size distribution is dominated by point sources, and the orientation is more uniformly distributed after the PSF correction.\\
\begin{figure}
\centering
\includegraphics[width=.25\textwidth]{./../plots/psf/ab.eps}~\hfill
\includegraphics[width=.25\textwidth]{./../plots/psf/ab_correc.eps}
\includegraphics[width=.25\textwidth]{./../plots/psf/map.eps}~\hfill
\includegraphics[width=.25\textwidth]{./../plots/psf/map_correc.eps}
\includegraphics[width=.25\textwidth]{./../plots/psf/e1.eps}~\hfill
\includegraphics[width=.25\textwidth]{./../plots/psf/e1_correc.eps}
\includegraphics[width=.25\textwidth]{./../plots/psf/e2.eps}~\hfill
\includegraphics[width=.25\textwidth]{./../plots/psf/e2_correc.eps}
\includegraphics[width=.25\textwidth]{./../plots/psf/Theta.eps}~\hfill
\includegraphics[width=.25\textwidth]{./../plots/psf/Theta_correc.eps}
\caption{PSF treatment applied to stars of one of the images of the DES
simulations: Semiaxis ($a$ cos $\theta$, $a$ sin $\theta$) before (\textit{left}) and after (\textit{right}) the
PSF deconvolution, which is more randomly distributed and considerably smaller (First row at top left corner
indicates the scale).}
\label{PSF}
\end{figure}
\subsection{Shear radial profiles}
Gravitational lensing maps the unlensed image, specified by coordinates $(x,y)$, to the lensed image $(x',y')$ using a matrix transformation:
\\
\begin{equation*}
\begin{pmatrix}
x\\
y
\end{pmatrix}
\begin{pmatrix}
1-g_{1} & -g_{2} \\
-g_{2} & 1+g_{1}
\end{pmatrix}
=
\begin{pmatrix}
x' \\
y' \\
\end{pmatrix}
\end{equation*}
\\
where $g_{1}$ and $g_{2}$ are the components of the reduced shear:
\begin{equation*}
g=\dfrac{\gamma}{1-\kappa}
\end{equation*}
If lensing is weak, the image of a circular source with ratio \textit{r}, appears elliptical, with axis given by
\begin{equation*}
a=\dfrac{r}{1-\kappa-\gamma} , b=\dfrac{r}{1-\kappa+\gamma}
\end{equation*}
Defining the ellipticity as
\begin{equation*}
e=\dfrac{a-b}{a+b}=\dfrac{\gamma}{1-\kappa}\approx\gamma
\end{equation*}
where \textit{g} becomes the normal shear, $\gamma$, since $\kappa \ll 1$, which generally holds in the weak lensing regime for clusters, and will be assumed henceforth here.\\
If the source has an intrinsic ellipticity $e_{s}$, the observed ellipticity in the weak lensing limit will be:
\begin{equation*}
e=e_{s}+\gamma
\end{equation*}
Assuming that unlensed galaxies are randomly oriented on the sky plane ($\langle e_{s} \rangle = 0$ ) and averaging over sufficiently many sources:
\begin{equation*}
\langle e \rangle=\langle \gamma \rangle
\end{equation*}
Hence, in the weaklensing approximation, we get an unbiased estimator of the reduced shear by averaging the shape of background galaxies in concentric annuli around the cluster center. Spherical symmetry also implies that the tangential component (E-mode) of the lensed galaxies traces the reduced shear, while the average of the component tilted at $\pi/4$ relative to the tangential component is the 'B-mode' and should be exactly zero for the case of perfect symmetry (see e.g. Sec. 4 of Bartelmann \& Schneider 2001 or Bernstein \& Jarvis 2002).\\
Because of the random orientation of the galaxies in the source plane, the error in the observed galaxy ellipticities and thus, on the estimated shear, will depend on
the number of galaxies averaged together to measure the shear (Schneider et al. 2000). The error can be estimated as:
\begin{equation*}
\sigma_{\gamma}\approx\dfrac{\sigma_{\epsilon}}{\sqrt{N}}
\end{equation*}
where $\sigma_{\epsilon}$ is the dispersion of the intrinsic ellipticity distribution ($\sigma_{\epsilon} \approx 0.3$) and $N$ is the number of objects in the annular bin.
\subsection{Background Galaxies selection and redshift distribution}
%
Background galaxies for the shear estimation, were selected as the galaxies with magnitude in the filter \textit{r'}, $m_{r}$, higher than $m_{L}$, and lower than $m_{max} + 1.0$, where $m_{L}$ is defined as the lowest magnitude $r'$ such that the probability that the galaxy is behind the cluster is higher than 0.7 and $m_{max}$ correspond to the peak of the magnitude distribution of the image in the filter \textit{r'} (both magnitudes for each cluster are listed in Table\,\ref{table:2}). The later cut in magnitude ensures that we are not taking into account galaxies that are too faint, given that they could have great errors in the shape measurements. Also, as we did for the simulated data, we discard galaxies with FWHM $<$ 5, with companions closer than 16 pixels and with $\sigma_{e} > 0.2$.\\
To compute $m_{L}$ we used the catalogue of photometric redshifts computed by Coupon et al. 2009, based on the public release Deep Field 1 of the Canada-France-Hawaii Telescope Legacy Survey. We estimated the fraction of galaxies with $z > z_{cluster}$ in magnitude bins for the \textit{r'} filter, and then we chose $m_{L}$ as the lowest magnitude for which the fraction of galaxies was greater than 0.7. Background galaxy density after the selection are listed for each cluster in Table 2.\\
Once we obtained a catalogue for the background galaxies, we average the two components of the ellipticities (E-mode and B-mode) in nonoverlapping logarithmic annuli, to conserve the S/N ratio in each bin. In order to take into account the contamination of foreground galaxies in the catalogue, we weighted the average value of the components of the galaxie's ellipticity with the probability that the galaxy was behind the cluster. We compute this probability using Coupon's catalogue, from the fraction of galaxies with $z > z_{cluster}$ for each bin in magnitude, \textit{r'}, and color (\textit{g' - r'} and \textit{r' - i'}) - Figure\,\ref{weigh}. Hence, given the magnitude and the color of each galaxy, we assigned to it a weigh, \textit{w}, as the fraction of galaxies with $z > z_{cluster}$ in that bin. \\
\begin{figure}
\centering
\includegraphics[scale=0.35]{./../plots/new-ones/distrib_z05.ps}\\
\caption{Fraction of galaxies with $z > 0.5$ ($n(z>0.5)/N_{tot}$), for different magnitudes in filter \textit{r'} and colors \textit{r'-i'}.}
\label{weigh}
\end{figure}
The average E-mode components correspond to the shear value which depends on the geometrical factor $ \beta = D_{LS}/ D_{L} $, where $D_{LS}$ is the angular diameter distance from the lens to the background source galaxy, and $D_{L}$ is the distance from the observer to the background galaxy. A galaxy at the same radial distance from the center of the cluster but at a different background redshift is sheared differently. To remove this variation, we divide the estimated shear by $\beta$, which results in the infinite redshift limit (z $\rightarrow \infty$) shear value, that is only function of the radial distance from the cluster center. To estimate $\beta$ we used once more Coupon's catalogue, since we can consider that this catalogue is complete up to our magnitude limit. We applied the photometric selection criteria to the catalogue and then we computed $\beta$ for the whole distribution of field galaxies, taking into account the contamination by foreground galaxies given our selection criteria by setting $\beta(z_{phot} < z_{cluster}) = 0$. The averaged geometrical factors for each cluster are given in Table\,\ref{table:2}.
\subsection{Fiting the profiles}
We finally estimate the cluster mass fitting the shear data with the singular isothermal
sphere (SIS) and the NFW profile (Navarro et al. 1997). We will explain briefly the lensing formulae for these two profiles:
\subsubsection{SIS profile}
The SIS mass model is the simplest one for describing a relaxed massive sphere with a constant and isotropic velocity dispersion, $\sigma_{V}$. This is mainly described by the density distribution:
\begin{equation*}
\rho(r) = \dfrac{\sigma_{V}^{2}}{2 \pi G r^{2}}
\end{equation*}
where $c$ is the speed of light and \textit{G} is the gravitational constant. From this ecuation, we can get the critical Einstein radius for source at $z \rightarrow \infty $ ($\beta =1$) as:
\begin{equation*}
\theta_{E} = \dfrac{4 \pi \sigma_{V}^{2}}{c^{2}}
\end{equation*}
In terms of wich one obtains:
\begin{equation*}
\kappa_{\theta} = \gamma_{\theta} = \dfrac{\theta_{E}}{2 \theta}
\end{equation*}
where $\theta$ is the distance to the cluster center. Hence, fitting the shear for different radius, we can estimate the Einstein radius, and from that, we can obtain an estimation of the mass within a radius $\theta$ as:
\begin{equation*}
M(\theta) = \dfrac{\theta_{E} c }{4G} \theta
\end{equation*}
\subsubsection{NFW profile}
The NFW profile is derived from fitting the density profile
of numerical simulations of cold dark matter halos (Navarro
et al. 1995, 1997). This profiles depend on two parameters, the virial radius, $R_{200}$, and a dimensionless concentration parameter:
\begin{equation*}
\rho(r) = \dfrac{\rho_{c} \delta_{c}}{(r/r_{s})(1+r/r_{s})^{2}}
\end{equation*}
where $r_{s}$ is the scale radius, $r_{s} = R_{200}/c$ and $\delta_{c}$ is the characteristic overdensity of the halo,
\begin{equation*}
\delta_{c} = \frac{200}{3} \dfrac{c^{3}}{\ln(1+c)+c/(1+c)}
\end{equation*}
The properties of the NFW profile in the context of gravitational lensing have been discussed by many authors (Bartelmann 1996, Wright \& Brainerd 2000). The radial dependence of the shear as function of $x:=\theta / r_{s}$ is given by:
\begin{equation*}
\gamma(x) = \dfrac{2r_{s}\delta_{c}\rho_{c}}{\Sigma_{c}} j(x)
\end{equation*}
where $\Sigma_{c}$ is the \textit{critical surface mass density} and $j(x)$ is given by: \begin{eqnarray*}
$j(x < 1) & = & \dfrac{4 \arctanh \sqrt{\frac{1-x}{1+x}}}{x^{2} \sqrt{1-x^{2}}} + \dfrac{2 \ln(\frac{x}{2})}{x^{2}} - \dfrac{1}{x^{2}-1} \\
& &+ \dfrac{2 \arctanh \sqrt{\frac{1-x}{1+x}}}{(x^{2}-1) \sqrt{1-x^{2}}} \\ $
\end{eqnarray*}
\begin{eqnarray*}
$j(x = 1) = 2 \ln \left(\dfrac{1}{2} \right) + \dfrac{5}{3}$\\
$j(x < 1) = \dfrac{4 \arctan \sqrt{\frac{1-x}{1+x}}}{x^{2} \sqrt{1-x^{2}}} + \dfrac{2 \ln(\frac{x}{2})}{x^{2}} - \dfrac{1}{x^{2}-1} + \dfrac{2 \arctan \sqrt{\frac{1-x}{1+x}}}{(x^{2}-1)^{3/2}}} $
\end{eqnarray*}
If we fit the shear for different radius we can have an estimation of the parameters $c$ and $R_{200}$. Once we obtain these parameters we can compute the mass within a radius $\theta$ as:
\begin{equation*}
M(\theta) = \int^{\theta}_0 2\pi \theta \Sigma(\theta) d\theta = \int^{\theta}_0 4\pi \theta \rho_{c} \delta_{c} r_{s} F(x) d\theta
\end{equation*}
where,
\begin{eqnarray*}
$F(x < 1) & = & \dfrac{1}{x^{2}-1} (1 - \dfrac{1}{\sqrt{1-x^{2}}} \arcosh \dfrac{1}{x} )$
$F(x = 1) = \dfrac{1}{3}$ \\
$F(x > 1) & = & \dfrac{1}{x^{2}-1} (1 - \dfrac{1}{\sqrt{x^{2}-1}} \arccos \dfrac{1}{x} )$
\end{eqnarray*}
There is a well-known degeneracy between the parameters $R_{200}$ and $c$ when fitting the shear profile in the weak lensing regime. This is due to the lack of information on the mass distribution near the cluster center and only a combination of strong and weak lensing can raise it and provide useful constraints on the concentration parameter. Since we do not have strong lensing modeling of the clusters in the sample,
we decided to fix the concentration parameter, $c = 5$, and fit the mass profile with only one free parameter, $R_{200}$.
\subsection{Testing the pipeline with simulated data}
%
To check the performance of our weak lensing analysis pipeline, we tested it on the DES cluster simulation images publically available \citep{2009arXiv0909.3856G}, with good results. This simulation consists of a sets of images, with different grades of difficulty, of sheared galaxies according to a SIS profile with a velocity dispersion of 1250\,km\,s$^{-1}$.\\
%
\begin{figure}{h!}
\centering
\includegraphics[scale=0.4]{./../plots/des/DES4_shearprofile.eps}\\
\includegraphics[scale=0.4]{./../plots/des/DES2_shearprofile.eps}\\
\includegraphics[scale=0.4]{./../plots/des/DES3_shearprofile.eps}\\
\caption{Shear blabla des.}
\label{DES}
\end{figure}
We applied our pipeline to three of the image files available, High Noise File, High Noise PSF Applied File and Low Noise PSF Applied File. For the PSF Applied files, we checked that our Im2shape implementation can recover point-like objects by applying the PSF correction to each star. The images contain only the sheared galaxies, hence all the galaxies detected were considered as background galaxies at $z=0.8$, which is the average of the redshifts of the galaxies. We cut the catalogue discarding the galaxies with FWHM $<$ 5, with companions closer than 16 pixels and with $\sigma_{e} > 0.2$, where $\sigma_{e}$ is defined as the quadratic sum of the errors $\sigma_{e1}$ and $\sigma_{e2}$ given by IM2SHAPE. These cuts are made in order to consider, for the shear measurements, insulated galaxies with good shape measurements. Shear profiles are shown in Figure\,\ref{DES}. For the most complex image that we treated, we obtained a deviation parameter of 1.6, defined as the number of $\sigma$ that the result is away from the input value of $\sigma_{V} = 1250 km/sec$, i.e. $\sigma = \dfrac{result - input}{error}$ .
\begin{table*}
\centering
\caption{Summary of the weak lensing analysis}\label{tab:esp}
\label{table:2}
\begin{tabular}{@{}ccccccccc@{}}
\hline
\hline
\rule{0pt}{1.05em}%
[VMF\,98] & Density of background galaxies & $m_{L}$ & $m_{max}$ & $\langle\beta\rangle$ & $\sigma_{V}$ & M$^{(1)}$($r < 0.5 Mpc$) & $R_{200}$ & M$^{(2)}$($r < 0.5 Mpc$) \\
Id. & (galaxies/arcmin$^{2}$) & & & & (km/sec) & $10^{14} M_{\odot} h_{70}^{-1}$ & Mpc & $10^{14} M_{\odot} h_{70}^{-1}$ \\
\hline
\rule{0pt}{1.05em}%
001 & 33 & 23.0 & 27.0 & 0.42 & - & - & - & - \\
022 & 15 & 20.7 & 26.0 & 0.63 & 600 $\pm$ 30 & 1.3 & 1.7 $\pm$ 0.2 & 1.2 \\
093 & 9 & 22.3 & 25.0 & 0.50 & - & - & - & - \\
097 & 31 & 23.0 & 27.0 & 0.44 & 580 $\pm$ 40 & 1.2 & 1.2 $\pm$ 0.4 & 0.83 \\
102 & 26 & 22.7 & 26.9 & 0.50 & 450 $\pm$ 40 & 0.75 & 1.1 $\pm$ 0.1 & 0.62 \\
119 & 18 & 24.0 $^{(3)}$ & 26.4 & 0.30 & 860 $\pm$ 30 & 2.7 & 1.4 $\pm$ 0.1 & 1.4 \\
124 & 26 & 19.5 & 26.7 & 0.72 & - & - & - & - \\
148 & 19 & 24.0 $^{(3)}$ & 26.9 & 0.27 & 860 $\pm$ 140 & 2.7 & 1.4 $\pm$ 0.7 & 1.4 \\
\hline
\end{tabular}
\\
\begin{tablenotes}{\begin{footnotesize}
\end{footnotesize}
\begin{small}
\begin{flushleft}
\item [] (1) From the fit of the SIS profile
\item [] (2) From the fit of the NFW profile
\item [] (3) Here $m_{L}$ is defined as the lowest magnitude $r'$ such that the probability that the galaxy is behind the cluster is higher than 0.6}
\end{flushleft}
\end{small}
\end{tablenote}
\end{table*}
% XXX Discuss the details of the shear recovered profiles and mention possible improvements
%
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Clusters Dark Matter Content:}
%
\label{sec:results}
From our weak lensing analysis, we estimated the mass of five of the eight clusters in the sample. For the clusters [VMF98]001, [VMF98]093 and [VMF98]124, the shear profiles obtained were compatible with a plane profile. In the case of [VMF98]093 the troubles in getting a good profile are mainly because of the low density of background galaxies given that $m_{max}$ is low for the images of this cluster. The drawbacks in the mass determination of [VMF98]001, are the presence of another group in the line of sight (Figure ...) and the possible presence of substructure in the cluster. Finally [VMF98]124 seems to be a relaxed cluster and there is no evidence of another group in the line of sight
%
\begin{figure*}
\centering
\includegraphics[width=1.0\textwidth,height=1.3\textwidth]{./../plots/shear-profiles/shear_profiles.ps}
\caption{Shear blabla c22 filter r.}
\label{fig5}
\end{figure*}
\ref{fig5}
%
% XXX How display this information in the best way XXX
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Conclusions:}
\label{sec:conclusions}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\citet{b9}.
\section*{Acknowledgments}
%
We thank the anonymous referee for their very useful comments that improved
the content and clarity of the manuscript.
%
This work was partially supported by the Consejo Nacional de Investigaciones Cient\'{\i}ficas y T\'ecnicas (CONICET, Argentina)
and the Secretar\'{\i}a de Ciencia y Tecnolog\'{\i}a de la Universidad Nacional de C\'ordoba (SeCyT-UNC, Argentina).
%
This research has made use of NASA's Astrophysics Data System and Cornell University arXiv repository.
%
\bibliographystyle{mn2e}
\bibliography{references}
% Produces the bibliography via % BibTeX.%
\appendix
\end{document}
