Geant4  10.02.p01
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 90228 2015-05-21 08:49:57Z gcosmo $
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 
62 void G4ChipsElasticModel::ModelDescription(std::ostream& outFile) const
63 {
64 
65  outFile << "The G4ChipsElasticModel model performs hadron-nucleus elastic\n"
66  << "scattering using the parameterized elastic cross sections\n"
67  << "of M. Kossov\n";
68 
69 }
70 
71 
72 G4double
74  G4double plab, G4int Z, G4int A)
75 {
76  G4int N = A - Z;
77  if(Z == 1 && N == 2) { N = 1; }
78  else if(Z == 2 && N == 1) { N = 2; }
79  G4int projPDG = p->GetPDGEncoding();
80  G4double cs = 0.;
81  if (projPDG==2212) { cs = pxsManager->GetChipsCrossSection(plab,Z,N,projPDG); }
82  else if(projPDG==2112) { cs = nxsManager->GetChipsCrossSection(plab,Z,N,projPDG); }
83  else if(projPDG==-2212){ cs = PBARxsManager->GetChipsCrossSection(plab,Z,N,projPDG); } //Pbar
84  else if(projPDG== 211) { cs = PIPxsManager->GetChipsCrossSection(plab,Z,N,projPDG); } // Pi+
85  else if(projPDG==-211) { cs = PIMxsManager->GetChipsCrossSection(plab,Z,N,projPDG); } // Pi-
86  else if(projPDG== 321) { cs = KPxsManager->GetChipsCrossSection(plab,Z,N,projPDG); } // K+
87  else if(projPDG==-321) { cs = KMxsManager->GetChipsCrossSection(plab,Z,N,projPDG); } // K-
88 
89  G4double t = 0.0;
90  if(cs > 0.0)
91  {
92  if (projPDG== 2212) { t = pxsManager->GetExchangeT(Z,N,projPDG); }
93  else if(projPDG== 2112) { t = nxsManager->GetExchangeT(Z,N,projPDG); }
94  else if(projPDG==-2212) { t = PBARxsManager->GetExchangeT(Z,N,projPDG); } // Pbar
95  else if(projPDG== 211) { t = PIPxsManager->GetExchangeT(Z,N,projPDG); } // Pi+
96  else if(projPDG== -211) { t = PIMxsManager->GetExchangeT(Z,N,projPDG); } // Pi-
97  else if(projPDG== 321) { t = KPxsManager->GetExchangeT(Z,N,projPDG); } // K+
98  else if(projPDG== -321) { t = KMxsManager->GetExchangeT(Z,N,projPDG); } // K-
99  }
100  else { t = G4HadronElastic::SampleInvariantT(p, plab, Z, A); }
101  return t;
102 }
103 
G4VCrossSectionDataSet * GetCrossSectionDataSet(const G4String &name, G4bool warning=true)
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)
virtual void ModelDescription(std::ostream &) const
G4ChipsPionMinusElasticXS * PIMxsManager
virtual G4double GetChipsCrossSection(G4double momentum, G4int Z, G4int N, G4int pdg)
int G4int
Definition: G4Types.hh:78
G4ChipsKaonMinusElasticXS * KMxsManager
double A(double temperature)
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)
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)