#include <G4AdjointPrimaryGenerator.hh>

Public Member Functions
void	ComputeAccumulatedDepthVectorAlongBackRay (G4ThreeVector glob_pos, G4ThreeVector direction, G4double ekin, G4ParticleDefinition *aPDef)

	G4AdjointPrimaryGenerator ()

	G4AdjointPrimaryGenerator (const G4AdjointPrimaryGenerator &)=delete

void	GenerateAdjointPrimaryVertex (G4Event anEvt, G4ParticleDefinition adj_part, G4double E1, G4double E2)

void	GenerateFwdPrimaryVertex (G4Event anEvt, G4ParticleDefinition adj_part, G4double E1, G4double E2)

G4AdjointPrimaryGenerator &	operator= (const G4AdjointPrimaryGenerator &)=delete

G4double	SampleDistanceAlongBackRayAndComputeWeightCorrection (G4double &weight_corr)

void	SetAdjointPrimarySourceOnAnExtSurfaceOfAVolume (const G4String &v_name)

void	SetSphericalAdjointPrimarySource (G4double radius, G4ThreeVector pos)

	~G4AdjointPrimaryGenerator ()

Private Attributes
G4ThreeVector	center_spherical_source

G4Navigator *	fLinearNavigator = nullptr

G4double	radius_spherical_source = 0.0

G4PhysicsFreeVector *	theAccumulatedDepthVector = nullptr

G4AdjointPosOnPhysVolGenerator *	theG4AdjointPosOnPhysVolGenerator = nullptr

G4SingleParticleSource *	theSingleParticleSource = nullptr

G4String	type_of_adjoint_source

Detailed Description

Definition at line 58 of file G4AdjointPrimaryGenerator.hh.

Constructor & Destructor Documentation

◆ G4AdjointPrimaryGenerator() [1/2]

G4AdjointPrimaryGenerator::G4AdjointPrimaryGenerator ( )

Definition at line 47 of file G4AdjointPrimaryGenerator.cc.

{
  center_spherical_source = G4ThreeVector(0.,0.,0.);
  type_of_adjoint_source="Spherical";
  theSingleParticleSource  = new G4SingleParticleSource();
 
  theSingleParticleSource->GetEneDist()->SetEnergyDisType("Pow");
  theSingleParticleSource->GetEneDist()->SetAlpha(-1.);
  theSingleParticleSource->GetPosDist()->SetPosDisType("Point");
  theSingleParticleSource->GetAngDist()->SetAngDistType("planar");
 
  theG4AdjointPosOnPhysVolGenerator = G4AdjointPosOnPhysVolGenerator::GetInstance();
}

References center_spherical_source, G4SingleParticleSource::GetAngDist(), G4SingleParticleSource::GetEneDist(), G4AdjointPosOnPhysVolGenerator::GetInstance(), G4SingleParticleSource::GetPosDist(), G4SPSEneDistribution::SetAlpha(), G4SPSAngDistribution::SetAngDistType(), G4SPSEneDistribution::SetEnergyDisType(), G4SPSPosDistribution::SetPosDisType(), theG4AdjointPosOnPhysVolGenerator, theSingleParticleSource, and type_of_adjoint_source.

◆ ~G4AdjointPrimaryGenerator()

G4AdjointPrimaryGenerator::~G4AdjointPrimaryGenerator ( )

Definition at line 63 of file G4AdjointPrimaryGenerator.cc.

{
  delete theSingleParticleSource;
}

References theSingleParticleSource.

◆ G4AdjointPrimaryGenerator() [2/2]

G4AdjointPrimaryGenerator::G4AdjointPrimaryGenerator ( const G4AdjointPrimaryGenerator & )

delete

Member Function Documentation

◆ ComputeAccumulatedDepthVectorAlongBackRay()

void G4AdjointPrimaryGenerator::ComputeAccumulatedDepthVectorAlongBackRay	(	G4ThreeVector	glob_pos,
		G4ThreeVector	direction,
		G4double	ekin,
		G4ParticleDefinition *	aPDef
	)

Definition at line 161 of file G4AdjointPrimaryGenerator.cc.

