TG4GeometryServices.h

Go to the documentation of this file.

#ifndef TG4_GEOMETRY_SERVICES_H
#define TG4_GEOMETRY_SERVICES_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 "TG4Globals.h"
#include "TG4MediumMap.h"
#include "TG4OpSurfaceMap.h"
#include "TG4Verbose.h"
 
#include <G4OpticalSurface.hh>
#include <G4SurfaceProperty.hh>
#include <G4Transform3D.hh>
#include <G4Version.hh>
#include <globals.hh>
 
#include <Rtypes.h>
#include <TMCOptical.h>
 
#include <map>
#include <vector>
 
class TG4MediumMap;
class TG4NameMap;
class TG4Limits;
class TG4G3CutVector;
class TG4G3ControlVector;
 
class G4Material;
class G4LogicalVolume;
class G4VPhysicalVolume;
class G4UserLimits;
class G4OpticalSurface;
 
class TGeoHMatrix;

struct TG4AssemblyLevels
{
  std::vector<G4String> fNames; 
  std::vector<G4int> fCopyNos;  
};

 
class TG4GeometryServices : public TG4Verbose
{
 public:
  TG4GeometryServices();
  virtual ~TG4GeometryServices();
 
  // static access method
  static TG4GeometryServices* Instance();
 
  // methods
  void BuildAssemblyLevels();
  // Snapshot reflection provenance on the master for worker name lookups.
  void BuildConstituentVolumes();
  const TG4AssemblyLevels& GetAssemblyLevels(const G4VPhysicalVolume* pv) const;
  // utilities
  G4double* CreateG4doubleArray(Float_t* array, G4int size, G4bool copyValues = true) const;
  G4double* CreateG4doubleArray(Double_t* array, G4int size, G4bool copyValues = true) const;
  G4String CutName(const char* name) const;
  G4String CutMaterialName(const char* name) const;
  G4String CutVolumePath(
    const G4String& volumePath, G4String& volName, G4int& copyNo) const;
  const G4String& UserVolumeName(const G4String& name) const;
  // Reflection provenance, without guessing from a name suffix.
  const G4String& GetConstituentVolumeName(G4LogicalVolume* lv) const;
  // VMC source name: constituent name followed by Gsposp normalization.
  G4String UserVolumeName(G4LogicalVolume* lv) const;
 
  G4OpticalSurfaceModel SurfaceModel(EMCOpSurfaceModel model) const;
  G4SurfaceType SurfaceType(EMCOpSurfaceType surfType) const;
  G4OpticalSurfaceFinish SurfaceFinish(EMCOpSurfaceFinish finish) const;
  void Convert(const G4Transform3D& transform, TGeoHMatrix& matrix) const;
 
  G4Material* MixMaterials(G4String name, G4double density,
    const TG4StringVector& matNames, const TG4doubleVector& matWeights);
  // printing
  void PrintLimits(const G4String& name) const;
  void PrintVolumeLimits(const G4String& volumeName) const;
  void PrintStatistics(G4bool open, G4bool close) const;
  void PrintLogicalVolumeStore() const;
  void PrintPhysicalVolumeStore() const;
  void PrintElementTable() const;
  void PrintMaterials() const;
  void PrintMaterialsProperties() const;
  void PrintMedia() const;
  void PrintCuts(const G4String& cutName) const;
  void PrintControls(const G4String& controlName) const;
 
  // set methods
  void SetWorld(G4VPhysicalVolume* world);
  void SetIsG3toG4(G4bool isG3toG4);
  void SetG3toG4Separator(char separator);
  void SetAccountAssemblyLevels(G4bool option);
 
  // get methods
  // volumes
  Int_t NofG3Volumes() const;
  Int_t NofG4LogicalVolumes() const;
  Int_t NofG4PhysicalVolumes() const;
  G4VPhysicalVolume* GetWorld() const;
 
  TG4Limits* GetLimits(G4UserLimits* limits) const;
  TG4Limits* GetLimits(G4UserLimits* limits, const TG4G3CutVector& cuts,
    const TG4G3ControlVector& controls) const;
 
  // All Geant4 representations of a VMC source volume.
  std::vector<G4LogicalVolume*> FindLogicalVolumes(
    const G4String& name, G4bool silent = false) const;
  // Representative only; use the plural form for assignments.
  G4LogicalVolume* FindLogicalVolume(
    const G4String& name, G4bool silent = false) const;
  G4VPhysicalVolume* FindPhysicalVolume(
    const G4String& name, G4int copyNo, G4bool silent = false) const;
  G4VPhysicalVolume* FindDaughter(const G4String& name, G4int copyNo,
    G4LogicalVolume* mlv, G4bool silent = false) const;
  TG4Limits* FindLimits(const G4String& name, G4bool silent = false) const;
  TG4Limits* FindLimits2(const G4String& name, G4bool silent = false) const;
  TG4Limits* FindLimits(const G4Material*, G4bool silent = false) const;
 
  // materials
  G4int GetMediumId(G4LogicalVolume* lv) const;
  G4double GetEffA(G4Material* material) const;
  G4double GetEffZ(G4Material* material) const;
  G4Material* FindMaterial(G4double a, G4double z, G4double density) const;
  G4Material* FindMaterial(G4double* a, G4double* z, G4double density,
    G4int nmat, G4double* wmat) const;
 
  TG4MediumMap* GetMediumMap() const;
  TG4OpSurfaceMap* GetOpSurfaceMap() const;
 
 private:
  TG4GeometryServices(const TG4GeometryServices& right);
  TG4GeometryServices& operator=(const TG4GeometryServices& right);
 
  // methods
  G4bool IsG3Volume(const G4String& lvName) const;
  G4bool CompareElement(G4double a, G4double z, const G4Element* elem) const;
  G4bool CompareMaterial(
    G4int nofElements, G4double density, const G4Material* material) const;
  G4double* ConvertAtomWeight(G4int nmat, G4double* a, G4double* wmat) const;
 
  // static data members
  static TG4GeometryServices* fgInstance;    
  static thread_local G4String fgBuffer;                  
  static const G4double fgkAZTolerance;      
  static const G4double fgkDensityTolerance; 
 
  //
  // data members

  G4bool fIsG3toG4;

  std::map<G4LogicalVolume*, G4LogicalVolume*> fConstituentVolumes;

  TG4MediumMap* fMediumMap;

  TG4OpSurfaceMap* fOpSurfaceMap;

  G4VPhysicalVolume* fWorld;

  G4bool fAccountAssemblyLevels;

  std::vector<TG4AssemblyLevels> fAssemblyLevels;

  G4bool fAssemblyLevelsBuilt;
};
 
// inline methods
 
inline TG4GeometryServices* TG4GeometryServices::Instance()
{
  return fgInstance;
}
 
inline void TG4GeometryServices::SetWorld(G4VPhysicalVolume* world)
{
  fWorld = world;
}
 
inline void TG4GeometryServices::SetIsG3toG4(G4bool isG3toG4)
{
  fIsG3toG4 = isG3toG4;
}
 
inline void TG4GeometryServices::SetAccountAssemblyLevels(G4bool option)
{
  fAccountAssemblyLevels = option;
}
 
inline G4VPhysicalVolume* TG4GeometryServices::GetWorld() const
{
  return fWorld;
}
 
inline TG4MediumMap* TG4GeometryServices::GetMediumMap() const
{
  return fMediumMap;
}
 
inline TG4OpSurfaceMap* TG4GeometryServices::GetOpSurfaceMap() const
{
  return fOpSurfaceMap;
}
 
#endif // TG4_GEOMETRY_SERVICES_H