[yt-users] Newbie: SFR for halos

[Nick Gnedin]

Folks,

Could someone advise me how to compute SFR for halos?

I have a dataset for which I computed SFR:

d = 
yt.load('/data/gnedin/REICMP/D/RT/CASE.REF2/rei10_a0.1003/rei10_a0.1003.art')
ad = d.all_data()
sfr = StarFormationRate(d,data_source=ad)

I also know how to add stellar masses to halos:

sm = d.add_deposited_particle_field(("STAR","MASS"),"sum")
h = yt.load('/data/gnedin/REICMP/D/RT/CASE.REF2/HC/a=0.1/hc.0.h5')
hc = 
HaloCatalog(data_ds=d,halos_ds=h,output_dir="/data/gnedin/REICMP/D/RT/CASE.REF2/HC/a=0.1")
hc.add_recipe("calculate_virial_quantities", ["radius", sm], 
overdensity_field=('artio', 'HVAR_GAS_DENSITY'))
hc.create()

But how do I add sfr.Msol_yr field from sfr as a deposited particle 
field for d?

Thank you so much,

Nick Gnedin