WLSPrimaryGeneratorAction Class Reference

#include <Doxymodules_optical.h>

Inheritance diagram for WLSPrimaryGeneratorAction:

Public Member Functions
	WLSPrimaryGeneratorAction (WLSDetectorConstruction *)
	~WLSPrimaryGeneratorAction () override
void	GeneratePrimaries (G4Event *) override
void	BuildEmissionSpectrum ()
void	SetOptPhotonPolar (G4double)
void	SetDecayTimeConstant (G4double)
void	SetUseSampledEnergy (G4bool v)

Protected Attributes
G4PhysicsTable *	fIntegralTable = nullptr

Private Member Functions
void	SetOptPhotonPolar ()
void	SetOptPhotonTime ()

Private Attributes
WLSDetectorConstruction *	fDetector = nullptr
G4GeneralParticleSource *	fParticleGun = nullptr
WLSPrimaryGeneratorMessenger *	fGunMessenger = nullptr
G4double	fTimeConstant = 0.
G4bool	fUseSampledEnergy = false

Static Private Attributes
static G4bool	fFirst = false

Detailed Description

Definition at line 96 of file Doxymodules_optical.h.

Constructor & Destructor Documentation

WLSPrimaryGeneratorAction()

WLSPrimaryGeneratorAction::WLSPrimaryGeneratorAction

(

WLSDetectorConstruction *

dc

)

Definition at line 57 of file WLSPrimaryGeneratorAction.cc.

{
  fDetector      = dc;
 
  fParticleGun = new G4GeneralParticleSource();
  fGunMessenger = new WLSPrimaryGeneratorMessenger(this);
}

~WLSPrimaryGeneratorAction()

WLSPrimaryGeneratorAction::~WLSPrimaryGeneratorAction ( )

override

Definition at line 68 of file WLSPrimaryGeneratorAction.cc.

{
  delete fParticleGun;
  delete fGunMessenger;
  if(fIntegralTable)
  {
    fIntegralTable->clearAndDestroy();
    delete fIntegralTable;
  }
}

Member Function Documentation

GeneratePrimaries()

void WLSPrimaryGeneratorAction::GeneratePrimaries ( G4Event *  anEvent )

override

Definition at line 142 of file WLSPrimaryGeneratorAction.cc.

{
  if(!fFirst)
  {
    fFirst = true;
    BuildEmissionSpectrum();
  }
 
  if(fUseSampledEnergy)
  {
    const G4MaterialTable* theMaterialTable = G4Material::GetMaterialTable();
 
    G4double sampledEnergy = 3. * eV;
 
    for(size_t j = 0; j < theMaterialTable->size(); ++j)
    {
      G4Material* fMaterial = (*theMaterialTable)[j];
      if(fMaterial->GetName() == "PMMA")
      {
        auto WLSIntensity =
          fMaterial->GetMaterialPropertiesTable()->GetProperty("WLSCOMPONENT");
 
        if(WLSIntensity)
        {
          auto WLSIntegral = (G4PhysicsFreeVector*) ((*fIntegralTable)(
              fMaterial->GetIndex()));
 
          G4double CIImax   = WLSIntegral->GetMaxValue();
          G4double CIIvalue = G4UniformRand() * CIImax;
 
          sampledEnergy = WLSIntegral->GetEnergy(CIIvalue);
        }
      }
    }
 
    // this does not work.
    G4String cmd = "/gun/energy " + G4UIcommand::ConvertToString(sampledEnergy / eV) + " eV";
    G4UImanager::GetUIpointer()->ApplyCommand(cmd);
  }
 
  // The code behind this line is not thread safe because polarization
  // and time are randomly selected and GPS properties are global
 
  G4AutoLock l(&gen_mutex);
  if(fParticleGun->GetParticleDefinition() == G4OpticalPhoton::Definition())
  {
    SetOptPhotonPolar();
    SetOptPhotonTime();
  }
 
  fParticleGun->GeneratePrimaryVertex(anEvent);
}

BuildEmissionSpectrum()

void WLSPrimaryGeneratorAction::BuildEmissionSpectrum

(

)

Definition at line 88 of file WLSPrimaryGeneratorAction.cc.

