Esempio n. 1
0
            continue
        fi = args.f.rstrip('/') + '/' + f

        # READ INPUT DATA FILE
        # --------------------
        this_run_time, this_run_x, this_run_y, this_run_d, this_run_d_err, this_run_t, this_run_h = read_data_file(fi, args.de, args.w)

        # CORRECT FOR NONLINEAR MOTION OF STAGES
        # --------------------------------------
        if args.nc:
            this_run_x, this_run_y = correct_stage_positions(F_CALX, F_CALY, this_run_x, this_run_y)

        # CORRECT FOR TEMPERATURE EXCURSIONS
        # ----------------------------------
        if args.nt:
            this_run_d = correct_for_temperature(F_CORT, this_run_d, this_run_t)

        TIME.append(np.asarray(this_run_time))
        x.append(np.mean(this_run_x))
        y.append(np.mean(this_run_y))
        d.append(np.mean(this_run_d))
        t.append(np.asarray(this_run_t))
        h.append(np.asarray(this_run_h))
        print
else:
    print "using file " + args.f

    # READ INPUT DATA FILE
    # --------------------
    this_run_time, this_run_x, this_run_y, this_run_d, this_run_d_err, this_run_t, this_run_h = read_data_file(args.f, args.de, args.w)
    # CORRECT FOR NONLINEAR MOTION OF STAGES
        print "ignoring file"
        continue
    fi = args.f.rstrip('/') + '/' + f
    # READ INPUT DATA FILE
    # --------------------
    this_time, x, y, d, d_err, t, h = read_data_file(fi, args.de, args.w)

    # CORRECT FOR NONLINEAR MOTION OF STAGES
    # --------------------------------------
    if args.nc:
        x, y = correct_stage_positions(F_CALX, F_CALY, x, y)

    # CORRECT FOR TEMPERATURE EXCURSIONS
    # ----------------------------------
    if args.nt:
        d = correct_for_temperature(F_CORT, d, t)

    run_data_x.append(np.mean(x))
    run_data_y.append(np.mean(y))
    run_data_d.append(np.mean(d))
    run_data_d_err.append(np.mean(d_err))
    print

with open("flat.dat", 'w') as out:
    for i in range(len(run_data_x)):
        if not np.isnan(run_data_x[i]):
            out.write(
                str(round(run_data_x[i], 2)) + '\t' +
                str(round(run_data_y[i], 2)) + '\t' +
                str(round(run_data_d[i], 2)) + '\t' +
                str(round(run_data_d_err[i], 2)) + '\n')
Esempio n. 3
0
        # READ INPUT DATA FILE
        # --------------------
        this_run_time, this_run_x, this_run_y, this_run_d, this_run_d_err, this_run_t, this_run_h = read_data_file(
            fi, args.de, args.w)

        # CORRECT FOR NONLINEAR MOTION OF STAGES
        # --------------------------------------
        if args.nc:
            this_run_x, this_run_y = correct_stage_positions(
                F_CALX, F_CALY, this_run_x, this_run_y)

        # CORRECT FOR TEMPERATURE EXCURSIONS
        # ----------------------------------
        if args.nt:
            this_run_d = correct_for_temperature(F_CORT, this_run_d,
                                                 this_run_t)

        TIME.append(np.asarray(this_run_time))
        x.append(np.mean(this_run_x))
        y.append(np.mean(this_run_y))
        d.append(np.mean(this_run_d))
        t.append(np.asarray(this_run_t))
        h.append(np.asarray(this_run_h))
        print
else:
    print "using file " + args.f

    # READ INPUT DATA FILE
    # --------------------
    this_run_time, this_run_x, this_run_y, this_run_d, this_run_d_err, this_run_t, this_run_h = read_data_file(
        args.f, args.de, args.w)
Esempio n. 4
0
        print "ignoring file"
        continue
    fi = args.f.rstrip('/') + '/' + f
    # READ INPUT DATA FILE
    # --------------------
    this_time, x, y, d, d_err, t, h = read_data_file(fi, args.de, args.w)

    # CORRECT FOR NONLINEAR MOTION OF STAGES
    # --------------------------------------
    if args.nc:
        x, y = correct_stage_positions(F_CALX, F_CALY, x, y)

    # CORRECT FOR TEMPERATURE EXCURSIONS
    # ----------------------------------
    if args.nt:
        d = correct_for_temperature(F_CORT, d, t)

    run_data_x.append(np.mean(x))
    run_data_y.append(np.mean(y))
    run_data_d.append(np.mean(d))
    run_data_d_err.append(np.mean(d_err))
    print

with open("flat.dat", 'w') as out:
    for i in range(len(run_data_x)):
        if not np.isnan(run_data_x[i]):
            out.write(str(round(run_data_x[i], 2)) + '\t' + str(round(run_data_y[i], 2)) + '\t' + str(round(run_data_d[i], 2)) + '\t' + str(round(run_data_d_err[i], 2)) + '\n')

if args.p:
    plt.scatter(run_data_x, run_data_y, c=run_data_d, s=50, linewidth=0, marker='s', cmap="autumn")
    plt.colorbar()