Running find_acm_on_ihmp on 20180929 UT data
oct2018_notes.html
Special Note: A lot of re-evaluation occurred
during this period. You can read about some of this
in my Oct2018 notes.
Here I review the run of find_ihmp_on_acm
on data taken the first night back on sky after Sep 2018 FPA take-down. 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/Sep29_2018).
Steven J. made the observations on 20180929 UT.
From Steven J.:
Remeasure 4 major ACAM positions
LRS2B - 735.9, 212.2
LRS2R - 738.2, 581.8
HPF regular setup - 145.7, 254.2
HPF direct setup - 127.9, 292.5
BIB - 454.8, 402.3
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
With the values above I can construct the Data.ALL file used as input
for find_ihmp_on_acm.
Note that for this run I have added the PREVIOUS ACAM VALUES to
the Data.ALL file so that we can consistently produce some "quick viusual checks"
before deriving the new parametrs. In other words, we have one table
where we can easily compare the new and old acm positions.
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 SJ (20180929 UT):
LOCATION : $scohtm/scocodes/het_eng/Work_HET/Sep29_2018
[sco@mcs shuf2]$ cat Data.ALL
# IHMP-acm positions after Sep2018 FPA take-down
# NOTE: I am changing find_ihmp_on_acm to read the PREVIOUS X,Y_acm values
# name X,Y_ihmp X,Y_acm X,Y_acm _prev
# data
555 3.3 -73.2 454.8 402.3 448.5 402.2
600 -35.0 10.1 145.7 254.2 140.8 257.7
056 -49.8 -148.7 735.9 212.2 731.6 212.0
066 49.8 -151.3 738.2 581.8 734.5 582.0
# Here I run the code and show SOME of the important output
[sco@mcs shuf2]$ find_ihmp_on_acm Data.ALL
**** Now Computing IHMP Origin on acm ****
# data
183.374 384.810
Plot the results? Y/N: Y
# I use the show() routine to adjust the plot and make a hardcopy.
# The plot file I create is: acm_20180929.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) = -0.3278 -+ 0.9629
# Mean dX (Y1-Y2) = 0.2390 -+ 0.5487
# dX sigma = 1.9258
# dY sigma = 1.0975
# dR (X,Y) sigma = 1.0975
# 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 454.800 402.300 453.273 402.491 1.527 -0.191 1.539 0.6943
2 145.700 254.200 148.731 254.881 -3.031 -0.681 3.107 1.4015
3 735.900 212.200 735.900 212.200 0.000 0.000 0.000 0.0000
4 738.200 581.800 738.007 579.972 0.193 1.828 1.838 0.8293
The mean residuals statistics were:
RMS of X residuals (file1-file2) = 1.926
RMS of Y residuals (file1-file2) = 1.097
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.833 0.0833 2.1612 3
268.917 0.0833 2.4767 3 # (Jul2018 Results)
268.807 0.3632 4.6793 3 # (Dec2017 results)
Hence, I would guess for NOW: acam = 270.00-268.83 = 1.17 -+ 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 (in 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 Sep2018 FPA take-down
# NOTE: I am changing find_ihmp_on_acm to read the PREVIOUS X,Y_acm values
# name X,Y_ihmp X,Y_acm X,Y_acm _prev
# data
555 3.3 -73.2 454.8 402.3 448.5 402.2
600 -35.0 10.1 145.7 254.2 140.8 257.7
056 -49.8 -148.7 735.9 212.2 731.6 212.0
066 49.8 -151.3 738.2 581.8 734.5 582.0
% find_ihmp_on_acm Data.ALL # the command line
# New shuffle parameters
Late July2018
acam_x_origin = 183.4 ± 1.8 (units of acm pixels) 179.3 ± 1.8 pixels
acam_y_origin = 384.8 ± 1.1 (units of acm pixels) 385.4 ± 1.4 pixels
acam = 1.17 ± 0.1 (degrees) 1.08 ± 0.1 degrees
|
The green values above were computed after the July 2018 FPA takedown
for comparison. The observed changes are fairly small.
UPDATE: New HPS ihmp coordinated should be used now!
Hi Steve,
Thanks a lot.
It looks like the coordinates in the fplane.txt file were not updated in February as I thought they were.
Could you check whether using the update coordinates I sent in previous email improves the calculation?
(Two relevant lines form the previous email copy pasted below)
600 -35.8 12.3 902 902 902 0.0 1.0
603 -25.9 17.5 905 905 905 0.0 1.0
Hence, my updataed file should be:
% cat Data.ALL # the input data file
# IHMP-acm positions after Sep2018 FPA take-down
# NOTE: I am changing find_ihmp_on_acm to read the PREVIOUS X,Y_acm values
# NOTE: Using Joe Ninan's Sep30,2018 reported ihmp values for HPF positions
# name X,Y_ihmp X,Y_acm X,Y_acm _prev
# data
555 3.3 -73.2 454.8 402.3 448.5 402.2
600 -35.8 12.3 145.7 254.2 140.8 257.7
056 -49.8 -148.7 735.9 212.2 731.6 212.0
066 49.8 -151.3 738.2 581.8 734.5 582.0
The position for 603 I should use in find_ihmp_on_acm: x,y_603 = -25.9,17.5
****After these changes I get:
**** Now Computing IHMP Origin on acm ****
183.809 386.198
**** Now Computing HPFacq_600 on acm ****
140.773 253.259
**** Now Computing HPFdirect_603 on acm ****
120.929 289.455
acam_x_origin = 183.809 -+ 2.5 acm pixels
acam_y_origin = 386.198 -+ 1.2 acm pixels
acam = 270 - 268.977 = 1.02 -+ 0.08 degrees
Predicted HPFacq_600 position from fit = 140.773 253.259 (acm pixels)
Predicted HPFacq_603 position from fit = 120.929 289.455 (acm pixels)
UPDATE: New data from SJ on Oct1,2018
I am now performing a new run with Data.ALL_oct01.
NOTE from SJ the morning of Oct01,2018:
BIB 453.0(1.0) 401.5(0.5)
LRS2-B 735.5(0.4) 212.6(0.4)
LRS2-R 737.9(0.2) 581.7(0.4)
HPF side CFB 217.9(0.2) 182.2(0.1)
HPFACQ 126.4(0.2) 292.2(0.3)
*** Looks like the HPFACQ value is the HPFdirect value!!!! This is a problem
For now I use the old (Sep29 value for HPFACQ=600)
Here is the new dile (/home/sco/sco/scohtm/scocodes/het_eng/Work_HET/Sep29_2018/S/Data.ALL_oct01)
# IHMP-acm positions after Sep2018 FPA take-down
# NOTE: Using Joe Ninan's Sep30,2018 reported ihmp values for HPF positions
# NOTE: These acm positions from SJ Oct01,2018
# name X,Y_ihmp X,Y_acm X,Y_acm _prev
# data
555 3.3 -73.2 453.0 401.5 448.5 402.2
056 -49.8 -148.7 735.4 212.6 731.6 212.0
066 49.8 -151.3 737.9 581.7 734.5 582.0
610 -56.6 -6.6 217.9 182.2 nan nan
# The old HPFacq line messes the solution up....
600 -35.8 12.3 145.7 254.2 140.8 257.7
Commands I run in ..../scocodes/het_eng/Work_HET/Sep29_2018/run_oct01
find_ihmp_on_acm Data.ALL_oct01
cat Final_Residuals.List
cat trs_rotatedet.explain
acam = 270.00-ang where ang is angle from trs_rotatedet.explain
NOTE: When I run with only the points: 555,056,066,0610
*** The 610 point never fits in very well!
*** May want to make find_ihmp_on_acm plot only one set of Text id.
The code that makes the plot is: trs_plot.py
Back to calling page