Ignore:
Timestamp:
08/26/04 14:13:36 (21 years ago)
Author:
wittek
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/mstarcam/MStarCamTrans.cc

    r4687 r4741  
    829829    //gLog << "mmtodeg = " << mmtodeg << endl;
    830830
     831    // aberr   is the ratio r_optaberr/r_ideal between
     832    //         the distance from the camera center with optical aberration and
     833    //         the distance from the camera center with an ideal imaging
     834    // the value 1.07 is valid if the expected position (with aberration)
     835    // in the camera is calculated as the average of the positions obtained
     836    // from the reflections at the individual mirrors
     837    Double_t aberr = 1.07;
     838
    831839    //--------------------------------------------------------------------
    832840
     
    937945        Loc0LocToCam(theta0, phi0, theta[j], phi[k],
    938946                           xx, yy);
    939         x[k] = xx * mmtodeg;
    940         y[k] = yy * mmtodeg;
     947        x[k] = xx * mmtodeg * aberr;
     948        y[k] = yy * mmtodeg * aberr;
    941949
    942950        //gLog << "theta0, phi0 = " << theta0 << ",  " << phi0
     
    10031011        Loc0LocToCam(theta0, phi0, theta1[j], phi1[k],
    10041012                           xx, yy);
    1005         x1[count] = xx * mmtodeg;
    1006         y1[count] = yy * mmtodeg;
     1013        x1[count] = xx * mmtodeg * aberr;
     1014        y1[count] = yy * mmtodeg * aberr;
    10071015
    10081016        //gLog << "theta0, phi0 = " << theta0 << ",  " << phi0
     
    11031111        Cel0CelToCam(dec0, hh0, dec[j], hx,
    11041112                           xx, yy);
    1105         xh[k] = xx * mmtodeg;
    1106         yh[k] = yy * mmtodeg;
     1113        xh[k] = xx * mmtodeg * aberr;
     1114        yh[k] = yy * mmtodeg * aberr;
    11071115
    11081116        //gLog << "dec0, h0 = " << dec0 << ",  " << h0
     
    11721180        Cel0CelToCam(dec0, hh0, dec1[j], hhx,
    11731181                           xx, yy);
    1174         xd[count] = xx * mmtodeg;
    1175         yd[count] = yy * mmtodeg;
     1182        xd[count] = xx * mmtodeg * aberr;
     1183        yd[count] = yy * mmtodeg * aberr;
    11761184
    11771185        //gLog << "dec0, h0 = " << dec0 << ",  " << h0
     
    12241232    return kTRUE;
    12251233}
     1234
Note: See TracChangeset for help on using the changeset viewer.