Geant4-11
G4ParticleHPPartial.hh
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//
27// 070618 Comment out unused private member leaking by T. Koi
28//
29// P. Arce, June-2014 Conversion neutron_hp to particle_hp
30//
31#ifndef G4ParticleHPPartial_h
32#define G4ParticleHPPartial_h 1
33
35
36#include "globals.hh"
37#include "G4ParticleHPVector.hh"
40
42{
43 public:
44
46 {
47 X = new G4double[n];
49 nData = n;
50 T=0;
51 }
52
54 {
55 T = new G4double[n2];
56 X = new G4double[n1];
57 data = new G4ParticleHPVector[n1];
58 nData = std::max(n1,n2);
59 }
60
61 void InitInterpolation(G4int i, std::istream & aDataFile)
62 {
63 data[i].InitInterpolation(aDataFile);
64 }
65
66 void InitInterpolation(std::istream & aDataFile)
67 {
68 theManager.Init(aDataFile);
69 }
70
71 void Init(std::istream & aDataFile)
72 {
73 G4int i;
74 G4double e;
75 for( i=0; i<nData; i++)
76 {
77 aDataFile >> e;
78 e *= CLHEP::eV;
79 SetX(i,e);
80 InitData(i, aDataFile, CLHEP::eV); // energy and probability for gammas
81 }
82 }
83
84 void InitData(G4int i, std::istream & aDataFile, G4double unit=1.)
85 {
86 G4int ii;
87 G4double eg, pg;
88 G4int neg;
89 aDataFile >> neg;
90 data[i].InitInterpolation(aDataFile);
91 for (ii=0; ii<neg; ii++)
92 {
93 aDataFile >> eg >> pg;
94 eg *= unit;
95 SetX(i,ii,eg);
96 SetY(i,ii,pg);
97 }
98 data[i].Hash();
99 }
100
102 {
103 delete [] X;
104 if(T!=0) delete [] T;
105 delete [] data;
106 }
107 inline G4int GetNumberOfEnergies() {return nData;}
108
109 inline void SetX(G4int i, G4double x) {X[i]=x;}
110 inline void SetT(G4int i, G4double x) {T[i]=x;}
111 inline void SetX(G4int i, G4int j, G4double x) {data[i].SetX(j,x);}
112 inline void SetY(G4int i, G4int j, G4double y) {data[i].SetY(j,y);}
113 void DoneSetXY( G4int i ) { data[i].Hash(); }
114
115 inline G4double GetX(G4int i) {return X[i];}
116 inline G4double GetT(G4int i) {return T[i];}
117 inline G4double GetX(G4int i, G4int j) {return data[i].GetX(j);}
118 inline G4double GetY(G4int i, G4int j) {return data[i].GetY(j);}
119 inline G4double GetY(G4int i, G4double e) {return data[i].GetY(e);}
120 inline G4int GetNEntries(G4int i) {return data[i].GetVectorLength();}
123
124 private:
125
129 // TKDB
130 //G4ParticleHPVector * theBuffer;
132 G4InterpolationManager theManager; // interpolate between different data[i]
134};
135
136#endif
static const G4double e1[44]
double G4double
Definition: G4Types.hh:83
int G4int
Definition: G4Types.hh:85
void Init(G4int aScheme, G4int aRange)
void Init(std::istream &aDataFile)
G4ParticleHPPartial(G4int n1, G4int n2)
void InitInterpolation(std::istream &aDataFile)
G4double Sample(G4double x)
void SetT(G4int i, G4double x)
G4double GetY(G4int i, G4int j)
G4double GetX(G4int i, G4int j)
G4ParticleHPInterpolator theInt
G4double GetT(G4int i)
G4InterpolationManager theManager
void SetX(G4int i, G4int j, G4double x)
void SetX(G4int i, G4double x)
void SetY(G4int i, G4int j, G4double y)
G4double GetY(G4int i, G4double e)
void InitData(G4int i, std::istream &aDataFile, G4double unit=1.)
G4ParticleHPVector * data
void InitInterpolation(G4int i, std::istream &aDataFile)
G4double GetX(G4int i)
void SetY(G4int i, G4double x)
void SetX(G4int i, G4double e)
G4double GetY(G4double x)
G4double GetX(G4int i) const
G4int GetVectorLength() const
void InitInterpolation(std::istream &aDataFile)
static constexpr double eV
T max(const T t1, const T t2)
brief Return the largest of the two arguments
pg
Definition: demo.py:37