"""

controlling function to run the script

input: none

output: /data/mta/www/mta_grat/<dir> where dir is <Mon><yy>

"""

#

#--- set data collecting range

#

(start, stop, dir) = find_time_interval()

#

#--- extract needed data fits files

#

run_arc4gl(start, stop)

#

#--- process fits files

#

run_idl(dir)

#

#--- update records

#

update_record(dir)

#

#--- move the data to an appropriate location

#

cmd = 'mv /data/mta4/Gratings/' + dir + ' /data/mta/www/mta_grat/.'

os.system(cmd)

#

#--- clean up

#

os.system('rm -rf *.fits.gz')

"""

extract acis and hrc evt1a.fits files using arc4gl

input: start --- start time in the format of 03/01/15

stop --- stop time

output: fits files (e.g., acisf17108_001N002_evt1a.fits.gz)

"""

#

#--- read a template and create the current command file

#

file = bin_dir + 'house_keeping/arc_template'

line = open(file, 'r').read()

line = line.replace('#START#', start)

line = line.replace('#STOP#', stop)

fo = open('./run_arc', 'w')

fo.write(line)

fo.close()

#

#--- run arc4gl

#

cmd1 = "/usr/bin/env PERL5LIB="

cmd2 = ' echo ' + hakama + ' |arc4gl -U' + dare + ' -Sarcocc -i./run_arc '

cmd = cmd1 + cmd2

bash(cmd, env=ascdsenv)

#

#--- remove unwanted fits files

#

"""

process fits files with an updated idl scripts

input: dir --- the name of output directory

output: dir/<stemp> --- a directory which contains the processed data

"""

#

#--- read a template and create the current command file

#

file = bin_dir + 'house_keeping/pl_template'

line = open(file, 'r').read()

line = line.replace('#DIR#', dir)

fo = open('./mk_idl_command.pl', 'w')

fo.write(line)

fo.close()

#

#--- make an output directory

#

cmd = 'mkdir /data/mta4/Gratings/' + dir

os.system(cmd)

#

#--- run a perl script to create an idl script

#

cmd1 = "/usr/bin/env PERL5LIB="

cmd2 = ' perl ./mk_idl_command.pl '

cmd = cmd1 + cmd2

bash(cmd, env=ascdsenv)

#

#--- run the idl script to process fits files

#

"""

create start and stop time and output directory name

input: test --- if it is test, yes. default: 'NO'

month --- test month input. default: 4. it is ignored if not test.

year --- test year input. default: 2015. it is ignored if not test.

output: start --- start time in the format of mm/01/yy (e.g. 03/01/15)

stop --- stop time

dir --- output directory name in the format of <Mon><yy> (e.g. Jan15)

"""

#

#--- for the case if you want to run the test, there is option...

#

if test == 'NO':

#

#--- this function gives you the current time. year in the format of 2015, and month in 4

#

[year, month, day, hours, min, sec, weekday, yday, dst] = tcnv.currentTime()

#

#--- we need start and stop times in the formats of 03/01/15, 04/01/15

#

year_start = year

year_end = year

this_month = str(month)

if month < 10:

this_month = '0' + this_month

temp_month = month -1

if temp_month < 1:

temp_month = 12

year_start -= 1

last_month = str(temp_month)

if temp_month < 10:

last_month = '0' + last_month

temp = str(year_start)

year_start = temp[2] + temp[3]

temp = str(year_end)

year_end = temp[2] + temp[3]

start = last_month + '/01/' + year_start

stop = this_month + '/01/' + year_end

#

#--- create the output directory name (format example: Apr15)

#

cmonth = tcnv.changeMonthFormat(temp_month)

dir = cmonth + year_start

"""

update the recormed of the process data

input: none but read from the current output directory

outpu: obslist in 'Angles', 'EdE', 'Focus', 'Zero' directories

"""

#

#--- find out which data are processed from the output directory

#

cmd = 'ls /data/mta4/Gratings/' + dir + '> ./ztemp'

os.system(cmd)

f = open('./ztemp', 'r')

data = [line.strip() for line in f.readlines()]

f.close()

os.system('rm ./ztemp')

line = ''

for each in data:

line = line + '../' + dir + '/' + each + '\n'

#

#--- there are four directories with "obslist"

#

for odir in ('Angles', 'EdE', 'Focus', 'Zero'):

file = '/data/mta/www/mta_grat/' + odir + '/obslist'

fo = open(file, 'a')

fo.write(line)

"""

testing functions

"""

#------------------------------------------------------------

def test_find_time_interval(self):

out = find_time_interval(test='YES')

comp = ('03/01/15', '04/01/15', 'Mar15')