XVCrystalIntegratedDensity Class Reference

#include <XVCrystalIntegratedDensity.hh>

Inheritance diagram for XVCrystalIntegratedDensity:

Public Member Functions
void	SetIntegrationPoints (unsigned int, unsigned int)
unsigned int	GetIntegrationPoints (unsigned int)
unsigned int	GetIntegrationPoints ()
void	SetNumberOfPoints (unsigned int)
unsigned int	GetNumberOfPoints ()
void	SetDensity (XVCrystalCharacteristic *)
XVCrystalCharacteristic *	GetDensity ()
void	SetPotential (XVCrystalCharacteristic *)
XVCrystalCharacteristic *	GetPotential ()
void	SetXPhysicalLattice (XPhysicalLattice *)
XPhysicalLattice *	GetXPhysicalLattice ()
void	SetParticleCharge (G4int)
G4int	GetParticleCharge ()
void	PrintOnFile (const G4String &)
void	ReadFromFile (const G4String &)
G4bool	HasBeenInitialized (XPhysicalLattice *, G4int)
G4double	GetIntegratedDensity (G4double, XPhysicalLattice *, G4int)
virtual void	InitializeTable ()
	XVCrystalIntegratedDensity ()
	~XVCrystalIntegratedDensity ()

Protected Member Functions
G4double	GetStep ()
virtual void	ComputePotentialParameters ()
virtual G4double	ComputeIntegratedDensity (G4double, G4int)
G4double	FindCatmullRomInterpolate (G4double &p0, G4double &p1, G4double &p2, G4double &p3, G4double &x)

Protected Attributes
G4double	fPotentialMinimum
G4double	fPotentialMaximum
G4double	fPotentialRange

Detailed Description

Definition at line 36 of file XVCrystalIntegratedDensity.hh.

Constructor & Destructor Documentation

XVCrystalIntegratedDensity::XVCrystalIntegratedDensity

(

)

Definition at line 32 of file XVCrystalIntegratedDensity.cc.

                                                      {
    fNumberOfPoints = 512;
    fIntegrationPoints[0] = 32;
    fIntegrationPoints[1] = 32;
    fIntegrationPoints[2] = 32;
}

XVCrystalIntegratedDensity::~XVCrystalIntegratedDensity

(

)

Definition at line 41 of file XVCrystalIntegratedDensity.cc.

                                                       {
}

Member Function Documentation

G4double XVCrystalIntegratedDensity::ComputeIntegratedDensity ( G4double  vPotentialInitial, G4int    )

protectedvirtual

Reimplemented in XCrystalIntegratedDensityPlanar.

