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
00035
00036 #include "G4QHadronElasticDataSet.hh"
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049 G4QHadronElasticDataSet::G4QHadronElasticDataSet(const G4String& dataSetName)
00050 : G4VCrossSectionDataSet(dataSetName)
00051 {
00052
00053 Description();
00054 }
00055
00056 void G4QHadronElasticDataSet::Description() const
00057 {
00058 char* dirName = getenv("G4PhysListDocDir");
00059 if (dirName) {
00060 std::ofstream outFile;
00061 G4String outFileName = GetName() + ".html";
00062 G4String pathName = G4String(dirName) + "/" + outFileName;
00063
00064 outFile.open(pathName);
00065 outFile << "<html>\n";
00066 outFile << "<head>\n";
00067
00068 outFile << "<title>Description of CHIPSElasticXS</title>\n";
00069 outFile << "</head>\n";
00070 outFile << "<body>\n";
00071
00072 outFile << "CHIPSElasticXS provides hadron-nuclear elastic cross\n"
00073 << "sections for all hadrons at all energies. These cross\n"
00074 << "sections represent parameterizations developed by M. Kossov.\n";
00075
00076 outFile << "</body>\n";
00077 outFile << "</html>\n";
00078 outFile.close();
00079 }
00080 }
00081
00082 G4bool G4QHadronElasticDataSet::IsIsoApplicable(const G4DynamicParticle* Pt, G4int, G4int,
00083 const G4Element*, const G4Material*)
00084 {
00085 G4ParticleDefinition* particle = Pt->GetDefinition();
00086 if (particle == G4Neutron::Neutron() ) return true;
00087 else if (particle == G4Proton::Proton() ) return true;
00088 else if (particle == G4PionMinus::PionMinus() ) return true;
00089 else if (particle == G4PionPlus::PionPlus() ) return true;
00090 else if (particle == G4KaonPlus::KaonPlus() ) return true;
00091 else if (particle == G4KaonMinus::KaonMinus() ) return true;
00092 else if (particle == G4KaonZeroLong::KaonZeroLong() ) return true;
00093 else if (particle == G4KaonZeroShort::KaonZeroShort() ) return true;
00094 else if (particle == G4Lambda::Lambda() ) return true;
00095 else if (particle == G4SigmaPlus::SigmaPlus() ) return true;
00096 else if (particle == G4SigmaMinus::SigmaMinus() ) return true;
00097 else if (particle == G4SigmaZero::SigmaZero() ) return true;
00098 else if (particle == G4XiMinus::XiMinus() ) return true;
00099 else if (particle == G4XiZero::XiZero() ) return true;
00100 else if (particle == G4OmegaMinus::OmegaMinus() ) return true;
00101 else if (particle == G4AntiNeutron::AntiNeutron() ) return true;
00102 else if (particle == G4AntiProton::AntiProton() ) return true;
00103 else if (particle == G4AntiLambda::AntiLambda() ) return true;
00104 else if (particle == G4AntiSigmaPlus::AntiSigmaPlus() ) return true;
00105 else if (particle == G4AntiSigmaMinus::AntiSigmaMinus() ) return true;
00106 else if (particle == G4AntiSigmaZero::AntiSigmaZero() ) return true;
00107 else if (particle == G4AntiXiMinus::AntiXiMinus() ) return true;
00108 else if (particle == G4AntiXiZero::AntiXiZero() ) return true;
00109 else if (particle == G4AntiOmegaMinus::AntiOmegaMinus() ) return true;
00110 return false;
00111 }
00112
00113
00114
00115
00116
00117
00118
00119
00120
00121
00122
00123
00124
00125
00126
00127
00128
00129
00130
00131
00132
00133
00134
00135
00136
00137
00138
00139
00140
00141
00142
00143
00144
00145
00146
00147
00148
00149
00150
00151
00152
00153
00154
00155
00156
00157
00158
00159
00160
00161
00162
00163
00164
00165
00166
00167
00168
00169
00170
00171
00172
00173
00174
00175
00176
00177
00178
00179
00180
00181
00182 G4double G4QHadronElasticDataSet::GetIsoCrossSection(const G4DynamicParticle* Pt, G4int Z,
00183 G4int A, const G4Isotope*,
00184 const G4Element*, const G4Material*)
00185 {
00186 G4ParticleDefinition* particle = Pt->GetDefinition();
00187 G4double Momentum=Pt->GetTotalMomentum();
00188 G4VQCrossSection* CSmanager=0;
00189 G4VQCrossSection* CSmanager2=0;
00190
00191 G4int pPDG=0;
00192 if(particle == G4Neutron::Neutron())
00193 {
00194 CSmanager=G4QNeutronElasticCrossSection::GetPointer();
00195 pPDG=2112;
00196 }
00197 else if(particle == G4Proton::Proton())
00198 {
00199 CSmanager=G4QProtonElasticCrossSection::GetPointer();
00200 pPDG=2212;
00201 }
00202 else if(particle == G4PionMinus::PionMinus())
00203 {
00204 CSmanager=G4QPionMinusElasticCrossSection::GetPointer();
00205 pPDG=-211;
00206 }
00207 else if(particle == G4PionPlus::PionPlus())
00208 {
00209 CSmanager=G4QPionPlusElasticCrossSection::GetPointer();
00210 pPDG=211;
00211 }
00212 else if(particle == G4KaonMinus::KaonMinus())
00213 {
00214 CSmanager=G4QKaonMinusElasticCrossSection::GetPointer();
00215 pPDG=-321;
00216 }
00217 else if(particle == G4KaonPlus::KaonPlus())
00218 {
00219 CSmanager=G4QKaonPlusElasticCrossSection::GetPointer();
00220 pPDG=321;
00221 }
00222 else if(particle == G4KaonZeroLong::KaonZeroLong() ||
00223 particle == G4KaonZeroShort::KaonZeroShort() ||
00224 particle == G4KaonZero::KaonZero() ||
00225 particle == G4AntiKaonZero::AntiKaonZero() )
00226 {
00227 CSmanager=G4QKaonMinusElasticCrossSection::GetPointer();
00228 CSmanager2=G4QKaonPlusElasticCrossSection::GetPointer();
00229 if(G4UniformRand() > 0.5) pPDG= 321;
00230 else pPDG=-321;
00231 }
00232 else if(particle == G4Lambda::Lambda())
00233 {
00234 CSmanager=G4QHyperonElasticCrossSection::GetPointer();
00235 pPDG=3122;
00236 }
00237 else if(particle == G4SigmaPlus::SigmaPlus())
00238 {
00239 CSmanager=G4QHyperonElasticCrossSection::GetPointer();
00240 pPDG=3222;
00241 }
00242 else if(particle == G4SigmaMinus::SigmaMinus())
00243 {
00244 CSmanager=G4QHyperonElasticCrossSection::GetPointer();
00245 pPDG=3112;
00246 }
00247 else if(particle == G4SigmaZero::SigmaZero())
00248 {
00249 CSmanager=G4QHyperonElasticCrossSection::GetPointer();
00250 pPDG=3212;
00251 }
00252 else if(particle == G4XiMinus::XiMinus())
00253 {
00254 CSmanager=G4QHyperonElasticCrossSection::GetPointer();
00255 pPDG=3312;
00256 }
00257 else if(particle == G4XiZero::XiZero())
00258 {
00259 CSmanager=G4QHyperonElasticCrossSection::GetPointer();
00260 pPDG=3322;
00261 }
00262 else if(particle == G4OmegaMinus::OmegaMinus())
00263 {
00264 CSmanager=G4QHyperonElasticCrossSection::GetPointer();
00265 pPDG=3334;
00266 }
00267 else if(particle == G4AntiNeutron::AntiNeutron())
00268 {
00269 CSmanager=G4QAntiBaryonElasticCrossSection::GetPointer();
00270 pPDG=-2112;
00271 }
00272 else if(particle == G4AntiProton::AntiProton())
00273 {
00274 CSmanager=G4QAntiBaryonElasticCrossSection::GetPointer();
00275 pPDG=-2212;
00276 }
00277 else if(particle == G4AntiLambda::AntiLambda())
00278 {
00279 CSmanager=G4QAntiBaryonElasticCrossSection::GetPointer();
00280 pPDG=-3122;
00281 }
00282 else if(particle == G4AntiSigmaPlus::AntiSigmaPlus())
00283 {
00284 CSmanager=G4QAntiBaryonElasticCrossSection::GetPointer();
00285 pPDG=-3222;
00286 }
00287 else if(particle == G4AntiSigmaMinus::AntiSigmaMinus())
00288 {
00289 CSmanager=G4QAntiBaryonElasticCrossSection::GetPointer();
00290 pPDG=-3112;
00291 }
00292 else if(particle == G4AntiSigmaZero::AntiSigmaZero())
00293 {
00294 CSmanager=G4QAntiBaryonElasticCrossSection::GetPointer();
00295 pPDG=-3212;
00296 }
00297 else if(particle == G4AntiXiMinus::AntiXiMinus())
00298 {
00299 CSmanager=G4QAntiBaryonElasticCrossSection::GetPointer();
00300 pPDG=-3312;
00301 }
00302 else if(particle == G4AntiXiZero::AntiXiZero())
00303 {
00304 CSmanager=G4QAntiBaryonElasticCrossSection::GetPointer();
00305 pPDG=-3322;
00306 }
00307 else if(particle == G4AntiOmegaMinus::AntiOmegaMinus())
00308 {
00309 CSmanager=G4QAntiBaryonElasticCrossSection::GetPointer();
00310 pPDG=-3334;
00311 }
00312 else
00313 {
00314 G4cerr << "-ERROR-G4QHadronElasticDataSet::GetIsoZACrossSection: PDG="
00315 << particle->GetPDGEncoding() << " isn't supported by CHIPS" << G4endl;
00316
00317
00318 return 0;
00319 }
00320 G4int N=A-Z;
00321 G4double CSI=CSmanager->GetCrossSection(true, Momentum, Z, N, pPDG);
00322 if(CSmanager2) CSI = (CSI + CSmanager2->GetCrossSection(true, Momentum, Z, N, pPDG))/2;
00323 return CSI;
00324 }