Geant4  10.02.p01
G4ParticleHPField.cc
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 // neutron_hp -- source file
27 // J.P. Wellisch, Nov-1996
28 // A prototype of the low energy neutron transport model.
29 //
30 // P. Arce, June-2014 Conversion neutron_hp to particle_hp
31 //
32 #include "G4ParticleHPField.hh"
33 #include "G4HadronicException.hh"
34 #include "G4ios.hh"
35 
36 
38  {
39  theData = new G4ParticleHPFieldPoint[100];
40  nPoints=100;
41  nEntries=0;
42  theData->SetData(0,0,0);
43  }
44 
46 
48  {
49  G4int low = 0;
50  G4int high = 0;
51  G4int i;
52  for (i=1; i<nEntries/10; i++)
53  {
54  if(theData[10*i].GetX()>e) break;
55  }
56  if(i==(nEntries/10))
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);
88  }
89 
91  {
93  for(G4int i=0; i<nEntries; i++)
94  {
95  G4cout << theData[i].GetX()<<" ";
96  for(G4int j=0; j<theData[i].GetDepth(); j++)
97  {
98  G4cout << theData[i].GetY(j)<<" ";
99  }
100  G4cout << G4endl;
101  }
102  }
103 
105  {
106  if(i>nEntries) throw G4HadronicException(__FILE__, __LINE__, "Skipped some index numbers in G4ParticleHPField");
107  if(i==nPoints)
108  {
109  nPoints += 50;
111 // G4cout << "copying 1"<<G4endl;
112  for (G4int j=0; j<nEntries; j++)
113  {
114  buff[j] = theData[j];
115  }
116 // G4cout << "copying 2"<<G4endl;
117  delete [] theData;
118  theData = buff;
119  }
120  if(i==nEntries) nEntries=i+1;
121  }
G4double GetY(G4int i, G4int j)
G4double GetX(G4int i)
int G4int
Definition: G4Types.hh:78
G4GLOB_DLL std::ostream G4cout
void SetData(G4double e, G4int i, G4double x)
const G4double x[NPOINTSGL]
#define G4endl
Definition: G4ios.hh:61
G4ParticleHPFieldPoint * theData
double G4double
Definition: G4Types.hh:76
G4double GetY(G4int i) const