I am trying to study the clustering of AGNs using the Halo Occupation Distribution formalism. I am trying to calculate the Mean Occupation Function of AGNs and their redshift evolution.
My basic workflow requirements are as follows :
Please let me know if anyone knows how to handle this. Also I would like to hear out any other suggestions that could help me.
Is your objective an observational comparison, or a theoretical study?
If you are intending primarily to compare to observations, I would encourage you to select "central versus satellite AGNs" in the same way as the data (this will have nothing to do with subhalos, flags, fuzz, etc).
If you are just doing a theory study, I would suggest to write down a reasonable definition of central vs satellite AGN and apply that to all PartType5 particles.
The simplest definition would be: for each FoF halo, call its "main" SMBH the central, and any other satellites. Any SMBHs outside of FoF halos could also be called centrals.
By "main" you could use (a) most massive, as you describe above, (b) most massive within the central subhalo, (c) most central, (d) etc.