G4LatticeLogical Class Reference

#include <G4LatticeLogical.hh>

Public Types
enum	{ MAXRES = 322 }

Public Member Functions
	G4LatticeLogical ()
virtual	~G4LatticeLogical ()=default
void	SetVerboseLevel (G4int vb)
G4bool	LoadMap (G4int, G4int, G4int, G4String)
G4bool	Load_NMap (G4int, G4int, G4int, G4String)
virtual G4double	MapKtoV (G4int, const G4ThreeVector &) const
virtual G4ThreeVector	MapKtoVDir (G4int, const G4ThreeVector &) const
void	Dump (std::ostream &os) const
void	DumpMap (std::ostream &os, G4int pol, const G4String &name) const
void	Dump_NMap (std::ostream &os, G4int pol, const G4String &name) const
void	SetDynamicalConstants (G4double Beta, G4double Gamma, G4double Lambda, G4double Mu)
void	SetScatteringConstant (G4double b)
void	SetAnhDecConstant (G4double a)
void	SetLDOS (G4double LDOS)
void	SetSTDOS (G4double STDOS)
void	SetFTDOS (G4double FTDOS)
G4double	GetBeta () const
G4double	GetGamma () const
G4double	GetLambda () const
G4double	GetMu () const
G4double	GetScatteringConstant () const
G4double	GetAnhDecConstant () const
G4double	GetLDOS () const
G4double	GetSTDOS () const
G4double	GetFTDOS () const

Detailed Description

Definition at line 40 of file G4LatticeLogical.hh.

Member Enumeration Documentation

anonymous enum

anonymous enum

Enumerator
MAXRES

Definition at line 43 of file G4LatticeLogical.hh.

  {
    MAXRES = 322
  };  // Maximum map resolution (bins)

Constructor & Destructor Documentation

G4LatticeLogical()

G4LatticeLogical::G4LatticeLogical

(

)

Definition at line 38 of file G4LatticeLogical.cc.

{
  for (G4int i = 0; i < 3; i++) {
    for (G4int j = 0; j < MAXRES; j++) {
      for (G4int k = 0; k < MAXRES; k++) {
        fMap[i][j][k] = 0.;
        fN_map[i][j][k].set(0., 0., 0.);
      }
    }
  }
}

~G4LatticeLogical()

virtual G4LatticeLogical::~G4LatticeLogical ( )

virtualdefault

Member Function Documentation

Dump()

void G4LatticeLogical::Dump

(

std::ostream &

os

)

const

Definition at line 199 of file G4LatticeLogical.cc.

{
  os << "dyn " << fBeta << " " << fGamma << " " << fLambda << " " << fMu << "\nscat " << fB
     << " decay " << fA << "\nLDOS " << fLDOS << " STDOS " << fSTDOS << " FTDOS " << fFTDOS
     << std::endl;
 
  Dump_NMap(os, 0, "LVec.ssv");
  Dump_NMap(os, 1, "FTVec.ssv");
  Dump_NMap(os, 2, "STVec.ssv");
 
  DumpMap(os, 0, "L.ssv");
  DumpMap(os, 1, "FT.ssv");
  DumpMap(os, 2, "ST.ssv");
}

Referenced by operator<<().

Dump_NMap()

void G4LatticeLogical::Dump_NMap	(	std::ostream &	os,
		G4int	pol,
		const G4String &	name ) const

Definition at line 230 of file G4LatticeLogical.cc.

{
  os << "VDir " << name << " "
     << (pol == 0    ? "L"
          : pol == 1 ? "FT"
          : pol == 2 ? "ST"
                     : "??")
     << " " << fDresTheta << " " << fDresPhi << std::endl;
 
  for (G4int iTheta = 0; iTheta < fDresTheta; iTheta++) {
    for (G4int iPhi = 0; iPhi < fDresPhi; iPhi++) {
      os << fN_map[pol][iTheta][iPhi].x() << " " << fN_map[pol][iTheta][iPhi].y() << " "
         << fN_map[pol][iTheta][iPhi].z() << std::endl;
    }
  }
}

Referenced by Dump().

DumpMap()

void G4LatticeLogical::DumpMap	(	std::ostream &	os,
		G4int	pol,
		const G4String &	name ) const

Definition at line 214 of file G4LatticeLogical.cc.

{
  os << "VG " << name << " "
     << (pol == 0    ? "L"
          : pol == 1 ? "FT"
          : pol == 2 ? "ST"
                     : "??")
     << " " << fVresTheta << " " << fVresPhi << std::endl;
 
  for (G4int iTheta = 0; iTheta < fVresTheta; iTheta++) {
    for (G4int iPhi = 0; iPhi < fVresPhi; iPhi++) {
      os << fMap[pol][iTheta][iPhi] << std::endl;
    }
  }
}

