def phat_data(inputs, object_mass=None):
vsfhs, vsfh_kws = prepare_vsfh_run(inputs)
nruns = len(vsfhs) * inputs.nsfhs
if nruns > 10:
os.system('ipcluster start -n=%i' % inputs.nprocs)
time.wait(45)
calibrate.main(vsfhs, vsfh_kws=vsfh_kws)
os.system('ipcluster stop')
opt_files, ir_files = load_data_files(inputs)
for i, vsfh in enumerate(vsfhs):
if nruns < 10:
if inputs.dry_run is False:
res_dict = vsfh.vary_the_SFH(object_mass=object_mass, dry_run=inputs.dry_run)
vsfh.write_results(res_dict)
pl = Plotting(vsfh)
gal_kw = {'filetype': inputs.data_ftype, 'angst': False, 'hla': False}
opt_gal = rsp.Galaxy(opt_files[i], filter1='F475W', filter2='F814W',
**gal_kw)
ir_gal = rsp.Galaxy(ir_files[i], filter1='F110W', filter2='F160W',
**gal_kw)
pl.compare_to_gal(opt_gal, ir_gal, 28, 28, narratio=False, ylim=(.1, 1e5))
