TG4ModelConfigurationManager.cxx

Go to the documentation of this file.

//------------------------------------------------
// The Geant4 Virtual Monte Carlo package
// Copyright (C) 2007 - 2015 Ivana Hrivnacova
// All rights reserved.
//
// For the licensing terms see geant4_vmc/LICENSE.
// Contact: root-vmc@cern.ch
//-------------------------------------------------
 
 
#include "TG4ModelConfigurationManager.h"
#include "TG4G3ControlVector.h"
#include "TG4GeometryServices.h"
#include "TG4Limits.h"
#include "TG4Medium.h"
#include "TG4MediumMap.h"
#include "TG4ModelConfiguration.h"
#include "TG4ModelConfigurationMessenger.h"
#include "TG4PhysicsManager.h"
 
#include <G4LogicalVolume.hh>
#include <G4LogicalVolumeStore.hh>
#include <G4Material.hh>
#include <G4ProductionCuts.hh>
#include <G4Region.hh>
#include <G4RegionStore.hh>
#include <G4TransportationManager.hh>
 
#ifdef G4MULTITHREADED
namespace
{
// Mutex to lock creating regions
G4Mutex setRegionsNamesMutex = G4MUTEX_INITIALIZER;
G4Mutex createRegionsMutex = G4MUTEX_INITIALIZER;
} // namespace
#endif
 
//_____________________________________________________________________________
TG4ModelConfigurationManager::TG4ModelConfigurationManager(
  const G4String& name, const G4String& availableModels)
  : TG4Verbose(G4String(name).append("ConfigurationManager")),
    // make the verbose name unique
    fMessenger(0),
    fName(name),
    fAvailableModels(availableModels),
    fVector(),
    fCreateRegionsDone(false)
{
 
  if (VerboseLevel() > 1) {
    G4cout << "TG4ModelConfigurationManager::TG4ModelConfigurationManager"
           << G4endl;
  }
 
  fMessenger = new TG4ModelConfigurationMessenger(this, availableModels);
}
 
//_____________________________________________________________________________
TG4ModelConfigurationManager::~TG4ModelConfigurationManager()
{
 
  delete fMessenger;
 
  ModelConfigurationVector::iterator it;
  for (it = fVector.begin(); it != fVector.end(); it++) {
    delete *it;
  }
}
 
//
// private methods
//
 
//_____________________________________________________________________________
void TG4ModelConfigurationManager::SetRegionsNames()
{
 
#ifdef G4MULTITHREADED
  G4AutoLock lm(&setRegionsNamesMutex);
#endif
  TG4MediumMap* mediumMap = TG4GeometryServices::Instance()->GetMediumMap();
 
  ModelConfigurationVector::iterator it;
  for (it = fVector.begin(); it != fVector.end(); it++) {
 
    const std::vector<G4String>& mediaNames = (*it)->GetRegionsMedia();
 
    // Process the vector and generate a new list with material names
    std::vector<G4String>::const_iterator itr;
    for (itr = mediaNames.begin(); itr != mediaNames.end(); itr++) {
 
      std::vector<TG4Medium*> media;
      G4bool isPattern = (*itr).find("*") != std::string::npos;
      if (isPattern) {
        // G4cout << "processing " << (*itr) << " as pattern." << G4endl;
        mediumMap->GetMedia(*itr, media);
      }
      else {
        // G4cout << "processing " << (*itr) << " as single medium." << G4endl;
        TG4Medium* medium = mediumMap->GetMedium(*itr);
        if (!medium) {
          // Add warning
          TString text = "Medium ";
          text += (*itr).data();
          text += " not found.";
          text += TG4Globals::Endl();
          text += "The special physics model will not be applied.";
          TG4Globals::Warning(
            "TG4ModelConfigurationManager", "SetRegionsNames", text);
          continue;
        }
        media.push_back(medium);
      }
 
      std::vector<TG4Medium*>::const_iterator itrm;
      for (itrm = media.begin(); itrm != media.end(); itrm++) {
        G4String materialName = (*itrm)->GetMaterial()->GetName();
        (*it)->SetOneRegion(materialName);
      }
    }
 
    if (VerboseLevel() > 1) {
      // Print the new regions names corresponding to materials
      G4cout << "TG4ModelConfigurationManager::SetRegionsNames: ";
      const std::vector<G4String>& regions = (*it)->GetRegions();
      std::vector<G4String>::const_iterator itm;
      for (itm = regions.begin(); itm != regions.end(); itm++) {
        G4cout << "\'" << (*itm) << "\' ";
      }
      G4cout << G4endl;
    }
  }
#ifdef G4MULTITHREADED
  lm.unlock();
#endif
}
 
//
// public methods
//
 
