3
_
@ s d Z dd ZdS )a<
******************************************************
*** Program for the convolution of line profile ***
*** with the instrumental response ***
******************************************************
This code is part of the LIFELINE program.
Copyright (C) 2020-2021 University of Liege (Belgium)
Enmanuelle Mossoux (STAR Institute)
This program is free software: you can redistribute it and/or modify it under the terms of the GNU General Public
License as published by the Free Software Foundation, either version 3 of the License, or any later version.
This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied
warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
details.
You should have received a copy of the GNU General Public License along with this program.
If not, see .
Run: from observed_line_profile_main import convolve
p=convolve(direct_LP, file_LP, direct_rmf_arf, RMF, ARF, distance)
# Parameters
# ==========
direct_LP - Directory where the line profile file is located
file_LP - file containing the theoretical line profile
direct_rmf_arf - Directory where the response matrix file (RMF) and ancillary response file (ARF) are located
RMF - Response matrix file
ARF - Ancillary response file
Distance - Distance to the observed binary [kpc].
For Athena, you can use the files computed for a XIFU mirror module radius Rmax=1190mm, a 2.3mm rib spacing and a on-axis case given in the main directory of the code.
RMF: XIFU_CC_BASELINECONF_THICKFILTER_2018_10_10.rmf.
ARF: XIFU_CC_BASELINECONF_THICKFILTER_2018_10_10.arf.
# Output
# ======
p - numpy array containing the convolved emission in each bin
# Versions
# ========
v1 - 03/03/2020
c B C s ddl m} dd l}dd l} dd l}
dd l} ddlm} dd lj }
dd l
}dd l }t | d | d}g }g }g }xp|j
}|sP |jd}|d dkrt|d }|d d
kr|jt|d |jt|d qr|j| qrW |j |j||j|j| }|j|j| }d | }|| d | d
| }d | }d| d| j |d | d d }|dd |d d8 d }|j|d |d g||d9 |d: gf}|j|}|j |d | } | d j} d} d} d| krd} d} | | jjd} | | jjd}!| | jjd}"| | jjd}#t|#j dkrP|#d d d }#| | jjd}$t|$j dkr|$d d d }$| | jjd}%| | jjd}&| | jjd}'|'|& d }(g })t | t ||d kst |!t ||d k r t d t d |
j xt!t"t|D ]}*||* }+|j#| |+|d k d },|j#|!|+|d kd }-|j$t"t|%t"|-d |,d; d f}.x(t!t"|-d |,d< d D ]}/|#t"|/|,d= }0|$t"|/|,d> }1|j%|1}2t&|1t'rxd}3t|1dkrB|j(|"t"|/|,d? d d }4|j#|2|4k d }3|3d d }3|0|3 }0|1|3 }1|"t"|/|,d@ |3|3|1 d }5n |"t"|/|,dA d d }5d d |5D |.|0|0|1 |/f< qW |j)|.dd }6||(j*d!|6j*d!d"d#}7|7|}8|)j|8 q W |d$ d% }9|d&kr~|j |d | }:|:d jjd} |:d jjd}!|:d jjd'};| |! d }(||(j*d!|;j*d!d"d#}7|7|}9g }xpt!t"t|D ]\}?|?t|d kr|>j/d.t.|d d/ t.t"|= d0 |>j/d1 P |>j/||? qHW x@t!t"t|D ],}@|>j/t.|dC d t.|<|@ d0 qW |>j |<|j-|< }Change the RMF file and run 'python observed_line_profile.py'.c S s g | ]}|qS r
).0xr
r
C/local/network/htdocs/Lifeline/codes3/observed_line_profile_main.py
s z convolve..)axisf8cubic)kindg g ? ZSPECRESPg DAz&The number of photons observed during zks is: z_convolved.datawz(# The number of photons observed during zks is
z3# tangential velocity (km/s) | spectrum (photon/s)
g?)r
z Energy (keV)z Normalized convolved emissionz_convolved.pdfzThe output file is r r r r r r r r r r r )5scipy.interpolater numpyossysmath
constantesr matplotlib.pyplotpyplotpandaspyfitsopenreadlinesplitfloatappendclosearrayargsortsortpi concatenatemedianheaderdatafieldlenshapeminmaxprintexitrangeintwherezeroscumsum
isinstancelistargmaxmeanastypeabssumnansumstrwriteplotylimxlabelylabelsavefig)BZ direct_LPZfile_LPdirect_rmf_arfRMFARFdistance expo_timer npr r r! r pltpdr& ZfLPvtangZemiss_thZphraseslinevecenergy
wavelengthwaveZ energy_thbinsZ
bin_lengthZrmfZ
rmf_headerZext_1Zext_2Z energy_loZ energy_himatrixZf_chanZn_chanchannelZe_minZe_maxZ energy_lo_hiZ matrix_totalibinZ
energy_binind1ind2Z
matrix_newZiindZ f_chan_hereZ n_chan_herecum_sumZ
ind_cumsumindZ matrix_hereZ mean_matrixZfinterpZ
matrix_interpZarf_binZarfZspecrespZLP_convZ
nbr_photonfprofileZiphraseiy_limitr
r
r convolve, s
& ,
4*&
" & $
.
$ *
,
ri N)__doc__ri r
r
r
r * s