XVCrystalPlanarAnalytical Class Reference

#include <XVCrystalPlanarAnalytical.hh>

Inheritance diagram for XVCrystalPlanarAnalytical:

Collaboration diagram for XVCrystalPlanarAnalytical:

Public Member Functions
void	SetNumberOfPlanes (G4int)
G4int	GetNumberOfPlanes ()
virtual G4double	ComputeECForSinglePlane (G4double, XPhysicalLattice *)=0
G4ThreeVector	ComputeEC (G4ThreeVector, XPhysicalLattice *)
G4ThreeVector	ComputeECFromVector (G4ThreeVector)
G4ThreeVector	ComputePositionInUnitCell (G4ThreeVector, XPhysicalLattice *)
virtual G4double	ComputeMaximum (XPhysicalLattice *)
virtual G4double	ComputeMinimum (XPhysicalLattice *)
virtual void	PrintOnFile (const G4String &, XPhysicalLattice *, G4double)
virtual void	ReadFromFile (const G4String &, XPhysicalLattice *, G4double=1)
virtual void	ReadFromECHARM (const G4String &, G4double=1)
void	InitializeVector ()
	XVCrystalPlanarAnalytical ()
	~XVCrystalPlanarAnalytical ()
Public Member Functions inherited from XVCrystalCharacteristic
XPhysicalLattice *	GetXPhysicalLattice (G4VPhysicalVolume *)
XUnitCell *	GetXUnitCell (G4VPhysicalVolume *)
XLogicalLattice *	GetLogicalLattice (G4VPhysicalVolume *)
void	InitializePhysicalLattice (XPhysicalLattice *)
G4ThreeVector	GetEC (G4ThreeVector, XPhysicalLattice *)
virtual G4double	ComputeTFScreeningRadius (XPhysicalLattice *)
virtual G4double	GetMaximum (XPhysicalLattice *)
virtual G4double	GetMinimum (XPhysicalLattice *)
G4bool	IsInitialized (XPhysicalLattice *)
	XVCrystalCharacteristic ()
	~XVCrystalCharacteristic ()

Additional Inherited Members
Protected Attributes inherited from XVCrystalCharacteristic
G4double	fMaximum
G4double	fMinimum
XPhysicalLattice *	fPhysicalLattice
G4PhysicsVector *	fVectorEC

Detailed Description

Definition at line 35 of file XVCrystalPlanarAnalytical.hh.

Constructor & Destructor Documentation

XVCrystalPlanarAnalytical::XVCrystalPlanarAnalytical

(

)

Definition at line 33 of file XVCrystalPlanarAnalytical.cc.

                                                    {
    fNumberOfPlanes = 4;
}

XVCrystalPlanarAnalytical::~XVCrystalPlanarAnalytical

(

)

Definition at line 39 of file XVCrystalPlanarAnalytical.cc.

                                                     {
}

Member Function Documentation

G4ThreeVector XVCrystalPlanarAnalytical::ComputeEC ( G4ThreeVector  vPositionVector, XPhysicalLattice *  vLattice  )

virtual

Implements XVCrystalCharacteristic.

Definition at line 57 of file XVCrystalPlanarAnalytical.cc.

                                     {
    G4VPhysicalVolume* vVolume = vLattice->GetVolume();
    
    G4double vInterplanarDistance =
        GetXPhysicalLattice(vVolume)->ComputeInterplanarPeriod();
    
    G4double vPosition = ComputePositionInUnitCell(vPositionVector,
                                                   vLattice).x();
    
    G4double vValue = 0.;
    for(int i=-int(GetNumberOfPlanes()/2);i<=+int(GetNumberOfPlanes()/2);i++){
        vValue += ComputeECForSinglePlane( ( vPosition + G4double(i) ) *
                                          vInterplanarDistance ,
                                          vLattice );
    }
    
    return G4ThreeVector(vValue,0.,0.);
}

Here is the call graph for this function:

Here is the caller graph for this function:

virtual G4double XVCrystalPlanarAnalytical::ComputeECForSinglePlane ( G4double  , XPhysicalLattice *    )

pure virtual

Implemented in XCrystalPlanarMoliereElectronDensity, XCrystalPlanarMoliereTempPotential, XCrystalPlanarMoliereElectricField, XCrystalPlanarMolierePotential, and XCrystalPlanarNucleiDensity.

Here is the caller graph for this function:

