Changeset 2761 for trunk/external-alma/atnf/PKSIO/NROReader.h

Timestamp:

02/06/13 13:13:20 (11 years ago)

Author:

Takeshi Nakazato

Message:

New Development: No

JIRA Issue: No

Ready for Test: Yes

Interface Changes: No

What Interface Changed: Please list interface changes

Test Programs: List test programs

Put in Release Notes: Yes/No?

Module(s): Module Names change impacts.

Description: Describe your changes here...

Added new option, freqref, to NRO filler. Posible values are:
1) 'rest' to import frequency in REST frame, which results in an exactly
same frequency label as NEWSTAR, and 2) 'vref' to import frequency
in the frame that source velocity refers, which results in the same
velocity label as NEWSTAR. The option must be given to scantable
constructor.

File:

• trunk/external-alma/atnf/PKSIO/NROReader.h (modified) (10 diffs)

trunk/external-alma/atnf/PKSIO/NROReader.h

@@ -42 +42 @@
 #include <casa/BasicSL/String.h>
 #include <measures/Measures/MPosition.h>
+#include <measures/Measures/MEpoch.h>
 #include <measures/Measures/MCPosition.h>
 #include <measures/Measures/MDirection.h>
@@ -56 +57 @@
 #include <atnf/PKSIO/NRODataset.h>
 #include <atnf/PKSIO/NRODataRecord.h>
-
-using namespace std ;
 
 // <summary>
@@ -103 +102 @@
   // Destructor.
   virtual ~NROReader() ;
+
+  // determine whether to import frequency as REST (frequency is same as 
+  // NEWSTAR/NOSTAR) or as is (velocity is same as NEWSTAR/NOSTAR)
+  void setFreqRefFromVREF( bool fromVREF ) ;
 
   // Read data header
@@ -133 +136 @@
                            uInt &beamno,
                            uInt &polno,
-                           vector<double> &freqs,   
+                           std::vector<double> &freqs,   
                            Vector<Double> &restfreq,  
                            uInt &refbeamno,
@@ -171 +174 @@
 
   // Get IF settings
-  virtual vector<Bool> getIFs() ;
+  virtual std::vector<Bool> getIFs() ;
 
   // Get Number of IFs
@@ -177 +180 @@
 
   // Get Beam settings
-  virtual vector<Bool> getBeams() ;
+  virtual std::vector<Bool> getBeams() ;
 
   // Get Number of Beams
@@ -188 +191 @@
 
   // Get spectrum
-  virtual vector< vector<double> > getSpectrum() ;
+  virtual std::vector< std::vector<double> > getSpectrum() ;
 
   // Get number of polarization
@@ -196 +199 @@
   virtual double getStartTime() ;
   virtual double getEndTime() ;
-  virtual vector<double> getStartIntTime() ;
+  virtual std::vector<double> getStartIntTime() ;
   //virtual double getStartIntTime( int i ) ;
 
   // Get Antenna Position in ITRF coordinate
-  virtual vector<double> getAntennaPosition() = 0 ;
+  virtual std::vector<double> getAntennaPosition() = 0 ;
 
   // Get SRCDIRECTION in RADEC(J2000)
@@ -209 +212 @@
   virtual void initConvert( int icoord, double t, char *epoch ) ;
 
+  // Shift frequency by given velocity with respect to specified 
+  // velocity reference
+  std::vector<double> shiftFrequency( const std::vector<double> &f, const double &v, const string &vref ) ; 
+
   // filename 
   string filename_ ;
@@ -222 +229 @@
   CountedPtr<MDirection::Convert> converter_ ;
   CountedPtr<MeasFrame> mf_ ;
+  MEpoch me_ ;
+  MPosition mp_ ;
   int coord_ ;
 
+  bool freqRefFromVREF_ ;
+
   // Logger
-  //LogIO os ;
+  LogIO os_ ;
 };