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 Backgrounds/85524
Full headers

From: jjung11@alaska.edu
Subject: evqpb spectrum doesn't match with sciece spectrum
Compose reply
Download message
Move To:
1 replies: 1
1 followups: 1

Private message: yes  no

Notes:

Notification:


Date: Thu, 8 Aug 2019 02:25:52 +0200
From: jjung11@alaska.edu
To: xmmhelp@sciops.esa.int
Subject: evqpb spectrum doesn't match with sciece spectrum
Full_Name: Jaewoong Jung
Submission from: (NULL) (193.147.152.102)


Dear XMM helpdesk team,

I'm sorry, I think I'm bothering you with this problem again and again, but I
still can't understand it, so your help would be appreciated.

I'm doing the background study, so I removed astrophysical sources from the
observation data.
And made a spectrum from large region, circle (-50,-180,16000), with both
science exposure and FWC exposure generated by evqpb. (I also removed same
region from the FWC files)
And I loaded two spectrum (science and FWC) into Xspec. I expected they should
be almost same, as I removed astrophysical sources, like in the figure 3 of (
https://www.cosmos.esa.int/web/xmm-newton/sas-thread-background  - here they
didn't removed sources, but in the high energy band, like 7.5-9.2 keV, science
and FWC should be exactly same.) But my result was not so.


The observation I used was 0150610101.

The SAS commands I used were

cifbuild
odfingest
emproc
epproc

#Proton flaring filtering
evselect table=0623_0150610101_EMOS1_0150610101_ImagingEvts.ds withrateset=Y
rateset=rateMOS1.fits maketimecolumn=Y timebinsize=100 makeratecolumn=Y
expression='#XMMEA_EM && (PI>10000) && (PATTERN==0)'
dsplot table=rateMOS1.fits x=TIME y=RATE
tabgtigen table=rateMOS1.fits expression='RATE<=.35' gtiset=MOS1gti.fits
evselect table=0623_0150610101_EMOS1_0150610101_ImagingEvts.ds withfilteredset=Y
filteredset=MOS1clean.fits destruct=Y keepfilteroutput=T expression="#XMMEA_EM
&& gti(MOS1gti.fits,TIME) && (PI>150)"
evselect table=MOS1clean.fits withrateset=Y rateset=rateMOS1clean.fits
maketimecolumn=Y timebinsize=100 makeratecolumn=Y expression='#XMMEA_EM
&&
(PI>10000) && (PATTERN==0)'
dsplot table=rateMOS1clean.fits x=TIME y=RATE

#Sources finding
evselect table=MOS1clean.fits imagebinning=binSize imageset=MOS1_image_full.fits
withimageset=y xcolumn=X ycolumn=Y ximagebinsize=40 yimagebinsize=40
expression='#XMMEA_EM && (PI in [200:12000]) && (PATTERN in
[0:12]) &&
(FLAG==0)'
edetect_chain imagesets=MOS1_image_full.fits eventsets=MOS1clean.fits
attitudeset=0623_0150610101_AttHk.ds pimin=200 pimax=12000 ecf=2.0
eboxl_list=MOS1_eboxlist_l.fits eboxm_list=MOS1_eboxlist_m.fits
esp_nsplinenodes=16 eml_list=MOS1_emllist.fits esen_mlmin=15
srcdisplay boxlistset=MOS1_emllist.fits imageset=MOS1_image_full.fits
sourceradius=.01

#Sources removal
region -w 0 -V 0 eventset=MOS1clean.fits operationstyle=global
radiusstyle=userfixed fixedradius=35 srclisttab=MOS1_emllist.fits
expression='ID_INST==0' outunit=detxy bkgregionset=MOS1reg.fits
evselect table=MOS1cleanf.fits imagebinning=binSize imageset=MOS1_bkg.fits
withimageset=y xcolumn=DETX ycolumn=DETY ximagebinsize=64 yimagebinsize=64
expression='#XMMEA_EM && (PI in [200:10000]) && (PATTERN in
[0:12])'
imgdisplay imagefile=MOS1_bkg.fits withimagefile=y