Definition at line 231 of file XVCrystalIntegratedDensity.cc.

                               {
    
    unsigned int i1,i2,i3;
    i1 = i2 = i3 = 0;
    
    G4ThreeVector vPositionTemp = G4ThreeVector(0.,0.,0.);
    G4double vDensity = 0.;
    
    G4ThreeVector vSize = fLattice->GetXUnitCell()->GetSize();
    while(i1<fIntegrationPoints[2]){
        vPositionTemp.setY(G4double(G4double(i3)/
                            G4double(fIntegrationPoints[2])*vSize.z()));
        while(i1<fIntegrationPoints[1]){
            vPositionTemp.setZ(G4double(G4double(i2)/
                                G4double(fIntegrationPoints[1])*vSize.y()));
            while(i1<fIntegrationPoints[0]){
                vPositionTemp.setX(G4double(G4double(i1)/
                                G4double(fIntegrationPoints[0])*vSize.x()));
                if(fPotential->GetEC(vPositionTemp,fLattice).x()
                   < vPotentialInitial){
                    vDensity += fDensity->GetEC(vPositionTemp,fLattice).x();
                }
                i1++;
            };
            i2++;
        };
        i3++;
    };
    
    vDensity *= fLattice->GetXUnitCell()->ComputeVolume();
    vDensity /= GetIntegrationPoints();
    
    return vDensity;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void XVCrystalIntegratedDensity::ComputePotentialParameters ( )

protectedvirtual

Definition at line 160 of file XVCrystalIntegratedDensity.cc.

                                                           {
    fPotentialMinimum = fPotential->GetMinimum(fLattice);
    if(fParticleCharge < 0.){
        fPotentialMinimum = - fPotential->GetMaximum(fLattice);
    }
    
    fPotentialMaximum = fPotential->GetMaximum(fLattice);
    if(fParticleCharge < 0.){
        fPotentialMaximum = - fPotential->GetMinimum(fLattice);
    }
    
    fPotentialRange = std::fabs(fPotentialMaximum - fPotentialMinimum);
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4double XVCrystalIntegratedDensity::FindCatmullRomInterpolate ( G4double &  p0, G4double &  p1, G4double &  p2, G4double &  p3, G4double &  x  )

protected

Definition at line 270 of file XVCrystalIntegratedDensity.cc.

{
    double a0, a1 , a2 , a3 , x2;
    
    x2 = x * x;
    a0 = -0.5 * p0 + 1.5 * p1 - 1.5 * p2 + 0.5 * p3;
    a1 = p0 - 2.5 * p1 + 2.0 * p2 - 0.5 * p3;
    a2 = -0.5 * p0 + 0.5 * p2;
    a3 = p1;
    
    return (a0 * x * x2 + a1 * x2 + a2 * x + a3);
}

XVCrystalCharacteristic * XVCrystalIntegratedDensity::GetDensity

(

)

Definition at line 109 of file XVCrystalIntegratedDensity.cc.

                                                               {
    return fDensity;
}

Here is the caller graph for this function:

G4double XVCrystalIntegratedDensity::GetIntegratedDensity	(	G4double	vPotential,
		XPhysicalLattice *	vLattice,
		G4int	vParticleCharge
	)

Definition at line 205 of file XVCrystalIntegratedDensity.cc.

                                           {
        
    G4double vPotentialModified = vPotential /
        std::fabs(G4double(vParticleCharge)) - fPotentialMinimum;
    
    // if the crystal has not been initialized return -1 -> Error!
    if(!HasBeenInitialized(vLattice,vParticleCharge))
        return -1.;
    // if the potential is higher than the maximum the average density
    // is equal to the one of the amorphous material
    else if(vPotentialModified >= std::fabs(fPotentialMaximum - fPotentialMinimum))
        return 1.;
    // if the value is less than zero (because of possible variation
    // due to centrifugal force) take the zero value
    else if(vPotentialModified < 0.) return fTableVector->Value(0.);
    //else if(vPotentialModified < 1. && vParticleCharge < 0.) return 1.;
    else{
        return fTableVector->Value(vPotentialModified);
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

unsigned int XVCrystalIntegratedDensity::GetIntegrationPoints

(

unsigned int

vIndex

)

Definition at line 59 of file XVCrystalIntegratedDensity.cc.

                                         {
    if(vIndex<3) {
        return fIntegrationPoints[vIndex];
    }
    else{
        return 0;
    }
}

unsigned int XVCrystalIntegratedDensity::GetIntegrationPoints

(

)

Definition at line 70 of file XVCrystalIntegratedDensity.cc.

                                                             {
    return fIntegrationPoints[0]*fIntegrationPoints[1]*fIntegrationPoints[2];
}

Here is the caller graph for this function:

unsigned int XVCrystalIntegratedDensity::GetNumberOfPoints

(

)

Definition at line 83 of file XVCrystalIntegratedDensity.cc.

                                                          {
    return fNumberOfPoints;
}

Here is the caller graph for this function:

G4int XVCrystalIntegratedDensity::GetParticleCharge

(

)

Definition at line 134 of file XVCrystalIntegratedDensity.cc.

                                                    {
    return fParticleCharge;
}

XVCrystalCharacteristic * XVCrystalIntegratedDensity::GetPotential

(

)

Definition at line 96 of file XVCrystalIntegratedDensity.cc.

                                                                 {
    return fPotential;
}

Here is the caller graph for this function:

G4double XVCrystalIntegratedDensity::GetStep ( )

protected

Definition at line 140 of file XVCrystalIntegratedDensity.cc.

                                            {
    return fPotentialRange / fNumberOfPoints;
}

Here is the caller graph for this function:

XPhysicalLattice * XVCrystalIntegratedDensity::GetXPhysicalLattice

(

)

Definition at line 122 of file XVCrystalIntegratedDensity.cc.

                                                                  {
    return fLattice;
}

Here is the caller graph for this function:

G4bool XVCrystalIntegratedDensity::HasBeenInitialized	(	XPhysicalLattice *	vLattice,
		G4int	vParticleCharge
	)

Definition at line 147 of file XVCrystalIntegratedDensity.cc.

                                         {
    //now it checks only of the table is initialized,
    //it does not check if the particular crystal is initialized.
    //To be changed in the future!
    if(fTableVector->GetVectorLength() == 0) return false;
    else if(vLattice!=fLattice) return false;
    else if(vParticleCharge!=fParticleCharge) return false;
    else return true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void XVCrystalIntegratedDensity::InitializeTable ( )

virtual

Definition at line 176 of file XVCrystalIntegratedDensity.cc.

                                                {
    
    ComputePotentialParameters();
    
    G4cout << "XVCrystalIntegratedDensity::InitializeTable()::";
    G4cout << "Potential Range =  " << fPotentialRange/CLHEP::eV;
    G4cout << " - Minimum = " << fPotentialMinimum / CLHEP::eV;
    G4cout << " - Maximum " << fPotentialMaximum / CLHEP::eV << G4endl;
    
    G4double vPotentialInitial = 0.;
    
    fTableVector =
    new G4PhysicsLinearVector(0.,
                               std::fabs(fPotentialMaximum - fPotentialMinimum),
                               GetNumberOfPoints());
    
    G4double vValue = 0;
    for(unsigned int i=0;i<GetNumberOfPoints();i++){
        vPotentialInitial = (fPotentialMinimum +
                             fPotentialRange * G4double(i+1)
                             / G4double(GetNumberOfPoints()));
        vValue = ComputeIntegratedDensity(vPotentialInitial,fParticleCharge);
        fTableVector->PutValue(i,vValue);
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

void XVCrystalIntegratedDensity::PrintOnFile

(

const G4String &

filename

)

Definition at line 289 of file XVCrystalIntegratedDensity.cc.

                                                                    {
    std::ofstream vFileOut;
    vFileOut.open(filename);
    G4double vStep = GetStep();
    
    vFileOut << "energy,dens" << std::endl;
    for(unsigned int i = 0;i < fNumberOfPoints;i++){
        vFileOut << i * vStep  / CLHEP::eV;
        vFileOut << " , ";
        vFileOut << GetIntegratedDensity((i+1) * vStep
                    + fPotentialMinimum,fLattice,fParticleCharge);
        vFileOut << std::endl;
    }
    vFileOut.close();
}

Here is the call graph for this function:

Here is the caller graph for this function:

void XVCrystalIntegratedDensity::ReadFromFile

(

const G4String &

filename

)

Definition at line 307 of file XVCrystalIntegratedDensity.cc.

                                                                     {
    std::ifstream vFileIn;
    vFileIn.open(filename);
    if(!vFileIn){
        G4cout << "XVCrystalIntegratedDensity::";
        G4cout << "ReadFromFile - ERROR READING FILE!!!!! ";
        G4cout << filename << G4endl;
    }
    
    fPotentialMinimum = +DBL_MAX;
    fPotentialMaximum = -DBL_MAX;
    G4double vDensity = 0.;
    
    vFileIn >> fPotentialMinimum;
    vFileIn >> fPotentialMaximum;
    
    fPotentialMinimum *= CLHEP::eV;
    fPotentialMaximum *= CLHEP::eV;

    std::vector<G4double> fTable;
    
    while(!vFileIn.eof()){
        vFileIn >> vDensity;
        if(vDensity < 1.E-2) {
            vDensity = 1.E-2;
        }
        fTable.push_back(vDensity);
    };
    
    fNumberOfPoints = fTable.size();
    
    fTableVector =
        new G4PhysicsLinearVector(0.,
                                  std::fabs(fPotentialMaximum - fPotentialMinimum),
                                  fNumberOfPoints);
    
    for(unsigned int i=0;i<fTable.size();i++){
        fTableVector->PutValue(i,fTable.at(i));
    }
    
    fPotentialRange = std::fabs(fPotentialMaximum - fPotentialMinimum);
    
    G4cout << "XVCrystalIntegratedDensity::InitializeTable()::";
    G4cout << "Potential Range =  " << fPotentialRange/CLHEP::eV;
    G4cout << " - Minimum = " << fPotentialMinimum / CLHEP::eV;
    G4cout << " - Maximum " << fPotentialMaximum / CLHEP::eV << G4endl;
    vFileIn.close();
}

Here is the call graph for this function:

Here is the caller graph for this function:

void XVCrystalIntegratedDensity::SetDensity

(

XVCrystalCharacteristic *

vDensity

)

Definition at line 103 of file XVCrystalIntegratedDensity.cc.

                                             {
    fDensity = vDensity;
}

Here is the caller graph for this function:

void XVCrystalIntegratedDensity::SetIntegrationPoints	(	unsigned int	vIndex,
		unsigned int	vIntegrationPoints
	)

Definition at line 47 of file XVCrystalIntegratedDensity.cc.

                                                     {
    if(vIndex<3) {
        if(vIntegrationPoints > 0){
            fIntegrationPoints[vIndex] = vIntegrationPoints;
        }
    }
}

Here is the caller graph for this function:

void XVCrystalIntegratedDensity::SetNumberOfPoints

(

unsigned int

vNumberOfPoints

)

Definition at line 77 of file XVCrystalIntegratedDensity.cc.

                                               {
    fNumberOfPoints = vNumberOfPoints;
}

Here is the caller graph for this function:

void XVCrystalIntegratedDensity::SetParticleCharge

(

G4int

vParticleCharge

)

Definition at line 128 of file XVCrystalIntegratedDensity.cc.

                                                                       {
    fParticleCharge = vParticleCharge;
}

Here is the caller graph for this function:

void XVCrystalIntegratedDensity::SetPotential

(

XVCrystalCharacteristic *

vPotential

)

Definition at line 90 of file XVCrystalIntegratedDensity.cc.

                                                 {
    fPotential = vPotential;
}

Here is the caller graph for this function:

void XVCrystalIntegratedDensity::SetXPhysicalLattice

(

XPhysicalLattice *

vLattice

)

Definition at line 116 of file XVCrystalIntegratedDensity.cc.

                                               {
    fLattice = vLattice;
}

Here is the caller graph for this function:

Member Data Documentation

G4double XVCrystalIntegratedDensity::fPotentialMaximum

protected

Definition at line 89 of file XVCrystalIntegratedDensity.hh.

G4double XVCrystalIntegratedDensity::fPotentialMinimum

protected

Definition at line 88 of file XVCrystalIntegratedDensity.hh.

G4double XVCrystalIntegratedDensity::fPotentialRange

protected

Definition at line 90 of file XVCrystalIntegratedDensity.hh.

---

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

• source/geant4.10.03.p02/examples/extended/exoticphysics/channeling/include/XVCrystalIntegratedDensity.hh
• source/geant4.10.03.p02/examples/extended/exoticphysics/channeling/src/XVCrystalIntegratedDensity.cc