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
csst-sims
csst_mci_sim
Commits
912f171d
Commit
912f171d
authored
Oct 25, 2024
by
Yan Zhaojun
Browse files
update
parent
7a2eca2f
Pipeline
#7087
failed with stage
in 0 seconds
Changes
5
Pipelines
1
Hide whitespace changes
Inline
Side-by-side
csst_mci_sim/earthshine.py
View file @
912f171d
...
@@ -57,7 +57,7 @@ def ill2flux(E, path):
...
@@ -57,7 +57,7 @@ def ill2flux(E, path):
# use template from sky_bkg (background_spec_hst.dat)
# use template from sky_bkg (background_spec_hst.dat)
filename
=
path
+
'MCI_inputData/refs/background_spec_hst.dat'
filename
=
path
+
'MCI_inputData/refs/background_spec_hst.dat'
cat_spec
=
pd
.
read_csv
(
filename
,
sep
=
'\s+'
,
header
=
None
,
comment
=
'#'
)
cat_spec
=
pd
.
read_csv
(
filename
,
sep
=
'
\
\
s+'
,
header
=
None
,
comment
=
'#'
)
wave0
=
cat_spec
[
0
].
values
# A
wave0
=
cat_spec
[
0
].
values
# A
spec0
=
cat_spec
[
2
].
values
# erg/s/cm^2/A/arcsec^2
spec0
=
cat_spec
[
2
].
values
# erg/s/cm^2/A/arcsec^2
...
...
csst_mci_sim/straylight.py
View file @
912f171d
...
@@ -57,7 +57,7 @@ def ill2flux(E, path):
...
@@ -57,7 +57,7 @@ def ill2flux(E, path):
# use template from sky_bkg (background_spec_hst.dat)
# use template from sky_bkg (background_spec_hst.dat)
filename
=
path
+
'MCI_inputData/refs/background_spec_hst.dat'
filename
=
path
+
'MCI_inputData/refs/background_spec_hst.dat'
cat_spec
=
pd
.
read_csv
(
filename
,
sep
=
'\s+'
,
header
=
None
,
comment
=
'#'
)
cat_spec
=
pd
.
read_csv
(
filename
,
sep
=
'
\
\
s+'
,
header
=
None
,
comment
=
'#'
)
wave0
=
cat_spec
[
0
].
values
# A
wave0
=
cat_spec
[
0
].
values
# A
spec0
=
cat_spec
[
2
].
values
# erg/s/cm^2/A/arcsec^2
spec0
=
cat_spec
[
2
].
values
# erg/s/cm^2/A/arcsec^2
...
...
csst_mci_sim/support/cosmicrays.py
View file @
912f171d
...
@@ -48,7 +48,7 @@ class cosmicrays():
...
@@ -48,7 +48,7 @@ class cosmicrays():
self
.
log
=
log
self
.
log
=
log
#image and size
#
image and size
self
.
image
=
image
.
copy
()
self
.
image
=
image
.
copy
()
self
.
ysize
,
self
.
xsize
=
self
.
image
.
shape
self
.
ysize
,
self
.
xsize
=
self
.
image
.
shape
...
@@ -180,9 +180,7 @@ class cosmicrays():
...
@@ -180,9 +180,7 @@ class cosmicrays():
return
crImage
return
crImage
##################################################
##################################################
############################################
############################################
############################################################################
#####################################
########################################
def
_drawEventsToCoveringFactor
(
self
,
coveringFraction
=
3.0
,
limit
=
1000
,
verbose
=
False
):
def
_drawEventsToCoveringFactor
(
self
,
coveringFraction
=
3.0
,
limit
=
1000
,
verbose
=
False
):
"""
"""
Generate cosmic ray events up to a covering fraction and include it to a cosmic ray map (self.cosmicrayMap).
Generate cosmic ray events up to a covering fraction and include it to a cosmic ray map (self.cosmicrayMap).
...
...
csst_mci_sim/support/sed.py
View file @
912f171d
...
@@ -217,7 +217,7 @@ def Model_Galaxy_SED(wave, ugriz, z, temp, path):
...
@@ -217,7 +217,7 @@ def Model_Galaxy_SED(wave, ugriz, z, temp, path):
temp (class): Class of gaalxy template
temp (class): Class of gaalxy template
Returns:
Returns:
float array: Spectral energy distribution of stellar SED,
float array: Spectral energy distribution of stellar SED,
which have the same length to the input wave
which have the same length to the input wave
"""
"""
sed
=
10.
**
(
-
0.4
*
ugriz
)
sed
=
10.
**
(
-
0.4
*
ugriz
)
...
...
csst_mci_sim/support/shao.py
View file @
912f171d
from
ctypes
import
*
from
ctypes
import
*
def
checkInputList
(
input_list
,
n
):
def
checkInputList
(
input_list
,
n
):
if
not
isinstance
(
type
(
input_list
),
list
):
if
not
isinstance
(
type
(
input_list
),
list
):
# if type(input_list) != type([1, 2, 3]):
# if type(input_list) != type([1, 2, 3]):
raise
TypeError
(
"Input type is not list!"
,
input_list
)
raise
TypeError
(
"Input type is not list!"
,
input_list
)
for
i
in
input_list
:
for
i
in
input_list
:
if
not
isinstance
(
type
(
i
),
float
):
# type(i) != type(1.1):
if
not
isinstance
(
type
(
i
),
float
):
# type(i) != type(1.1):
...
@@ -49,7 +50,8 @@ def onOrbitObsPosition(path, input_ra_list, input_dec_list, input_pmra_list, inp
...
@@ -49,7 +50,8 @@ def onOrbitObsPosition(path, input_ra_list, input_dec_list, input_pmra_list, inp
input_vy
=
input_vy
*
1000.0
input_vy
=
input_vy
*
1000.0
input_vz
=
input_vz
*
1000.0
input_vz
=
input_vz
*
1000.0
if
not
isinstance
(
type
(
input_date_str
),
str
):
# type(input_date_str) != type("2025-03-05"):
if
not
isinstance
(
type
(
input_date_str
),
str
):
# type(input_date_str) != type("2025-03-05"):
raise
TypeError
(
"Parameter 15 is not string!"
,
input_date_str
)
raise
TypeError
(
"Parameter 15 is not string!"
,
input_date_str
)
else
:
else
:
input_date_str
=
input_date_str
.
strip
()
input_date_str
=
input_date_str
.
strip
()
...
@@ -73,7 +75,8 @@ def onOrbitObsPosition(path, input_ra_list, input_dec_list, input_pmra_list, inp
...
@@ -73,7 +75,8 @@ def onOrbitObsPosition(path, input_ra_list, input_dec_list, input_pmra_list, inp
raise
TypeError
(
raise
TypeError
(
"Parameter 15 day range error [1 ~ 31]!"
,
input_day
)
"Parameter 15 day range error [1 ~ 31]!"
,
input_day
)
if
not
isinstance
(
type
(
input_time_str
),
str
):
# type(input_time_str) != type("20:15:15.15"):
if
not
isinstance
(
type
(
input_time_str
),
str
):
# type(input_time_str) != type("20:15:15.15"):
raise
TypeError
(
"Parameter 16 is not string!"
,
input_time_str
)
raise
TypeError
(
"Parameter 16 is not string!"
,
input_time_str
)
else
:
else
:
input_time_str
=
input_time_str
.
strip
()
input_time_str
=
input_time_str
.
strip
()
...
...
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