Документ взят из кэша поисковой машины. Адрес оригинального документа : http://hea-www.harvard.edu/AstroStat/Stat310_1516/vlk_IACHEC2016.pdf Дата изменения: Mon Feb 29 09:54:16 2016 Дата индексирования: Sun Apr 10 12:38:10 2016 Кодировка: Поисковые слова: р п р п р п р п р п р п р п р п р п р п р п р п

Calibration Uncertainties Working Group IACHEC 2016

Systematic errors in calibration are important, and must be dealt with, either by working to eliminate them, or by providing people with means to deal with them: these are the two main goals of this WG.

1


Schedule
·

Mar 1, WG meeting, 9:00am-10:45am IST
· · · · ·

Intro to WG and pyBLoCXS, Vinay Kashyap Intro to Cal Concordance, Herman Marshall Updates to XSPEC, Keith Arnaud [skype] Status of Cal Concordance Project, Yang Chen/Xufei Wang/Xiao-Li Meng [skype] Discussion

·

Mar 2, Improving Cross-Calibration Status, 9:45am-12:45pm IST
· · ·

Monte Carlo constraints on instrument calibration, Jeremy Drake NuSTAR and PyBlocks, Kristin Madsen [skype] Panel Discussion: what next?
2


Calibration has Uncertainties
·

The fundamental equation of observational astronomy C(i,j,k1,k2,tf,t;) = dt dxdy dE f(x,y,E,t;) R(t,tf) PSF(x,y,E;t) RMF(E,k;x,y,t) ARF(E;x,y,t)

·

Calibration analysis inverts the usual analysis method
·

Given ARF, RMF, PSF, evaluate expected model spectrum to compare with observed counts Given known model spectrum, compare with observed counts to evaluate ARF, RMF, PSF
3

·


Calibration has Uncertainties
·

How to find the uncertainties?

·

Once known, how to account for them?

·

And then how to minimize them?

4


Calibration has Uncertainties
·

How to find and tabulate the uncertainties?
·

MCCal

·

Once known, how to account for them?
·

pyBLoCXS

·

And then how to minimize them?
·

Concordance
5


pyBLoCXS
Vinay Kashyap (CXC/CfA)
David van Dyk, Hyunsook Lee, Jin Xu, Jeremy Drake, Pete Ratzlaff, Alanna Connors, et al.

MCMC scheme to incorporate defined calibration uncertainty into analysis

Simulated = Nominal + Bias + randomized components + residuals

6


fitting to simulated data f(;) = ­ e­ ()

7

-- Jin Xu and Shandong Min


fitting to simulated data f(;) = ­ e­ ()

p(|D,A)

7

-- Jin Xu and Shandong Min


fitting to simulated data f(;) = ­ e­ ()

p(|D,A) p(|D,Ai)
7
-- Jin Xu and Shandong Min


fitting to simulated data f(;) = ­ e­ ()

p(|D,A) p(|D,Ai)

p(A) p(|D,A)

7

-- Jin Xu and Shandong Min


fitting to simulated data f(;) = ­ e­ ()

p(|D,A) p(|D,Ai)

p(A) p(|D,A)

p(A,|D)

7

-- Jin Xu and Shandong Min


fitting to simulated data f(;) = ­ e­ ()

p(|D,A) p(|D,Ai)

p(A) p(|D,A)

p(A,|D) p(A('),|D)

7

-- Jin Xu and Shandong Min


pyBLoCXS resources
·

Lee et al. 2011, Accounting for Calibration Uncertainties in X-ray Analysis: Effective Areas in Spectral Fitting, ApJ 731, 126 [2011ApJ...731..126L] Xu et al. 2014, A Fully Bayesian Method for Jointly Fitting Instrumental Calibration and X-Ray Spectral Models, ApJ 794 97X [2014ApJ...794...97X] Sherpa (PragBayes): http://cxc.harvard.edu/sherpa/ahelp/ pyblocxs.html github (FullBayes): https://github.com/astrostat/pyblocxs tutorial from IACHEC 2014: http://hea-www.harvard.edu/AstroStat/ Demo/pyBLoCXS/IACHEC2014/
8

·

·

· ·