def hrc_dose_plot_monthly_report(indir='NA', outdir='NA'):
"""
read hrc database, and plot history of exposure.
input: indir --- data directory path
outdir --- output directory path
output: <outdir>/hrc_max_exp.gif
"""
#
#--- setting indir and outdir if not given
#
if indir == 'NA':
indir = data_out
if outdir == 'NA':
outdir = img_dir
#
#--- read HRC I data
#
out = expf.readExpData(indir, 'hrci')
idate = out[0]
imax_acc = out[7]
#
#--- read HRC S data
#
out = expf.readExpData(indir, 'hrcs')
smax_acc = out[7]
#
#--- plot data
#
plot_max_dose(idate, imax_acc, smax_acc)
#
#--- move the plot to img directory
#
cmd = 'mv hrc_max_exp.gif ' + outdir
os.system(cmd)
def hrc_dose_make_data_html(indir='NA', outdir='NA'):
"""
read hrc database, and create html page:
input: indir --- data directory
outdir --- output direcotry
"""
if indir == 'NA':
indir = data_out
if outdir == 'NA':
outdir = data_out
for hrc in ('hrci', 'hrcs'):
#
#--- just in a case, clear up the files before reading them
#
expf.clean_data(indir)
#
#--- read HRC histrical data
#
data = expf.readExpData(indir, hrc)
#
#--- create a HTML page to display histrical data
#
print_html_page(indir, outdir, hrc, data)
def acis_dose_make_data_html(indir='NA', outdir='NA'):
"""
read data and create html pages
input: indir --- input data directory
outdir --- html directory
"""
#
#--- setting indir and outdir if not given
#
if indir == 'NA':
indir = data_out
if outdir == 'NA':
outdir = data_out
#
#--- read data
#
for ccd in ('i_2', 'i_3', 's_2', 's_3'):
for sec in range(0, 4):
inst = ccd + '_n_' + str(sec)
data = expf.readExpData(indir, inst)
#
#--- write html page
#
outfile = outdir + inst + '.html'
write_html(ccd, sec, data, outfile)
#
#--- update top html page (this is done in HRC)
#
# line = house_keeping + 'exposure.html'
# with open(line, 'r') as f:
# data = f.read()
#
# now = time.strftime("%m %d %Y", time.gmtime())
# [syear, smon, day] = mcf.today_date()
# smon -= 1
#
# if smon < 1:
# smon = 12
# syear -= 1
# lyear = str(syear)
# lmon = mcf.add_leading_zero(smon)
#
# data = data.replace("#DATE#", now)
# data = data.replace("#YEAR#", lyear)
# data = data.replace("#MON#", lmon)
#
# ofile = web_dir + 'exposure.html'
# with open(ofile, 'w') as fo:
# fo.write(data)
#
#--- update plot page htmls
#
update_plt_html_date()
def hrc_dose_plot_exposure_stat(indir='NA', outdir='NA', clean='NA'):
"""
read hrc database, and plot history of exposure.
input: indir --- data directory path
outdir --- output directory path
output: outdir/<hrc>.png
"""
#
#--- setting indir and outdir if not given
#
if indir == 'NA':
indir = data_out
if outdir == 'NA':
outdir = plot_dir
#
#--- clean up the data sets before reading
#
if clean != 'NA':
expf.clean_data(indir)
#
#--- start plotting data
#
for detector in ('hrci', 'hrcs'):
#
#--- read data
#
idata = expf.readExpData(indir, detector)
#
#--- plot data
#
try:
plot_hrc_dose(idata)
#
#--- trim the edge of the plot and move to the plot directory
#
outfile = detector + '.png'
cmd = 'convert hrc.png -trim ' + plot_dir + detector + '.png'
os.system(cmd)
mcf.rm_files('hrc.png')
except:
pass
def acis_dose_plot_exposure_stat(indir='NA', outdir='NA', clean='NA'):
"""
read acis database, and plot history of exposure.
input: indir --- input directory
outdir --- output directory
clean --- if not 'NA', clean up the data sets before plotting the data
output: <outdir>/<inst>.png
"""
#
#--- setting indir and outdir if not given
#
if indir == 'NA':
indir = data_out
if outdir == 'NA':
outdir = plot_dir
#
#--- clean up the data sets before reading
#
if clean != 'NA':
expf.clean_data(indir)
#
#--- start plotting data
#
for ccd in ('i_2', 'i_3', 's_2', 's_3'):
for sec in range(0, 4):
inst = ccd + '_n_' + str(sec)
idata = expf.readExpData(indir, inst)
#
#--- plot data
#
plot_acis_dose(idata)
#
#--- trim the edge and move to the plot directory
#
cmd = 'convert acis.png -trim ' + outdir + inst + '.png'
os.system(cmd)
mcf.rm_files('rm acis.png')
def hrc_dose_make_data_html(indir = 'NA', outdir = 'NA', indir2 = 'NA', outdir2 = 'NA', comp_test = 'NA'):
'read hrc database, and create html page: input: indir (data directory)'
#
#--- setting indir and outdir if not given
#
if comp_test == 'test':
indir = test_data_out
outdir = test_data_out
indir2 = test_data_out_hrc
outdir2 = test_data_out_hrc
else:
if indir == 'NA':
indir = data_out
if outdir == 'NA':
outdir = data_out
if indir2 == 'NA':
indir2 = data_out_hrc
if outdir2 == 'NA':
outdir2 = data_out_hrc
for hrc in ('hrci', 'hrcs'):
#
#--- center exposure map --------
#
#
#--- just in a case, clear up the files before reading them
#
expf.clean_data(indir)
#
#--- read HRC histrical data
#
[date,year,month,mean_acc,std_acc,min_acc,min_apos, max_acc,max_apos,asig1, asig2, asig3, mean_dff,std_dff, \
min_dff, min_dpos,max_dff,max_dpos,dsig1, dsig2, dsig3] = expf.readExpData(indir,hrc)
#
#--- create a HTML page to display histrical data
#
printHtml(indir, outdir, hrc, date, year,month,mean_acc,std_acc,min_acc,min_apos,max_acc,max_apos, asig1, asig2, asig3, \
mean_dff,std_dff,min_dff,min_dpos,max_dff,max_dpos, dsig1, dsig2, dsig3)
#
#--- full exposure map --------
#
#
#--- just in a case, clear up the files before reading them
#
# expf.clean_data(indir2)
for i in range(0,10):
if hrc == 'hrci' and i == 9:
break
hrc2 = hrc + '_' + str(i) #--- naming is slightly different from the center exposure maps
#
#--- read HRC histrical data --- sec='full' indicating reading full exposure maps
#
[date,year,month,mean_acc,std_acc,min_acc,min_apos, max_acc,max_apos,asig1, asig2, asig3, mean_dff,std_dff, \
min_dff, min_dpos,max_dff,max_dpos,dsig1, dsig2, dsig3] = expf.readExpData(indir2, hrc2)
#
#--- create a HTML page to display histrical data
#
printHtml(indir2, outdir2, hrc2, date, year,month,mean_acc,std_acc,min_acc,min_apos,max_acc,max_apos, asig1, asig2, asig3, \
mean_dff,std_dff,min_dff,min_dpos,max_dff,max_dpos, dsig1, dsig2, dsig3)
def acis_dose_make_data_html(indir = 'NA', outdir = 'NA', comp_test = 'NA'):
"""
read data and create html pages
input: indir, outdir, both are optional
"""
#
#--- setting indir and outdir if not given
#
if indir == 'NA':
if comp_test == 'test':
indir = test_data_out
else:
indir = data_out
if outdir == 'NA':
if comp_test == 'test':
outdir = test_data_out
else:
outdir = data_out
#
#--- read data
#
for ccd in ('i_2', 'i_3', 's_2', 's_3'):
for sec in range(0, 4):
inst = ccd + '_n_' + str(sec)
date = []
year = []
month = []
mean_acc = []
std_acc = []
min_acc = []
min_apos = []
max_acc = []
max_apos = []
asig1 = []
asig2 = []
asig3 = []
mean_dff = []
std_dff = []
min_dff = []
min_dpos = []
max_dff = []
max_dpos = []
dsig1 = []
dsig2 = []
dsig3 = []
[date,year,month,mean_acc,std_acc,min_acc,min_apos, max_acc,max_apos,asig1, asig2, asig3, mean_dff,std_dff, \
min_dff, min_dpos,max_dff,max_dpos,dsig1, dsig2, dsig3] = expf.readExpData(indir, inst)
#
#--- write html page
#
write_html(ccd, sec, year,month,mean_acc,std_acc,min_acc,min_apos, max_acc,max_apos,asig1, asig2, asig3, \
mean_dff,std_dff,min_dff, min_dpos,max_dff,max_dpos,dsig1, dsig2, dsig3)
out_name = outdir + inst + '.html'
cmd = 'mv acis.html ' + out_name
os.system(cmd)
#
#--- update top html page
#
line = house_keeping + 'exposure.html'
f = open(line, 'r')
data = f.read()
f.close()
now = tcnv.currentTime("Display")
ndate = 'Last Update: ' + now
ndata = data.replace("Last Update:", ndate)
line = web_dir + 'exposure.html'
f = open(line, 'w')
f.write(ndata)
f.close()
#
#--- update plot page htmls
#
update_plt_html_date(comp_test)
def hrc_dose_plot_monthly_report(indir='NA', outdir='NA', comp_test='NA'):
'read hrc database, and plot history of exposure. input: data directory path, output directory path '
#
#--- setting indir and outdir if not given
#
if comp_test == 'test':
indir = data_out
outdir = test_img_dir
else:
if indir == 'NA':
indir = data_out
if outdir == 'NA':
outdir = img_dir
#
#--- read HRC I data
#
[idate, iyear,imonth,imean_acc,istd_acc,imin_acc,imin_apos, imax_acc,imax_apos,iasig1, iasig2, iasig3, \
imean_dff,istd_dff,imin_dff, imin_dpos,imax_dff,imax_dpos,idsig1, idsig2, idsig3] = expf.readExpData(indir,'hrci')
#
#--- read HRC S data
#
[sdate, syear,smonth,smean_acc,sstd_acc,smin_acc,smin_apos, smax_acc,smax_apos,sasig1, sasig2, sasig3, \
smean_dff,sstd_dff,smin_dff, smin_dpos,smax_dff,smax_dpos,sdsig1, sdsig2,sdsig3] = expf.readExpData(indir,'hrcs')
#
#--- plot data
#
plot_max_dose(idate, imax_acc, smax_acc)
#
#--- move the plot to img directory
#
cmd = 'mv hrc_max_exp.gif ' + outdir
os.system(cmd)
def acis_dose_plot_exposure_stat(indir='NA',
outdir='NA',
indir2='NA',
outdir2='NA',
clean='NA'):
'read acis database, and plot history of exposure. input: data directory path, output directory path '
#
#--- setting indir and outdir if not given
#
if indir == 'NA':
indir = data_out
if outdir == 'NA':
outdir = plot_dir
#
#--- clean up the data sets before reading
#
if clean != 'NA':
expf.clean_data(indir)
#
#--- start plotting data
#
for ccd in ('i_2', 'i_3', 's_2', 's_3'):
for sec in range(0, 4):
inst = ccd + '_n_' + str(sec)
date = []
year = []
month = []
mean_acc = []
std_acc = []
min_acc = []
min_apos = []
max_acc = []
max_apos = []
m10_acc = []
m10_apos = []
mean_dff = []
std_dff = []
min_dff = []
min_dpos = []
max_dff = []
max_dpos = []
m10_dff = []
m10_dpos = []
expf.readExpData(indir, inst, date, year,month,mean_acc,std_acc,min_acc,min_apos, max_acc,max_apos,m10_acc, \
m10_apos,mean_dff,std_dff,min_dff, min_dpos,max_dff,max_dpos,m10_dff,m10_dpos)
#
#--- plot data
#
plot_acis_dose(date, mean_dff, max_dff, m10_dff, mean_acc, max_acc,
m10_acc)
#
#--- move to the plot directory
#
outfile = inst + '.gif'
cmd = 'convert acis.png ' + plot_dir + outfile
os.system(cmd)
os.system('rm acis.png')
def acis_dose_make_data_html(indir='NA', outdir='NA', comp_test='NA'):
"""
read data and create html pages
input: indir, outdir, both are optional
"""
#
#--- setting indir and outdir if not given
#
if indir == 'NA':
if comp_test == 'test':
indir = test_data_out
else:
indir = data_out
if outdir == 'NA':
if comp_test == 'test':
outdir = test_data_out
else:
outdir = data_out
#
#--- read data
#
for ccd in ('i_2', 'i_3', 's_2', 's_3'):
for sec in range(0, 4):
inst = ccd + '_n_' + str(sec)
date = []
year = []
month = []
mean_acc = []
std_acc = []
min_acc = []
min_apos = []
max_acc = []
max_apos = []
asig1 = []
asig2 = []
asig3 = []
mean_dff = []
std_dff = []
min_dff = []
min_dpos = []
max_dff = []
max_dpos = []
dsig1 = []
dsig2 = []
dsig3 = []
[date,year,month,mean_acc,std_acc,min_acc,min_apos, max_acc,max_apos,asig1, asig2, asig3, mean_dff,std_dff, \
min_dff, min_dpos,max_dff,max_dpos,dsig1, dsig2, dsig3] = expf.readExpData(indir, inst)
#
#--- write html page
#
write_html(ccd, sec, year,month,mean_acc,std_acc,min_acc,min_apos, max_acc,max_apos,asig1, asig2, asig3, \
mean_dff,std_dff,min_dff, min_dpos,max_dff,max_dpos,dsig1, dsig2, dsig3)
out_name = outdir + inst + '.html'
cmd = 'mv acis.html ' + out_name
os.system(cmd)
#
#--- update top html page
#
line = house_keeping + 'exposure.html'
f = open(line, 'r')
data = f.read()
f.close()
now = time.strftime("%m %d %Y", time.gmtime())
smon = int(time.strftime("%m", time.gmtime())) - 1
syear = int(time.strftime("%Y", time.gmtime()))
if smon < 1:
smon = 12
syear -= 1
lmon = str(smon)
if smon < 10:
lmon = '0' + lmon
lyear = str(syear)
data = data.replace("#DATE#", now)
data = data.replace("#YEAR#", lyear)
data = data.replace("#MON#", lmon)
line = web_dir + 'exposure.html'
f = open(line, 'w')
f.write(data)
f.close()
#
#--- update plot page htmls
#
update_plt_html_date(comp_test)
def acis_dose_make_data_html(indir='NA', outdir='NA', comp_test='NA'):
"""
read data and create html pages
input: indir, outdir, both are optional
"""
#
#--- setting indir and outdir if not given
#
if indir == 'NA':
if comp_test == 'test':
indir = test_data_out
else:
indir = data_out
if outdir == 'NA':
if comp_test == 'test':
outdir = test_data_out
else:
outdir = data_out
#
#--- read data
#
for ccd in ('i_2', 'i_3', 's_2', 's_3'):
for sec in range(0, 4):
inst = ccd + '_n_' + str(sec)
date = []
year = []
month = []
mean_acc = []
std_acc = []
min_acc = []
min_apos = []
max_acc = []
max_apos = []
asig1 = []
asig2 = []
asig3 = []
mean_dff = []
std_dff = []
min_dff = []
min_dpos = []
max_dff = []
max_dpos = []
dsig1 = []
dsig2 = []
dsig3 = []
[date,year,month,mean_acc,std_acc,min_acc,min_apos, max_acc,max_apos,asig1, asig2, asig3, mean_dff,std_dff, \
min_dff, min_dpos,max_dff,max_dpos,dsig1, dsig2, dsig3] = expf.readExpData(indir, inst)
#
#--- write html page
#
write_html(ccd, sec, year,month,mean_acc,std_acc,min_acc,min_apos, max_acc,max_apos,asig1, asig2, asig3, \
mean_dff,std_dff,min_dff, min_dpos,max_dff,max_dpos,dsig1, dsig2, dsig3)
out_name = outdir + inst + '.html'
cmd = 'mv acis.html ' + out_name
os.system(cmd)
#
#--- update top html page
#
line = house_keeping + 'exposure.html'
f = open(line, 'r')
data = f.read()
f.close()
now = tcnv.currentTime("Display")
ndate = 'Last Update: ' + now
ndata = data.replace("Last Update:", ndate)
line = web_dir + 'exposure.html'
f = open(line, 'w')
f.write(ndata)
f.close()
#
#--- update plot page htmls
#
update_plt_html_date(comp_test)
def hrc_dose_make_data_html(indir='NA',
outdir='NA',
indir2='NA',
outdir2='NA',
comp_test='NA'):
'read hrc database, and create html page: input: indir (data directory)'
#
#--- setting indir and outdir if not given
#
if comp_test == 'test':
indir = test_data_out
outdir = test_data_out
indir2 = test_data_out_hrc
outdir2 = test_data_out_hrc
else:
if indir == 'NA':
indir = data_out
if outdir == 'NA':
outdir = data_out
if indir2 == 'NA':
indir2 = data_out_hrc
if outdir2 == 'NA':
outdir2 = data_out_hrc
for hrc in ('hrci', 'hrcs'):
#
#--- center exposure map --------
#
#
#--- just in a case, clear up the files before reading them
#
expf.clean_data(indir)
#
#--- read HRC histrical data
#
[date,year,month,mean_acc,std_acc,min_acc,min_apos, max_acc,max_apos,asig1, asig2, asig3, mean_dff,std_dff, \
min_dff, min_dpos,max_dff,max_dpos,dsig1, dsig2, dsig3] = expf.readExpData(indir,hrc)
#
#--- create a HTML page to display histrical data
#
printHtml(indir, outdir, hrc, date, year,month,mean_acc,std_acc,min_acc,min_apos,max_acc,max_apos, asig1, asig2, asig3, \
mean_dff,std_dff,min_dff,min_dpos,max_dff,max_dpos, dsig1, dsig2, dsig3)
#
#--- full exposure map --------
#
#
#--- just in a case, clear up the files before reading them
#
# expf.clean_data(indir2)
for i in range(0, 10):
if hrc == 'hrci' and i == 9:
break
hrc2 = hrc + '_' + str(
i
) #--- naming is slightly different from the center exposure maps
#
#--- read HRC histrical data --- sec='full' indicating reading full exposure maps
#
[date,year,month,mean_acc,std_acc,min_acc,min_apos, max_acc,max_apos,asig1, asig2, asig3, mean_dff,std_dff, \
min_dff, min_dpos,max_dff,max_dpos,dsig1, dsig2, dsig3] = expf.readExpData(indir2, hrc2)
#
#--- create a HTML page to display histrical data
#
printHtml(indir2, outdir2, hrc2, date, year,month,mean_acc,std_acc,min_acc,min_apos,max_acc,max_apos, asig1, asig2, asig3, \
mean_dff,std_dff,min_dff,min_dpos,max_dff,max_dpos, dsig1, dsig2, dsig3)
def hrc_dose_plot_monthly_report(indir = 'NA', outdir = 'NA', comp_test = 'NA'):
'read hrc database, and plot history of exposure. input: data directory path, output directory path '
#
#--- setting indir and outdir if not given
#
if comp_test == 'test':
indir = data_out
outdir = test_img_dir
else:
if indir == 'NA':
indir = data_out
if outdir == 'NA':
outdir = img_dir
#
#--- read HRC I data
#
[idate, iyear,imonth,imean_acc,istd_acc,imin_acc,imin_apos, imax_acc,imax_apos,iasig1, iasig2, iasig3, \
imean_dff,istd_dff,imin_dff, imin_dpos,imax_dff,imax_dpos,idsig1, idsig2, idsig3] = expf.readExpData(indir,'hrci')
#
#--- read HRC S data
#
[sdate, syear,smonth,smean_acc,sstd_acc,smin_acc,smin_apos, smax_acc,smax_apos,sasig1, sasig2, sasig3, \
smean_dff,sstd_dff,smin_dff, smin_dpos,smax_dff,smax_dpos,sdsig1, sdsig2,sdsig3] = expf.readExpData(indir,'hrcs')
#
#--- plot data
#
plot_max_dose(idate, imax_acc, smax_acc)
#
#--- move the plot to img directory
#
cmd = 'mv hrc_max_exp.gif ' + outdir
os.system(cmd)
def acis_dose_plot_exposure_stat(indir = 'NA', outdir = 'NA', indir2 = 'NA', outdir2 = 'NA', clean = 'NA'):
'read acis database, and plot history of exposure. input: data directory path, output directory path '
#
#--- setting indir and outdir if not given
#
if indir == 'NA':
indir = data_out
if outdir == 'NA':
outdir = plot_dir
#
#--- clean up the data sets before reading
#
if clean != 'NA':
expf.clean_data(indir)
#
#--- start plotting data
#
for ccd in ('i_2', 'i_3', 's_2', 's_3'):
for sec in range(0, 4):
inst = ccd + '_n_' + str(sec)
date = []
year = []
month = []
mean_acc = []
std_acc = []
min_acc = []
min_apos = []
max_acc = []
max_apos = []
m10_acc = []
m10_apos = []
mean_dff = []
std_dff = []
min_dff = []
min_dpos = []
max_dff = []
max_dpos = []
m10_dff = []
m10_dpos = []
expf.readExpData(indir, inst, date, year,month,mean_acc,std_acc,min_acc,min_apos, max_acc,max_apos,m10_acc, \
m10_apos,mean_dff,std_dff,min_dff, min_dpos,max_dff,max_dpos,m10_dff,m10_dpos)
#
#--- plot data
#
plot_acis_dose(date, mean_dff, max_dff, m10_dff, mean_acc, max_acc, m10_acc)
#
#--- move to the plot directory
#
outfile = inst + '.gif'
cmd = 'convert acis.png ' + plot_dir + outfile
os.system(cmd)
os.system('rm acis.png')