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 G4NeutronHPKallbachMannSyst_h
00030 #define G4NeutronHPKallbachMannSyst_h 1
00031
00032 #include "globals.hh"
00033
00034 class G4NeutronHPKallbachMannSyst
00035 {
00036 public:
00037
00038 G4NeutronHPKallbachMannSyst(G4double aCompoundFraction,
00039 G4double anIncidentEnergy, G4double anIncidentMass,
00040 G4double aProductEnergy, G4double aProductMass,
00041 G4double aResidualMass, G4int aResidualA, G4int aResidualZ,
00042 G4double aTargetMass, G4int aTargetA, G4int aTargetZ)
00043 {
00044 theCompoundFraction = aCompoundFraction;
00045 theIncidentEnergy = anIncidentEnergy;
00046 theIncidentMass = anIncidentMass;
00047 theProductEnergy = aProductEnergy;
00048 theProductMass = aProductMass;
00049 theResidualMass = aResidualMass;
00050 theResidualA = aResidualA;
00051 theResidualZ = aResidualZ;
00052 theTargetMass = aTargetMass;
00053 theTargetA = aTargetA;
00054 theTargetZ = aTargetZ;
00055 }
00056
00057 ~G4NeutronHPKallbachMannSyst() {};
00058
00059 G4double Sample(G4double anEnergy);
00060
00061 G4double Kallbach(G4double cosTh, G4double anEnergy);
00062
00063 G4double GetKallbachZero(G4double anEnergy);
00064
00065 G4double A(G4double anEnergy);
00066
00067 G4double SeparationEnergy(G4int Ac, G4int Nc, G4int AA, G4int ZA);
00068
00069 private:
00070
00071 G4double theCompoundFraction;
00072 G4double theIncidentEnergy;
00073 G4double theIncidentMass;
00074 G4double theProductEnergy;
00075 G4double theProductMass;
00076 G4double theResidualMass;
00077 G4double theTargetMass;
00078 G4int theResidualA;
00079 G4int theResidualZ;
00080 G4int theTargetA;
00081 G4int theTargetZ;
00082 };
00083
00084 #endif