The primary generator action class with particle gun. More...

#include <Doxymodules_medical.h>

Inheritance diagram for Dicom2PrimaryGeneratorAction:

Public Member Functions
	Dicom2PrimaryGeneratorAction ()

virtual	~Dicom2PrimaryGeneratorAction ()

virtual void	GeneratePrimaries (G4Event *)

const G4ParticleGun *	GetParticleGun () const

Private Attributes
G4ParticleGun *	fParticleGun

G4Box *	fEnvelopeBox

G4VPhysicalVolume *	fEnvelopeVol

G4ThreeVector	fPosCenter

G4ThreeVector	fPosDelta

G4double	fGeomFactor

Detailed Description

The primary generator action class with particle gun.

The default kinematic is a 6 MeV gamma, randomly distribued in front of the phantom across 80% of the (X,Y) phantom size.

Definition at line 71 of file Doxymodules_medical.h.

Constructor & Destructor Documentation

◆ Dicom2PrimaryGeneratorAction()

Dicom2PrimaryGeneratorAction::Dicom2PrimaryGeneratorAction ( )

Definition at line 49 of file Dicom2PrimaryGeneratorAction.cc.

: G4VUserPrimaryGeneratorAction(),
  fParticleGun(new G4ParticleGun(1)),
  fEnvelopeBox(nullptr),
  fEnvelopeVol(nullptr),
  fPosCenter(G4ThreeVector(0.0)),
  fPosDelta(G4ThreeVector(0.0)),
  fGeomFactor(0.8)
{
    // default particle kinematic
    G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
    G4String particleName;
    G4ParticleDefinition* particle
            = particleTable->FindParticle(particleName="e-");
    fParticleGun->SetParticleDefinition(particle);
    fParticleGun->SetParticleMomentumDirection(G4ThreeVector(0.,0.,1.));
    fParticleGun->SetParticleEnergy(10.*MeV);
}

◆ ~Dicom2PrimaryGeneratorAction()

Dicom2PrimaryGeneratorAction::~Dicom2PrimaryGeneratorAction ( )

virtual

Definition at line 70 of file Dicom2PrimaryGeneratorAction.cc.

{
    delete fParticleGun;
}

Member Function Documentation

◆ GeneratePrimaries()

void Dicom2PrimaryGeneratorAction::GeneratePrimaries ( G4Event * anEvent )

virtual

Definition at line 77 of file Dicom2PrimaryGeneratorAction.cc.

{
    //this function is called at the begining of ecah event
    //
 
    // In order to avoid dependence of PrimaryGeneratorAction
    // on DetectorConstruction class we get phantomContainer volume
    // from G4LogicalVolumeStore.
    if (!fEnvelopeBox)
    {
        G4LogicalVolume* envLV
                = G4LogicalVolumeStore::GetInstance()->GetVolume("phantomContainer");
        if(envLV)
            fEnvelopeBox = dynamic_cast<G4Box*>(envLV->GetSolid());
    }
 
    if(!fEnvelopeVol)
    {
        fEnvelopeVol = G4PhysicalVolumeStore::GetInstance()->GetVolume("phantomContainer");
    }
 
    // get the center of the phantom
    if(fEnvelopeVol)
    {
        fPosCenter = fEnvelopeVol->GetObjectTranslation();
    }
    else
    {
        G4ExceptionDescription msg;
        msg << "Envelope physical volume not found.\n";
        msg << "The gun will be place at the center.";
        G4Exception("Dicom2PrimaryGeneratorAction::GeneratePrimaries()",
                    "DICOM20002",JustWarning, msg);
    }
 
    // get the size of the phantom
    if(fEnvelopeBox)
    {
        fPosDelta = G4ThreeVector(2.0 * fEnvelopeBox->GetXHalfLength(),
                                  2.0 * fEnvelopeBox->GetYHalfLength(),
                                  2.0 * fEnvelopeBox->GetZHalfLength());
    }
    else
    {
        G4ExceptionDescription msg;
        msg << "Envelope volume of box shape not found.\n";
        msg << "The gun will be place at the center.";
        G4Exception("Dicom2PrimaryGeneratorAction::GeneratePrimaries()",
                    "DICOM20003",JustWarning, msg);
    }
 
    // start at center
    G4ThreeVector pos = fPosCenter;
    // get a random displacement
    G4ThreeVector delta(fPosDelta.x() * (G4UniformRand()-0.5),
                        fPosDelta.y() * (G4UniformRand()-0.5),
                        fPosDelta.z() * (G4UniformRand()-0.5));
    // scale by geom factor
    delta *= fGeomFactor;
 
    //G4cout << "Starting position: " << G4BestUnit(pos + delta, "Length")
    //       << ", pos = " << G4BestUnit(fPosCenter, "Length")
    //       << ", size = " << G4BestUnit(fPosDelta, "Length")
    //       << ", delta = " << G4BestUnit(delta, "Length")
    //       << G4endl;
 
    fParticleGun->SetParticlePosition(pos + delta);
    fParticleGun->SetParticleMomentumDirection(G4RandomDirection());
    fParticleGun->GeneratePrimaryVertex(anEvent);
}