{
  if(fIntegralTable)
    return;
 
  const G4MaterialTable* theMaterialTable = G4Material::GetMaterialTable();
  G4int numOfMaterials                    = G4Material::GetNumberOfMaterials();
 
  if(!fIntegralTable)
    fIntegralTable = new G4PhysicsTable(numOfMaterials);
 
  for(G4int i = 0; i < numOfMaterials; ++i)
  {
    auto vec = new G4PhysicsFreeVector();
 
    G4MaterialPropertiesTable* MPT =
      (*theMaterialTable)[i]->GetMaterialPropertiesTable();
 
    if(MPT)
    {
      G4MaterialPropertyVector* theWLSVector = MPT->GetProperty("WLSCOMPONENT");
 
      if(theWLSVector)
      {
        G4double currentIN = (*theWLSVector)[0];
        if(currentIN >= 0.0)
        {
          G4double currentPM  = theWLSVector->Energy(0);
          G4double currentCII = 0.0;
          vec->InsertValues(currentPM, currentCII);
          G4double prevPM  = currentPM;
          G4double prevCII = currentCII;
          G4double prevIN  = currentIN;
 
          for(size_t j = 1; j < theWLSVector->GetVectorLength(); ++j)
          {
            currentPM  = theWLSVector->Energy(j);
            currentIN  = (*theWLSVector)[j];
            currentCII = 0.5 * (prevIN + currentIN);
            currentCII = prevCII + (currentPM - prevPM) * currentCII;
            vec->InsertValues(currentPM, currentCII);
            prevPM  = currentPM;
            prevCII = currentCII;
            prevIN  = currentIN;
          }
        }
      }
    }
    fIntegralTable->insertAt(i, vec);
  }
}

SetOptPhotonPolar() [1/2]

void WLSPrimaryGeneratorAction::SetOptPhotonPolar

(

G4double

angle

)

Definition at line 205 of file WLSPrimaryGeneratorAction.cc.

{
  if(fParticleGun->GetParticleDefinition()->GetParticleName() !=
     "opticalphoton")
  {
    G4cout << "-> warning from WLSPrimaryGeneratorAction::SetOptPhotonPolar()"
           << ":  the ParticleGun is not an opticalphoton" << G4endl;
    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);
}

SetDecayTimeConstant()

void WLSPrimaryGeneratorAction::SetDecayTimeConstant

(

G4double

time

)

Definition at line 81 of file WLSPrimaryGeneratorAction.cc.

{
  fTimeConstant = time;
}

SetUseSampledEnergy()

void WLSPrimaryGeneratorAction::SetUseSampledEnergy ( G4bool  v )

inline

Definition at line 59 of file WLSPrimaryGeneratorAction.hh.

{ fUseSampledEnergy = v; }

SetOptPhotonPolar() [2/2]

void WLSPrimaryGeneratorAction::SetOptPhotonPolar ( )

private

Definition at line 197 of file WLSPrimaryGeneratorAction.cc.

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

SetOptPhotonTime()

void WLSPrimaryGeneratorAction::SetOptPhotonTime ( )

private

Definition at line 232 of file WLSPrimaryGeneratorAction.cc.

{
  G4double time = -std::log(G4UniformRand()) * fTimeConstant;
  fParticleGun->SetParticleTime(time);
}

Member Data Documentation

fIntegralTable

G4PhysicsTable* WLSPrimaryGeneratorAction::fIntegralTable = nullptr

protected

Definition at line 62 of file WLSPrimaryGeneratorAction.hh.

fDetector

WLSDetectorConstruction* WLSPrimaryGeneratorAction::fDetector = nullptr

private

Definition at line 68 of file WLSPrimaryGeneratorAction.hh.

fParticleGun

G4GeneralParticleSource* WLSPrimaryGeneratorAction::fParticleGun = nullptr

private

Definition at line 69 of file WLSPrimaryGeneratorAction.hh.

fGunMessenger

WLSPrimaryGeneratorMessenger* WLSPrimaryGeneratorAction::fGunMessenger = nullptr

private

Definition at line 70 of file WLSPrimaryGeneratorAction.hh.

fFirst

G4bool WLSPrimaryGeneratorAction::fFirst = false

staticprivate

Definition at line 72 of file WLSPrimaryGeneratorAction.hh.

fTimeConstant

G4double WLSPrimaryGeneratorAction::fTimeConstant = 0.

private

Definition at line 73 of file WLSPrimaryGeneratorAction.hh.

fUseSampledEnergy

G4bool WLSPrimaryGeneratorAction::fUseSampledEnergy = false

private

Definition at line 74 of file WLSPrimaryGeneratorAction.hh.

---

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

• .doxygen/Doxymodules_optical.h
• extended/optical/wls/include/WLSPrimaryGeneratorAction.hh
• extended/optical/wls/src/WLSPrimaryGeneratorAction.cc