//_____________________________________________________________________________
void TG4ModelConfigurationManager::CreateRegions()
{
 
  if (VerboseLevel() > 1) {
    G4cout << "TG4ModelConfigurationManager::CreateRegions" << G4endl;
  }
 
  // Return if no models are registered
  if (!fVector.size()) return;
 
  // Return if regions were already created
  if (fCreateRegionsDone) return;
 
  // Generate new regions names based on material names
  SetRegionsNames();
 
#ifdef G4MULTITHREADED
  G4AutoLock lm(&createRegionsMutex);
  if (!fCreateRegionsDone) {
#endif
    // Loop over logical volumes
    G4LogicalVolumeStore* lvStore = G4LogicalVolumeStore::GetInstance();
    for (G4int i = 0; i < G4int(lvStore->size()); i++) {
      G4LogicalVolume* lv = (*lvStore)[i];
      G4String materialName = lv->GetMaterial()->GetName();
 
      if (VerboseLevel() > 2) {
        G4cout << "Processing volume " << lv->GetName() << ", material "
               << materialName << G4endl;
      }
 
      // Skip volumes with materials which are not in the regions list
      G4bool isModelRegion = false;
      ModelConfigurationVector::const_iterator it;
      for (it = fVector.begin(); it != fVector.end(); it++) {
        if ((*it)->HasRegion(materialName)) isModelRegion = true;
      }
 
      if (!isModelRegion) {
        if (VerboseLevel() > 2) {
          G4cout << "Material " << materialName << " is not in selection"
                 << G4endl;
        }
        continue;
      }
 
      // If region already exists, only add the logical volume
      // and continue the loop
      G4Region* region =
        G4RegionStore::GetInstance()->GetRegion(materialName, false);
 
      if (!region) {
        region = new G4Region(materialName);
 
        if (VerboseLevel() > 1) {
          G4cout << "Created region " << region->GetName() << region << G4endl;
        }
      }
 
      if (VerboseLevel() > 1) {
        G4cout << "Adding region " << region->GetName() << " to lv "
               << lv->GetName() << G4endl;
      }
      region->AddRootLogicalVolume(lv);
      lv->SetRegion(region);
    }
    fCreateRegionsDone = true;
#ifdef G4MULTITHREADED
    lm.unlock();
  }
#endif
}
 
//_____________________________________________________________________________
void TG4ModelConfigurationManager::SetModel(const G4String& modelName)
{
 
  if (VerboseLevel() > 1) {
    G4cout << "TG4ModelConfigurationManager::SetModel: " << modelName << G4endl;
  }
 
  if (!GetModelConfiguration(modelName, false)) {
    fVector.push_back(new TG4ModelConfiguration(modelName));
  }
  else {
    TString text = "Cannot create model ";
    text += modelName.data();
    text += " twice.";
    TG4Globals::Warning("TG4ModelConfigurationManager", "SetModel", text);
  }
}
 
//_____________________________________________________________________________
void TG4ModelConfigurationManager::SetModelParticles(
  const G4String& modelName, const G4String& particles)
{
 
  TG4ModelConfiguration* modelConfiguration =
    GetModelConfiguration(modelName, "SetModelParticles");
 
  if (!modelConfiguration) {
    TString text = "The model configuration ";
    text += modelName.data();
    text += " is not defined.";
    TG4Globals::Warning("TG4ModelConfigurationManager", "SetModelParticles",
      text + TG4Globals::Endl() + TString("Setting will be ignored."));
    return;
  }
 
  modelConfiguration->SetParticles(particles);
}
 
//_____________________________________________________________________________
void TG4ModelConfigurationManager::SetModelRegions(
  const G4String& modelName, const G4String& regionsMedia)
{
 
  TG4ModelConfiguration* modelConfiguration =
    GetModelConfiguration(modelName, "SetModelRegions");
 
  if (!modelConfiguration) {
    TString text = "The model configuration ";
    text += modelName.data();
    text += " is not defined.";
    TG4Globals::Warning("TG4ModelConfigurationManager", "SetModelRegions",
      text + TG4Globals::Endl() + TString("Setting will be ignored."));
    return;
  }
 
  modelConfiguration->SetRegionsMedia(regionsMedia);
}
 
//_____________________________________________________________________________
void TG4ModelConfigurationManager::SetOneModelRegion(
  const G4String& modelName, const G4String& regionMedium)
{
 
  TG4ModelConfiguration* modelConfiguration =
    GetModelConfiguration(modelName, "SetModelRegions");
 
  if (!modelConfiguration) {
    TString text = "The model configuration ";
    text += modelName.data();
    text += " is not defined.";
    TG4Globals::Warning("TG4ModelConfigurationManager", "SetModelRegions",
      text + TG4Globals::Endl() + TString("Setting will be ignored."));
    return;
  }
  modelConfiguration->SetOneRegionMedium(regionMedium);
}
 
//_____________________________________________________________________________
TG4ModelConfiguration* TG4ModelConfigurationManager::GetModelConfiguration(
  const G4String& modelName, G4bool warn) const
{
 
  ModelConfigurationVector::const_iterator it;
  for (it = fVector.begin(); it != fVector.end(); it++) {
    if ((*it)->GetModelName() == modelName) return *it;
  }
 
  if (warn) {
    TString text = "Model configuration ";
    text += modelName.data();
    text += " does not exist";
    TG4Globals::Warning(
      "TG4ModelConfigurationManager", "GetModelConfiguration", text);
  }
 
  return 0;
}