SOF_profiles.cc

Go to the documentation of this file.

#include "SOF_profiles.h"
#include "TCanvas.h"
#include "TFile.h"
#include "TProfile.h"
#include "TH1F.h"
#include "TDirectory.h"
#include "TLine.h"
#include "TGaxis.h"
#include "TLegend.h"
#include <iostream>

TCanvas* printSOF(TFile* file, const int run, const int firstLS, const int zoom) {
  TCanvas* canout = nullptr;

  char rname[400];
  sprintf(rname, "run_%d", run);

  if (file == nullptr)
    return canout;

  TProfile* badmod = nullptr;
  TH1F* evtanydcs = nullptr;
  TH1F* evtdcson = nullptr;
  TH1F* evtnostrip = nullptr;

  char dname[422];

  sprintf(dname, "ssqhistory/%s", rname);
  if (file->cd(dname)) {
    badmod = (TProfile*)gDirectory->Get("badmodrun_HVDCS");
  }

  sprintf(dname, "eventtimedistranydcs/%s", rname);
  if (file->cd(dname)) {
    evtanydcs = (TH1F*)gDirectory->Get("orbit");
  }

  sprintf(dname, "eventtimedistribution/%s", rname);
  if (file->cd(dname)) {
    evtdcson = (TH1F*)gDirectory->Get("orbit");
  }

  sprintf(dname, "eventtimedistrnostrip/%s", rname);
  if (file->cd(dname)) {
    evtnostrip = (TH1F*)gDirectory->Get("orbit");
  }

  if (badmod && evtanydcs && evtdcson && evtnostrip) {
    badmod->SetStats(kFALSE);
    badmod->SetLineColor(kGreen);
    badmod->SetMarkerColor(kGreen);
    evtanydcs->SetStats(kFALSE);
    evtanydcs->SetLineColor(kBlue);
    evtdcson->SetStats(kFALSE);
    evtdcson->SetLineColor(kRed);
    evtnostrip->SetStats(kFALSE);
    evtnostrip->SetLineColor(kRed);
    evtnostrip->SetFillColor(kRed);
    evtnostrip->SetFillStyle(1000);

    canout = new TCanvas;

    badmod->Draw();
    badmod->GetYaxis()->SetRangeUser(1, 100000);
    badmod->GetYaxis()->SetTitle("");
    badmod->SetTitle(rname);
    badmod->GetXaxis()->SetRangeUser((firstLS - zoom) * 262144, (firstLS + zoom) * 262144);
    evtanydcs->Draw("same");
    evtdcson->Draw("same");
    evtnostrip->Rebin(int(badmod->GetBinWidth(1) / evtnostrip->GetBinWidth(1)));
    std::cout << "rebin " << int(badmod->GetBinWidth(1) / evtnostrip->GetBinWidth(1)) << std::endl;
    evtnostrip->Draw("same");

    TGaxis* lsaxis = new TGaxis((firstLS - zoom) * 262144,
                                100000,
                                (firstLS + zoom) * 262144,
                                100000,
                                firstLS - zoom + 1,
                                firstLS + zoom + 1,
                                2 * zoom,
                                "-SM");
    //    TGaxis* lsaxis = new TGaxis(badmod->GetXaxis()->GetXmin(),100000,
    //              badmod->GetXaxis()->GetXmax(),100000,
    //              badmod->GetXaxis()->GetXmin()/262144+1,badmod->GetXaxis()->GetXmax()/262144+1,50,"-SM");
    lsaxis->Draw();

    TLine* line = new TLine((firstLS - 1) * 262144, 1, (firstLS - 1) * 262144, 100000);
    line->SetLineWidth(2);
    line->SetLineStyle(2);
    line->Draw();

    TLegend* leg = new TLegend(.5, .65, .9, .85, "");
    leg->AddEntry(badmod, "Modules with HV off", "l");
    leg->AddEntry(evtanydcs, "Events with any DCS bit", "l");
    leg->AddEntry(evtdcson, "Events with DCS bit ON", "l");
    leg->AddEntry(evtnostrip, "DCS bit ON No strip clus (masked FED)", "f");
    leg->AddEntry(line, "first good LS (DCSonly JSON)", "l");

    leg->SetFillStyle(0);
    leg->Draw();

    canout->SetLogy(1);
  }

  return canout;
}