Merge pull request #116 from lsst/tickets/DM-8052

Load LST into VisitInfo to approximate ERA tickets/DM-8052
lsst · Nov 18, 2016 · 22f008f · 22f008f
2 parents f354919 + 64cdf2b
commit 22f008f
Show file tree

Hide file tree

Showing 3 changed files with 225 additions and 270 deletions.
diff --git a/include/lsst/afw/image/VisitInfo.h b/include/lsst/afw/image/VisitInfo.h
@@ -177,6 +177,12 @@ class VisitInfo : public table::io::PersistableFacade<VisitInfo>, public table::
 
     bool isPersistable() const { return true; }
 
+    // get the local sidereal time on the meridian (equivalent, but not equal, to Local Mean Sidereal Time)
+    geom::Angle getLocalEra() const;
+
+    // get hour angle at the boresight
+    geom::Angle getBoresightHourAngle() const;
+
 protected:
 
     virtual std::string getPersistenceName() const;

diff --git a/src/image/VisitInfo.cc b/src/image/VisitInfo.cc
@@ -447,5 +447,12 @@ void VisitInfo::write(OutputArchiveHandle & handle) const {
     handle.saveCatalog(cat);
 }
 
+geom::Angle VisitInfo::getLocalEra() const {
+    return getEra() + getObservatory().getLongitude();
+}
+
+geom::Angle VisitInfo::getBoresightHourAngle() const {
+    return getLocalEra() - getBoresightRaDec()[0];
+}
 
 }}}  // namespace