Raw Visibilities

Following the steps described in the previous chapters we will calculate one raw visibility curve for the calibrator HD 10380 and one for the AGN of NGC 1068. Although the thresholds have been determined interactively I set it here with set_thresholds() to get reproducable results. Here it is more important to show the basic steps than getting optimized results. Thus improved visibilities are probably easy to find!

HD 10380

  • calfile=midigui(dir='<path>')
  • select files with index 1, 2, and 3
  • cal = xmdv(calfile)
  • confirm masks
  • cal -> gui,/true
  • cal -> set_thresholds,15523.9,0.00000,3.89843e-05

NGC 1068

  • scifile=midigui(dir='<path>')
  • select files with index 5, 6, and 7
  • sci = xmdv(scifile)
  • confirm masks
  • sci -> gui,/true
  • sci -> set_thresholds,3172.32,0.00000,4.00620e-05

Instrumental Visibility

In a next step we want to determine the instrumental visibility, i.e. the visibility of an unresolved point source, from a calibrator object of known diameter with the routine

  insvis = x->instruvisi(diameter,NAME=name)

where diameter is the diameter of the calibrator in mas. If the object is listed in the calibrator database, i.e. in 'heidelberg/', the diameter can be omitted. The name that is used to find the object in the database is either read from the MIDI data files (can be checked with object = cal -> get_objectname()) or can be given with the parameter NAME. The instrumental visibility is calculated by dividing the raw visibility with the expected visibility of the calibrator. When /VISPLOT is specified one gets a plot of the instrumental visibility, i.e. insvis[2,*] vs. insvis[0,*]. The parameter /PLOT launches the display for the wavelength-binned spectral power of the fringes.

HD 10380 is in the database. Thus the output of the command

  insvis = cal -> instruvisi()


  Looking up hd10380 in calibrator database
  Found:   2.85000 mas
  object:   hd10380
  Baseline:   U1 - U3
  ESO's baseline len.:   72.7130
  Date of observation:   2003-06-15T09:37:19.000
  Local Sid Time:   80914.915
  Right Ascension:   25.357838
  Declination:   5.4878200
  Our baseline length:   72.7233, p.a.:   7.63431
  Diameter 2.85000mas = 1.38172e-08rad

The last column is the expected visibility of the calibrator and thus the instrumental visibility is close to the raw visibility of HD 10380.

Calibrated Visibility

To calculate the calibrated visibility of the science object, i.e. to divide the raw visibiliy of the science object by the instrumental visibility, please use the routine

  calvis = sci->calibratedvisi(cal,diameter,NAME=name,/PLOT,/VISPLOT)

where diameter is the diameter of the calibrator in mas. Alternatively, the diameter can be read from the database ('heidelberg/') by using the name of the calibrator that is either set automatically from the fits header (can be checked with object = cal -> get_objectname()) or manually by using the paramter NAME. The advantage of the routine calibratedvisi() over calling instruvisi() manually and dividing the raw visibility of the science target by the result is the adjustment of the lambda-binning for the visibility of the calibrator, i.e. calibratedvisi() ensures that the same lambda-binning used for the science target is applied to the calibrator. The option /VISPLOT leads to a plot of the calibrated visibility. /PLOT launches the display for the wavelength-binned spectral power of the fringes for both the raw visibility of the science target and the now lambda-bins-adjusted instrumental visibility.

  calvis = sci->calibratedvisi(cal,/VISPLOT)

leads to

The wavelength and the calibrated visibility is stored in calvis[0,*] and calvis[2,*].