temp = None
print("\t ... done!")
# observed coordinates
else:
print("Observed coordinates not supported yet...")
exit(1)
trcat.compute_catalogueProperties(param.findDouble('boxside'))
##########################################################
# Generate chain-mesh of the input tracers catalogue
ChM = cbl.ChainMesh3D(2. * trcat.mps(), trcat.var(cbl.Var__X_),
trcat.var(cbl.Var__Y_), trcat.var(cbl.Var__Z_),
vdcat.Max(cbl.Var__Radius_))
##########################################################
# Finally building the cleaned catalogue
# sets the radius if not read from file:
if not param.findBool('Radius'):
limit = param.findVectorDouble('delta_r')
radii = [delta_r[1] for ii in range(vdcat.nObjects())]
vdcat.set_var(cbl.Var__Radius_, radii)
# sets the central density if not read from file:
if not param.findBool('centralDensity'):
vdcat.compute_centralDensity(trcat, ChM, trcat.numdensity(),
param.findDouble('ratio'))
trcat = temp
temp = None
print "\t ... done!"
# observed coordinates
else:
print "Observed coordinates not supported yet..."
exit(1)
trcat.compute_catalogueProperties(param.findDouble('boxside'))
##########################################################
# Generate chain-mesh of the input tracers catalogue
ChM = cbl.ChainMesh3D(2. * trcat.mps(), trcat.var(et._X_), trcat.var(et._Y_),
trcat.var(et._Z_), vdcat.Max(et._Radius_))
##########################################################
# Finally building the cleaned catalogue
# sets the radius if not read from file:
if not param.findBool('Radius'):
limit = param.findVectorDouble('delta_r')
radii = [delta_r[1] for ii in range(vdcat.nObjects())]
vdcat.set_var(et._Radius_, radii)
# sets the central density if not read from file:
if not param.findBool('centralDensity'):
vdcat.compute_centralDensity(trcat, ChM, trcat.numdensity(),
param.findDouble('ratio'))