ProduceBadRunsByFailureMode.C

Go to the documentation of this file.

void ProduceBadRunsByFailureMode{

  TFile file("MetricsTree_All.root","read");
  TTree *tree = (TTree*)file.Get("mytree");

  int run, subrun;
  int pass_all, pass_hits,pass_db,pass_time,pass_slc,pass_trk,pass_runlen,pass_empty, ngooddb, ngoodcdb;

  tree->SetBranchAddress("run",&run);
  tree->SetBranchAddress("subrun",&subrun);
  tree->SetBranchAddress("ngooddb",&ngooddb);
  tree->SetBranchAddress("ngoodcdb",&ngoodcdb);
  tree->SetBranchAddress("pass_all",&pass_all);
  tree->SetBranchAddress("pass_hits",&pass_hits);
  tree->SetBranchAddress("pass_db",&pass_db);
  tree->SetBranchAddress("pass_time",&pass_time);
  tree->SetBranchAddress("pass_slc",&pass_slc);
  tree->SetBranchAddress("pass_trk",&pass_trk);
  tree->SetBranchAddress("pass_runlen",&pass_runlen);
  tree->SetBranchAddress("pass_empty",&pass_empty);

  ofstream procfile;
  procfile.open("FD_BRL_v0.txt");

  int nentries = tree->GetEntries();

  TString rowstr;

  ifstream fileListItr("ManBadSubruns.txt");

  int badrun, badsubrun;
  //int removedreason;
  set<int> cleanruns;
  set<int> alreadyused; 
  while(fileListItr >> badrun >> badsubrun){
    cleanruns.insert((100*badrun)+badsubrun);
    cout << "Manually added Run/Subrun: " << badrun <<"/"<< badsubrun<< endl;
  }

  for(int i = 0; i < nentries; i++){

    if(i%(nentries/10)==0) cout << int(100.0*i/nentries+0.5) << "% done" << endl;
    int removedreason = 0;
    tree->GetEntry(i);
    //putting in run numbers is optional, mostly used to conform to epoch periolds.
    if(cleanruns.count((run*100)+subrun)){
      removedreason = 8;
      rowstr = TString::Format("%d %d %d %d", run, subrun, removedreason, ngooddb);
      procfile << rowstr.Data() << endl;
      alreadyused.insert((run*100)+subrun);
    }
    if(!pass_all){
    if(!pass_hits){
      removedreason = 1; 
      rowstr = TString::Format("%d %d %d %d", run, subrun, removedreason, ngooddb);
      procfile << rowstr.Data() << endl;
    } 
    if(!pass_db || ngoodcdb < 4){
      removedreason = 2;
      rowstr = TString::Format("%d %d %d %d", run, subrun, removedreason, ngooddb);
      procfile << rowstr.Data() << endl;
    }
    if(!pass_time){
      removedreason = 3;
      rowstr = TString::Format("%d %d %d %d", run, subrun, removedreason, ngooddb);
      procfile << rowstr.Data() << endl;
    }
    if(!pass_slc){
      removedreason = 4;
      rowstr = TString::Format("%d %d %d %d", run, subrun, removedreason, ngooddb);
      procfile << rowstr.Data() << endl;
    }
    if(!pass_trk){
      removedreason = 5;
      rowstr = TString::Format("%d %d %d %d", run, subrun, removedreason, ngooddb);
      procfile << rowstr.Data() << endl;
    }
    if(!pass_runlen){
      removedreason = 6;
      rowstr = TString::Format("%d %d %d %d", run, subrun, removedreason, ngooddb);
      procfile << rowstr.Data() << endl;
    }
    if(!pass_empty){
      removedreason = 7;
      rowstr = TString::Format("%d %d %d %d", run, subrun, removedreason, ngooddb);
      procfile << rowstr.Data() << endl;
    }
    //if(cleanruns.count((run*100)+subrun)){
    //  removedreason = 9;
    //  rowstr = TString::Format("%d %d %d", run, subrun, removedreason);
    //  procfile << rowstr.Data() << endl;
   // }
    }
    if(pass_all && ngooddb < 14 && ngoodcdb > 3){
      removedreason = 0;
      rowstr = TString::Format("%d %d %d %d", run, subrun, removedreason, ngooddb);
      procfile << rowstr.Data() << endl;
    }
  }

   set<int>::iterator it; 
  for(it = cleanruns.begin(); it != cleanruns.end(); it++){
    int runsubrun = *it;
    if(!alreadyused.count(runsubrun)){
    std::cout <<"Manual run not in regular dataset... adding... " << std::endl;
    int temprun = runsubrun/100;
    int tempsubrun = runsubrun%100; 
    rowstr = TString::Format("%d %d %d %d", temprun, tempsubrun, 8, -1);
    procfile << rowstr.Data() << endl;
    }
  }
  procfile.close();

}