G4AdjointPrimaryGeneratorAction Class Reference

#include <G4AdjointPrimaryGeneratorAction.hh>

Inheritance diagram for G4AdjointPrimaryGeneratorAction:

Public Member Functions
	G4AdjointPrimaryGeneratorAction ()
	~G4AdjointPrimaryGeneratorAction ()
void	GeneratePrimaries (G4Event *)
void	SetRndmFlag (const G4String &val)
void	SetEmin (G4double val)
void	SetEmax (G4double val)
void	SetEminIon (G4double val)
void	SetEmaxIon (G4double val)
void	SetSphericalAdjointPrimarySource (G4double radius, G4ThreeVector pos)
void	SetAdjointPrimarySourceOnAnExtSurfaceOfAVolume (const G4String &volume_name)
void	ConsiderParticleAsPrimary (const G4String &particle_name)
void	NeglectParticleAsPrimary (const G4String &particle_name)
void	SetPrimaryIon (G4ParticleDefinition *adjointIon, G4ParticleDefinition *fwdIon)
void	UpdateListOfPrimaryParticles ()
size_t	GetNbOfAdjointPrimaryTypes ()
std::vector< G4ParticleDefinition * >	GetListOfPrimaryFwdParticles ()
const G4String &	GetPrimaryIonName ()
Public Member Functions inherited from G4VUserPrimaryGeneratorAction
	G4VUserPrimaryGeneratorAction ()
virtual	~G4VUserPrimaryGeneratorAction ()
virtual void	GeneratePrimaries (G4Event *anEvent)=0

Detailed Description

Definition at line 78 of file G4AdjointPrimaryGeneratorAction.hh.

Constructor & Destructor Documentation

G4AdjointPrimaryGeneratorAction()

G4AdjointPrimaryGeneratorAction::G4AdjointPrimaryGeneratorAction

(

)

Definition at line 47 of file G4AdjointPrimaryGeneratorAction.cc.

  : Emin(0.), Emax(0.), EminIon(0.), EmaxIon(0.), NbOfAdjointPrimaryTypes(0),
    index_particle(100000), last_generated_part_was_adjoint(false),
    radius_spherical_source(0.), fwd_ion(0), adj_ion(0), 
    ion_name("not_defined")
{
  theAdjointPrimaryGenerator= new G4AdjointPrimaryGenerator();
 
  PrimariesConsideredInAdjointSim[G4String("e-")]=false;
  PrimariesConsideredInAdjointSim[G4String("gamma")]=false;
  PrimariesConsideredInAdjointSim[G4String("proton")]=false;
  PrimariesConsideredInAdjointSim[G4String("ion")]=false;
 
  ListOfPrimaryFwdParticles.clear();
  ListOfPrimaryAdjParticles.clear();
}

~G4AdjointPrimaryGeneratorAction()

G4AdjointPrimaryGeneratorAction::~G4AdjointPrimaryGeneratorAction

(

)

Definition at line 65 of file G4AdjointPrimaryGeneratorAction.cc.

{
  delete theAdjointPrimaryGenerator;
}

Member Function Documentation

ConsiderParticleAsPrimary()

void G4AdjointPrimaryGeneratorAction::ConsiderParticleAsPrimary

(

const G4String &

particle_name

)

Definition at line 194 of file G4AdjointPrimaryGeneratorAction.cc.

{
  if (PrimariesConsideredInAdjointSim.find(particle_name) != PrimariesConsideredInAdjointSim.end()){
    PrimariesConsideredInAdjointSim[particle_name]=true;
  }
  UpdateListOfPrimaryParticles();
}

Referenced by G4AdjointSimManager::ConsiderParticleAsPrimary().

GeneratePrimaries()

void G4AdjointPrimaryGeneratorAction::GeneratePrimaries ( G4Event *  anEvent )

virtual

Implements G4VUserPrimaryGeneratorAction.

Definition at line 71 of file G4AdjointPrimaryGeneratorAction.cc.

