| 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 <<
34 #include "G4HadronicException.hh" 33 #include "G4HadronicException.hh"
35 #include "G4ios.hh" 34 #include "G4ios.hh"
36 35
37 G4ParticleHPField::G4ParticleHPField() <<
38 { <<
39 theData = new G4ParticleHPFieldPoint[100]; <<
40 nPoints = 100; <<
41 nEntries = 0; <<
42 theData->SetData(0, 0, 0); <<
43 } <<
44 <<
45 G4ParticleHPField::~G4ParticleHPField() <<
46 { <<
47 delete[] theData; <<
48 } <<
49 36
50 G4double G4ParticleHPField::GetY(G4double e, G << 37 G4ParticleHPField::G4ParticleHPField()
51 { << 38 {
52 G4int low = 0; << 39 theData = new G4ParticleHPFieldPoint[100];
53 G4int high = 0; << 40 nPoints=100;
54 G4int i; << 41 nEntries=0;
55 for (i = 1; i < nEntries / 10; i++) { << 42 theData->SetData(0,0,0);
56 if (theData[10 * i].GetX() > e) break; <<
57 } 43 }
58 if (i == (nEntries / 10)) { << 44
59 i = 10 * i; << 45 G4ParticleHPField::~G4ParticleHPField(){ delete [] theData;}
60 while (i < nEntries) // Loop checking, 11 << 46
>> 47 G4double G4ParticleHPField::GetY(G4double e, G4int j)
>> 48 {
>> 49 G4int low = 0;
>> 50 G4int high = 0;
>> 51 G4int i;
>> 52 for (i=1; i<nEntries/10; i++)
61 { 53 {
62 if (theData[i++].GetX() > e) break; << 54 if(theData[10*i].GetX()>e) break;
63 } <<
64 if (i == nEntries) { <<
65 low = nEntries - 1; <<
66 high = nEntries - 2; <<
67 } <<
68 else { <<
69 low = i - 1; <<
70 high = i; <<
71 } <<
72 } <<
73 else { <<
74 for (G4int jj = 0; jj < 10; jj++) { <<
75 if (theData[i].GetX() < e) break; <<
76 i--; <<
77 } 55 }
78 low = i; << 56 if(i==(nEntries/10))
79 high = i + 1; << 57 {
>> 58 i=10*i;
>> 59 while (i<nEntries) // Loop checking, 11.05.2015, T. Koi
>> 60 {
>> 61 if(theData[i++].GetX()>e) break;
>> 62 }
>> 63 if (i==nEntries)
>> 64 {
>> 65 low = nEntries-1;
>> 66 high = nEntries-2;
>> 67 }else{
>> 68 low = i-1;
>> 69 high = i;
>> 70 }
>> 71 }else{
>> 72 for (G4int jj=0; jj<10; jj++)
>> 73 {
>> 74 if(theData[i].GetX()<e) break;
>> 75 i--;
>> 76 }
>> 77 low = i;
>> 78 high = i+1;
>> 79 }
>> 80 G4double x1, x2, y1, y2, x, y;
>> 81 x = e;
>> 82 x1 = theData[low] .GetX();
>> 83 x2 = theData[high].GetX();
>> 84 y1 = theData[low] .GetY(j);
>> 85 y2 = theData[high].GetY(j);
>> 86 y = x*(y2-y1)/(x2-x1);
>> 87 return y += y2-x2*(y2-y1)/(x2-x1);
80 } 88 }
81 G4double x1, x2, y1, y2, x, y; <<
82 x = e; <<
83 x1 = theData[low].GetX(); <<
84 x2 = theData[high].GetX(); <<
85 y1 = theData[low].GetY(j); <<
86 y2 = theData[high].GetY(j); <<
87 y = x * (y2 - y1) / (x2 - x1); <<
88 return y += y2 - x2 * (y2 - y1) / (x2 - x1); <<
89 } <<
90 89
91 void G4ParticleHPField::Dump() << 90 void G4ParticleHPField::Dump()
92 { << 91 {
93 G4cout << nEntries << G4endl; << 92 G4cout << nEntries<<G4endl;
94 for (G4int i = 0; i < nEntries; ++i) { << 93 for(G4int i=0; i<nEntries; i++)
95 G4cout << theData[i].GetX() << " "; << 94 {
96 for (G4int j = 0; j < theData[i].GetDepth( << 95 G4cout << theData[i].GetX()<<" ";
97 G4cout << theData[i].GetY(j) << " "; << 96 for(G4int j=0; j<theData[i].GetDepth(); j++)
>> 97 {
>> 98 G4cout << theData[i].GetY(j)<<" ";
>> 99 }
>> 100 G4cout << G4endl;
98 } 101 }
99 G4cout << G4endl; <<
100 } 102 }
101 } << 103
102 << 104 void G4ParticleHPField::Check(G4int i)
103 void G4ParticleHPField::Check(G4int i) << 105 {
104 { << 106 if(i>nEntries) throw G4HadronicException(__FILE__, __LINE__, "Skipped some index numbers in G4ParticleHPField");
105 if (i > nEntries) << 107 if(i==nPoints)
106 throw G4HadronicException(__FILE__, __LINE << 108 {
107 "Skipped some in << 109 nPoints += 50;
108 if (i == nPoints) { << 110 G4ParticleHPFieldPoint * buff = new G4ParticleHPFieldPoint[nPoints];
109 nPoints += 50; << 111 // G4cout << "copying 1"<<G4endl;
110 const std::size_t fsize = nPoints > 0 ? nP << 112 for (G4int j=0; j<nEntries; j++)
111 auto buff = new G4ParticleHPFieldPoint[fsi << 113 {
112 for (G4int j = 0; j < nEntries; ++j) { << 114 buff[j] = theData[j];
113 buff[j] = theData[j]; << 115 }
>> 116 // G4cout << "copying 2"<<G4endl;
>> 117 delete [] theData;
>> 118 theData = buff;
114 } 119 }
115 delete[] theData; << 120 if(i==nEntries) nEntries=i+1;
116 theData = buff; <<
117 } 121 }
118 if (i == nEntries) nEntries = i + 1; <<
119 } <<
120 122