#spectrum generation
evselect table=MOS1cleanf.fits withspectrumset=y spectrumset=MOS1spec.fits
energycolumn=PI spectralbinsize=5 withspecranges=y specchannelmin=0
specchannelmax=11999 expression='#XMMEA_EM && (PATTERN<=12)
&& ((DETX,DETY) in
circle (-50,-180,16000))' withfilteredset=y filteredset=MOS1cleanfs.fits
backscale spectrumset=MOS1spec.fits badpixlocation=MOS1cleanfs.fits
rmfgen spectrumset=MOS1spec.fits rmfset=MOS1.rmf
arfgen spectrumset=MOS1spec.fits arfset=MOS1.arf withrmfset=y rmfset=MOS1.rmf
badpixlocation=MOS1cleanfs.fits detmaptype=flat extendedsource=y

#FWC file and spectrum generation
evqpb table=0623_0150610101_EMOS1_0150610101_ImagingEvts.ds exposurefactor=2.0
attfile=0623_0150610101_AttHk.ds outset=evqpb_MOS1.fits
evselect table=evqpb_MOS1.fits withfilteredset=y filteredset=MOS1_qpb_clean.fits
destruct=y keepfilteroutput=T expression="gti(MOS1gti.fits,TIME)"
evselect table=MOS1_qpb_cleanf.fits imagebinning=binSize
imageset=MOS1_bkg_fwc.fits withimageset=y xcolumn=DETX ycolumn=DETY
ximagebinsize=64 yimagebinsize=64 expression='#XMMEA_EM && (PI in
[200:10000])
&& (PATTERN in [0:12])'
imgdisplay imagefile=MOS1_bkg_fwc.fits withimagefile=y withcommands=t
evselect table=MOS1_qpb_cleanf.fits withspectrumset=y
spectrumset=MOS1fwc_spec.fits energycolumn=PI spectralbinsize=5 withspecranges=y
specchannelmin=0 specchannelmax=11999 expression='#XMMEA_EM &&
(PATTERN<=12) &&
((DETX,DETY) in circle (-50,-180,16000))' withfilteredset=y
filteredset=MOS1_qpb_cleanfs.fits zcolumn=EWEIGHT
backscale spectrumset=MOS1fwc_spec.fits badpixlocation=MOS1_qpb_cleanfs.fits
fparkey none MOS1spec.fits BACKFILE add=yes
fparkey none MOS1spec.fits ANCRFILE add=yes
fparkey none MOS1fwc_spec.fits BACKFILE add=yes
fparkey none MOS1fwc_spec.fits ANCRFILE add=yes

Then I used Xspec and typed

     1  data MOS1spec.fits
     2  resp MOS1.rmf
     3  arf MOS1.arf
     4  back MOS1fwc_spec.fits
     5  setplot energy
     6  setplot command log x off
     7  setplot background
     8  setplot command color

Message of length 5537 truncated

Followup 1

Compose reply
Download message
From: Jaewoong Jung <jjung11@alaska.edu>
Date: Wed, 7 Aug 2019 16:43:41 -0800
Subject: Re: evqpb spectrum doesn't match with sciece spectrum (PR#85524)
To: Simone Migliari <xmmhelp@sciops.esa.int>
--00000000000016f279058f905a63
Content-Type: multipart/alternative; boundary="00000000000016f274058f905a61"

--00000000000016f274058f905a61
Content-Type: text/plain; charset="UTF-8"

This is the plot between 7.5-9.2 keV for science-FWC (black) spectrum.
(Don't mind red spectrum. It's raw FWC.)
It's not quite zero, and actually 10% of the total science spectrum as I
wrote in the request.

On Wed, Aug 7, 2019 at 4:26 PM <xmmhelp@sciops.esa.int> wrote:

> European  Space  Agency    -    XMM-Newton  Science  Operations  Centre
> -----------------------------------------------------------------------
>
> Dear colleague,
>
> Thank you for submitting your request. We shall deal with it as quickly
> as possible.
>
> Please  note that both your incoming mail and  the answer to it will be
> visible via  the  XMM-Newton Helpdesk  web interface, unless you marked
> your request as PRIVATE.
>
> You can enquire about its status by following this URL:
>
>
> http://xmm.esac.esa.int/xmmhelp?findid=85524
>
> Please note that  your feedback is extremely  useful  for us to improve
> the service  provided to the  Scientific  Community. When  you consider
> that  your  problem  or  question  has  been answered, please visit the
> following web address and send us your comments:
>
>
>
> https://xmm-tools.cosmos.esa.int/external/xmm_user_support/user_feedback/ufeedback.php?service=helpdesk&name_ref=85524
>
>
> Thank you very much for your attention and time.
>
>
> 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).
>
>
>

