G4AdjointPrimaryGenerator Class Reference

#include <G4AdjointPrimaryGenerator.hh>

Public Member Functions
	G4AdjointPrimaryGenerator ()
	~G4AdjointPrimaryGenerator ()
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 &volume_name)
void	ComputeAccumulatedDepthVectorAlongBackRay (G4ThreeVector glob_pos, G4ThreeVector direction, G4double ekin, G4ParticleDefinition *aPartDef)
G4double	SampleDistanceAlongBackRayAndComputeWeightCorrection (G4double &weight_corr)

Detailed Description

Definition at line 67 of file G4AdjointPrimaryGenerator.hh.

Constructor & Destructor Documentation

G4AdjointPrimaryGenerator::G4AdjointPrimaryGenerator

(

)

Definition at line 53 of file G4AdjointPrimaryGenerator.cc.

  : radius_spherical_source(0.),fLinearNavigator(0),theAccumulatedDepthVector(0)
{
  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();

}

Here is the call graph for this function:

G4AdjointPrimaryGenerator::~G4AdjointPrimaryGenerator

(

)

Definition at line 70 of file G4AdjointPrimaryGenerator.cc.

{
  delete theSingleParticleSource;
}

Member Function Documentation

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

Definition at line 154 of file G4AdjointPrimaryGenerator.cc.

{ if (!fLinearNavigator) 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) delete theAccumulatedDepthVector;
 theAccumulatedDepthVector = new G4PhysicsOrderedFreeVector();
 //if (theAccumulatedCSDepthVector) delete theAccumulatedCSDepthVector;
 //theAccumulatedCSDepthVector = new G4PhysicsOrderedFreeVector();

 G4double acc_depth=0.;
 G4double acc_length=0.;
 //G4double acc_cs_depth=0.;
 //theAccumulatedCSDepthVector->InsertValues(acc_cs_depth, acc_length);
 theAccumulatedDepthVector->InsertValues(acc_length,acc_depth);

 while (newStep > 0. && thePhysVolume) {
   acc_length+=newStep;
   /*
   const G4MaterialCutsCouple* theMatCutsCouple=
            thePhysVolume->GetLogicalVolume()->GetMaterialCutsCouple();


   acc_cs_depth+=newStep*G4AdjointCSManager::GetAdjointCSManager()->GetTotalAdjointCS(aPartDef,
                                                             ekin,
                                                        theMatCutsCouple);
   theAccumulatedCSDepthVector->InsertValues(acc_cs_depth, acc_length);*/

   acc_depth+=newStep*thePhysVolume->GetLogicalVolume()->GetMaterial()->GetDensity();
   theAccumulatedDepthVector->InsertValues(acc_length,acc_depth);
   position=position+newStep*direction;
   thePhysVolume =
          fLinearNavigator->LocateGlobalPointAndSetup(position,0,false);
   newStep =fLinearNavigator->ComputeStep(position,direction,1.e50,
                     safety);
 }


}

Here is the call graph for this function:

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

Definition at line 76 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);



}

Here is the call graph for this function:

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

Definition at line 105 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);
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4AdjointPrimaryGenerator::SampleDistanceAlongBackRayAndComputeWeightCorrection

(

G4double &

weight_corr

)

Definition at line 204 of file G4AdjointPrimaryGenerator.cc.

{G4double rand = G4UniformRand();
 G4double distance = theAccumulatedDepthVector->FindLinearEnergy(rand);
 /*
 G4double acc_cs_depth=theAccumulatedCSDepthVector->GetEnergy(distance);
 weight_corr=std::exp(-acc_cs_depth);*/
 weight_corr=1.;
 return distance;
}

Here is the call graph for this function:

void G4AdjointPrimaryGenerator::SetAdjointPrimarySourceOnAnExtSurfaceOfAVolume

(

const G4String &

volume_name

)

Definition at line 144 of file G4AdjointPrimaryGenerator.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 129 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);
}

Here is the call graph for this function:

Here is the caller graph for this function:

---

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

• geant4.10.03.p01/source/event/include/G4AdjointPrimaryGenerator.hh
• geant4.10.03.p01/source/event/src/G4AdjointPrimaryGenerator.cc