Referenced by Dump().

GetAnhDecConstant()

G4double G4LatticeLogical::GetAnhDecConstant ( ) const

inline

Definition at line 87 of file G4LatticeLogical.hh.

{ return fA; }

Referenced by G4LatticePhysical::GetAnhDecConstant().

GetBeta()

G4double G4LatticeLogical::GetBeta ( ) const

inline

Definition at line 82 of file G4LatticeLogical.hh.

{ return fBeta; }

Referenced by G4LatticePhysical::GetBeta().

GetFTDOS()

G4double G4LatticeLogical::GetFTDOS ( ) const

inline

Definition at line 90 of file G4LatticeLogical.hh.

{ return fFTDOS; }

Referenced by G4LatticePhysical::GetFTDOS().

GetGamma()

G4double G4LatticeLogical::GetGamma ( ) const

inline

Definition at line 83 of file G4LatticeLogical.hh.

{ return fGamma; }

Referenced by G4LatticePhysical::GetGamma().

GetLambda()

G4double G4LatticeLogical::GetLambda ( ) const

inline

Definition at line 84 of file G4LatticeLogical.hh.

{ return fLambda; }

Referenced by G4LatticePhysical::GetLambda().

GetLDOS()

G4double G4LatticeLogical::GetLDOS ( ) const

inline

Definition at line 88 of file G4LatticeLogical.hh.

{ return fLDOS; }

Referenced by G4LatticePhysical::GetLDOS().

GetMu()

G4double G4LatticeLogical::GetMu ( ) const

inline

Definition at line 85 of file G4LatticeLogical.hh.

{ return fMu; }

Referenced by G4LatticePhysical::GetMu().

GetScatteringConstant()

G4double G4LatticeLogical::GetScatteringConstant ( ) const

inline

Definition at line 86 of file G4LatticeLogical.hh.

{ return fB; }

Referenced by G4LatticePhysical::GetScatteringConstant().

GetSTDOS()

G4double G4LatticeLogical::GetSTDOS ( ) const

inline

Definition at line 89 of file G4LatticeLogical.hh.

{ return fSTDOS; }

Referenced by G4LatticePhysical::GetSTDOS().

Load_NMap()

G4bool G4LatticeLogical::Load_NMap	(	G4int	tRes,
		G4int	pRes,
		G4int	polarizationState,
		G4String	map )

Definition at line 89 of file G4LatticeLogical.cc.

{
  if (tRes > MAXRES || pRes > MAXRES) {
    G4cerr << "G4LatticeLogical::LoadMap exceeds maximum resolution of " << MAXRES << " by "
           << MAXRES << ". terminating." << G4endl;
    return false;  // terminate if resolution out of bounds.
  }
 
  std::ifstream fMapFile(map.data());
  if (! fMapFile.is_open()) {
    return false;
  }
 
  G4double x, y, z;  // Buffers to read coordinates from file
  G4ThreeVector dir;
  for (G4int theta = 0; theta < tRes; theta++) {
    for (G4int phi = 0; phi < pRes; phi++) {
      fMapFile >> x >> y >> z;
      dir.set(x, y, z);
      fN_map[polarizationState][theta][phi] = dir.unit();  // Enforce unity
    }
  }
 
  if (verboseLevel != 0) {
    G4cout << "\nG4LatticeLogical::Load_NMap(" << map << ") successful"
           << " (Vdir " << tRes << " x " << pRes << " for polarization " << polarizationState
           << ")." << G4endl;
  }
 
  fDresTheta = tRes;  // store map dimensions
  fDresPhi = pRes;
  return true;
}

Referenced by G4LatticeReader::ProcessNMap().

LoadMap()

G4bool G4LatticeLogical::LoadMap	(	G4int	tRes,
		G4int	pRes,
		G4int	polarizationState,
		G4String	map )

Definition at line 53 of file G4LatticeLogical.cc.

{
  if (tRes > MAXRES || pRes > MAXRES) {
    G4cerr << "G4LatticeLogical::LoadMap exceeds maximum resolution of " << MAXRES << " by "
           << MAXRES << ". terminating." << G4endl;
    return false;  // terminate if resolution out of bounds.
  }
 
  std::ifstream fMapFile(map.data());
  if (! fMapFile.is_open()) {
    return false;
  }
 
  G4double vgrp = 0.;
  for (G4int theta = 0; theta < tRes; theta++) {
    for (G4int phi = 0; phi < pRes; phi++) {
      fMapFile >> vgrp;
      fMap[polarizationState][theta][phi] = vgrp * (m / s);
    }
  }
 
  if (verboseLevel != 0) {
    G4cout << "\nG4LatticeLogical::LoadMap(" << map << ") successful"
           << " (Vg scalars " << tRes << " x " << pRes << " for polarization " << polarizationState
           << ")." << G4endl;
  }
 
  fVresTheta = tRes;  // store map dimensions
  fVresPhi = pRes;
  return true;
}

