54 result->
SetData(running++, x, y);
58 else if( left.
GetX(i)+right.
GetX(j) == 0
59 || std::abs((right.
GetX(j)-left.
GetX(i))/(left.
GetX(i)+right.
GetX(j))) > 0.001 )
63 result->
SetData(running++, x, y);
75 result->
SetData(running++, x, y);
126 if(&right ==
this)
return *
this;
232 G4int s_tmp = 0,
n=0, m_tmp=0;
249 && std::abs(std::abs(xp-xa)/xa) < 0.0000001
254 active = passive; a=p;
259 active = passive; a=p;
268 anX = passive->
GetXsec(p)-deltaX;
295 G4int count = 0, current = 2, start = 1;
303 x1=aBuff[count].
GetX();
304 y1=aBuff[count].
GetY();
307 for(
G4int j=start; j<current; j++)
310 if(x1-x2 == 0) y = (y2+y1)/2.;
314 aBuff[++count] =
theData[current-1];
337 std::vector<G4double>::iterator i;
341 if(std::abs(aX-aBlock) < 0.1*
MeV)
405 if ( ibin < 0 )
G4cout <<
"TKDB 080807 " << rand <<
G4endl;
421 value = rand * ( x2 - x1 ) + x1;
G4double GetEnergy(G4int i) const
G4double Get15percentBorder()
G4int GetVectorLength() const
G4double the15percentBorderCash
G4double Lin(G4double x, G4double x1, G4double x2, G4double y1, G4double y2)
void SetData(G4int i, G4double x, G4double y)
G4ParticleHPInterpolator theInt
G4ParticleHPVector & operator=(const G4ParticleHPVector &right)
G4double GetXsec(G4int i)
G4InterpolationManager theManager
void SetY(G4int i, G4double x)
void AppendScheme(G4int aPoint, const G4InterpolationScheme &aScheme)
G4ParticleHPVector & operator+(G4ParticleHPVector &left, G4ParticleHPVector &right)
G4GLOB_DLL std::ostream G4cout
G4ParticleHPDataPoint * theData
G4InterpolationScheme GetScheme(G4int index) const
void Merge(G4ParticleHPVector *active, G4ParticleHPVector *passive)
G4double GetX(G4int i) const
void SetPoint(G4int i, const G4ParticleHPDataPoint &it)
G4double GetY(G4double x)
void IntegrateAndNormalise()
G4double Interpolate(G4InterpolationScheme aScheme, G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const
G4bool IsBlocked(G4double aX)
G4double the50percentBorderCash
T max(const T t1, const T t2)
brief Return the largest of the two arguments
std::vector< G4double > theBlocked
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
G4ParticleHPInterpolator theLin
std::vector< G4double > theBuffered
G4int GetMinIndex(G4double e) const
G4InterpolationScheme GetScheme(G4int anIndex)
const G4ParticleHPDataPoint & GetPoint(G4int i) const
void ThinOut(G4double precision)
G4double Get50percentBorder()