XLogicalLattice Class Reference

#include <XLogicalLattice.hh>

Collaboration diagram for XLogicalLattice:

Public Member Functions
	XLogicalLattice ()
	~XLogicalLattice ()
void	SetDynamicalConstants (double, double, double, double)
void	SetScatteringConstant (G4double)
void	SetAnhDecConstant (G4double)
void	SetLDOS (double)
void	SetSTDOS (double)
void	SetFTDOS (double)
double	GetBeta ()
double	GetGamma ()
double	GetLambda ()
double	GetMu ()
G4double	GetScatteringConstant ()
G4double	GetAnhDecConstant ()
double	GetLDOS ()
double	GetSTDOS ()
double	GetFTDOS ()
bool	LoadMap (int, int, int, string)
bool	Load_NMap (int, int, int, string)
double	MapKtoV (int, G4ThreeVector)
G4ThreeVector	MapKtoVDir (int, G4ThreeVector)

Private Attributes
int	fVresTheta
int	fVresPhi
int	fDresTheta
int	fDresPhi
double	fMap [3][MAXRES][MAXRES]
double	fN_map [3][MAXRES][MAXRES][3]
double	fA
double	fB
double	fDosL
double	fDosST
double	fDosFT
double	fBeta
double	fGamma
double	fLambda
double	fMu
ifstream	fMapFile
int	fThetaRes
int	fPhiRes

Detailed Description

Definition at line 40 of file XLogicalLattice.hh.

Constructor & Destructor Documentation

XLogicalLattice()

XLogicalLattice::XLogicalLattice

(

)

Definition at line 32 of file XLogicalLattice.cc.

                                {
  fVresTheta=0;
  fVresPhi=0;
  fDresTheta=0;
  fDresPhi=0;
  fA=0;
  fB=0;
  fDosL=0;
  fDosST=0;
  fDosFT=0;
  fBeta=0;
  fLambda=0;
  fGamma=0;
  fMu=0;
}

~XLogicalLattice()

XLogicalLattice::~XLogicalLattice

(

)

Definition at line 51 of file XLogicalLattice.cc.

{;}

Member Function Documentation

GetAnhDecConstant()

G4double XLogicalLattice::GetAnhDecConstant

(

)

Definition at line 155 of file XLogicalLattice.cc.

{
  return fA;
}

Here is the caller graph for this function:

GetBeta()

double XLogicalLattice::GetBeta

(

)

Definition at line 115 of file XLogicalLattice.cc.

{
  return fBeta;
}

Here is the caller graph for this function:

GetFTDOS()

double XLogicalLattice::GetFTDOS

(

)

Definition at line 179 of file XLogicalLattice.cc.

{
  return fDosFT;
}

Here is the caller graph for this function:

GetGamma()

double XLogicalLattice::GetGamma

(

)

Definition at line 123 of file XLogicalLattice.cc.

{
  return fGamma;
}

Here is the caller graph for this function:

GetLambda()

double XLogicalLattice::GetLambda

(

)

Definition at line 131 of file XLogicalLattice.cc.

{
  return fLambda;
}

Here is the caller graph for this function:

GetLDOS()

double XLogicalLattice::GetLDOS

(

)

Definition at line 163 of file XLogicalLattice.cc.

{
  return fDosL;
}

Here is the caller graph for this function:

GetMu()

double XLogicalLattice::GetMu

(

)

Definition at line 139 of file XLogicalLattice.cc.

{
  return fMu;
}

Here is the caller graph for this function:

GetScatteringConstant()

G4double XLogicalLattice::GetScatteringConstant

(

)

Definition at line 147 of file XLogicalLattice.cc.

{
  return fB;
}

Here is the caller graph for this function:

GetSTDOS()

double XLogicalLattice::GetSTDOS

(

)

Definition at line 171 of file XLogicalLattice.cc.

{
  return fDosST;
}

Here is the caller graph for this function:

Load_NMap()

bool XLogicalLattice::Load_NMap	(	int	tRes,
		int	pRes,
		int	polarizationState,
		string	m
	)

Definition at line 227 of file XLogicalLattice.cc.

