| Geant4 Cross Reference |
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 // Hadronic Process: Nuclear Preequilibrium
28 // by V. Lara
29 //
30 // Modified:
31 // 27.08.2010 V.Ivanchenko moved constructor and destructor to source,
32 // simplify run time computations making inlined
33 //
34
35 #include "G4PreCompoundFragmentVector.hh"
36
37 G4PreCompoundFragmentVector::G4PreCompoundFragmentVector(pcfvector * avector)
38 : theChannels(nullptr), nChannels(0)
39 {
40 SetVector(avector);
41 }
42
43 void G4PreCompoundFragmentVector::SetVector(pcfvector * avector)
44 {
45 if(avector != theChannels) {
46 delete theChannels;
47 theChannels = avector;
48 }
49 if(theChannels) {
50 nChannels = (G4int)theChannels->size();
51 probabilities.resize(nChannels, 0.0);
52 } else {
53 nChannels = 0;
54 probabilities.clear();
55 }
56 }
57
58 //for inverse cross section choice
59 void G4PreCompoundFragmentVector::SetOPTxs(G4int opt)
60 {
61 for (G4int i=0; i<nChannels; ++i) {
62 (*theChannels)[i]->SetOPTxs(opt);
63 }
64 }
65
66 //for superimposed Coulomb Barrier for inverse cross sections
67 void G4PreCompoundFragmentVector::UseSICB(G4bool use)
68 {
69 for (G4int i=0; i< nChannels; ++i) {
70 (*theChannels)[i]->UseSICB(use);
71 }
72 }
73
74 G4double G4PreCompoundFragmentVector::CalculateProbabilities(
75 const G4Fragment & aFragment)
76 {
77 //G4cout << "## G4PreCompoundFragmentVector::CalculateProbabilities nCh= "
78 // << nChannels << G4endl;
79 G4double probtot = 0.0;
80 for (G4int i=0; i<nChannels; ++i) {
81 if ((*theChannels)[i]->Initialize(aFragment)) {
82 G4double prob = (*theChannels)[i]->CalcEmissionProbability(aFragment);
83 probtot += prob;
84 }
85 probabilities[i] = probtot;
86 //G4cout<< " probtot= " << probtot << " for "<< i << "-th channel" <<G4endl;
87 }
88 return probtot;
89 }
90
91 G4VPreCompoundFragment* G4PreCompoundFragmentVector::ChooseFragment()
92 {
93 //G4cout << "## G4PreCompoundFragmentVector::ChooseFragment nCh= "
94 // << nChannels << G4endl;
95 G4double x = probabilities[nChannels-1]*G4UniformRand();
96 G4int i=0;
97 for (; i<nChannels; ++i) {
98 if(x <= probabilities[i]) { break; }
99 }
100 return (*theChannels)[i];
101 }
102
103