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 "G4PiMinusStopCo.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 G4PiMinusStopCo::npRatio = 2.5;
00055
00056
00057
00058
00059
00060 G4double G4PiMinusStopCo::nFinalNucleons = 1.38;
00061
00062
00063
00064
00065
00066 G4int G4PiMinusStopCo::eKinEntries = 11;
00067
00068 G4double G4PiMinusStopCo::eKinData[11] = {0.085, 0.09, 0.09, 0.15,
00069 0.1, 0.09, 0.08,
00070 0.04, 0.03, 0.02, 0.01};
00071
00072 G4double G4PiMinusStopCo::eKin[12] = { 15., 17.5, 25., 33.,
00073 42., 52., 62.,
00074 75., 85., 95., 105. };
00075
00076
00077
00078
00079
00080 G4int G4PiMinusStopCo::angleEntries = 7;
00081
00082 G4double G4PiMinusStopCo::angleData[7] =
00083 {6., 8., 9., 10., 25., 40., 45. };
00084
00085 G4double G4PiMinusStopCo::angle[8] = { 0.5, 0.7, 0.87, 1.4, 2.1,
00086 2.44, 2.8, 3.1415927 };
00087
00088
00089
00090
00091
00092 G4PiMinusStopCo::G4PiMinusStopCo()
00093
00094 {
00095
00096
00097 _clusterSize = 2;
00098
00099
00100 theR = 1. / (1. + npRatio);
00101
00102 _definitions = new std::vector<G4ParticleDefinition*>();
00103 _momenta = new std::vector<G4LorentzVector*>();
00104
00105 std::vector<double> eKinVector;
00106 std::vector<double> eKinDataVector;
00107 int i;
00108 for (i=0; i<eKinEntries; i++)
00109 {
00110 eKinVector.push_back(eKin[i]);
00111 eKinDataVector.push_back(eKinData[i]);
00112 }
00113 eKinVector.push_back(eKin[eKinEntries]);
00114 _distributionE = new G4DistributionGenerator(eKinVector,eKinDataVector);
00115
00116 std::vector<double> angleVector;
00117 std::vector<double> angleDataVector;
00118 for (i=0; i<angleEntries; i++)
00119 {
00120 angleVector.push_back(angle[i]);
00121 angleDataVector.push_back(angleData[i]);
00122 }
00123 angleVector.push_back(angle[angleEntries]);
00124 _distributionAngle = new G4DistributionGenerator(angleVector,angleDataVector);
00125 }
00126
00127
00128
00129
00130 G4PiMinusStopCo::~G4PiMinusStopCo()
00131 {}
00132
00133 G4double G4PiMinusStopCo::FinalNucleons()
00134 {
00135 return nFinalNucleons;
00136 }
00137