temet.collab package

Submodules

temet.collab.bartlett module

James Bartlett / MAGIC mission proposal

magicCGMEmissionMaps(single=False, subhaloID=None)

Emission maps (single, or in stacked M* bins) for MAGIC-II proposal.

Parameters:

single (bool) – render a single subhalo per mass bin for visualization (no stacking).
subhaloID (int) – if not None, just visualize this specific subhalo (single must be True).

magicCGMEmissionTrends(): Emission summary statisics (auxCat-based) as a function of galaxy properties, for MAGIC-II proposal.

temet.collab.explore module

Misc exploration plots and testing, checks for others.

amyDIGzProfiles(): Use some projections to create the SB(em lines) vs z plot.

martinSubboxProj3DGrid(): Compute (i) 2D histo projection, (ii) 2D sphMap projection, (iii) 3D sphMap grid then projection (todo), of subbox gas.

auroraVoyage2050WhitePaper(): Create plots for Aurora’s ESA Voyage 2050 white paper.

smitaXMMproposal(): Dependence of OVII on sSFR at fixed mass.

nachoAngularQuenchingDens(): Variation of CGM gas density with azimuthal angle (for Martin Navarro+20).

nachoAngularQuenchingImage(): Images of delta rho/rho (for Martin Navarro+20).

omega_metals_z(metal_mass=True, hih2=False, mstar=False, mstarZ=False, hot=False, higal=False): Compute Omega_Q(z) for various components (Q). Rob Yates paper 2021.

abhijeetMgIISurfDens(): Test for Anand+ (2022).

xenoSNevo_profiles(): Xeno idealized SNe runs: density profiles vs time.

xenoSNevo_movie(conf=1): Xeno idealized SNe runs: render visualization frames for movie.

arjenMasses5kpc(): Explore Mtot_5kpc vs M*_5kpc.

yenting_vis_sample(redshift=1.0): For the raw TNG-Cluster halos (not in the virtual box), render some views of RIZ stellar composite and SFR, to identify rings like Yen-Ting is after.

benedetta_vis_sample(): For all TNG300-1 centrals at z=1, Mhalo > 5e13, plot stellar RIZ (observed-frame) composites and SFR maps, in a few projections. One plot per halo.

erica_tng50_sfrmaps(): Render some SFR surface density maps of TNG50 galaxies for Nelson, E.+2021 vs. 3D-HST paper.

temet.collab.peroux module

Celine Peroux

writeH2CDDFBand(): Use H2 CDDFs with many variations (TNG100) to derive an envelope band, f(N_H2) vs. N_H2, and write a text file. https://arxiv.org/abs/1909.08624 (Figures 6, 7, 9).

galaxyImageH2(): Metallicity distribution in CGM image: Klitsch+ (2019) paper figure. https://arxiv.org/abs/1909.08624 (Figure 8).

galaxyImage(): MgII emission image for Celine proposal 2022.

radialProfilesHIH2(): Compute stacked radial profiles of N_HI(b) and N_H2(b). https://arxiv.org/abs/2011.01935 (Figure 1).

numdensHIVsColumn(): Re-create Rahmati+ (2013) Fig 2. https://arxiv.org/abs/2011.01935 (Figure 1).

galaxyRotationAngles(): For Daniela filaments project.

temet.collab.peroux_azimuthal module

Azimuthal angle dependence of CGM properties (Peroux+, TNG50) https://arxiv.org/abs/2009.07809

singleHaloImage(sP, subhaloInd=440839, conf=0): Metallicity distribution in CGM image.

metallicityVsVradProjected(sP, shIDs=[440839], directCells=False, clean=False, distBin=None, ylim=None): Plot correlation of gas metallicity and vrad (i.e. mass flow rate) in projection by using the images directly.

metallicityVsTheta(sPs, dataField, massBins, distBins, min_NHI=[None], ptRestrictions=None, fullbox=False, nThetaBins=90, addObs=False, addEagle=False, sizefac=1.0, ylim=None, percs=[38, 50, 62], distRvir=False): Use some projections to create the Z_gas vs. theta plot. dataField : what to plot on y-axis massBins: list of 2-tuples, one or more M* bins distBins: list of 2-tuples, one or more impact parameter bins ptRestrictions: pass to gridBox(), e.g. {‘NeutralHydrogenAbundance’:[‘gt’,1e-3]} min_NHI: list of scalars, enforce minimum N_HI column of pixels to consider (log cm-2) fullbox: do global projections out to larger distance, otherwise fof-local scope distRvir: distBins are in units of rvir

