Geant4_10
G4HadronElastic.hh
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 // $Id: G4HadronElastic.hh 66892 2013-01-17 10:57:59Z gunter $
28 //
29 // Geant4 Header : G4HadronElastic
30 //
31 // Author : V.Ivanchenko 29 June 2009 (redesign old elastic model)
32 //
33 // Modified:
34 //
35 // Class Description
36 // Default model for elastic scattering; GHEISHA algorithm is used
37 // Class Description - End
38 
39 #ifndef G4HadronElastic_h
40 #define G4HadronElastic_h 1
41 
42 #include "globals.hh"
43 #include "G4HadronicInteraction.hh"
44 #include "G4HadProjectile.hh"
45 #include "G4Nucleus.hh"
46 #include "G4NucleiProperties.hh"
47 
49 
51 {
52 public:
53 
54  G4HadronElastic(const G4String& name = "hElasticLHEP");
55 
56  virtual ~G4HadronElastic();
57 
58  // implementation of the G4HadronicInteraction interface
59  virtual G4HadFinalState * ApplyYourself(const G4HadProjectile & aTrack,
60  G4Nucleus & targetNucleus);
61 
62  // sample momentum transfer using Lab. momentum
64  G4double plab,
65  G4int Z, G4int A);
66 
67  inline void SetLowestEnergyLimit(G4double value);
68 
69  inline G4double LowestEnergyLimit() const;
70 
72  G4double plab, G4int Z, G4int A);
73 
74  virtual void Description() const;
75 
76 private:
77 
78  G4ParticleDefinition* theProton;
79  G4ParticleDefinition* theNeutron;
80  G4ParticleDefinition* theDeuteron;
81  G4ParticleDefinition* theAlpha;
82 
83  G4double lowestEnergyLimit;
84 
85 };
86 
88 {
89  lowestEnergyLimit = value;
90 }
91 
93 {
94  return lowestEnergyLimit;
95 }
96 
97 inline G4double
99  G4double plab, G4int Z, G4int A)
100 {
101  G4double m1 = p->GetPDGMass();
102  G4double m12= m1*m1;
104  return plab*mass2/std::sqrt(m12 + mass2*mass2 + 2.*mass2*std::sqrt(m12 + plab*plab));
105 }
106 
107 #endif
virtual G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
static G4double GetNuclearMass(const G4double A, const G4double Z)
virtual G4double SampleInvariantT(const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A)
G4HadronElastic(const G4String &name="hElasticLHEP")
const char * p
Definition: xmltok.h:285
const XML_Char * name
Definition: expat.h:151
int G4int
Definition: G4Types.hh:78
virtual void Description() const
Float_t Z
Definition: plot.C:39
G4double LowestEnergyLimit() const
G4double ComputeMomentumCMS(const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A)
G4double GetPDGMass() const
void SetLowestEnergyLimit(G4double value)
virtual ~G4HadronElastic()
const XML_Char int const XML_Char * value
Definition: expat.h:331
double G4double
Definition: G4Types.hh:76