# cosmology parameters
H0      = 100.0
omega_m = 0.30

# optional, nside for healpix, default 64
healpix_nside = 64

# masking style, for quadrant cuts. "none", "sdss", "equatorial"
mask_style = "equatorial"

# shear style
#  "reduced": normal reduced shear shapes, source catalog rows are like
#      ra dec g1 g2 weight ...
#  "lensfit": for lensfit with sensitivities
#      ra dec g1 g2 g1sens g2sens weight ...
#  "metacal": for metacal with response tensor g1sens=R11, g2sens=R22, g12sens=(R12+R21)/2 
#      ra dec g1 g2 g1sens g2sens g12sens weight ...

shear_style = "reduced"

# source id style
#    "none": no source id in first column (default)
#   "index": integer source id in first column
#   "somcell": integer source id in first column, integer SOMCELL id in second column

sourceid_style = "index"

# sigma crit style
#   "point": point z for sources. Implies the last column in source cat is z
#  "interp": Interpolate 1/sigma_crit calculated from full P(z)).
#            Implies last N columns in source cat are 1/sigma_crit(zlens)_i
#  "sample": random sample of p(z). Implies the second-to-last column in
#            source cat is a mean-z estimate to be used for weighting and
#            last column in source cat contains a random sample from the p(z). 

sigmacrit_style = "point"

# units of shear, "shear", "deltasig" or "both".  Defaults to "deltasig". For pipeline
# mode with variable shear styles, "both" is recommended (in which case output format is
# fixed and allows calculating all necessary quantities, see below).
shear_units = "deltasig"

# weight style 
# only used in shear_units = "both"
#
# "uniform": calculate mean shear and mean sigma_crit^-1, weighting sources by the 
#            weight column in the input catalog only
# "optimal": weight sources by the product of weight colum in the input catalog and 
#            expectation value of sigma_crit^-1 for a minimum variance estimate (default)
#
# The output allows you to calculate (weighted) mean shear and DeltaSigma in both styles. 

weight_style = "optimal"


# number of logarithmically spaced radial bins to use
nbin = 21

# index of innermost and outermost radius bin for which xshear should print pairs to logfile 
# (optional, default is 0 for both which means no printing)
pairlog_rmin = 0
pairlog_rmax = 0
# maximum number of pairs logged in each radial bin 
# first pairs in that bin are printed - make sure to shuffle your shape catalog 
# if you need a representative sample; 0 (default) to print all pairs
pairlog_nmax = 0

# min and max radius (units default to Mpc, see below)
rmin = 0.02
rmax = 35.15

# units of radius (Mpc or arcmin). If not set defaults to Mpc; alternative is "arcmin"
r_units = "Mpc"

# demand zs > zl + zdiff_min
# optional, used for sigmacrit_style "point" and "sample" (on z_mean for the latter)
zdiff_min       = 0.2
# note that if using weight_style "uniform" and sigma_crit_style "sample" the
#zdiff_min is not enforced in selecting source galaxies

sigmacrit_style = "interp"

# zlens values for the 1/sigma_crit(zlens) values tabulated for each source
# note the elements of arrays can be separated by either spaces or commas
zlvals = [0.02 0.035 0.05 0.065 0.08 0.095 0.11 0.125 0.14 0.155 0.17 0.185 0.2 0.215 0.23 0.245 0.26 0.275 0.29 0.305 0.32 0.335 0.35 0.365 0.38 0.395 0.41]

r_units     = "arcmin"
shear_units = "shear"

r_units     = "arcmin"
shear_units = "deltasig"

shear_units = "both"