moduleOccupancyTrend.cc

Go to the documentation of this file.

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

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

  gROOT->SetStyle("Plain");

  char* filelist;
  char* modulelist;

  filelist= argv[1];
  modulelist= argv[2];

  std::string detid;
  std::string hn;

  TLegend leg(0.1,0.7,0.2,0.9);
  leg.SetFillStyle(0);

  std::ifstream inmodules(modulelist);

  while(1){
    
    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(1){

      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(0);
      histo->SetStats(kFALSE);

      if (hn.find("Summary")==std::string::npos) histo->Scale(1/EvtNum);

      double numberPerEvent= histo->Integral();

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