source: trunk/source/processes/hadronic/models/de_excitation/fission/include/G4CompetitiveFission.hh @ 1340

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

update ti head

File size: 5.4 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: G4CompetitiveFission.hh,v 1.3 2006/06/29 20:13:19 gunter Exp $
28// GEANT4 tag $Name: geant4-09-03-ref-09 $
29//
30// Hadronic Process: Nuclear De-excitations
31// by V. Lara (Oct 1998)
32
33#ifndef G4CompetitiveFission_h
34#define G4CompetitiveFission_h 1
35
36#include "G4VEvaporationChannel.hh"
37#include "G4Fragment.hh"
38#include "G4VFissionBarrier.hh"
39#include "G4FissionBarrier.hh"
40#include "G4VEmissionProbability.hh"
41#include "G4FissionProbability.hh"
42#include "G4VLevelDensityParameter.hh"
43#include "G4FissionLevelDensityParameter.hh"
44#include "G4FissionParameters.hh"
45#include "G4ParticleTable.hh"
46#include "G4IonTable.hh"
47#include "Randomize.hh"
48
49//#define debug
50
51class G4CompetitiveFission : public G4VEvaporationChannel
52{
53public:
54 
55  G4CompetitiveFission();
56  virtual ~G4CompetitiveFission();
57
58private:
59  G4CompetitiveFission(const G4CompetitiveFission &right);
60
61  const G4CompetitiveFission & operator=(const G4CompetitiveFission &right);
62public:
63  G4bool operator==(const G4CompetitiveFission &right) const;
64  G4bool operator!=(const G4CompetitiveFission &right) const;
65
66public:
67  G4FragmentVector * BreakUp(const G4Fragment &theNucleus);
68
69  void Initialize(const G4Fragment & fragment);
70
71  inline void SetFissionBarrier(G4VFissionBarrier * aBarrier)
72  {
73    if (MyOwnFissionBarrier) delete theFissionBarrierPtr;
74    theFissionBarrierPtr = aBarrier;
75    MyOwnFissionBarrier = false;
76  }
77
78  inline void SetEmissionStrategy(G4VEmissionProbability * aFissionProb)
79  {
80    if (MyOwnFissionProbability) delete theFissionProbabilityPtr;
81    theFissionProbabilityPtr = aFissionProb;
82    MyOwnFissionProbability = false;
83  }
84
85
86  inline void SetLevelDensityParameter(G4VLevelDensityParameter * aLevelDensity)
87  { 
88    if (MyOwnLevelDensity) delete theLevelDensityPtr;
89    theLevelDensityPtr = aLevelDensity;
90    MyOwnLevelDensity = false;
91  }
92
93
94  inline G4double GetFissionBarrier(void) const { return FissionBarrier; }
95
96  inline G4double GetEmissionProbability(void) const { return FissionProbability; }
97
98  inline G4double GetLevelDensityParameter(void) const { return LevelDensityParameter; }
99
100  inline G4double GetMaximalKineticEnergy(void) const { return MaximalKineticEnergy; }
101private:
102
103  // Maximal Kinetic Energy that can be carried by fragment
104  G4double MaximalKineticEnergy;
105
106
107  // For Fission barrier
108  G4VFissionBarrier * theFissionBarrierPtr;
109  G4double FissionBarrier;
110  G4bool MyOwnFissionBarrier;
111
112  // For Fission probability emission
113  G4VEmissionProbability * theFissionProbabilityPtr;
114  G4double FissionProbability;
115  G4bool MyOwnFissionProbability;
116
117
118  // For Level Density calculation
119  G4bool MyOwnLevelDensity;
120  G4VLevelDensityParameter * theLevelDensityPtr;
121  G4double LevelDensityParameter;
122
123
124
125
126  //  --------------------
127
128
129  // Sample AtomicNumber of Fission products
130  G4int FissionAtomicNumber(const G4int A, const G4FissionParameters & theParam);
131  G4double MassDistribution(const G4double x, const G4double A, const G4FissionParameters & theParam);
132
133
134  // Sample Charge of fission products
135  G4int FissionCharge(const G4double A, const G4double Z, const G4double Af);
136
137
138  // Sample Kinetic energy of fission products
139  G4double FissionKineticEnergy(const G4double A, const G4double Z,
140                                const G4double Af1, const G4double Zf1,
141                                const G4double Af2, const G4double Zf2,
142                                const G4double U, const G4double Tmax,
143                                const G4FissionParameters & theParam);
144   
145
146
147  G4double Ratio(const G4double A,const G4double A11,const G4double B1,const G4double A00);
148  G4double SymmetricRatio(const G4double A,const G4double A11);
149  G4double AsymmetricRatio(const G4double A,const G4double A11);
150
151
152
153  G4ThreeVector IsotropicVector(const G4double Magnitude = 1.0);
154
155
156#ifdef debug
157  void CheckConservation(const G4Fragment & theInitialState,
158                         G4FragmentVector * Result) const;
159#endif
160
161
162};
163
164
165
166#endif
167
168
Note: See TracBrowser for help on using the repository browser.