Changeset 1779 for branches/mergetest/src/STMath.h

Timestamp:

07/29/10 19:13:46 (14 years ago)

Author:

Kana Sugimoto

Message:

New Development: Yes

JIRA Issue: No (test merging alma branch)

Ready for Test: Yes

Interface Changes: Yes

What Interface Changed:

Test Programs:

Put in Release Notes: No

Module(s):

Description:

Location:

branches/mergetest

Files:

• . (modified) (1 prop)
• src/STMath.h (modified) (9 diffs)

branches/mergetest/src/STMath.h

@@ -119 +119 @@
                   const std::string& mode, bool tsys=false );
 
+  // array operation
+  casa::CountedPtr<Scantable>
+    arrayOperate( const casa::CountedPtr<Scantable>& in,
+                  const std::vector<float> val,
+                  const std::string& mode,
+                  const std::string& opmode="channel",  
+                  bool tsys=false );
+
+  // channel operation
+  casa::CountedPtr<Scantable>
+    arrayOperateChannel( const casa::CountedPtr<Scantable>& in, 
+                         const std::vector<float> val,
+                         const std::string& mode, bool tsys=false );
+
+  // row operation
+  casa::CountedPtr<Scantable>
+    arrayOperateRow( const casa::CountedPtr<Scantable>& in, 
+                     const std::vector<float> val,
+                     const std::string& mode, bool tsys=false );
+
+  // 2d array operation
+  casa::CountedPtr<Scantable>
+    array2dOperate( const casa::CountedPtr<Scantable>& in,
+                  const std::vector< std::vector<float> > val,
+                  const std::string& mode, bool tsys=false );
+
   casa::CountedPtr<Scantable>
     binaryOperate( const casa::CountedPtr<Scantable>& left,
@@ -134 +160 @@
   /**
     * Calibrate total power scans (translated from GBTIDL)
-    * @param calon uncalibrated Scantable with CAL noise signal
+    * @param calon uncalibrated Scantable with CAL noise signal 
     * @param caloff uncalibrated Scantable with no CAL signal
     * @param tcal optional scalar Tcal, CAL temperature (K)
-    * @return casa::CountedPtr<Scantable> which holds a calibrated Scantable
+    * @return casa::CountedPtr<Scantable> which holds a calibrated Scantable 
     * (spectrum - average of the two CAL on and off spectra;
     * tsys - mean Tsys = <caloff>*Tcal/<calon-caloff> + Tcal/2)
-    */
+    */              
   casa::CountedPtr<Scantable> dototalpower( const casa::CountedPtr<Scantable>& calon,
                                             const casa::CountedPtr<Scantable>& caloff,
@@ -151 +177 @@
     * @param smoothref optional Boxcar smooth width of the reference scans
     * default: no smoothing (=1)
-    * @param tsysv optional scalar Tsys value at the zenith, required to
-    * set tau, as well
+    * @param tsysv optional scalar Tsys value at the zenith, required to 
+    * set tau, as well 
     * @param tau optional scalar Tau value
     * @return casa::CountedPtr<Scantable> which holds combined scans
@@ -163 +189 @@
                                         casa::Float tau=0.0 );
 
- /**
+  /**
     * Calibrate GBT Nod scan pairs (translated from GBTIDL)
     * @param s Scantable which contains Nod scans
@@ -170 +196 @@
     * @param tsysv optional scalar Tsys value at the zenith, required to
     * set tau, as well
-    * @param tau optional scalar Tau value
+    * @param tau optional scalar Tau value 
     * @param tcal optional scalar Tcal, CAL temperature (K)
     * @return casa::CountedPtr<Scantable> which holds calibrated scans
@@ -199 +225 @@
                                     casa::Float tcal=0.0 );
 
+  /**
+   * 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
+   * @param ref
+   * @param choffset
+   **/
+  casa::CountedPtr<Scantable> dofold( const casa::CountedPtr<Scantable> &sig,
+                                      const casa::CountedPtr<Scantable> &ref,
+                                      casa::Double choffset,
+                                      casa::Double choffset = 0.0 );
+
+  /**
+   * ALMA calibration
+   **/
+  casa::CountedPtr<Scantable> almacal( const casa::CountedPtr<Scantable>& s, 
+                                       const casa::String calmode ) ;
+  casa::CountedPtr<Scantable> almacalfs( const casa::CountedPtr<Scantable>& s ) ;
 
   casa::CountedPtr<Scantable>
@@ -206 +277 @@
                                const std::vector<bool>& mask,
                                const std::string& which);
+
+  std::vector< int > minMaxChan(const casa::CountedPtr<Scantable>& in,
+                                const std::vector<bool>& mask,
+                                const std::string& which);
 
   casa::CountedPtr<Scantable> bin( const casa::CountedPtr<Scantable>& in,
@@ -266 +341 @@
              double end, const std::string& mode="frequency");
 
+  // test for average spectra with different channel/resolution
+  casa::CountedPtr<Scantable>
+    new_average( const std::vector<casa::CountedPtr<Scantable> >& in,
+                 const bool& compel, 
+                 const std::vector<bool>& mask = std::vector<bool>(),
+                 const std::string& weight = "NONE",
+                 const std::string& avmode = "SCAN" )
+    throw (casa::AipsError) ;
+
 private:
   casa::CountedPtr<Scantable>  applyToPol( const casa::CountedPtr<Scantable>& in,
@@ -301 +385 @@
     maskedArray( const casa::Vector<casa::Float>& s,
                  const casa::Vector<casa::uChar>& f );
+  casa::MaskedArray<casa::Double>
+    maskedArray( const casa::Vector<casa::Double>& s,
+                 const casa::Vector<casa::uChar>& f );
   casa::Vector<casa::uChar>
     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<float> getTsysFromTime( string reftime, casa::CountedPtr<Scantable>& s, string mode="before" ) ;
+  vector<int> getRowIdFromTime( string reftime, casa::CountedPtr<Scantable>& s ) ;
+
+  // Chopper-Wheel type calibration
+  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 ) ;
+  // Tsys * (ON-OFF)/OFF
+  vector<float> getCalibratedSpectra( casa::CountedPtr<Scantable>& on,
+                                      casa::CountedPtr<Scantable>& off,
+                                      int index ) ;
+  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_;
 };