00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029 #ifndef G4DNADingfelderChargeDecreaseModel_h
00030 #define G4DNADingfelderChargeDecreaseModel_h 1
00031
00032 #include "G4VEmModel.hh"
00033 #include "G4ParticleChangeForGamma.hh"
00034 #include "G4ProductionCutsTable.hh"
00035
00036 #include "G4Proton.hh"
00037 #include "G4DNAGenericIonsManager.hh"
00038 #include "G4NistManager.hh"
00039
00040 class G4DNADingfelderChargeDecreaseModel : public G4VEmModel
00041 {
00042
00043 public:
00044
00045 G4DNADingfelderChargeDecreaseModel(const G4ParticleDefinition* p = 0,
00046 const G4String& nam = "DNADingfelderChargeDecreaseModel");
00047
00048 virtual ~G4DNADingfelderChargeDecreaseModel();
00049
00050 virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&);
00051
00052 virtual G4double CrossSectionPerVolume(const G4Material* material,
00053 const G4ParticleDefinition* p,
00054 G4double ekin,
00055 G4double emin,
00056 G4double emax);
00057
00058 virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
00059 const G4MaterialCutsCouple*,
00060 const G4DynamicParticle*,
00061 G4double tmin,
00062 G4double maxEnergy);
00063
00064 protected:
00065
00066 G4ParticleChangeForGamma* fParticleChangeForGamma;
00067
00068 private:
00069
00070 const std::vector<G4double>* fpMolWaterDensity;
00071
00072 std::map<G4String,G4double,std::less<G4String> > lowEnergyLimit;
00073 std::map<G4String,G4double,std::less<G4String> > highEnergyLimit;
00074
00075 G4bool isInitialised;
00076 G4int verboseLevel;
00077
00078
00079
00080 G4double PartialCrossSection(G4double energy, G4int level, const G4ParticleDefinition* particle);
00081
00082 G4double Sum(G4double energy, const G4ParticleDefinition* particle);
00083
00084 G4int RandomSelect(G4double energy, const G4ParticleDefinition* particle);
00085
00086 G4int numberOfPartialCrossSections[3];
00087
00088 G4double f0[2][3];
00089 G4double a0[2][3];
00090 G4double a1[2][3];
00091 G4double b0[2][3];
00092 G4double b1[2][3];
00093 G4double c0[2][3];
00094 G4double d0[2][3];
00095 G4double x0[2][3];
00096 G4double x1[2][3];
00097
00098
00099
00100 G4int NumberOfFinalStates(G4ParticleDefinition* particleDefinition, G4int finalStateIndex);
00101
00102 G4ParticleDefinition* OutgoingParticleDefinition(G4ParticleDefinition* particleDefinition, G4int finalStateIndex);
00103
00104 G4double WaterBindingEnergyConstant(G4ParticleDefinition* particleDefinition, G4int finalStateIndex);
00105
00106 G4double OutgoingParticleBindingEnergyConstant(G4ParticleDefinition* particleDefinition, G4int finalStateIndex);
00107
00108
00109
00110 G4DNADingfelderChargeDecreaseModel & operator=(const G4DNADingfelderChargeDecreaseModel &right);
00111 G4DNADingfelderChargeDecreaseModel(const G4DNADingfelderChargeDecreaseModel&);
00112
00113 };
00114
00115
00116
00117 #endif