Changeset 1388 for branches/alma


Ignore:
Timestamp:
07/27/07 02:02:13 (17 years ago)
Author:
TakTsutsumi
Message:

merged from NRAO version of ASAP2.1 with ALMA specific modifications

Location:
branches/alma/src
Files:
8 edited

Legend:

Unmodified
Added
Removed
  • branches/alma/src/STDefs.h

    r1103 r1388  
    4141
    4242  enum Instrument {UNKNOWNINST=0,
     43                   ALMA,
    4344                   ATPKSMB,
    4445                   ATPKSHOH,
     
    4647                   TIDBINBILLA,
    4748                   CEDUNA,
     49                   GBT,
    4850                   HOBART,
    4951                   N_INSTRUMENTS};
  • branches/alma/src/STFiller.h

    r1353 r1388  
    9999  casa::String filename_;
    100100  casa::CountedPtr< Scantable > table_;
    101   casa::Int nIF_, nBeam_, nPol_, nChan_;
     101  casa::Int nIF_, nBeam_, nPol_, nChan_, nInDataRow;
    102102  casa::uInt ifOffset_, beamOffset_;
    103103  casa::Vector<casa::Bool> haveXPol_;
  • branches/alma/src/STFitter.h

    r1386 r1388  
    6464  void reset();
    6565  bool fit();
     66  // Fit via linear method
     67  bool lfit();
    6668  bool computeEstimate();
    6769
  • branches/alma/src/STMath.h

    r1374 r1388  
    132132                                        bool preserve = true );
    133133
     134  /**
     135    * Calibrate total power scans (translated from GBTIDL)
     136    * @param calon uncalibrated Scantable with CAL noise signal
     137    * @param caloff uncalibrated Scantable with no CAL signal
     138    * @param tcal optional scalar Tcal, CAL temperature (K)
     139    * @return casa::CountedPtr<Scantable> which holds a calibrated Scantable
     140    * (spectrum - average of the two CAL on and off spectra;
     141    * tsys - mean Tsys = <caloff>*Tcal/<calon-caloff> + Tcal/2)
     142    */
     143  casa::CountedPtr<Scantable> dototalpower( const casa::CountedPtr<Scantable>& calon,
     144                                            const casa::CountedPtr<Scantable>& caloff,
     145                                            casa::Float tcal=1.0 );
     146
     147  /**
     148    * Combine signal and reference scans (translated from GBTIDL)
     149    * @param sig Scantable which contains signal scans
     150    * @param ref Scantable which contains reference scans
     151    * @param smoothref optional Boxcar smooth width of the reference scans
     152    * default: no smoothing (=1)
     153    * @param tsysv optional scalar Tsys value at the zenith, required to
     154    * set tau, as well
     155    * @param tau optional scalar Tau value
     156    * @return casa::CountedPtr<Scantable> which holds combined scans
     157    * (spectrum = (sig-ref)/ref * Tsys )
     158    */
     159  casa::CountedPtr<Scantable> dosigref( const casa::CountedPtr<Scantable>& sig,
     160                                        const casa::CountedPtr<Scantable>& ref,
     161                                        int smoothref=1,
     162                                        casa::Float tsysv=0.0,
     163                                        casa::Float tau=0.0 );
     164
     165 /**
     166    * Calibrate GBT Nod scan pairs (translated from GBTIDL)
     167    * @param s Scantable which contains Nod scans
     168    * @param scans Vector of scan numbers
     169    * @param smoothref optional Boxcar smooth width of the reference scans
     170    * @param tsysv optional scalar Tsys value at the zenith, required to
     171    * set tau, as well
     172    * @param tau optional scalar Tau value
     173    * @param tcal optional scalar Tcal, CAL temperature (K)
     174    * @return casa::CountedPtr<Scantable> which holds calibrated scans
     175    */
     176  casa::CountedPtr<Scantable> donod( const casa::CountedPtr<Scantable>& s,
     177                                     const std::vector<int>& scans,
     178                                     int smoothref=1,
     179                                     casa::Float tsysv=0.0,
     180                                     casa::Float tau=0.0,
     181                                     casa::Float tcal=0.0 );
     182
     183  /**
     184    * Calibrate frequency switched scans (translated from GBTIDL)
     185    * @param s Scantable which contains frequency switched  scans
     186    * @param scans Vector of scan numbers
     187    * @param smoothref optional Boxcar smooth width of the reference scans
     188    * @param tsysv optional scalar Tsys value at the zenith, required to
     189    * set tau, as well
     190    * @param tau optional scalar Tau value
     191    * @param tcal optional scalar Tcal, CAL temperature (K)
     192    * @return casa::CountedPtr<Scantable> which holds calibrated scans
     193    */
     194  casa::CountedPtr<Scantable> dofs( const casa::CountedPtr<Scantable>& s,
     195                                    const std::vector<int>& scans,
     196                                    int smoothref=1,
     197                                    casa::Float tsysv=0.0,
     198                                    casa::Float tau=0.0,
     199                                    casa::Float tcal=0.0 );
     200
     201
    134202  casa::CountedPtr<Scantable>
    135203    freqSwitch( const casa::CountedPtr<Scantable>& in );
  • branches/alma/src/STMathWrapper.h

    r1353 r1388  
    9191                                               preserve ) ); }
    9292
     93  ScantableWrapper dototalpower( const ScantableWrapper& calon,
     94                             const ScantableWrapper& caloff, casa::Float tcal= 0 )
     95  { return ScantableWrapper( STMath::dototalpower( calon.getCP(), caloff.getCP(), tcal ) ); }
     96
     97  ScantableWrapper dosigref( const ScantableWrapper& sig,
     98                             const ScantableWrapper& ref,
     99                             int smoothref = 0, casa::Float tsysv=0.0, casa::Float tau=0.0)
     100  { return ScantableWrapper( STMath::dosigref( sig.getCP(), ref.getCP(), smoothref, tsysv, tau ) ); }
     101
     102  ScantableWrapper donod( const ScantableWrapper& s,
     103                          const std::vector<int>& scans,
     104                          int smoothref = 0,
     105                          casa::Float tsysv=0.0, casa::Float tau=0.0, casa::Float tcal=0.0 )
     106  { return ScantableWrapper( STMath::donod( s.getCP(), scans, smoothref, tsysv, tau, tcal ) ); }
     107
     108  ScantableWrapper dofs( const ScantableWrapper& s,
     109                         const std::vector<int>& scans,
     110                         int smoothref = 0,
     111                         casa::Float tsysv=0.0, casa::Float tau=0.0, casa::Float tcal=0.0 )
     112  { return ScantableWrapper( STMath::dofs( s.getCP(), scans, smoothref, tsysv, tau, tcal ) ); }
     113
    93114  ScantableWrapper
    94115    freqSwitch( const ScantableWrapper& in )
  • branches/alma/src/STWriter.h

    r1353 r1388  
    8181                      casa::Vector<casa::Complex>& xpol,
    8282                      const casa::Table& tab);
     83
     84  void replacePtTab(const casa::Table& tab, const std::string& fname);
     85
    8386  std::string     format_;
    8487  PKSwriter* writer_;
  • branches/alma/src/Scantable.h

    r1385 r1388  
    381381    { STFitEntry fe; fitTable_.getEntry(fe, mfitidCol_(row)); return fe; }
    382382
     383  //Added by TT
     384  /**
     385   * Get the antenna name
     386   * @return antenna name string
     387   */
     388  std::string getAntennaName() const;
     389
     390  /**
     391   * For GBT MS data only. check a scan list
     392   * against the information found in GBT_GO table for
     393   * scan number orders to get correct pairs.
     394   * @param[in] scan list
     395   * @return status
     396   */
     397  int checkScanInfo(const std::vector<int>& scanlist) const;
     398
     399  /**
     400   * Get the direction as a vector, for a specific row
     401   * @param[in] whichrow the row numbyyer
     402   * @return the direction in a vector
     403   */
     404  std::vector<double> getDirectionVector(int whichrow) const;
     405
    383406private:
    384407
  • branches/alma/src/ScantableWrapper.h

    r1385 r1388  
    199199    { return table_->columnNames(); }
    200200
     201  std::string getAntennaName() const
     202    { return table_->getAntennaName(); }
     203
     204  int checkScanInfo(const vector<int>& scanlist) const
     205    { return table_->checkScanInfo(scanlist); }
     206
     207  std::vector<double>  getDirectionVector(int whichrow) const
     208    { return table_->getDirectionVector(whichrow); }
     209
    201210private:
    202211  casa::CountedPtr<Scantable> table_;
Note: See TracChangeset for help on using the changeset viewer.