/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_13_patch3/src/Validation/Geometry/src/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( MaterialBudgetData* data, const std::string& fileName ): MaterialBudgetFormat( data )
{
  const char * fnamechar = fileName.c_str();
  theFile = new std::ofstream(fnamechar, std::ios::out);
  std::cout <<"Dumping  Material Budget to " << fileName << std::endl;
  if (theFile->fail()){
    std::cerr <<" Error opening file" << fileName << std::endl;
  }
}


MaterialBudgetTxt::~MaterialBudgetTxt()
{
  theFile->close();
}


void MaterialBudgetTxt::fillStartTrack()
{
   std::cout << " Track "<< G4EventManager::GetEventManager()->GetConstCurrentEvent()->GetEventID() << " " << theData->getEta() << " " << theData->getPhi() << std::endl;
  (*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;
}