ReadChunkList.C

Go to the documentation of this file.

/*
# Copyright (C) 2019 Gabriele Vedovato
#
# This program is free software: you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program.  If not, see <https://www.gnu.org/licenses/>.
*/
 
 
#define CHUNK_MAX_SIZE     100
 
 
int ReadChunkList(TString ifile, int* chunk=NULL, double* start=NULL, double* stop=NULL) {
 
  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;exit(1);}
  if(isize>CHUNK_MAX_SIZE) {cout << "Error : File " << ifile << " > " << CHUNK_MAX_SIZE << endl;exit(1);}
 
//  int   chunk[CHUNK_MAX_SIZE];
//  double start[CHUNK_MAX_SIZE];
//  double stop[CHUNK_MAX_SIZE];
  if(chunk==NULL) chunk = new int[CHUNK_MAX_SIZE];
  if(start==NULL) start = new double[CHUNK_MAX_SIZE];
  if(stop==NULL)  stop  = new double[CHUNK_MAX_SIZE];
 
  char   note[CHUNK_MAX_SIZE][256];
  double days[CHUNK_MAX_SIZE];
  char   bmon[CHUNK_MAX_SIZE][256];
  int    bday[CHUNK_MAX_SIZE];
  char   emon[CHUNK_MAX_SIZE][256];
  int    eday[CHUNK_MAX_SIZE];
 
  char   sdummy[256];
 
  double LH_coincident_days = 0;
 
  cout.precision(10); 
 
  int k=0;
  while(true) {
    //   O2              22 1187312718 1187740818 3.5  days      Aug 21 00:00 - Aug 26 00:00
    //in >> note[k] >> chunk[k] >> start[k] >> stop[k] >> days[k] >> sdummy >> bmon[k] >> bday[k] >> sdummy >> emon[k] >> eday[k];
    in >> note[k] >> chunk[k] >> start[k] >> stop[k] >> days[k] >> sdummy >> bmon[k] >> bday[k] >> sdummy >> sdummy >> emon[k] >> eday[k] >> sdummy;
    if(!in.good()) break;
    if(note[k][0]=='#') continue;
    cout << "\t" << chunk[k] << "\t" << start[k] << "\t" << stop[k] << "\t" << days[k] << "\t" 
         << " days\t" << bmon[k] << " " << bday[k] << " - " << emon[k] << " " << eday[k] << "\t\t" << note[k] << endl;
    LH_coincident_days += days[k];
    k++;
  }
  in.close();
 
  cout << endl;
  cout << "LH_coincident_days = " << LH_coincident_days << endl;
  cout << endl;
 
//  exit(0);
 
  return k;
}