BXA in practice (with pyxspec)

• Set up spectrum & model
• Set up parameter space
• Run!

s = Spectrum('example-file.fak')

s.ignore("**"); s.notice("0.2-8.0")

m = Model("pow")

m.powerlaw.norm.values = ",,1e-10,1e-10,1e1,1e1" # 10^-10 .. 10

m.powerlaw.PhoIndex.values = ",,1,1,3,3" # 1 .. 3

# define prior

transformations = [

bxa.create_uniform_prior_for( m, m.powerlaw.PhoIndex),

bxa.create_jeffreys_prior_for(m, m.powerlaw.norm),

]

bxa.standard_analysis(transformations, outputfiles_basename = 'simple-example')