esetIFAR.C

Go to the documentation of this file.

{
  // Definitions of the search bins used for the GW150914 analysis
  #include "../../cwb/postproduction/O1/GW150914_search_bins.hh"

  // Merge Label
  TString merge_label = "M1.V_cat2LH_hvetoLH_cat3LH";

  // Definition of the FAR file names
  TString far_dir          = "../ER8b_12Sep20Oct_C01_BKG_LF_rMRA_run1ato10a/report/dump/";
  TString far_eunmodeled   = far_dir+"FARvsRHO_ER8b_12Sep20Oct_C01_BKG_LF_rMRA_run1ato10a_eunmodeled.txt";
  TString far_econstrained = far_dir+"FARvsRHO_ER8b_12Sep20Oct_C01_BKG_LF_rMRA_run1ato10a_econstrained.txt";
  TString far_chirp        = far_dir+"FARvsRHO_ER8b_12Sep20Oct_C01_BKG_LF_rMRA_run1ato10a_chirp.txt";

  char cmd[1024];

  // add ifar parameter to the wave root file for the eunmodeled bin
  sprintf(cmd,"${CWB_SCRIPTS}/cwb_setifar.csh %s \"%s\" %s eunmodeled exclusive",
          merge_label.Data(),eunmodeled.GetTitle(),far_eunmodeled.Data());
  cout << cmd << endl;
  gSystem->Exec(cmd);

  // add ifar parameter to the wave root file for the econstrained bin
  sprintf(cmd,"${CWB_SCRIPTS}/cwb_setifar.csh %s.S_eunmodeled \"%s\" %s econstrained exclusive",
          merge_label.Data(),econstrained.GetTitle(),far_econstrained.Data());
  cout << cmd << endl;
  gSystem->Exec(cmd);

  // add ifar parameter to the wave root file for the chirp bin
  sprintf(cmd,"${CWB_SCRIPTS}/cwb_setifar.csh %s.S_eunmodeled.S_econstrained \"%s\" %s chirp exclusive",
          merge_label.Data(),chirp.GetTitle(),far_chirp.Data());
  cout << cmd << endl;
  gSystem->Exec(cmd);
  
  exit(0);
}