Double_t TPointStar::GetResidual

(

Double_t *

err = 0

)

const

Definition at line 55 of file TPointStar.cc.

References fRawAz, fRawEl, fStarAz, and fStarEl.

Referenced by TPointGui::Fit(), and GetMag().

{
    const Double_t del = fRawEl-fStarEl;
    const Double_t daz = fRawAz-fStarAz;

    const double x = cos(fRawEl) * cos(fStarEl) * cos(fStarAz-fRawAz);
    const double y = sin(fRawEl) * sin(fStarEl);

    const Double_t d = x + y;

    if (err)
    {
        // Error of one pixel in the CCD
        const Double_t e1 = 45./3600*TMath::DegToRad()  /4 * 0.5;

        // Error of the SE readout
        const Double_t e2 = 360./16384*TMath::DegToRad()/4 * 0.5;

        const Double_t e11 =  sin(del)+cos(fRawEl)*sin(fStarEl)*(1-cos(daz));
        const Double_t e12 =  cos(fRawEl)*cos(fStarEl)*sin(daz);

        const Double_t e21 = -sin(del)+sin(fRawEl)*cos(fStarEl)*(1-cos(daz));
        const Double_t e22 = -cos(fRawEl)*cos(fStarEl)*sin(daz);

        const Double_t err1  = sqrt(1-d*d);
        const Double_t err2  = (e11*e11 + e12*e12)*e1*e1;
        const Double_t err3  = (e21*e21 + e22*e22)*e2*e2;

        *err = sqrt(err2+err3)/err1 * TMath::RadToDeg();
    }

    const Double_t dist = acos(d);
    return dist * TMath::RadToDeg();
}

Here is the caller graph for this function: