File indexing completed on 2025-01-18 09:58:48
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
0047
0048
0049
0050
0051
0052 #ifndef G4PAIPhotData_h
0053 #define G4PAIPhotData_h 1
0054
0055 #include <vector>
0056 #include "globals.hh"
0057 #include "G4PAIxSection.hh"
0058 #include "G4SandiaTable.hh"
0059
0060 class G4PhysicsLogVector;
0061 class G4PhysicsTable;
0062 class G4MaterialCutsCouple;
0063 class G4PAIPhotModel;
0064
0065 class G4PAIPhotData
0066 {
0067
0068 public:
0069
0070 explicit G4PAIPhotData(G4double tmin, G4double tmax, G4int verbose);
0071
0072 ~G4PAIPhotData();
0073
0074 void Initialise(const G4MaterialCutsCouple*, G4double cut, G4PAIPhotModel*);
0075
0076 G4double DEDXPerVolume(G4int coupleIndex, G4double scaledTkin,
0077 G4double cut) const;
0078
0079 G4double CrossSectionPerVolume(G4int coupleIndex, G4double scaledTkin,
0080 G4double tcut, G4double tmax) const;
0081
0082 G4double GetPlasmonRatio( G4int coupleIndex, G4double scaledTkin ) const;
0083
0084 G4double SampleAlongStepTransfer(G4int coupleIndex, G4double kinEnergy,
0085 G4double scaledTkin,
0086 G4double stepFactor) const;
0087 G4double SampleAlongStepPhotonTransfer(G4int coupleIndex, G4double kinEnergy,
0088 G4double scaledTkin,
0089 G4double stepFactor) const;
0090 G4double SampleAlongStepPlasmonTransfer(G4int coupleIndex, G4double kinEnergy,
0091 G4double scaledTkin,
0092 G4double stepFactor) const;
0093
0094 G4double SamplePostStepTransfer(G4int coupleIndex,
0095 G4double scaledTkin) const;
0096 G4double SamplePostStepPhotonTransfer(G4int coupleIndex,
0097 G4double scaledTkin) const;
0098 G4double SamplePostStepPlasmonTransfer(G4int coupleIndex,
0099 G4double scaledTkin) const;
0100
0101
0102 G4PAIPhotData & operator=(const G4PAIPhotData &right) = delete;
0103 G4PAIPhotData(const G4PAIPhotData&) = delete;
0104
0105 private:
0106
0107 G4double GetEnergyTransfer(G4int coupleIndex, size_t iPlace,
0108 G4double position) const;
0109 G4double GetEnergyPhotonTransfer(G4int coupleIndex, size_t iPlace,
0110 G4double position) const;
0111 G4double GetEnergyPlasmonTransfer(G4int coupleIndex, size_t iPlace,
0112 G4double position) const;
0113
0114 G4int fTotBin;
0115 G4double fLowestKineticEnergy;
0116 G4double fHighestKineticEnergy;
0117
0118 G4PhysicsLogVector* fParticleEnergyVector;
0119
0120 G4PAIxSection fPAIxSection;
0121 G4SandiaTable fSandia;
0122
0123 std::vector<G4PhysicsTable*> fPAIxscBank;
0124 std::vector<G4PhysicsTable*> fPAIphotonBank;
0125 std::vector<G4PhysicsTable*> fPAIplasmonBank;
0126
0127 std::vector<G4PhysicsTable*> fPAIdEdxBank;
0128 std::vector<G4PhysicsLogVector*> fdEdxTable;
0129
0130 std::vector<G4PhysicsLogVector*> fdNdxCutTable;
0131 std::vector<G4PhysicsLogVector*> fdNdxCutPhotonTable;
0132 std::vector<G4PhysicsLogVector*> fdNdxCutPlasmonTable;
0133
0134 std::vector<G4PhysicsLogVector*> fdEdxCutTable;
0135 };
0136
0137 #endif
0138
0139
0140
0141
0142
0143
0144