40 #include "PrimaryGeneratorAction.hh"
62 G4double x0,y0,
z0,theta,phi,xMom0,yMom0,zMom0,e0;
71 while (! (std::sqrt(x0*x0+y0*y0)<= sizeMax) )
86 while(theta>1.1e-6*
rad);
90 xMom0=std::sin(theta)*std::cos(phi);
91 yMom0=std::sin(theta)*std::sin(phi);
92 zMom0=std::cos(theta);
97 <<
"-> Event # " << numEvent
ThreeVector shoot(const G4int Ap, const G4int Af)
G4ParticleDefinition * FindParticle(G4int PDGEncoding)
static constexpr double mm
CLHEP::Hep3Vector G4ThreeVector
static constexpr double rad
void SetParticleMomentumDirection(G4ParticleMomentum aMomentumDirection)
virtual void GeneratePrimaryVertex(G4Event *evt)
void SetParticlePosition(G4ThreeVector aPosition)
G4GLOB_DLL std::ostream G4cout
void GeneratePrimaries(G4Event *)
void SetParticleEnergy(G4double aKineticEnergy)
static G4ParticleTable * GetParticleTable()
static constexpr double MeV
~PrimaryGeneratorAction()
static constexpr double micrometer
static constexpr double twopi
void SetParticleDefinition(G4ParticleDefinition *aParticleDefinition)