CreateChunkCuts.C

Go to the documentation of this file.

#define MACRO_READ_CHUNK        "/home/waveburst/O1/MACROS/ReadChunkList.C"

#define CHUNK_MAX_SIZE  100

#define EXECUTE

{

  CWB::Toolbox::checkFile(MACRO_READ_CHUNK);

  // load macro
  gROOT->LoadMacro(gSystem->ExpandPathName(MACRO_READ_CHUNK));

  int    chunk[CHUNK_MAX_SIZE];
  double start[CHUNK_MAX_SIZE];
  double stop[CHUNK_MAX_SIZE];

  int nChunks = ReadChunkList(chunk,start,stop);

  char tcut[256];

  char ofile[1024];
  sprintf(ofile,"/home/waveburst/O1/CHUNKS/O1_Chunks_Cuts.hh");
  cout << "output file list : " << ofile << endl;
#ifdef EXECUTE
  ofstream out;
  out.open(ofile,ios::out);
#endif
  for(int k=0;k<nChunks;k++) {
//    cout << "\t" << chunk[k] << "\t" << (int)start[k] << "\t" << (int)stop[k] << endl;

    sprintf(tcut,"TCut O1_K%02d_cut(\"O1_K%02d_cut\",\"time[0]>%d && time[0]<=%d\");",chunk[k],chunk[k],start[k],stop[k]);
    
#ifdef EXECUTE
    out << tcut << endl;
#else
    cout << tcut << endl;
#endif

  }
#ifdef EXECUTE
  out.close();
#endif

  exit(0);

}