Running find_acm_on_ihmp on 20180727 UT data
Here I review the run of find_ihmp_on_acm
on data taken the first night back on sky after July 2018 FPA take-down.
This follows the description of the Dec2017 analysis,
except I have eliminated some of the more verbose detail.
The orignial data file (Data.ALL) and output files from my last run of the
find_acm_on_ihmp can be viewed in the source directory for this web doc
($scohtm/het_eng/Work_HET/Jul27_2018). Sergey made the observations and Matt
collated some of the results using the fplabe.txt file used on 20180727 UT.
Here are my notes from each:
From Sergey:
Remeasure 4 major ACAM positions
BIB 449 402
HPF regular 141 258
LRS2-B 732 212
LRS2-R 735 582
From Matt:
ACAMX ACAMY FPLANE X FPLANEY
448.5 402.2 3.3 -73.2
140.8 257.7 -35.0 10.1
731.6 212.0 -49.8 -148.7
734.5 582.0 49.8 -151.3
If you want the individual measurements then you can use this:
ACAMX ACAMY FPLANE X FPLANE Y
447.97 400.77 3.3 -73.2
448.67 402.9 3.3 -73.2
448.97 402.99 3.3 -73.2
139.36 259.12 -35 10.1
140.69 257.87 -35 10.1
142.41 255.96 -35 10.1
732.95 212.25 -49.8 -148.7
731.16 212.25 -49.8 -148.7
730.57 211.65 -49.8 -148.7
735.49 584.98 49.8 -151.3
733.71 581.41 49.8 -151.3
734.3 579.63 49.8 -151.3
With the values above I can construct the Data.ALL file used as input
for find_ihmp_on_acm.
For our purposes here, the important acm parameters in shuffle.cfg are:
# From the shuffle.cfg file:
# ACAM coordinates for clarity and use in defining candidate stars that fall within the acam
acam_x_origin = 189.1
acam_y_origin = 386.8
acam_pix_scale = 0.2709
acam_x_length = 775
acam_y_length = 773
# section with the angle offsets in degrees for the various parts of the focal
# plane assembly. The offset should be 0 degree for all of them. The overall 90
# degrees rotation is already accounted for in the code
acam = 1.07
I ran find_ihmp_on_acm on the new Data.ALL
file made with the values from Sergey and Matt:
LOCATION : /home/sco/sco/scohtm/scocodes/het_eng/Work_HET/Jul27_2018
[sco@mcs shuf2]$ cat Data.ALL
# IHMP-acm positions after Jul2018 FPA take-down
# name X,Y_ihmp X,Y_acm
# data
555 3.3 -73.2 448.5 402.2
600 -35.0 10.1 140.8 257.7
056 -49.8 -148.7 731.6 212.0
066 49.8 -151.3 734.5 582.0
[sco@mcs shuf2]$ find_ihmp_on_acm Data.ALL
**** Now Computing IHMP Origin on acm ****
# data
179.328 385.420
Plot the results? Y/N: Y
# I use the show() routine to adjust the plot and make a hardcopy.
# The plot file I creat is: acm_20180227.png
I can see all residuals in: Final_Residuals.List
---------------------------------------------------------------------------------
Below I show the two files I inspect to get the acm center errors and
to derive the new "acam" value (the rotation ange used by shuffle).
% cat Final_Residuals.List
Residuals List, coordinates are in units of file: acm.xy
X1,Y1 = original coordinates from acm.xy
X2,Y2 = coordinates from ihmp.xy transformed to the acm.xy system
# Output from xy_tranrot_res.sh run:
# Mean dX (X1-X2) = -1.0480 -+ 0.9075
# Mean dX (Y1-Y2) = 0.9710 -+ 0.7140
# dX sigma = 1.8150
# dY sigma = 1.4280
# dR (X,Y) sigma = 1.4280
# Number of residuals = 4
# List of XY residuals from xy_tranrot_res.sh
line X1 Y1 X2 Y2 dX dY Res Res_norm
# data
1 448.500 402.200 449.252 402.705 -0.752 -0.505 0.906 0.3922
2 140.800 257.700 144.494 255.542 -3.694 2.158 4.278 1.8525
3 731.600 212.000 731.600 212.000 0.000 0.000 0.000 0.0000
4 734.500 582.000 734.246 579.769 0.254 2.231 2.245 0.9723
The mean residuals statistics were:
RMS of X residuals (file1-file2) = 1.815
RMS of Y residuals (file1-file2) = 1.428
Number of residuals points used = 4
We could also use the rotation angle measure:
[sco@mcs shuf2]$ cat trs_rotatedet.explain
Final values (Theta_deg,err,MeanRes_pixels,Npoints):
268.917 0.0833 2.4767 3
268.807 0.3632 4.6793 3 (Dec2017 results)
Hence, I would guess: acam = 270.00-268.92 = 1.08 -+ 0.1 degrees
The plot of the original acm coordinated (F1 Original = the small blue squares)
and the transformed IHMP coordinates (F2 Final = the large red squares) is shown
below.
 |
The IHMP to acm fit made with the positions (acm coordinate frame in pixel units) for
056 (LRS2-B), 066(LRS2-R), 555(BIB), and 600(HPFacq). In this case the small blue
squares are the original acm positions. The larger red squares are the IHMP positions
transformed to the acm system using the results from our trs_solve_2 run. These results
are contained in the local file named TRS.final. Below are the (1) input file to the
processing script find_ihmp_on_acm, (2) the command line call to make the run, and
(3) the new acm parameters for shuffle:
% cat Data.ALL # the input data file
# IHMP-acm positions after Jul2018 FPA take-down
# name X,Y_ihmp X,Y_acm
# data
555 3.3 -73.2 448.5 402.2
600 -35.0 10.1 140.8 257.7
056 -49.8 -148.7 731.6 212.0
066 49.8 -151.3 734.5 582.0
% find_ihmp_on_acm Data.ALL # the command line
# New shuffle parameters
Early Dec2017
acam_x_origin = 179.3 -+1.8 (units of acm pixels) 182.7 pixels
acam_y_origin = 385.4 -+1.4 (units of acm pixels) 384.8 pixels
acam = 1.08 -+0.1 (degrees) 1.20 degrees
|
The green values above were compute in early Dec2017 just prior that engineering
run. The revised Dec2017 shuffle values presented only small changes, as do the
newly derived paramters for out Jul2018 data.
Back to calling page