52 result->
SetData(running++, x, y);
56 else if( left.
GetX(i)+right.
GetX(j) == 0
57 || std::abs((right.
GetX(j)-left.
GetX(i))/(left.
GetX(i)+right.
GetX(j))) > 0.001 )
61 result->
SetData(running++, x, y);
73 result->
SetData(running++, x, y);
124 if(&right ==
this)
return *
this;
230 G4int s_tmp = 0,
n=0, m_tmp=0;
247 && std::abs(std::abs(xp-xa)/xa) < 0.0000001
252 active = passive; a=p;
257 active = passive; a=p;
266 anX = passive->
GetXsec(p)-deltaX;
293 G4int count = 0, current = 2, start = 1;
301 x1=aBuff[count].
GetX();
302 y1=aBuff[count].
GetY();
305 for(
G4int j=start; j<current; j++)
308 if(x1-x2 == 0) y = (y2+y1)/2.;
312 aBuff[++count] =
theData[current-1];
335 std::vector<G4double>::iterator i;
339 if(std::abs(aX-aBlock) < 0.1*
MeV)
403 if ( ibin < 0 )
G4cout <<
"TKDB 080807 " << rand <<
G4endl;
419 value = rand * ( x2 - x1 ) + x1;
G4double the50percentBorderCash
G4double GetEnergy(G4int i) const
G4double Get50percentBorder()
G4double GetY(G4double x)
G4int GetVectorLength() const
G4double the15percentBorderCash
void SetPoint(G4int i, const G4NeutronHPDataPoint &it)
G4int GetMinIndex(G4double e) const
void ThinOut(G4double precision)
void Merge(G4NeutronHPVector *active, G4NeutronHPVector *passive)
std::vector< G4double > theBlocked
G4double GetX(G4int i) const
void SetData(G4int i, G4double x, G4double y)
G4NeutronHPInterpolator theInt
G4NeutronHPVector & operator+(G4NeutronHPVector &left, G4NeutronHPVector &right)
G4double Interpolate(G4InterpolationScheme aScheme, G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const
G4bool IsBlocked(G4double aX)
G4InterpolationScheme GetScheme(G4int anIndex)
void AppendScheme(G4int aPoint, const G4InterpolationScheme &aScheme)
void SetY(G4int i, G4double x)
G4GLOB_DLL std::ostream G4cout
G4double Get15percentBorder()
G4InterpolationScheme GetScheme(G4int index) const
G4InterpolationManager theManager
G4double GetXsec(G4int i)
const G4NeutronHPDataPoint & GetPoint(G4int i) const
T max(const T t1, const T t2)
brief Return the largest of the two arguments
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
void IntegrateAndNormalise()
G4NeutronHPInterpolator theLin
G4double Lin(G4double x, G4double x1, G4double x2, G4double y1, G4double y2)
std::vector< G4double > theBuffered
G4NeutronHPVector & operator=(const G4NeutronHPVector &right)
G4NeutronHPDataPoint * theData