NasmythPointingModel

java.lang.Object
- stella.util.NasmythPointingModel

```
public class NasmythPointingModel
extends Object
```
Use this class to analyse differential pointing data obtained during a WiFSIP pointing campaign.

Author:

wooz

Nested Class Summary

Nested Classes
Modifier and Type	Class and Description
`static class`	`NasmythPointingModel.FileData` An access class to a pointing model that reads data from a file.
`private static class`	`NasmythPointingModel.Harmonics` A pure harmonic model
`private static class`	`NasmythPointingModel.ZonalHarmonicFunction` A pure harmonic function

Field Summary

Fields
Modifier and Type	Field and Description
`private Vector2D[]`	`azalt` The telescope nautic azimuth and altitude in radians.
`private static double`	`DEFLATITUDE` Default site, Teide observatory.
`private double`	`phi` The latitude of the site, degrees.
`private double[]`	`xaux` The measured x-coordinate on the auxiliary guider CCD.
`private double[]`	`yaux` The measured y-coordinate on the auxiliary guider CCD.

Constructor Summary

Constructors
Constructor and Description

NasmythPointingModel()

NasmythPointingModel(double latdeg)

Constructors
Constructor and Description
`NasmythPointingModel()`
`NasmythPointingModel(double latdeg)`

Method Summary

All Methods Static Methods Instance Methods Concrete Methods
Modifier and Type	Method and Description
`private static VectorG[]`	`convertAzAlt(Vector2D[] azalt)` We convert the argument vector2D array into to vectors, one with all azimuth data in radians, the other with cosine of zenith distance.
`DataModel`	`getAuxiliaryXModel()` We return the pure x-model for the guiding case.
`DataModel`	`getAuxiliaryYModel()` We return the pure y-model for the guiding case.
`static double`	`getAuxModelX(VectorG params, VectorG etael)` This method calculates the x-coordinate on the auxiliary CCD using the provided parameters with the following meaning, all in aux-ccd pixel, if not angles.
`static double`	`getAuxModelY(VectorG params, VectorG etael)` This method calculates the x-coordinate on the auxiliary CCD using the provided parameters with the following meaning, all in aux-ccd pixel, if not angles.
`Vector2D[]`	`getAzAlt()`
`static Vector2D[]`	`getAzCosz(Vector2D[] azalt)` Converts for the harmonic model the azimuth and height to azimuth (plain copy) and cosine of zenith distance.
`static Vector2D[]`	`getEtaElevation(Vector2D[] azalt, double phi)` From the measured nautical azimuth and elevation in radians, we calculate the internal coordinates, namely parallactic angle and elevation.
`GeneralLinearRegression`	`getHarmonicXFit(SphericalHarmonics.LM[] lm, double hlow, double zenith)` We return a pure harmonic model.
`DataModel`	`getHarmonicXModel(SphericalHarmonics.LM[] lm, double hlow, double zenith)` We return a pure harmonic model.
`GeneralLinearRegression`	`getHarmonicYFit(SphericalHarmonics.LM[] lm, double hlow, double zenith)` We return a pure harmonic model.
`DataModel`	`getHarmonicYModel(SphericalHarmonics.LM[] lm, double hlow, double zenith)` We return a pure harmonic model.
`static Vector2D`	`getOrigin(VectorG params, VectorG etael)` This method calculates the origin of the unperturbed Az/El system as seen from the science CCD.
`double[]`	`getXaux()`
`double[]`	`getYaux()`
`void`	`setAzAlt(Vector2D[] azalt)`
`void`	`setXaux(double[] xaux)`
`void`	`setYaux(double[] yaux)`
`VectorG`	`solve(DataModel dm, GeneralLinearRegression solution)`
`VectorG`	`solveY(DataModel dm)` We solve for the y model.

Methods inherited from class java.lang.Object
clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

- Field Detail
  - DEFLATITUDE
```
private static final double DEFLATITUDE
```
    Default site, Teide observatory.
  - phi
```
private double phi
```
    The latitude of the site, degrees.
  - azalt
```
private Vector2D[] azalt
```
    The telescope nautic azimuth and altitude in radians.
  - xaux
```
private double[] xaux
```
    The measured x-coordinate on the auxiliary guider CCD.
  - yaux
```
private double[] yaux
```
    The measured y-coordinate on the auxiliary guider CCD.
- Constructor Detail
  - NasmythPointingModel
```
public NasmythPointingModel()
```
  - NasmythPointingModel
```
public NasmythPointingModel(double latdeg)
```
- Method Detail
  - getAzAlt
```
public Vector2D[] getAzAlt()
```
    Returns:
    
    the azalt
  - setAzAlt
```
public void setAzAlt(Vector2D[] azalt)
```
    Parameters:
    
    azalt - the azalt to set
  - getXaux
```
public double[] getXaux()
```
    Returns:
    
    the xaux
  - setXaux
```
public void setXaux(double[] xaux)
```
    Parameters:
    
    xaux - the xaux to set
  - getYaux
```
public double[] getYaux()
```
    Returns:
    
    the yaux
  - setYaux
```
public void setYaux(double[] yaux)
```
    Parameters:
    
    yaux - the yaux to set
  - getAuxModelX