G4ThreeVector XVCrystalPlanarAnalytical::ComputeECFromVector ( G4ThreeVector  vPosition )

virtual

Implements XVCrystalCharacteristic.

Definition at line 80 of file XVCrystalPlanarAnalytical.cc.

                                            {
    G4double vInterplanarPeriod = fPhysicalLattice->ComputeInterplanarPeriod();
    G4double vX = vPosition.x();
        if (vX < 0.0) {
            vX += ((int( - vX / vInterplanarPeriod) + 1.0 )
                  * vInterplanarPeriod);
        }
        else if( vX > vInterplanarPeriod ){
            vX -= ( int( vX / vInterplanarPeriod) * vInterplanarPeriod );
        }
        return G4ThreeVector(fVectorEC->Value(vX),0.,0.);
}

Here is the call graph for this function:

G4double XVCrystalPlanarAnalytical::ComputeMaximum ( XPhysicalLattice *  vLattice )

virtual

Reimplemented from XVCrystalCharacteristic.

Reimplemented in XCrystalPlanarMoliereTempPotential, and XCrystalPlanarMolierePotential.

Definition at line 121 of file XVCrystalPlanarAnalytical.cc.

                                                                            {
    unsigned int vPrecision = 1024;
    G4VPhysicalVolume* vVolume = vLattice->GetVolume();
    G4double vStep = GetXPhysicalLattice(vVolume)->ComputeInterplanarPeriod()
        / vPrecision;
    
    G4double vMaximum = -DBL_MAX;
    G4double vValue;
    
    for(unsigned int i=0;i<vPrecision;i++){
        if( (vValue = GetEC(G4ThreeVector(vStep * i,0.,0.),vLattice).x() )
           > vMaximum) {vMaximum = vValue;}
    }
    return vMaximum;
}

Here is the call graph for this function:

G4double XVCrystalPlanarAnalytical::ComputeMinimum ( XPhysicalLattice *  vLattice )

virtual

Reimplemented from XVCrystalCharacteristic.

Reimplemented in XCrystalPlanarMoliereTempPotential, and XCrystalPlanarMolierePotential.

Definition at line 139 of file XVCrystalPlanarAnalytical.cc.

                                                                            {
    unsigned int vPrecision = 1024;
    G4VPhysicalVolume* vVolume = vLattice->GetVolume();
    G4double vStep = GetXPhysicalLattice(vVolume)->ComputeInterplanarPeriod()
        / vPrecision;
    
    G4double vMinimum = +DBL_MAX;
    G4double vValue;
    
    for(unsigned int i=0;i<vPrecision;i++){
        if( (vValue = GetEC(G4ThreeVector(vStep * i,0.,0.),vLattice).x() )
           < vMinimum) {vMinimum = vValue;}
    }
    return vMinimum;
}

Here is the call graph for this function:

G4ThreeVector XVCrystalPlanarAnalytical::ComputePositionInUnitCell ( G4ThreeVector  vPosition, XPhysicalLattice *  vLattice  )

virtual

Reimplemented from XVCrystalCharacteristic.

