#include <Doxymodules_field.h>

Inheritance diagram for G4BlinePrimaryGeneratorAction:

Public Member Functions
	G4BlinePrimaryGeneratorAction ()=default

	~G4BlinePrimaryGeneratorAction () override=default

void	GeneratePrimaries (G4Event *anEvent) override

void	SetUserPrimaryAction (G4VUserPrimaryGeneratorAction *anAction)

Private Attributes
G4VUserPrimaryGeneratorAction *	fUserPrimaryAction = nullptr

G4bool	fFirstPartOfBline = true

G4ThreeVector	fBlineStartPosition

G4double	fT0 = 0.

Detailed Description

Definition at line 20 of file Doxymodules_field.h.

Constructor & Destructor Documentation

◆ G4BlinePrimaryGeneratorAction()

G4BlinePrimaryGeneratorAction::G4BlinePrimaryGeneratorAction ( )

default

◆ ~G4BlinePrimaryGeneratorAction()

G4BlinePrimaryGeneratorAction::~G4BlinePrimaryGeneratorAction ( )

overridedefault

Member Function Documentation

◆ GeneratePrimaries()

void G4BlinePrimaryGeneratorAction::GeneratePrimaries ( G4Event * anEvent )

override

Definition at line 50 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
 
    auto  tmpEvent = new G4Event();
    fUserPrimaryAction->GeneratePrimaries(tmpEvent);
    fBlineStartPosition = tmpEvent->GetPrimaryVertex()->GetPosition();
    fT0 = tmpEvent->GetPrimaryVertex()->GetT0();
    delete tmpEvent;
  }
  fFirstPartOfBline = false;
 
  auto  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;
 
  auto  particle = new G4PrimaryParticle(pdef,px,py,pz);
  particle->SetMass( mass );
  particle->SetCharge(pdef->GetPDGCharge());
  primary_vertex->SetPrimary( particle );
 
  anEvent->AddPrimaryVertex( primary_vertex );
}