December 2007 track map

Steps required to use Kim's Track Map processor in Matlab

  • Identify the start time and stop time of the total dataset (CW+CCW scans) in HH:MM:SS and in MJD
  • For example, for 2007-12-10, these values are:
    • 04:56:55 to 09:48:29
    • 54444.3073032407 to 54444.4086689815
  • Assemble the accelerometer, inclinometer and antenna position FITS files from /home/gbtlogs into subdirectories:
    • The files must cover the time range of the entire track scan.
    • You can get all the files using a few wildcards: 2007_12_10_04* 2007_12_10_05* 2007_12_10_06*, etc.
    • Important: Be sure to include _03* files if the track scan begins before the first line of the 04* file.
    • Note that for the inclinometer temperature data, each file covers several days, so you will need to identify which one (or possibly two) covers your data.
    • Note that on May 2008, the directory ServoMonitor-ServoMonitor-Az_El_50Hz was renamed to ServoMonitor-ServoMonitor-Az_El_FastRate, but the track map processor shall continue to use the old name. Just copy the files from the new location.
on target machine (that has Matlab), e.g. ssh ptcs@ptcs:
cd gbt/pointing/10_28_07/ExampleData
mkdir TPTCSPNT_071210_Track
cd TPTCSPNT_071210_Track
mkdir Accel1  Accel2  Accel3  AzEl50Hz  IncX1  IncX1Temp  IncX2  IncX2Temp  IncY1  IncY1Temp  IncY2  IncY2Temp

