39 #include "PrimaryGeneratorAction.hh"
119 boundYHalfLength*boundYHalfLength+
120 boundZHalfLength*boundZHalfLength);
122 G4double sinTheta = std::sqrt(1.-cosTheta*cosTheta);
125 G4ThreeVector positionStart(boundPos.
x()+radius*sinTheta*std::cos(phi),
126 boundPos.
y()+radius*sinTheta*std::sin(phi),
127 boundPos.
z()+radius*cosTheta);
140 (positionDir-positionStart).unit());