CbmLitRadLengthQaReport.cxx

Go to the documentation of this file.

 
#include "CbmLitRadLengthQaReport.h"
#include "CbmDrawHist.h"
#include "CbmHistManager.h"
#include "CbmReportElement.h"
 
#include "TProfile2D.h"
 
CbmLitRadLengthQaReport::CbmLitRadLengthQaReport() : CbmSimulationReport() {
  SetReportName("radlength_qa");
}
 
CbmLitRadLengthQaReport::~CbmLitRadLengthQaReport() {}
 
void CbmLitRadLengthQaReport::Create() {
  Out().precision(3);
  Out() << R()->DocumentBegin();
  Out() << R()->Title(0, GetTitle());
 
  PrintCanvases();
 
  Out() << R()->DocumentEnd();
}
 
void CbmLitRadLengthQaReport::Draw() {
  DrawDetector("Total");
  DrawDetector("Mvd");
  DrawDetector("Sts");
  DrawDetector("Rich");
  DrawDetector("Trd");
  DrawDetector("Much");
  DrawDetector("Tof");
 
  DrawDetectorStation("Mvd");
  DrawDetectorStation("Sts");
  DrawDetectorStation("Trd");
  DrawDetectorStation("Much");
  DrawDetectorStation("MuchAbsorber");
}
 
void CbmLitRadLengthQaReport::DrawDetector(const string& detName) {
  if (HM()->Exists("hrl_Thickness_" + detName + "_H1")) {
    TCanvas* canvas1 =
      CreateCanvas(string("hrl_Thickness_" + detName + "_H1").c_str(),
                   string("hrl_Thickness_" + detName + "_H1").c_str(),
                   1200,
                   1000);
    DrawH1(HM()->H1("hrl_Thickness_" + detName + "_H1"));
  }
 
  if (HM()->Exists("hrl_Thickness_" + detName + "_P2")) {
    TCanvas* canvas2 =
      CreateCanvas(string("hrl_Thickness_" + detName + "_P2").c_str(),
                   string("hrl_Thickness_" + detName + "_P2").c_str(),
                   1200,
                   1000);
    DrawH2(HM()->P2("hrl_Thickness_" + detName + "_P2"));
  }
 
  if (HM()->Exists("hrl_RadThickness_" + detName + "_H1")) {
    TCanvas* canvas3 =
      CreateCanvas(string("hrl_RadThickness_" + detName + "_H1").c_str(),
                   string("hrl_RadThickness_" + detName + "_H1").c_str(),
                   1200,
                   1000);
    DrawH1(HM()->H1("hrl_RadThickness_" + detName + "_H1"));
  }
 
  if (HM()->Exists("hrl_RadThickness_" + detName + "_P2")) {
    TCanvas* canvas4 =
      CreateCanvas(string("hrl_RadThickness_" + detName + "_P2").c_str(),
                   string("hrl_RadThickness_" + detName + "_P2").c_str(),
                   1200,
                   1000);
    DrawH2(HM()->P2("hrl_RadThickness_" + detName + "_P2"));
  }
 
  if (HM()->Exists("hrl_ThicknessSilicon_" + detName + "_H1")) {
    TCanvas* canvas5 =
      CreateCanvas(string("hrl_ThicknessSilicon_" + detName + "_H1").c_str(),
                   string("hrl_ThicknessSilicon_" + detName + "_H1").c_str(),
                   1200,
                   1000);
    DrawH1(HM()->H1("hrl_ThicknessSilicon_" + detName + "_H1"));
  }
 
  if (HM()->Exists("hrl_ThicknessSilicon_" + detName + "_P2")) {
    TCanvas* canvas6 =
      CreateCanvas(string("hrl_ThicknessSilicon_" + detName + "_P2").c_str(),
                   string("hrl_ThicknessSilicon_" + detName + "_P2").c_str(),
                   1200,
                   1000);
    DrawH2(HM()->P2("hrl_ThicknessSilicon_" + detName + "_P2"));
  }
}
 
void CbmLitRadLengthQaReport::DrawDetectorStation(const string& detName) {
  DrawH1ByPattern("hrl_RadThickness_" + detName + "_.+_H1");
  DrawH2ByPattern("hrl_RadThickness_" + detName + "_.+_P2");
  DrawH1ByPattern("hrl_Thickness_" + detName + "_.+_H1");
  DrawH2ByPattern("hrl_Thickness_" + detName + "_.+_P2");
  DrawH1ByPattern("hrl_ThicknessSilicon_" + detName + "_.+_H1");
  DrawH2ByPattern("hrl_ThicknessSilicon_" + detName + "_.+_P2");
}
 
ClassImp(CbmLitRadLengthQaReport)