```
public static final double getAuxModelX(VectorG params,
                                        VectorG etael)
```
    This method calculates the x-coordinate on the auxiliary CCD using the provided parameters with the following meaning, all in aux-ccd pixel, if not angles. Angles in radians.
    - #0: Distance of M3 point of origin seen from science CCD
    - #1: Angle of M3 POO seen from science CCD at derot=0.
    - #2: Distance of M2-M0 point of origin seen from science CCD.
    - #3: Angle of M2-M0 POO from science CCD at E=0 minus angle at index 1.
    - #4: Constant offset
    The argument is already in parallactic angle and elevation, both radians
  - getAuxModelY
```
public static final double getAuxModelY(VectorG params,
                                        VectorG etael)
```
    This method calculates the x-coordinate on the auxiliary CCD using the provided parameters with the following meaning, all in aux-ccd pixel, if not angles. Angles in radians.
    - #0: Distance of M3 point of origin seen from science CCD
    - #1: Angle of M3 POO seen from science CCD at derot=0.
    - #2: Distance of M2-M0 point of origin seen from science CCD.
    - #3: Angle of M2-M0 POO from science CCD at E=0 minus angle at index 1.
    - #4: Constant offset
    - #5: Differential tube flexure.
    The argument is already in parallactic angle and elevation, both radians
  - getOrigin
```
public static final Vector2D getOrigin(VectorG params,
                                       VectorG etael)
```
    This method calculates the origin of the unperturbed Az/El system as seen from the science CCD. The parameters used in the argument Vector are the first four parameters, others are ignored.
    - #0: Distance of M3 point of origin seen from science CCD
    - #1: Angle of M3 POO seen from science CCD at derot=0.
    - #2: Distance of M2-M0 point of origin seen from science CCD.
    - #3: Angle of M2-M0 POO from science CCD at E=0 minus angle at index 1.
    The argument is already in parallactic angle and elevation, both radians
  - getEtaElevation
```
public static final Vector2D[] getEtaElevation(Vector2D[] azalt,
                                               double phi)
```
    From the measured nautical azimuth and elevation in radians, we calculate the internal coordinates, namely parallactic angle and elevation. Formular for the parallactica angle η can be written as
```
 tan(η) = -sin(az)/(cos(E)tan(φ)-sin(E)*cos(az))
 
```
    with φ the geographic latitude of the observer.
    Parameters:
    
    azalt - The nautical azimuth and elevation of the measures, radians
    
    phi - Latitude of observer, radians.
  - getAzCosz
```
public static final Vector2D[] getAzCosz(Vector2D[] azalt)
```
    Converts for the harmonic model the azimuth and height to azimuth (plain copy) and cosine of zenith distance.
  - getAuxiliaryXModel
```
public DataModel getAuxiliaryXModel()
```
    We return the pure x-model for the guiding case.
    
    Returns:
    
    An independent model for the auxiliary x coordinate
  - getAuxiliaryYModel
```
public DataModel getAuxiliaryYModel()
```
    We return the pure y-model for the guiding case.
    
    Returns:
    
    An independent model for the auxiliary y coordinate
  - getHarmonicXModel
```
public DataModel getHarmonicXModel(SphericalHarmonics.LM[] lm,
                                   double hlow,
                                   double zenith)
```
    We return a pure harmonic model. The model fit is actually a residual fit to the RMS to a constant (the average in x)
  - getHarmonicYModel
```
public DataModel getHarmonicYModel(SphericalHarmonics.LM[] lm,
                                   double hlow,
                                   double zenith)
```
    We return a pure harmonic model. The model fit is actually a residual fit to y0+y1*cos(h).
  - getHarmonicXFit
```
public GeneralLinearRegression getHarmonicXFit(SphericalHarmonics.LM[] lm,
                                               double hlow,
                                               double zenith)
```
    We return a pure harmonic model. The model fit is actually a residual fit to y0+y1*cos(h).
  - getHarmonicYFit
```
public GeneralLinearRegression getHarmonicYFit(SphericalHarmonics.LM[] lm,
                                               double hlow,
                                               double zenith)
```
    We return a pure harmonic model. The model fit is actually a residual fit to y0+y1*cos(h).
  - convertAzAlt
```
private static VectorG[] convertAzAlt(Vector2D[] azalt)
```
    We convert the argument vector2D array into to vectors, one with all azimuth data in radians, the other with cosine of zenith distance.
    
    Parameters:
    
    azalt -
    
    Returns:
    
    A two-dim VectorG array, at first index azimuth, cos(z) on second
  - solve
```
public VectorG solve(DataModel dm,
                     GeneralLinearRegression solution)
```
  - solveY
```
public VectorG solveY(DataModel dm)
```
    We solve for the y model. The starting values for the Amoeba are estimated and hard-coded.

Class NasmythPointingModel

Nested Class Summary

Field Summary

Constructor Summary

Method Summary

Methods inherited from class java.lang.Object

Field Detail

DEFLATITUDE

phi

azalt

xaux

yaux

Constructor Detail

NasmythPointingModel

NasmythPointingModel

Method Detail

getAzAlt

setAzAlt

getXaux

setXaux

getYaux

setYaux

getAuxModelX

getAuxModelY

getOrigin

getEtaElevation

getAzCosz

getAuxiliaryXModel

getAuxiliaryYModel

getHarmonicXModel

getHarmonicYModel

getHarmonicXFit

getHarmonicYFit

convertAzAlt

solve

solveY