G4RegionReporter.cc

Go to the documentation of this file.

#include "SimG4Core/Application/interface/G4RegionReporter.h"

#include "G4Region.hh"
#include "G4RegionStore.hh"
#include "G4LogicalVolume.hh"
#include "G4ProductionCuts.hh"
#include "G4SystemOfUnits.hh"

#include <iostream>
#include <iomanip>
#include <fstream>

G4RegionReporter::G4RegionReporter()
{}
  
G4RegionReporter::~G4RegionReporter()
{}

void G4RegionReporter::ReportRegions(const std::string& ss)
{
  std::ofstream fout(ss.c_str(), std::ios::out);
  if(fout.fail()){
    std::cout << "SimG4CoreApplication ReportRegions WARNING : "
          << "error opening file <" << ss << ">" << std::endl;
    return;
  } 
  G4RegionStore* regStore =  G4RegionStore::GetInstance();  

  unsigned int numRegions= regStore->size();

  unsigned int i;

  fout << std::endl;
  fout << "#---------------------------------------------------------------------";
  fout << "------------------------------------" << std::endl;
  fout << "## List of Regions, root logical volumes and cuts. " << std::endl;
  fout << "## Number of regions = " << numRegions << std::endl;

  //  Banner
  fout << "# " << std::setw(24) << " Region, " 
         << std::setw(38) << " LogicalVolume, " 
         << " Cuts:Gamma, Electron, Positron, Proton, Units" << std::endl;
  fout << "#---------------------------------------------------------------------";  
  fout << "------------------------------------" << std::endl;

  for( i=0; i<numRegions; ++i)
  {
     G4Region* region = regStore->at(i);
     G4ProductionCuts* prodCuts = region->GetProductionCuts(); 

     G4LogicalVolume* lv;

     G4double lengthUnit    = CLHEP::mm; 
     G4String lengthUnitName= "mm";
     unsigned int pmax= 4;  // g, e-, e+, proton

     std::vector<G4LogicalVolume*>::iterator rootLVItr
                         = region->GetRootLogicalVolumeIterator();
     size_t numRootLV = region->GetNumberOfRootVolumes();

     for(size_t iLV=0; iLV < numRootLV; ++iLV, ++rootLVItr )
     {
        // Cover each root logical volume in this region

        //Set the couple to the proper logical volumes in that region
        lv = *rootLVItr;

        // fout << " Region=" << region->GetName()
        //     << " Logical-Volume = " << lv->GetName();
        char quote='"';
        std::ostringstream regName, lvName;
        regName << quote << region->GetName() << quote;
        lvName << quote << lv->GetName()     << quote; 
        fout << " " << std::setw(26) << regName.str() << " ,"; 
        fout << " " << std::setw(36) <<  lvName.str() << " ,"; 

        unsigned int ic;
        for( ic= 0; ic< pmax; ++ic)
        {
           G4double cutLength = prodCuts->GetProductionCut(ic); 
           fout << " " << std::setw(5) << cutLength / lengthUnit;
           if (ic < pmax-1)       {
         fout << " , "; 
           } else {
         fout << " ,   " << lengthUnitName;
           }
        }
        fout << std::endl;
     }
  }
  fout << "#---------------------------------------------------------------------";  
  fout << "------------------------------------" << std::endl;
  fout << std::endl;
  fout.close();
}