def calculate_normalized_nonthermal_temperature_profile(input_data) :
d = input_data
normalized_T = {}
for hid in d['halo_ids'] :
Tnt = calculate_T_nonthermal(sigr=d['sigr'][hid],vr=d['vr'][hid],
sigt=d['sigt'][hid],vt=d['vt'][hid])
integrated_halo_props = CHP(physical_radial_profile=d['Rmid'][hid],
mass_enclosed_profile=d['Mtot'][hid],
delta=d['delta'],
aexp=d['aexp'])
Tdelta = calculate_T_normalization(Mvir=integrated_halo_props.calculate_virial_mass(),
delta=d['delta'],
aexp=d['aexp'])
Rscaled = d['Rmid'][hid]/integrated_halo_props.calculate_virial_radius()
normalized_T[hid] = Tnt/Tdelta
normalized_T[hid] = make_profile(x=Rscaled,y=normalized_T[hid])
return normalized_T
def calculate_normalized_nonthermal_temperature_profile(input_data):
d = input_data
normalized_T = {}
for hid in d['halo_ids']:
Tnt = calculate_T_nonthermal(sigr=d['sigr'][hid],
vr=d['vr'][hid],
sigt=d['sigt'][hid],
vt=d['vt'][hid])
integrated_halo_props = CHP(physical_radial_profile=d['Rmid'][hid],
mass_enclosed_profile=d['Mtot'][hid],
delta=d['delta'],
aexp=d['aexp'])
Tdelta = calculate_T_normalization(
Mvir=integrated_halo_props.calculate_virial_mass(),
delta=d['delta'],
aexp=d['aexp'])
Rscaled = d['Rmid'][
hid] / integrated_halo_props.calculate_virial_radius()
normalized_T[hid] = Tnt / Tdelta
normalized_T[hid] = make_profile(x=Rscaled, y=normalized_T[hid])
return normalized_T
def calculate_bulk_all_ratio_profile(input_data) :
d = input_data
mass_ratio = {}
for hid in d['halo_ids'] :
mass_ratio[hid] = d['M_gas_bulk'][hid]/d['M_gas_all'][hid]
mass_ratio[hid] = make_profile(x=d['Rscaled'][hid],y=mass_ratio[hid])
return mass_ratio
def calculate_circular_velocity_squared(input_data):
d = input_data
vcirc_squared = {}
for hid in d['halo_ids']:
vcirc_squared[hid] = calculate_vcirc2(mass_Msun=d['M_tot'][hid],
radius_kpc=d['Rmid'][hid])
vcirc_squared[hid] = make_profile(x=d['Rscaled'][hid],
y=vcirc_squared[hid])
return vcirc_squared
def calculate_circular_velocity_squared(input_data) :
d = input_data
vcirc_squared = {}
for hid in d['halo_ids'] :
vcirc_squared[hid] = calculate_vcirc2(mass_Msun=d['M_tot'][hid],
radius_kpc=d['Rmid'][hid])
vcirc_squared[hid] = make_profile(x=d['Rscaled'][hid],
y=vcirc_squared[hid])
return vcirc_squared
def calculate_normalized_density_profile(input_data) :
d = input_data
normalized_rho = {}
rho_delta = calculate_rho_normalization(delta=d['delta'],
aexp=d['aexp'])
for hid in d['halo_ids'] :
dv = calculate_dv(rin=d['Rin'][hid],rout=d['Rout'][hid])
rho = d['M_gas'][hid]/dv #Msun/kpc^3
normalized_rho[hid] = rho*d['units']/rho_delta
normalized_rho[hid] = make_profile(x=d['Rscaled'][hid],y=normalized_rho[hid])
return normalized_rho
def calculate_normalized_temperature_profile(input_data) :
d = input_data
normalized_T = {}
for hid in d['halo_ids'] :
Tdelta = calculate_T_normalization(Mvir=d['Mvir'][hid],
delta=d['delta'],
aexp=d['aexp'])
normalized_T[hid] = d['T_mw'][hid]*d['units']/Tdelta
normalized_T[hid] = make_profile(x=d['Rscaled'][hid],y=normalized_T[hid])
return normalized_T
def calculate_normalized_entropy_profile(input_data):
d = input_data
normalized_S = {}
for hid in d['halo_ids']:
Sdelta = calculate_S_normalization(Mvir=d['Mvir'][hid],
delta=d['delta'],
aexp=d['aexp'])
normalized_S[hid] = d['S_mw'][hid] * d['units'] / Sdelta
normalized_S[hid] = make_profile(x=d['Rscaled'][hid],
y=normalized_S[hid])
return normalized_S
def calculate_normalized_temperature_profile(input_data):
d = input_data
normalized_T = {}
for hid in d['halo_ids']:
Tdelta = calculate_T_normalization(Mvir=d['Mvir'][hid],
delta=d['delta'],
aexp=d['aexp'])
normalized_T[hid] = d['T_mw'][hid] * d['units'] / Tdelta
normalized_T[hid] = make_profile(x=d['Rscaled'][hid],
y=normalized_T[hid])
return normalized_T
def calculate_total_temperature_profile(input_data) :
d = input_data
T_tot = {}
for hid in d['halo_ids'] :
Ttot = calculate_Ttot(sigr=d['sigr'][hid],vr=d['vr'][hid],
sigt=d['sigt'][hid],vt=d['vt'][hid],
Tmw=d['T_mw'][hid])
T_tot_normalized[hid] = Ttot
T_tot_normalized[hid] = make_profile(x=d['Rscaled'][hid],
y=T_tot_normalized[hid])
return T_tot
def calculate_nonthermal_fraction_profile(input_data) :
d = input_data
fTnt = {}
for hid in d['halo_ids'] :
Tnt = calculate_T_nonthermal(sigr=d['sigr'][hid],vr=d['vr'][hid],
sigt=d['sigt'][hid],vt=d['vt'][hid])
Ttot = calculate_Ttot(sigr=d['sigr'][hid],vr=d['vr'][hid],
sigt=d['sigt'][hid],vt=d['vt'][hid],
Tmw=d['T_mw'][hid])
fTnt[hid] = Tnt/Ttot
fTnt[hid] = make_profile(x=d['Rscaled'][hid],
y=fTnt[hid])
return fTnt
def calculate_total_temperature_profile(input_data):
d = input_data
T_tot = {}
for hid in d['halo_ids']:
Ttot = calculate_Ttot(sigr=d['sigr'][hid],
vr=d['vr'][hid],
sigt=d['sigt'][hid],
vt=d['vt'][hid],
Tmw=d['T_mw'][hid])
T_tot_normalized[hid] = Ttot
T_tot_normalized[hid] = make_profile(x=d['Rscaled'][hid],
y=T_tot_normalized[hid])
return T_tot
def calculate_normalized_total_temperature_profile(input_data) :
d = input_data
T_tot_normalized = {}
for hid in d['halo_ids'] :
Tdelta = calculate_T_normalization(Mvir=d['Mvir'][hid],
delta=d['delta'],
aexp=d['aexp'])
Ttot = calculate_Ttot(sigr=d['sigr'][hid],vr=d['vr'][hid],
sigt=d['sigt'][hid],vt=d['vt'][hid],
Tmw=d['T_mw'][hid])
T_tot_normalized[hid] = Ttot/Tdelta
T_tot_normalized[hid] = make_profile(x=d['Rscaled'][hid],
y=T_tot_normalized[hid])
return T_tot_normalized
def calculate_normalized_nonthermal_temperature_profile(input_data):
d = input_data
normalized_Tnt = {}
for hid in d['halo_ids']:
Tnt = calculate_T_nonthermal(sigr=d['sigr'][hid],
vr=d['vr'][hid],
sigt=d['sigt'][hid],
vt=d['vt'][hid])
Tdelta = calculate_T_normalization(Mvir=d['Mvir'][hid],
delta=d['delta'],
aexp=d['aexp'])
normalized_Tnt[hid] = Tnt / Tdelta
normalized_Tnt[hid] = make_profile(x=d['Rscaled'][hid],
y=normalized_Tnt[hid])
return normalized_Tnt
def calculate_normalized_temperature_profile(input_data) :
d = input_data
normalized_T = {}
for hid in d['halo_ids'] :
integrated_halo_props = CHP(physical_radial_profile=d['Rmid'][hid],
mass_enclosed_profile=d['Mtot'][hid],
delta=d['delta'],
aexp=d['aexp'])
Tdelta = calculate_T_normalization(Mvir=integrated_halo_props.calculate_virial_mass(),
delta=d['delta'],
aexp=d['aexp'])
Rscaled = d['Rmid'][hid]/integrated_halo_props.calculate_virial_radius()
normalized_T[hid] = d['T_mw'][hid]*d['units']/Tdelta
normalized_T[hid] = make_profile(x=Rscaled,y=normalized_T[hid])
return normalized_T
def calculate_normalized_total_temperature_profile(input_data):
d = input_data
T_tot_normalized = {}
for hid in d['halo_ids']:
Tdelta = calculate_T_normalization(Mvir=d['Mvir'][hid],
delta=d['delta'],
aexp=d['aexp'])
Ttot = calculate_Ttot(sigr=d['sigr'][hid],
vr=d['vr'][hid],
sigt=d['sigt'][hid],
vt=d['vt'][hid],
Tmw=d['T_mw'][hid])
T_tot_normalized[hid] = Ttot / Tdelta
T_tot_normalized[hid] = make_profile(x=d['Rscaled'][hid],
y=T_tot_normalized[hid])
return T_tot_normalized