Category of the project |
II |
---|
Using 3 years of Fermi-LAT data we have analyzed the region of HESS J1857, a PWN detected at TeV energy by HESS, MAGIC and Veritas.
The corresponding pulsar PSR J1856+0245, found in Arecibo PALFA survey, remains undetected in Fermi-data.
An excess of emission was observed (TS~38) in Fermi data consistent with a point source at the position of HESS J1857+026. Its spectrum is well fitted by a hard power law (gamma~1.5) and the SED obtained using Fermi data is consistent with those of HESS and MAGIC.
We are aiming for a A&A letter with a first draft before the end of the next week.
=Task completed
= To do
Required Tasks |
Status |
SG.Coord Approval |
Internal Referee/s Approval |
Pub. Board Approval |
---|---|---|---|---|
First presentation to the group |
|
|
|
|
Paper category |
|
REQUIRED |
|
|
External Authors |
|
REQUIRED |
|
REQUIRED |
LAT Internal Technical review* |
|
REQUIRED |
|
|
Final Draft to internal referee |
|
REQUIRED |
|
|
Revised draft and sign-up (2 weeks) |
|
REQUIRED |
REQUIRED |
|
Walkthrough/Runtrough |
|
REQUIRED |
REQUIRED |
|
Draft with author list/ack. |
|
|
|
|
Request to submit |
REQUIRED |
REQUIRED |
REQUIRED |
|
Request to submit on ArXiv |
REQUIRED |
|
REQUIRED |
|
Request to resubmit after Journal referee comments |
REQUIRED |
REQUIRED |
REQUIRED |
*Presentation of the LAT data analysis to the group.
LAT Contact Authors:
Name |
contribution to this project |
---|---|
R. Rousseau |
Analysis of Fermi data |
M.-H. Grondin |
Search for pulsed emission |
A. Van Etten |
Broad-Band modelling |
M. Lemoine-Goumard |
Overall coordination |
Other LAT Contributors:
Name |
contribution to this project |
---|---|
B. Stappers |
Search for pulsed emission (ephemerid not covered by the MoU) (Affiliated member) |
A. Lyne |
Search for pulsed emission (ephemerid not covered by the MoU) (PTC) |
C. Espinoza |
Search for pulsed emission (ephemerid not covered by the MoU) (Affiliated member) |
T. Reposeur |
Search for pulsed emission |
D.-A. Smith |
Search for pulsed emission |
S. Johnston |
Search for pulsed emission (Affiliated member) |
J. Hessels |
X-Ray upper limit (PSC) |
V. Kaspi |
X-Ray upper limit (PTC) |
F. Camillo |
X-Ray upper limit (PTC) |
External Authors (Requirements: SG Coord Approval,Pub-board Approval)
Name |
contribution to this project |
|||
---|---|---|---|---|
|
|
|
||
|
|
|||
|
|
Agenda of the day
- Standards for High Level Analisys
- LAT Statistics Board
- Interactive LAT Source Catalog
- Analysis Users forum
- LAT Analysis Page (Workbook)
General Information (this is just an example, please update it)
Source list (if applicable)
NAME |
1/2FGL NAME |
RA |
DEC |
Notes |
---|---|---|---|---|
|
|
|
|
|
Data Set |
Pass 7 |
Event class |
Source (evclsmin/max = 2) |
Energy range |
100 MeV - 300 GeV |
Time interval |
UTC_Start -UTC_stop (MET 239557417 - ) |
ROI size |
10° |
Zenith angle (applied also to gtltcube?) |
< 100° |
Time cuts filter |
DATA_QUAL==1 && LAT_CONFIG==1 && ABS(ROCK_ANGLE)<52 |
Science Tools version |
v9r21p0 |
IRFs P7_V6_Source |
|
Diffuse emission |
ring_2years_P76_v0.fits, isotrop_2year_P76_source_v0.txt and limb_2year_P76_source_v0_smooth.txt |
Optimizer and tolerance |
Minuit (1e-3 ABS) |
Catalog/s |
Preliminary 2-year (gll_psc24month_v2.fits) |
Spatial and Spectral Analysis
The whole analysis is summarized here https://confluence.slac.stanford.edu/pages/viewpage.action?pageId=100515585.
HESS data has shown an emission located at RA=284.3°, DEC=2.68°. It remained unidentified untill the detection of PSR J1856+026 in Arecibo PALFA survey (Hessels et al 2008).
Here we tested the source for extension using pointlike then gtlike and found no significant extension.
Fig. 1 Residual TS map obtained between 10GeV and 300GeV. Here HESS J1857 is not included in the model.
The green contours are those obtained using HESS data (Aharonian et al.,2008).
The SED shown a hard spectrum consistent with those of HESS and MAGIC. The best fit between 300MeV and 300GeV yields to the following parameters :
Int Flux(100MeV-100GeV) = (5.79 ± 0.75 ± 3.11) X 10^{-9} MeV/cm^2/s
Gamma = 1.52 ± 0.16 ± 0.55
With a TS of 38.7.
Systematic Errors Analysis
The systematics were estimated using the bracketting IRFS, fitting the source galactic background level at +/- 6% of the value of the best fit and the systematics on the shape of the source were computed using a template consistent with the extension derived using HESS data.
Broad-Band modelling
Other Analysis
HESS 1857 is powered by a young and energetic pulsar, with a gamma-ray efficiency of 3.7% consistent with those observed using HESS data (3.1%) and consistent with the other values observed for the PWNe detected by Fermi.
An upper limit on the DC emission pulsar yield to a luminosity < 7.47 X 10^34 erg/s.
version |
---|
01/11/2011 : Second draft (v1)
20/09/2011 : First draft (v0) (one column , two column )
Figures
Title |
N. |
Plot/image |
Details |
---|---|---|---|
|
|
|
|
Tables
Title |
N. |
link |
Details |
---|---|---|---|
|
|
|
|
Internal referee report 1 (A. Caliandro)
It is important to describe the shape and the size of the ROI.
In section 3.2 you write that sources for the xml model are extracted in a region of 10deg around HESS J1857. I guess you used the same ROI for the data selection. If it is the case (that you used 10deg ROI for the data selection) I suggest to use for the xml model a larger extraction radius (for example 15deg). In this way you avoid possible issues at the borders of your region, and the diffuse emission will be better constrained. Plausibly, this will lead to a TS value for your source slightly higher.
To evaluate the systematics due to the Galactic diffuse emission, you chnge the normalization by +-6%. This is fine, but 6% refers to the uncertainty of the model we used with P6V3. The new model you are using to analyze P7V6 data most probably and hopefully is better than the previous one.
So the systematics you evaluated are over-estimated. YOu ´can have a feeling on how large is the over-estimation comparing the systematics due to the diffuse with those due to IRFS.
Sincerely, I do not know whether a study equivalent to those performed for W51C and W49 is already done for the new Galactic diffuse model.
Since this paper is a CatII, I do not think that this kind of work should be included in it. But in the text should be clarify what I wrote above.
Anyway, this is a general problem that we should discuss within the collaboration, involving the diffuse group.
#------------------------------------------------------------------------------------------------#
I will comment section by section here below.
-- Introduction.
I found quite confusing the third paragraph of the introduction (L36)
In particular the period:
In such sources, TeV radiation can be ex-
plained by Inverse Compton (IC) scattering of an old population of leptons produced earlier in the pulsar’s life on ambient photon ?lds (CMB, IR, ...) or by ? decay from the interaction of accelerated hadrons with nuclei of the interstellar medium.
What do you mean exactly with 'an old population of leptons'?
Then, in the first paragraph you say that 'The dissipation of the rotational energy of a pulsar leads to the creation of a relativistic wind made of electron/positron pairs.'.
So why you expect pi0 decay in PWN? Maybe a more detailed explanation, and some reference will improve.
-- Timing analysis
This is a request mainly for the radio people. It is possible to evaluate the significance of the last frequency derivative in the ephemeris?
In other words, how much the fourth derivative improve the fit respect to a fit with just three derivatives?
-- Spatial and spectral analysis
For the description of the diffuse models you could also reference the 2FGL paper.
The description of pointlike is quite poor. Further also gtlike is used in binned mode here (I guess). I would avoid to enter in the details of the pointlike description in tis paper. I suggest to just state that pointlike is optimized for the search of the extension of a source, while it is not the case for gtlike. I would cite also the paper in preparation about search for spatially extended sources, in addition to Kerr 2011.
-- Shape and position of the source
Since you finisc the previous chapter talking about the refitting of W44, to avoid confusion I would rename this one as ' Shape and position of the HESS J1857 counterpart'
Since you cite Fig1 (b) as first, I would exchange the positions of the maps, so that the TSmap for E>10GeV is in panel (a).
The TSmap for E>10GeV is not a residual map, as you write in the caption. Otherwise we should see no source, but just statistical fluctuations.
The ellipse errors on the localization of the source should be specified. you state that the localization is consistent with the position found in HESS. What about the consistency with the position of the radio pulsar?
-- Spectral analysis
I already commented on the systematic uncertainties.
The analysis concerning the evaluation of the pulsar DC emission (described in the last paragraph) is a little bit tricky. I ask here the same question than above. Is the radio position of the pulsar compatible with the pointlike localization? If so, you can not disentangle the pulsar DC emission from the PWN emission.
Authors response to referee 1
*"It is important to describe the shape and the size of the ROI.*
*In section 3.2 you write that sources for the xml model are extracted in a region of 10deg around HESS J1857. I guess you used the same ROI for the data selection. If it is the case (that you used 10deg ROI for the data selection) I suggest to use for the xml model a larger extraction radius (for example 15deg). In this way you avoid possible issues at the borders of your region, and the diffuse emission will be better constrained. Plausibly, this will lead to a TS value for your source slightly higher."*
This is a typo, we used a 15° extraction region to build the model. We fixed it in the text.
"
*I found quite confusing the third paragraph of the introduction (L36)*
*In particular the period:* *In such sources, TeV radiation can be explained by Inverse Compton (IC) scattering of an old population of leptons produced earlier in the pulsar’s life on ambient photon ?lds (CMB, IR, ...) or by ?decay from the interaction of accelerated hadrons with nuclei of the interstellar medium.* *What do you mean exactly with 'an old population of leptons'?* *Then, in the first paragraph you say that 'The dissipation of the rotational energy of a pulsar leads to the creation of a relativistic wind made of electron/positron pairs.'. ** **So why you expect pi0 decay in PWN? Maybe a more detailed explanation, and some reference will improve."
We removed these part of the text in the new version to be able to fit in the 4-page limit of A&A letters.
*-- Timing analysis*
*This is a request mainly for the radio people. It is possible to evaluate the significance of the last frequency derivative in the ephemeris?* *In other words, how much the fourth derivative improve the fit respect to a fit with just three derivatives?*
By including a fourth derivative the RMS of the timing residuals decreases in a bout 30%.
It goes from 1551.106 us to 1081.297 us, which is equivalent to a decrease from 19 to 13 milliperiods.
The reduced chi-square of the fit also decreases significantly, changing from 22 to 5.
Even though they have no physical meaning, the use of higher order derivatives helps to flat the residuals,
allowing a more accurate folding of the gamma-ray photons on the rotational period.
*The ellipse errors on the localization of the source should be specified. you state that the localization is consistent with the position found in HESS. What about the consistency with the position of the radio pulsar?*
We added the average error on the position.
*Is the radio position of the pulsar compatible with the pointlike localization? If so, you can not disentangle the pulsar DC emission from the PWN emission.*
The radio position of the pulsar is consistent with those of the GeV source. Nevertheless, it is clear that the hard spectrum that we obtain at very high energy does not ressemble the one of a pulsar. Then, to prevent any contamination of the upper limit on the pulsar by the PWN, we added the PWN best fit parameters in our source model. This way, we can be sure that the derived upper limit is robust (as was previously done for HESS J1825).
internal referee report 2 (name)
Authors response to referee 2
Journal referee report
Authors response to the referee
- SCIENCE GROUP PAGES
- LAT Speaker Bureau
- Templates for Presentations
- OTHER LINKS