/**
 * \file CbmLitRadLengthQaReport.cxx
 * \author Andrey Lebedev <andrey.lebedev@gsi.de>
 * \date 2013
 */
#include "CbmLitRadLengthQaReport.h"
#include "CbmReportElement.h"
#include "CbmHistManager.h"
#include "CbmDrawHist.h"

#include "TH1.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()
{
   TCanvas* canvas1 = CreateCanvas("hrl_Thickness_Total_H1", "hrl_Thickness_Total_H1", 1000, 1000);
   DrawH1(HM()->H1("hrl_Thickness_Total_H1"));

   TCanvas* canvas2 = CreateCanvas("hrl_Thickness_Total_P2", "hrl_Thickness_Total_P2", 1000, 1000);
   DrawH2(HM()->P2("hrl_Thickness_Total_P2"));

   TCanvas* canvas3 = CreateCanvas("hrl_RadThickness_Total_H1", "hrl_RadThickness_Total_H1", 1000, 1000);
   DrawH1(HM()->H1("hrl_RadThickness_Total_H1"));

   TCanvas* canvas4 = CreateCanvas("hrl_RadThickness_Total_P2", "hrl_RadThickness_Total_P2", 1000, 1000);
   DrawH2(HM()->P2("hrl_RadThickness_Total_P2"));

   TCanvas* canvas5 = CreateCanvas("hrl_ThicknessSilicon_Total_H1", "hrl_ThicknessSilicon_Total_H1", 1000, 1000);
   DrawH1(HM()->H1("hrl_ThicknessSilicon_Total_H1"));

   TCanvas* canvas6 = CreateCanvas("hrl_ThicknessSilicon_Total_P2", "hrl_ThicknessSilicon_Total_P2", 1000, 1000);
   DrawH2(HM()->P2("hrl_ThicknessSilicon_Total_P2"));

   DrawH1ByPattern("hrl_RadThickness_Trd_.+_.+_H1");
   DrawP2ByPattern("hrl_RadThickness_Trd_.+_.+_P2");
   DrawH1ByPattern("hrl_Thickness_Trd_.+_.+_H1");
   DrawP2ByPattern("hrl_Thickness_Trd_.+_.+_P2");
   DrawH1ByPattern("hrl_ThicknessSilicon_Trd_.+_.+_H1");
   DrawP2ByPattern("hrl_ThicknessSilicon_Trd_.+_.+_P2");
}

void CbmLitRadLengthQaReport::DrawH1ByPattern(
      const string& pattern)
{
   std::vector<TH1*> histos = HM()->H1Vector(pattern);
   Int_t nofHistos = histos.size();
   for (Int_t i = 0; i < nofHistos; i++) {
      TH1* hist = histos[i];
      TCanvas* canvas = CreateCanvas(hist->GetName(), hist->GetName(), 1000, 1000);
      DrawH1(hist);
   }
}

void CbmLitRadLengthQaReport::DrawP2ByPattern(
      const string& pattern)
{
   std::vector<TProfile2D*> profiles = HM()->P2Vector(pattern);
   Int_t nofProfiles = profiles.size();
   for (Int_t i = 0; i < nofProfiles; i++) {
      TProfile2D* profile = profiles[i];
      TCanvas* canvas = CreateCanvas(profile->GetName(), profile->GetName(), 1000, 1000);
      DrawH2(profile);
   }
}

ClassImp(CbmLitRadLengthQaReport)