{
  if((tRes>MAXRES)||(pRes>MAXRES)){
    G4cout<<"\nk-V map exceeds maximum resolution of "<<MAXRES<<
      " by "<<MAXRES<<" terminating."<<endl;
    return false; //terminate if resolution out of bounds.
  }


  fMapFile.clear();

  fThetaRes=tRes;
  fPhiRes=pRes;
  fMapFile.open(m.data());
  if(!fMapFile.is_open()) return false;
  
  for(int theta = 0; theta<fThetaRes; theta++){
    for(int phi = 0; phi<fPhiRes; phi++){
      for(int coord = 0; coord<3; coord++){
        fMapFile>>fN_map[polarizationState][theta][phi][coord];
      }
    }
  }
  G4cout<<"\n";
  fMapFile.close();
  G4cout<<"\nXLogicalLattice::Load_NMap() sucessful"<<endl;
  fDresTheta=tRes; //store map dimesnions
  fDresPhi=pRes;
  return true;
}

LoadMap()

bool XLogicalLattice::LoadMap	(	int	tRes,
		int	pRes,
		int	polarizationState,
		string	m
	)

Definition at line 190 of file XLogicalLattice.cc.

{
  if((tRes>MAXRES)||(pRes>MAXRES)){
    G4cout<<"\nk-V fMap exceeds maximum resolution of "<<MAXRES<<
      " by "<<MAXRES<<". terminating."<<endl;
    return false; //terminate if resolution out of bounds.
  }


  fMapFile.clear();

  fThetaRes=tRes;
  fPhiRes=pRes;
  fMapFile.open(m.data());
  if(!fMapFile.is_open()) return false;
  
  for(int theta = 0; theta<fThetaRes; theta++){
    for(int phi = 0; phi<fPhiRes; phi++){
        fMapFile>>fMap[polarizationState][theta][phi];        
    }
  }
  fMapFile.close();
  G4cout<<"\nXLogicalLattice::LoadMap() sucessful (Vg scalars)."<<endl;
  fVresTheta=tRes; //store map dimensions
  fVresPhi=pRes;
  return true;
}

MapKtoV()

double XLogicalLattice::MapKtoV	(	int	polarizationState,
		G4ThreeVector	k
	)

Definition at line 265 of file XLogicalLattice.cc.

{
  //Given the phonon wave vector k, phonon physical volume Vol 
  //and polarizationState(0=LON, 1=FT, 2=ST), 
  //returns phonon velocity in m/s

  double theta, phi, tRes, pRes;

  tRes=pi/(fVresTheta);//The summant "-1" is required:index=[0:array length-1]
  pRes=2*pi/(fVresPhi);

  
  theta=k.getTheta();
  phi=k.getPhi();

  if(phi<0) phi = phi + 2*pi;
  if(theta>pi) theta=theta-pi;
  //phi=[0 to 2 pi] in accordance with DMC //if(phi>pi/2) phi=phi-pi/2;
  if(fMap[polarizationState][int(theta/tRes)][int(phi/pRes)]==0){
      G4cout<<"\nFound v=0 for polarization "<<polarizationState
      <<" theta "<<theta<<" phi "<<phi<< " translating to map coords "
      << "theta "<< int(theta/tRes) << " phi " << int(phi/pRes)<<endl;
  }

  return fMap[polarizationState][int(theta/tRes)][int(phi/pRes)];
  
}

Here is the call graph for this function:

Here is the caller graph for this function:

MapKtoVDir()

G4ThreeVector XLogicalLattice::MapKtoVDir	(	int	polarizationState,
		G4ThreeVector	k
	)

Definition at line 296 of file XLogicalLattice.cc.

