29 #include "CLHEP/Random/Stat.h"
40 for(
unsigned int i=0;i<3;i++) {
62 for(
unsigned int i=0;i<3;i++){
63 vTau[i] = (pow( vTVA / aTF *
fBeta[i] , 2. ) / 2.0);
69 for(
unsigned int i=0;i<3;i++){
71 vTemp += ( exp(-vX/ aTF *
fBeta[i] ) *
72 (1.-CLHEP::HepStat::erfQ((vTVA / aTF *
fBeta[i] - vX/ vTVA) / pow(2.,0.5))) );
73 vTemp += ( exp( vX/ aTF *
fBeta[i] ) *
74 (1.-CLHEP::HepStat::erfQ((vTVA / aTF *
fBeta[i] + vX/ vTVA) / pow(2.,0.5))) );
75 vValueForSinglePlane += ( vTemp *
fGamma[i] * exp( vTau[i] ) /2.0);
84 vValueForSinglePlane *= aTF;
86 vValueForSinglePlane *= (CLHEP::elm_coupling);
90 return vValueForSinglePlane;
CLHEP::Hep3Vector G4ThreeVector
XPhysicalLattice * GetXPhysicalLattice(G4VPhysicalVolume *)
XCrystalPlanarMoliereTempPotential()
G4double ComputeMaximum(XPhysicalLattice *)
virtual G4double ComputeTFScreeningRadius(XPhysicalLattice *)
G4double ComputeAtomVolumeDensity()
G4double GetThermalVibrationAmplitude()
G4ThreeVector GetEC(G4ThreeVector, XPhysicalLattice *)
G4double ComputeECForSinglePlane(G4double, XPhysicalLattice *)
G4double ComputeMinimum(XPhysicalLattice *)
~XCrystalPlanarMoliereTempPotential()
G4double ComputeDirectPeriod(G4int, G4int, G4int)
XUnitCell * GetXUnitCell(G4VPhysicalVolume *)
G4VPhysicalVolume * GetVolume()