Changeset 18629


Ignore:
Timestamp:
09/18/16 14:55:45 (8 years ago)
Author:
tbretz
Message:
Disentengle from root and Mars were appropriate, some tidying.
Location:
trunk/FACT++/drive
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/FACT++/drive/TPointStar.cc

    r18618 r18629  
    11#include "TPointStar.h"
    22
     3#include <iostream>
    34#include <TMath.h>
    45
    5 #include "MLog.h"
    6 #include "MLogManip.h"
    7 
    86#include "MPointing.h"
    9 
    10 ClassImp(TPointStar);
    117
    128using namespace std;
     
    4339AltAz TPointStar::GetRawAltAz() const   { return AltAz(fRawEl, fRawAz); }
    4440
    45 void TPointStar::AdjustEl(Double_t del) { fStarEl += del*TMath::DegToRad(); }
    46 void TPointStar::AdjustAz(Double_t daz) { fStarAz += daz*TMath::DegToRad(); }
    47 
    4841void TPointStar::Adjust(const MPointing &bend)
    4942{
     
    7366    {
    7467        // Error of one pixel in the CCD
    75         const Double_t e1 = 32./3600*TMath::DegToRad()  * 0.5;
     68        const Double_t e1 = 45./3600*TMath::DegToRad()  /4 * 0.5;
    7669
    77         // Error of one SE unit
    78         const Double_t e2 = 360./16384*TMath::DegToRad() * 0.5;
     70        // Error of the SE readout
     71        const Double_t e2 = 360./16384*TMath::DegToRad()/4 * 0.5;
    7972
    8073        const Double_t e11 =  sin(del)+cos(fRawEl)*sin(fStarEl)*(1-cos(daz));
     
    109102    if (n<4)
    110103    {
    111         gLog << err << "Read: ERROR - Not enough numbers" << endl;
     104        cout << "Read: ERROR - Not enough numbers" << endl;
    112105        return fin;
    113106    }
     
    120113    set.fRawEl  = v[3]*TMath::DegToRad();
    121114
     115
     116
    122117    if (fin)
    123118    {
    124119        Double_t res, err;
    125120        res = set.GetResidual(&err);
    126         gLog << inf << "Read: " << v[0] << " " << v[1] << "  :  " << v[2] << " " << v[3] << "  :  " << v[2]-v[0] << " " << v[3]-v[1] << "  :  " << res << " " << err << " " << err/res << endl;
     121        cout << "Read: " << v[0] << " " << v[1] << "  :  " << v[2] << " " << v[3] << "  :  " << v[2]-v[0] << " " << v[3]-v[1] << "  :  " << res << " " << err << " " << err/res << endl;
    127122    }
    128123
  • trunk/FACT++/drive/TPointStar.h

    r18618 r18629  
    55#include <TNamed.h>
    66#endif
    7 
    8 //class istream;
    9 //class ostream;
    107
    118class ZdAz;
     
    4239    Double_t GetMag() const { return fMag; }
    4340    Double_t GetResidual(Double_t *err=0) const;
    44     //Double_t GetResidual() const;
    4541
    4642    Double_t GetDEl() const;//     { return (fRawEl-fStarEl)*TMath::RadToDeg(); }
     
    6056    AltAz GetRawAltAz() const;//   { return AltAz(fRawEl, fRawAz); }
    6157
    62     void AdjustEl(Double_t del);// { fStarEl += del*TMath::DegToRad(); }
    63     void AdjustAz(Double_t daz);// { fStarAz += daz*TMath::DegToRad(); }
    64 
    65     void Adjust(const MPointing &bend);/*
    66     {
    67         AltAz p = bend(GetStarAltAz());
    68         fStarEl = p.Alt();
    69         fStarAz = p.Az();
    70     }*/
    71     void AdjustBack(const MPointing &bend);/*
    72     {
    73         AltAz p = bend.CorrectBack(GetRawAltAz());
    74         fRawEl = p.Alt();
    75         fRawAz = p.Az();
    76     }*/
    77     ClassDef(TPointStar, 1)
     58    void Adjust(const MPointing &bend);
     59    void AdjustBack(const MPointing &bend);
    7860};
    79 
    80 /*
    81 class TPoint : public TPointStar
    82 {
    83 public:
    84     TPoint(Double_t sel=0, Double_t saz=0, Double_t rel=0, Double_t raz=0)
    85         : TPointStar(sel, saz, rel, raz)
    86     {
    87     }
    88     ClassDef(TPoint, 1)
    89 };
    90 */
    9161
    9262std::istream &operator>>(std::istream &fin, TPointStar &set);
Note: See TracChangeset for help on using the changeset viewer.