modulediff.cc

Go to the documentation of this file.

#include <Riostream.h>
#include <TDirectory.h>
#include <TFile.h>
#include <TROOT.h>
#include <TStyle.h>
#include <TKey.h>
#include <TH1.h>
#include <TH2.h>
#include <TH2D.h>
#include <TCanvas.h>
#include <TGraph.h>
#include <TPaveStats.h>
#include <TText.h>
#include <TLegend.h>
#include <cstring>
#include <utility>
#include <vector>
#include <sstream>
#include <algorithm>
#include <TString.h>
#include <TColor.h>
#include <dirent.h>
 
using namespace std;
 
int read_badmodlist(int, string, vector<int>&);
void get_difference(vector<int>, vector<int>, vector<int>&);
int get_filename(int, string, string&);
int search_closest_run(int, string);
void modulediff(int run_2, string repro_run2);
 
int main(int argc, char* argv[]) {
  if (argc == 3) {
    char* crun2 = argv[1];
    char* repro_run2 = argv[2];
 
    int run2 = 0;
    sscanf(crun2, "%d", &run2);
 
    std::cout << "ready to run modulediff " << run2 << " repro run " << repro_run2 << std::endl;
 
    modulediff(run2, repro_run2);
 
  } else {
    std::cout << "Too few arguments: " << argc << std::endl;
    return -1;
  }
  return 0;
}
 
void modulediff(int run_2, string repro_run2) {
  vector<int> badmodlist_run1;
  vector<int> badmodlist_run2;
  vector<int> modules_recovered;
  vector<int> modules_malformed;
  int run_1;
  string repro_run1;
 
  int res = search_closest_run(run_2, repro_run2);
  if (res > 0) {
    run_1 = res;
    repro_run1 = repro_run2;
  } else {
    printf("closest run not found, exiting.\n");
    return;
  }
 
  int flag1 = read_badmodlist(run_1, repro_run1, badmodlist_run1);
  int flag2 = read_badmodlist(run_2, repro_run2, badmodlist_run2);
 
  if (flag1 < 0 || flag2 < 0) {
    cout << "Error: file not found." << endl;
    return;
  }
 
  get_difference(badmodlist_run1, badmodlist_run2, modules_recovered);
  get_difference(badmodlist_run2, badmodlist_run1, modules_malformed);
 
  //save into file
  std::ofstream outfile;
  string namefile = "modulediff_emailbody.txt";
  outfile.open(namefile.c_str());
 
  outfile << "Recovered modules in run " << run_2 << ": " << (int)modules_recovered.size() << endl;
  outfile << "New bad modules in run   " << run_2 << ": " << (int)modules_malformed.size() << endl;
  outfile << "Using reference run " << run_1 << endl << endl;
 
  outfile << "Recovered modules in run " << run_2 << ":" << endl;
  if (modules_recovered.empty())
    outfile << " -" << endl;
  for (unsigned int i = 0; i < modules_recovered.size(); i++)
    outfile << " " << modules_recovered[i] << endl;
 
  outfile << "New bad modules that appeared in run " << run_2 << ":" << endl;
  if (modules_malformed.empty())
    outfile << " -" << endl;
  for (unsigned int i = 0; i < modules_malformed.size(); i++)
    outfile << " " << modules_malformed[i] << endl;
 
  outfile.close();
 
  //create two flat files to run the locatemodule script on later
 
  if (!modules_recovered.empty()) {
    std::ofstream outfile_good;
    outfile_good.open("modulediff_good.txt");
    for (unsigned int i = 0; i < modules_recovered.size(); i++)
      outfile_good << " " << modules_recovered[i] << endl;
    outfile_good.close();
  }
 
  if (!modules_malformed.empty()) {
    std::ofstream outfile_bad;
    outfile_bad.open("modulediff_bad.txt");
    for (unsigned int i = 0; i < modules_malformed.size(); i++)
      outfile_bad << " " << modules_malformed[i] << endl;
    outfile_bad.close();
  }
}
 
