Logged in as guest
Viewing Calibration/5635 Full headers
Private message: yes no
Notes: Notification:
Date: Mon, 5 Aug 2002 17:15:38 +0100 (BST) From: Igor Antokhin <igor@astro.gla.ac.uk> To: xmmhelp@xmm.vilspa.esa.es Subject: XMM PSF
Dear Sir/Madam, for my current XMM data analysis I need to know the encircled energy or PSF for MOS1, MOS2 and PN at various off-axis angles. Please note that it is NOT about correcting for the out-of-the-aperture contribution. I need to know the numbers (the fraction of the total PSF within a particular set of circular apertures). The xmmsas_20020605_1701/doc/arfgen/node21.html states that "For the common case of a circular or annular region centred on the source, a tabulation of enclosed energy as a function of energy, off-axis angle and box radius is used." I'd tried to find this table in the CCF directory but failed. Could you please point me to this file (if available) or to the way I can compute the numbers I need? Thanks, Igor Antokhin
From: Michel G. Breitfellner <xmmhelp@xmm.vilspa.esa.es> To: igor@astro.gla.ac.uk Subject: Re: XMM PSF (PR#5635) Date: Wed Aug 7 11:11:15 2002
Dear Dr. Antokhin, Please see the XMM-Newton Current Calibration File Release Note http://xmm.vilspa.esa.es/ccf/releasenotes/ XMM-CCF-REL-69 Encircled Energy Function of the X-ray telescopes 3 Apr 2001 http://xmm.vilspa.esa.es/docs/documents/CAL-SRN-0069-1-0.ps.gz With best regards, Michel G. Breitfellner XMM-Newton User Support Group
Date: Wed, 7 Aug 2002 14:31:41 +0100 (BST) From: Igor Antokhin <igor@astro.gla.ac.uk> To: "Michel G. Breitfellner" <xmmhelp@xmm.vilspa.esa.es> Subject: Re: XMM PSF (PR#5635)
This message is in MIME format. The first part should be readable text, while the remaining parts are likely unreadable without MIME-aware tools. Send mail to mime@docserver.cac.washington.edu for more info. ---779972051-8496021-1028727101=:1763 Content-Type: TEXT/PLAIN; charset=US-ASCII Dear Michel, thanks for your answer. Still it did not provide me with the means to calculate the encircled fraction of the PSF. I cannot access the CAL routines directly and the CAL manual (version 2.1, see the relevant extract from it attached to this message) does not really explain the algorithm of the subroutine CAL_getEncircledEnergy (which I hoped it would so I could reproduce it myself). In fact, it is puzzling. E.g, the formula for Rarcsec seems to be irrelevant as it is not explained how this parameter is used in the algorithm. In the files XRT[23]_XENCIREN_0003.CCF the parameters cor1 and cor2 are equal to zero. Does this mean that for MOS2 and PN no correction can be made for any off-axis angles? Even if so, I would still like to be able to calculate the numbers for MOS1. Thanks, Igor Antokhin On Wed, 7 Aug 2002, Michel G. Breitfellner wrote: > Dear Dr. Antokhin, > > Please see the XMM-Newton Current Calibration File Release Note > > http://xmm.vilspa.esa.es/ccf/releasenotes/ > > XMM-CCF-REL-69 Encircled Energy Function of the X-ray telescopes 3 Apr 2001 > > http://xmm.vilspa.esa.es/docs/documents/CAL-SRN-0069-1-0.ps.gz > > With best regards, > Michel G. Breitfellner > XMM-Newton User Support Group > > -------------------------------------------- Igor Antokhin Astronomy and Astrophysics Group University of Glasgow Kelvin Building, G12 8QQ Glasgow, UK E-mail: igor@astro.gla.ac.uk ---779972051-8496021-1028727101=:1763 Content-Type: APPLICATION/postscript; name="CAL_man_p15.ps" Content-Transfer-Encoding: BASE64 Content-ID: <Pine.LNX.4.44.0208071431410.1763@oberon.astro.gla.ac.uk> Content-Description: Content-Disposition: attachment; filename="CAL_man_p15.ps" JSFQUy1BZG9iZS0yLjAKJSVDcmVhdG9yOiBkdmlwc2sgNS41OGYgQ29weXJp Z2h0IDE5ODYsIDE5OTQgUmFkaWNhbCBFeWUgU29mdHdhcmUKJSVUaXRsZTog Y2FsaGIuZHZpCiUlUGFnZXM6IDEKJSVQYWdlT3JkZXI6IEFzY2VuZAolJUJv dW5kaW5nQm94OiAwIDAgNTk2IDg0MgolJURvY3VtZW50UGFwZXJTaXplczog QTQKJSVFbmRDb21tZW50cwolRFZJUFNDb21tYW5kTGluZTogZHZpcHMgLW8g Y2FsaGIucHMgY2FsaGIKJURWSVBTUGFyYW1ldGVyczogZHBpPTYwMCwgY29t bWVudHMgcmVtb3ZlZAolRFZJUFNTb3VyY2U6ICBUZVggb3V0cHV0IDIwMDEu MDcuMTY6MTAwNgolJUJlZ2luUHJvY1NldDogdGV4LnBybwovVGVYRGljdCAy NTAgZGljdCBkZWYgVGVYRGljdCBiZWdpbiAvTntkZWZ9ZGVmIC9Ce2JpbmQg ZGVmfU4gL1N7ZXhjaH1OCi9Ye1MgTn1CIC9UUnt0cmFuc2xhdGV9TiAvaXNs cyBmYWxzZSBOIC92c2l6ZSAxMSA3MiBtdWwgTiAvaHNpemUgOC41IDcyCm11 bCBOIC9sYW5kcGx1czkwe2ZhbHNlfWRlZiAvQHJpZ2lue2lzbHN7WzAgbGFu ZHBsdXM5MHsxIC0xfXstMSAxfQppZmVsc2UgMCAwIDBdY29uY2F0fWlmIDcy IFJlc29sdXRpb24gZGl2IDcyIFZSZXNvbHV0aW9uIGRpdiBuZWcgc2NhbGUK aXNsc3tsYW5kcGx1czkwe1ZSZXNvbHV0aW9uIDcyIGRpdiB2c2l6ZSBtdWwg MCBleGNofXtSZXNvbHV0aW9uIC03MiBkaXYKaHNpemUgbXVsIDB9aWZlbHNl IFRSfWlmIFJlc29sdXRpb24gVlJlc29sdXRpb24gdnNpemUgLTcyIGRpdiAx IGFkZCBtdWwKVFJbbWF0cml4IGN1cnJlbnRtYXRyaXh7ZHVwIGR1cCByb3Vu ZCBzdWIgYWJzIDAuMDAwMDEgbHR7cm91bmR9aWZ9CmZvcmFsbCByb3VuZCBl eGNoIHJvdW5kIGV4Y2hdc2V0bWF0cml4fU4gL0BsYW5kc2NhcGV7L2lzbHMg dHJ1ZSBOfUIKL0BtYW51YWxmZWVke3N0YXR1c2RpY3QgL21hbnVhbGZlZWQg dHJ1ZSBwdXR9QiAvQGNvcGllc3svI2NvcGllcyBYfUIKL0ZNYXRbMSAwIDAg LTEgMCAwXU4gL0ZCQlswIDAgMCAwXU4gL25uIDAgTiAvSUUgMCBOIC9jdHIg MCBOIC9kZi10YWlsewovbm4gOCBkaWN0IE4gbm4gYmVnaW4gL0ZvbnRUeXBl IDMgTiAvRm9udE1hdHJpeCBmbnRyeCBOIC9Gb250QkJveCBGQkIgTgpzdHJp bmcgL2Jhc2UgWCBhcnJheSAvQml0TWFwcyBYIC9CdWlsZENoYXJ7Q2hhckJ1 aWxkZXJ9TiAvRW5jb2RpbmcgSUUgTgplbmQgZHVwey9mb28gc2V0Zm9udH0y IGFycmF5IGNvcHkgY3Z4IE4gbG9hZCAwIG5uIHB1dCAvY3RyIDAgTlt9QiAv ZGZ7Ci9zZiAxIE4gL2ZudHJ4IEZNYXQgTiBkZi10YWlsfUIgL2Rmc3tkaXYg L3NmIFggL2ZudHJ4W3NmIDAgMCBzZiBuZWcgMCAwXQpOIGRmLXRhaWx9QiAv RXtwb3Agbm4gZHVwIGRlZmluZWZvbnQgc2V0Zm9udH1CIC9jaC13aWR0aHtj aC1kYXRhIGR1cApsZW5ndGggNSBzdWIgZ2V0fUIgL2NoLWhlaWdodHtjaC1k YXRhIGR1cCBsZW5ndGggNCBzdWIgZ2V0fUIgL2NoLXhvZmZ7CjEyOCBjaC1k YXRhIGR1cCBsZW5ndGggMyBzdWIgZ2V0IHN1Yn1CIC9jaC15b2Zme2NoLWRh dGEgZHVwIGxlbmd0aCAyIHN1YgpnZXQgMTI3IHN1Yn1CIC9jaC1keHtjaC1k YXRhIGR1cCBsZW5ndGggMSBzdWIgZ2V0fUIgL2NoLWltYWdle2NoLWRhdGEK ZHVwIHR5cGUgL3N0cmluZ3R5cGUgbmV7Y3RyIGdldCAvY3RyIGN0ciAxIGFk ZCBOfWlmfUIgL2lkIDAgTiAvcncgMCBOCi9yYyAwIE4gL2dwIDAgTiAvY3Ag MCBOIC9HIDAgTiAvc2YgMCBOIC9DaGFyQnVpbGRlcntzYXZlIDMgMSByb2xs IFMgZHVwCi9iYXNlIGdldCAyIGluZGV4IGdldCBTIC9CaXRNYXBzIGdldCBT IGdldCAvY2gtZGF0YSBYIHBvcCAvY3RyIDAgTiBjaC1keAowIGNoLXhvZmYg Y2gteW9mZiBjaC1oZWlnaHQgc3ViIGNoLXhvZmYgY2gtd2lkdGggYWRkIGNo LXlvZmYKc2V0Y2FjaGVkZXZpY2UgY2gtd2lkdGggY2gtaGVpZ2h0IHRydWVb MSAwIDAgLTEgLS4xIGNoLXhvZmYgc3ViIGNoLXlvZmYKLjEgc3ViXXtjaC1p bWFnZX1pbWFnZW1hc2sgcmVzdG9yZX1CIC9Eey9jYyBYIGR1cCB0eXBlIC9z dHJpbmd0eXBlIG5le119CmlmIG5uIC9iYXNlIGdldCBjYyBjdHIgcHV0IG5u IC9CaXRNYXBzIGdldCBTIGN0ciBTIHNmIDEgbmV7ZHVwIGR1cApsZW5ndGgg MSBzdWIgZHVwIDIgaW5kZXggUyBnZXQgc2YgZGl2IHB1dH1pZiBwdXQgL2N0
From: Michel G. Breitfellner <xmmhelp@xmm.vilspa.esa.es> To: igor@astro.gla.ac.uk Subject: Re: XMM PSF (PR#5635) Date: Wed Aug 7 16:07:48 2002
Dear Dr. Antokhin, Your question was forwarded to our PSF specialist. We will inform you as soon as possible. With best regards, Michel G. Breitfellner XMM-Newton User Support Group
From: Michel G. Breitfellner <xmmhelp@xmm.vilspa.esa.es> To: igor@astro.gla.ac.uk Subject: Re: XMM PSF (PR#5635) Date: Mon Aug 12 12:47:18 2002
Dear Igor, Please find enclosed the answer of our specialist. With best regards, Michel G. Breitfellner XMM-Newton User Support Group -------------------------------------------- Dear Igor, the latest formulation of the encircled energy is contained in the CCF files XRT1_XPSF_0005.CCF, XRT2_XPSF_0005.CCF and XRT3_XPSF_0004.CCF for the MOS-1, MOS-2 and PN cameras respectively. In these files you will find an extension called KING_PARAMS. This contains the columns 4 BINTABLE KING_PARAMS 8 16(3) 80 0 1 Column Name Format Dims Units TLMIN 1 ENERGY E eV 2 THETA E rad 3 PARAMS 2E Where the PARAMS column contains the King function parameters core radius (in arcseconds) and index (alpha). These parameters, tabulated against photon energy and off-axis angle (THETA in radians) can be converted to an encircled energy fraction using the formula given in the attached postscript file. In the SAS we interpolate in energy and theta using a spline but this isn't critical and a linear interpolation should be adequate. I hope this helps, Richard ------------------------------------
Date: Mon, 12 Aug 2002 14:10:04 +0100 (BST) From: Igor Antokhin <igor@astro.gla.ac.uk> To: "Michel G. Breitfellner" <xmmhelp@xmm.vilspa.esa.es> Subject: Re: XMM PSF (PR#5635)
On Mon, 12 Aug 2002, Michel G. Breitfellner wrote: > Dear Igor, > > Please find enclosed the answer of our specialist. ... > ... These parameters, > tabulated against photon energy and off-axis angle (THETA in radians) > can be converted to an encircled energy fraction using the formula > given in the attached postscript file. Thanks for the answer. Could you please send me the ps file which was NOT attached to your original message? Thanks again, Igor Antokhin
Date: Mon, 12 Aug 2002 16:33:43 +0200 From: "Michel G. Breitfellner" <mbreitfe@xmm.vilspa.esa.es> To: igor@astro.gla.ac.uk CC: xmmhelp@xmm.vilspa.esa.es Subject: XMM Helpdesk request 5635: XMM PSF (PR#5635)
This is a multi-part message in MIME format. --------------2E618022144926D997C8FD0C Content-Type: text/plain; charset=us-ascii Content-Transfer-Encoding: 7bit Dear Igor, Sorry, I forgot that our helpdesk system has difficulties handling attachments. Therefore I send you the file from my personal email account. Best regards, Michel G. Breitfellner XMM-Newton User Support Group -------------------------------------------- Dear Igor, the latest formulation of the encircled energy is contained in the CCF files XRT1_XPSF_0005.CCF, XRT2_XPSF_0005.CCF and XRT3_XPSF_0004.CCF for the MOS-1, MOS-2 and PN cameras respectively. In these files you will find an extension called KING_PARAMS. This contains the columns 4 BINTABLE KING_PARAMS 8 16(3) 80 0 1 Column Name Format Dims Units TLMIN 1 ENERGY E eV 2 THETA E rad 3 PARAMS 2E Where the PARAMS column contains the King function parameters core radius (in arcseconds) and index (alpha). These parameters, tabulated against photon energy and off-axis angle (THETA in radians) can be converted to an encircled energy fraction using the formula given in the attached postscript file. In the SAS we interpolate in energy and theta using a spline but this isn't critical and a linear interpolation should be adequate. I hope this helps, Richard ------------------------------------ --------------2E618022144926D997C8FD0C Content-Type: APPLICATION/octet-stream; name="calhb_eef.ps.gz" Content-Transfer-Encoding: BASE64 Content-Description: Extract from the next version of the CAL handbook Content-Disposition: attachment; filename="calhb_eef.ps.gz" H4sICO02Uj0AA2NhbGhiX2VlZi5wcwCkWm1v4ziS/s5fwcPdALPAJhYlx5Ib c4PuKNKigZ1cb7uvkUNjcFDbcuxrx/bIdjpZrf/71VNVlOWXZJKZJCIpslis N1YVqfzwbx8GZ+9Gi6/lWXgemB9+SKuyWC+qN3Z0P12uvtmL84tkbNPF8rGa 3k7W1vWT3l+p7Hftx2I0HRYzmz2WdrAYr78XVUkYPk3Xs/KNpZHJ13PCQl0f itty9cY6bf5XNSpphXerYTkfUd/lYjMfTee3l4uHNzag34t+zybdkIauFsPN XTlffyiWZTWY/hNo3nVpIJuP0sUdhlbmh6vP7z8M8FrMR3+fzkulvk0Dg3wo quKuXJcVYRktp/8ZBcFf7VDR2Kq8W9yXI4UdLDbVkDDZT+WNXWzWy83ahkHg zp07D6M3rtftg/Tydjr/UC2Gg3L9xq7Lh/NltTAdmnM1HdKEi8CO0BiVY+s7 v2KO7VzX1LnFQOey/jqdjwC0vbadQV0+DCfba9O5qQf2entpO58+1uuqmK9m xboEyHQ1W9lxMVuVlt7uVyQY65yNQ3u3maFrwl3J+QX1Ge2bkXSWs82qH9Q8 VdZ+S3qdzmtgrL8EtgXk7Jnb1mfOuq2ZjkssBuUEvw4X82Gx3k7HWPBjuVrM NuvpYk6c3qPn80HXvLy1K1JFaXiR1gptSJoIYOEFFAeWxVC3YM4EyEx2QFul 7NNH0NOCbeMWpDrbOluMRphtPn38clesq+mDHW6qioxA3urRZmnxVDBMu9p8 tcXXlQ3OA/pxdraueYAW3JrxoipmM4UEwa3mr6tSMUJnb8E4yWFZ1qK/dbUp oV7TeUuGuylm47Ic1at1sd6s2Gg6u24BJhuENbwdLpbTclV3/l0a9gZI8l+K 9RfHKjrj6ldaNL+8/BKo2uh1PqcW1e8zqYfrShqj8dm6mM5qA5BEbPbaUluN NV/M158el6WNAI23X0Rw4/m6evB9l7SDLS1pr82KRue3tvO1IOXc2KKqikcy 9CnNo415Q83NdDZKJ0VVo+C3soKcsvlwAW9gichrU2JfbJZ1Z7xYWJLnmJbZ hoqP2KfiHuvPFsWIWCGKsVE9Y18grtGY2FqNSe+EXeiFrOhNuRagVQ3jAOCN gn2hNkRHFWwYMjR7k7J6uVhiTTaYckyuB/Q1dBLIcHL2fTpaT2pqjIp1AVAz K+e364m9YNO6LT3gpISLbUNaheweQD4sxsSUCxN7AjjywKh1wiMmnIANAWMA 68K4BT96OAXtDsiY3pE394AGgGvYSEeUz+05jZMusAIrBY9swGvsH8IzHYkJ Vt9Rm041lPfbpdqo1n+TinVCdctq6lVxD8t0tPdoLw5YxmJ5WDe05FzLB24P djYor0KYsnpjoc9mV7AcTGBV4FblaBtVsThao6xpZk5BDVnCsBhOylF5Px2W O5gdCmzs3bY9O3cNQkXOeM5Z9L/WXuhbLu8KitBVuaKAXUKUV3VnOCQuntLE r3DksNedcAge7GLToH9fOBgZWN46NFu9otkzBnR6AYs0CRwbCb6KFmt2Y1vx ROr72tDS0nGFjq+LZb1ZlZV4Pno7mywW3+y3+eL7vPavwNgZvLes8GsrsctA dgjd64XtfLanPDqT6Zg8tBA52FcH7S4DWs77fTj4zj8+b+vOx88+vLB7QbyE gdADR1USxUSLyt964g0G2sSXLeJXk8X3JfSHgECuvlq3mOZ3BjUyc9fBooRt tsMa+am9N0fhiWXfoQyG7MPeGHj68KKn/vKaGHb0flFT9wAaFC8NAQS0gUSN UNnwnh3plsKb6V3EiTuPe4JaYumN3e+ViHwDtpY1mETr435EMhyCLq8WtCV7 gV8bO38zo2UD33yU5n2tbzce4sZ+hqV8NjWqj59tK1ZqnKLsa7Sh1++TsirZ O/uelhuLSet1YGP4vemc0sF7ymJp/Mer6Wo5Kx7/Uv7mzaNrWzA/Xpc3nzC4 qLaMW9IoTY62mlWpxVDgqutbNtRPH3lf0x+1HMSObMgKU8IjTzUQWA0Btbb3 rd/f25Po2kgELS1gJDV8Hp3SSZL8x2fFPC/JOElAnGuS6u/MybzGThsAv++8 CqsZxUB6D5QkipvG9wmMRFLfN56Sw75tubCiFowck8sZOWU2BkRcwx5Ie2+A ZAf6C2XKCMEYwjb2Efor9/u3YX3Wtb9sDXDXZxE1qVXWZ6G0xkh2uXVL5sGN Sa0901qB/q/Wed9q05XW93qPllm9tLTOd7Tv0I6kPUc7lPYCbSft32oD58Tt qoY75eaKmjpzTU1F+EBrDdqrPRI9oQS+pSnEk1JqOBiok2QWqYc6dmKWcxed ofTocnSSWS3L4bSY8Wnm4NwykJNNEPgscdAefTuQmVfluNjM1kTJZGV7LkRk v1/ZuB/y8WShe1zryWoo6dm9bxTzW43zf3//wScI00FZzk1z7qkm3LE7CM3K Ne3UGr5nvliX3Ci6XM3K22JGLU63eaNs5tM1XCYmvp1wV90aUbc2RHJNAPeH AEa9oQdgD8jc3ijRToT/9l5G7o9HhrPpknJT7gu1b0LhnrKGGjLSpX3Pve8h 4RAtLCPpINFQOqjOn17YZ4q4mmMGwCbTuoGaCNRkH4qkM5s91B0qFPds9ojX R33dVDRKRfP6iFcdJTcjJz1Gh3iJ/L3z9pcCZvXf7Qh/N9LwuNeFcIxsxXTQ xXlL+dsORIyM4HwOEshZjv6KB+0j3w5Md6N2shIKICPhUwPO2lMhjyDVeE/Z jxoOJf6dxahePywXK/hYyjjWD8ACK/dJxs4tDizlk/POEo5KvaqBttm93hXV t7qBrY+8Kbn6k6PiM7dmW/fgSHjP70Zf2FebaXtch1/aR1mVUkj524xEAYa2 W5R6EvYxZEjhYr1egFfNPgrJkzq3Q/swknSGtBH1KX9DykG6HP+T5Dpe0WJ5 51//yz/pYlNNywrnTEpP7stqPZ4Vt/WHgoLy5awYIjPCVYkopl4+rOzycaWh cLmsbMGpm6RtVT1fjEkLK47eA45FA0RRvOK6w/A8GMYD4Owjl3RuP5zgEnL5 C0o/Sp7XWi+CDRu0JN/k2w9KlFGjN+RetIKjcbnEaC1tyAxPLc/07kho1rWn 1w2MYBZ2BH+DeW/GPl08a9sI7dNHo+mvDeMT
Date: Tue, 13 Aug 2002 18:09:11 +0100 (BST) From: Igor Antokhin <igor@astro.gla.ac.uk> To: "Michel G. Breitfellner" <mbreitfe@xmm.vilspa.esa.es> cc: xmmhelp@xmm.vilspa.esa.es Subject: Re: XMM Helpdesk request 5635: XMM PSF (PR#5635)
Dear Michel, thanks for your response. Could you please clarify a couple of points so that I will not make any mistakes using the info you provided earlier? 1.In your message and in the postscipt description you say that the function used is the King function. In XRT1_XPSF_0005.CCF, the extention containing the parameters ENERGY, THETA and PARAMS is called GLE_GAUSS_PARAM (I could not find the KING_PARAMS extension). Is it still ok to use the formula provided in your postscript document (which I honestly doubt)? 2.The size of PARAMS is 6E, not 2E so there are 6 parameters for every offset angle and energy, not just r_c and alpha. What is their meaning and how can I use them? 3.How about the XRT[23]_XENCIREN_0003.CCF files? I has been told before that it is these files which are being used to correct for the fraction of the PSF outside of one's aperture. Still, cor1 and cor2 from those files are equal to 0 for MOS2 anf PN. This is not directly related to my current question but I would like to know whether my psf is being corrected if I use e.g. a 30arcsec aperture within 5arcmin from the center of the FOV and use appropriate task switches to make the correction. Thanks again, Igor Antokhin On Mon, 12 Aug 2002, Michel G. Breitfellner wrote: > Dear Igor, > > Sorry, I forgot that our helpdesk system has difficulties handling > attachments. Therefore I send you the file from my personal email > account. > > Best regards, > Michel G. Breitfellner > XMM-Newton User Support Group > > -------------------------------------------- > Dear Igor, > > the latest formulation of the encircled energy is contained in the > CCF files XRT1_XPSF_0005.CCF, XRT2_XPSF_0005.CCF and XRT3_XPSF_0004.CCF > for the MOS-1, MOS-2 and PN cameras respectively. > In these files you will find an extension called KING_PARAMS. > This contains the columns > > 4 BINTABLE KING_PARAMS 8 16(3) 80 0 1 > > Column Name Format Dims Units TLMIN > 1 ENERGY E eV > 2 THETA E rad > 3 PARAMS 2E > > Where the PARAMS column contains the King function parameters > core radius (in arcseconds) and index (alpha). These parameters, > tabulated against photon energy and off-axis angle (THETA in radians) > can be converted to an encircled energy fraction using the formula > given in the attached postscript file. > In the SAS we interpolate in energy and theta using a spline but this > isn't critical and a linear interpolation should be adequate. > > I hope this helps, > > Richard > ------------------------------------ -------------------------------------------- Igor Antokhin Astronomy and Astrophysics Group University of Glasgow Kelvin Building, G12 8QQ Glasgow, UK E-mail: igor@astro.gla.ac.uk
From: Matthias Ehle <xmmhelp@xmm.vilspa.esa.es> To: igor@astro.gla.ac.uk Subject: Re: XMM PSF (PR#5635) Date: Wed Aug 14 07:17:34 2002
Dear Igor, (Michel is on leave for 4 weeks now but I hope that me too, I can help you a bit with your questions on the PSF...) you wrote: > 1.In your message and in the postscipt description you say that the > function used is the King function. In XRT1_XPSF_0005.CCF, the extention > containing the parameters ENERGY, THETA and PARAMS is called > GLE_GAUSS_PARAM (I could not find the KING_PARAMS extension). Is it still > ok to use the formula provided in your postscript document (which I > honestly doubt)? hmmm, I just had a look at the XRT1_XPSF_0005.CCF file myself and I see the extension called KING_PARAMS (and no GLE_GAUSS_PARAM extension)! Please make sure you are looking at the correct file - you can download it from ftp://xmm.vilspa.esa.es/pub/ccf/constituents > 2.The size of PARAMS is 6E, not 2E so there are 6 parameters for every > offset angle and energy, not just r_c and alpha. What is their meaning and > how can I use them? in the correct XRT1_XPSF_0005.CCF file, the size is 2E in fact > 3.How about the XRT[23]_XENCIREN_0003.CCF files? I has been told before > that it is these files which are being used to correct for the fraction of > the PSF outside of one's aperture. Still, cor1 and cor2 from those files > are equal to 0 for MOS2 anf PN. This is not directly related to my current > question but I would like to know whether my psf is being corrected if I > use e.g. a 30arcsec aperture within 5arcmin from the center of the FOV and > use appropriate task switches to make the correction. Please, have a look at the CCF release notes available at http://xmm.vilspa.esa.es/ccf/releasenotes/ccfindex.html There you should find a reasonable description of what the different CCFs are and gow they are related. Please, let me know if still some info is needed here - in that case I will need to contact our CCF specialists. Cheers, Matthias. Matthias Ehle XMM-Newton SOC User Support Group
Date: Wed, 14 Aug 2002 15:56:30 +0100 (BST) From: Igor Antokhin <igor@astro.gla.ac.uk> To: Matthias Ehle <xmmhelp@xmm.vilspa.esa.es> Subject: Re: XMM PSF (PR#5635)
On Wed, 14 Aug 2002, Matthias Ehle wrote: > Dear Igor, > > (Michel is on leave for 4 weeks now but I hope that me too, I can help you > a bit with your questions on the PSF...) > > you wrote: > > 1.In your message and in the postscipt description you say that the > > function used is the King function. In XRT1_XPSF_0005.CCF, the extention > > containing the parameters ENERGY, THETA and PARAMS is called > > GLE_GAUSS_PARAM (I could not find the KING_PARAMS extension). Is it still > > ok to use the formula provided in your postscript document (which I > > honestly doubt)? > > hmmm, I just had a look at the XRT1_XPSF_0005.CCF file myself and I see > the extension called KING_PARAMS (and no GLE_GAUSS_PARAM extension)! > Please make sure you are looking at the correct file - you can download it from > ftp://xmm.vilspa.esa.es/pub/ccf/constituents Hi Matthias, I have just downloaded the file you had mentined. Attached is the main screen as seen by the "fv program" (XRT1_XPSF_0005_CCF_1.gif). No "KING_PARAMS" extention. The "GLE_GAUSS_PARAM" extention is shown in XRT1_XPSF_0005_CCF_1.gif from which it is seen that PARAMS have 6 elements. I am lost. Do I look into the wrong place? Thank you for your help. Igor Antokhin
Date: Fri, 16 Aug 2002 11:16:01 +0100 (BST) From: Igor Antokhin <igor@astro.gla.ac.uk> To: Matthias Ehle <xmmhelp@xmm.vilspa.esa.es> Subject: Re: XMM PSF (PR#5635)
This message is in MIME format. The first part should be readable text, while the remaining parts are likely unreadable without MIME-aware tools. Send mail to mime@docserver.cac.washington.edu for more info. ---779972051-132576782-1029492961=:6873 Content-Type: TEXT/PLAIN; charset=US-ASCII Sorry Matthias, I forgot to attach the plots to my last message. Best, Igor > On Wed, 14 Aug 2002, Matthias Ehle wrote: > > > Dear Igor, > > > > (Michel is on leave for 4 weeks now but I hope that me too, I can help you > > a bit with your questions on the PSF...) > > > > you wrote: > > > 1.In your message and in the postscipt description you say that the > > > function used is the King function. In XRT1_XPSF_0005.CCF, the extention > > > containing the parameters ENERGY, THETA and PARAMS is called > > > GLE_GAUSS_PARAM (I could not find the KING_PARAMS extension). Is it still > > > ok to use the formula provided in your postscript document (which I > > > honestly doubt)? > > > > hmmm, I just had a look at the XRT1_XPSF_0005.CCF file myself and I see > > the extension called KING_PARAMS (and no GLE_GAUSS_PARAM extension)! > > Please make sure you are looking at the correct file - you can download it from > > ftp://xmm.vilspa.esa.es/pub/ccf/constituents > > Hi Matthias, > > I have just downloaded the file you had mentined. Attached is the main > screen as seen by the "fv program" (XRT1_XPSF_0005_CCF_1.gif). No > "KING_PARAMS" extention. The "GLE_GAUSS_PARAM" extention is shown in > XRT1_XPSF_0005_CCF_1.gif from which it is seen that PARAMS have 6 > elements. I am lost. Do I look into the wrong place? Thank you for your > help. > > Igor Antokhin ---779972051-132576782-1029492961=:6873 Content-Type: IMAGE/gif; name="XRT1_XPSF_0005_CCF_1.gif" Content-Transfer-Encoding: BASE64 Content-ID: <Pine.LNX.4.44.0208161116010.6873@oberon.astro.gla.ac.uk> Content-Description: Content-Disposition: attachment; filename="XRT1_XPSF_0005_CCF_1.gif" R0lGODdhawITAvcAAAQCBFyCtCxCXKzC1JSSlERijBQiLExKTNTi7HyaxJSq zGRiZCQyRDxSdMzKzISChKyqrExynDQyNGyOtERCRCQmJNTW1HRydLS2tDRK ZCQ6TAwWHMTCxJyenFRWVPz+/IymxJy21DxafFR6rNza3GyKtExqlBwqPOTi 5IyOjBQeJGSKtJyyzGxqbLSytFRypHSWvCQuRNTe7Hx6fLy6vDROdCw6VLzK 3Ozq7GSCtDRGZJyanExmlBwiNFRSVISixJyuzCQ2TDxWdNTS1IyKjLSurDw6 PHSSvCwuLMzW5DRObBQWJKSmpFxeXIyqzERehFR2pCw+VCxGXLTC3JSWlERm jExOTISexJSuzGRmZCQyTMzOzISGhKyurGyOvERGRCwqLHR2dKy+1DRKbMTG xKSipFxaXIymzKS61DxahFx+rNze3ExunBwuPOTm5HyWvBQeLGxubHx+fLy+ vLzO3Ozu7GSGtBwmNDxWfDw+PNTa7BQaJFR2rCw+XAMDB0BACAgIgOTkvhUV EQgICCUYmAIkvoAWEQIICBik2N7s2/8V/78Iv/9+wXeuhwgDBwhACK3/cM8j ZQwAEwgACACgUAAq4QAq/wBAvxgA+N4A3v8A/78Av8FEANqNAAUPAAhAAAgE QOQAEBUAFUwIUALk4YAV/wIIv+AI+N3d3v///7+/vwnccQCujgADBQBACAAI ZgDkAAAVAAAIAAEYAAAkAAAWAAAMAQDrAKh+AN6uAP8DAL9AAAAEAAAAAAgc BMDtAEYVAEAIAACgCAHd5AD/FQC/CHxEAK2NAAcPAAiAAMAAAEYAAEAAAAEI AGs4AALdAAD/AAC/AADiUAC74QAD/6QIAN7kAP8VAL8IAAgYAMAkAEYWAFhE Ad6NAP8PALzmZq0aAAcHAAgDA2tQBALqABNcCAIC5AAAFQAACEhW596dq/8e A79AQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA AAAAAAAAAAAAAAAAAAAAACwAAAAAawITAgAI/gABCBxIsKDBgwgTKlzIsKHD hxAjSpxIsaLFixgzatzIsaPHjyBDirxIBYCMG2hYqFzJsqXLlzBjypxJs6bN mzhz6tzJs6fPn0CDCqWZoMYcMhaSrrFAwsLSpk+dMpUKdWrUq1azVt1KtStW rl+9ahULdmzYs2bTll1Lti1atm/dqpW7ZeANGG9gwDiy94jfv3zzCsZLeLDh wogPK07MeLHjxpAfS45MebLlypgva87MebPnzqA/i4ZcQEcXMiTWuEGxevWa 17Bjy55Nu7bt27hz695dG8Vr32uAC/9NPHjx4caTC/ft2jUK38+js46++nn1 6dapR1fNHLp31cqP/t9eXlcgmjd/vZRYYae9+/YrSngJTX+0/fr47+vPz3+/ //4A/nfZEwIwwQEJbuCAQx0MLqigggxGKOGEFFZo4YUYZqjhhhx26OGHIF7o 2hpNTWVVcA264cZrJDTVooupudEgDsC9aGOLqsk4x0BY+FWCFH3YYIMGGgwZ xJBCahAAYEw26eSTUEYp5ZRUVmnllVhmqeWWXHbp5ZdghulXGkEQ4AJqwSXo Go7BLcjgBwvCWYecdMZp55x31onnnnr2meeffALqZ6CEyhnioYhm6MYFWWTR xKMLQAppFhY0eEELjj6qqaSUNhgHppFuqukClKKAAw0DKaCeHVHoMQAW/meA cEass4LghRZsHOHFrrtO0OuvvAYr7LDB+uqFscYSq+ywyAK77LPANgvttNQe G221y0qbLLbCasstsd5+2y2wH4g7rrXo7lrut+ESuy6zIsSggQ5K4CGCA0NY UEUBIohgbxr4SifwwAQXbPDBCCes8MLS0ejww6ZCrGDEFE9sscQYVyxxgqau hoPHIJuKggVZfGDyySh/4IAPNKxGcsowb8GyyyXDjPLKGDhgQRcDgbCeGjaI EV8JdhC9HtEN3FFACUw3bTLTQzct9dRUl4By00NHXfXWU2cNNddgS73C00dP /YHUZ6fMtM1Obx3f2UenffXaJ9M9d9lRk+00/tx77z231nbzPTfMX1fNd9iF D5024lx7XfbaaDPuOOCRcy2rEwqUMEYXHHDwhB5TsODz5hxscePpqKeu+uqs t+7666gzBfvsrMv+ou20t0iGGR88F9x3KMhBQQe68+778dAJv4PuTfQOfPLD 0+AAEz3nEMAIQZyRQw4A2JGD9wB8P0YPIlhvRwA5fOD9B9tvH8D58L8vf/zw s++9HR+gb/3298/vP/3/a9/+vve9/xkQgOlbn/sK+D31oQ9/+CNgBO+nPgma 7HwFxGAD9ZdA+zVwgh/0YP8WyD/7OXCCJzRhAyWowf6BEIIiVJ8INeg++jkw gwd83wDT90EcAvCH/jpsHwfX154HVhCI8BPgEB+owyMykYHt2Z4IotCBLmBA CGLQ4RR3gAEy5CspYAyjGMdIxjKa8YxoTKMazTiEEoXxi3BMShwtMMc6yvGO dHRjHsX4xTXOgXe4gxEJ5GCEFCTljx8IpO0ImYI2IlKRL2LkaXYwkCsEQA0R 0EICAjAQNaihk2OAQw0CcElSfkANpSRlKVHJyvyxcpWXdGUsrafKV9oSlre8 pSpxyctcwvKUqeRl/gLgSlmqwZjIJCYvVznMYjJTma2EJix3GUtSojJ/w7ym MkuJzWkK85nJNCY1e0lOX+7Smb5MpzdrWU1UlnOZ1PQlMN1JT2t6MgBC/tBC CphQBCUo4JJqwEMQUgCBOeBrCAjNl0LpyNCFOrShEH2oRCN
From: Matthias Ehle <xmmhelp@xmm.vilspa.esa.es> To: igor@astro.gla.ac.uk Subject: Re: XMM PSF (PR#5635) Date: Fri Aug 16 10:46:00 2002
Dear Igor, OK. now I see what's happening: There are much more file extensions in XRT1_XPSF_0005.CCF then you show in your plot! You need to use the scroll bar on the right hand side of fv to display the other extension (in total there are 71). Extension 70 is the one called KING_PARAMS. BTW, the fv window sometimes does not display the complete extension names: GLE_GAUSS_PARAM is indeed an extension called SINGLE_GAUSS_PARAMS (you see the complete name if you display the extension header with fv, the name is given in the EXTNAME keyword) I am glad that we could solve this 'mystery'. Cheers, Matthias,
Date: Fri, 16 Aug 2002 11:57:45 +0100 (BST) From: Igor Antokhin <igor@astro.gla.ac.uk> To: Matthias Ehle <xmmhelp@xmm.vilspa.esa.es> Subject: Re: XMM PSF (PR#5635)
Dear Matthias, thanks for your reply. Stupid of me to not check the rest of the list. It is because after a few first extenstions everything else seemed to be the PSF image examples so I thought they would just go till the end of the list. Stupid of me. Thanks a lot, Igor On Fri, 16 Aug 2002, Matthias Ehle wrote: > Dear Igor, > > OK. now I see what's happening: > > There are much more file extensions in XRT1_XPSF_0005.CCF then you show > in your plot! You need to use the scroll bar on the right hand side of fv > to display the other extension (in total there are 71). > Extension 70 is the one called KING_PARAMS. > > BTW, the fv window sometimes does not display the complete extension names: > GLE_GAUSS_PARAM is indeed an extension called SINGLE_GAUSS_PARAMS > (you see the complete name if you display the extension header with fv, > the name is given in the EXTNAME keyword) > > I am glad that we could solve this 'mystery'. > > Cheers, > Matthias, > -------------------------------------------- Igor Antokhin Astronomy and Astrophysics Group University of Glasgow Kelvin Building, G12 8QQ Glasgow, UK E-mail: igor@astro.gla.ac.uk
Date: Thu, 22 Aug 2002 11:03:59 +0100 (BST) From: Igor Antokhin <igor@astro.gla.ac.uk> To: Matthias Ehle <xmmhelp@xmm.vilspa.esa.es> Subject: Re: XMM PSF (PR#5635)
Dear Matthias, I have successfully used the info on the encircled PSF energy you provided. Thanks. Just one more question: the King's formula from XMM-PS-GM-20 (2.2), p.16, has 5arcmin term in the divisor. It looks like I can use this formula up to 5arcmin off-axis angles only. I recall reading somewhere that XMM PSF is defined within this distance from the center of the FOV. Yet the CCFs (XRT1_XPSF_0005.CCF and others) provide the interpolated r_c and alpha up to ~15arcmin. Does this mean that I could use the King's formula beyond 5arcmin? One of my apertures is 7armin off-axis. Thank you for your help, Igor -------------------------------------------- Igor Antokhin Astronomy and Astrophysics Group University of Glasgow Kelvin Building, G12 8QQ Glasgow, UK E-mail: igor@astro.gla.ac.uk
From: Matthias Ehle <xmmhelp@xmm.vilspa.esa.es> To: igor@astro.gla.ac.uk Subject: Re: XMM PSF (PR#5635) Date: Thu Aug 22 10:47:16 2002
Dear Igor, I have forwarded your question to our XRT experts and will let you know their reply asap. Please note, however, that due to the vacation season, I am not sure if I can get a quick answer... Cheers, Matthias.
From: Matthias Ehle <xmmhelp@xmm.vilspa.esa.es> To: igor@astro.gla.ac.uk Subject: Re: XMM PSF (PR#5635) Date: Tue Aug 27 15:42:54 2002
Dear Igor, please find attached now the expert's reply to your question: The cited formula including the 5arcmin term, is for the EEF. 5arcmin is not the off-axis angle of the source, but a normalization radius. The EEF provides the enrcicled energy fraction enclosed within a radius R. We decided to normalize it at 5 arcmin. In practise, we assume that all the energy of the source is enclosed within a radius of 5 arcmin from the source center (whatever its off-axis position is), or in other words, the EEF evaluates the fraction of the energy within a radius of 5 arcmin. This is an approximation, and a discussion on this point can be found in the report XMM-SOC-CAL-TN-0022 (at the XMM vilspa page). Hence, the formula can be used for any off-axis position of the source. However, let me recall that the data allowed a good calibration only for low off-axis angles (say < about 4 arcmin). For larger off-axis angles only low energies are well calibrated (see the report for the "range of application" of the PSF/EEF model) I hope this info will help you with your further analysis. Cheers, Matthias.
Date: Wed, 28 Aug 2002 11:32:24 +0100 (BST) From: Igor Antokhin <igor@astro.gla.ac.uk> To: Matthias Ehle <xmmhelp@xmm.vilspa.esa.es> Subject: Re: XMM PSF (PR#5635)
Dear Matthias, thanks for your reply which was very useful. Still, while answering my query, it raised other question... Could you please ask your expert to answer just one question regarding my particular case? I have 4 XMM observations of a source (a WR star), in 2 of them (rev 284,285) it is located in the center of the FOV, in two others (rev 115,116) - 7arcmin off-axis. I understand that the 5arcmin term in the EEF formula is not an off-axis angle but a normalization radius. 5armin appeared in my message, because I saw this value as the off-axis limit (where PSF was still defined), in some documentation. I actually just checked the latest release notes and learned that currently the EEF is at least 5% accurate at energies below 4.5keV within the central 10armin. Now, what I see is that the EPIC flux of my source at the 7armin off-axis position is about 15% HIGHER that the flux in the on-axis position, the difference is the same in two energy bands 0.5-2.4keV and 2.4-10.0keV. In both cases, I used the extraction aperture R=1armin and did NOT correct the extracted spectra for the EEF. According to the XRT PSF CCF, the difference between the fluxes due to this correction would be less that 1 per cent. I computed the fluxes in xspec, accounting for the rmf and arf matrixes. I understand that at least for the soft band the difference exceeds possible EEF errors. However, you also mentioned that "...the data allowed a good calibration only for low off-axis angles (say < about 4 arcmin). For larger off-axis angles only low energies are well calibrated (see the report for the "range of application" of the PSF/EEF model)". So the question is: did you mean the above 4.5keV limit? Are there any other calibration uncertainties (maybe apart from the PSF/EEF issue) which could account for the 15% difference of fluxes? I would like to know whether the difference can be considered as real variations or not... Thanks for your help, Igor Antokhin On Tue, 27 Aug 2002, Matthias Ehle wrote: > Dear Igor, > > please find attached now the expert's reply to your question: > > The cited formula including the 5arcmin term, is for the EEF. 5arcmin is not > the off-axis angle of the source, but a normalization radius. ... > Hence, the formula can be used for any off-axis position of the source. > However, let me recall that the data allowed a good calibration only for > low off-axis angles (say < about 4 arcmin). For larger off-axis angles only low > energies are well calibrated (see the report for the "range of > application" of the PSF/EEF model) ... > Cheers, > Matthias.
From: Matthias Ehle <xmmhelp@xmm.vilspa.esa.es> To: igor@astro.gla.ac.uk Subject: Re: XMM PSF (PR#5635) Date: Fri Sep 20 08:43:27 2002
Dear Igor, sorry for my late reply on your PSF followup question. Somehow the system failed informing me that there was a new question raised by you. I have forwarded your report to our experts and will let you know their reply asap. Cheers and sorry again for the 'hickup' in our helpdesk system. Cheers, Matthias.
From: Matthias Ehle <xmmhelp@xmm.vilspa.esa.es> To: igor@astro.gla.ac.uk Subject: Re: XMM PSF (PR#5635) Date: Tue Sep 24 08:33:15 2002
Dear Igor, please, find below an answer from our calibration experts on the subject of flux differences between on- and off-axis source position: you wrote: > Are there any other > calibration uncertainties (maybe apart from the PSF/EEF issue) which > could account for the 15% difference of fluxes? I would like to know > whether the difference can be considered as real variations or not... This is almost certainly due to a calibration error in the vignetting. An uncertainty in the position of the optical-axis used within the SAS is currently giving an error, which we estimate to be about MOS-1 3% MOS-2 10% PN 8% at 7 arcmins off-axis for an energy of 4.5 keV. There is a small energy dependence on the above values with lower energies being affected less. At the moment you should use these numbers as a systematic error on the 7' off-axis fluxes. There is also a 2-4% error introduced by the same effect on the on-axis measurements. We are actively working on this problem and have just taken a set of calibration measurements which should resolve the issue within the next 1-2 months. This may not be enough to solve the full 15% difference that you see and after the vignetting correction becomes available we may need to revisit other sources of potential error to investigate the remaining 5% discrepancy. It is very unlikely that uncertainties in the PSF are causing more than 1% of this particular problem for such a large extraction radius. Best regards, Matthias.