Changeset 1633 for branches/alma/src/STMath.h

Timestamp:

09/16/09 17:11:16 (15 years ago)

Author:

Takeshi Nakazato

Message:

New Development: No

JIRA Issue: Yes CAS-1422

Ready to Release: Yes

Interface Changes: Yes

What Interface Changed: Defined calibrate() and almacal() in python/asapmath.py

Defined cwcal() in STMath class

Test Programs: List test programs

Put in Release Notes: Yes

Module(s): task_sdaverage

Description: Describe your changes here...

asapmath.py is changed to be able to calibrate data for APEX.
This modification is tentatively considered a calibration of ALMA
single-dish data. However, it must be updated in the future since
I don't know how raw ALMA data are provided and I have to change
code along with read ALMA data.

The calibrate() function takes calibration mode from its argument and
looks antenna name of the input scantable, and calls appropriate calibration
function depending on the calibration mode and antenna name.
If antenna name include 'APEX' or 'ALMA', newly defined calibration function
apexcal() is called. For other antenna name, one of the existing calibration
functions (calps, calnod, calfs, auto_quotient) is called.

File:

• branches/alma/src/STMath.h (modified) (3 diffs)

branches/alma/src/STMath.h

@@ -200 +200 @@
 
   /**
+   * Calibrate data with Chopper-Wheel like calibration method 
+   * which adopts position switching by antenna motion, 
+   * wobbler (nutator) switching and On-The-Fly observation.
+   * 
+   * The method is applicable to APEX, and other telescopes other than GBT.
+   *
+   * @param a Scantable which contains ON and OFF scans
+   * @param a string that indicates calibration mode 
+   * @param a string that indicates antenna name
+   **/
+  casa::CountedPtr<Scantable> cwcal( const casa::CountedPtr<Scantable>& s,
+                                       const casa::String calmode, 
+                                       const casa::String antname );
+
+  /**
+   * Calibrate frequency switched scans with Chopper-Wheel like 
+   * calibration method.
+   *
+   * The method is applicable to APEX, and other telescopes other than GBT.
+   * 
+   * @param a Scantable which contains ON and OFF scans
+   * @param a string that indicates antenna name
+   **/
+  casa::CountedPtr<Scantable> cwcalfs( const casa::CountedPtr<Scantable>& s,
+                                       const casa::String antname );
+
+
+  /**
    * Folding frequency-switch data
    * @param sig
@@ -207 +235 @@
   casa::CountedPtr<Scantable> dofold( const casa::CountedPtr<Scantable> &sig,
                                       const casa::CountedPtr<Scantable> &ref,
-                                      casa::Int choffset );
+                                      casa::Double choffset,
+                                      casa::Double choffset = 0.0 );
 
   casa::CountedPtr<Scantable>
@@ -330 +359 @@
     flagsFromMA(const casa::MaskedArray<casa::Float>& ma);
 
+  vector<float> getSpectrumFromTime( string reftime, casa::CountedPtr<Scantable>& s, string mode = "before" ) ;
+  vector<float> getTcalFromTime( string reftime, casa::CountedPtr<Scantable>& s, string mode="before" ) ;
+  vector<int> getRowIdFromTime( string reftime, casa::CountedPtr<Scantable>& s ) ;
+  vector<float> getCalibratedSpectra( casa::CountedPtr<Scantable>& on,
+                                      casa::CountedPtr<Scantable>& off,
+                                      casa::CountedPtr<Scantable>& sky,
+                                      casa::CountedPtr<Scantable>& hot,
+                                      casa::CountedPtr<Scantable>& cold,
+                                      int index,
+                                      string antname ) ;
+  vector<float> getFSCalibratedSpectra( casa::CountedPtr<Scantable>& sig,
+                                        casa::CountedPtr<Scantable>& ref,
+                                        casa::CountedPtr<Scantable>& sky,
+                                        casa::CountedPtr<Scantable>& hot,
+                                        casa::CountedPtr<Scantable>& cold,
+                                        int index ) ;
+  vector<float> getFSCalibratedSpectra( casa::CountedPtr<Scantable>& sig,
+                                        casa::CountedPtr<Scantable>& ref,
+                                        vector< casa::CountedPtr<Scantable> >& sky,
+                                        vector< casa::CountedPtr<Scantable> >& hot,
+                                        vector< casa::CountedPtr<Scantable> >& cold,
+                                        int index ) ;
+  double getMJD( string strtime ) ;
+
   bool insitu_;
 };