#include <G4BlinePrimaryGeneratorAction.hh>

Inheritance diagram for G4BlinePrimaryGeneratorAction:

Collaboration diagram for G4BlinePrimaryGeneratorAction:

Public Member Functions
	G4BlinePrimaryGeneratorAction ()

virtual	~G4BlinePrimaryGeneratorAction ()

virtual void	GeneratePrimaries (G4Event *anEvent)

void	SetUserPrimaryAction (G4VUserPrimaryGeneratorAction *anAction)

Public Member Functions inherited from G4VUserPrimaryGeneratorAction
	G4VUserPrimaryGeneratorAction ()

virtual	~G4VUserPrimaryGeneratorAction ()

Private Attributes
G4VUserPrimaryGeneratorAction *	fUserPrimaryAction

G4bool	fFirstPartOfBline

G4ThreeVector	fBlineStartPosition

G4double	fT0

Detailed Description

Definition at line 60 of file G4BlinePrimaryGeneratorAction.hh.

Constructor & Destructor Documentation

◆ G4BlinePrimaryGeneratorAction()

G4BlinePrimaryGeneratorAction::G4BlinePrimaryGeneratorAction ( )

Definition at line 51 of file G4BlinePrimaryGeneratorAction.cc.

 {
   fUserPrimaryAction = 0;
   fFirstPartOfBline = true;
 }

◆ ~G4BlinePrimaryGeneratorAction()

G4BlinePrimaryGeneratorAction::~G4BlinePrimaryGeneratorAction ( )

virtual

Definition at line 59 of file G4BlinePrimaryGeneratorAction.cc.

60 {

61 }

Member Function Documentation

◆ GeneratePrimaries()

void G4BlinePrimaryGeneratorAction::GeneratePrimaries ( G4Event * anEvent )

virtual

Implements G4VUserPrimaryGeneratorAction.

Definition at line 65 of file G4BlinePrimaryGeneratorAction.cc.

 {
   if (!fUserPrimaryAction)
   {
     G4Exception("G4BlinePrimaryGeneratorAction::GeneratePrimaries()",
                 "NullPointer", JustWarning,
                 "Primary generator action not defined !");
     return;
   }
 
   // For the first part of a bline the start position and time are defined
   // by using the USER primary action while for the second part the previous
   // values are taken.
   
   if (fFirstPartOfBline)
   {
     // set the position and time defined by using the USER primary action
  
     G4Event* tmpEvent = new G4Event();    
     fUserPrimaryAction->GeneratePrimaries(tmpEvent);
     fBlineStartPosition = tmpEvent->GetPrimaryVertex()->GetPosition();
     fT0 = tmpEvent->GetPrimaryVertex()->GetT0();
     delete tmpEvent;
   }
   fFirstPartOfBline = false;   
 
   G4PrimaryVertex* primary_vertex = 
     new G4PrimaryVertex(fBlineStartPosition, fT0);
 
   // Define the particle to be tracked as Charged Geantino
     
   G4ChargedGeantino* pdef = G4ChargedGeantino::ChargedGeantino();
    
   G4double mass =  pdef->GetPDGMass();
   G4double energy = 10000.*MeV + mass;
   G4double pmom = std::sqrt(energy*energy-mass*mass);
 
   // The momentum direction and energy do not have an effect in tracing of 
   // bline but still need to be defined.
    
   G4double px = 0.;
   G4double py = 0.;
   G4double pz = pmom;
 
   G4PrimaryParticle* particle = new G4PrimaryParticle(pdef,px,py,pz);
   particle->SetMass( mass );
   particle->SetCharge(pdef->GetPDGCharge());
   primary_vertex->SetPrimary( particle );
  
   anEvent->AddPrimaryVertex( primary_vertex );
 }