- Timestamp:
- 09/18/16 14:55:45 (8 years ago)
- Location:
- trunk/FACT++/drive
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/FACT++/drive/TPointStar.cc
r18618 r18629 1 1 #include "TPointStar.h" 2 2 3 #include <iostream> 3 4 #include <TMath.h> 4 5 5 #include "MLog.h"6 #include "MLogManip.h"7 8 6 #include "MPointing.h" 9 10 ClassImp(TPointStar);11 7 12 8 using namespace std; … … 43 39 AltAz TPointStar::GetRawAltAz() const { return AltAz(fRawEl, fRawAz); } 44 40 45 void TPointStar::AdjustEl(Double_t del) { fStarEl += del*TMath::DegToRad(); }46 void TPointStar::AdjustAz(Double_t daz) { fStarAz += daz*TMath::DegToRad(); }47 48 41 void TPointStar::Adjust(const MPointing &bend) 49 42 { … … 73 66 { 74 67 // 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; 76 69 77 // Error of one SE unit78 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; 79 72 80 73 const Double_t e11 = sin(del)+cos(fRawEl)*sin(fStarEl)*(1-cos(daz)); … … 109 102 if (n<4) 110 103 { 111 gLog << err<< "Read: ERROR - Not enough numbers" << endl;104 cout << "Read: ERROR - Not enough numbers" << endl; 112 105 return fin; 113 106 } … … 120 113 set.fRawEl = v[3]*TMath::DegToRad(); 121 114 115 116 122 117 if (fin) 123 118 { 124 119 Double_t res, err; 125 120 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; 127 122 } 128 123 -
trunk/FACT++/drive/TPointStar.h
r18618 r18629 5 5 #include <TNamed.h> 6 6 #endif 7 8 //class istream;9 //class ostream;10 7 11 8 class ZdAz; … … 42 39 Double_t GetMag() const { return fMag; } 43 40 Double_t GetResidual(Double_t *err=0) const; 44 //Double_t GetResidual() const;45 41 46 42 Double_t GetDEl() const;// { return (fRawEl-fStarEl)*TMath::RadToDeg(); } … … 60 56 AltAz GetRawAltAz() const;// { return AltAz(fRawEl, fRawAz); } 61 57 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); 78 60 }; 79 80 /*81 class TPoint : public TPointStar82 {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 */91 61 92 62 std::istream &operator>>(std::istream &fin, TPointStar &set);
Note:
See TracChangeset
for help on using the changeset viewer.