#include <Doxymodules_optical.h>

Inheritance diagram for OpNovicePrimaryGeneratorAction:

Public Member Functions
	OpNovicePrimaryGeneratorAction ()

	~OpNovicePrimaryGeneratorAction () override

void	GeneratePrimaries (G4Event *) override

void	SetOptPhotonPolar ()

void	SetOptPhotonPolar (G4double)

G4ParticleGun *	GetParticleGun ()

Private Attributes
G4ParticleGun *	fParticleGun = nullptr

OpNovicePrimaryGeneratorMessenger *	fGunMessenger = nullptr

Detailed Description

Definition at line 22 of file Doxymodules_optical.h.

Constructor & Destructor Documentation

◆ OpNovicePrimaryGeneratorAction()

OpNovicePrimaryGeneratorAction::OpNovicePrimaryGeneratorAction ( )

Definition at line 44 of file OpNovicePrimaryGeneratorAction.cc.

  : G4VUserPrimaryGeneratorAction()
  , fParticleGun(nullptr)
{
  G4int n_particle = 1;
  fParticleGun     = new G4ParticleGun(n_particle);
  // create a messenger for this class
  fGunMessenger = new OpNovicePrimaryGeneratorMessenger(this);
  // default kinematic
  //
  G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
  G4ParticleDefinition* particle = particleTable->FindParticle("e+");
 
  fParticleGun->SetParticleDefinition(particle);
  fParticleGun->SetParticleTime(0.0 * ns);
  fParticleGun->SetParticlePosition(G4ThreeVector(0., 0., 0.));
  fParticleGun->SetParticleMomentumDirection(G4ThreeVector(1., 0., 0.));
  fParticleGun->SetParticleEnergy(500.0 * keV);
}

◆ ~OpNovicePrimaryGeneratorAction()

OpNovicePrimaryGeneratorAction::~OpNovicePrimaryGeneratorAction ( )

override

Definition at line 65 of file OpNovicePrimaryGeneratorAction.cc.

{
  delete fParticleGun;
  delete fGunMessenger;
}

Member Function Documentation

◆ GeneratePrimaries()

void OpNovicePrimaryGeneratorAction::GeneratePrimaries ( G4Event * anEvent )

override

Definition at line 72 of file OpNovicePrimaryGeneratorAction.cc.

{
  fParticleGun->GeneratePrimaryVertex(anEvent);
}

◆ SetOptPhotonPolar() [1/2]

void OpNovicePrimaryGeneratorAction::SetOptPhotonPolar ( )

Definition at line 78 of file OpNovicePrimaryGeneratorAction.cc.

{
  G4double angle = G4UniformRand() * 360.0 * deg;
  SetOptPhotonPolar(angle);
}

◆ SetOptPhotonPolar() [2/2]

void OpNovicePrimaryGeneratorAction::SetOptPhotonPolar ( G4double angle )

Definition at line 85 of file OpNovicePrimaryGeneratorAction.cc.

{
  if(fParticleGun->GetParticleDefinition()->GetParticleName() !=
     "opticalphoton")
  {
    G4ExceptionDescription ed;
    ed << "Warning: the particleGun is not an opticalphoton";
    G4Exception("OpNovicePrimaryGeneratorAction::SetOptPhotonPolar()",
                "OpNovice_010", JustWarning, ed);
    return;
  }
 
  G4ThreeVector normal(1., 0., 0.);
  G4ThreeVector kphoton = fParticleGun->GetParticleMomentumDirection();
  G4ThreeVector product = normal.cross(kphoton);
  G4double modul2       = product * product;
 
  G4ThreeVector e_perpend(0., 0., 1.);
  if(modul2 > 0.)
    e_perpend = (1. / std::sqrt(modul2)) * product;
  G4ThreeVector e_paralle = e_perpend.cross(kphoton);
 
  G4ThreeVector polar =
    std::cos(angle) * e_paralle + std::sin(angle) * e_perpend;
  fParticleGun->SetParticlePolarization(polar);
}