예제 #1
0
def products(output):
    """Standardized output products for any XSPEC fit"""
    pdf(output + ".pdf", cmd="ldata delchi")
    wdata(output + ".qdp")
    xs_utils.dump_fit_log(output + ".log")
    xs_utils.dump_dict(xs_utils.fit_dict(), output + ".json")
예제 #2
0
def products(output):
    """Standardized output products for any XSPEC fit"""
    pdf(output + ".pdf", cmd="ldata delchi")
    wdata(output + ".qdp")
    xs_utils.dump_fit_log(output + ".log")
    xs_utils.dump_dict(xs_utils.fit_dict(), output + ".json")
예제 #3
0
# Fit and dump fit information/plot
# ---------------------------------

xs.Fit.renorm()
xs.Fit.perform()

if args.free_en:
    for cname, en in zip(instr.componentNames, lines):
        comp = instr.__getattribute__(cname)
        comp.LineE.frozen = False
        comp.Sigma.frozen = False
    xs.Fit.perform()


xs.Plot.xAxis = "keV"
xs.Plot.xLog = True
xs.Plot.yLog = True
xs.Plot.addCommand("co 2 on 2")
xs.Plot.addCommand("rescale y 1e-4 1")

xs.Plot.device = f_plot + "/cps"
xs.Plot("ldata resid delchi")
call(["ps2pdf", f_plot, f_plot + ".pdf"])
call(["rm", f_plot])

# Save current XSPEC state
dump_dict(fit_dict(), f_fit)
dump_fit_log(f_log)