{  
  //Given the phonon wave vector k, phonon physical volume Vol 
  //and polarizationState(0=LON, 1=FT, 2=ST), 
  //returns phonon propagation direction as dimensionless unit vector

  double 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;

  G4int iTheta = int(theta/tRes+0.5);
  G4int iPhi = int(phi/pRes+0.5);

  
  G4ThreeVector v(fN_map[polarizationState][iTheta][iPhi][0],
                  fN_map[polarizationState][iTheta][iPhi][1],
                  fN_map[polarizationState][iTheta][iPhi][2]);


  //v.Set(sin(theta)*cos(phi),sin(theta)*sin(phi),cos(theta));
  //v.Set(k.getX(), k.getY(), k.getZ());

  return v.unit();
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetAnhDecConstant()

void XLogicalLattice::SetAnhDecConstant

(

G4double

a

)

Definition at line 79 of file XLogicalLattice.cc.

{
  //Constant governing rate of anharmonic down conversion of L-phonons,
  //for use with XPhononDownconversionProcess
  fA=a;
}

SetDynamicalConstants()

void XLogicalLattice::SetDynamicalConstants	(	double	Beta,
		double	Gamma,
		double	Lambda,
		double	Mu
	)

Definition at line 54 of file XLogicalLattice.cc.

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

}

SetFTDOS()

void XLogicalLattice::SetFTDOS

(

double

FTDOS

)

Definition at line 106 of file XLogicalLattice.cc.

{
  //Fast transverse phonon density of states
  fDosFT=FTDOS;
}

SetLDOS()

void XLogicalLattice::SetLDOS

(

double

LDOS

)

Definition at line 88 of file XLogicalLattice.cc.

{
  //Longitudinal phonon density of states
  fDosL=LDOS;
}

SetScatteringConstant()

void XLogicalLattice::SetScatteringConstant

(

G4double

b

)

Definition at line 69 of file XLogicalLattice.cc.

{
  //Constant governing the rate of isotope scattering, for use with
  //XPhononScatteringProcess
  fB=b;
}

Here is the caller graph for this function:

SetSTDOS()

void XLogicalLattice::SetSTDOS

(

double

STDOS

)

Definition at line 97 of file XLogicalLattice.cc.

{
  //Slow transverse phonon density of states
  fDosST=STDOS;
}

Member Data Documentation

fA

double XLogicalLattice::fA

private

Definition at line 53 of file XLogicalLattice.hh.

fB

double XLogicalLattice::fB

private

Definition at line 54 of file XLogicalLattice.hh.

fBeta

double XLogicalLattice::fBeta

private

Definition at line 58 of file XLogicalLattice.hh.

fDosFT

double XLogicalLattice::fDosFT

private

Definition at line 57 of file XLogicalLattice.hh.

fDosL

double XLogicalLattice::fDosL

private

Definition at line 55 of file XLogicalLattice.hh.

fDosST

double XLogicalLattice::fDosST

private

Definition at line 56 of file XLogicalLattice.hh.

fDresPhi

int XLogicalLattice::fDresPhi

private

Definition at line 47 of file XLogicalLattice.hh.

fDresTheta

int XLogicalLattice::fDresTheta

private

Definition at line 46 of file XLogicalLattice.hh.

fGamma

double XLogicalLattice::fGamma

private

Definition at line 58 of file XLogicalLattice.hh.

fLambda

double XLogicalLattice::fLambda

private

Definition at line 58 of file XLogicalLattice.hh.

fMap

double XLogicalLattice::fMap[3][MAXRES][MAXRES]

private

Definition at line 49 of file XLogicalLattice.hh.

fMapFile

ifstream XLogicalLattice::fMapFile

private

Definition at line 61 of file XLogicalLattice.hh.

fMu

double XLogicalLattice::fMu

private

Definition at line 58 of file XLogicalLattice.hh.

fN_map

double XLogicalLattice::fN_map[3][MAXRES][MAXRES][3]

private

Definition at line 50 of file XLogicalLattice.hh.

fPhiRes

int XLogicalLattice::fPhiRes

private

Definition at line 62 of file XLogicalLattice.hh.

fThetaRes

int XLogicalLattice::fThetaRes

private

Definition at line 62 of file XLogicalLattice.hh.

fVresPhi

int XLogicalLattice::fVresPhi

private

Definition at line 45 of file XLogicalLattice.hh.

fVresTheta

int XLogicalLattice::fVresTheta

private

Definition at line 44 of file XLogicalLattice.hh.

---

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

• XLogicalLattice.hh
• XLogicalLattice.cc