fit_planet_cal - Fit absolute calibration model

This script lives in analysis.det_sens. Its purpose is to analyze planet peak height results and produce a calibration model for a particular data set (e.g. a single season of a single array and band). This code operates on a list of planet peak heights, along with a solid angle estimate, and does not process planet maps directly.

Running the fitter

After setting up the configuration file, run:

moby2 fit_planet_cal fit_planet_cal.in

This will call the driver program in main().

The input file will need modification to specify appropriate data selection and cuts.

Configuration options

The configuration file for fit_planet_cal is fairly dense. An example can be found here: fit_planet_cal.in.

You might want to refer to the docs/source for load_amplitudes()

Auto-doc

moby2.analysis.det_sens.fit_planet_cal.main(args_in)

moby2.analysis.det_sens.fit_planet_cal.load_amplitudes(params, info_dict)

Read planet peak heights from some kind of file. The info_dict will be used to format the filename according to what data subset is being analyzed. The params argument is a dict that looks like one of these:

{'type': 'quick_pick',
 'filename': 'quick_beam.pik',
 'field': 'gaussian_amp',
 'masks': ...,      #optional
 'also_load': ...,  #optional
}

{'type': 'column_file',
 'filename': 'peaks_{fcode}.txt',
 'columns': [0,1],  #optional
}

{'type': 'solid_angle',
 'filename': 'solid_angle/table.fits'
}