stackedImageProjection(): Testing.

paperPlots(): Driver.

paperPlotsB(): Azimuthal angle variation of B field strength.

temet.collab.pillepich module

Annalisa Pillepich

stellarMergerContribution(sP): Analysis routine for TNG flagship paper on stellar mass content.

stellarMergerContributionPlot(): Driver.

tngMethods2_windPatterns(conf=1, pageNum=0): Plot gas streamlines (galaxy wind patterns), 4x2, top four from L25n512_0000 and bottom four from L25n512_0010 (Illustris model), matched. https://arxiv.org/abs/1703.02970 (Figure 5).

loop_patterns()

annalisa_tng50_presentation(setNum=0, stars=False): TNG50 presentation paper: face-on + edge-on combination, 5x5 systems.

gjoshi_clustermaps(conf=0, haloID=0): Author: A. Pillepich

apillepich_TNG50MWM31s_bubbles_top30(setType='top30', partType='gas', partField='P_gas', rotation='edge-on'): Author: A. Pillepich

pgalan_tng50_fornax(setType='fornax10', partType='stars', partField='coldens_msunkpc2'): Author: A. Pillepich

cengler_tng50_MWM31satellites(setNum=0, setType='MWM31satellites_selection', setMaps='stellarmass'): Author: A. Pillepich

apillepich_TNG50MWM31s_maps(setType='TNG50MWM31s', setMaps='stellarlight'): Author: A. Pillepich

temet.collab.rubin module

Kate Rubin / HST MST 2024 Proposal

hubbleMCT_gibleVis(conf=1): Visualization of CGM emission from a GIBLE or TNG50 halo.

hubbleMCT_emissionTrends(simname='tng50-1', cQuant=None): Hubble MST Proposal 2024 of Kate Rubin, and MAGIC-2 proposal.

hubbleMCT_emissionTrendsVsSim(): Combine results from above into a summary plot.

temet.collab.truong_xrayangular module

Angular dependence/anisotropy of x-ray emission of the CGM/ICM (Troung+21, TNG50). https://arxiv.org/abs/2109.06884 OVII and OVIII emission (Troung+22). https://arxiv.org/abs/2307.01277

stackedHaloImage(sP, mStarBin, conf=0, renderIndiv=False, median=True, rvirUnits=False, depthFac=1.0, stack2Dmaps=False): Stacked halo-scale image: delta rho/rho (for Martin Navarro+21) and x-ray SB (for Truong+21). Orient all galaxies edge-on, and remove average radial profile, to highlight angular variation.

stackedPropVsTheta(sP, mStarBin, distBins, conf=0, depthFac=1.0, stack2Dmaps=True, distRvir=False, nThetaBins=45): Stacked plot of quantity vs azimuthal angle.

singleHaloImage(): Quick test.

paperPlots(): Plots for Truong+21 x-ray emission angular dependence paper.

fullboxEmissionO8(): Create fullbox emission figure for Truong+22.

temet.collab.zooming2 module

“Zooming in on accretion” paper series (II) - Suresh+ 2019. http://arxiv.org/abs/1811.01949

check_box(snap)

visualize_halo(conf=1, quadrant=False, snap=None): Visualize single final halo of h2_L11_12_FP (boosted, sims.zooms2) at z=2.25.

visualize_compare_vs_normal(conf=1): Visualize single final halo of h2_L11_12_FP (boosted, sims.zooms2) vs h2_L11_FP (unboosted) at z=2.25.

phase_diagram_ovi()

phase_diagram_coolingtime()

phase_diagram_vs_L11()

phase_diagram_ovi_tng50_comparison()

figure1_res_statistics(conf=0): Figure 1: resolution statistics in mass/size for gas cells, comparing runs.

tracer_ambient_hot_halo(): Check the existence of an ambient/pre-existing hot halo at r<0.25rvir, as opposed to the possibility that all hot gas is arising from wind.

gas_components_time_evo(): Plot redshift evolution of total mass in different halo gas components (Fig 4).

gas_components_radial_profiles(): Compare CGM cool-dense and non-cool-dense radial profiles between runs.

mgii_radial_profile(): Compare MgII column density profiles.

hi_covering_frac(): Plot radial profiles of HI covering fractions.

Module contents

Collaborations / work for other people.