#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>

Functions
void	get_difference (vector< int >, vector< int >, vector< int > &)

int	get_filename (int, string, string &)

int	main (int argc, char *argv[])

void	modulediff (int run_2, string repro_run2)

int	read_badmodlist (int, string, vector< int > &)

int	search_closest_run (int, string)

Function Documentation

◆ get_difference()

void get_difference	(	vector< int >	badlist1,
		vector< int >	badlist2,
		vector< int > &	difflist
	)

Definition at line 121 of file modulediff.cc.

References testProducerWithPsetDescEmpty_cfi::i1, and testProducerWithPsetDescEmpty_cfi::i2.

Referenced by modulediff().

                                                                                        {
   //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]);
   }
 }

◆ get_filename()

int get_filename	(	int	run,
		string	repro_type,
		string &	filename
	)

Definition at line 182 of file modulediff.cc.

References gather_cfg::cout, dqmPostProcessing_online::DIR, Calorimetry_cff::dp, fileCollector::dqmfile, corrVsCorr::filename, writedatasetfile::run, fileinputsource_cfi::runstr, AlCaHLTBitMon_QueryRunRegistry::string, and compareTotals::thisdir.

Referenced by read_badmodlist(), and search_closest_run().

                                                                {
   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;
 }

◆ main()

int main	(	int	argc,
		char *	argv[]
	)

===============================================================================================================================================================================================

variant2: for each run define phi-averaged A for normalization channel (Dref,16) and then, divide Rijk on it, i.e. get RRijk

eta=27

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=19

eta=17

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=20

eta=19

eta=18

eta=27 L1=1

eta=25 L1=1

eta=23 L1=1

eta=22 L1=1

eta=21 L1=1

eta=29 L1=1

eta=26 L1=1

eta=24 L1=1

eta=20 L1=1

eta=19 L1=1

eta=18 L1=1

eta=17 L1=1

eta=28 L7=1

eta=27 L7=1

eta=25 L7=1

eta=23 L7=1

eta=22 L7=1

eta=21 L7=1

eta=26 L7=1

eta=24 L7=1

eta=20 L7=1

eta=19 L7=1

eta=18 L7=1

eta=17 L7=1

eta=27

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=19

eta=17

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=20

eta=19

eta=18

eta=27 L1=1

eta=25 L1=1

eta=23 L1=1

eta=22 L1=1

eta=21 L1=1

eta=26 L1=1

eta=24 L1=1

eta=20 L1=1

eta=19 L1=1

eta=18 L1=1

eta=17 L1=1

eta=28 L7=1

eta=27 L7=1

eta=25 L7=1

eta=23 L7=1

eta=22 L7=1

eta=21 L7=1

eta=26 L7=1

eta=24 L7=1

eta=20 L7=1

eta=19 L7=1

eta=18 L7=1

eta=17 L7=1

eta=27

eta=28

errA with average Amplitudes

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Sticking ADC counts in neughbouring TS for HB:

Sticking ADC counts in neughbouring TS for HE:

RBX:

errA with average Amplitudes

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

RBX:

errA with average Amplitudes

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

RBX:

Prepare maps of good/bad channels:

Definition at line 32 of file modulediff.cc.

References dir2webdir::argc, GCPpyPlots::argv, gather_cfg::cout, modulediff(), and compare_using_db::run2.

                                  {
   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;
 }

◆ modulediff()

void modulediff	(	int	run_2,
		string	repro_run2
	)

Definition at line 51 of file modulediff.cc.

References gather_cfg::cout, get_difference(), mps_fire::i, createfilelist::int, PixelBaryCentreAnalyzer_cfg::outfile, read_badmodlist(), and search_closest_run().

Referenced by main().

                                               {
   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();
   }
 }

◆ read_badmodlist()

int read_badmodlist	(	int	run,
		string	repro_type,
		vector< int > &	badlist
	)

Definition at line 135 of file modulediff.cc.

References gather_cfg::cout, ALCARECOPPSCalTrackBasedSel_cff::detid, fileCollector::dqmfile, corrVsCorr::filename, RemoveAddSevLevel::flag, get_filename(), mps_fire::i, submitPVResolutionJobs::key, GetRecoTauVFromDQM_MC_cff::next, and writedatasetfile::run.

Referenced by modulediff().

                                                                       {
   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;
 }

◆ search_closest_run()

int search_closest_run	(	int	thisrun,
		string	repro_type
	)

Definition at line 226 of file modulediff.cc.

References corrVsCorr::filename, and get_filename().

Referenced by modulediff().

                                                        {
   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;
 }

Functions

Function Documentation

◆ get_difference()

◆ get_filename()

◆ main()

◆ modulediff()

◆ read_badmodlist()

◆ search_closest_run()