Referenced by G4LatticeReader::ProcessMap().

MapKtoV()

G4double G4LatticeLogical::MapKtoV ( G4int polarizationState, const G4ThreeVector & k ) const

virtual

Definition at line 128 of file G4LatticeLogical.cc.

{
  G4double theta, phi, tRes, pRes;
 
  tRes = pi / fVresTheta;
  pRes = twopi / fVresPhi;
 
  theta = k.getTheta();
  phi = k.getPhi();
 
  if (phi < 0) {
    phi = phi + twopi;
  }
  if (theta > pi) {
    theta = theta - pi;
  }
 
  G4double Vg = fMap[polarizationState][int(theta / tRes)][int(phi / pRes)];
 
  if (Vg == 0) {
    G4cout << "\nFound v=0 for polarization " << polarizationState << " theta " << theta << " phi "
           << phi << " translating to map coords "
           << "theta " << int(theta / tRes) << " phi " << int(phi / pRes) << G4endl;
  }
 
  if (verboseLevel > 1) {
    G4cout << "G4LatticeLogical::MapKtoV theta,phi=" << theta << " " << phi << " : ith,iph "
           << int(theta / tRes) << " " << int(phi / pRes) << " : V " << Vg << G4endl;
  }
 
  return Vg;
}

Referenced by G4LatticePhysical::MapKtoV().

MapKtoVDir()

G4ThreeVector G4LatticeLogical::MapKtoVDir ( G4int polarizationState, const G4ThreeVector & k ) const

virtual

Definition at line 166 of file G4LatticeLogical.cc.

{
  G4double theta, phi, tRes, pRes;
 
  tRes = pi / (fDresTheta - 1);  // The summant "-1" is required:index=[0:array length-1]
  pRes = 2 * pi / (fDresPhi - 1);
 
  theta = k.getTheta();
  phi = k.getPhi();
 
  if (theta > pi) {
    theta = theta - pi;
  }
  // phi=[0 to 2 pi] in accordance with DMC //if(phi>pi/2) phi=phi-pi/2;
  if (phi < 0) {
    phi = phi + 2 * pi;
  }
 
  auto iTheta = int(theta / tRes + 0.5);
  auto iPhi = int(phi / pRes + 0.5);
 
  if (verboseLevel > 1) {
    G4cout << "G4LatticeLogical::MapKtoVDir theta,phi=" << theta << " " << phi << " : ith,iph "
           << iTheta << " " << iPhi << " : dir " << fN_map[polarizationState][iTheta][iPhi]
           << G4endl;
  }
 
  return fN_map[polarizationState][iTheta][iPhi];
}

Referenced by G4LatticePhysical::MapKtoVDir().

SetAnhDecConstant()

void G4LatticeLogical::SetAnhDecConstant ( G4double a )

inline

Definition at line 77 of file G4LatticeLogical.hh.

{ fA = a; }

Referenced by G4LatticeReader::ProcessValue().

SetDynamicalConstants()

void G4LatticeLogical::SetDynamicalConstants ( G4double Beta, G4double Gamma, G4double Lambda, G4double Mu )

inline

Definition at line 68 of file G4LatticeLogical.hh.

  {
    fBeta = Beta;
    fGamma = Gamma;
    fLambda = Lambda;
    fMu = Mu;
  }

Referenced by G4LatticeReader::ProcessConstants().

SetFTDOS()

void G4LatticeLogical::SetFTDOS ( G4double FTDOS )

inline

Definition at line 80 of file G4LatticeLogical.hh.

{ fFTDOS = FTDOS; }

Referenced by G4LatticeReader::ProcessValue().

SetLDOS()

void G4LatticeLogical::SetLDOS ( G4double LDOS )

inline

Definition at line 78 of file G4LatticeLogical.hh.

{ fLDOS = LDOS; }

Referenced by G4LatticeReader::ProcessValue().

SetScatteringConstant()

void G4LatticeLogical::SetScatteringConstant ( G4double b )

inline

Definition at line 76 of file G4LatticeLogical.hh.

{ fB = b; }

Referenced by G4LatticeReader::ProcessValue().

SetSTDOS()

void G4LatticeLogical::SetSTDOS ( G4double STDOS )

inline

Definition at line 79 of file G4LatticeLogical.hh.

{ fSTDOS = STDOS; }

Referenced by G4LatticeReader::ProcessValue().

SetVerboseLevel()

void G4LatticeLogical::SetVerboseLevel ( G4int vb )

inline

Definition at line 52 of file G4LatticeLogical.hh.

{ verboseLevel = vb; }

---

The documentation for this class was generated from the following files:

• G4LatticeLogical.hh
• G4LatticeLogical.cc