--00000000000016f274058f905a61
Content-Type: text/html; charset="UTF-8"
Content-Transfer-Encoding: quoted-printable

<div dir=3D"ltr">This is the plot between 7.5-9.2 keV for science-FWC
(blac=
k) spectrum. (Don&#39;t mind red spectrum. It&#39;s raw
FWC.)<div>It&#39;s =
not quite zero, and actually 10% of the total science spectrum as I wrote i=
n the request.=C2=A0=C2=A0</div></div><br><div
class=3D"gmail_quote"><div d=
ir=3D"ltr" class=3D"gmail_attr">On Wed, Aug 7, 2019 at 4:26 PM &lt;<a
href=
=3D"mailto:xmmhelp@sciops.esa.int">xmmhelp@sciops.esa.int</a>&gt;
wrote:<br=
></div><blockquote class=3D"gmail_quote" style=3D"margin:0px 0px 0px
0.8ex;=
border-left:1px solid rgb(204,204,204);padding-left:1ex">European=C2=A0 Spa=
ce=C2=A0 Agency=C2=A0 =C2=A0 -=C2=A0 =C2=A0 XMM-Newton=C2=A0 Science=C2=A0 =
Operations=C2=A0 Centre<br>
----------------------------------------------------------------------- <br=
>
<br>
Dear colleague,<br>
<br>
Thank you for submitting your request. We shall deal with it as
quickly<br>
as possible. <br>
<br>
Please=C2=A0 note that both your incoming mail and=C2=A0 the answer to it w=
ill be <br>
visible via=C2=A0 the=C2=A0 XMM-Newton Helpdesk=C2=A0 web interface, unless=
 you marked <br>
your request as PRIVATE.<br>
<br>
You can enquire about its status by following this URL:<br>
<br>
<br>
<a href=3D"http://xmm.esac.esa.int/xmmhelp?findid=3D85524" rel=3D"noreferre=
r" target=3D"_blank">http://xmm.esac.esa.int/xmmhelp?findid=3D85524</a>=C2=
=A0 <br>
<br>
Please note that=C2=A0 your feedback is extremely=C2=A0 useful=C2=A0 for us=
 to improve<br>
the service=C2=A0 provided to the=C2=A0 Scientific=C2=A0 Community. When=C2=
=A0 you consider <br>
that=C2=A0 your=C2=A0 problem=C2=A0 or=C2=A0 question=C2=A0 has=C2=A0 been =
answered, please visit the <br>
following web address and send us your comments:<br>
<br>
<br>
<a href=3D"https://xmm-tools.cosmos.esa.int/external/xmm_user_support/user_=
feedback/ufeedback.php?service=3Dhelpdesk&amp;name_ref=3D85524" rel=3D"nore=
ferrer" target=3D"_blank">https://xmm-tools.cosmos.esa.int/external/xmm_use=
r_support/user_feedback/ufeedback.php?service=3Dhelpdesk&amp;name_ref=3D855=
24</a>=C2=A0 <br>
<br>
Thank you very much for your attention and time. <br>
<br>
<br>
This message is intended only for the recipient(s) named above. It may cont=
ain proprietary information and/or<br>
protected content. Any unauthorised disclosure, use, retention or dissemina=
tion is prohibited. If you have received<br>
this e-mail in error, please notify the sender immediately. ESA applies app=
ropriate organisational measures to protect<br

Message of length 72468 truncated


Reply 1

Resend
From: Jacobo Ebrero <xmmhelp@sciops.esa.int>
To: jjung11@alaska.edu
Subject: Re: evqpb spectrum doesn't match with sciece spectrum (%2%555%2)
Date: Thu Aug  8 10:09:06 2019
Dear Jaewoong,

Your analysis seems to be right and the result looks reasonable too. The QPB is
typically more or less constant and it is expected (although it may vary from
one observation to another) to be of the order of 10%. Estimations vary but you
can find some here:

https://xmm-tools.cosmos.esa.int/external/xmm_user_support/documentation/uhb/epicintbkgd.html

or here:

https://www.cosmos.esa.int/web/xmm-newton/epic-background-components

Best regards,
Jacobo

----------------------
Dr. Jacobo Ebrero
XMM-Newton SOC
Community Support and Scientific Planning Team

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