# This notebook has been added in 2021.
# I am using old CLUES data to estimate the center of mass peculiar velocity.
import numpy as np
import matplotlib.pyplot as plt
import generateInfo as gi
sim_id = ["10909", "16953", "2710"]
data_path = "../data/"
for sim in sim_id:
FOF_file=data_path+"IC{}_halos_catalog.dat".format(sim)
ID_file=data_path+"index_LG_IC{}".format(sim)
output_file=data_path+"pair_physical_vweb_IC{}.dat".format(sim)
gi.generate_info(FOF_file, ID_file, output_file=output_file, hubble=0.73)
physical_data = np.loadtxt(output_file)
vcm = physical_data[2]
print('VCM', vcm)
1 894 674 selected in total 1 VCM 357.821 1 1106 888 selected in total 1 VCM 448.458949 1 830 642 selected in total 1 VCM 356.634919
357.821