G4ParticleHPVector.cc File Reference

#include "G4ParticleHPVector.hh"
#include "G4SystemOfUnits.hh"
#include "G4Threading.hh"

Go to the source code of this file.

Functions
G4ParticleHPVector &	operator+ (G4ParticleHPVector &left, G4ParticleHPVector &right)

Function Documentation

operator+()

G4ParticleHPVector & operator+	(	G4ParticleHPVector &	left,
		G4ParticleHPVector &	right )

Definition at line 41 of file G4ParticleHPVector.cc.

{
  auto result = new G4ParticleHPVector;
  G4int j = 0;
  G4double x;
  G4double y;
  G4int running = 0;
  for (G4int i = 0; i < left.GetVectorLength(); i++) {
    while (j < right.GetVectorLength())  // Loop checking, 11.05.2015, T. Koi
    {
      if (right.GetX(j) < left.GetX(i) * 1.001) {
        x = right.GetX(j);
        y = right.GetY(j) + left.GetY(x);
        result->SetData(running++, x, y);
        j++;
      }
      // else if(std::abs((right.GetX(j)-left.GetX(i))/(left.GetX(i)+right.GetX(j)))>0.001)
      else if (left.GetX(i) + right.GetX(j) == 0
               || std::abs((right.GetX(j) - left.GetX(i)) / (left.GetX(i) + right.GetX(j))) > 0.001)
      {
        x = left.GetX(i);
        y = left.GetY(i) + right.GetY(x);
        result->SetData(running++, x, y);
        break;
      }
      else {
        break;
      }
    }
    if (j == right.GetVectorLength()) {
      x = left.GetX(i);
      y = left.GetY(i) + right.GetY(x);
      result->SetData(running++, x, y);
    }
  }
  result->ThinOut(0.02);
  return *result;
}