Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
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 100378 2016-10-19 15:03:27Z 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 
56  explicit G4PreCompoundFragmentVector(pcfvector * avector);
57 
59 
60  void SetVector(pcfvector * avector);
61 
62  void SetOPTxs(G4int);
63 
64  void UseSICB(G4bool);
65 
66  inline G4double CalculateProbabilities(const G4Fragment & aFragment);
67 
69 
70 private:
71 
74  operator=(const G4PreCompoundFragmentVector &right) = delete;
75  G4bool operator==(const G4PreCompoundFragmentVector &right) const = delete;
76  G4bool operator!=(const G4PreCompoundFragmentVector &right) const = delete;
77 
78  pcfvector * theChannels;
79  G4DataVector probabilities;
80 
81  G4int nChannels;
82 };
83 
84 inline G4double
86 {
87  //G4cout << "## G4PreCompoundFragmentVector::CalculateProbabilities" << G4endl;
88  G4double probtot = 0.0;
89  G4double prob;
90  for (G4int i=0; i< nChannels; ++i) {
91  (*theChannels)[i]->Initialize(aFragment);
92  prob = 0.0;
93  if ((*theChannels)[i]->IsItPossible(aFragment)) {
94  prob = (*theChannels)[i]->CalcEmissionProbability(aFragment);
95  }
96  probtot += prob;
97  probabilities[i] = probtot;
98  //G4cout<<" prob= "<<prob<<" for "<<(*theChannels)[i]->GetName()<<G4endl;
99  }
100  return probtot;
101 }
102 
104 {
105  G4double x = probabilities[nChannels-1]*G4UniformRand();
106  G4int i=0;
107  for (; i<nChannels; ++i) {
108  if(x <= probabilities[i]) { break; }
109  }
110  return (*theChannels)[i];
111 }
112 
113 #endif
114 
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:97
bool G4bool
Definition: G4Types.hh:79
G4VPreCompoundFragment * ChooseFragment()
double G4double
Definition: G4Types.hh:76