Geant4_10
G4PreCompoundFragmentVector.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 // $Id: G4PreCompoundFragmentVector.hh 68028 2013-03-13 13:48:15Z gcosmo $
27 //
28 // Hadronic Process: Nuclear Preequilibrium
29 // by V. Lara
30 //
31 // Modified:
32 // 03.09.2008 by J. M. Quesada for external choice of inverse
33 // cross section option
34 // 06.09.2008 JMQ Also external choice has been added for:
35 // - superimposed Coulomb barrier (if useSICB=true)
36 // 27.08.2010 V.Ivanchenko simplify and make more efficient by adding extra
37 // vector of probabilities, moved constructor and destructor to source,
38 // simplify run time computations making inlined
39 //
40 
41 #ifndef G4PreCompoundFragmentVector_h
42 #define G4PreCompoundFragmentVector_h 1
43 
45 #include "G4DataVector.hh"
46 #include "Randomize.hh"
47 #include "globals.hh"
48 #include <vector>
49 
50 typedef std::vector<G4VPreCompoundFragment*> pcfvector;
51 
53 {
54 public:
55 
57 
59 
60  void SetVector(pcfvector * avector);
61 
62  void SetOPTxs(G4int);
63 
64  void UseSICB(G4bool);
65 
66  inline void Initialize(const G4Fragment & aFragment);
67 
68  inline G4double CalculateProbabilities(const G4Fragment & aFragment);
69 
71 
72 private:
73 
76  operator=(const G4PreCompoundFragmentVector &right);
77  G4bool operator==(const G4PreCompoundFragmentVector &right) const;
78  G4bool operator!=(const G4PreCompoundFragmentVector &right) const;
79 
80  pcfvector * theChannels;
81  G4DataVector probabilities;
82 
83  G4int nChannels;
84 
85 };
86 
87 inline void
89 {
90  for (G4int i=0; i< nChannels; ++i) {
91  (*theChannels)[i]->Initialize(aFragment);
92  }
93 }
94 
95 inline G4double
97 {
98  //G4cout << "## G4PreCompoundFragmentVector::CalculateProbabilities" << G4endl;
99  G4double probtot = 0.0;
100  for (G4int i=0; i< nChannels; ++i) {
101  G4double prob = (*theChannels)[i]->CalcEmissionProbability(aFragment);
102  probtot += prob;
103  probabilities[i] = probtot;
104  //G4cout<<" prob= "<<prob<<" for "<<(*theChannels)[i]->GetName()<<G4endl;
105  }
106  return probtot;
107 }
108 
110 {
111  G4double x = probabilities[nChannels-1]*G4UniformRand();
112  G4int i=0;
113  for (; i<nChannels; ++i) {
114  if(x <= probabilities[i]) { break; }
115  }
116  return (*theChannels)[i];
117 }
118 
119 #endif
120 
G4double CalculateProbabilities(const G4Fragment &aFragment)
tuple x
Definition: test.py:50
int G4int
Definition: G4Types.hh:78
std::vector< G4VPreCompoundFragment * > pcfvector
#define G4UniformRand()
Definition: Randomize.hh:87
bool G4bool
Definition: G4Types.hh:79
G4VPreCompoundFragment * ChooseFragment()
double G4double
Definition: G4Types.hh:76
void Initialize(const G4Fragment &aFragment)