#include <iostream>
#include <cstdio>
#include <TDirectory.h>
#include <TFile.h>
#include <TKey.h>
#include <TH1.h>
#include <Riostream.h>
#include <string>
#include <vector>
#include <set>
#include <fstream>
#include <sstream>
#include <cstdint>
#include <cstdlib>
#include "listbadmodule.h"

Functions
void	listbadmodule (std::string filename, std::string pclfilename)

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

Function Documentation

◆ listbadmodule()

void listbadmodule	(	std::string	filename,
		std::string	pclfilename
	)

Definition at line 38 of file listbadmodule.cc.

References gather_cfg::cout, debug, corrVsCorr::filename, RemoveAddSevLevel::flag, mps_fire::i, crabWrapper::key, mps_splice::line, GetRecoTauVFromDQM_MC_cff::next, timingPdfMaker::outfile, AlCaHLTBitMon_QueryRunRegistry::string, and testProducerWithPsetDescEmpty_cfi::x1.

Referenced by main().

                                                               {
   int debug = 1;
 
   // extract fully bad modules from PCLBadComponents txt file
 
   std::set<unsigned int> pclbadmods;
 
   std::ifstream pclfile(pclfilename);
 
   char line[400];
   unsigned int pcldetid;
   char sixapvs[] = "1 1 1 1 1 1";
   char fourapvs[] = "1 1 x x 1 1";
 
   while (pclfile.getline(line, 400)) {
     if (strstr(line, sixapvs) || strstr(line, fourapvs)) {
       std::stringstream linestream;
       linestream << line;
       linestream >> pcldetid;
       //      std::cout << pcldetid << std::endl;
       pclbadmods.insert(pcldetid);
     }
   }
 
   std::vector<std::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");
 
   std::string nrun = filename.substr(filename.find("_R000") + 5, 6);
   int fileNum = atoi(nrun.c_str());
   std::cout << " ------   Run " << fileNum << std::endl;
 
   std::ofstream outfile;
   std::string namefile;
   namefile = "QualityTestOBSOLETE_run" + nrun + ".txt";
   outfile.open(namefile.c_str());
 
   TFile *myfile = TFile::Open(filename.c_str());
   if (debug == 1) {
     std::cout << " Opened " << filename << std::endl;
   }
   std::string topdir = "DQMData/Run " + nrun + "/SiStrip/Run summary/MechanicalView";
   gDirectory->cd(topdir.c_str());
   TDirectory *mec1 = gDirectory;
 
   //get the summary first
   std::vector<int> nbadmod;
   for (unsigned int i = 0; i < subdet.size(); i++) {
     int nbad = 0;
     std::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()))) {
         std::string sflag = key->GetName();
         if (sflag.empty())
           continue;
         nbad++;
       }
     }
     nbadmod.push_back(nbad);
     mec1->cd();
   }
 
   outfile << "Number of bad modules in total excluding PCL-only bad modules:" << std::endl;
   outfile << "--------------------------------------------------------------" << std::endl;
   outfile << subdet.at(0) << ": " << nbadmod.at(0) << std::endl;
   outfile << subdet.at(1) << ": " << nbadmod.at(1) << std::endl;
   outfile << subdet.at(2) << ": " << nbadmod.at(2) << std::endl;
   outfile << subdet.at(3) << ": " << nbadmod.at(3) << std::endl;
   outfile << subdet.at(4) << ": " << nbadmod.at(4) << std::endl;
   outfile << subdet.at(5) << ": " << nbadmod.at(5) << std::endl;
   outfile << "-------------------------------" << std::endl;
 
   outfile << std::endl << "List of bad modules per partition:" << std::endl;
   outfile << "----------------------------------" << std::endl;
 
   std::set<unsigned int>::const_iterator pclbadmod = pclbadmods.begin();
 
   for (unsigned int i = 0; i < subdet.size(); i++) {
     std::string badmodule_dir = subdet[i] + "/BadModuleList";
     outfile << " " << std::endl;
     outfile << "SubDetector " << subdet[i] << std::endl;
     outfile << " " << std::endl;
     std::cout << badmodule_dir.c_str() << std::endl;
     if (gDirectory->cd(badmodule_dir.c_str())) {
       //
       // Loop to find bad module for each partition
       //
       TIter next(gDirectory->GetListOfKeys());
       TKey *key;
 
       while ((key = dynamic_cast<TKey *>(next()))) {
         std::string sflag = key->GetName();
         if (sflag.empty())
           continue;
         std::string detid = sflag.substr(sflag.find('<') + 1, 9);
         size_t pos1 = sflag.find('/');
         sflag = sflag.substr(sflag.find('<') + 13, pos1 - 2);
         int flag = atoi(sflag.c_str());
         sscanf(detid.c_str(), "%u", &pcldetid);
         // the following loop add modules which are bad only for the PCL to the list.
         // It requires that the bad modules from QT are sorted as the one of the PCL
         while (pclbadmod != pclbadmods.end() && pcldetid > *pclbadmod) {
           outfile << "Module " << *pclbadmod << " PCLBadModule " << std::endl;
           pclbadmod++;
         }
         std::string message;
         message = "Module " + detid;
         if (((flag >> 0) & 0x1) > 0)
           message += " Fed BadChannel : ";
         if (((flag >> 1) & 0x1) > 0)
           message += " # of Digi : ";
         if (((flag >> 2) & 0x1) > 0)
           message += " # of Clusters :";
         if (((flag >> 3) & 0x1) > 0)
           message += " Excluded FED Channel ";
         if (((flag >> 4) & 0x1) > 0)
           message += " DCSError ";
         if (pclbadmods.find(pcldetid) != pclbadmods.end()) {
           message += " PCLBadModule ";
           pclbadmod = pclbadmods.find(pcldetid);
           pclbadmod++;
         }
         outfile << message.c_str() << std::endl;
       }
     }
     mec1->cd();
   }
   myfile->Close();
   outfile.close();
 }

◆ 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:

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 22 of file listbadmodule.cc.

References dir2webdir::argc, GCPpyPlots::argv, gather_cfg::cout, corrVsCorr::filename, and listbadmodule().

                                  {
   if (argc == 3) {
     char *filename = argv[1];
     char *pclfilename = argv[2];
 
     std::cout << "ready to prepare list of bad modules " << filename << std::endl;
 
     listbadmodule(filename, pclfilename);
 
   } else {
     std::cout << "Too few arguments: " << argc << std::endl;
     return -1;
   }
   return 0;
 }

Functions

Function Documentation

◆ listbadmodule()

◆ main()