void get_difference(vector<int> badlist1, vector<int> badlist2, vector<int>& difflist) {
  //check if any element of badlist1 is in badlist2
  for (unsigned int i1 = 0; i1 < badlist1.size(); i1++) {
    bool thisrecovered = true;
    for (unsigned int i2 = 0; i2 < badlist2.size(); i2++)
      if (badlist1[i1] == badlist2[i2]) {
        thisrecovered = false;
        break;
      }
    if (thisrecovered)
      difflist.push_back(badlist1[i1]);
  }
}
 
int read_badmodlist(int run, string repro_type, vector<int>& badlist) {
  string filename;
  int flag = get_filename(run, repro_type, filename);
 
  if (flag < 0) {
    cout << "reading problem" << endl;
    return -1;
  }
 
  vector<string> subdet;
  subdet.push_back("TIB");
  subdet.push_back("TID/MINUS");
  subdet.push_back("TID/PLUS");
  subdet.push_back("TOB");
  subdet.push_back("TEC/MINUS");
  subdet.push_back("TEC/PLUS");
 
  string nrun = filename.substr(filename.find("_R000") + 5, 6);
 
  TFile* dqmfile = TFile::Open(filename.c_str(), "READ");
  string topdir = "DQMData/Run " + nrun + "/SiStrip/Run summary/MechanicalView";
  gDirectory->cd(topdir.c_str());
  TDirectory* mec1 = gDirectory;
 
  for (unsigned int i = 0; i < subdet.size(); i++) {
    string badmodule_dir = subdet[i] + "/BadModuleList";
    if (gDirectory->cd(badmodule_dir.c_str())) {
      TIter next(gDirectory->GetListOfKeys());
      TKey* key;
 
      while ((key = dynamic_cast<TKey*>(next()))) {
        string sflag = key->GetName();
        if (sflag.empty())
          continue;
 
        string detid = sflag.substr(sflag.find('<') + 1, 9);
        badlist.push_back(atoi(detid.c_str()));
      }
    } else {
      //cout << "no dir " << badmodule_dir << " in filename " << filename << endl;
    }
    mec1->cd();
  }
  dqmfile->Close();
  return 0;
}
 
int get_filename(int run, string repro_type, string& filename) {
  stringstream runstr;
  runstr << run;
 
  stringstream rundirprefix;
  rundirprefix << "000" << run / 100 << "xx/";
 
  stringstream thisdir;
  //thisdir << "/afs/cern.ch/cms/CAF/CMSCOMM/COMM_DQM/data/OfflineData/Run2011/" << repro_type.c_str() << "/" << rundirprefix.str();
  thisdir << "/afs/cern.ch/cms/CAF/CMSCOMM/COMM_DQM/data/OfflineData/HIRun2011/" << repro_type.c_str() << "/"
          << rundirprefix.str();
 
  string thisdir2 = thisdir.str();
  DIR* dp;
 
  if ((dp = opendir(thisdir2.c_str())) == nullptr) {
    cout << "dir " << thisdir2.c_str() << " not found" << endl;
    return -1;
  }
 
  struct dirent* dirp;
 
  string dqmfile;
 
  while ((dirp = readdir(dp)) != nullptr) {
    string dirfile = string(dirp->d_name);
    if (dirfile.find("__DQM") != string::npos && dirfile.find(runstr.str()) != string::npos) {
      dqmfile = dirfile;
      break;
    }
  }
 
  closedir(dp);
 
  if (dqmfile.size() < 10) {
    //cout << "file " << dqmfile << " not found" << endl;
    return -1;
  }
 
  filename = thisdir.str() + dqmfile;
 
  return 0;
}
 
int search_closest_run(int thisrun, string repro_type) {
  string filename;
 
  for (int test_run = thisrun - 1; test_run > thisrun - 1000; test_run--) {
    int res = get_filename(test_run, repro_type, filename);
    if (res == 0)
      return test_run;
  }
 
  return -1;
}