log into prospero and run:
cd /home/gbtlogs/Accelerometer-Accelerometer1-AccelerometerData
scp ./2007_12_10_04* ./2007_12_10_05* ./2007_12_10_06* ./2007_12_10_07* ./2007_12_10_08* ./2007_12_10_09* serenity:gbt/pointing/10_28_07/ExampleData/TPTCSPNT_071210_Track/Accel1
cd ../Accelerometer-Accelerometer2-AccelerometerData
scp ./2007_12_10_04* ./2007_12_10_05* ./2007_12_10_06* ./2007_12_10_07* ./2007_12_10_08* ./2007_12_10_09* serenity:gbt/pointing/10_28_07/ExampleData/TPTCSPNT_071210_Track/Accel2
cd ../Accelerometer-Accelerometer3-AccelerometerData
scp ./2007_12_10_04* ./2007_12_10_05* ./2007_12_10_06* ./2007_12_10_07* ./2007_12_10_08* ./2007_12_10_09* serenity:gbt/pointing/10_28_07/ExampleData/TPTCSPNT_071210_Track/Accel3
cd ../ServoMonitor-ServoMonitor-Az_El_50Hz
scp ./2007_12_10_04* ./2007_12_10_05* ./2007_12_10_06* ./2007_12_10_07* ./2007_12_10_08* ./2007_12_10_09* serenity:gbt/pointing/10_28_07/ExampleData/TPTCSPNT_071210_Track/AzEl50Hz
cd ../Inclinometer-Inclinometer-x1data/
scp ./2007_12_10_04* ./2007_12_10_05* ./2007_12_10_06* ./2007_12_10_07* ./2007_12_10_08* ./2007_12_10_09* serenity:gbt/pointing/10_28_07/ExampleData/TPTCSPNT_071210_Track/IncX1
cd ../Inclinometer-Inclinometer-x1temperature/
scp ./2007_12_09* serenity:gbt/pointing/10_28_07/ExampleData/TPTCSPNT_071210_Track/IncX1Temp
cd ../Inclinometer-Inclinometer-x2data/
scp ./2007_12_10_04* ./2007_12_10_05* ./2007_12_10_06* ./2007_12_10_07* ./2007_12_10_08* ./2007_12_10_09* serenity:gbt/pointing/10_28_07/ExampleData/TPTCSPNT_071210_Track/IncX2
cd ../Inclinometer-Inclinometer-x2temperature/
scp ./2007_12_09* serenity:gbt/pointing/10_28_07/ExampleData/TPTCSPNT_071210_Track/IncX2Temp
cd ../Inclinometer-Inclinometer-y1data/
scp ./2007_12_10_04* ./2007_12_10_05* ./2007_12_10_06* ./2007_12_10_07* ./2007_12_10_08* ./2007_12_10_09* serenity:gbt/pointing/10_28_07/ExampleData/TPTCSPNT_071210_Track/IncY1
cd ../Inclinometer-Inclinometer-y1temperature/
scp ./2007_12_09* serenity:gbt/pointing/10_28_07/ExampleData/TPTCSPNT_071210_Track/IncY1Temp
cd ../Inclinometer-Inclinometer-y2data/
scp ./2007_12_10_04* ./2007_12_10_05* ./2007_12_10_06* ./2007_12_10_07* ./2007_12_10_08* ./2007_12_10_09* serenity:gbt/pointing/10_28_07/ExampleData/TPTCSPNT_071210_Track/IncY2
cd ../Inclinometer-Inclinometer-y2temperature/
scp ./2007_12_09* serenity:gbt/pointing/10_28_07/ExampleData/TPTCSPNT_071210_Track/IncY2Temp

  • If the track was run CW, then CCW, use the file mainTrackMap2.m as a script template (located in ../10_28_07/ModelScript/trackMap)
  • If the track was run CCW, then CW, use the file mainTrackMap2_dec2007.m as a script template
  • Changes required to the template:
    • line 11: set dataDir
    • line 12: set startTime
    • line 13: set endTime
    • various lines: change 10dec2007 to the date of the new data (in order to get output files with the correct names).
    • After making these changes, save the script as: mainTrackMap2_new_date.m
  • Run matlab
    • Execute the script: mainTrackMap2_new_date.m
      • The elements of the data structure is as follows:
        1. time (seconds)
        2. Accel 1 X
        3. Accel 1 Y
        4. Accel 1 Z
        5. Accel 2 X
        6. Accel 2 Y
        7. Accel 2 Z
        8. Accel 3 X
        9. Accel 3 Y
        10. Accel 3 Z
        11. X1 inclinometer angle
        12. X2 inclinometer angle
        13. Y1 inclinometer angle
        14. Y2 inclinometer angle
        15. Az (50Hz, actually 10Hz)
        16. El (50Hz, actually 10Hz)
        17. AzRate (50Hz, actually 10Hz)
        18. ElRate (50Hz, actually 10Hz)
        19. AzAccel (50Hz, actually 10Hz)
        20. ElAccel (50Hz, actually 10Hz)
        21. X1 inclinometer temperature
        22. X2 inclinometer tempreature
        23. Y1 inclinometer temperature
        24. Y2 inclinometer temperature
    • If you have also recorded dwells at quantized positions, you can use median_15apr2008.m to compute median filtered results.
  • Here is another example for 04Sep2007
  • Here is another example for January 31, 2008
  • Here is another example for April 15, 2008
  • Here is another example for September 01, 2008


Revised, simplified version

This version eliminates the need to copy many log files to local directories.
  • Revised script to fetch the fits log files from /home/gbtlogs, rather than copying them all to a special directory.

First, log into the "ptcs" machine: "ssh ptcs@ptcs"

  • new version called "mainTrackMap3_28Mar2009.m"
  • in directory ~ptcs/pointing/10_28_07/ModelScript/trackMap/

  • modified routines :
    • fitsTemplate3.m --> fitsTemplate3_logs.m
    • injestFits.m --> injestFits3_logs.m

  • in "mainTrackMap3"
    • set dataDir='/home/gbtlogs';
    • set "startTime" and "endTime" to mjd of start and end of data
    • set namepat = 'yyyy_mm_dd' to find the FITS files for the correct day.

  • You probably still will have to fiddle with the azOffset and c0 and C360 to partition in to CW and CCW sections
    • this part still needs some work.


--
ToddHunter - 10 Dec 2007
Topic revision: r11 - 2010-06-28, FrankGhigo
This site is powered by FoswikiCopyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding NRAO Public Wiki? Send feedback