XPhysicalLattice Class Reference

#include <XPhysicalLattice.hh>

Collaboration diagram for XPhysicalLattice:

Public Member Functions
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 ()
	XPhysicalLattice ()
	XPhysicalLattice (G4VPhysicalVolume *, XLogicalLattice *)
	~XPhysicalLattice ()
double	MapKtoV (int, G4ThreeVector)
G4ThreeVector	MapKtoVDir (int, G4ThreeVector)
G4VPhysicalVolume *	GetVolume ()
void	SetPhysicalVolume (G4VPhysicalVolume *)
void	SetXLogicalLattice (XLogicalLattice *)
void	SetLatticeOrientation (G4double, G4double)
void	SetLatticeOrientation (G4double, G4double, G4double)
void	SetMillerOrientation (int, int, int)
void	SetUnitCell (XUnitCell *)
G4ThreeVector	ProjectMomentumVectorFromWorldToLattice (G4ThreeVector &, G4ThreeVector &)
G4ThreeVector	ProjectMomentumVectorFromLatticeToWorld (G4ThreeVector &, G4ThreeVector &)
G4ThreeVector	GetLatticeDirection (G4ThreeVector &)
XUnitCell *	GetXUnitCell ()
XLogicalLattice *	GetLogicalLattice ()
G4int	GetMiller (G4int)
G4ThreeVector	GetLatticeAngles ()
G4ThreeVector	GetCurvatureRadius ()
void	SetCurvatureRadius (G4ThreeVector)
G4ThreeVector	ComputeBendingAngle (G4ThreeVector &)
G4bool	IsBent ()
G4double	ComputeInterplanarPeriod ()
void	SetThermalVibrationAmplitude (G4double)
G4double	GetThermalVibrationAmplitude ()

Public Attributes
G4AffineTransform	fLocalToGlobal
G4AffineTransform	fGlobalToLocal

Detailed Description

Definition at line 40 of file XPhysicalLattice.hh.

Constructor & Destructor Documentation

XPhysicalLattice::XPhysicalLattice

(

)

Definition at line 38 of file XPhysicalLattice.cc.

                                  {
    fLattice=NULL;
    fVolume=NULL;
    fTheta=0.;
    fPhi=0.;
    fOmega=0.;

    
    fCurvatureRadius = G4ThreeVector(0.,0.,0.); // if cr = 0 == no bending
    fThermalVibrationAmplitude = 0.1 * angstrom; // no physical meaning
    fMillerOrientation[0] = 2;
    fMillerOrientation[1] = 2;
    fMillerOrientation[2] = 0;
}

XPhysicalLattice::XPhysicalLattice	(	G4VPhysicalVolume *	Vol,
		XLogicalLattice *	Lat
	)

Definition at line 55 of file XPhysicalLattice.cc.

                                                        {
    fLattice=Lat;
    fVolume=Vol;
    fA=fLattice->GetAnhDecConstant();
    fB=fLattice->GetScatteringConstant();
    fDosL=fLattice->GetLDOS();
    fDosST=fLattice->GetSTDOS();
    fDosFT=fLattice->GetFTDOS();
    fBeta=fLattice->GetBeta();
    fGamma=fLattice->GetGamma();
    fLambda=fLattice->GetLambda();
    fMu=fLattice->GetMu();
    
    G4RotationMatrix *rot = fVolume->GetObjectRotation();
    
    fGlobalToLocal = G4AffineTransform(*rot);
    fLocalToGlobal = fGlobalToLocal.Invert();
}

Here is the call graph for this function:

XPhysicalLattice::~XPhysicalLattice

(

)

Definition at line 77 of file XPhysicalLattice.cc.

{;}

Member Function Documentation

G4ThreeVector XPhysicalLattice::ComputeBendingAngle

(

G4ThreeVector &

vPosition

)

