#!/usr/bin/env python
import argparse
from absoluteHeightTask import absoluteHeightTask
parser = argparse.ArgumentParser(description="Perform absolute height analysis on OGP data.")
parser.add_argument("ogp_file", help="OGP data file of xyz tuples containing gauge block and sensor metrology points")
parser.add_argument("sensor_id", help="LSST ID number of sensor")
parser.add_argument("--zoffset", type=float, default=-1.2, help="z-offset")
args = parser.parse_args()
absoluteHeightTask(args.sensor_id, args.ogp_file, zoffset=args.zoffset)
#!/usr/bin/env python
import sys
import siteUtils
import metUtils
from absoluteHeightTask import absoluteHeightTask
sensor_id = siteUtils.getUnitId()
ccd_vendor = siteUtils.getCcdVendor()
if ccd_vendor == 'e2v':
print "Absolute height analysis cannot not be performed with e2v vendor data. Exiting producer script."
sys.exit(0)
patterns = dict(e2v='*CT100*.csv', ITL='*.txt')
#
# Find the vendor metrology scan data from the Data Catalog. Sort the
# results in case of multiple vendor ingests so that we can use the
# most recent one based on SR-RCV-1 job id, assuming the filenames
# are the same for each delivery.
#
met_file = metUtils.get_met_scan_data(sensor_id, patterns[ccd_vendor],
sort=True)[-1]
absoluteHeightTask(sensor_id, met_file, dtype=ccd_vendor,
pickle_file='abs_height.pickle')
direcs = siteUtils.dependency_glob('*',
jobname=siteUtils.getProcessName('OGP_Dimensional_Metrology_Upload'),description='OGP Absolute Height scan:')
print "direcs:"
print direcs
print "------------------------------"
dimdir = ""
for line in direcs :
if "DimensionalMetrology" in line :
dimdir = line
print "dimensional metrology directory = %s" % dimdir
met_file = ""
for root, dirs, met_files in os.walk(dimdir):
for file in met_files:
if file.endswith(".DAT"):
met_file=os.path.join(root, file)
#met_file = met_file[0] # siteUtils returns met_file as a list with one member;
# here take the first (and only) member
print "met_file = %s" % met_file
# Below dtype is always 'OGP', i.e., for this task the source of the data
# is always the OGP scanner at BNL
absoluteHeightTask(sensor_id, met_file, dtype='OGP',
pickle_file='abs_height.pickle')
