| 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 // G4DecayTable class implementation
27 //
28 // Author: H.Kurashige, 7 July 1996
29 // --------------------------------------------------------------------
30
31 #include "G4DecayTable.hh"
32
33 #include "Randomize.hh"
34 #include "globals.hh"
35
36 G4DecayTable::G4DecayTable()
37 {
38 channels = new G4VDecayChannelVector;
39 }
40
41 G4DecayTable::~G4DecayTable()
42 {
43 // remove and delete all contents
44 for (const auto channel : *channels) {
45 delete channel;
46 }
47 channels->clear();
48 delete channels;
49 channels = nullptr;
50 parent = nullptr;
51 }
52
53 void G4DecayTable::Insert(G4VDecayChannel* aChannel)
54 {
55 if (parent == nullptr) {
56 parent = (G4ParticleDefinition*)(aChannel->GetParent());
57 }
58 if (parent != aChannel->GetParent()) {
59 #ifdef G4VERBOSE
60 G4cout << " G4DecayTable::Insert :: bad G4VDecayChannel (mismatch parent) "
61 << " " << parent->GetParticleName()
62 << " input:" << aChannel->GetParent()->GetParticleName() << G4endl;
63 #endif
64 }
65 else {
66 G4double br = aChannel->GetBR();
67 for (auto iCh = channels->cbegin(); iCh != channels->cend(); ++iCh) {
68 if (br > (*iCh)->GetBR()) {
69 channels->insert(iCh, aChannel);
70 return;
71 }
72 }
73 channels->push_back(aChannel);
74 }
75 }
76
77 G4VDecayChannel* G4DecayTable::SelectADecayChannel(G4double parentMass)
78 {
79 // check if contents exist
80 if (channels->empty()) return nullptr;
81
82 if (parentMass < 0.) parentMass = parent->GetPDGMass();
83
84 G4double sumBR = 0.;
85 for (const auto channel : *channels) {
86 if (!(channel->IsOKWithParentMass(parentMass))) continue;
87 sumBR += channel->GetBR();
88 }
89 if (sumBR <= 0.0) {
90 #ifdef G4VERBOSE
91 G4cout << " G4DecayTable::SelectADecayChannel :: no possible DecayChannel"
92 << " " << parent->GetParticleName() << G4endl;
93 #endif
94 return nullptr;
95 }
96
97 const std::size_t MAX_LOOP = 10000;
98 for (std::size_t loop_counter = 0; loop_counter < MAX_LOOP; ++loop_counter) {
99 G4double sum = 0.0;
100 G4double br = sumBR * G4UniformRand();
101 // select decay channel
102 for (const auto channel : *channels) {
103 sum += channel->GetBR();
104 if (!(channel->IsOKWithParentMass(parentMass))) continue;
105 if (br < sum) return channel;
106 }
107 }
108 return nullptr;
109 }
110
111 void G4DecayTable::DumpInfo() const
112 {
113 G4cout << "G4DecayTable: " << parent->GetParticleName() << G4endl;
114 G4int index = 0;
115 for (const auto channel : *channels) {
116 G4cout << index << ": ";
117 channel->DumpInfo();
118 index += 1;
119 }
120 G4cout << G4endl;
121 }
122