| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // neutron_hp -- source file 26 // neutron_hp -- source file
27 // J.P. Wellisch, Nov-1996 27 // J.P. Wellisch, Nov-1996
28 // A prototype of the low energy neutron trans 28 // A prototype of the low energy neutron transport model.
29 // 29 //
30 // P. Arce, June-2014 Conversion neutron_hp to 30 // P. Arce, June-2014 Conversion neutron_hp to particle_hp
31 // 31 //
32 #include "G4ParticleHPField.hh" 32 #include "G4ParticleHPField.hh"
33 33
34 #include "G4HadronicException.hh" 34 #include "G4HadronicException.hh"
35 #include "G4ios.hh" 35 #include "G4ios.hh"
36 36
37 G4ParticleHPField::G4ParticleHPField() 37 G4ParticleHPField::G4ParticleHPField()
38 { 38 {
39 theData = new G4ParticleHPFieldPoint[100]; 39 theData = new G4ParticleHPFieldPoint[100];
40 nPoints = 100; 40 nPoints = 100;
41 nEntries = 0; 41 nEntries = 0;
42 theData->SetData(0, 0, 0); 42 theData->SetData(0, 0, 0);
43 } 43 }
44 44
45 G4ParticleHPField::~G4ParticleHPField() 45 G4ParticleHPField::~G4ParticleHPField()
46 { 46 {
47 delete[] theData; 47 delete[] theData;
48 } 48 }
49 49
50 G4double G4ParticleHPField::GetY(G4double e, G 50 G4double G4ParticleHPField::GetY(G4double e, G4int j)
51 { 51 {
52 G4int low = 0; 52 G4int low = 0;
53 G4int high = 0; 53 G4int high = 0;
54 G4int i; 54 G4int i;
55 for (i = 1; i < nEntries / 10; i++) { 55 for (i = 1; i < nEntries / 10; i++) {
56 if (theData[10 * i].GetX() > e) break; 56 if (theData[10 * i].GetX() > e) break;
57 } 57 }
58 if (i == (nEntries / 10)) { 58 if (i == (nEntries / 10)) {
59 i = 10 * i; 59 i = 10 * i;
60 while (i < nEntries) // Loop checking, 11 60 while (i < nEntries) // Loop checking, 11.05.2015, T. Koi
61 { 61 {
62 if (theData[i++].GetX() > e) break; 62 if (theData[i++].GetX() > e) break;
63 } 63 }
64 if (i == nEntries) { 64 if (i == nEntries) {
65 low = nEntries - 1; 65 low = nEntries - 1;
66 high = nEntries - 2; 66 high = nEntries - 2;
67 } 67 }
68 else { 68 else {
69 low = i - 1; 69 low = i - 1;
70 high = i; 70 high = i;
71 } 71 }
72 } 72 }
73 else { 73 else {
74 for (G4int jj = 0; jj < 10; jj++) { 74 for (G4int jj = 0; jj < 10; jj++) {
75 if (theData[i].GetX() < e) break; 75 if (theData[i].GetX() < e) break;
76 i--; 76 i--;
77 } 77 }
78 low = i; 78 low = i;
79 high = i + 1; 79 high = i + 1;
80 } 80 }
81 G4double x1, x2, y1, y2, x, y; 81 G4double x1, x2, y1, y2, x, y;
82 x = e; 82 x = e;
83 x1 = theData[low].GetX(); 83 x1 = theData[low].GetX();
84 x2 = theData[high].GetX(); 84 x2 = theData[high].GetX();
85 y1 = theData[low].GetY(j); 85 y1 = theData[low].GetY(j);
86 y2 = theData[high].GetY(j); 86 y2 = theData[high].GetY(j);
87 y = x * (y2 - y1) / (x2 - x1); 87 y = x * (y2 - y1) / (x2 - x1);
88 return y += y2 - x2 * (y2 - y1) / (x2 - x1); 88 return y += y2 - x2 * (y2 - y1) / (x2 - x1);
89 } 89 }
90 90
91 void G4ParticleHPField::Dump() 91 void G4ParticleHPField::Dump()
92 { 92 {
93 G4cout << nEntries << G4endl; 93 G4cout << nEntries << G4endl;
94 for (G4int i = 0; i < nEntries; ++i) { << 94 for (G4int i = 0; i < nEntries; i++) {
95 G4cout << theData[i].GetX() << " "; 95 G4cout << theData[i].GetX() << " ";
96 for (G4int j = 0; j < theData[i].GetDepth( << 96 for (G4int j = 0; j < theData[i].GetDepth(); j++) {
97 G4cout << theData[i].GetY(j) << " "; 97 G4cout << theData[i].GetY(j) << " ";
98 } 98 }
99 G4cout << G4endl; 99 G4cout << G4endl;
100 } 100 }
101 } 101 }
102 102
103 void G4ParticleHPField::Check(G4int i) 103 void G4ParticleHPField::Check(G4int i)
104 { 104 {
105 if (i > nEntries) 105 if (i > nEntries)
106 throw G4HadronicException(__FILE__, __LINE 106 throw G4HadronicException(__FILE__, __LINE__,
107 "Skipped some in 107 "Skipped some index numbers in G4ParticleHPField");
108 if (i == nPoints) { 108 if (i == nPoints) {
109 nPoints += 50; 109 nPoints += 50;
110 const std::size_t fsize = nPoints > 0 ? nP << 110 auto buff = new G4ParticleHPFieldPoint[nPoints];
111 auto buff = new G4ParticleHPFieldPoint[fsi << 111 // G4cout << "copying 1"<<G4endl;
112 for (G4int j = 0; j < nEntries; ++j) { << 112 for (G4int j = 0; j < nEntries; j++) {
113 buff[j] = theData[j]; 113 buff[j] = theData[j];
114 } 114 }
>> 115 // G4cout << "copying 2"<<G4endl;
115 delete[] theData; 116 delete[] theData;
116 theData = buff; 117 theData = buff;
117 } 118 }
118 if (i == nEntries) nEntries = i + 1; 119 if (i == nEntries) nEntries = i + 1;
119 } 120 }
120 121