G4RegionReporter Class Reference

#include <G4RegionReporter.h>

Public Member Functions
	G4RegionReporter ()
void	ReportRegions (const std::string &ss)
	~G4RegionReporter ()

Detailed Description

Definition at line 6 of file G4RegionReporter.h.

Constructor & Destructor Documentation

G4RegionReporter::G4RegionReporter

(

)

Definition at line 13 of file G4RegionReporter.cc.

{}

G4RegionReporter::~G4RegionReporter

(

)

Definition at line 16 of file G4RegionReporter.cc.

{}

Member Function Documentation

void G4RegionReporter::ReportRegions

(

const std::string &

ss

)

Definition at line 19 of file G4RegionReporter.cc.

References gather_cfg::cout, groupFilesInBlocks::fout, i, and MillePedeFileConverter_cfg::out.

Referenced by RunManagerMT::initG4(), and RunManager::initG4().

{
  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();
}