55 result->
SetData(running++, x, y);
59 else if( left.
GetX(i)+right.
GetX(j) == 0
60 || std::abs((right.
GetX(j)-left.
GetX(i))/(left.
GetX(i)+right.
GetX(j))) > 0.001 )
64 result->
SetData(running++, x, y);
76 result->
SetData(running++, x, y);
127 if(&right ==
this)
return *
this;
240 G4int s_tmp = 0,
n=0, m_tmp=0;
257 && std::abs(std::abs(xp-xa)/xa) < 0.0000001
262 active = passive; a=
p;
267 active = passive; a=
p;
276 anX = passive->
GetXsec(p)-deltaX;
303 G4int count = 0, current = 2, start = 1;
311 x1=aBuff[count].
GetX();
312 y1=aBuff[count].
GetY();
318 for (
G4int j=start; j<current; j++ ) {
321 aBuff[++count] =
theData[current-1];
327 for(
G4int j=start; j<current; j++)
330 if(x1-x2 == 0) y = (y2+
y1)/2.;
334 aBuff[++count] =
theData[current-1];
358 std::vector<G4double>::iterator i;
362 if(std::abs(aX-aBlock) < 0.1*
MeV)
395 G4int icounter_max=1024;
399 if ( icounter > icounter_max ) {
400 G4cout <<
"Loop-counter exceeded the threshold value at " << __LINE__ <<
"th line of " << __FILE__ <<
"." <<
G4endl;
422 G4int jcounter_max=1024;
426 if ( jcounter > jcounter_max ) {
427 G4cout <<
"Loop-counter exceeded the threshold value at " << __LINE__ <<
"th line of " << __FILE__ <<
"." <<
G4endl;
440 if ( ibin < 0 )
G4cout <<
"TKDB 080807 " << rand <<
G4endl;
456 value = rand * ( x2 -
x1 ) + x1;
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4ParticleHPInterpolator theLin
G4double Get15percentBorder()
std::vector< ExP01TrackerHit * > a
G4bool IsBlocked(G4double aX)
static constexpr double MeV
Float_t y1[n_points_granero]
G4InterpolationScheme GetScheme(G4int anIndex)
G4double Interpolate(G4InterpolationScheme aScheme, G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const
G4double GetX(G4int i) const
Float_t x1[n_points_granero]
std::vector< G4double > theBlocked
G4double Get50percentBorder()
G4int GetMinIndex(G4double e) const
G4InterpolationScheme GetScheme(G4int index) const
G4InterpolationManager theManager
G4double GetXsec(G4int i)
Float_t y2[n_points_geant4]
G4double GetY(G4double x)
const XML_Char int const XML_Char * value
void ThinOut(G4double precision)
BasicVector3D< float > operator+(const BasicVector3D< float > &v)
void SetData(G4int i, G4double x, G4double y)
G4ParticleHPVector & operator=(const G4ParticleHPVector &right)
G4double G4ParticleHPJENDLHEData::G4double result
void SetPoint(G4int i, const G4ParticleHPDataPoint &it)
const G4ParticleHPDataPoint & GetPoint(G4int i) const
void SetY(G4int i, G4double x)
G4double the50percentBorderCash
G4double Lin(G4double x, G4double x1, G4double x2, G4double y1, G4double y2)
G4GLOB_DLL std::ostream G4cout
G4int GetVectorLength() const
G4ParticleHPDataPoint * theData
G4double the15percentBorderCash
void AppendScheme(G4int aPoint, const G4InterpolationScheme &aScheme)
Float_t x2[n_points_geant4]
std::vector< G4double > theBuffered
G4ParticleHPInterpolator theInt
void Merge(G4ParticleHPVector *active, G4ParticleHPVector *passive)
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
G4double GetEnergy(G4int i) const
void IntegrateAndNormalise()