Geant4_10
G4XMesonBaryonElastic.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 //
27 
28 #include "globals.hh"
29 #include "G4ios.hh"
30 #include "G4XMesonBaryonElastic.hh"
31 #include "G4KineticTrack.hh"
32 #include "G4Gamma.hh"
33 #include "G4PionPlus.hh"
34 #include "G4Proton.hh"
35 #include "G4XAqmElastic.hh"
36 #include "G4XPDGElastic.hh"
37 
39 {
40  // As a first approximation the model is assumed to be valid over
41  // the entire energy range
42  lowLimit = 0.;
43  highLimit = DBL_MAX;
44 }
45 
46 
48 { }
49 
50 
52 {
53  return (this == (G4XMesonBaryonElastic *) &right);
54 }
55 
56 
58 {
59  return (this != (G4XMesonBaryonElastic *) &right);
60 }
61 
62 
64 {
65  G4double sigma;
66 
67  // No gamma-baryon elastic scattering
68  const G4ParticleDefinition* defLight = FindLightParticle(trk1,trk2);
69  if (defLight == G4Gamma::GammaDefinition())
70  {
71  sigma = 0;
72  }
73  else
74  {
75 
76  G4LorentzVector p41 = trk1.Get4Momentum();
77  G4ThreeVector p3 = trk1.GetPosition();
79 
80  G4KineticTrack piTrk(def,
81  trk1.GetFormationTime(),
82  p3,
83  (G4LorentzVector&)p41);
84 
85  G4LorentzVector p42 = trk2.Get4Momentum();
87  trk2.GetFormationTime(),
88  (G4ThreeVector)trk2.GetPosition(),
89  (G4LorentzVector&)p42);
90 
91  G4XAqmElastic aqm;
92  G4double xAqmDummy = aqm.CrossSection(piTrk,pTrk);
93  G4double xAqm = aqm.CrossSection(trk1,trk2);
94  G4double factor = 1.;
95  if (xAqmDummy != 0.0)
96  {
97  factor = xAqm / xAqmDummy;
98  }
99  G4XPDGElastic pdg;
100 
101  sigma = pdg.CrossSection(piTrk,pTrk);
102  sigma = sigma * factor;
103  }
104 
105  return sigma;
106 }
107 
108 
110 {
111  G4String name("MesonBaryonElasticCrossSection");
112  return name;
113 }
114 
115 
116 
118 {
119  G4bool answer = InLimits(e,lowLimit,highLimit);
120 
121  return answer;
122 }
G4bool operator==(const G4XMesonBaryonElastic &right) const
CLHEP::Hep3Vector G4ThreeVector
const G4ThreeVector & GetPosition() const
static G4Proton * ProtonDefinition()
Definition: G4Proton.cc:88
const XML_Char * name
Definition: expat.h:151
virtual G4double CrossSection(const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
virtual G4String Name() const
G4double GetFormationTime() const
bool G4bool
Definition: G4Types.hh:79
static G4PionPlus * PionPlusDefinition()
Definition: G4PionPlus.cc:93
virtual G4double CrossSection(const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
G4bool operator!=(const G4XMesonBaryonElastic &right) const
virtual G4bool IsValid(G4double e) const
double G4double
Definition: G4Types.hh:76
virtual G4double CrossSection(const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
const G4LorentzVector & Get4Momentum() const
#define DBL_MAX
Definition: templates.hh:83
const G4ParticleDefinition * FindLightParticle(const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
G4bool InLimits(G4double e, G4double eLow, G4double eHigh) const
static G4Gamma * GammaDefinition()
Definition: G4Gamma.cc:81