Changeset 7202 for trunk/MagicSoft/Mars/mastro

Timestamp:

07/21/05 11:44:05 (19 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mastro

Files:

• MAstro.cc (modified) (1 diff)
• MAstro.h (modified) (1 diff)

trunk/MagicSoft/Mars/mastro/MAstro.cc

@@ -601 +601 @@
     return v1.XYvector(); //TVector2(v1.Y(), -v1.X());//v1.XYvector();
 }
+
+// --------------------------------------------------------------------------
+//
+// Calculate the absolute misspointing from the nominal zenith angle nomzd
+// and the deviations in zd (devzd) and az (devaz).
+// All values given in deg, the return value, too.
+//
+Double_t MAstro::GetDevAbs(Double_t nomzd, Double_t devzd, Double_t devaz)
+{
+    const Double_t pzd = nomzd * TMath::DegToRad();
+    const Double_t azd = devzd * TMath::DegToRad();
+    const Double_t aaz = devaz * TMath::DegToRad();
+
+    const double el = TMath::Pi()/2-pzd;
+
+    const double dphi2 = aaz/2.;
+    const double cos2  = TMath::Cos(dphi2)*TMath::Cos(dphi2);
+    const double sin2  = TMath::Sin(dphi2)*TMath::Sin(dphi2);
+    const double d     = TMath::Cos(azd)*cos2 - TMath::Cos(2*el)*sin2;
+
+    return TMath::ACos(d)*TMath::RadToDeg();
+}

trunk/MagicSoft/Mars/mastro/MAstro.h

@@ -70 +70 @@
     static Int_t    GetMagicPeriod(Double_t mjd);
 
+    // Get distance between v1 and v0 in a plain perpendicular to v0 in distance dist
     static TVector2 GetDistOnPlain(const TVector3 &v0, TVector3 v1, Double_t dist=1);
+
+    // Calculate absolute misspointing from error in zd and error in az
+    static Double_t GetDevAbs(Double_t nomzd, Double_t devzd, Double_t devaz);
 
     ClassDef(MAstro, 0)