G4AdjointPrimaryGenerator Class Reference

#include <G4AdjointPrimaryGenerator.hh>

Public Member Functions
	G4AdjointPrimaryGenerator ()
	~G4AdjointPrimaryGenerator ()
	G4AdjointPrimaryGenerator (const G4AdjointPrimaryGenerator &)=delete
G4AdjointPrimaryGenerator &	operator= (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)
void	SetSphericalAdjointPrimarySource (G4double radius, G4ThreeVector pos)
void	SetAdjointPrimarySourceOnAnExtSurfaceOfAVolume (const G4String &v_name)
void	ComputeAccumulatedDepthVectorAlongBackRay (G4ThreeVector glob_pos, G4ThreeVector direction, G4double ekin, G4ParticleDefinition *aPDef)
G4double	SampleDistanceAlongBackRayAndComputeWeightCorrection (G4double &weight_corr)

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();
}

~G4AdjointPrimaryGenerator()

G4AdjointPrimaryGenerator::~G4AdjointPrimaryGenerator

(

)

Definition at line 63 of file G4AdjointPrimaryGenerator.cc.

{
  delete 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 G4PhysicsOrderedFreeVector();
 
  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);
  }
}

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);
}

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);
}

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;
}

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"); 
}

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);
}

Referenced by G4AdjointPrimaryGeneratorAction::SetSphericalAdjointPrimarySource().

---

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

• G4AdjointPrimaryGenerator.hh
• G4AdjointPrimaryGenerator.cc