Changeset 8324 for trunk/MagicSoft/Mars/mastro

Timestamp:

02/16/07 09:21:13 (18 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mastro

Files:

• MAstro.cc (modified) (1 diff)
• MAstro.h (modified) (2 diffs)
• MObservatory.cc (modified) (1 diff)

trunk/MagicSoft/Mars/mastro/MAstro.cc

@@ -446 +446 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Estimates the time at which a source culminates.
+//
+// ra: right ascension [rad]
+// elong: observers longitude [rad]
+// mjd: modified julian date (utc)
+//
+// return time in [-12;12]
+//
+Double_t MAstro::EstimateCulminationTime(Double_t mjd, Double_t elong, Double_t ra)
+{
+    // startime at 1.1.2000 for greenwich 0h
+    const Double_t gmt0 = 6.664520;
+
+    // difference of startime for greenwich for two calendar days [h]
+    const Double_t d0 = 0.06570982224;
+
+    // mjd of greenwich 1.1.2000 0h
+    const Double_t mjd0 = 51544;
+
+    // mjd today
+    const Double_t mjd1 = TMath::Floor(mjd);
+
+    // scale between star-time and sun-time
+    const Double_t scale = 1;//1.00273790926;
+
+    const Double_t UT = (ra-elong)*RadToHor() - (gmt0 + d0 * (mjd1-mjd0))/scale;
+
+    return fmod(2412 + UT, 24) - 12;
+}
 
 // --------------------------------------------------------------------------

trunk/MagicSoft/Mars/mastro/MAstro.h

@@ -31 +31 @@
     static Double_t HorToDeg()        { return 1./15; }
 
-    static Double_t RadToHor()        { return 15/TMath::TwoPi(); }
-    static Double_t HorToRad()        { return TMath::TwoPi()/15; }
+    static Double_t RadToHor()        { return 24/TMath::TwoPi(); }
+    static Double_t HorToRad()        { return TMath::TwoPi()/24; }
 
     // Angle/Time conversion functions
@@ -86 +86 @@
     static Int_t    GetMagicPeriod(Double_t mjd);
 
+    static Double_t EstimateCulminationTime(Double_t mjd, Double_t longit, Double_t ra);
+
     // Estimate some parameters around the sun
     static Double_t GetSunRaDec(Double_t mjd, Double_t &ra, Double_t &dec);

trunk/MagicSoft/Mars/mastro/MObservatory.cc

@@ -110 +110 @@
     *fLog << all;
     *fLog << underline << fObservatoryName << ":" << endl;
-    *fLog << "Latitude:  " << TMath::Abs(fLatitude*kRad2Deg)  << " deg " << (fLatitude  > 0 ? "W" : "E") << endl;
-    *fLog << "Longitude: " << TMath::Abs(fLongitude*kRad2Deg) << " deg " << (fLongitude < 0 ? "N" : "S") << endl;
+    *fLog << "Latitude:  " << TMath::Abs(fLatitude*kRad2Deg)  << " deg " << (fLatitude  > 0 ? "N" : "S") << endl;
+    *fLog << "Longitude: " << TMath::Abs(fLongitude*kRad2Deg) << " deg " << (fLongitude < 0 ? "E" : "W") << endl;
     *fLog << "Height:    " << fHeight << "m" << endl;
 }