source: trunk/source/processes/hadronic/models/coherent_elastic/include/G4HadronElastic.hh @ 1315

Last change on this file since 1315 was 1315, checked in by garnier, 14 years ago

update geant4-09-04-beta-cand-01 interfaces-V09-03-09 vis-V09-03-08

File size: 4.8 KB
Line 
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,v 1.32 2010/01/13 15:42:06 mkossov Exp $
28// GEANT4 tag $Name: geant4-09-04-beta-cand-01 $
29//
30//
31// G4 Model: Low energy elastic scattering with 4-momentum balance
32// Derived fron G4LElastic of F.W. Jones, TRIUMF, 04-JUN-96
33// Uses  G4ElasticHadrNucleusHE and G4VQCrossSection
34// 
35// Modified:
36// 14-Dec-05 V.Ivanchenko rename the class
37// 13-Apr-06 V.Ivanchenko move to coherent_elastic
38// 25-Jul-06 V.Ivanchenko add 19 MeV low energy, below which S-wave is sampled
39// 20-Oct-06 V.Ivanchenko default ekinhigh = GeV (use HE model)
40// 16-Nov-06 V.Ivanchenko remove definition
41// 16-Nov-06 V.Ivanchenko default ekinhigh = 0.4*GeV
42// 16-Nov-06 V.Ivanchenko cleanup and rename Set methods and variables
43// 28-Mar-07 V.Ivanchenko add NIST manager
44// 11-May-07 V.Ivanchenko remove unused method Defs1
45// 13.01.10: M.Kosov: Use G4Q(Pr/Neut)ElasticCS instead of G4QElasticCS
46//
47//
48// Class Description
49// Final state production model for hadron nuclear elastic scattering;
50// Class Description - End
51
52
53#ifndef G4HadronElastic_h
54#define G4HadronElastic_h 1
55 
56#include "globals.hh"
57#include "G4HadronicInteraction.hh"
58#include "G4HadProjectile.hh"
59#include "G4Nucleus.hh"
60
61enum G4ElasticGenerator
62{
63  fLElastic = 0,
64  fHElastic,
65  fQElastic,
66  fSWave
67};
68
69class G4ParticleDefinition;
70class G4VQCrossSection;
71class G4ElasticHadrNucleusHE;
72
73class G4HadronElastic : public G4HadronicInteraction
74{
75public:
76
77  G4HadronElastic(G4ElasticHadrNucleusHE* HModel = 0);
78
79  virtual ~G4HadronElastic();
80 
81  G4HadFinalState * ApplyYourself(const G4HadProjectile & aTrack, 
82                                  G4Nucleus & targetNucleus);
83
84  G4VQCrossSection* GetCS();
85
86  G4ElasticHadrNucleusHE* GetHElastic();
87
88  void SetPlabLowLimit(G4double value);
89
90  void SetHEModelLowLimit(G4double value);
91
92  void SetQModelLowLimit(G4double value);
93
94  void SetLowestEnergyLimit(G4double value);
95
96  void SetRecoilKinEnergyLimit(G4double value);
97
98  G4double SampleT(G4double p, G4double m1, G4double m2, G4double A);
99
100private:
101
102  G4int Rtmi(G4double* x, G4double xli, G4double xri, G4double eps, 
103             G4int iend,
104             G4double aa, G4double bb, G4double cc, G4double dd, 
105             G4double rr);
106
107  G4double Fctcos(G4double t, 
108                  G4double aa, G4double bb, G4double cc, G4double dd, 
109                  G4double rr);
110
111  static G4VQCrossSection* pCManager;
112  static G4VQCrossSection* nCManager;
113
114  G4ElasticHadrNucleusHE*     hElastic;
115
116  G4ParticleDefinition* theProton;
117  G4ParticleDefinition* theNeutron;
118  G4ParticleDefinition* theDeuteron;
119  G4ParticleDefinition* theAlpha;
120  const G4ParticleDefinition* thePionPlus;
121  const G4ParticleDefinition* thePionMinus;
122
123  G4double lowEnergyRecoilLimit; 
124  G4double lowEnergyLimitHE; 
125  G4double lowEnergyLimitQ; 
126  G4double lowestEnergyLimit; 
127  G4double plabLowLimit;
128
129};
130
131inline void G4HadronElastic::SetRecoilKinEnergyLimit(G4double value)
132{
133  lowEnergyRecoilLimit = value;
134}
135
136inline void G4HadronElastic::SetPlabLowLimit(G4double value)
137{
138  plabLowLimit = value;
139}
140
141inline void G4HadronElastic::SetHEModelLowLimit(G4double value)
142{
143  lowEnergyLimitHE = value;
144}
145
146inline void G4HadronElastic::SetQModelLowLimit(G4double value)
147{
148  lowEnergyLimitQ = value;
149}
150
151inline void G4HadronElastic::SetLowestEnergyLimit(G4double value)
152{
153  lowestEnergyLimit = value;
154}
155
156#endif
Note: See TracBrowser for help on using the repository browser.