{
  if (fLinearNavigator == nullptr)
  {
    fLinearNavigator = G4TransportationManager::GetTransportationManager()
                     ->GetNavigatorForTracking();
  }
  G4ThreeVector position = glob_pos;
  G4double safety=1.;
  G4VPhysicalVolume* thePhysVolume =
        fLinearNavigator->LocateGlobalPointAndSetup(position);
  G4double newStep = fLinearNavigator->ComputeStep(position,direction,1.e50,
                                                   safety);
  if (theAccumulatedDepthVector != nullptr) {delete theAccumulatedDepthVector;}
  theAccumulatedDepthVector = new G4PhysicsFreeVector();
 
  G4double acc_depth=0.;
  G4double acc_length=0.;
  theAccumulatedDepthVector->InsertValues(acc_length,acc_depth);
 
  while (newStep > 0. && thePhysVolume != nullptr)
  {
    acc_length+=newStep;
    acc_depth+=newStep*thePhysVolume->GetLogicalVolume()
                                    ->GetMaterial()->GetDensity();
    theAccumulatedDepthVector->InsertValues(acc_length,acc_depth);
    position=position+newStep*direction;
    thePhysVolume = fLinearNavigator
                  ->LocateGlobalPointAndSetup(position,nullptr,false);
    newStep = fLinearNavigator->ComputeStep(position,direction,1.e50,safety);
  }
}

References G4Navigator::ComputeStep(), fLinearNavigator, G4Material::GetDensity(), G4VPhysicalVolume::GetLogicalVolume(), G4LogicalVolume::GetMaterial(), G4TransportationManager::GetNavigatorForTracking(), G4TransportationManager::GetTransportationManager(), G4PhysicsFreeVector::InsertValues(), G4Navigator::LocateGlobalPointAndSetup(), and theAccumulatedDepthVector.

◆ GenerateAdjointPrimaryVertex()

void G4AdjointPrimaryGenerator::GenerateAdjointPrimaryVertex	(	G4Event *	anEvt,
		G4ParticleDefinition *	adj_part,
		G4double	E1,
		G4double	E2
	)

Definition at line 70 of file G4AdjointPrimaryGenerator.cc.

{
  if (type_of_adjoint_source == "ExternalSurfaceOfAVolume")
  {
    // Generate position and direction relative to the external surface
    // of sensitive volume
 
    G4double costh_to_normal=1.;
    G4ThreeVector pos =G4ThreeVector(0.,0.,0.);
    G4ThreeVector direction = G4ThreeVector(0.,0.,1.);
    theG4AdjointPosOnPhysVolGenerator
      ->GenerateAPositionOnTheExtSurfaceOfThePhysicalVolume(pos, direction,
                                                            costh_to_normal);
    if (costh_to_normal <1.e-4)  { costh_to_normal = 1.e-4; }
 
    // compute now the position along the ray backward direction
    //
    theSingleParticleSource->GetAngDist()
      ->SetParticleMomentumDirection(-direction);
    theSingleParticleSource->GetPosDist()->SetCentreCoords(pos);
  }
 
  theSingleParticleSource->GetEneDist()->SetEmin(E1); 
  theSingleParticleSource->GetEneDist()->SetEmax(E2); 
 
  theSingleParticleSource->SetParticleDefinition(adj_part);
  theSingleParticleSource->GeneratePrimaryVertex(anEvent);
}

References G4AdjointPosOnPhysVolGenerator::GenerateAPositionOnTheExtSurfaceOfThePhysicalVolume(), G4SingleParticleSource::GeneratePrimaryVertex(), G4SingleParticleSource::GetAngDist(), G4SingleParticleSource::GetEneDist(), G4SingleParticleSource::GetPosDist(), pos, G4SPSPosDistribution::SetCentreCoords(), G4SPSEneDistribution::SetEmax(), G4SPSEneDistribution::SetEmin(), G4SingleParticleSource::SetParticleDefinition(), G4SPSAngDistribution::SetParticleMomentumDirection(), theG4AdjointPosOnPhysVolGenerator, theSingleParticleSource, and type_of_adjoint_source.

◆ GenerateFwdPrimaryVertex()

void G4AdjointPrimaryGenerator::GenerateFwdPrimaryVertex	(	G4Event *	anEvt,
		G4ParticleDefinition *	adj_part,
		G4double	E1,
		G4double	E2
	)

Definition at line 103 of file G4AdjointPrimaryGenerator.cc.

