#!/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')