Changeset 1779 for branches/mergetest/src/STFiller.cpp

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/STFiller.cpp (modified) (18 diffs)

branches/mergetest/src/STFiller.cpp

@@ -25 +25 @@
 #include <tables/Tables/TableRow.h>
 
+#include <measures/Measures/MDirection.h>
+#include <measures/Measures/MeasConvert.h>
+
 #include <atnf/PKSIO/PKSrecord.h>
 #include <atnf/PKSIO/PKSreader.h>
@@ -30 +33 @@
  #include <casa/System/ProgressMeter.h>
 #endif
+#include <casa/System/ProgressMeter.h>
+#include <atnf/PKSIO/NROReader.h>
+#include <casa/Logging/LogIO.h>
+
+#include <time.h>
+
 
 #include "STDefs.h"
@@ -45 +54 @@
   header_(0),
   table_(0),
-  refRx_(".*(e|w|_R)$")
+  refRx_(".*(e|w|_R)$"),
+  nreader_(0)
 {
 }
@@ -53 +63 @@
   header_(0),
   table_(stbl),
-  refRx_(".*(e|w|_R)$")
+  refRx_(".*(e|w|_R)$"),
+  nreader_(0)
 {
 }
@@ -61 +72 @@
   header_(0),
   table_(0),
-  refRx_(".*(e|w|_R)$")
+  refRx_(".*(e|w|_R)$"),
+  nreader_(0)
 {
-  open(filename, whichIF, whichBeam);
+  open(filename, "", whichIF, whichBeam);
 }
 
@@ -70 +82 @@
   close();
 }
