CreateChunkFiles.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/>.
*/
 
 
#include "ReadChunkList.C"
 
#define CHUNK_FILE_LIST         "Chunk_List.txt"
#define CHUNK_MAX_SIZE          100
 
#define EXECUTE
 
void CreateChunkFiles(TString run, TString search="") {
 
  // get CWB_CONFIG
  char cwb_config_env[1024] = "";
  if(gSystem->Getenv("CWB_CONFIG")!=NULL) {
    strcpy(cwb_config_env,TString(gSystem->Getenv("CWB_CONFIG")).Data());
  }
 
  char chunk_file_list[1024];
  sprintf(chunk_file_list,"%s/%s/CHUNKS/%s/%s",cwb_config_env,run.Data(),search.Data(),CHUNK_FILE_LIST);
  cout << chunk_file_list << endl;
 
  int    chunk[CHUNK_MAX_SIZE];
  double start[CHUNK_MAX_SIZE];
  double stop[CHUNK_MAX_SIZE];
 
  int nChunks = ReadChunkList(chunk_file_list, chunk, start, stop);
 
  char ofile[1024];
  for(int k=0;k<nChunks;k++) {
//    cout << "\t" << chunk[k] << "\t" << (int)start[k] << "\t" << (int)stop[k] << endl;
 
    sprintf(ofile,"%s/%s/CHUNKS/%s/K%02d.period",cwb_config_env,run.Data(),search.Data(),chunk[k]);
 
#ifdef EXECUTE
    ofstream out;
    out.open(ofile,ios::out);
    out.precision(16);
    out << start[k] << "\t" << stop[k] << endl;
    out.close();
#else
    cout.precision(16);
    cout << start[k] << "\t" << stop[k] << endl;
#endif
 
    cout << "output file list : " << ofile << endl;
  }
 
  // create chunk99 : full period
  sprintf(ofile,"%s/%s/CHUNKS/%s/K%02d.period",cwb_config_env,run.Data(),search.Data(),99);
#ifdef EXECUTE
  ofstream out;
  out.open(ofile,ios::out);
  out.precision(16);
  out << 0 << "\t" << 2147483647 << endl;
  out.close();
#else
  cout.precision(16);
  cout << 0 << "\t" << 2147483647 << endl;
#endif
  cout << "output file list : " << ofile << endl;
 
  exit(0);
 
}