{
   if ( !last_generated_part_was_adjoint ) {
     
     index_particle++;
     if (index_particle >= ListOfPrimaryAdjParticles.size()) index_particle =0;
     
    
     G4double E1=Emin;
     G4double E2=Emax;
     if (!ListOfPrimaryAdjParticles[index_particle]) UpdateListOfPrimaryParticles();//ion has not been created yet
    
     if (ListOfPrimaryAdjParticles[index_particle]->GetParticleName() == "adj_proton") {
        E1=EminIon;
        E2=EmaxIon;
     }
     if (ListOfPrimaryAdjParticles[index_particle]->GetParticleType() == "adjoint_nucleus") {
        G4int A= ListOfPrimaryAdjParticles[index_particle]->GetAtomicMass();
        E1=EminIon*A;
        E2=EmaxIon*A;
     }
     theAdjointPrimaryGenerator->GenerateAdjointPrimaryVertex(anEvent,
                                  ListOfPrimaryAdjParticles[index_particle],
                                  E1,E2);
     G4PrimaryVertex* aPrimVertex = anEvent->GetPrimaryVertex();
  
  
      p=aPrimVertex->GetPrimary()->GetMomentum();
      pos=aPrimVertex->GetPosition();
      G4double pmag=p.mag();
      
      G4double m0=ListOfPrimaryAdjParticles[index_particle]->GetPDGMass();
      G4double ekin=std::sqrt( m0*m0 + pmag*pmag) -m0;
    
  
      //The factor pi is to normalise the weight to the directional flux
      G4double adjoint_source_area = G4AdjointSimManager::GetInstance()->GetAdjointSourceArea();
      G4double adjoint_weight = ComputeEnergyDistWeight(ekin,E1,E2)*adjoint_source_area*pi;
 
      aPrimVertex->SetWeight(adjoint_weight);
 
      last_generated_part_was_adjoint =true;
      G4AdjointSimManager::GetInstance()->SetAdjointTrackingMode(true);
      G4AdjointSimManager::GetInstance()->RegisterAdjointPrimaryWeight(adjoint_weight);
   }
   else {
          //fwd particle equivalent to the last generated adjoint particle ios generated    
      G4PrimaryVertex* aPrimVertex = new G4PrimaryVertex();
      aPrimVertex->SetPosition(pos.x(),pos.y(),pos.z());
      aPrimVertex->SetT0(0.);
      G4PrimaryParticle* aPrimParticle = new G4PrimaryParticle(ListOfPrimaryFwdParticles[index_particle],
                               -p.x(),-p.y(),-p.z()); 
    
      aPrimVertex->SetPrimary(aPrimParticle);
      anEvent->AddPrimaryVertex(aPrimVertex);                          
      last_generated_part_was_adjoint =false;
      G4AdjointSimManager::GetInstance()->SetAdjointTrackingMode(false);
   }           
}

GetListOfPrimaryFwdParticles()

std::vector< G4ParticleDefinition * > G4AdjointPrimaryGeneratorAction::GetListOfPrimaryFwdParticles ( )

inline

Definition at line 100 of file G4AdjointPrimaryGeneratorAction.hh.

{return ListOfPrimaryFwdParticles;}

Referenced by G4AdjointSimManager::GetListOfPrimaryFwdParticles(), and G4AdjointSimManager::RegisterAtEndOfAdjointTrack().

GetNbOfAdjointPrimaryTypes()

size_t G4AdjointPrimaryGeneratorAction::GetNbOfAdjointPrimaryTypes ( )

inline

Definition at line 99 of file G4AdjointPrimaryGeneratorAction.hh.

{return ListOfPrimaryAdjParticles.size();}

Referenced by G4AdjointSimManager::RunAdjointSimulation().

GetPrimaryIonName()

const G4String & G4AdjointPrimaryGeneratorAction::GetPrimaryIonName ( )

inline

Definition at line 101 of file G4AdjointPrimaryGeneratorAction.hh.

{return ion_name;}

Referenced by G4AdjointSimManager::GetPrimaryIonName().

NeglectParticleAsPrimary()

void G4AdjointPrimaryGeneratorAction::NeglectParticleAsPrimary

(

const G4String &

particle_name

)

Definition at line 203 of file G4AdjointPrimaryGeneratorAction.cc.

{
  if (PrimariesConsideredInAdjointSim.find(particle_name) != PrimariesConsideredInAdjointSim.end()){
    PrimariesConsideredInAdjointSim[particle_name]= false;
  }
  UpdateListOfPrimaryParticles();
}

Referenced by G4AdjointSimManager::NeglectParticleAsPrimary().

SetAdjointPrimarySourceOnAnExtSurfaceOfAVolume()

void G4AdjointPrimaryGeneratorAction::SetAdjointPrimarySourceOnAnExtSurfaceOfAVolume

(

const G4String &

volume_name

)

Definition at line 187 of file G4AdjointPrimaryGeneratorAction.cc.

{
  type_of_adjoint_source ="ExternalSurfaceOfAVolume";
  theAdjointPrimaryGenerator->SetAdjointPrimarySourceOnAnExtSurfaceOfAVolume(volume_name);
}

