ReadLoudestList.C

Go to the documentation of this file.

#define LOUDEST_MAX_SIZE   1000

/*
  TString run[LOUDEST_MAX_SIZE];
  int     chunk[LOUDEST_MAX_SIZE];
  double  gps[LOUDEST_MAX_SIZE];
  TString bbh_name[LOUDEST_MAX_SIZE];
  double  ifar_sec[LOUDEST_MAX_SIZE];
  double  ifar_year[LOUDEST_MAX_SIZE];
  int     obs_time_sec[LOUDEST_MAX_SIZE];
  double  obs_time_day[LOUDEST_MAX_SIZE];
  double  expected[LOUDEST_MAX_SIZE];
  int     observed[LOUDEST_MAX_SIZE];
  double  cumul_FAP[LOUDEST_MAX_SIZE];
  double  sigma[LOUDEST_MAX_SIZE];

  int nLoudest = ReadLoudestList(ifile, run, chunk, gps, bbh_name, ifar_sec, ifar_year,
                                 obs_time_sec, obs_time_day, expected, observed, cumul_FAP, sigma); 
  
*/

int ReadLoudestList(TString ifile, TString* run, int* chunk, double* gps, TString* bbh_name, double* ifar_sec, double* ifar_year,
                    int* obs_time_sec, double* obs_time_day, double* expected, int* observed, double* cumul_FAP, double* sigma) {

  CWB::Toolbox::checkFile(ifile);

  // Open chunk list
  ifstream in;
  in.open(ifile.Data(),ios::in);
  if (!in.good()) {cout << "Error Opening File : " << ifile << endl;exit(1);}

  int isize=0;
  char str[1024];
  int fpos=0;
  while(true) {
    in.getline(str,1024);
    if (!in.good()) break;
    if(str[0] != '#') isize++;
  }
//  cout << "size " << isize << endl;
  in.clear(ios::goodbit);
  in.seekg(0, ios::beg);
  if(isize==0)  {cout << "Error : File " << ifile << " is empty" << endl;return 0;}
  if(isize>LOUDEST_MAX_SIZE) {cout << "Error : File " << ifile << " > " << LOUDEST_MAX_SIZE << endl;exit(1);}
 
  char srun[256];
  char sbbh[256];

  int k=0;
  while(true) {
    fpos=in.tellg();
    in.getline(str,1024);
    if (!in.good()) break;
    if(str[0] == '#' || str[0]=='\0') continue;
    in.seekg(fpos, ios::beg);
    in >> srun >> chunk[k] >> gps[k] >> sbbh >> ifar_sec[k] >> ifar_year[k] >> obs_time_sec[k] 
       >> obs_time_day[k] >> expected[k] >> observed[k] >> cumul_FAP[k] >> sigma[k];
    if(!in.good()) break;
    run[k]=srun;
    bbh_name[k]=sbbh;
    //cout << "\t" << run[k] << "\t" << chunk[k] << "\t" << gps[k] << "\t" << bbh_name[k] << "\t" << ifar_sec[k] << "\t" << ifar_year[k] << "\t" << obs_time_sec[k] 
    //     << "\t" << obs_time_day[k] << "\t" << expected[k] << "\t" << observed[k] << "\t" << cumul_FAP[k] << "\t" << sigma[k] << endl;
    k++;
    if(k>=LOUDEST_MAX_SIZE) {
      cout << "WARNING: loudest events exceed LOUDEST_MAX_SIZE = " << LOUDEST_MAX_SIZE << endl; 
      break;
    }
  }
  in.close();

  return k;
}