I am writing to ask some questions on looking for the mock image for a given halo.
I am using sunpy to extract fits files from illustris database, so as is shown in the instruction, I need to input the subhaloID to specify the image. Say here is a Halo (haloID = 2), how could I find the FOF image of this halo? I assume that the subhaloID of the this image should be the subhaloID of the first subhalo of this halo, so I use the code below to read out the subhaloID :
GroupFirstSub = il.groupcat.loadHalos(basePath,135,fields=['GroupFirstSub'])
When HaloID=2, GroupFirstSub=1030. But sunpy can not find the image of subhalo-1030. I think I must misunderstand something.
Could you please instruct me to figure it out?
PS: One more native question, "the SubfindDensity in dm class shows the the local total comoving mass density", the "total" means all the components ( dm+stars/wind+gas+.....), right? If I want to calculate the total density map of a halo, I can just load in the SubfindDensity then perform linear interpolation, do I understand right?
If you could help me out about these questions, I really appreciate it.
Thank you so much!
All the best,
Hi Nan Li,
Let me know if I can clarify anything further
Thank you so much for your reply, it is pretty clear to me now, and extremely helpful. really appreciate it.