Subversion Repositories f9daq

#include <stdio.h>
#include <time.h>
#include <stdlib.h>

#include <TROOT.h>
#include <TLatex.h>
#include <TH1F.h>
#include <TH2D.h>
#include <TH3D.h>
#include <TGraph.h>
#include <TCanvas.h>
#include <TLegend.h>
#include <TStyle.h>
#include <TSystem.h>
#include <TFile.h>
#include <TTree.h>
#include <TDirectory.h>
#include <TPaveText.h>

#include "base.h"
#include "info.h"
#include "savetoroot.h"

int info(int runNumber, int save) {

  int HAPDnumber; //Number of HAPD: 0, 1, 2 or 3!

  char text[256];
  char pdfname[256];
  char buf[256];
  char path[128];
  int vrstaPisave = 82;

  const char * serialNumber;
  std::string serialNumberTemp;
  std::string HAPDserialNumber,FEBserialNumber;

  unsigned int id, len, cas;
  int vmon[24], vset[24];
  int imon[24];
  //int status[6];
  //int min_imon[6], max_imon[6], min_vmon[6], max_vmon[6], min_status[6], max_status[6];
  int min_chipBias = 0, max_chipBias = 0;
  int sizeOfTree;
  time_t itime;
  TDatime dh(1995,01,01,00,00,00); // to je za SetTimeOffset

  char * viri[6] = {"Guard", "Chip A", "Chip B", "Chip C", "Chip D", "HV"};

  TFile * data[3];
  TCanvas * c;
  TPad * pad1, * pad2, * pad3, * pad4, * pad5;
  TPaveText  * sn, * timeStamp, * imeMeritve, * zacetekMeritve;
  TGraph * gr[10];
  TLegend * leg[2];
  TTree * tree;

  /********** GLOBALNE **********/
  gStyle->SetOptStat(0);
  //gStyle->SetOptTitle(1);
  gStyle->SetTextFont(vrstaPisave);
  gStyle->SetTitleFont(vrstaPisave);
  gStyle->SetTitleFontSize(0.08);

  gStyle->SetPalette(55);

  //gROOT->SetBatch(kTRUE); //Does not show canvases

  sprintf(path,"../modules/%04d/%04d_1_Treshold_Waveform.root",runNumber,runNumber);
  data[0] = new TFile(path);
  if (!file_exists(path)){
    printf("File '%s' does not exist.\n",path);
    return -1;
  }
  serialNumber = getSN(runNumber);
  printf("base.c return: %s\n", serialNumber);
  serialNumberTemp = serialNumber;
 
  sprintf(path,"../modules/%04d/%04d_2_2DX.root",runNumber,runNumber);
  data[1] = new TFile(path);
  if (!file_exists(path)){
    printf("File '%s' does not exist.\n",path);
    return -1;
  }
  sprintf(path,"../modules/%04d/%04d_3_2DY.root",runNumber,runNumber);
  data[2] = new TFile(path);
  if (!file_exists(path)){
    printf("File '%s' does not exist.\n",path);
    return -1;
  }

  for(HAPDnumber=0;HAPDnumber<4;HAPDnumber++){
         
          SaveToRootInit();

   //std::cout << serialNumberTemp << std::endl;
    FEBserialNumber = serialNumberTemp.substr(0,serialNumberTemp.find(","));
    serialNumberTemp.erase(0,FEBserialNumber.length()+1);
    if(!FEBserialNumber.compare("noserial")) continue;
    HAPDserialNumber = FEBserialNumber.substr(0,FEBserialNumber.find("/"));
    FEBserialNumber.erase(0,FEBserialNumber.find("/")+1);
 
    printf("%s \t %s \n",FEBserialNumber.c_str(),HAPDserialNumber.c_str());

    sprintf(buf,"%04d_%s_Monitor",runNumber,HAPDserialNumber.c_str());
    c = new TCanvas(buf,buf,0,0,1500,1200);
    /******* PRVI PAD *******/
    sprintf(buf,"%04d_%s_Monitor_pad1",runNumber,HAPDserialNumber.c_str());
    pad1 = new TPad(buf,"SN",0,0.95,1,1,0,0);
      pad1->SetFillStyle(4000);
      pad1->Draw();
      sn = new TPaveText(0.2,0,0.8,1,"ndc");
        sn->SetBorderSize(0);
        sn->SetFillColor(4000);
        sn->SetTextSize(.8);
        sn->SetTextAlign(22);
        sn->SetTextFont(vrstaPisave);

    /******* DRUGI PAD *******/
    sprintf(buf,"%04d_%s_Monitor_pad2",runNumber,HAPDserialNumber.c_str());
    pad2 = new TPad(buf,"Date",0.05,0.9,0.95,0.95,0,0);
      pad2->SetFillStyle(4000);
      pad2->Draw();
      timeStamp = new TPaveText(0.2,0,0.8,1,"ndc");
        timeStamp->SetBorderSize(0);
        timeStamp->SetFillColor(4000);
        timeStamp->SetTextSize(.4);
        timeStamp->SetTextAlign(22);
        timeStamp->SetTextFont(vrstaPisave);

    /******* TRETJI PAD *******/
    sprintf(buf,"%04d_%s_Monitor_pad3",runNumber,HAPDserialNumber.c_str());
    pad3 = new TPad(buf,"Napetosti",0,0,1,0.9,0,0);
      pad3->SetFillStyle(4000);
      pad3->Draw();
      pad3->cd();
      pad3->Divide(0,3);

    for (int k=0;k<3;k++){
      pad3->cd(k+1);

      sprintf(buf,"%04d_%s_Monitor_pad4",runNumber,HAPDserialNumber.c_str());
      pad4 = new TPad(buf,"Naslov",0.05,0.9,0.95,1,0,0);
        pad4->SetFillStyle(4000);
        pad4->Draw();
        imeMeritve = new TPaveText(0,0,1,1,"ndc");
          imeMeritve->SetBorderSize(0);
          imeMeritve->SetFillColor(kSpring+1);
          imeMeritve->SetTextSize(.6);
          imeMeritve->SetTextAlign(22);
          imeMeritve->SetTextFont(vrstaPisave);

        zacetekMeritve = new TPaveText(0.8,0,1,1,"ndc");
          zacetekMeritve->SetBorderSize(0);
          zacetekMeritve->SetFillColor(kSpring+1);
          zacetekMeritve->SetTextSize(.4);
          zacetekMeritve->SetTextAlign(22);
          zacetekMeritve->SetTextFont(vrstaPisave);

      switch (k){
        case 0: sprintf(text, "Treshold and Waveform");break;
        case 1: sprintf(text, "2D in X");break;
        case 2: sprintf(text, "2D in Y");break;
      }
     
      pad4->cd();
        imeMeritve->AddText(text);
        imeMeritve->Draw();

      pad3->cd(k+1);

      sprintf(buf,"%04d_%s_Monitor_pad5",runNumber,HAPDserialNumber.c_str());
      pad5 = new TPad(buf,"Grafi",0,0,1,.9,0,0);
        pad5->SetFillStyle(4000);
        pad5->Draw();
        pad5->cd();
        pad5->Divide(4,0);

      tree = (TTree*) data[k]->Get("Monitor");
        tree->SetBranchAddress("id",&id);
        tree->SetBranchAddress("len",&len);
        tree->SetBranchAddress("time",&cas);
        tree->SetBranchAddress("imon",&imon);
        tree->SetBranchAddress("vmon",&vmon);
        tree->SetBranchAddress("vset",&vset);
        //tree->SetBranchAddress("status",&status);

      sizeOfTree = tree->GetEntries();

      gr[0] = new TGraph(sizeOfTree);
      gr[1] = new TGraph(sizeOfTree);
      gr[2] = new TGraph(sizeOfTree);
      gr[3] = new TGraph(sizeOfTree);
      gr[4] = new TGraph(sizeOfTree);
      gr[5] = new TGraph(sizeOfTree);
      gr[6] = new TGraph(sizeOfTree);
      gr[7] = new TGraph(sizeOfTree);
      gr[8] = new TGraph(sizeOfTree);
      gr[9] = new TGraph(sizeOfTree);

      for (int i=0;i<sizeOfTree;i++){
        tree->GetEntry(i);

        cas = cas-2208988800; //(2208988800)CVI zapiše čas v sekundah od 00:00:00 1.1.1900 naprej, UNIX čas je od 00:00:00 1.1.1970
        // cas-u morš še nekaj odštet, kr kasneje ko delaš TimeOffset, le ta vodi čas od 00:00:00 1.1.1995 naprej.
        gr[0]->SetPoint(i,cas-788918400,vmon[HAPDnumber*6+0]);
        gr[1]->SetPoint(i,cas-788918400,vmon[HAPDnumber*6+1]-vset[HAPDnumber*6+1]);
        gr[2]->SetPoint(i,cas-788918400,vmon[HAPDnumber*6+2]-vset[HAPDnumber*6+2]);
        gr[3]->SetPoint(i,cas-788918400,vmon[HAPDnumber*6+3]-vset[HAPDnumber*6+3]);
        gr[4]->SetPoint(i,cas-788918400,vmon[HAPDnumber*6+4]-vset[HAPDnumber*6+4]);
        gr[5]->SetPoint(i,cas-788918400,vmon[HAPDnumber*6+5]);
        gr[6]->SetPoint(i,cas-788918400,imon[HAPDnumber*6+1]/1000.);
        gr[7]->SetPoint(i,cas-788918400,imon[HAPDnumber*6+2]/1000.);
        gr[8]->SetPoint(i,cas-788918400,imon[HAPDnumber*6+3]/1000.);
        gr[9]->SetPoint(i,cas-788918400,imon[HAPDnumber*6+4]/1000.);
        if (!i){
          itime = cas;
          zacetekMeritve->AddText(ctime(&itime));
        } else if (k==2 && i==sizeOfTree-1){
          itime = cas;
          sprintf(buf,"End of measuring: %s",ctime(&itime));
          timeStamp->AddText(buf);
        }
       
        for (int j=1;j<5;j++){
          if (!i && j==1) {
            min_chipBias = max_chipBias = vmon[HAPDnumber*6+j]-vset[HAPDnumber*6+j];
          } else {
            if (min_chipBias > vmon[HAPDnumber*6+j]-vset[HAPDnumber*6+j]) min_chipBias = vmon[HAPDnumber*6+j]-vset[HAPDnumber*6+j];
            else if (max_chipBias < vmon[HAPDnumber*6+j]-vset[HAPDnumber*6+j]) max_chipBias = vmon[HAPDnumber*6+j]-vset[HAPDnumber*6+j];
          }
        }
      }

      leg[0] = new TLegend(0.55,0.1,.9,0.35);
      leg[1] = new TLegend(0.70,0.70,.9,0.9);
      leg[0]->SetFillColor(0);
      leg[1]->SetFillColor(0);
      for (int i=0;i<10;i++){
        gr[i]->GetYaxis()->CenterTitle();
        gr[i]->GetYaxis()->SetTitleFont(vrstaPisave);
        gr[i]->GetYaxis()->SetLabelFont(vrstaPisave);
        gr[i]->GetYaxis()->SetTitleOffset(1.6);


        gr[i]->GetXaxis()->SetTitleFont(vrstaPisave);
        gr[i]->GetXaxis()->SetLabelFont(vrstaPisave);
        //gr[i]->GetXaxis()->CenterTitle();
        gr[i]->GetXaxis()->SetTimeDisplay(1);
        gr[i]->GetXaxis()->SetTimeOffset(dh.Convert());
        gr[i]->GetXaxis()->SetTimeFormat("#splitline{%H:%M}{%d-%m}");
        gr[i]->GetXaxis()->SetLabelOffset(0.02);
        gr[i]->GetXaxis()->SetTitleOffset(1.5);


        if(!i){
          pad5->cd(1);
         
          sprintf(text, "Guard bias (RMS = %.0f mV);Time;Bias [mV]", gr[i]->GetRMS(2));
          gr[i]->SetTitle(text);
          gr[i]->GetYaxis()->SetRangeUser(170000,180000);

          gr[i]->Draw("AL");
         
          sprintf(text, "Guard_%s_T_%d", HAPDserialNumber.c_str(), k);
          SaveToRootAddGraph(gr[i], text);
        } else if (i<5 && i!=0){
          pad5->cd(2);
          sprintf(text, "%s (RMS = %.0f mV)", viri[i], gr[i]->GetRMS(2));
          leg[0]->AddEntry(gr[i],text,"L");


          //gr[i]->GetYaxis()->SetRangeUser(min_chipBias-50,max_chipBias+50);
          gr[i]->GetYaxis()->SetRangeUser(-1000,1000);

          gr[i]->SetTitle("Chip (vmon - vset) difference;Time;Bias difference [mV]");
          if (i==1) gr[i]->Draw("AL");

          gr[i]->SetLineColor(i);
          gr[i]->SetLineWidth(1);

          gr[i]->Draw("LSAME");
         
          sprintf(text, "BiasDifference_%s_CHIP_%d_T_%d", HAPDserialNumber.c_str(), i-1, k);
                  SaveToRootAddGraph(gr[i], text);
        } else if (i>5){
          pad5->cd(3);
          sprintf(text, "%s", viri[i-5]);
          leg[1]->AddEntry(gr[i],text,"L");

          gr[i]->GetYaxis()->SetRangeUser(-1,11);

          gr[i]->SetTitle("Chip leakage current;Time;Leakage current [mA]");
          if (i==6) gr[i]->Draw("AL");

          gr[i]->SetLineColor(i);
          gr[i]->SetLineWidth(1);

          gr[i]->Draw("LSAME");
         
          sprintf(text, "LeakageCurrent_%s_CHIP_%d_T_%d", HAPDserialNumber.c_str(), i-6, k);
                  SaveToRootAddGraph(gr[i], text);
        } else if (i==5){
          pad5->cd(4);
          sprintf(text, "HV (RMS = %.0f mV);Time;Bias [mV]", gr[i]->GetRMS(2));
          gr[i]->GetYaxis()->SetRangeUser(6900000,7100000);

          gr[i]->SetTitle(text);

          gr[i]->Draw("AL");
         
          sprintf(text, "HV_%s_T_%d", HAPDserialNumber.c_str(), k);
                  SaveToRootAddGraph(gr[i], text);
        }
      }

      pad5->cd(2);
        leg[0]->Draw();
      pad5->cd(3);
        leg[1]->Draw();
     
      pad4->cd();
        zacetekMeritve->Draw();
    }

    pad1->cd();
      sprintf(buf,"RUN#: %04d (pos %d)   FEB#: %s   HAPD#: %s",runNumber, HAPDnumber, FEBserialNumber.c_str(),HAPDserialNumber.c_str());
      sn->AddText(buf);
      sn->Draw();
   
    pad2->cd();
      timeStamp->Draw();

    c->Modified();
    c->Update();

    if (save) {
      switch (save) {
        case 1: sprintf(pdfname,"../modules/%04d/%04d_%s_0_Info.pdf",runNumber,runNumber,HAPDserialNumber.c_str());break;
        case 2: sprintf(pdfname,"../modules/%04d/%04d_%s.pdf(",runNumber,runNumber,HAPDserialNumber.c_str());break;
       
      }
      c->SaveAs(pdfname,"pdf");
    }
    sprintf(pdfname,"../modules/%04d/%04d_%s_out.root",runNumber,runNumber,HAPDserialNumber.c_str());
    SaveToRootWrite(pdfname, HAPDserialNumber.c_str());
  }
  for(int k=0;k<3;k++) delete data[k];

  return 0;
}