Skip to content
GitLab
Projects
Groups
Snippets
/
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
Andrew Strong
galplot
Commits
46993239
Commit
46993239
authored
Jan 14, 2020
by
Andrew Strong
Browse files
Sources.cc: 4FGL but spectra, profiles seem problem: source spectra parameters interpretation?
parent
58371361
Changes
1
Hide whitespace changes
Inline
Side-by-side
source/Sources.cc
View file @
46993239
...
...
@@ -692,24 +692,29 @@ void Sources::readFromFile(const std::string & fileName){
table
.
column
(
"SPECTRAL_INDEX"
).
read
(
index
,
1
,
nRows
);
//Set the pivot point to 100 MeV
pivot
=
100
;
}
catch
(
CCfits
::
Table
::
NoSuchColumn
){
}
catch
(
CCfits
::
Table
::
NoSuchColumn
){
try
{
std
::
cout
<<
" trying
xNickName"
<<
std
::
endl
;
// NickName is in 3FGL, so avoid this test, want next block
std
::
cout
<<
" trying
Source_Name for 4FGL"
<<
std
::
endl
;
// 4FGL AWS20200113
//Read in source names
table
.
column
(
"
xNick
Name"
).
read
(
names
,
1
,
nRows
);
table
.
column
(
"
Source_
Name"
).
read
(
names
,
1
,
nRows
);
//And then longitude and latitude
table
.
column
(
"GLON"
).
read
(
l
,
1
,
nRows
);
table
.
column
(
"GLAT"
).
read
(
b
,
1
,
nRows
);
//flux > 100 MeV, prefactor and spectral index next
table
.
column
(
"Flux100"
).
read
(
flux100
,
1
,
nRows
);
//AWS20140109 NB not in final 2FGL. Set later from energy ranges.
//
table.column("Flux100").read(flux100,1,nRows); //AWS20140109 NB not in final 2FGL. Set later from energy ranges.
table
.
column
(
"Flux1000"
).
read
(
flux1000
,
1
,
nRows
);
//AWS20110806 2FGL
std
::
cout
<<
"read flux1000: first entry="
<<
flux1000
[
1
]
<<
std
::
endl
;
table
.
column
(
"Flux_Density"
).
read
(
prefactor
,
1
,
nRows
);
table
.
column
(
"Spectral_Index"
).
read
(
index
,
1
,
nRows
);
table
.
column
(
"Pivot_Energy"
).
read
(
pivot
,
1
,
nRows
);
table
.
column
(
"PL_Flux_Density"
).
read
(
prefactor
,
1
,
nRows
);
//4FGL
table
.
column
(
"PL_Index"
).
read
(
index
,
1
,
nRows
);
//4FGL
table
.
column
(
"Pivot_Energy"
).
read
(
pivot
,
1
,
nRows
);
//4FGL
// table.column("Flux_Density").read(prefactor,1,nRows); //3FGL
// table.column("Spectral_Index").read(index,1,nRows);
// table.column("Pivot_Energy").read(pivot,1,nRows);
}
catch
(
CCfits
::
Table
::
NoSuchColumn
){
try
{
std
::
cout
<<
" trying Source_Name"
<<
std
::
endl
;
// final official FSSC 2FGL
std
::
cout
<<
" trying Source_Name
for 3FGL
"
<<
std
::
endl
;
// final official FSSC 2
,3,4
FGL
//Read in source names
table
.
column
(
"Source_Name"
).
read
(
names
,
1
,
nRows
);
//And then longitude and latitude
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment