#include <iostream>
#include <fstream>
#include <sstream>
#include <string>
#include <TROOT.h>
#include <TStyle.h>
#include "TFile.h"
#include "THStack.h"
#include "TH1D.h"
#include "TDirectoryFile.h"
#include "TCanvas.h"
#include "TLine.h"
#include "TLegend.h"
#include "TRint.h"

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

Function Documentation

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

RBX:

errA with average Amplitudes

Summed Amplitude Plots:

RBX:

errA with average Amplitudes

Summed Amplitude Plots:

RBX:

Prepare maps of good/bad channels:

Definition at line 17 of file moduleOccupancyTrend.cc.

                                  {
  
   gROOT->SetStyle("Plain");
  
   char* filelist;
   char* modulelist;
  
   filelist= argv[1];
   modulelist= argv[2];
  
   int stripmin = 999;
   int stripmax = 999;
  
   if (argv[3] && argv[4]){
     stripmin=atoi(argv[3]);
     stripmax=atoi(argv[4]);
   }
  
   std::string detid;
   std::string hn;
  
   TLegend leg(0.1,0.7,0.2,0.9);
   leg.SetFillStyle(0);
  
   std::ifstream inmodules(modulelist);
  
   while(true){
     
     inmodules >> detid;
     if (!inmodules.good()) break;
  
     hn="ClusterDigiPosition__det__" + detid;//std::to_string(detid);
     
     TCanvas c1("c1","c1",1600,900);
   
     c1.SetBatch(kTRUE);
     c1.SetLogy(1);
     c1.SetGridy(1);
  
     //cout << detid << endl;
     std::ifstream fileToCountLines(filelist);
     std::size_t lines_count =0;
     std::string line;
  
     while (std::getline(fileToCountLines , line))
         ++lines_count;
  
     const float dim=lines_count;
  
     std::ifstream filesin(filelist);
  
     std::string filename;
     int k=0;
  
     TH1D* trend=new TH1D("trend","trend",int(dim),0.5,dim+0.5);
     trend->SetMarkerSize(3);
     trend->SetMarkerStyle(8);
     trend->SetMarkerColor(4);
  
     std::string ttitle=hn+"_trend";
     trend->SetTitle(ttitle.c_str());
  
     double max=-1;
     double min=1000000;
  
     leg.Clear();
     TFile* fin;
  
     THStack *hs = new THStack("hs","");
  
  
     while(true){
  
       filesin >> filename;
       if (!filesin.good()) break;      
  
       std::string runNum= filename.substr(filename.find("run_")+4, 6);
       //std::cout << runNum << std::endl;
      
       fin=TFile::Open(filename.c_str());
  
       if (!fin) {std::cout << "Cannot open file " << filename.c_str() << std::endl;
     return 0;
       }
       
       TH1D* Events=(TH1D*)fin->Get("TotEvents");
  
       double EvtNum=Events->GetBinContent(1);
  
       TH1D* histo=(TH1D*) fin->Get(hn.c_str());
       
       if (!histo) {std::cout << "Cannot open histo " << hn.c_str() << std::endl;
     return 0;
       }
  
       histo->SetDirectory(nullptr);
       histo->SetStats(kFALSE);
  
       if (hn.find("Summary")==std::string::npos) histo->Scale(1/EvtNum);
  
       double numberPerEvent;
       if (stripmin==999 && stripmax==999) numberPerEvent=histo->Integral();
       else numberPerEvent=histo->Integral(stripmin,stripmax);
  
       
       if (max<=histo->GetBinContent(histo->GetMaximumBin())) max=histo->GetBinContent(histo->GetMaximumBin());
       if (min>histo->GetBinContent(histo->GetMinimumBin())) min=histo->GetBinContent(histo->GetMinimumBin());
  
       histo->SetLineColor(k+1);
       histo->SetMarkerStyle(9);
       histo->SetMarkerColor(k+1);
            
       trend->SetBinContent(k+1,numberPerEvent);
       trend->GetXaxis()->SetBinLabel(k+1,runNum.c_str());
  
       leg.AddEntry(histo, runNum.c_str(), "L");
       hs->Add(histo);
       k++;
  
       fin->Close();
     }
  
     if (min==0) min=1.e-6;
  
     max=max*10;// in a way one can read the legend
  
     hs->SetMaximum(max);
     hs->SetMinimum(min);
     hs->SetTitle(hn.c_str());
  
     hs->Draw("nostack");
  
     TLine l;
     
     l.SetLineColor(4);
     l.SetLineStyle(2);
     l.SetLineWidth(3);
  
     l.DrawLine(128,min,128,max);
     l.DrawLine(256,min,256,max);
     l.DrawLine(384,min,384,max);
     l.DrawLine(384,min,384,max);
     l.DrawLine(512,min,512,max);
     l.DrawLine(640,min,640,max);
  
     leg.Draw();
     std::string outname= hn+"_Super.png";
     c1.SaveAs(outname.c_str());
     
     c1.SetGridx(1);
     
     double mintrend=0;
     if (trend->GetMinimum()==0) mintrend=1.e-4;
     else mintrend=trend->GetMinimum();
  
     trend->SetStats(kFALSE);
  
     trend->GetYaxis()->SetRangeUser(mintrend*0.5,trend->GetMaximum()*2);
     trend->Draw("P");
     outname=hn+"_Trend.png";
     c1.SaveAs(outname.c_str());
     c1.Clear();
  
     delete trend;      
     
   }  
   return 0;
 }

References cmsBatch::argv, alignmentValidation::c1, gather_cfg::cout, MillePedeFileConverter_cfg::e, harvestHGCalValidationPlots::filelist, corrVsCorr::filename, groupFilesInBlocks::fin, timingPdfMaker::histo, dqmdumpme::k, cmsLHEtoEOSManager::l, mps_splice::line, SiStripPI::max, min(), DeadROCCounter_Phase1::outname, DeadROC_duringRun::runNum, and AlCaHLTBitMon_QueryRunRegistry::string.

Functions

Function Documentation

◆ main()