CreateChunkCuts.C

Go to the documentation of this file.

#define CHUNK_FILE_LIST         "Chunk_List.txt"
#define MACRO_READ_CHUNK        "ReadChunkList.C"
#define CHUNK_MAX_SIZE          100

#define OFILE_CHUNK_CUTS   "Chunks_Cuts.hh"

//#define EXECUTE

void CreateChunkCuts(TString run) {

  // 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",cwb_config_env,run.Data(),CHUNK_FILE_LIST);
  cout << chunk_file_list << endl;

  char macro_read_chunk_path[1024];
  sprintf(macro_read_chunk_path,"%s/MACROS/%s",cwb_config_env,MACRO_READ_CHUNK);

  CWB::Toolbox::checkFile(macro_read_chunk_path);

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

  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 tcut[256];

  char ofile[1024];
  sprintf(ofile,"%s/%s/CHUNKS/%s",cwb_config_env,run.Data(),OFILE_CHUNK_CUTS);
  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 O2_K%02d_cut(\"O2_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

  }

  // create cut chunk99 : full period
  sprintf(tcut,"TCut O2_K%02d_cut(\"O2_K%02d_cut\",\"time[0]>%d && time[0]<=%d\");",99,99,start[0],stop[nChunks-1]);
#ifdef EXECUTE
  out << tcut << endl;
#else
  cout << tcut << endl;
#endif


#ifdef EXECUTE
  out.close();
#endif

  exit(0);

}