Up to top level
AO15   AO16   AO17   AO18   AO19   Backgrounds   Calibration   Conference   Data   Docs   EPICMOS   EPICpn   Feedback   Gallery   Misc   OM   Pending   PhD_Theses   Publications   RGS   RadMonitor   SAS_Hardware   SAS_WS   SASv16.0   SASv16.0_Installation   SASv16.1   SASv16.1_Installation   SASv17.0   SASv17.0_Installation   SASv18.0   SASv18.0_Installation   SciSim   Simulators_other   Suggestions   Trash   Visibility   XMM-bouncing   XMM-news   XRPS   XSA   esas   incoming  

Logged in as guest

Viewing SASv17.0/84825
Full headers

From: p.zhou@uva.nl
Subject: incompatibility between mask-filter and arfgen?
Compose reply
Download message
Move To:
6 replies: 1 2 3 4 5 6
3 followups: 1 2 3

Private message: yes  no

Notes:

Notification:


Date: Tue, 2 Apr 2019 23:35:41 +0200
From: p.zhou@uva.nl
To: xmmhelp@sciops.esa.int
Subject: incompatibility between mask-filter and arfgen?
Full_Name: Ping Zhou
Submission from: (NULL) (193.147.152.102)


Dear colleague,

I came across a problem when extracting the .arf file for a spectrum. Because
the spectral extraction area is not a simply region, I used a mask image instead
of a region, so the expression for filtering is like
"mask(mask_image.fits,0,0,X,Y) && FLAG == 0 && PATTERN <=4.. 

The extracted image, detector image all look fine, and the pha file seems to be
OK. But the .arf file is problematic, showing SPECRESP values (<=3 cm^2) over
100 times less than the expected values (<=800 cm^2). On the other hand, the
backscale values are not correct. If I choose to correct the bad pixels using
.withbadpixcorr=true badpixlocation=${evt2}., >99% of the good pixels are
filtered out for the value of backscale. 

There are no errors reported, but there is an unusual warning saying: 
** arfgen: warning (detmapXBoundsExceeded), The DSS x-axis extent
[-52658.83059,50476.01687] exceeds the detmap x-axis bounds [-18283.5,13996.5].
This may produce an inappropriate ARF
** arfgen: warning (detmapYBoundsExceeded), The DSS y-axis extent
[-50091.07462,53043.77679] exceeds the detmap y-axis bounds [-17527.5,15352.5].
This may produce an inappropriate ARF

I am wondering whether there is an incompatibility between the arfgen and the
mask-filter. How can I solve this problem?

The SAS version: sasversion:- sasversion (sasversion-1.2) 
[xmmsas_20170719_1539-16.1.0]
HESAOFT version: heasoft-6.24/x86_64-apple-darwin16.7.0

The followings are the commands that I used&#65306;

# environmental variables
evt2=`echo ../${instr}*-clean.fits`
      pi_bin=5
      imgbin=82
      pi_max=20479
      filter='(FLAG == 0) && (PATTERN <= 4)'
# spectral extraction

      express="mask(mask_bin_${i}.fits,0,0,X,Y) && ${filter}"
      evselect table=$evt2 withspectrumset=yes spectrumset=spec_${i}.pi \
        energycolumn=PI spectralbinsize=${pi_bin} withspecranges=yes
specchannelmin=0 specchannelmax=${pi_max} \
        writedss=yes expression="${express}"
      evselect table=$evt2 imagebinning=binSize imageset=img_${i}.fits
withimageset=yes \
        xcolumn=X ycolumn=Y ximagebinsize=${imgbin} yimagebinsize=${imgbin}
expression="${express}"
     #generate detector map
      evselect table="${evt2}:EVENTS" destruct=false withfilteredset=true
withimageset=true imageset=detmap.fits xcolumn=DETX ycolumn=DETY
imagebinning=binSize ximagebinsize=120 yimagebinsize=120 expression="${express}"
writedss=true updateexposure=true
       # generate rmf with rmfgen
      rmfgen spectrumset=spec_${i}.pi rmfset=spec_${i}.rmf detmaptype=dataset