{
  if (type_of_adjoint_source == "ExternalSurfaceOfAVolume")
  {
    // Generate position and direction relative to the external surface
    // of sensitive volume
 
    G4double costh_to_normal=1.;
    G4ThreeVector pos =G4ThreeVector(0.,0.,0.);
    G4ThreeVector direction = G4ThreeVector(0.,0.,1.);
    theG4AdjointPosOnPhysVolGenerator
      ->GenerateAPositionOnTheExtSurfaceOfThePhysicalVolume(pos, direction,
                                                            costh_to_normal);
    if (costh_to_normal <1.e-4)  { costh_to_normal =1.e-4; }
    theSingleParticleSource->GetAngDist()
      ->SetParticleMomentumDirection(direction);
    theSingleParticleSource->GetPosDist()->SetCentreCoords(pos);
  }
 
  theSingleParticleSource->GetEneDist()->SetEmin(E1);
  theSingleParticleSource->GetEneDist()->SetEmax(E2);
 
  theSingleParticleSource->SetParticleDefinition(fwd_part);
  theSingleParticleSource->GeneratePrimaryVertex(anEvent);
}

References G4AdjointPosOnPhysVolGenerator::GenerateAPositionOnTheExtSurfaceOfThePhysicalVolume(), G4SingleParticleSource::GeneratePrimaryVertex(), G4SingleParticleSource::GetAngDist(), G4SingleParticleSource::GetEneDist(), G4SingleParticleSource::GetPosDist(), pos, G4SPSPosDistribution::SetCentreCoords(), G4SPSEneDistribution::SetEmax(), G4SPSEneDistribution::SetEmin(), G4SingleParticleSource::SetParticleDefinition(), G4SPSAngDistribution::SetParticleMomentumDirection(), theG4AdjointPosOnPhysVolGenerator, theSingleParticleSource, and type_of_adjoint_source.

Referenced by G4AdjointPrimaryGeneratorAction::GeneratePrimaries().

◆ operator=()

G4AdjointPrimaryGenerator & G4AdjointPrimaryGenerator::operator= ( const G4AdjointPrimaryGenerator & )

delete

◆ SampleDistanceAlongBackRayAndComputeWeightCorrection()

G4double G4AdjointPrimaryGenerator::SampleDistanceAlongBackRayAndComputeWeightCorrection ( G4double & weight_corr )

Definition at line 199 of file G4AdjointPrimaryGenerator.cc.

{
  G4double rand = G4UniformRand();
  G4double distance = theAccumulatedDepthVector->FindLinearEnergy(rand);
  weight_corr=1.;
  return distance;
}

References G4PhysicsVector::FindLinearEnergy(), G4UniformRand, and theAccumulatedDepthVector.

◆ SetAdjointPrimarySourceOnAnExtSurfaceOfAVolume()

void G4AdjointPrimaryGenerator::SetAdjointPrimarySourceOnAnExtSurfaceOfAVolume ( const G4String & v_name )

Definition at line 150 of file G4AdjointPrimaryGenerator.cc.

{
  theG4AdjointPosOnPhysVolGenerator->DefinePhysicalVolume1(volume_name);
  type_of_adjoint_source ="ExternalSurfaceOfAVolume";
  theSingleParticleSource->GetPosDist()->SetPosDisType("Point");
  theSingleParticleSource->GetAngDist()->SetAngDistType("planar"); 
}

References G4AdjointPosOnPhysVolGenerator::DefinePhysicalVolume1(), G4SingleParticleSource::GetAngDist(), G4SingleParticleSource::GetPosDist(), G4SPSAngDistribution::SetAngDistType(), G4SPSPosDistribution::SetPosDisType(), theG4AdjointPosOnPhysVolGenerator, theSingleParticleSource, and type_of_adjoint_source.

Referenced by G4AdjointPrimaryGeneratorAction::SetAdjointPrimarySourceOnAnExtSurfaceOfAVolume().

◆ SetSphericalAdjointPrimarySource()

void G4AdjointPrimaryGenerator::SetSphericalAdjointPrimarySource	(	G4double	radius,
		G4ThreeVector	pos
	)

Definition at line 133 of file G4AdjointPrimaryGenerator.cc.

{ 
  radius_spherical_source = radius;
  center_spherical_source = center_pos;
  type_of_adjoint_source = "Spherical"; 
  theSingleParticleSource->GetPosDist()->SetPosDisType("Surface");
  theSingleParticleSource->GetPosDist()->SetPosDisShape("Sphere");
  theSingleParticleSource->GetPosDist()->SetCentreCoords(center_pos);
  theSingleParticleSource->GetPosDist()->SetRadius(radius);
  theSingleParticleSource->GetAngDist()->SetAngDistType("cos");
  theSingleParticleSource->GetAngDist()->SetMaxTheta(pi);
  theSingleParticleSource->GetAngDist()->SetMinTheta(halfpi);
}