PlotDB_Web.C

Go to the documentation of this file.

#include "SetupHistos.C"
#include "TMath.h"

/********************/
#include <algorithm>
#include <limits>
#include <cmath>
#include <vector>
#include <iostream>
/*********************/






void UpdateStamp(int ptype, int maxt){

  TDatime tnow(maxt);
  TString lasttime = TString::Format("Last updated on: %s", tnow.AsString());

  TLatex *ltx = MiscText(0.94,0.9,0.03,lasttime);

  ltx->SetTextAngle(-90);

}

void DrawConfig(int dbcfg){

  for(int db=1; db<=14; db++){
    TPave *box = new TPave(0.26+0.02*db,0.945,0.28+0.02*db,0.975,4,"BR NDC");
    box->SetLineWidth(2);
    box->SetShadowColor(0);
    box->SetFillStyle(1001);
    //box->DrawPave(0.17+0.03*db,0.95,0.20+0.03*db,0.98,4,"NDC BR");
//    box->SetFillColor(10);
    int bit = pow(2,db-1);
    int isgood = ( (dbcfg - dbcfg % bit) / bit ) % 2;
    if(isgood) box->SetFillColor(kGreen);
    else       box->SetFillColor(kWhite);
    box->Draw();
  }

}



void PlotDB_Web(int ptype=2){




  
//  TCanvas* ct=new TCanvas("ct","ct",500,500);
  
  gStyle->SetTimeOffset(0);

  TDatime t;
  int GMToffset = (t.Convert(true) - t.Convert());

  gSystem->Setenv("TZ","UTC");

  TString stdP1 = "run>11495&&par==1";
  t.Set();
  int maxt = t.Convert();
  int mint = TDatime(20130729,0).Convert();

//  TDatime(maxt).Print();

  int dt = 3600*24;
  if(ptype==1) dt *= 7;
  if(ptype==0) dt *= 30;

  TString tsel = TString::Format("&& lastsec > %d && lastsec < %d",maxt - dt + 2*GMToffset, maxt + 2*GMToffset);
  stdP1 += tsel;

  int tdiv = 706;
  int tbin = 3600*24;

  TString tlab = "Date (dd/mm/yy - Central Time)";
  TString tform = "%d/%m/%y";
  TString suf = "Month";

  if(ptype==1){
    tbin = 3600*6;
    tdiv = 408;
    tlab = "Date (dd/mm - Central Time)";
    tform = "%d/%m";
    suf = "Week";
  }
  if(ptype==2){
    tbin = 3600;
    tdiv = 309;
    tlab = "Central Time (hour:min)";
    tform = "%H:%M";
    suf = "Day";
  }

  TChain *mytree = new TChain("mytree");
  mytree->Add("MetricsTree_All.root");

//  MakeLongCanvas();

  TCanvas *c1 = new TCanvas("c1","c1",1000,600);
  c1->SetBottomMargin(0.12);
  c1->SetLeftMargin(0.12);
  c1->SetRightMargin(0.07);
  c1->SetTopMargin(0.07);



  

  int dbencoded, dbaencoded, run, subrun, par, pass_all, lastsec;
  int preliminary;

  mytree->SetBranchAddress("dbencoded",&dbencoded);
  mytree->SetBranchAddress("dbaencoded",&dbaencoded);
  mytree->SetBranchAddress("lastsec",&lastsec);
  mytree->SetBranchAddress("run",&run);
  mytree->SetBranchAddress("subrun",&subrun);
  mytree->SetBranchAddress("par",&par);
  mytree->SetBranchAddress("pass_all",&pass_all);
  mytree->SetBranchAddress("preliminary",&preliminary);

//  preliminary = false;

  int nentries = mytree->GetEntries();

  TH2D *h2 = new TH2D("","",dt/10,maxt-dt,maxt,14,1,15);

  vector<int> dbconfigs;

  for(int i=0; i<nentries; i++){

    //if(i%(nentries/1000)==0) cout << i*100.0/nentries << "% done" << endl;

    mytree->GetEntry(i);
//    if(par!=1 || !pass_all || run<16540) continue;
    if(par!=1 || lastsec - 2*GMToffset < maxt - dt){
      continue;
    }
//    if(run<16540) continue;
    //if( mytree2->GetEntries(TString::Format("run==%d && subrun==%d && pass_all",run,subrun)) ) continue;
    //if( !IsGood(run,subrun) ) continue;

    dbconfigs.push_back(dbencoded);
    
    for(int db=1; db<=14; db++){
      int bit = pow(2,db-1);
      int isgood = ( (dbencoded - dbencoded % bit) / bit ) % 2;
      int isactive = ( (dbaencoded - dbaencoded % bit) / bit ) % 2;
      if(!preliminary){
        if(isgood)        h2->Fill(lastsec - 2*GMToffset,db,2);
        else if(isactive) h2->Fill(lastsec - 2*GMToffset,db,1);
      }
      else{
        if(isgood)        h2->Fill(lastsec - 2*GMToffset,db,1.75);
        else if(isactive) h2->Fill(lastsec - 2*GMToffset,db,1.25);
      }
    }

  }

  if(dbconfigs.size()==0) dbconfigs.push_back(0);

  int *idx = new int[dbconfigs.size()];

   TMath::Sort(int(dbconfigs.size()),&dbconfigs[0],idx,kFALSE);
//  std::sort(dbconfigs.begin(),dbconfigs.end());
  
  int napp = 0;
  int maxapp = 0;
  int tempcfg = -1;
  int maxcfg = -1;
  
  for(int i=0; i<int(dbconfigs.size()); i++){
    if(tempcfg != dbconfigs[idx[i]]){
      if(napp > maxapp){
        maxapp = napp;
        maxcfg = tempcfg;
      }
      napp = 0;
      tempcfg = dbconfigs[idx[i]];
    }
    napp++;
  }
  if(napp > maxapp){
    maxapp = napp;
    maxcfg = tempcfg;
  }

  delete[] idx;
  
  SetHist(h2);

//  h3->GetXaxis()->SetRangeUser(11495,16300);
  h2->GetYaxis()->SetNdivisions(14);
  for(int i=1; i<=14; i++) h2->GetYaxis()->SetBinLabel(i,TString::Format("%d",i));
    
  h2->SetTitle(";"+tlab+";Diblock Number;");
  h2->GetXaxis()->SetTimeDisplay(1);
  h2->GetXaxis()->SetTimeFormat(tform);
  h2->GetXaxis()->SetNdivisions(tdiv);
  h2->SetStats(kFALSE);
  h2->GetXaxis()->CenterTitle();
  h2->GetYaxis()->CenterTitle();
  h2->GetYaxis()->SetLabelSize(0.06);

  // Get moodier colours for colz
  const Int_t NRGBs = 5;
  const Int_t NCont = 255;
  Double_t stops[NRGBs] = { 0.00, 0.25, 0.50, 0.75, 1.00 };
  Double_t red[NRGBs]   = { 1.00, 1.00, 1.00, 0.60, 0.00 };
  Double_t green[NRGBs] = { 0.00, 0.70, 1.00, 0.80, 0.80 };
  Double_t blue[NRGBs]  = { 0.00, 0.70, 1.00, 0.60, 0.00 };
  TColor::CreateGradientColorTable(NRGBs, stops, red, green, blue, NCont);
  gStyle->SetNumberContours(NCont);

  h2->GetZaxis()->SetRangeUser(1,2);

  h2->Draw("col");


// ******************************************************************
//  ct->cd();
// 
// *****************************************************************

  c1->SetGridy();

  if(dbconfigs.size()==1 && dbconfigs[0]==0) maxapp = 0;
  
  TString bestcfg = TString::Format("# of Subruns = %d; Percentage = %2.1f%%",
                                    maxapp,100.*maxapp/dbconfigs.size());

  MiscText(0.13,0.95,0.035,"Most Common = ");
  MiscText(0.58,0.95,0.035,bestcfg);
  DrawConfig(maxcfg);

  UpdateStamp(ptype,maxt);
  c1->SaveAs("pngfiles/FarDet-t02-P1GoodDataSelDetConfig"+suf+".png");
c1->Close();
  
}

TString GetConfig(int dbcfg){

  TString out = "";

  for(int db=1; db<=14; db++){
    int bit = pow(2,db-1);
    int isgood = ( (dbcfg - dbcfg % bit) / bit ) % 2;
    if(isgood) out += "1";
    else       out += "0";
  }
  
  return out;
  
}