detmaparray=detmap.fits
       # generate arf with arfgen
      arfgen spectrumset=spec_${i}.pi arfset=spec_${i}.arf detmaptype=dataset
detmaparray=detmap.fits extendedsource=yes withrmfset=true rmfset=spec_${i}.rmf 
setbackscale=yes withbadpixcorr=true badpixlocation=${evt2}

Best regards,
Ping


Reply 1

Resend
From: Nora Loiseau <xmmhelp@sciops.esa.int>
To: p.zhou@uva.nl
Subject: Re: incompatibility between mask-filter and arfgen? (PR#84825)
Date: Wed Apr  3 12:55:11 2019
Dear Ping Zhou,
> 
> I came across a problem when extracting the .arf file for a spectrum.
Because
> the spectral extraction area is not a simply region, I used a mask image
instead
> of a region, so the expression for filtering is like
> "mask(mask_image.fits,0,0,X,Y) && FLAG == 0 && PATTERN
<=4.. 

I consulted an expert and he advises to extract the spectrum with the command:

"(X,Y) in circle(25000,25000,15000) && mask(mask_image.fits,0,0,X,Y)
&& FLAG ==
0 && PATTERN <=4.. 

we imagine that you are interested in an extended source, then arfgen should be
run with extendesource=yes  

Please contact us again if you encounter more problems,

Best regards,

Nora


----
Dr. Nora Loiseau
XMM-Newton User Support Group


Followup 1

Compose reply
Download message
From: Ping Zhou <p.zhou@uva.nl>
Subject: Re: incompatibility between mask-filter and arfgen? (PR#84825)
Date: Thu, 4 Apr 2019 11:56:06 +0200
To: Nora Loiseau <xmmhelp@sciops.esa.int>
--Apple-Mail=_7E000771-5D5F-4DE5-A1B5-CE2EDED4C94E
Content-Transfer-Encoding: quoted-printable
Content-Type: text/plain;
	charset=utf-8

Dear Nora,

Thanks a lot for the reply and the suggestion.=20

Yes, I=E2=80=99m interested in an extended source and used =
extendedsource=3Dyes in arfgen.

Following your suggestion, I tested a few regions at the beginning of =
the expression.
It seems that the SPECRESP values in arf file strongly depend on the =
region selection
(e.g, circle(25000,25000,15000), circle(25000,25000,20000)) rather than =
the mask
file. The larger region I choose, the smaller SPECRESP values I obtain.

For comparison, I made two arf files using
1. expression=3D"(X,Y) in circle(25000,25000,15000) && mask
=E2=80=A6=E2=80=
=9D
and 2. expression=3D"box(25261.721,23168.657,15460,15900,49.998953, X,Y) =
 &&=20
mask =E2=80=A6=E2=80=9D respectively. The 2nd region is selected to =
cover the whole pn camera.

I expected to see the arf files look similar given the same mask-filter =
(~1 arcmin^2 area),=20
but they turned out to be different. Please check the attached image =
(left=E2=80=94expression 1,=20
right=E2=80=94 expression 2).

Do I have to use the region (X,Y) in circle(25000,25000,15000)? Although =
it covers the=20
mask-region that I just tested, it doesn=E2=80=99t cover some other =
regions at the CCD edges.=20
So there could be a limitation of using it.

Best regards,
Ping





> On 3 Apr 2019, at 12:55 PM, Nora Loiseau <xmmhelp@sciops.esa.int> =
wrote:
>=20
> Dear Ping Zhou,
>>=20
>> I came across a problem when extracting the .arf file for a spectrum. =
Because
>> the spectral extraction area is not a simply region, I used a mask =
image
> instead
>> of a region, so the expression for filtering is like
>> "mask(mask_image.fits,0,0,X,Y) && FLAG =3D=3D 0 &&
PATTERN <=3D4=C2=94.=
=20
>=20
> I consulted an expert and he advises to extract the spectrum with the =
command:
>=20
> "(X,Y) in circle(25000,25000,15000) &&
mask(mask_image.fits,0,0,X,Y) =
&& FLAG =3D=3D
> 0 && PATTERN <=3D4..=20
>=20
> we imagine that you are interested in an extended source, then arfgen =
should be
> run with extendesource=3Dyes =20
>=20
> Please contact us again if you encounter more problems,
>=20
> Best regards,
>=20
> Nora
>=20
>=20
> ----
> Dr. Nora Loiseau
> XMM-Newton User Support Group
>=20
> This message is intended only for the recipient(s) named above. It may =
contain proprietary information and/or
> protected content. Any unauthorised disclosure, use, retention or =
dissemination is prohibited. If you have received
> this e-mail in error, please notify the sender immediately. ESA =
applies appropriate organisational measures to protect
> personal data, in case of data privacy queries, please contact the ESA =
Data Protection Officer (dpo@esa.int).
>=20
>=20


--Apple-Mail=_7E000771-5D5F-4DE5-A1B5-CE2EDED4C94E
Content-Type: multipart/related;
	type="text/html";
	boundary="Apple-Mail=_96CD3E02-27FE-4CF5-8737-452EC8424031"


--Apple-Mail=_96CD3E02-27FE-4CF5-8737-452EC8424031
Content-Transfer-Encoding: quoted-printable
Content-Type: text/html;
	charset=utf-8

<html><head><meta http-equiv=3D"Content-Type"
content=3D"text/html =
charset=3Dutf-8"></head><body style=3D"word-wrap: break-word; =
-webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" =
class=3D"">Dear Nora,<div class=3D""><br
class=3D""></div><div =
class=3D"">Thanks a lot for the reply and the =
suggestion.&nbsp;</div><div class=3D""><br
class=3D""></div><div =
class=3D"">Yes, I=E2=80=99m interested in an extended source and used =
extendedsource=3Dyes in arfgen.</div><div class=3D""><br =
class=3D""></div><div class=3D"">Following your suggestion, I
tested a =
few regions at the beginning of the expression.</div><div
class=3D"">It =
seems that the SPECRESP values in arf file strongly depend on the region =
selection</div><div class=3D"">(e.g, circle(25000,25000,15000), =
circle(25000,25000,20000)) rather than the mask</div><div
class=3D"">file.=
 The larger region I choose, the smaller SPECRESP values I =
obtain.</div><div class=3D""><br
class=3D""></div><div class=3D"">For =
comparison, I made two arf files using</div><div class=3D"">1. =
expression=3D"(X,Y) in circle(25000,25000,15000) &amp;&amp; mask =
=E2=80=A6=E2=80=9D</div><div class=3D"">and 2. =
expression=3D"box(25261.721,23168.657,15460,15900,49.998953, X,Y)&nbsp; =
&amp;&amp;&nbsp;</div><div class=3D"">mask
=E2=80=A6=E2=80=9D =
respectively. The 2nd region is selected to cover the whole pn =
camera.</div><div class=3D""><br
class=3D

Message of length 205378 truncated


Reply 2

Resend
From: Nora Loiseau <xmmhelp@sciops.esa.int>
To: p.zhou@uva.nl
Subject: Re: incompatibility between mask-filter and arfgen? (PR#84825)
Date: Thu Apr  4 18:21:24 2019
Dear Ping,

I am consulting a SAS expert and I will come back to you asap.

Best regards,

Nora


> 
> Yes, I=E2=80=99m interested in an extended source and used =
> extendedsource=3Dyes in arfgen.
> 
> Following your suggestion, I tested a few regions at the beginning of =
> the expression.
> It seems that the SPECRESP values in arf file strongly depend on the =
> region selection
> (e.g, circle(25000,25000,15000), circle(25000,25000,20000)) rather than =
> the mask
> file. The larger region I choose, the smaller SPECRESP values I obtain.
> 
> For comparison, I made two arf files using
> 1. expression=3D"(X,Y) in circle(25000,25000,15000) && mask
=E2=80=A6=E2=80=
> =9D
> and 2. expression=3D"box(25261.721,23168.657,15460,15900,49.998953, X,Y) =
>  &&=20
> mask =E2=80=A6=E2=80=9D respectively. The 2nd region is selected to =
> cover the whole pn camera.
> 
> I expected to see the arf files look similar given the same mask-filter =
----
Dr. Nora Loiseau
XMM-Newton User Support Group


Reply 3

Resend
From: Nora Loiseau <xmmhelp@sciops.esa.int>
To: p.zhou@uva.nl
Subject: Re: incompatibility between mask-filter and arfgen? (PR#84825)
Date: Fri Apr  5 12:49:51 2019
Dear Ping,

Could you please send us the mask_image.fits file and the command that you use
when
you call arfgen ?

Regards,

Nora

> Thanks a lot for the reply and the suggestion.=20
> 
> Yes, I=E2=80=99m interested in an extended source and used =
> extendedsource=3Dyes in arfgen.
> 
> Following your suggestion, I tested a few regions at the beginning of =
> the expression.
> It seems that the SPECRESP values in arf file strongly depend on the =
> region selection
> (e.g, circle(25000,25000,15000), circle(25000,25000,20000)) rather than =
> the mask
> file. The larger region I choose, the smaller SPECRESP values I obtain.
> 
> For comparison, I made two arf files using
> 1. expression=3D"(X,Y) in circle(25000,25000,15000) && mask
=E2=80=A6=E2=80=
> =9D
> and 2. expression=3D"box(25261.721,23168.657,15460,15900,49.998953, X,Y) =
>  &&=20
> mask =E2=80=A6=E2=80=9D respectively. The 2nd region is selected to =
> cover the whole pn camera.
> 
> I expected to see the arf files look similar given the same mask-filter =
----
Dr. Nora Loiseau
XMM-Newton User Support Group


Followup 2

Compose reply
Download message
From: Ping Zhou <p.zhou@uva.nl>
Subject: Re: incompatibility between mask-filter and arfgen? (PR#84825)
Date: Fri, 5 Apr 2019 12:57:44 +0200
To: Nora Loiseau <xmmhelp@sciops.esa.int>
--Apple-Mail=_D23F9AB5-A27D-43B0-B7AB-38B7E1D3EBAA
Content-Transfer-Encoding: quoted-printable
Content-Type: text/plain;
	charset=us-ascii

Dear Nora,

Attached is the mask image.=20

This is the arfgen command:

arfgen spectrumset=3Dspec_${i}.pi arfset=3Dspec_${i}.arf =
detmaptype=3Ddataset detmaparray=3Ddetmap.fits extendedsource=3Dyes =
withrmfset=3Dtrue rmfset=3Dspec_${i}.rmf setbackscale=3Dyes =
withbadpixcorr=3Dtrue badpixlocation=3D${evt2}=20

For the pha, detmap and rmf extraction, I used:
     express=3D"box(25261.721,23168.657,15460,15900,49.998953, X,Y) && =
mask(mask_bin_${i}.fits,0,0,X,Y) && (FLAG =3D=3D 0) && (PATTERN
<=3D 4)"
    =20
     evselect table=3D$evt2 withspectrumset=3Dyes =
spectrumset=3Dspec_${i}.pi \
        energycolumn=3DPI spectralbinsize=3D5 withspecranges=3Dyes =
specchannelmin=3D0 specchannelmax=3D20479 \
        writedss=3Dyes expression=3D"${express}"

     evselect table=3D"${evt2}:EVENTS" destruct=3Dfalse =
withfilteredset=3Dtrue withimageset=3Dtrue imageset=3Ddetmap.fits =
xcolumn=3DDETX ycolumn=3DDETY imagebinning=3DbinSize ximagebinsize=3D120 =
yimagebinsize=3D120 expression=3D"${express}" writedss=3Dtrue =
updateexposure=3Dtrue
    =20
     rmfgen spectrumset=3Dspec_${i}.pi rmfset=3Dspec_${i}.rmf =
detmaptype=3Ddataset detmaparray=3Ddetmap.fits


Thanks,
Ping


> On 5 Apr 2019, at 12:49 PM, Nora Loiseau <xmmhelp@sciops.esa.int> =
wrote:
>=20
> Dear Ping,
>=20
> Could you please send us the mask_image.fits file and the command that =
you use
> when
> you call arfgen ?
>=20
> Regards,
>=20
> Nora
>=20
>> Thanks a lot for the reply and the suggestion.=3D20
>>=20
>> Yes, I=3DE2=3D80=3D99m interested in an extended source and used =3D
>> extendedsource=3D3Dyes in arfgen.
>>=20
>> Following your suggestion, I tested a few regions at the beginning of =
=3D
>> the expression.
>> It seems that the SPECRESP values in arf file strongly depend on the =
=3D
>> region selection
>> (e.g, circle(25000,25000,15000), circle(25000,25000,20000)) rather =
than =3D
>> the mask
>> file. The larger region I choose, the smaller SPECRESP values I =
obtain.
>>=20
>> For comparison, I made two arf files using
>> 1. expression=3D3D"(X,Y) in circle(25000,25000,15000) && mask =
=3DE2=3D80=3DA6=3DE2=3D80=3D
>> =3D9D
>> and 2. expression=3D3D"box(25261.721,23168.657,15460,15900,49.998953, =
X,Y) =3D
>> &&=3D20
>> mask =3DE2=3D80=3DA6=3DE2=3D80=3D9D respectively. The 2nd region is =
selected to =3D
>> cover the whole pn camera.
>>=20
>> I expected to see the arf files look similar given the same =
mask-filter =3D
> ----
> Dr. Nora Loiseau
> XMM-Newton User Support Group
>=20
> This message is intended only for the recipient(s) named above. It may =
contain proprietary information and/or
> protected content. Any unauthorised disclosure, use, retention or =
dissemination is prohibited. If you have received
> this e-mail in error, please notify the sender immediately. ESA =
applies appropriate organisational measures to protect
> personal data, in case of data privacy queries, please contact the ESA =
Data Protection Officer (dpo@esa.int).
>=20
>=20


--Apple-Mail=_D23F9AB5-A27D-43B0-B7AB-38B7E1D3EBAA
Content-Type: multipart/mixed;
	boundary="Apple-Mail=_A21C4978-020B-473D-9AD3-92E45392F202"


--Apple-Mail=_A21C4978-020B-473D-9AD3-92E45392F202
Content-Transfer-Encoding: quoted-printable
Content-Type: text/html;
	charset=us-ascii

<html><head><meta http-equiv=3D"Content-Type"
content=3D"text/html =
charset=3Dus-ascii"></head><body style=3D"word-wrap: break-word; =
-webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" =
class=3D"">Dear Nora,<div class=3D""><br
class=3D""></div><div =
class=3D"">Attached is the mask image.&nbsp;</div><div
class=3D""><br =
class=3D""></div><div class=3D"">This is the arfgen
command:</div><div =
class=3D""><br class=3D""></div><div class=3D"">arfgen =
spectrumset=3Dspec_${i}.pi arfset=3Dspec_${i}.arf detmaptype=3Ddataset =
detmaparray=3Ddetmap.fits extendedsource=3Dyes withrmfset=3Dtrue =
rmfset=3Dspec_${i}.rmf setbackscale=3Dyes withbadpixcorr=3Dtrue =
badpixlocation=3D${evt2}&nbsp;</div><div class=3D""><br =
class=3D""></div><div class=3D""><font color=3D"#919191"
class=3D"">For =
the pha, detmap and rmf extraction, I used:</font></div><div =
class=3D""><div class=3D""><font color=3D"#919191"
class=3D"">&nbsp; =
&nbsp; &nbsp;express=3D"box(25261.721,23168.657,15460,15900,49.998953, =
X,Y) &amp;&amp; mask(mask_bin_${i}.fits,0,0,X,Y) &amp;&amp;
(FLAG =3D=3D =
0) &amp;&amp; (PATTERN &lt;=3D 4)"&l

Message of length 4414833 truncated


Reply 4

Resend
From: Nora Loiseau <xmmhelp@sciops.esa.int>
To: p.zhou@uva.nl
Subject: Re: incompatibility between mask-filter and arfgen? (PR#84825)
Date: Thu Apr 11 16:22:39 2019
Dear Ping,
 I contacted the SAS developer and he informs me that the problem is not so
simple.
We will come back to you as soon as he has an answer.

Best regards,

Nora

>
> 
> Dear Nora,
> 
> Attached is the mask image.=20
> 
----
Dr. Nora Loiseau
XMM-Newton User Support Group


Reply 5

Resend
From: Nora Loiseau <xmmhelp@sciops.esa.int>
To: p.zhou@uva.nl
Subject: Re: incompatibility between mask-filter and arfgen? (PR#84825)
Date: Thu Apr 25 20:03:30 2019
Dear Ping,

I am very sorry for the delay in coming back to you. 
The SAS expert says that the way of using a MASK to include data 
is not well supported within the SAS. If you want to use the MASK 
to include the region then it is better to use a polygon to select 
that region. Then you will be able to extract the spectrum from 
the polygon and run arfgen with the extendedsource=yes 
option.
If what you want it to exclude that region then please get back in touch with
us.

Best regards,

Nora

> 
> Attached is the mask image.=20
> 
> This is the arfgen command:
> 
> arfgen spectrumset=3Dspec_${i}.pi arfset=3Dspec_${i}.arf =
> detmaptype=3Ddataset detmaparray=3Ddetmap.fits extendedsource=3Dyes =
> withrmfset=3Dtrue rmfset=3Dspec_${i}.rmf setbackscale=3Dyes =
> withbadpixcorr=3Dtrue badpixlocation=3D${evt2}=20
> 
> For the pha, detmap and rmf extraction, I used:
>      express=3D"box(25261.721,23168.657,15460,15900,49.998953, X,Y)
&& =
> mask(mask_bin_${i}.fits,0,0,X,Y) && (FLAG =3D=3D 0) &&
(PATTERN <=3D 4)"
>     =20
>      evselect table=3D$evt2 withspectrumset=3Dyes =
> spectrumset=3Dspec_${i}.pi \
>         energycolumn=3DPI spectralbinsize=3D5 withspecranges=3Dyes =
> specchannelmin=3D0 specchannelmax=3D20479 \
>         writedss=3Dyes expression=3D"${express}"
> 
>      evselect table=3D"${evt2}:EVENTS" destruct=3Dfalse =
> withfilteredset=3Dtrue withimageset=3Dtrue imageset=3Ddetmap.fits =
> xcolumn=3DDETX ycolumn=3DDETY imagebinning=3DbinSize ximagebinsize=3D120 =
> yimagebinsize=3D120 expression=3D"${express}" writedss=3Dtrue =
> updateexposure=3Dtrue
>     =20
>      rmfgen spectrumset=3Dspec_${i}.pi rmfset=3Dspec_${i}.rmf =
> detmaptype=3Ddataset detmaparray=3Ddetmap.fits
> 
> 
> Thanks,
> Ping
> 
> 
>> On 5 Apr 2019, at 12:49 PM, Nora Loiseau <xmmhelp@sciops.esa.int>
=
> wrote:
>>=20
>> Dear Ping,
>>=20
>> Could you please send us the mask_image.fits file and the command that
=
> you use
>> when
>> you call arfgen ?
>>=20
>> Regards,
>>=20
>> Nora
>>=20
>>> Thanks a lot for the reply and the suggestion.=3D20
>>>=20
>>> Yes, I=3DE2=3D80=3D99m interested in an extended source and used
=3D
>>> extendedsource=3D3Dyes in arfgen.
>>>=20
>>> Following your suggestion, I tested a few regions at the beginning
of =
> =3D
>>> the expression.
>>> It seems that the SPECRESP values in arf file strongly depend on
the =
> =3D
>>> region selection
>>> (e.g, circle(25000,25000,15000), circle(25000,25000,20000)) rather
=
> than =3D
>>> the mask
>>> file. The larger region I choose, the smaller SPECRESP values I =
> obtain.
>>>=20
>>> For comparison, I made two arf files using
>>> 1. expression=3D3D"(X,Y) in circle(25000,25000,15000) &&
mask =
> =3DE2=3D80=3DA6=3DE2=3D80=3D
>>> =3D9D
>>> and 2. expression=3D3D"box(25261.721,23168.657,15460,15900,49.998953,
=
> X,Y) =3D
>>> &&=3D20
>>> mask =3DE2=3D80=3DA6=3DE2=3D80=3D9D respectively. The 2nd region is
=
> selected to =3D
>>> cover the whole pn camera.
>>>=20
>>> I expected to see the arf files look similar given the same =
> mask-filter =3D
>> ----
>> Dr. Nora Loiseau
>> XMM-Newton User Support Group
>>=20
>> This message is intended only for the recipient(s) named above. It may
=
> contain proprietary information and/or
>> protected content. Any unauthorised disclosure, use, retention or =
> dissemination is prohibited. If you have received
>> this e-mail in error, please notify the sender immediately. ESA =
> applies appropriate organisational measures to protect
>> personal data, in case of data privacy queries, please contact the ESA
=
> Data Protection Officer (dpo@esa.int).
>>=20
>>=20
> 
> 
> --Apple-Mail=_D23F9AB5-A27D-43B0-B7AB-38B7E1D3EBAA
> Content-Type: multipart/mixed;
> 	boundary="Apple-Mail=_A21C4978-020B-473D-9AD3-92E45392F202"
> 
> 
> --Apple-Mail=_A21C4978-020B-473D-9AD3-92E45392F202
> Content-Transfer-Encoding: quoted-printable
> Content-Type: text/html;
> 	charset=us-ascii
> 
> <html><head><meta http-equiv=3D"Content-Type"
content=3D"text/html =
> charset=3Dus-ascii"></head><body style=3D"word-wrap:
break-word; =
> -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" =
> class=3D"">Dear Nora,<div class=3D""><br
class=3D""></div><div =
> class=3D"">Attached is the mask image..</div><div
c

Message of length 4759330 truncated


Reply 6

Resend
From: Nora Loiseau <xmmhelp@sciops.esa.int>
To: p.zhou@uva.nl
Subject: Re: incompatibility between mask-filter and arfgen? (PR#84825)
Date: Thu Apr 25 20:10:06 2019
Dear Ping,

I am very sorry for the delay in coming back to you.
The SAS expert says that the way of using a MASK to include data
is not well supported within the SAS. If you want to use the MASK
to include the region then it is better to use a polygon to select
that region. Then you will be able to extract the spectrum from
the polygon and run arfgen with the extendedsource=yes
option.
If what you want it to exclude that region then please get back in touch with
us.

Best regards,

Nora


> 
> Dear Nora,
> 
> Attached is the mask image.=20
> 
> This is the arfgen command:
> 
> arfgen spectrumset=3Dspec_${i}.pi arfset=3Dspec_${i}.arf =
> detmaptype=3Ddataset detmaparray=3Ddetmap.fits extendedsource=3Dyes =
> withrmfset=3Dtrue rmfset=3Dspec_${i}.rmf setbackscale=3Dyes =
> withbadpixcorr=3Dtrue badpixlocation=3D${evt2}=20
> 
> For the pha, detmap and rmf extraction, I used:
>      express=3D"box(25261.721,23168.657,15460,15900,49.998953, X,Y)
&& =
> mask(mask_bin_${i}.fits,0,0,X,Y) && (FLAG =3D=3D 0) &&
(PATTERN <=3D 4)"
----
Dr. Nora Loiseau
XMM-Newton User Support Group


Followup 3

Compose reply
Download message
From: Ping Zhou <p.zhou@uva.nl>
Subject: Re: incompatibility between mask-filter and arfgen? (PR#84825)
Date: Fri, 26 Apr 2019 11:03:29 +0200
To: Nora Loiseau <xmmhelp@sciops.esa.int>
Dear Nora,

Thank you very much for helping me with this question. And thank your =
SAS colleague as well.=20

Then I know I need to switch to polygon.

Best,
Ping
=20

> On 25 Apr 2019, at 8:10 PM, Nora Loiseau <xmmhelp@sciops.esa.int> =
wrote:
>=20
> Dear Ping,
>=20
> I am very sorry for the delay in coming back to you.
> The SAS expert says that the way of using a MASK to include data
> is not well supported within the SAS. If you want to use the MASK
> to include the region then it is better to use a polygon to select
> that region. Then you will be able to extract the spectrum from
> the polygon and run arfgen with the extendedsource=3Dyes
> option.
> If what you want it to exclude that region then please get back in =
touch with
> us.
>=20
> Best regards,
>=20
> Nora
>=20
>=20
>>=20
>> Dear Nora,
>>=20
>> Attached is the mask image.=3D20
>>=20
>> This is the arfgen command:
>>=20
>> arfgen spectrumset=3D3Dspec_${i}.pi arfset=3D3Dspec_${i}.arf =3D
>> detmaptype=3D3Ddataset detmaparray=3D3Ddetmap.fits =
extendedsource=3D3Dyes =3D
>> withrmfset=3D3Dtrue rmfset=3D3Dspec_${i}.rmf setbackscale=3D3Dyes =3D
>> withbadpixcorr=3D3Dtrue badpixlocation=3D3D${evt2}=3D20
>>=20
>> For the pha, detmap and rmf extraction, I used:
>>     express=3D3D"box(25261.721,23168.657,15460,15900,49.998953, X,Y) =
&& =3D
>> mask(mask_bin_${i}.fits,0,0,X,Y) && (FLAG =3D3D=3D3D 0)
&& (PATTERN =
<=3D3D 4)"
> ----
> Dr. Nora Loiseau
> XMM-Newton User Support Group
>=20
> This message is intended only for the recipient(s) named above. It may =
contain proprietary information and/or
> protected content. Any unauthorised disclosure, use, retention or =
dissemination is prohibited. If you have received
> this e-mail in error, please notify the sender immediately. ESA =
applies appropriate organisational measures to protect
> personal data, in case of data privacy queries, please contact the ESA =
Data Protection Officer (dpo@esa.int).
>=20
>=20


Up to top level
AO15   AO16   AO17   AO18   AO19   Backgrounds   Calibration   Conference   Data   Docs   EPICMOS   EPICpn   Feedback   Gallery   Misc   OM   Pending   PhD_Theses   Publications   RGS   RadMonitor   SAS_Hardware   SAS_WS   SASv16.0   SASv16.0_Installation   SASv16.1   SASv16.1_Installation   SASv17.0   SASv17.0_Installation   SASv18.0   SASv18.0_Installation   SciSim   Simulators_other   Suggestions   Trash   Visibility   XMM-bouncing   XMM-news   XRPS   XSA   esas   incoming  

Logged in as guest


Please make your (short) question the subject of your request!


Web interface using JitterBug ... back to the XMM home page