Definition at line 96 of file XVCrystalPlanarAnalytical.cc.

                                                     {
    G4VPhysicalVolume* vVolume = vLattice->GetVolume();

    G4double vInterplanarPeriod =
        GetXPhysicalLattice(vVolume)->ComputeInterplanarPeriod();
    
    G4double vPositionX = vPosition.x();
    
    if((vPositionX >= 0.) &&
       (vPositionX < vInterplanarPeriod)){
        return G4ThreeVector(vPositionX/vInterplanarPeriod,0.,0.);
    }
    else if(vPositionX == vInterplanarPeriod){
        return G4ThreeVector(0.,0.,0.);
    }
    else{
        vPositionX -= std::fmod(vPosition.x(),vInterplanarPeriod)
            * vInterplanarPeriod;
        return G4ThreeVector(vPositionX/vInterplanarPeriod,0.,0.);
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4int XVCrystalPlanarAnalytical::GetNumberOfPlanes

(

)

Definition at line 50 of file XVCrystalPlanarAnalytical.cc.

                                                  {
    return fNumberOfPlanes;
}

Here is the caller graph for this function:

void XVCrystalPlanarAnalytical::InitializeVector ( )

virtual

Implements XVCrystalCharacteristic.

Definition at line 251 of file XVCrystalPlanarAnalytical.cc.

                                                {
    G4int imax = 4096;
    G4double vXposition = 0.;
    G4double vInterplanarPeriod = fPhysicalLattice->ComputeInterplanarPeriod();
    
    fVectorEC = new G4PhysicsLinearVector(0,
                                          vInterplanarPeriod*(imax-1)/imax,
                                          imax);
    for(G4int i = 0;i<imax;i++){
        vXposition = double(i) / double(imax) * vInterplanarPeriod;
        fVectorEC->PutValue(i,ComputeEC(G4ThreeVector(vXposition,0.,0.),
                                        fPhysicalLattice).x());
    }
}

Here is the call graph for this function:

void XVCrystalPlanarAnalytical::PrintOnFile ( const G4String &  filename, XPhysicalLattice *  vLattice, G4double  vUnit  )

virtual

Implements XVCrystalCharacteristic.

Definition at line 157 of file XVCrystalPlanarAnalytical.cc.

                                                           {
    std::ofstream vFileOut;
    vFileOut.open(filename);
    vFileOut << "pos val" << std::endl;
    
    G4int imax = 8192;
    G4double vXposition = 0.;
    G4double vInterplanarPeriod = vLattice->ComputeInterplanarPeriod();
    
    for(G4int i = 0;i<imax;i++){
        vXposition = double(i) / double(imax) * vInterplanarPeriod;
        vFileOut << vXposition / CLHEP::angstrom << " ";
        vFileOut << GetEC(G4ThreeVector(vXposition,0.,0.),vLattice).x() / vUnit;
        vFileOut << std::endl;
    }
    
    vFileOut.close();
}

Here is the call graph for this function:

void XVCrystalPlanarAnalytical::ReadFromECHARM ( const G4String &  filename, G4double  vConversion = 1  )

virtual

Implements XVCrystalCharacteristic.

Definition at line 214 of file XVCrystalPlanarAnalytical.cc.

                                                                    {
    std::ifstream vFileIn;
    vFileIn.open(filename);
    
    G4int imax;
    G4double xmax;
    
    vFileIn >> imax;
    vFileIn >> xmax;

    xmax *= CLHEP::meter;
    fMaximum = -DBL_MAX;
    fMinimum = +DBL_MAX;

    fVectorEC = new G4PhysicsLinearVector(0,xmax,imax);
    
    for(G4int i=0;i<imax; i++){ 
        double vTempX;
        vFileIn >> vTempX;

        vTempX *= vConversion;
        if(vTempX > fMaximum) {fMaximum = vTempX;}
        if(vTempX < fMinimum) {fMinimum = vTempX;}
        fVectorEC->PutValue(i,vTempX);
    }

    G4cout << "XVCrystalPlanarAnalytical::ReadFromECHARM() - " <<
     vConversion << " " << imax << " " << xmax << " " << 
     fMinimum << " " << fMaximum << G4endl;
   
    vFileIn.close();
    
}

Here is the call graph for this function:

void XVCrystalPlanarAnalytical::ReadFromFile ( const G4String &  filename, XPhysicalLattice *  , G4double  = 1  )

virtual

Implements XVCrystalCharacteristic.

Definition at line 180 of file XVCrystalPlanarAnalytical.cc.

                                                                        {
    std::ifstream vFileIn;
    vFileIn.open(filename);
    std::string vTempString;
    vFileIn >> vTempString >> vTempString;
 
    double vTempX;
    double vTempVal;
    std::vector<double> vTempXvector;
    std::vector<double> vTempValvector;
    
    while(!vFileIn.eof()){ 
        vFileIn >> vTempX;
        vFileIn >> vTempVal;
        vTempXvector.push_back(vTempX);
        vTempValvector.push_back(vTempVal);
    }
    
    vFileIn.close();

    G4int imax = vTempXvector.size();
    G4double vInterplanarPeriod = vTempXvector.at(imax-1);

    fVectorEC = new G4PhysicsLinearVector(0,
                                          vInterplanarPeriod*(imax-1)/imax,
                                          imax);
    for(G4int i = 0;i<imax;i++){
        fVectorEC->PutValue(i,vTempValvector.at(i));
    }
}

Here is the call graph for this function:

void XVCrystalPlanarAnalytical::SetNumberOfPlanes

(

G4int

vNumberOfPlanes

)

Definition at line 44 of file XVCrystalPlanarAnalytical.cc.

                                                                      {
    fNumberOfPlanes = vNumberOfPlanes;
}

---

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

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