CMSG4RegionReporter Class Reference

#include <CMSG4RegionReporter.h>

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

Detailed Description

Definition at line 6 of file CMSG4RegionReporter.h.

Constructor & Destructor Documentation

CMSG4RegionReporter::CMSG4RegionReporter

(

)

Definition at line 14 of file CMSG4RegionReporter.cc.

{}

CMSG4RegionReporter::~CMSG4RegionReporter

(

)

Definition at line 16 of file CMSG4RegionReporter.cc.

{}

Member Function Documentation

void CMSG4RegionReporter::ReportRegions

(

const std::string &

ss

)

Definition at line 18 of file CMSG4RegionReporter.cc.

References groupFilesInBlocks::fout, mps_fire::i, MillePedeFileConverter_cfg::out, and HLT_2018_cff::region.

Referenced by CMSG4CheckOverlap::CMSG4CheckOverlap().

                                                           {
  std::ofstream fout(ss.c_str(), std::ios::out);
  if (fout.fail()) {
    edm::LogWarning("SimG4CoreGeometry") << "CMSG4RegionReporter: file <" << ss
                                         << "> is not opened - no report provided";
    return;
  }
  G4RegionStore* regStore = G4RegionStore::GetInstance();

  unsigned int numRegions = regStore->size();

  unsigned int i;

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

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

  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 << "\n";
    }
  }
  fout << "#---------------------------------------------------------------------";
  fout << "------------------------------------"
       << "\n";
  fout << "\n";
  fout.close();
}