Changeset 3872 in Sophya for trunk/AddOn/TAcq/branap.cc


Ignore:
Timestamp:
Aug 28, 2010, 12:43:04 AM (15 years ago)
Author:
ansari
Message:

Ajout classe BRFFTCalculator et programme specmfib.cc, Reza 28/08/2010

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/AddOn/TAcq/branap.cc

    r3776 r3872  
    2525  nzones_=nzon;
    2626  npaqinzone_=npaqz;
     27  fgdatafft_=false;  fgsinglechannel_=false;
    2728  prtlevel_=0;
    2829  nbcalgrp_=1;
     
    150151  char flnm[1024];
    151152  sprintf(flnm,"%s/signal%d.fits",dirlist_[0].c_str(),imin_);
    152   return DecodeMiniFitsHeader(flnm,paqsize_, npaqinfile_);
     153  return DecodeMiniFitsHeader(flnm,paqsize_, npaqinfile_,fgdatafft_, fgsinglechannel_);
    153154}
    154155
     
    165166  cout << " PaqSize=" << paqsize_ << "  -  NZones=" << nzones_ << " NPaqZone=" << npaqinzone_
    166167       << " PrtLevel=" << prtlevel_ << endl;
     168  cout << " AcqMode: " << ((fgdatafft_)?" Data_FFT " : " Data_Raw " )
     169       << ((fgsinglechannel_)?" SingleChannel " : " TwoChannels " ) << endl;
    167170  cout << " NbVisibCalculator in Group: " <<  nbcalgrp_ << " with N//threads: " << nthreads_ << endl;
    168171
     
    172175
    173176/* --Fonction-- */
    174 int BRAnaParam::DecodeMiniFitsHeader(const char* filename, uint_4& paqsz, uint_4& npaq)
     177int BRAnaParam::DecodeMiniFitsHeader(const char* filename, uint_4& paqsz, uint_4& npaq,
     178                                     bool& fgdatafft, bool& fgsinglechannel)
    175179{
    176180  cout << " DecodeMiniFitsHeader - Opening file: " << filename << endl;
    177181  MiniFITSFile mff(filename, MF_Read);
     182  string acqmode=mff.GetKeyValue("ACQMODE");
    178183  cout << "DecodeMiniFitsHeader()... Type=" << mff.DataTypeToString() << " NAxis1=" << mff.NAxis1()
    179        << " NAxis2=" << mff.NAxis2() << endl;
     184       << " NAxis2=" << mff.NAxis2() << " AcqMode=" << acqmode << endl;
    180185  paqsz = mff.NAxis1();
    181186  npaq = mff.NAxis2();
     187  if (acqmode.substr(0,3)=="fft") fgdatafft=true;
     188  if (acqmode.find("1c") < acqmode.length()) fgsinglechannel=true;
    182189  return 0;
    183190}
Note: See TracChangeset for help on using the changeset viewer.