XCrystalPlanarMoliereTempPotential Class Reference

#include <XCrystalPlanarMoliereTempPotential.hh>

Inheritance diagram for XCrystalPlanarMoliereTempPotential:

Collaboration diagram for XCrystalPlanarMoliereTempPotential:

Public Member Functions
G4double	ComputeECForSinglePlane (G4double, XPhysicalLattice *)
G4double	ComputeMaximum (XPhysicalLattice *)
G4double	ComputeMinimum (XPhysicalLattice *)
	XCrystalPlanarMoliereTempPotential ()
	~XCrystalPlanarMoliereTempPotential ()
Public Member Functions inherited from XVCrystalPlanarAnalytical
void	SetNumberOfPlanes (G4int)
G4int	GetNumberOfPlanes ()
G4ThreeVector	ComputeEC (G4ThreeVector, XPhysicalLattice *)
G4ThreeVector	ComputeECFromVector (G4ThreeVector)
G4ThreeVector	ComputePositionInUnitCell (G4ThreeVector, 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 36 of file XCrystalPlanarMoliereTempPotential.hh.

Constructor & Destructor Documentation

XCrystalPlanarMoliereTempPotential::XCrystalPlanarMoliereTempPotential

(

)

Definition at line 34 of file XCrystalPlanarMoliereTempPotential.cc.

                                                                      {
    fAlfa[0] = 0.1;
    fAlfa[1] = 0.55;
    fAlfa[2] = 0.35;
    
    fBeta[0] = 6.0;
    fBeta[1] = 1.2;
    fBeta[2] = 0.3;
    
    for(unsigned int i=0;i<3;i++) {
        fGamma[i] = fAlfa[i]/fBeta[i];
    }
}

XCrystalPlanarMoliereTempPotential::~XCrystalPlanarMoliereTempPotential

(

)

Definition at line 50 of file XCrystalPlanarMoliereTempPotential.cc.

                                                                       {
}

Member Function Documentation

G4double XCrystalPlanarMoliereTempPotential::ComputeECForSinglePlane ( G4double  vX, XPhysicalLattice *  vLattice  )

virtual

Implements XVCrystalPlanarAnalytical.

Definition at line 56 of file XCrystalPlanarMoliereTempPotential.cc.

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

    G4double aTF = ComputeTFScreeningRadius(vLattice);
    G4double vTVA = vLattice->GetThermalVibrationAmplitude();
    
    G4double vTau[3];
    for(unsigned int i=0;i<3;i++){
        vTau[i] = (std::pow( vTVA / aTF * fBeta[i] , 2. ) / 2.0);
    }

    
    G4double vValueForSinglePlane = 0.;
    
    for(unsigned int i=0;i<3;i++){
        G4double vTemp = 0.;
        vTemp += ( std::exp(-vX/ aTF * fBeta[i] ) *
                 (1.-CLHEP::HepStat::erfQ((vTVA / aTF * fBeta[i] - vX/ vTVA)
                     / std::pow(2.,0.5))) );
        vTemp += ( std::exp( vX/ aTF * fBeta[i] ) *
                 (1.-CLHEP::HepStat::erfQ((vTVA / aTF * fBeta[i] + vX/ vTVA)
                     / std::pow(2.,0.5))) );
        vValueForSinglePlane += ( vTemp * fGamma[i] * std::exp( vTau[i] ) /2.0);
    }
    
    vValueForSinglePlane *= 2. * CLHEP::pi;
    vValueForSinglePlane *= GetXUnitCell(vVolume)->ComputeDirectPeriod(
                                GetXPhysicalLattice(vVolume)->GetMiller(0),
                                GetXPhysicalLattice(vVolume)->GetMiller(1),
                                GetXPhysicalLattice(vVolume)->GetMiller(2));
    
    vValueForSinglePlane *= aTF;

    vValueForSinglePlane *= (CLHEP::elm_coupling);

    vValueForSinglePlane *= (GetXUnitCell(vVolume)->ComputeAtomVolumeDensity());

    return vValueForSinglePlane;
}

Here is the call graph for this function:

G4double XCrystalPlanarMoliereTempPotential::ComputeMaximum ( XPhysicalLattice *  vLattice )

virtual

Reimplemented from XVCrystalPlanarAnalytical.

Definition at line 101 of file XCrystalPlanarMoliereTempPotential.cc.

                                          {

    G4double vMaximum = GetEC(G4ThreeVector(0.,0.,0.),
                              vLattice).x();
    
    return vMaximum;
}

Here is the call graph for this function:

G4double XCrystalPlanarMoliereTempPotential::ComputeMinimum ( XPhysicalLattice *  vLattice )

virtual

Reimplemented from XVCrystalPlanarAnalytical.

Definition at line 112 of file XCrystalPlanarMoliereTempPotential.cc.

                                          {
    G4VPhysicalVolume* vVolume = vLattice->GetVolume();
    G4double vInterplanarDistance =
        GetXUnitCell(vVolume)->ComputeDirectPeriod(
                            GetXPhysicalLattice(vVolume)->GetMiller(0),
                            GetXPhysicalLattice(vVolume)->GetMiller(1),
                            GetXPhysicalLattice(vVolume)->GetMiller(2));
    
    G4double vMinimum = GetEC(G4ThreeVector(vInterplanarDistance/2.,0.,0.),
                              vLattice).x();
    
    return vMinimum;
}

Here is the call graph for this function:

---

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

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