Changeset 8915 for trunk


Ignore:
Timestamp:
06/02/08 16:23:28 (16 years ago)
Author:
tbretz
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars/mpointing
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • trunk/MagicSoft/Mars/mpointing/MPointing.cc

    r8829 r8915  
    101101#define DEBUG(txt)
    102102
     103void ZdAz::Round()
     104{
     105    fX = TMath::Nint(fX);
     106    fY = TMath::Nint(fY);
     107}
     108
     109void ZdAz::Abs()
     110{
     111    fX = TMath::Abs(fX);
     112    fY = TMath::Abs(fY);
     113}
     114
    103115void MPointing::Init(const char *name, const char *title)
    104116{
    … …  
    357369}
    358370
     371ZdAz MPointing::AddOffsets(const ZdAz &zdaz) const
     372{
     373    AltAz p(TMath::Pi()/2-zdaz.Zd(), zdaz.Az());
     374
     375    AltAz c = AddOffsets(p);
     376
     377    return ZdAz(TMath::Pi()/2-c.Alt(), c.Az());
     378}
     379
     380TVector3 MPointing::AddOffsets(const TVector3 &v) const
     381{
     382    AltAz p(TMath::Pi()/2-v.Theta(), v.Phi());
     383    AltAz c = AddOffsets(p);
     384
     385    TVector3 rc;
     386    rc.SetMagThetaPhi(1, TMath::Pi()/2-c.Alt(), c.Az());
     387    return rc;
     388}
     389
    359390AltAz MPointing::SubtractOffsets(const AltAz &aa) const
    360391{
    … …  
    367398
    368399    return p;
     400}
     401
     402ZdAz MPointing::SubtractOffsets(const ZdAz &zdaz) const
     403{
     404    AltAz p(TMath::Pi()/2-zdaz.Zd(), zdaz.Az());
     405
     406    AltAz c = SubtractOffsets(p);
     407
     408    return ZdAz(TMath::Pi()/2-c.Alt(), c.Az());
     409}
     410
     411TVector3 MPointing::SubtractOffsets(const TVector3 &v) const
     412{
     413    AltAz p(TMath::Pi()/2-v.Theta(), v.Phi());
     414    AltAz c = SubtractOffsets(p);
     415
     416    TVector3 rc;
     417    rc.SetMagThetaPhi(1, TMath::Pi()/2-c.Alt(), c.Az());
     418    return rc;
    369419}
    370420
    … …  
    707757}
    708758
     759void MPointing::GetError(TArrayD &par) const
     760{
     761    par = fError;
     762    for (int i=0; i<kNumPar; i++)
     763        par[i] *= TMath::RadToDeg();
     764}
     765
     766TVector2 MPointing::GetDxy() const
     767{
     768    return TVector2(fDx, fDy)*TMath::RadToDeg();
     769}
     770
     771Double_t MPointing::GetPx() const
     772{
     773    return fPx*TMath::RadToDeg();
     774}
     775
     776Double_t MPointing::GetPy() const
     777{
     778    return fPy*TMath::RadToDeg();
     779}
     780
    709781void MPointing::SetMinuitParameters(TMinuit &m, Int_t n) const
    710782{

  • trunk/MagicSoft/Mars/mpointing/MPointing.h

    r8877 r8915  
    9191    // MSlewing only?!?
    9292    double Ratio() const { return fX/fY; }
    93     void Round()         { fX=(int)(floor(fX+.5)); fY=(int)(floor(fY+.5)); }
    94     void Abs()           { fX=fabs(fX); fY=fabs(fY); }
     93    void Round();
     94    void Abs();
    9595
    9696    ClassDef(ZdAz, 0)
    … …  
    250250    }
    251251
    252     AltAz  AddOffsets(const AltAz &aa) const;
    253     ZdAz   AddOffsets(const ZdAz &zdaz) const
    254     {
    255         AltAz p(TMath::Pi()/2-zdaz.Zd(), zdaz.Az());
    256         AltAz c = AddOffsets(p);
    257         return ZdAz(TMath::Pi()/2-c.Alt(), c.Az());
    258     }
    259     TVector3 AddOffsets(const TVector3 &v) const
    260     {
    261         AltAz p(TMath::Pi()/2-v.Theta(), v.Phi());
    262         AltAz c = AddOffsets(p);
    263         TVector3 rc;
    264         rc.SetMagThetaPhi(1, TMath::Pi()/2-c.Alt(), c.Az());
    265         return rc;
    266     }
    267 
    268     AltAz  SubtractOffsets(const AltAz &aa) const;
    269     ZdAz   SubtractOffsets(const ZdAz &zdaz) const
    270     {
    271         AltAz p(TMath::Pi()/2-zdaz.Zd(), zdaz.Az());
    272         AltAz c = SubtractOffsets(p);
    273         return ZdAz(TMath::Pi()/2-c.Alt(), c.Az());
    274     }
    275     TVector3 SubtractOffsets(const TVector3 &v) const
    276     {
    277         AltAz p(TMath::Pi()/2-v.Theta(), v.Phi());
    278         AltAz c = SubtractOffsets(p);
    279         TVector3 rc;
    280         rc.SetMagThetaPhi(1, TMath::Pi()/2-c.Alt(), c.Az());
    281         return rc;
    282     }
     252    AltAz    AddOffsets(const AltAz &aa) const;
     253    ZdAz     AddOffsets(const ZdAz &zdaz) const;
     254    TVector3 AddOffsets(const TVector3 &v) const;
     255
     256    AltAz    SubtractOffsets(const AltAz &aa) const;
     257    ZdAz     SubtractOffsets(const ZdAz &zdaz) const;
     258    TVector3 SubtractOffsets(const TVector3 &v) const;
    283259
    284260    void SetParameters(const Double_t *par, Int_t n=kNumPar);
    … …  
    294270        GetParameters(par.GetArray());
    295271    }
    296     void GetError(TArrayD &par) const
    297     {
    298         par = fError;
    299         for (int i=0; i<kNumPar; i++)
    300             par[i] *= TMath::RadToDeg();
    301     }
     272    void GetError(TArrayD &par) const;
    302273
    303274    Double_t &operator[](UInt_t i) { return *fCoeff[i]; }
    … …  
    327298     */
    328299
    329     TVector2 GetDxy() const { return TVector2(fDx, fDy)*TMath::RadToDeg(); }
    330 
    331     Double_t GetPx() const { return fPx*TMath::RadToDeg(); }
    332     Double_t GetPy() const { return fPy*TMath::RadToDeg(); }
     300    TVector2 GetDxy() const;// { return TVector2(fDx, fDy)*TMath::RadToDeg(); }
     301
     302    Double_t GetPx() const;// { return fPx*TMath::RadToDeg(); }
     303    Double_t GetPy() const;// { return fPy*TMath::RadToDeg(); }
    333304
    334305    Bool_t IsPxValid() const { return fError[kPX]>0; }

  • trunk/MagicSoft/Mars/mpointing/MPointingPos.cc

    r7887 r8915  
    5757const TString MPointingPos::gsDefName  = "MPointingPos";
    5858const TString MPointingPos::gsDefTitle = "Container storing the (corrected) telescope pointing position";
     59
     60Double_t MPointingPos::GetZdRad() const
     61{
     62    return fZd*TMath::DegToRad();
     63}
     64
     65Double_t MPointingPos::GetAzRad() const
     66{
     67    return fAz*TMath::DegToRad();
     68}
     69
     70Double_t MPointingPos::GetRaRad() const
     71{
     72    return fRa*TMath::DegToRad()*15;
     73}
     74
     75Double_t MPointingPos::GetDecRad() const
     76{
     77    return fDec*TMath::DegToRad();
     78}
    5979
    6080// --------------------------------------------------------------------------

  • trunk/MagicSoft/Mars/mpointing/MPointingPos.h

    r7550 r8915  
    5151    Double_t GetAz() const  { return fAz; }
    5252
    53     Double_t GetZdRad() const  { return fZd*TMath::DegToRad(); }
    54     Double_t GetAzRad() const  { return fAz*TMath::DegToRad(); }
     53    Double_t GetZdRad() const;//  { return fZd*TMath::DegToRad(); }
     54    Double_t GetAzRad() const;//  { return fAz*TMath::DegToRad(); }
    5555
    5656    ZdAz GetZdAz() const { return ZdAz(fZd, fAz); }
    … …  
    6060    Double_t GetDec() const { return fDec; }
    6161
    62     Double_t GetRaRad() const  { return fRa*TMath::DegToRad()*15; }
    63     Double_t GetDecRad() const { return fDec*TMath::DegToRad(); }
     62    Double_t GetRaRad() const;//  { return fRa*TMath::DegToRad()*15; }
     63    Double_t GetDecRad() const;// { return fDec*TMath::DegToRad(); }
    6464
    6565    Double_t RotationAngle(const MObservatory &o) const;
Note: See TracChangeset for help on using the changeset viewer.