|
|
|
#include <fstream>#include <iostream>#include <vector>#include <map>#include <string>#include <sstream>#include "TFile.h"#include "TVectorT.h"#include "TGraph.h"#include "TGaxis.h"#include "TCanvas.h"#include "TH1F.h"Go to the source code of this file.
Functions | |
| int | main (int argc, char *argv[]) |
| void | makeTKTrend (const char *inFileName, const char *outFileName, std::string subDetName, std::string partName, const unsigned int partNumber) |
| int main | ( | int | argc, |
| char * | argv[] | ||
| ) |
Definition at line 18 of file makeTKTrend.cc.
References gather_cfg::cout, makeTKTrend(), and partName.
{
if(argc==6) {
char* inFileName = argv[1];
char* outFileName = argv[2];
char* subDetName = argv[3];
char* partName = argv[4];
int partNumber = atoi(argv[5]);
std::cout << "ready to make trend plots from "
<< inFileName << " to " << outFileName << " for "
<< subDetName << " " << partName << " " << partNumber << std::endl;
int returncode = 0;
makeTKTrend(inFileName,outFileName,subDetName,partName,partNumber);
return returncode;
}
else {std::cout << "Too few arguments: " << argc << std::endl; return -1; }
return -9;
}
| void makeTKTrend | ( | const char * | inFileName, |
| const char * | outFileName, | ||
| std::string | subDetName, | ||
| std::string | partName, | ||
| const unsigned int | partNumber | ||
| ) |
Definition at line 44 of file makeTKTrend.cc.
References j, GetRecoTauVFromDQM_MC_cff::outFile, partName, and makeHLTPrescaleTable::values.
Referenced by main().
{
// Maps <Run number, nBad>
std::map<unsigned int, unsigned int> badModulesTK;
std::map<unsigned int, unsigned int> badFibersTK;
std::map<unsigned int, unsigned int> badAPVsTK;
std::map<unsigned int, unsigned int> badStripsTK;
std::map<unsigned int, unsigned int> badStripsFromAPVsTK;
std::map<unsigned int, unsigned int> allBadStripsTK;
std::ostringstream oss;
if(partName=="All") partName="";
// // Number of modules, fibers, APVs, strips for each tracker part
// std::vector< std::string > tkParts;
// tkParts.push_back("Tracker");
// tkParts.push_back("TIB");
// tkParts.push_back("TID");
// tkParts.push_back("TOB");
// tkParts.push_back("TEC");
// tkParts.push_back("TIB Layer 1");
// tkParts.push_back("TIB Layer 2");
// tkParts.push_back("TIB Layer 3");
// tkParts.push_back("TIB Layer 4");
// tkParts.push_back("TID+ Disk 1");
// tkParts.push_back("TID+ Disk 2");
// tkParts.push_back("TID+ Disk 3");
// tkParts.push_back("TID- Disk 1");
// tkParts.push_back("TID- Disk 2");
// tkParts.push_back("TID- Disk 3");
// tkParts.push_back("TOB Layer 1");
// tkParts.push_back("TOB Layer 2");
// tkParts.push_back("TOB Layer 3");
// tkParts.push_back("TOB Layer 4");
// tkParts.push_back("TOB Layer 5");
// tkParts.push_back("TOB Layer 6");
// tkParts.push_back("TEC+ Disk 1");
// tkParts.push_back("TEC+ Disk 2");
// tkParts.push_back("TEC+ Disk 3");
// tkParts.push_back("TEC+ Disk 4");
// tkParts.push_back("TEC+ Disk 5");
// tkParts.push_back("TEC+ Disk 6");
// tkParts.push_back("TEC+ Disk 7");
// tkParts.push_back("TEC+ Disk 8");
// tkParts.push_back("TEC+ Disk 9");
// tkParts.push_back("TEC- Disk 1");
// tkParts.push_back("TEC- Disk 2");
// tkParts.push_back("TEC- Disk 3");
// tkParts.push_back("TEC- Disk 4");
// tkParts.push_back("TEC- Disk 5");
// tkParts.push_back("TEC- Disk 6");
// tkParts.push_back("TEC- Disk 7");
// tkParts.push_back("TEC- Disk 8");
// tkParts.push_back("TEC- Disk 9");
//
// std::vector<unsigned int> nModules;
// nModules.push_back(15148);
// nModules.push_back(2724);
// nModules.push_back(816);
// nModules.push_back(5208);
// nModules.push_back(6400);
// nModules.push_back(672);
// nModules.push_back(864);
// nModules.push_back(540);
// nModules.push_back(648);
// nModules.push_back(136);
// nModules.push_back(136);
// nModules.push_back(136);
// nModules.push_back(136);
// nModules.push_back(136);
// nModules.push_back(136);
// nModules.push_back(1008);
// nModules.push_back(1152);
// nModules.push_back(648);
// nModules.push_back(720);
// nModules.push_back(792);
// nModules.push_back(888);
// nModules.push_back(408);
// nModules.push_back(408);
// nModules.push_back(408);
// nModules.push_back(360);
// nModules.push_back(360);
// nModules.push_back(360);
// nModules.push_back(312);
// nModules.push_back(312);
// nModules.push_back(272);
// nModules.push_back(408);
// nModules.push_back(408);
// nModules.push_back(408);
// nModules.push_back(360);
// nModules.push_back(360);
// nModules.push_back(360);
// nModules.push_back(312);
// nModules.push_back(312);
// nModules.push_back(272);
//
// std::vector<unsigned int> nFibers;
// nFibers.push_back(36392);
// nFibers.push_back(6984);
// nFibers.push_back(2208);
// nFibers.push_back(12096);
// nFibers.push_back(15104);
// nFibers.push_back(2016);
// nFibers.push_back(2592);
// nFibers.push_back(1080);
// nFibers.push_back(1296);
// nFibers.push_back(368);
// nFibers.push_back(368);
// nFibers.push_back(368);
// nFibers.push_back(368);
// nFibers.push_back(368);
// nFibers.push_back(368);
// nFibers.push_back(2016);
// nFibers.push_back(2304);
// nFibers.push_back(1296);
// nFibers.push_back(1440);
// nFibers.push_back(2376);
// nFibers.push_back(2664);
// nFibers.push_back(992);
// nFibers.push_back(992);
// nFibers.push_back(992);
// nFibers.push_back(848);
// nFibers.push_back(848);
// nFibers.push_back(848);
// nFibers.push_back(704);
// nFibers.push_back(704);
// nFibers.push_back(624);
// nFibers.push_back(992);
// nFibers.push_back(992);
// nFibers.push_back(992);
// nFibers.push_back(848);
// nFibers.push_back(848);
// nFibers.push_back(848);
// nFibers.push_back(704);
// nFibers.push_back(704);
// nFibers.push_back(624);
//
// std::vector<unsigned int> nAPVs;
// nAPVs.push_back(72784);
// nAPVs.push_back(13968);
// nAPVs.push_back(4416);
// nAPVs.push_back(24192);
// nAPVs.push_back(30208);
// nAPVs.push_back(4032);
// nAPVs.push_back(5184);
// nAPVs.push_back(2160);
// nAPVs.push_back(2592);
// nAPVs.push_back(736);
// nAPVs.push_back(736);
// nAPVs.push_back(736);
// nAPVs.push_back(736);
// nAPVs.push_back(736);
// nAPVs.push_back(736);
// nAPVs.push_back(4032);
// nAPVs.push_back(4608);
// nAPVs.push_back(2592);
// nAPVs.push_back(2880);
// nAPVs.push_back(4752);
// nAPVs.push_back(5328);
// nAPVs.push_back(1984);
// nAPVs.push_back(1984);
// nAPVs.push_back(1984);
// nAPVs.push_back(1696);
// nAPVs.push_back(1696);
// nAPVs.push_back(1696);
// nAPVs.push_back(1408);
// nAPVs.push_back(1408);
// nAPVs.push_back(1248);
// nAPVs.push_back(1984);
// nAPVs.push_back(1984);
// nAPVs.push_back(1984);
// nAPVs.push_back(1696);
// nAPVs.push_back(1696);
// nAPVs.push_back(1696);
// nAPVs.push_back(1408);
// nAPVs.push_back(1408);
// nAPVs.push_back(1248);
//
// std::vector<unsigned int> nStrips;
// nStrips.push_back(9316352);
// nStrips.push_back(1787904);
// nStrips.push_back(565248);
// nStrips.push_back(3096576);
// nStrips.push_back(3866624);
// nStrips.push_back(516096);
// nStrips.push_back(663552);
// nStrips.push_back(276480);
// nStrips.push_back(331776);
// nStrips.push_back(94208);
// nStrips.push_back(94208);
// nStrips.push_back(94208);
// nStrips.push_back(94208);
// nStrips.push_back(94208);
// nStrips.push_back(94208);
// nStrips.push_back(516096);
// nStrips.push_back(589824);
// nStrips.push_back(331776);
// nStrips.push_back(368640);
// nStrips.push_back(608256);
// nStrips.push_back(681984);
// nStrips.push_back(253952);
// nStrips.push_back(253952);
// nStrips.push_back(253952);
// nStrips.push_back(217088);
// nStrips.push_back(217088);
// nStrips.push_back(217088);
// nStrips.push_back(180224);
// nStrips.push_back(180224);
// nStrips.push_back(159744);
// nStrips.push_back(253952);
// nStrips.push_back(253952);
// nStrips.push_back(253952);
// nStrips.push_back(217088);
// nStrips.push_back(217088);
// nStrips.push_back(217088);
// nStrips.push_back(180224);
// nStrips.push_back(180224);
// nStrips.push_back(159744);
//
// // Map with <name of tracker part, count of channels in the part>
// std::map<std::string, unsigned int> allModulesTK;
// std::map<std::string, unsigned int> allFibersTK;
// std::map<std::string, unsigned int> allAPVsTK;
// std::map<std::string, unsigned int> allStripsTK;
// for(unsigned int i = 0; i < tkParts.size(); i++)
// {
// allModulesTK[tkParts[i].c_str()] = nModules[i];
// allFibersTK[tkParts[i].c_str()] = nFibers[i];
// allAPVsTK[tkParts[i].c_str()] = nAPVs[i];
// allStripsTK[tkParts[i].c_str()] = nStrips[i];
// }
//
// // for(std::map< std::string, unsigned int>::iterator it = allStripsTK.begin(); it != allStripsTK.end(); it++)
// // {
// // std::cout << it->first.c_str() << " " << it->second << std::endl;
// // }
//
// // Finds number of channels from above map
// std::string completePartName = subDetName;
// if(partName.compare("") != 0)
// completePartName += " " + partName;
// if(partNumber != 0)
// {
// oss.str("");
// oss << partNumber;
// completePartName += " " + oss.str();
// }
// // Total number of channels in currently processed map
// const unsigned int nModulesInPart = allModulesTK[completePartName.c_str()];
// const unsigned int nFibersInPart = allFibersTK[completePartName.c_str()];
// const unsigned int nAPVsInPart = allAPVsTK[completePartName.c_str()];
// const unsigned int nStripsInPart = allStripsTK[completePartName.c_str()];
// Read input file
std::ifstream resultsFile(inFileName);
unsigned int runIOV;
unsigned int values[6];
do
{
resultsFile >> runIOV;
resultsFile >> values[0] >> values[1] >> values[2] >> values[3] >> values[4] >> values[5];
// std::cout << runIOV << " " << values[0] << " " << values[1] << " " << values[2] << " " << values[3] << " " << values[4] << std::endl;
badModulesTK[runIOV]=values[0];
badFibersTK[runIOV]=values[1];
badAPVsTK[runIOV]=values[2];
badStripsTK[runIOV]=values[3];
badStripsFromAPVsTK[runIOV]=values[4];
allBadStripsTK[runIOV]=values[5];
}
while(!resultsFile.eof());
const unsigned int IOVSize = badStripsTK.size();
// Create histograms
std::string histoName;
std::string histoTitle;
oss.str("");
histoName = "hBadModules" + subDetName + partName;
if(partNumber!=0)
{
oss << partNumber;
histoName += oss.str();
}
oss.str("");
histoTitle = "Bad modules in " + subDetName;
if(partName!="")
{
histoTitle += " " + partName;
}
if(partNumber!=0)
{
oss << partNumber;
histoTitle += " " + oss.str();
}
TH1F* hBadModulesTK = new TH1F(histoName.c_str(), histoTitle.c_str(), IOVSize, 0.5, IOVSize+0.5);
oss.str("");
histoName = "hBadFibers" + subDetName + partName;
if(partNumber!=0)
{
oss << partNumber;
histoName += oss.str();
}
oss.str("");
histoTitle = "Bad fibers in " + subDetName;
if(partName!="")
{
histoTitle += " " + partName;
}
if(partNumber!=0)
{
oss << partNumber;
histoTitle += " " + oss.str();
}
TH1F* hBadFibersTK = new TH1F(histoName.c_str(), histoTitle.c_str(), IOVSize, 0.5, IOVSize+0.5);
oss.str("");
histoName = "hBadAPVs" + subDetName + partName;
if(partNumber!=0)
{
oss << partNumber;
histoName += oss.str();
}
oss.str("");
histoTitle = "Bad APVs in " + subDetName;
if(partName!="")
{
histoTitle += " " + partName;
}
if(partNumber!=0)
{
oss << partNumber;
histoTitle += " " + oss.str();
}
TH1F* hBadAPVsTK = new TH1F(histoName.c_str(), histoTitle.c_str(), IOVSize, 0.5, IOVSize+0.5);
oss.str("");
histoName = "hBadStrips" + subDetName + partName;
if(partNumber!=0)
{
oss << partNumber;
histoName += oss.str();
}
oss.str("");
histoTitle = "Bad strips in " + subDetName;
if(partName!="")
{
histoTitle += " " + partName;
}
if(partNumber!=0)
{
oss << partNumber;
histoTitle += " " + oss.str();
}
TH1F* hBadStripsTK = new TH1F(histoName.c_str(), histoTitle.c_str(), IOVSize, 0.5, IOVSize+0.5);
oss.str("");
histoName = "hBadStripsFromAPVs" + subDetName + partName;
if(partNumber!=0)
{
oss << partNumber;
histoName += oss.str();
}
oss.str("");
histoTitle = "Bad strips from APVs in " + subDetName;
if(partName!="")
{
histoTitle += " " + partName;
}
if(partNumber!=0)
{
oss << partNumber;
histoTitle += " " + oss.str();
}
TH1F* hBadStripsFromAPVsTK = new TH1F(histoName.c_str(), histoTitle.c_str(), IOVSize, 0.5, IOVSize+0.5);
oss.str("");
histoName = "hAllBadStrips" + subDetName + partName;
if(partNumber!=0)
{
oss << partNumber;
histoName += oss.str();
}
oss.str("");
histoTitle = "All bad strips in " + subDetName;
if(partName!="")
{
histoTitle += " " + partName;
}
if(partNumber!=0)
{
oss << partNumber;
histoTitle += " " + oss.str();
}
TH1F* hAllBadStripsTK = new TH1F(histoName.c_str(), histoTitle.c_str(), IOVSize, 0.5, IOVSize+0.5);
unsigned int j = 0;
for(std::map<unsigned int, unsigned int>::iterator iMap=badModulesTK.begin(); iMap!=badModulesTK.end(); iMap++ )
{
hBadModulesTK->SetBinContent(++j,/*(double)*/iMap->second/*/(double)nModulesInPart*/);
oss.str("");
oss << iMap->first;
hBadModulesTK->GetXaxis()->SetBinLabel(j,oss.str().c_str());
// std::cout << hBadModulesTK->GetBinContent(j) << std::endl;
}
TCanvas* cBadModulesTK = new TCanvas();
hBadModulesTK->Draw();
cBadModulesTK->Update();
j = 0;
for(std::map<unsigned int, unsigned int>::iterator iMap=badFibersTK.begin(); iMap!=badFibersTK.end(); iMap++ )
{
hBadFibersTK->SetBinContent(++j,/*(double)*/iMap->second/*/(double)nFibersInPart*/);
oss.str("");
oss << iMap->first;
hBadFibersTK->GetXaxis()->SetBinLabel(j,oss.str().c_str());
// std::cout << hBadFibersTK->GetBinContent(j) << std::endl;
}
TCanvas* cBadFibersTK = new TCanvas();
hBadFibersTK->Draw();
cBadFibersTK->Update();
j = 0;
for(std::map<unsigned int, unsigned int>::iterator iMap=badAPVsTK.begin(); iMap!=badAPVsTK.end(); iMap++ )
{
hBadAPVsTK->SetBinContent(++j,/*(double)*/iMap->second/*/(double)nAPVsInPart*/);
oss.str("");
oss << iMap->first;
hBadAPVsTK->GetXaxis()->SetBinLabel(j,oss.str().c_str());
// std::cout << hBadAPVsTK->GetBinContent(j) << std::endl;
}
TCanvas* cBadAPVsTK = new TCanvas();
hBadAPVsTK->Draw();
cBadAPVsTK->Update();
j = 0;
for(std::map<unsigned int, unsigned int>::iterator iMap=badStripsTK.begin(); iMap!=badStripsTK.end(); iMap++ )
{
hBadStripsTK->SetBinContent(++j,/*(double)*/iMap->second/*/(double)nStripsInPart*/);
oss.str("");
oss << iMap->first;
hBadStripsTK->GetXaxis()->SetBinLabel(j,oss.str().c_str());
// std::cout << hBadStripsTK->GetBinContent(j) << std::endl;
}
TCanvas* cBadStripsTK = new TCanvas();
hBadStripsTK->Draw();
cBadStripsTK->Update();
j = 0;
for(std::map<unsigned int, unsigned int>::iterator iMap=badStripsFromAPVsTK.begin(); iMap!=badStripsFromAPVsTK.end(); iMap++ )
{
hBadStripsFromAPVsTK->SetBinContent(++j,/*(double)*/iMap->second/*/(double)nStripsInPart*/);
oss.str("");
oss << iMap->first;
hBadStripsFromAPVsTK->GetXaxis()->SetBinLabel(j,oss.str().c_str());
// std::cout << hBadStripsTK->GetBinContent(j) << std::endl;
}
TCanvas* cBadStripsFromAPVsTK = new TCanvas();
hBadStripsFromAPVsTK->Draw();
cBadStripsFromAPVsTK->Update();
j = 0;
for(std::map<unsigned int, unsigned int>::iterator iMap=allBadStripsTK.begin(); iMap!=allBadStripsTK.end(); iMap++ )
{
hAllBadStripsTK->SetBinContent(++j,/*(double)*/iMap->second/*/(double)nStripsInPart*/);
oss.str("");
oss << iMap->first;
hAllBadStripsTK->GetXaxis()->SetBinLabel(j,oss.str().c_str());
// std::cout << hAllBadStripsTK->GetBinContent(j) << std::endl;
}
TCanvas* cAllBadStripsTK = new TCanvas();
hAllBadStripsTK->Draw();
cAllBadStripsTK->Update();
// Write histograms to output file
TFile* outFile = new TFile(outFileName, "UPDATE");
outFile->cd();
hBadModulesTK->Write();
hBadFibersTK->Write();
hBadAPVsTK->Write();
hBadStripsTK->Write();
hBadStripsFromAPVsTK->Write();
hAllBadStripsTK->Write();
}
1.7.3