Write code to calculate eta_mb as a function of time

aakepley commented 7 years ago

I've put together some code to do the flux calibration at /home/scratch/akepley/gbt_hcn/degas/dcr_test/peak_fit.py. You'll also need caldb.py in same directory.

I've run this all the way through and it works and delivers appropriate values for the eta_a and eta_mb. Some things that might be improved:

Right now I only fit the first az scan. Fitting both az scans and both el scans separately and using the average of each for the eta_mb calculation would improve the calculation of the beam size. This would probably require some error checking if one of the scans is bad.
We're using a very beta web data base query for the flux calibration. After looking at the available options, this seemed the easiest to do rather than go through NAASC-only CASA routines. However, there may be bugs and bad values, so good to have a check and not use blindly.
We also don't have a good flux cal for the first session due to a vane issue. Might be good to have a way to put eta_a and eta_mb in manually.

aakepley commented 6 years ago

Amanda will flesh out strategy. Jialu (@JialLi) has done some work. Also need to involve David Frayer.

low-sky commented 6 years ago

Run this as a PR model

aakepley commented 5 years ago

Based on talk with Erik: I need to take the logs, parse them, figure out what scans are flux calibration scans, and then fit a Gaussian to them to get eta_mb. I have the Gaussian fitting code in fit_peak. Need to add the driver code and clean things up a bit.

aakepley commented 5 years ago

Need to compare with David's manual fits.

GBTSpectroscopy / degas

Write code to calculate eta_mb as a function of time #14