Referenced by G4AdjointSimManager::DefineAdjointSourceOnTheExtSurfaceOfAVolume().

SetEmax()

void G4AdjointPrimaryGeneratorAction::SetEmax

(

G4double

val

)

Definition at line 139 of file G4AdjointPrimaryGeneratorAction.cc.

{
  Emax=val;
  EmaxIon=val;
}

Referenced by G4AdjointSimManager::SetAdjointSourceEmax().

SetEmaxIon()

void G4AdjointPrimaryGeneratorAction::SetEmaxIon

(

G4double

val

)

Definition at line 152 of file G4AdjointPrimaryGeneratorAction.cc.

{
  EmaxIon=val;
}

SetEmin()

void G4AdjointPrimaryGeneratorAction::SetEmin

(

G4double

val

)

Definition at line 132 of file G4AdjointPrimaryGeneratorAction.cc.

{
  Emin=val;
  EminIon=val;
}

Referenced by G4AdjointSimManager::SetAdjointSourceEmin().

SetEminIon()

void G4AdjointPrimaryGeneratorAction::SetEminIon

(

G4double

val

)

Definition at line 146 of file G4AdjointPrimaryGeneratorAction.cc.

{
  EminIon=val;
}

SetPrimaryIon()

void G4AdjointPrimaryGeneratorAction::SetPrimaryIon	(	G4ParticleDefinition *	adjointIon,
		G4ParticleDefinition *	fwdIon
	)

Definition at line 244 of file G4AdjointPrimaryGeneratorAction.cc.

{
  fwd_ion = fwdIon;
  adj_ion = adjointIon;
  UpdateListOfPrimaryParticles();
}

Referenced by G4AdjointSimManager::SetPrimaryIon().

SetRndmFlag()

void G4AdjointPrimaryGeneratorAction::SetRndmFlag ( const G4String &  val )

inline

Definition at line 88 of file G4AdjointPrimaryGeneratorAction.hh.

{ rndmFlag = val;}

SetSphericalAdjointPrimarySource()

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

Definition at line 178 of file G4AdjointPrimaryGeneratorAction.cc.

{ 
  radius_spherical_source = radius;
  center_spherical_source = center_pos;
  type_of_adjoint_source ="Spherical";
  theAdjointPrimaryGenerator->SetSphericalAdjointPrimarySource(radius,center_pos);
}

Referenced by G4AdjointSimManager::DefineSphericalAdjointSource(), and G4AdjointSimManager::DefineSphericalAdjointSourceWithCentreAtTheCentreOfAVolume().

UpdateListOfPrimaryParticles()

void G4AdjointPrimaryGeneratorAction::UpdateListOfPrimaryParticles

(

)

Definition at line 212 of file G4AdjointPrimaryGeneratorAction.cc.

{
    G4ParticleTable* theParticleTable = G4ParticleTable::GetParticleTable();
    ListOfPrimaryFwdParticles.clear();
    ListOfPrimaryAdjParticles.clear();
    std::map<G4String, G4bool>::iterator iter;
    for( iter = PrimariesConsideredInAdjointSim.begin(); iter != PrimariesConsideredInAdjointSim.end(); ++iter ) {
    if(iter->second) {
        G4String fwd_particle_name = iter->first;
        if ( fwd_particle_name != "ion") {
            G4String adj_particle_name = G4String("adj_") + fwd_particle_name;
            ListOfPrimaryFwdParticles.push_back(theParticleTable->FindParticle(fwd_particle_name));
            ListOfPrimaryAdjParticles.push_back(theParticleTable->FindParticle(adj_particle_name));
        }
        else {
            if (fwd_ion ){
                ion_name=fwd_ion->GetParticleName();
                G4String adj_ion_name=G4String("adj_") +ion_name;
                ListOfPrimaryFwdParticles.push_back(fwd_ion);
                ListOfPrimaryAdjParticles.push_back(adj_ion);
            }
            else {
                ListOfPrimaryFwdParticles.push_back(0);
                ListOfPrimaryAdjParticles.push_back(0);
                
            }   
        }
    }   
   }
}

Referenced by ConsiderParticleAsPrimary(), GeneratePrimaries(), NeglectParticleAsPrimary(), G4AdjointSimManager::RunAdjointSimulation(), and SetPrimaryIon().

---

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

• G4AdjointPrimaryGeneratorAction.hh
• G4AdjointPrimaryGeneratorAction.cc