MaterialBudgetHistos.cc

Go to the documentation of this file.

#include "Validation/Geometry/interface/MaterialBudgetHistos.h"
#include "Validation/Geometry/interface/MaterialBudgetData.h"

MaterialBudgetHistos::MaterialBudgetHistos(std::shared_ptr<MaterialBudgetData> data,
                                           std::shared_ptr<TestHistoMgr> mgr,
                                           const std::string& fileName)
    : MaterialBudgetFormat(data), hmgr(mgr) {
  theFileName = fileName;
  book();
}

void MaterialBudgetHistos::book() {
  edm::LogInfo("MaterialBudget") << " MaterialBudgetHistos: Booking Histos";
  hmgr->addHistoProf1(new TProfile("10", "MB prof Eta ", 250, -5., 5.));
  hmgr->addHisto1(new TH1F("11", "Eta ", 501, -5., 5.));
  hmgr->addHistoProf1(new TProfile("20", "MB prof Phi ", 180, -3.1416, 3.1416));
  hmgr->addHisto1(new TH1F("21", "Phi ", 360, -3.1416, 3.1416));
  hmgr->addHistoProf2(new TProfile2D("30", "MB prof Eta  Phi ", 250, -5., 5., 180, -3.1416, 3.1416));
  hmgr->addHisto2(new TH2F("31", "Eta vs Phi ", 501, -5., 5., 180, -3.1416, 3.1416));
}

void MaterialBudgetHistos::fillStartTrack() {}

void MaterialBudgetHistos::fillPerStep() {}

void MaterialBudgetHistos::fillEndTrack() {
  hmgr->getHisto1(11)->Fill(theData->getEta());
  hmgr->getHisto1(21)->Fill(theData->getPhi());
  hmgr->getHisto2(31)->Fill(theData->getEta(), theData->getPhi());

  hmgr->getHistoProf1(10)->Fill(theData->getEta(), theData->getTotalMB());
  hmgr->getHistoProf1(20)->Fill(theData->getPhi(), theData->getTotalMB());
  hmgr->getHistoProf2(30)->Fill(theData->getEta(), theData->getPhi(), theData->getTotalMB());
}

void MaterialBudgetHistos::endOfRun() {
  edm::LogInfo("MaterialBudget") << "MaterialBudgetHistos: Writing Histos ROOT file to" << theFileName;
  hmgr->save(theFileName);
}