------- Forwarded Message
Date: Wed, 16 Jul 2008 06:50:37 -0400
From: owner-chandra-users@head.cfa.harvard.edu
To: chandra-users-approval@head.cfa.harvard.edu
Subject: BOUNCE chandra-users@head.cfa.harvard.edu: Non-member submission from
["Aneta Siemiginowska" <asiemiginowska@cfa.harvard.edu>]
>From owner-chandra-users@head.cfa.harvard.edu Wed Jul 16 06:50:36 2008
Received: from psmtp.com (exprod7mx202.postini.com [64.18.2.42])
by head.cfa.harvard.edu (w/d) with SMTP id m6GAoY5b011200
for <chandra-users@head.cfa.harvard.edu>; Wed, 16 Jul 2008 06:50:35 -04
00 (EDT)
Received: from source ([209.85.198.249]) by exprod7mx202.postini.com ([64.18.6.
14]) with SMTP;
Wed, 16 Jul 2008 06:50:34 EDT
Received: by rv-out-0708.google.com with SMTP id c5so4949112rvf.24
for <chandra-users@head.cfa.harvard.edu>; Wed, 16 Jul 2008 03:50:33 -07
00 (PDT)
DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed;
d=gmail.com; s=gamma;
h=domainkey-signature:received:received:message-id:date:from:sender
:to:subject:cc:in-reply-to:mime-version:content-type
:content-transfer-encoding:content-disposition:references
:x-google-sender-auth;
bh=JGjtr9W2OmoKNT0QiUxn23TJrhhY9xolojD1hW96l3Y=;
b=xLodl0Fkzefrg8NjdXtcBX/Rn9iTvbGNu/FCzwVyRwVRv0MmJTIyw+2r0jWzA4Y/I2
qyC6VtMSgd4JrQNzk0tgg/cxp8IBa6DPWupTATk+W2N5L/wJd5tec3qF+78qrXUdTI9p
rsJ2BgoDz7c68UX/KoQ076hzwEuimO5kqTeMY=
DomainKey-Signature: a=rsa-sha1; c=nofws;
d=gmail.com; s=gamma;
h=message-id:date:from:sender:to:subject:cc:in-reply-to:mime-version
:content-type:content-transfer-encoding:content-disposition
:references:x-google-sender-auth;
b=Y6+iNc18hDFQuLMM5QRwuNg86voCIVXTUdQsWQrME87H4sz0vggBnfvjJbVLN7vvuM
BvO6QUDP+hVqHx2Te9WIaa8WQVacdkngU2vSt83zGm4R+ezKm8dkaMpiMmc7ECLaTu2E
kf5IkhIOreZLTBBDmx8d84MqDKvOiDF1Euj4w=
Received: by 10.141.137.6 with SMTP id p6mr8163965rvn.279.1216205433817;
Wed, 16 Jul 2008 03:50:33 -0700 (PDT)
Received: by 10.141.195.11 with HTTP; Wed, 16 Jul 2008 03:50:33 -0700 (PDT)
Message-ID: <337019db0807160350v53b1d5a0lb8d7ec28acb56a32@mail.gmail.com>
Date: Wed, 16 Jul 2008 06:50:33 -0400
From: "Aneta Siemiginowska" <asiemiginowska@cfa.harvard.edu>
Sender: anetasie@gmail.com
To: "John ZuHone" <zuhone@flash.uchicago.edu>
Subject: Re: Convolving radial profile with PSF in sherpa
Cc: chandra-users@head.cfa.harvard.edu
In-Reply-To: <53C10725-9769-4DDA-B615-7ACA958E8711@flash.uchicago.edu>
MIME-Version: 1.0
Content-Type: text/plain; charset=ISO-8859-1
Content-Transfer-Encoding: 7bit
Content-Disposition: inline
References: <66321179-C04F-4778-8BB2-B2CF78EC7A7D@flash.uchicago.edu>
<1216178271.15258.109.camel@waimea>
<53C10725-9769-4DDA-B615-7ACA958E8711@flash.uchicago.edu>
X-Google-Sender-Auth: ef64e4a95511a5e9
X-pstn-neptune: 0/0/0.00/0
X-pstn-levels: (S:99.90000/99.90000 CV:99.0000 R:95.9108 P:95.9108 M:97.028
2 C:98.6951 )
X-pstn-settings: 3 (1.0000:1.0000) s cv gt3 gt2 gt1 r p m c
X-pstn-addresses: from <asiemiginowska@cfa.harvard.edu> [1157/51]
On Tue, Jul 15, 2008 at 11:29 PM, John ZuHone <zuhone@flash.uchicago.edu> wrote
:
> Aneta,
>
> Thanks for the helpful comments and code. Just one quick question:
>
> Would I create the 1D PSF by creating a 2D PSF (with, e.g. mkpsf) and then
> using dmextract to create a radial profile?
>
You should use CHART to create a PSF rays, then MARX to make a PSF
event file and then extract radial profile. Check the Thread on CIAO
web pages:
http://cxc.harvard.edu/chart/threads/index.html
Good Luck,
Aneta
> Thanks again,
>
> John Z
>
> On Jul 15, 2008, at 10:17 PM, Aneta Siemiginowska wrote:
>
>>
>> Hello,
>>
>> There are two options depending on whether you are using CIAO 3.4 or
>> CIAO 4 version of Sherpa.
>>
>> In CIAO 3.4 the psf data file can be input into "fpsf1d" instrument
>> model. You should see an example on the command line when you type
>> "ahelp fpsf1d"
>>
>> CIAO4 contains Sherpa Beta, but you should be able to work with radial
>> profile in CIAO4 too. There is no thread example. I include a list of
>> the commands in CIAO4 Python Sherpa Beta below.
>>
>> Hope this helps
>> Aneta Siemiginowska
>>
>>
>>
>>
>> ---------------------------------
>> In CIAO4 Sherpa beta - Python version use the following commands:
>>
>> load_table("rprofile_rmid.fits[cols RMID, SUR_BRI, SUR_BRI_ERR]")
>>
>> # check the data
>> get_data()
>>
>> # plot and modify plots
>>
>> plot_data()
>> log_scale(XY_AXIS)
>> set_plot_xlabel("RMID [ pixel ]")
>> set_plot_ylabel("SUR BRI [ cnts/pixel^{2} ]")
>>
>>
>> # read psf 1d data file
>>
>> load_table("psf","psf_rprofile_rmids.fits[cols RMID, SUR_BRI]")
>> print (get_data("psf"))
>>
>>
>> # plot psf
>>
>> get_data_plot_prefs()["yerrorbars"]=False
>> plot_data("psf")
>> log_scale(XY_AXIS)
>> set_plot_xlabel("RMID [ pixel ]")
>> set_plot_ylabel("SUR BRI [ cnts/pixel^{2} ]")
>> get_data_plot_prefs()["yerrorbars"]=True
>>
>> #check statistics
>> print get_stat_name()
>>
>> # check optimization method
>> print(get_method_name())
>>
>> # change optimization method
>>
>> set_method("neldermead")
>>
>> # set beta model
>>
>> set_model(beta1d.src)
>> src.r0 = 105
>> src.beta = 4
>> src.xpos = 0
>> src.ampl = 0.01
>> src.ampl.max = 10
>> freeze(src.xpos)
>>
>> # check the model
>> print get_model()
>>
>>
>> # set psf data file as a 1d psf model
>>
>> set_psf_model(get_data("psf"), psf1d.pf)
>> pf.radial = 1
>> pf.xsize = 36
>> pf.xcenter = 18
>>
>> # check model expression with the psf
>>
>> print("sherpa> print get_full_model()")
>> print(get_full_model())
>>
>>
>> fit()
>>
>> # plot and modify the plot labels, axes etc.
>>
>> plot_fit()
>>
>> print('sherpa> log_scale XY_AXIS')
>> print('sherpa> set_plot_xlabel "RMID [ pixel ]"')
>> print('sherpa> set_plot_ylabel "SUR BRI [ cnts/pixel^{2} ]"')
>> log_scale(XY_AXIS)
>> set_plot_xlabel("RMID [ pixel ]")
>> set_plot_ylabel("SUR BRI [ cnts/pixel^{2} ]")
>>
>>
>> # You should be careful of the soft parameter limits
>>
>>
>> On Tue, 2008-07-15 at 18:33 -0500, John ZuHone wrote:
>>>
>>> Hello,
>>>
>>> I was wondering how to convolve a 1d radial profile model (such as
>>> beta1d) in sherpa with a file-based PSF and fit it to a radial profile.
>>>
>>> Thanks,
>>>
>>> John Z
>>>
>
>
------- End of Forwarded Message
This archive was generated by hypermail 2b29 : Tue Feb 12 2013 - 01:00:12 EST