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
00030
00031
00032
00033
00034 #include "G4ios.hh"
00035
00036 #include "G4PiMinusStopLi.hh"
00037
00038 #include <vector>
00039
00040 #include "globals.hh"
00041 #include "Randomize.hh"
00042 #include "G4Proton.hh"
00043 #include "G4Neutron.hh"
00044 #include "G4ParticleTypes.hh"
00045 #include "G4ReactionKinematics.hh"
00046 #include "G4DynamicParticleVector.hh"
00047 #include "G4LorentzVector.hh"
00048 #include "G4PiMinusStopMaterial.hh"
00049 #include "G4DistributionGenerator.hh"
00050
00051
00052
00053
00054 G4double G4PiMinusStopLi::npRatio = 6.3;
00055
00056
00057
00058 G4double G4PiMinusStopLi::nFinalNucleons = 1.9;
00059
00060
00061
00062
00063
00064 G4int G4PiMinusStopLi::eKinEntries = 21;
00065
00066 G4double G4PiMinusStopLi::eKinData[21] = { 0.0018, 0.0025, 0.003, 0.045,
00067 0.007, 0.014, 0.023,
00068 0.4, 0.09, 0.18,
00069 0.25, 0.3, 0.25, 0.2,
00070 0.18, 0.08, 0.05,
00071 0.023, 0.012, 0.007, 0.02};
00072
00073 G4double G4PiMinusStopLi::eKin[22] = { 15., 17.5, 22.5, 27.5,
00074 32.5, 37.5, 42.5,
00075 47.5, 52.5, 57.5,
00076 62.5, 67.5, 72.5, 77.5,
00077 82.5, 87.5, 92.5,
00078 97.5, 102.5, 105. };
00079
00080
00081
00082
00083
00084 G4int G4PiMinusStopLi::angleEntries = 7;
00085
00086 G4double G4PiMinusStopLi::angleData[7] =
00087 { 0.17, 0.4, 0.7, 1.1, 1.3, 20., 70. };
00088
00089 G4double G4PiMinusStopLi::angle[8] = { 1.308997, 1.570796, 1.832596, 2.094395,
00090 2.356194, 2.617994, 2.967060, 3.1415927 };
00091
00092
00093
00094
00095
00096 G4PiMinusStopLi::G4PiMinusStopLi()
00097
00098 {
00099
00100
00101 _clusterSize = 2;
00102
00103
00104 theR = 1. / (1. + npRatio);
00105
00106 _definitions = new std::vector<G4ParticleDefinition*>();
00107 _momenta = new std::vector<G4LorentzVector*>();
00108
00109 std::vector<double> eKinVector;
00110 std::vector<double> eKinDataVector;
00111 int i;
00112 for (i=0; i<eKinEntries; i++)
00113 {
00114 eKinVector.push_back(eKin[i]);
00115 eKinDataVector.push_back(eKinData[i]);
00116 }
00117 eKinVector.push_back(eKin[eKinEntries]);
00118 _distributionE = new G4DistributionGenerator(eKinVector,eKinDataVector);
00119
00120 std::vector<double> angleVector;
00121 std::vector<double> angleDataVector;
00122 for (i=0; i<angleEntries; i++)
00123 {
00124 angleVector.push_back(angle[i]);
00125 angleDataVector.push_back(angleData[i]);
00126 }
00127 angleVector.push_back(angle[angleEntries]);
00128 _distributionAngle = new G4DistributionGenerator(angleVector,angleDataVector);
00129 }
00130
00131
00132
00133
00134 G4PiMinusStopLi::~G4PiMinusStopLi()
00135 {}
00136
00137 G4double G4PiMinusStopLi::FinalNucleons()
00138 {
00139 return nFinalNucleons;
00140 }
00141