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'))