Geant4  10.01
G4ChipsElasticModel.cc
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 // $Id: G4ChipsElasticModel.cc 66892 2013-01-17 10:57:59Z gunter $
27 //
28 //---------------------------------------------------------------------
29 //
30 // Geant4 Class : G4ChipsElasticModel
31 //
32 // Author : V.Ivanchenko 29 June 2009
33 //
34 // Modified:
35 // 13.01.10: M.Kosov: Use G4Q(Pr/Neut)ElasticCS instead of G4QElasticCS
36 //
37 //---------------------------------------------------------------------
38 // CHIPS model of hadron elastic scattering
39 //
40 
41 #include "G4ChipsElasticModel.hh"
42 #include "G4ParticleDefinition.hh"
43 #include <iostream>
44 
46 
48 {
56  //Description();
57 }
58 
60 {}
61 
63 {
64  char* dirName = getenv("G4PhysListDocDir");
65  if (dirName) {
66  std::ofstream outFile;
67  G4String outFileName = GetModelName() + ".html";
68  G4String pathName = G4String(dirName) + "/" + outFileName;
69  outFile.open(pathName);
70  outFile << "<html>\n";
71  outFile << "<head>\n";
72 
73  outFile << "<title>Description of G4ChipsElasticModel</title>\n";
74  outFile << "</head>\n";
75  outFile << "<body>\n";
76 
77  outFile << "The G4ChipsElasticModel model performs hadron-nucleus elastic\n"
78  << "scattering using the parameterized elastic cross sections\n"
79  << "of M. Kossov\n";
80 
81  outFile << "</body>\n";
82  outFile << "</html>\n";
83  outFile.close();
84  }
85 }
86 
87 
88 G4double
90  G4double plab, G4int Z, G4int A)
91 {
92  G4int N = A - Z;
93  if(Z == 1 && N == 2) { N = 1; }
94  else if(Z == 2 && N == 1) { N = 2; }
95  G4int projPDG = p->GetPDGEncoding();
96  G4double cs = 0.;
97  if (projPDG==2212) { cs = pxsManager->GetChipsCrossSection(plab,Z,N,projPDG); }
98  else if(projPDG==2112) { cs = nxsManager->GetChipsCrossSection(plab,Z,N,projPDG); }
99  else if(projPDG==-2212){ cs = PBARxsManager->GetChipsCrossSection(plab,Z,N,projPDG); } //Pbar
100  else if(projPDG== 211) { cs = PIPxsManager->GetChipsCrossSection(plab,Z,N,projPDG); } // Pi+
101  else if(projPDG==-211) { cs = PIMxsManager->GetChipsCrossSection(plab,Z,N,projPDG); } // Pi-
102  else if(projPDG== 321) { cs = KPxsManager->GetChipsCrossSection(plab,Z,N,projPDG); } // K+
103  else if(projPDG==-321) { cs = KMxsManager->GetChipsCrossSection(plab,Z,N,projPDG); } // K-
104 
105  G4double t = 0.0;
106  if(cs > 0.0)
107  {
108  if (projPDG== 2212) { t = pxsManager->GetExchangeT(Z,N,projPDG); }
109  else if(projPDG== 2112) { t = nxsManager->GetExchangeT(Z,N,projPDG); }
110  else if(projPDG==-2212) { t = PBARxsManager->GetExchangeT(Z,N,projPDG); } // Pbar
111  else if(projPDG== 211) { t = PIPxsManager->GetExchangeT(Z,N,projPDG); } // Pi+
112  else if(projPDG== -211) { t = PIMxsManager->GetExchangeT(Z,N,projPDG); } // Pi-
113  else if(projPDG== 321) { t = KPxsManager->GetExchangeT(Z,N,projPDG); } // K+
114  else if(projPDG== -321) { t = KMxsManager->GetExchangeT(Z,N,projPDG); } // K-
115  }
116  else { t = G4HadronElastic::SampleInvariantT(p, plab, Z, A); }
117  return t;
118 }
119 
G4VCrossSectionDataSet * GetCrossSectionDataSet(const G4String &name, G4bool warning=true)
virtual void Description() const
virtual G4double SampleInvariantT(const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A)
virtual G4double GetChipsCrossSection(G4double momentum, G4int Z, G4int N, G4int pdg)
virtual G4double SampleInvariantT(const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A)
G4ChipsKaonPlusElasticXS * KPxsManager
G4double GetExchangeT(G4int tZ, G4int tN, G4int pPDG)
G4ChipsNeutronElasticXS * nxsManager
G4double GetExchangeT(G4int tZ, G4int tN, G4int pPDG)
G4ChipsPionMinusElasticXS * PIMxsManager
virtual G4double GetChipsCrossSection(G4double momentum, G4int Z, G4int N, G4int pdg)
const G4String & GetModelName() const
int G4int
Definition: G4Types.hh:78
G4ChipsKaonMinusElasticXS * KMxsManager
G4ChipsPionPlusElasticXS * PIPxsManager
G4ChipsAntiBaryonElasticXS * PBARxsManager
G4double GetExchangeT(G4int tZ, G4int tN, G4int pPDG)
static const char * Default_Name()
static G4CrossSectionDataSetRegistry * Instance()
static const char * Default_Name()
virtual G4double GetChipsCrossSection(G4double momentum, G4int Z, G4int N, G4int pdg)
static const G4double A[nN]
virtual G4double GetChipsCrossSection(G4double momentum, G4int Z, G4int N, G4int pdg)
G4double GetExchangeT(G4int tZ, G4int tN, G4int pPDG)
virtual G4double GetChipsCrossSection(G4double momentum, G4int Z, G4int N, G4int pdg)
virtual G4double GetChipsCrossSection(G4double momentum, G4int Z, G4int N, G4int pdg)
static const char * Default_Name()
G4ChipsProtonElasticXS * pxsManager
double G4double
Definition: G4Types.hh:76
G4double GetExchangeT(G4int tZ, G4int tN, G4int pPDG)
virtual G4double GetChipsCrossSection(G4double momentum, G4int Z, G4int N, G4int pdg)
G4double GetExchangeT(G4int tZ, G4int tN, G4int pPDG)
static const char * Default_Name()
G4double GetExchangeT(G4int tZ, G4int tN, G4int pPDG)