Definition at line 379 of file XPhysicalLattice.cc.

                                                                           {

    G4double vAngleX = 0.;
    G4double vAngleY = 0.;
    
    if(GetCurvatureRadius().x() != 0){
        vAngleX = vPosition.phi();
    }
    
    return G4ThreeVector(vAngleX,0.,vAngleY);
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4double XPhysicalLattice::ComputeInterplanarPeriod

(

)

Definition at line 413 of file XPhysicalLattice.cc.

                                                   {
    G4double vInterplanarPeriod =
        GetXUnitCell()->ComputeDirectPeriod(GetMiller(0),
                                            GetMiller(1),
                                            GetMiller(2));
    return vInterplanarPeriod;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4double XPhysicalLattice::GetAnhDecConstant

(

)

Definition at line 165 of file XPhysicalLattice.cc.

{
    return fA;
}

double XPhysicalLattice::GetBeta

(

)

Definition at line 129 of file XPhysicalLattice.cc.

                                {
    return fBeta;
}

G4ThreeVector XPhysicalLattice::GetCurvatureRadius

(

)

Definition at line 362 of file XPhysicalLattice.cc.

                                                  {
    return fCurvatureRadius;
}

Here is the caller graph for this function:

double XPhysicalLattice::GetFTDOS

(

)

Definition at line 187 of file XPhysicalLattice.cc.

{
    return fDosFT;
}

double XPhysicalLattice::GetGamma

(

)

Definition at line 136 of file XPhysicalLattice.cc.

                                 {
    return fGamma;
}

double XPhysicalLattice::GetLambda

(

)

Definition at line 142 of file XPhysicalLattice.cc.

                                  {
    return fLambda;
}

G4ThreeVector XPhysicalLattice::GetLatticeAngles

(

)

Definition at line 407 of file XPhysicalLattice.cc.

                                                {
    return G4ThreeVector(fPhi,fTheta,fOmega);
}

G4ThreeVector XPhysicalLattice::GetLatticeDirection

(

G4ThreeVector &

vPosition

)

Definition at line 336 of file XPhysicalLattice.cc.

                                                                           {
        G4ThreeVector dir = G4ThreeVector(0.,0.,1.);
    return ProjectMomentumVectorFromLatticeToWorld(dir,
                                                   vPosition);
}

Here is the call graph for this function:

double XPhysicalLattice::GetLDOS

(

)

Definition at line 173 of file XPhysicalLattice.cc.

                                {
    return fDosL;
}

XLogicalLattice * XPhysicalLattice::GetLogicalLattice

(

)

Definition at line 393 of file XPhysicalLattice.cc.

                                                    {
    return fLattice;
}

Here is the caller graph for this function:

G4int XPhysicalLattice::GetMiller

(

G4int

vIndex

)

Definition at line 399 of file XPhysicalLattice.cc.

                                             {
    if(vIndex<3 && vIndex>=0)
        return fMillerOrientation[vIndex];
    else return -1;
}

Here is the caller graph for this function:

double XPhysicalLattice::GetMu

(

)

Definition at line 149 of file XPhysicalLattice.cc.

{
    return fMu;
}

G4double XPhysicalLattice::GetScatteringConstant

(

)

Definition at line 157 of file XPhysicalLattice.cc.

{
    return fB;
}

double XPhysicalLattice::GetSTDOS

(

)

Definition at line 180 of file XPhysicalLattice.cc.

                                 {
    return fDosST;
}

G4double XPhysicalLattice::GetThermalVibrationAmplitude

(

)

Definition at line 294 of file XPhysicalLattice.cc.

                                                       {
    return fThermalVibrationAmplitude;
}

Here is the caller graph for this function:

G4VPhysicalVolume * XPhysicalLattice::GetVolume

(

)

Definition at line 231 of file XPhysicalLattice.cc.

                                              {
    return fVolume;
}

Here is the caller graph for this function:

XUnitCell * XPhysicalLattice::GetXUnitCell

(

)

Definition at line 350 of file XPhysicalLattice.cc.

                                         {
    return fUnitCell;
}

Here is the caller graph for this function:

G4bool XPhysicalLattice::IsBent

(

)

Definition at line 368 of file XPhysicalLattice.cc.

                               {
    if(fCurvatureRadius.x() != 0.) {
        return true;
    }
    else {
        return false;
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 199 of file XPhysicalLattice.cc.

                                                                      {
    double groupVelocity;
    
    k.rotate(G4ThreeVector(0,1,0), fTheta).rotate(G4ThreeVector(0,0,1), fPhi);
    groupVelocity = fLattice->MapKtoV(polarizationState, k);
    k.rotate(G4ThreeVector(0,0,1), -fPhi).rotate(G4ThreeVector(0,1,0), -fTheta);
    
    return groupVelocity;
}

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 216 of file XPhysicalLattice.cc.

                                                           {
    
    G4ThreeVector GroupVelocity;
    
    k=k.rotate(G4ThreeVector(0,1,0), fTheta).rotate(G4ThreeVector(0,0,1), fPhi);
    GroupVelocity = fLattice->MapKtoVDir(polarizationState, k);
    
    return GroupVelocity.rotate(G4ThreeVector(0,0,1), -fPhi)
        .rotate(G4ThreeVector(0,1,0), -fTheta).unit();
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4ThreeVector XPhysicalLattice::ProjectMomentumVectorFromLatticeToWorld	(	G4ThreeVector &	vMomentum,
		G4ThreeVector &	vPosition
	)

Definition at line 319 of file XPhysicalLattice.cc.

                                                                 {
    vMomentum.rotate(G4ThreeVector(0.,0.,1.), -fPhi)
        .rotate(G4ThreeVector(0.,1.,0.), -fTheta)
        .rotate(G4ThreeVector(1.,0.,0.), fOmega);
    
    if(IsBent() ){
        G4ThreeVector vBendingAngle = ComputeBendingAngle(vPosition);
        vMomentum.rotate(G4ThreeVector(0.,1.,0.), -vBendingAngle.x())
            .rotate(G4ThreeVector(1.,0.,0.), -vBendingAngle.z());
    }

    return vMomentum;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4ThreeVector XPhysicalLattice::ProjectMomentumVectorFromWorldToLattice	(	G4ThreeVector &	vMomentum,
		G4ThreeVector &	vPosition
	)

Definition at line 301 of file XPhysicalLattice.cc.

                                                                 {
    vMomentum.rotate(G4ThreeVector(1.,0.,0.),fOmega)
        .rotate(G4ThreeVector(0.,1.,0.), fTheta)
        .rotate(G4ThreeVector(0.,0.,1.), fPhi);

    if(IsBent() ){
        G4ThreeVector vBendingAngle = ComputeBendingAngle(vPosition);
        vMomentum.rotate(G4ThreeVector(1.,0.,0.), vBendingAngle.z())
            .rotate(G4ThreeVector(0.,1.,0.),vBendingAngle.x());
    }

    return vMomentum;
}

Here is the call graph for this function:

void XPhysicalLattice::SetAnhDecConstant

(

G4double

b

)

Definition at line 102 of file XPhysicalLattice.cc.

                                                  {
    fB=b;
}

void XPhysicalLattice::SetCurvatureRadius

(

G4ThreeVector

cr

)

Definition at line 356 of file XPhysicalLattice.cc.

                                                         {
    fCurvatureRadius =  cr;
}

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

Definition at line 81 of file XPhysicalLattice.cc.

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

void XPhysicalLattice::SetFTDOS

(

double

FTDOS

)

Definition at line 122 of file XPhysicalLattice.cc.

                                           {
    fDosFT = FTDOS;
}

void XPhysicalLattice::SetLatticeOrientation	(	G4double	t_rot,
		G4double	p_rot
	)

Definition at line 245 of file XPhysicalLattice.cc.

                                                                          {
    fTheta=t_rot;
    fPhi= p_rot;
}

void XPhysicalLattice::SetLatticeOrientation	(	G4double	t_rot,
		G4double	o_rot,
		G4double	p_rot
	)

Definition at line 252 of file XPhysicalLattice.cc.

                                                            {
    fTheta = t_rot;
    fOmega = o_rot;
    fPhi = p_rot;
}

void XPhysicalLattice::SetLDOS

(

double

LDOS

)

Definition at line 109 of file XPhysicalLattice.cc.

                                         {
    fDosL=LDOS;
}

void XPhysicalLattice::SetMillerOrientation	(	int	l,
		int	k,
		int	n
	)

Definition at line 262 of file XPhysicalLattice.cc.

                                                              {
    //fTheta=pi/2-std::atan2(n+0.000001,l+0.000001)*rad;
    //fPhi=pi/2-std::atan2(l+0.000001,k+0.000001)*rad;
    
    // // // added for channeling // // //
    fMillerOrientation[0]=l;
    fMillerOrientation[1]=k;
    fMillerOrientation[2]=n;
    // // // // // // // // // // // // //
}

void XPhysicalLattice::SetPhysicalVolume

(

G4VPhysicalVolume *

Vol

)

Definition at line 238 of file XPhysicalLattice.cc.

                                                              {
    fVolume=Vol;
}

void XPhysicalLattice::SetScatteringConstant

(

G4double

a

)

Definition at line 95 of file XPhysicalLattice.cc.

                                                      {
    fA=a;
}

void XPhysicalLattice::SetSTDOS

(

double

STDOS

)

Definition at line 115 of file XPhysicalLattice.cc.

{
    fDosST = STDOS;
}

void XPhysicalLattice::SetThermalVibrationAmplitude

(

G4double

vThermalVibrationAmplitude

)

Definition at line 288 of file XPhysicalLattice.cc.

                                                                 {
    fThermalVibrationAmplitude = vThermalVibrationAmplitude;
}

void XPhysicalLattice::SetUnitCell

(

XUnitCell *

cell

)

Definition at line 344 of file XPhysicalLattice.cc.

                                                 {
    fUnitCell = cell;
}

void XPhysicalLattice::SetXLogicalLattice

(

XLogicalLattice *

Lat

)

Definition at line 276 of file XPhysicalLattice.cc.

                                                             {
    fLattice=Lat;
}

Member Data Documentation

G4AffineTransform XPhysicalLattice::fGlobalToLocal

Definition at line 57 of file XPhysicalLattice.hh.

G4AffineTransform XPhysicalLattice::fLocalToGlobal

Definition at line 56 of file XPhysicalLattice.hh.

---

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

• source/geant4.10.03.p03/examples/extended/exoticphysics/channeling/include/XPhysicalLattice.hh
• source/geant4.10.03.p03/examples/extended/exoticphysics/channeling/src/XPhysicalLattice.cc