TG4RegionsManager.h

Go to the documentation of this file.

#ifndef TG4_REGIONS_MANAGER_H
#define TG4_REGIONS_MANAGER_H
 
//------------------------------------------------
// The Geant4 Virtual Monte Carlo package
// Copyright (C) 2007 - 2014 Ivana Hrivnacova
// All rights reserved.
//
// For the licensing terms see geant4_vmc/LICENSE.
// Contact: root-vmc@cern.ch
//-------------------------------------------------
 
 
#include "TG4VRegionsManager.h"
 
#include <globals.hh>
 
#include <map>
 
class TG4Limits;
 
class G4Region;
class G4Material;
class G4VRangeToEnergyConverter;
class G4MaterialCutsCouple;
 
 
class TG4RegionsManager : public TG4VRegionsManager
{
 public:
  TG4RegionsManager();
  ~TG4RegionsManager() override = default;
 
  // methods
  void DefineRegions() override;;
  void UpdateProductionCutsTable() override {}
 
  void CheckRegions() const override;
  void PrintRegions(std::ostream& output) const override;
  void LoadRegions();
  void DumpRegion(const G4String& volName) const;
 
  // set methods
  void SetRangePrecision(G4int precision);
  void SetEnergyTolerance(G4double tolerance);
  void SetLoad(G4bool isLoad);
  void SetFromG4Table(G4bool isG4Table);
 
  // get methods
  G4int GetRangePrecision() const;
  G4double GetEnergyTolerance() const;
  G4String GetFileName() const;
  G4bool IsG4Table() const;
  G4bool IsLoad() const;
 
 private:
  using TG4RegionData = std::array<G4double, fgkValuesSize>;
 
  TG4RegionsManager(const TG4RegionsManager& right) = delete;
  TG4RegionsManager& operator=(const TG4RegionsManager& right) = delete;
 
  // methods
  G4bool Iterate(G4double energyCut, G4double& lowerCut, G4double& higherCut,
    G4double defaultRangeCut, G4double lowEdgeEnergy, G4double highEdgeEnergy,
    G4int nbin, std::map<G4double, G4double>& map,
    std::map<G4double, G4double>::const_iterator& it, G4Material* material,
    G4VRangeToEnergyConverter& converter) const;
 
  std::pair<G4double,G4double>
    ConvertEnergyToRange(G4double energyCut, G4Material* material,
    G4VRangeToEnergyConverter& converter, G4double defaultRangeValue) const;
 
  std::pair<G4double,G4double>
    GetRangeCut(G4double energyCut, G4Material* material,
    G4VRangeToEnergyConverter& converter, G4double defaultRangeValue) const;
 
  void CheckRegionsRanges() const;
  void PrintFromMap(std::ostream& output) const;
 
  //
  // static data members
 
  static constexpr G4int fgkDefaultRangePrecision = 5;
  static constexpr G4double fgkDefaultEnergyTolerance = 0.01;
  static constexpr G4int fgkNofBins = 10;
  static constexpr G4int fgkMinRangeOrder = -3;
  static constexpr G4int fgkMaxRangeOrder = 6;
 
  //
  // data members
 
  TG4RegionsMessenger fMessenger;
  G4int fRangePrecision = fgkDefaultRangePrecision;
  G4double fEnergyTolerance = fgkDefaultEnergyTolerance;
  G4bool fIsG4Table = false;
  G4bool fIsLoad = false;
  std::map<G4String, TG4RegionData> fRegionData;
};
 
inline void TG4RegionsManager::SetRangePrecision(G4int precision)
{
  G4cout << "### New precision: " << precision << G4endl;
  fRangePrecision = precision;
}
 
inline void TG4RegionsManager::SetEnergyTolerance(G4double tolerance)
{
  G4cout << "### New tolerance: " << tolerance << G4endl;
  fEnergyTolerance = tolerance;
}
 
inline void TG4RegionsManager::SetFromG4Table(G4bool isG4Table)
{
  fIsG4Table = isG4Table;
}
 
inline G4int TG4RegionsManager::GetRangePrecision() const
{
  return fRangePrecision;
}
 
inline G4double TG4RegionsManager::GetEnergyTolerance() const
{
  return fEnergyTolerance;
}
 
inline G4bool TG4RegionsManager::IsG4Table() const
{ return fIsG4Table; }
 
inline G4bool TG4RegionsManager::IsLoad() const { return fIsLoad; }
 
#endif // TG4_REGIONS_MANAGER_H