Determining which halos are clusters

Jasmine Fitzsimons
  • 13 Feb

Hello! I have written a code to extract which halos are categorized as clusters:

h = 0.6774
m200c_msun = halos["Group_M_Crit200"] * 1e10 / h

mask = m200c_msun > 1e14
cluster_IDs = np.where(mask)[0]

Am I utilizing Group_M_Crit200 correctly? Also I am using a mass threshold of 1e14 to categorize a cluster - is this common practice?

Dylan Nelson
  • 13 Feb

Your question is a bit too vague for me to help. It depends on what you want to do. It's more of a science question, so I think you should ask your supervisor/collaborator.

Your code (i.e. unit conversions) are correct.

The use of M200c (vs other options), and 1e14 (vs other options) is up to you. These are, you're right, common choices.

  • Page 1 of 1