PrimaryGenerator Class Reference

#include <PrimaryGenerator.hh>

Inheritance diagram for PrimaryGenerator:

Collaboration diagram for PrimaryGenerator:

Public Member Functions
	PrimaryGenerator ()
	~PrimaryGenerator ()
virtual void	GeneratePrimaryVertex (G4Event *)
Public Member Functions inherited from G4VPrimaryGenerator
	G4VPrimaryGenerator ()
virtual	~G4VPrimaryGenerator ()
G4ThreeVector	GetParticlePosition ()
G4double	GetParticleTime ()
void	SetParticlePosition (G4ThreeVector aPosition)
void	SetParticleTime (G4double aTime)

Additional Inherited Members
Static Public Member Functions inherited from G4VPrimaryGenerator
static G4bool	CheckVertexInsideWorld (const G4ThreeVector &pos)
Protected Attributes inherited from G4VPrimaryGenerator
G4ThreeVector	particle_position
G4double	particle_time

Detailed Description

Definition at line 44 of file PrimaryGenerator.hh.

Constructor & Destructor Documentation

PrimaryGenerator()

PrimaryGenerator::PrimaryGenerator

(

)

Definition at line 48 of file PrimaryGenerator.cc.

: G4VPrimaryGenerator()
{ }

~PrimaryGenerator()

PrimaryGenerator::~PrimaryGenerator

(

)

Definition at line 54 of file PrimaryGenerator.cc.

{ }

Member Function Documentation

GeneratePrimaryVertex()

void PrimaryGenerator::GeneratePrimaryVertex ( G4Event *  event )

virtual

Implements G4VPrimaryGenerator.

Definition at line 59 of file PrimaryGenerator.cc.

{
  //vertex A uniform on a cylinder
  //
  const G4double r = 2*mm;
  const G4double zmax = 8*mm;
  //
  G4double alpha = twopi*G4UniformRand();     //alpha uniform in (0, 2*pi)
  G4double ux = std::cos(alpha);
  G4double uy = std::sin(alpha);
  G4double z = zmax*(2*G4UniformRand() - 1);  //z uniform in (-zmax, +zmax)
  G4ThreeVector positionA(r*ux,r*uy,z);
  G4double timeA = 0*s;
  // 
  G4PrimaryVertex* vertexA = new G4PrimaryVertex(positionA, timeA);
  
  //particle 1 at vertex A
  //
  G4ParticleDefinition* particleDefinition
           = G4ParticleTable::GetParticleTable()->FindParticle("geantino");
  G4PrimaryParticle* particle1 = new G4PrimaryParticle(particleDefinition);
  particle1->SetMomentumDirection(G4ThreeVector(ux,uy,0));    
  particle1->SetKineticEnergy(1*MeV);
  //
  vertexA->SetPrimary(particle1);
  event->AddPrimaryVertex(vertexA);

  //vertex (B) symetric to vertex A
  //
  alpha += pi;
  ux = std::cos(alpha);
  uy = std::sin(alpha);
  G4ThreeVector positionB(r*ux,r*uy,z);
  G4double timeB = 1*s;
  // 
  G4PrimaryVertex* vertexB = new G4PrimaryVertex(positionB, timeB);
  
  //particles 2 and 3 at vertex B
  //
  const G4double dalpha = 10*deg;
  ux = std::cos(alpha + dalpha);
  uy = std::sin(alpha + dalpha);        
  G4PrimaryParticle* particle2 = new G4PrimaryParticle(particleDefinition);
  particle2->SetMomentumDirection(G4ThreeVector(ux,uy,0));    
  particle2->SetKineticEnergy(1*keV);
  //
  ux = std::cos(alpha - dalpha);
  uy = std::sin(alpha - dalpha);        
  G4PrimaryParticle* particle3 = new G4PrimaryParticle(particleDefinition);
  particle3->SetMomentumDirection(G4ThreeVector(ux,uy,0));    
  particle3->SetKineticEnergy(1*GeV);
  //
  vertexB->SetPrimary(particle2);
  vertexB->SetPrimary(particle3);  
  event->AddPrimaryVertex(vertexB);
}

Here is the call graph for this function:

---

The documentation for this class was generated from the following files:

• PrimaryGenerator.hh
• PrimaryGenerator.cc