integral/ridge
1
0
Fork 0
generated from erosita/uds
Code Issues Pull Requests Packages Projects Releases Wiki Activity

IKI

Browse Source
This commit is contained in:
Roman Krivonos 2024-08-27 14:02:54 +03:00
parent 5292290a35
commit 6781f091f9
5 changed files with 19 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

BIN
GalProf/detcnts.E01.ALL.resid.gal.fits
View File

Binary file not shown.

BIN
GalProf/detcnts.E01.ALL.resid.galprof.fits
View File

Binary file not shown.

8
scripts/03_grxe_galprof.py
View File

@ -58,14 +58,14 @@ maxiters=10
with open(ignored_rev_file, 'rb') as fp: with open(ignored_rev_file, 'rb') as fp:
ignored_rev = pickle.load(fp) ignored_rev = pickle.load(fp)
print(ignored_rev) #print(ignored_rev)
print("{} orbits ignored".format(len(ignored_rev))) print("{} orbits ignored".format(len(ignored_rev)))
ign=ignored_rev.tolist() ign=ignored_rev.tolist()
glon, step = np.linspace(-lon_max, lon_max, num=lon_nbin, endpoint=False,retstep=True) glon, step = np.linspace(-lon_max, lon_max, num=lon_nbin, endpoint=False,retstep=True)
print(glon,step) #print(glon,step)
# #
@ -99,6 +99,8 @@ lat_arr=[]
for i in range(lon_nbin): for i in range(lon_nbin):
dkey="{:04d}".format(i) dkey="{:04d}".format(i)
#if(i != 5):
# continue
df0 = df.query('LON > {} & LON <= {} & LAT > {} & LAT < {} & REV != @ign'.format(glon[i],glon[i]+step,-lat_max,lat_max)) df0 = df.query('LON > {} & LON <= {} & LAT > {} & LAT < {} & REV != @ign'.format(glon[i],glon[i]+step,-lat_max,lat_max))
@ -124,7 +126,7 @@ for i in range(lon_nbin):
lon_arr.append(row['LON']) lon_arr.append(row['LON'])
lat_arr.append(row['LAT']) lat_arr.append(row['LAT'])
#sys.exit()
if not os.path.exists(profdir): if not os.path.exists(profdir):

2
scripts/03_grxe_galprof.sh
View File

@ -1,4 +1,4 @@
./03_grxe_galprof.py E01 71 ./03_grxe_galprof.py E01 51
./03_grxe_galprof.py E13 17 ./03_grxe_galprof.py E13 17
./03_grxe_galprof.py E14 21 ./03_grxe_galprof.py E14 21

18
scripts/03_grxe_galprof_plot.py
View File

@ -4,7 +4,7 @@ __author__ = "Roman Krivonos"
__copyright__ = "Space Research Institute (IKI)" __copyright__ = "Space Research Institute (IKI)"
from astropy.table import Table, Column from astropy.table import Table, Column
from matplotlib import ticker
import matplotlib.pyplot as plt import matplotlib.pyplot as plt
import sys import sys
@ -51,30 +51,38 @@ ax1.ticklabel_format(style='sci', axis='y', scilimits=(0,0))
ax2.ticklabel_format(style='sci', axis='y', scilimits=(0,0)) ax2.ticklabel_format(style='sci', axis='y', scilimits=(0,0))
ax3.ticklabel_format(style='sci', axis='y', scilimits=(0,0)) ax3.ticklabel_format(style='sci', axis='y', scilimits=(0,0))
ax1.xaxis.set_minor_locator(ticker.MultipleLocator(10))
ax2.xaxis.set_minor_locator(ticker.MultipleLocator(10))
ax3.xaxis.set_minor_locator(ticker.MultipleLocator(10))
ax1.yaxis.set_minor_locator(ticker.MultipleLocator(0.2))
ax2.yaxis.set_minor_locator(ticker.MultipleLocator(0.2))
ax3.yaxis.set_minor_locator(ticker.MultipleLocator(0.2))
lon=(df1['LON1']+df1['LON2'])/2 lon=(df1['LON1']+df1['LON2'])/2
ax1.errorbar(lon, df1['GRXE_SIM_FLUX']/scale, ax1.errorbar(lon, df1['GRXE_SIM_FLUX']/scale,
yerr=df1['GRXE_SIM_ERROR']/scale, yerr=df1['GRXE_SIM_ERROR']/scale,
xerr=(df1['LON2']-df1['LON1'])/2, xerr=(df1['LON2']-df1['LON1'])/2,
fmt='o' ) fmt='o', color='black')
ax1.grid(visible=True) ax1.grid(visible=True)
lon=(df2['LON1']+df2['LON2'])/2 lon=(df2['LON1']+df2['LON2'])/2
ax2.errorbar(lon, df2['GRXE_SIM_FLUX']/scale, ax2.errorbar(lon, df2['GRXE_SIM_FLUX']/scale,
yerr=df2['GRXE_SIM_ERROR']/scale, yerr=df2['GRXE_SIM_ERROR']/scale,
xerr=(df2['LON2']-df2['LON1'])/2,fmt='o' ) xerr=(df2['LON2']-df2['LON1'])/2,fmt='o' , color='black')
ax2.grid(visible=True) ax2.grid(visible=True)
lon=(df3['LON1']+df3['LON2'])/2 lon=(df3['LON1']+df3['LON2'])/2
ax3.errorbar(lon, df3['GRXE_SIM_FLUX']/scale, ax3.errorbar(lon, df3['GRXE_SIM_FLUX']/scale,
yerr=df3['GRXE_SIM_ERROR']/scale, yerr=df3['GRXE_SIM_ERROR']/scale,
xerr=(df3['LON2']-df3['LON1'])/2, fmt='o' ) xerr=(df3['LON2']-df3['LON1'])/2, fmt='o' , color='black')
ax3.grid(visible=True) ax3.grid(visible=True)
#x = np.arange(-10, 10, 0.001) #x = np.arange(-10, 10, 0.001)
#plot normal distribution with mean 0 and standard deviation 1 #plot normal distribution with mean 0 and standard deviation 1
#plt.plot(x, norm.pdf(x, 0, 1), color='red', linewidth=2) #plt.plot(x, norm.pdf(x, 0, 1), color='red', linewidth=2)
plt.xlabel('Galactic Longitude',fontsize=14, fontweight='normal') plt.xlabel('Galactic Longitude, deg.',fontsize=14, fontweight='normal')
ax2.set_ylabel('GRXE flux, x100 mCrab',fontsize=14, fontweight='normal') ax2.set_ylabel('GRXE flux, x100 mCrab',fontsize=14, fontweight='normal')
#plt.xscale('linear') #plt.xscale('linear')