def calc_profile(base_id, base, snap_num, sub_id , interact=0): import os from read_halo import read_halo from center_halo import center_halo #from J_align_halo import J_align_halo from extract_galaxy import extract_galaxy from rotate_halo import rotate_halo #from visualize_halo import visualize_halo #from J_profile_halo import J_profile_halo from J_profile_halo2 import J_profile_halo from gal_halo_comparison import gal_halo_comparison from bdd_halo import bdd_halo from save_halo import save_halo save_dir = base_id+"/snap"+str(snap_num).zfill(3)+"/" halo = Halo() halo.base_id = base_id halo.base = base halo.snap_num = snap_num halo.sub_id = sub_id halo.save_dir = save_dir halo.interact = interact # Minimum threshold for # of particles type_thresh = 25 halo.type_thresh = type_thresh # Error flag halo.err_flag = 0 if False: read_halo(halo) center_halo(halo) extract_galaxy(halo) rotate_halo(halo) gal_halo_comparison(halo) #visualize_halo(halo) J_profile_halo(halo) bdd_halo(halo) save_halo(halo) if True: read_halo(halo) center_halo(halo) extract_galaxy(halo) rotate_halo(halo) J_profile_halo(halo) #print "Did rotation fail? ",halo.rotate_flag print "J_profile test done" return halo
def calc_profile(base_id, base, snap_num, sub_id , interact=0): import os from read_halo import read_halo from center_halo import center_halo #from J_align_halo import J_align_halo from extract_galaxy import extract_galaxy #from rotate_halo import rotate_halo #from visualize_halo import visualize_halo from J_profile_halo2 import J_profile_halo save_dir = base_id+"/" halo = Halo() halo.base_id = base_id halo.base = base halo.snap_num = snap_num halo.sub_id = sub_id halo.save_dir = save_dir halo.interact = interact # Minimum threshold for # of particles type_thresh = 25 halo.type_thresh = type_thresh # Error flag halo.err_flag = 0 read_halo(halo) center_halo(halo) extract_galaxy(halo) rotate_halo(halo) J_profile_halo2(halo) #print "Did rotation fail? ",halo.rotate_flag print "J_profile test done" return halo