MaterialBudgetTxt.cc

Go to the documentation of this file.

#include "Validation/Geometry/interface/MaterialBudgetTxt.h"
#include "Validation/Geometry/interface/MaterialBudgetData.h"
#include "G4EventManager.hh"
#include "G4Event.hh"

MaterialBudgetTxt::MaterialBudgetTxt(std::shared_ptr<MaterialBudgetData> data, const std::string& fileName)
    : MaterialBudgetFormat(data) {
  const char* fnamechar = fileName.c_str();
  theFile = new std::ofstream(fnamechar, std::ios::out);
  edm::LogInfo("MaterialBudget") << "MaterialBudgetTxt: Dumping  Material Budget to " << fileName;
  if (theFile->fail()) {
    edm::LogError("MaterialBudget") << "MaterialBudgetTxt: Error opening file" << fileName;
  }
}

MaterialBudgetTxt::~MaterialBudgetTxt() {}

void MaterialBudgetTxt::fillStartTrack() {
  (*theFile) << " Track " << G4EventManager::GetEventManager()->GetConstCurrentEvent()->GetEventID() << " "
             << theData->getEta() << " " << theData->getPhi() << std::endl;
  // + 1 was GEANT3 notation  (*theFile)<< " Track "<< G4EventManager::GetEventManager()->GetConstCurrentEvent()->GetEventID() + 1<< " " << theData->getEta() << " " << theData->getPhi() << std::endl;
}

void MaterialBudgetTxt::fillPerStep() {
  (*theFile) << "step " << theData->getTrkLen() << " " << theData->getPVname() << " " << theData->getPVcopyNo() << " "
             << theData->getTotalMB() << " " << theData->getRadLen() << std::endl;
  //-    std::cout << "step "<< theData->getTrkLen() << " " << theData->getPVname() << " " << theData->getPVcopyNo()  << " " << theData->getTotalMB() << " " << theData->getRadLen() << std::endl;
}

void MaterialBudgetTxt::fillEndTrack() {
  (*theFile) << G4EventManager::GetEventManager()->GetConstCurrentEvent()->GetEventID() << " "
             << "finalTrkMB " << theData->getTotalMB() << std::endl;
}

void MaterialBudgetTxt::endOfRun() { theFile->close(); }