-
-void STFiller::open( const std::string& filename, int whichIF, int whichBeam )
+  
+void STFiller::open( const std::string& filename, const std::string& antenna, int whichIF, int whichBeam, casa::Bool getPt )
 {
   if (table_.null())  {
@@ -93 +105 @@
   Vector<Bool> beams, ifs;
   Vector<uInt> nchans,npols;
-  if ( (reader_ = getPKSreader(inName, 0, 0, format, beams, ifs,
+
+  //
+  // if isNRO_ is true, try NROReader
+  //
+  // 2008/11/11 Takeshi Nakazato
+  isNRO_ = fileCheck() ;
+  if ( isNRO_ ) {
+    if ( (nreader_ = getNROReader( inName, format )) == 0 ) {
+      throw(AipsError("Creation of NROReader failed")) ;
+    }
+    else {
+      openNRO( whichIF, whichBeam ) ;
+      return ;
+    }
+  }
+  //
+
+  if ( (reader_ = getPKSreader(inName, antenna, 0, 0, format, beams, ifs,
                               nchans, npols, haveXPol_,haveBase, haveSpectra
                               )) == 0 )  {
@@ -197 +226 @@
   Vector<Int> start(nIF_, 1);
   Vector<Int> end(nIF_, 0);
-  reader_->select(beams, ifs, start, end, ref, True, haveXPol_[0]);
+  reader_->select(beams, ifs, start, end, ref, True, haveXPol_[0], False, getPt);
   table_->setHeader(*header_);
   //For MS, add the location of POINTING of the input MS so one get
   //pointing data from there, if necessary.
-  //Also find nrow in MS
+  //Also find nrow in MS 
   nInDataRow = 0;
   if (format == "MS2") {
-    Path datapath(inName);
+    Path datapath(inName); 
     String ptTabPath = datapath.absoluteName();
     Table inMS(ptTabPath);
@@ -216 +245 @@
     }
   }
+  String freqFrame = header_->freqref;
+  //translate frequency reference frame back to 
+  //MS style (as PKSMS2reader converts the original frame 
+  //in FITS standard style)
+  if (freqFrame == "TOPOCENT") {
+    freqFrame = "TOPO";
+  } else if (freqFrame == "GEOCENER") {
+    freqFrame = "GEO";
+  } else if (freqFrame == "BARYCENT") {
+    freqFrame = "BARY";
+  } else if (freqFrame == "GALACTOC") {
+    freqFrame = "GALACTO";
+  } else if (freqFrame == "LOCALGRP") {
+    freqFrame = "LGROUP";
+  } else if (freqFrame == "CMBDIPOL") {
+    freqFrame = "CMB";
+  } else if (freqFrame == "SOURCE") {
+    freqFrame = "REST";
+  }
+  // set both "FRAME" and "BASEFRAME"
+  table_->frequencies().setFrame(freqFrame, false);
+  table_->frequencies().setFrame(freqFrame,true);
   //table_->focus().setParallactify(true);
 }
@@ -222 +273 @@
 {
   delete reader_;reader_=0;
+  delete nreader_;nreader_=0;
   delete header_;header_=0;
   table_ = 0;
@@ -229 +281 @@
 {
   int status = 0;
+
+  // 
+  // for NRO data
+  //
+  // 2008/11/12 Takeshi Nakazato
+  if ( isNRO_ ) {
+    status = readNRO() ;
+    return status ;
+  }
+  //
+
+/**
+  Int    beamNo, IFno, refBeam, scanNo, cycleNo;
+  Float  azimuth, elevation, focusAxi, focusRot, focusTan,
+    humidity, parAngle, pressure, temperature, windAz, windSpeed;
+  Double bandwidth, freqInc, interval, mjd, refFreq, srcVel;
+  String          fieldName, srcName, tcalTime, obsType;
+  Vector<Float>   calFctr, sigma, tcal, tsys;
+  Matrix<Float>   baseLin, baseSub;
+  Vector<Double>  direction(2), scanRate(2), srcDir(2), srcPM(2), restFreq(1);
+  Matrix<Float>   spectra;
+  Matrix<uChar>   flagtra;
+  Complex         xCalFctr;
+  Vector<Complex> xPol;
+**/
 
   Double min = 0.0;
@@ -236 +313 @@
 #endif
   PKSrecord pksrec;
+  pksrec.srcType=-1;
   int n = 0;
+  bool isGBTFITS = false ;
+  if ((header_->antennaname.find( "GBT" ) != String::npos) && File(filename_).isRegular()) {
+    FILE *fp = fopen( filename_.c_str(), "r" ) ;
+    fseek( fp, 640, SEEK_SET ) ;
+    char buf[81] ;
+    fread( buf, 80, 1, fp ) ;
+    buf[80] = '\0' ;
+    if ( strstr( buf, "NRAO_GBT" ) != NULL ) {
+      isGBTFITS = true ;
+    }
+    fclose( fp ) ;
+  } 
   while ( status == 0 ) {
     status = reader_->read(pksrec);
@@ -288 +378 @@
     //*srcnCol = pksrec.srcName;//.before(rx2);
     *srctCol = match;
+    if ( pksrec.srcType != -1 ) {
+      *srctCol = pksrec.srcType ;
+    }
     RecordFieldPtr<uInt> beamCol(rec, "BEAMNO");
     *beamCol = pksrec.beamNo-beamOffset_-1;
@@ -298 +391 @@
     uInt id;
     /// @todo this has to change when nchan isn't global anymore
-    id = table_->frequencies().addEntry(Double(header_->nchan/2),
-                                        pksrec.refFreq, pksrec.freqInc);
+    //id = table_->frequencies().addEntry(Double(header_->nchan/2),
+    //                                    pksrec.refFreq, pksrec.freqInc);
+    if ( pksrec.nchan == 1 ) {
+      id = table_->frequencies().addEntry(Double(0),
+                                          pksrec.refFreq, pksrec.freqInc);
+    }
+    else {
+      id = table_->frequencies().addEntry(Double(pksrec.nchan/2),
+                                          pksrec.refFreq, pksrec.freqInc);
+    }
     RecordFieldPtr<uInt> mfreqidCol(rec, "FREQ_ID");
     *mfreqidCol = id;
+    //*ifCol = id;
 
     id = table_->molecules().addEntry(pksrec.restFreq);
@@ -317 +419 @@
 
     RecordFieldPtr<uInt> mfocusidCol(rec, "FOCUS_ID");
-    id = table_->focus().addEntry(pksrec.parAngle, pksrec.focusAxi,
+    id = table_->focus().addEntry(pksrec.parAngle, pksrec.focusAxi, 
                                   pksrec.focusTan, pksrec.focusRot);
     *mfocusidCol = id;
@@ -335 +437 @@
     // into 2-4 rows in the scantable
     Vector<Float> tsysvec(1);
-    // Why is pksrec.spectra.ncolumn() == 3 for haveXPol_ == True
+    // Why is spectra.ncolumn() == 3 for haveXPol_ == True
     uInt npol = (pksrec.spectra.ncolumn()==1 ? 1: 2);
     for ( uInt i=0; i< npol; ++i ) {
       tsysvec = pksrec.tsys(i);
       *tsysCol = tsysvec;
-      *polnoCol = i;
+      if (isGBTFITS)
+        *polnoCol = pksrec.polNo ;
+      else
+        *polnoCol = i;
 
       *specCol = pksrec.spectra.column(i);
@@ -347 +452 @@
       row.put(table_->table().nrow()-1, rec);
     }
+
+    RecordFieldPtr< uInt > flagrowCol(rec, "FLAGROW");
+    *flagrowCol = pksrec.flagrow;
+
     if ( haveXPol_[0] ) {
       // no tsys given for xpol, so emulate it
@@ -381 +490 @@
 }
 
+/**
+ * For NRO data
+ *
+ * 2008/11/11 Takeshi Nakazato
+ **/
+void STFiller::openNRO( int whichIF, int whichBeam )
+{
+  // open file
+  // DEBUG
+  time_t t0 ;
+  time( &t0 ) ;
+  tm *ttm = localtime( &t0 ) ;
+  LogIO os( LogOrigin( "STFiller", "openNRO()", WHERE ) ) ;
+//   cout << "STFiller::openNRO()  Start time = " << t0 
+//        << " (" 
+//        << ttm->tm_year + 1900 << "/" << ttm->tm_mon + 1 << "/" << ttm->tm_mday 
+//        << " " 
+//        << ttm->tm_hour << ":" << ttm->tm_min << ":" << ttm->tm_sec 
+//        << ")" << endl ;
+  os << "Start time = " << t0 
+     << " (" 
+     << ttm->tm_year + 1900 << "/" << ttm->tm_mon + 1 << "/" << ttm->tm_mday 
+     << " " 
+     << ttm->tm_hour << ":" << ttm->tm_min << ":" << ttm->tm_sec 
+     << ")" << LogIO::POST ;
+
+  // fill STHeader
+  header_ = new STHeader() ;
+
+  if ( nreader_->getHeaderInfo( header_->nchan,
+                                header_->npol,
+                                nIF_,
+                                nBeam_,
+                                header_->observer,
+                                header_->project,
+                                header_->obstype,
+                                header_->antennaname,
+                                header_->antennaposition,
+                                header_->equinox,
+                                header_->freqref,
+                                header_->reffreq,
+                                header_->bandwidth,
+                                header_->utc,
+                                header_->fluxunit,
+                                header_->epoch,
+                                header_->poltype ) ) {
+//     cout << "STFiller::openNRO()  Failed to get header information." << endl ;
+//     return ;
+    throw( AipsError("Failed to get header information.") ) ;
+  }
+
+  // set FRAME and BASEFRAME keyword of FREQUENCIES table
+  if ( header_->freqref != "TOPO" ) {
+    table_->frequencies().setFrame( header_->freqref, false ) ;
+    table_->frequencies().setFrame( header_->freqref, true ) ;
+  }
+
+  ifOffset_ = 0;
+  vector<Bool> ifs = nreader_->getIFs() ;
+  if ( whichIF >= 0 ) {
+    if ( whichIF >= 0 && whichIF < nIF_ ) {
+      for ( int i = 0 ; i < nIF_ ; i++ ) 
+        ifs[i] = False ;
+      ifs[whichIF] = True ;
+      header_->nif = 1;
+      nIF_ = 1;
+      ifOffset_ = whichIF;
+    } else {
+      delete reader_;
+      reader_ = 0;
+      delete header_;
+      header_ = 0;
+      throw(AipsError("Illegal IF selection"));
+    }
+  }
+
+  beamOffset_ = 0;
+  vector<Bool> beams = nreader_->getBeams() ;
+  if (whichBeam>=0) {
+    if (whichBeam>=0 && whichBeam<nBeam_) {
+      for ( int i = 0 ; i < nBeam_ ; i++ ) 
+        beams[i] = False ;
+      beams[whichBeam] = True;
+      header_->nbeam = 1;
+      nBeam_ = 1;
+      beamOffset_ = whichBeam;
+    } else {
+      delete reader_;
+      reader_ = 0;
+      delete header_;
+      header_ = 0;
+      throw(AipsError("Illegal Beam selection"));
+    }
+  }
+
+  header_->nbeam = nBeam_ ;
+  header_->nif = nIF_ ;
+
+  // set header
+  table_->setHeader( *header_ ) ;
+
+  // DEBUG
+  time_t t1 ;
+  time( &t1 ) ;
+  ttm = localtime( &t1 ) ;
+//   cout << "STFiller::openNRO()  End time = " << t1 
+//        << " (" 
+//        << ttm->tm_year + 1900 << "/" << ttm->tm_mon + 1 << "/" << ttm->tm_mday 
+//        << " " 
+//        << ttm->tm_hour << ":" << ttm->tm_min << ":" << ttm->tm_sec 
+//        << ")" << endl ;
+//   cout << "STFiller::openNRO()  Elapsed time = " << t1 - t0 << " sec" << endl ;
+  os << "End time = " << t1 
+     << " (" 
+     << ttm->tm_year + 1900 << "/" << ttm->tm_mon + 1 << "/" << ttm->tm_mday 
+     << " " 
+     << ttm->tm_hour << ":" << ttm->tm_min << ":" << ttm->tm_sec 
+     << ")" << endl ;
+  os << "Elapsed time = " << t1 - t0 << " sec" << endl ;
+  os.post() ;
+  //
+
+  return ;
+}
+
+int STFiller::readNRO()
+{
+  // DEBUG
+  time_t t0 ;
+  time( &t0 ) ;
+  tm *ttm = localtime( &t0 ) ;
+  LogIO os( LogOrigin( "STFiller", "readNRO()", WHERE ) ) ;
+//   cout << "STFiller::readNRO()  Start time = " << t0 
+//        << " (" 
+//        << ttm->tm_year + 1900 << "/" << ttm->tm_mon + 1 << "/" << ttm->tm_mday 
+//        << " " 
+//        << ttm->tm_hour << ":" << ttm->tm_min << ":" << ttm->tm_sec 
+//        << ")" << endl ;
+  os << "Start time = " << t0 
+     << " (" 
+     << ttm->tm_year + 1900 << "/" << ttm->tm_mon + 1 << "/" << ttm->tm_mday 
+     << " " 
+     << ttm->tm_hour << ":" << ttm->tm_min << ":" << ttm->tm_sec 
+     << ")" << LogIO::POST ;
+  //
+
+  // fill row
+  uInt id ;
+  uInt imax = nreader_->getRowNum() ;
+  vector< vector<double > > freqs ;
+  uInt i = 0 ;
+  int count = 0 ;
+  uInt scanno ;
+  uInt cycleno ;
+  uInt beamno ;
+  uInt polno ;
+  vector<double> fqs ;
+  Vector<Double> restfreq ;
+  uInt refbeamno ;
+  Double scantime ;
+  Double interval ;
+  String srcname ;
+  String fieldname ;
+  Array<Float> spectra ;
+  Array<uChar> flagtra ;
+  Array<Float> tsys ;
+  Array<Double> direction ;
+  Float azimuth ;
+  Float elevation ;
+  Float parangle ;
+  Float opacity ;
+  uInt tcalid ;
+  Int fitid ;
+  uInt focusid ;
+  Float temperature ;
+  Float pressure ;
+  Float humidity ;
+  Float windvel ;
+  Float winddir ;
+  Double srcvel ;
+  Array<Double> propermotion ;
+  Vector<Double> srcdir ;
+  Array<Double> scanrate ;
+  for ( i = 0 ; i < imax ; i++ ) {
+    string scanType = nreader_->getScanType( i ) ;
+    Int srcType = -1 ;
+    if ( scanType.compare( 0, 2, "ON") == 0 ) {
+      // os << "ON srcType: " << i << LogIO::POST ;
+      srcType = 0 ;
+    }
+    else if ( scanType.compare( 0, 3, "OFF" ) == 0 ) {
+      //os << "OFF srcType: " << i << LogIO::POST ;
+      srcType = 1 ;
+    }
+    else if ( scanType.compare( 0, 4, "ZERO" ) == 0 ) {
+      //os << "ZERO srcType: " << i << LogIO::POST ;
+      srcType = 2 ;
+    }
+    else {
+      //os << "Undefined srcType: " << i << LogIO::POST ;
+      srcType = 3 ;
+    }
+ 
+    // if srcType is 2 (ZERO scan), ignore scan
+    if ( srcType != 2 && srcType != -1 && srcType != 3 ) { 
+      TableRow row( table_->table() ) ;
+      TableRecord& rec = row.record();
+
+      if ( nreader_->getScanInfo( i,
+                                  scanno,
+                                  cycleno,
+                                  beamno,
+                                  polno,
+                                  fqs,
+                                  restfreq,
+                                  refbeamno,
+                                  scantime,
+                                  interval,
+                                  srcname,
+                                  fieldname,
+                                  spectra,
+                                  flagtra,
+                                  tsys,
+                                  direction,
+                                  azimuth,
+                                  elevation,
+                                  parangle,
+                                  opacity,
+                                  tcalid,
+                                  fitid,
+                                  focusid,
+                                  temperature,
+                                  pressure,
+                                  humidity,
+                                  windvel,
+                                  winddir,
+                                  srcvel,
+                                  propermotion,
+                                  srcdir,
+                                  scanrate ) ) {
+//         cerr << "STFiller::readNRO()  Failed to get scan information." << endl ;
+//         return 1 ;
+        throw( AipsError("Failed to get scan information.") ) ;
+      }
+
+      RecordFieldPtr<uInt> scannoCol( rec, "SCANNO" ) ;
+      *scannoCol = scanno ;
+      RecordFieldPtr<uInt> cyclenoCol(rec, "CYCLENO") ;
+      *cyclenoCol = cycleno ;
+      RecordFieldPtr<uInt> beamCol(rec, "BEAMNO") ;
+      *beamCol = beamno ;
+      RecordFieldPtr<uInt> ifCol(rec, "IFNO") ;
+      RecordFieldPtr< uInt > polnoCol(rec, "POLNO") ;
+      *polnoCol = polno ;
+      RecordFieldPtr<uInt> mfreqidCol(rec, "FREQ_ID") ;
+      if ( freqs.size() == 0 ) {
+        id = table_->frequencies().addEntry( Double( fqs[0] ),
+                                             Double( fqs[1] ),
+                                             Double( fqs[2] ) ) ;
+        *mfreqidCol = id ;
+        *ifCol = id ;
+        freqs.push_back( fqs ) ;
+      }
+      else {
+        int iadd = -1 ;
+        for ( uInt iif = 0 ; iif < freqs.size() ; iif++ ) {
+          //os << "freqs[" << iif << "][1] = " << freqs[iif][1] << LogIO::POST ;
+          double fdiff = abs( freqs[iif][1] - fqs[1] ) / freqs[iif][1] ;
+          //os << "fdiff = " << fdiff << LogIO::POST ;
+          if ( fdiff < 1.0e-8 ) {
+            iadd = iif ;
+            break ;
+          }
+        }
+        if ( iadd == -1 ) {
+          id = table_->frequencies().addEntry( Double( fqs[0] ),
+                                               Double( fqs[1] ),
+                                               Double( fqs[2] ) ) ;
+          *mfreqidCol = id ;
+          *ifCol = id ;
+          freqs.push_back( fqs ) ;
+        }
+        else {
+          *mfreqidCol = iadd ;
+          *ifCol = iadd ;
+        }
+      }
+      RecordFieldPtr<uInt> molidCol(rec, "MOLECULE_ID") ;
+      id = table_->molecules().addEntry( restfreq ) ;
+      *molidCol = id ;
+      RecordFieldPtr<Int> rbCol(rec, "REFBEAMNO") ;
+      *rbCol = refbeamno ;
+      RecordFieldPtr<Double> mjdCol( rec, "TIME" ) ;
+      *mjdCol = scantime ;
+      RecordFieldPtr<Double> intervalCol( rec, "INTERVAL" ) ;
+      *intervalCol = interval ;
+      RecordFieldPtr<String> srcnCol(rec, "SRCNAME") ;
+      *srcnCol = srcname ;
+      RecordFieldPtr<Int> srctCol(rec, "SRCTYPE") ;
+      *srctCol = srcType ;     
+      RecordFieldPtr<String> fieldnCol(rec, "FIELDNAME"); 
+      *fieldnCol = fieldname ;
+      RecordFieldPtr< Array<Float> > specCol(rec, "SPECTRA") ;
+      *specCol = spectra ;
+      RecordFieldPtr< Array<uChar> > flagCol(rec, "FLAGTRA") ;
+      *flagCol = flagtra ;
+      RecordFieldPtr< Array<Float> > tsysCol(rec, "TSYS") ;
+      *tsysCol = tsys ;
+      RecordFieldPtr< Array<Double> > dirCol(rec, "DIRECTION") ;
+      *dirCol = direction ;
+      RecordFieldPtr<Float> azCol(rec, "AZIMUTH") ;
+      *azCol = azimuth ;
+      RecordFieldPtr<Float> elCol(rec, "ELEVATION") ;
+      *elCol = elevation ;
+      RecordFieldPtr<Float> parCol(rec, "PARANGLE") ;
+      *parCol = parangle ;
+      RecordFieldPtr<Float> tauCol(rec, "OPACITY") ;
+      *tauCol = opacity ;
+      RecordFieldPtr<uInt> mcalidCol(rec, "TCAL_ID") ;
+      *mcalidCol = tcalid ;
+      RecordFieldPtr<Int> fitCol(rec, "FIT_ID") ;
+      *fitCol = fitid ;
+      RecordFieldPtr<uInt> mfocusidCol(rec, "FOCUS_ID") ;
+      *mfocusidCol = focusid ;
+      RecordFieldPtr<uInt> mweatheridCol(rec, "WEATHER_ID") ;
+      id = table_->weather().addEntry( temperature,
+                                       pressure,
+                                       humidity,
+                                       windvel,
+                                       winddir ) ;
+      *mweatheridCol = id ;          
+      RecordFieldPtr<Double> svelCol(rec, "SRCVELOCITY") ;
+      *svelCol = srcvel ;
+      RecordFieldPtr<Array<Double> > spmCol(rec, "SRCPROPERMOTION") ;
+      *spmCol = propermotion ;
+      RecordFieldPtr<Array<Double> > sdirCol(rec, "SRCDIRECTION") ;
+      *sdirCol = srcdir ;
+      RecordFieldPtr<Array<Double> > srateCol(rec, "SCANRATE");
+      *srateCol = scanrate ;
+
+      table_->table().addRow() ;
+      row.put(table_->table().nrow()-1, rec) ;
+    }
+    else {
+      count++ ;
+    }
+    // DEBUG
+    //int rownum = nreader_->getRowNum() ;
+    //os << "Finished row " << i << "/" << rownum << LogIO::POST ;
+    //
+  }
+
+  // DEBUG
+  time_t t1 ;
+  time( &t1 ) ;
+  ttm = localtime( &t1 ) ;
+//   cout << "STFiller::readNRO()  Processed " << i << " rows" << endl ;
+//   cout << "STFiller::readNRO()  Added " << i - count << " rows (ignored " 
+//        << count << " \"ZERO\" scans)" << endl ;
+//   cout << "STFiller::readNRO()  End time = " << t1 
+//        << " (" 
+//        << ttm->tm_year + 1900 << "/" << ttm->tm_mon + 1 << "/" << ttm->tm_mday 
+//        << " " 
+//        << ttm->tm_hour << ":" << ttm->tm_min << ":" << ttm->tm_sec 
+//        << ")" << endl ;
+//   cout << "STFiller::readNRO()  Elapsed time = " << t1 - t0 << " sec" << endl ;
+  os << "Processed " << i << " rows" << endl ;
+  os << "Added " << i - count << " rows (ignored " 
+     << count << " \"ZERO\" scans)" << endl ;
+  os.post() ;
+  os << "End time = " << t1 
+     << " (" 
+     << ttm->tm_year + 1900 << "/" << ttm->tm_mon + 1 << "/" << ttm->tm_mday 
+     << " " 
+     << ttm->tm_hour << ":" << ttm->tm_min << ":" << ttm->tm_sec 
+     << ")" << endl ;
+  os << "Elapsed time = " << t1 - t0 << " sec" << endl ;
+  os.post() ;
+  //
+
+  return 0 ;
+} 
+
+Bool STFiller::fileCheck()
+{
+  bool bval = false ;
+
+  // if filename_ is directory, return false
+  File inFile( filename_ ) ;
+  if ( inFile.isDirectory() )
+    return bval ;
+  
+  // if beginning of header data is "RW", return true
+  // otherwise, return false ; 
+  FILE *fp = fopen( filename_.c_str(), "r" ) ;
+  char buf[9] ;
+  char buf2[80] ;
+  fread( buf, 4, 1, fp ) ;
+  buf[4] = '\0' ;
+  fseek( fp, 640, SEEK_SET ) ;
+  fread( buf2, 80, 1, fp ) ;
+  if ( ( strncmp( buf, "RW", 2 ) == 0 ) || ( strstr( buf2, "NRO45M" ) != NULL ) ) {
+    bval = true ;
+  }
+  fclose( fp ) ;
+  return bval ;
+}
+
 }//namespace asap