File indexing completed on 2025-01-31 09:22:05
0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023
0024
0025
0026
0027
0028
0029
0030
0031
0032
0033
0034
0035
0036
0037
0038
0039
0040
0041
0042
0043
0044
0045
0046 #ifndef G4RDPhotoElectricAngularGeneratorPolarized_h
0047 #define G4RDPhotoElectricAngularGeneratorPolarized_h 1
0048
0049 #include "G4RDVPhotoElectricAngularDistribution.hh"
0050 #include "G4ios.hh"
0051 #include "globals.hh"
0052 #include "G4RotationMatrix.hh"
0053
0054 class G4RDPhotoElectricAngularGeneratorPolarized : public G4RDVPhotoElectricAngularDistribution
0055 {
0056
0057 public:
0058
0059 G4RDPhotoElectricAngularGeneratorPolarized(const G4String& name);
0060
0061 ~G4RDPhotoElectricAngularGeneratorPolarized();
0062
0063 G4ThreeVector GetPhotoElectronDirection(const G4ThreeVector& direction,
0064 const G4double kineticEnergy,
0065 const G4ThreeVector& polarization,
0066 const G4int shellId) const;
0067
0068 void PrintGeneratorInformation() const;
0069
0070 protected:
0071 G4ThreeVector SetPerpendicularVector(const G4ThreeVector& a) const;
0072
0073 private:
0074
0075
0076 G4RDPhotoElectricAngularGeneratorPolarized & operator=(const G4RDPhotoElectricAngularGeneratorPolarized &right);
0077 G4RDPhotoElectricAngularGeneratorPolarized(const G4RDPhotoElectricAngularGeneratorPolarized&);
0078
0079 void PhotoElectronGetMajorantSurfaceAandCParameters(const G4int shellLevel,
0080 const G4double beta,
0081 G4double *majorantSurfaceParameterA,
0082 G4double *majorantSurfaceParameterC) const;
0083
0084 void PhotoElectronGeneratePhiAndTheta(const G4int shellLevel,
0085 const G4double beta,
0086 const G4double aBeta,
0087 const G4double cBeta,
0088 G4double *pphi,
0089 G4double *ptheta) const;
0090
0091 G4ThreeVector PhotoElectronComputeFinalDirection(const G4RotationMatrix& rotation,
0092 const G4double theta,
0093 const G4double phi) const;
0094
0095 G4RotationMatrix PhotoElectronRotationMatrix(const G4ThreeVector& direction,
0096 const G4ThreeVector& polarization) const;
0097
0098 G4double GetMax(const G4double arg1, const G4double arg2) const;
0099
0100 G4double CrossSectionMajorantFunction(const G4double theta, const G4double cBeta) const;
0101 G4double DSigmaKshellGavrila1959(const G4double beta, const G4double theta, const G4double phi) const;
0102 G4double DSigmaL1shellGavrila(const G4double beta, const G4double theta, const G4double phi) const;
0103
0104 G4double betaArray[3];
0105 G4double aMajorantSurfaceParameterTable[980][2], cMajorantSurfaceParameterTable[980][2];
0106 };
0107
0108 #endif
0109