#include <G4AdjointPhotoElectricModel.hh>

Inheritance diagram for G4AdjointPhotoElectricModel:

Collaboration diagram for G4AdjointPhotoElectricModel:

Public Member Functions
	G4AdjointPhotoElectricModel ()

	~G4AdjointPhotoElectricModel ()

virtual void	SampleSecondaries (const G4Track &aTrack, G4bool IsScatProjToProjCase, G4ParticleChange *fParticleChange)

virtual G4double	AdjointCrossSection (const G4MaterialCutsCouple *aCouple, G4double primEnergy, G4bool IsScatProjToProjCase)

virtual G4double	GetAdjointCrossSection (const G4MaterialCutsCouple *aCouple, G4double primEnergy, G4bool IsScatProjToProjCase)

G4double	AdjointCrossSectionPerAtom (const G4Element *anElement, G4double electronEnergy)

void	SetTheDirectPEEffectModel (G4PEEffectFluoModel *aModel)

virtual void	CorrectPostStepWeight (G4ParticleChange *fParticleChange, G4double old_weight, G4double adjointPrimKinEnergy, G4double projectileKinEnergy, G4bool IsScatProjToProjCase)

Public Member Functions inherited from G4VEmAdjointModel
	G4VEmAdjointModel (const G4String &nam)

virtual	~G4VEmAdjointModel ()

virtual G4double	DiffCrossSectionPerAtomPrimToSecond (G4double kinEnergyProj, G4double kinEnergyProd, G4double Z, G4double A=0.)

virtual G4double	DiffCrossSectionPerAtomPrimToScatPrim (G4double kinEnergyProj, G4double kinEnergyScatProj, G4double Z, G4double A=0.)

virtual G4double	DiffCrossSectionPerVolumePrimToSecond (const G4Material *aMaterial, G4double kinEnergyProj, G4double kinEnergyProd)

virtual G4double	DiffCrossSectionPerVolumePrimToScatPrim (const G4Material *aMaterial, G4double kinEnergyProj, G4double kinEnergyScatProj)

virtual G4double	GetSecondAdjEnergyMaxForScatProjToProjCase (G4double PrimAdjEnergy)

virtual G4double	GetSecondAdjEnergyMinForScatProjToProjCase (G4double PrimAdjEnergy, G4double Tcut=0)

virtual G4double	GetSecondAdjEnergyMaxForProdToProjCase (G4double PrimAdjEnergy)

virtual G4double	GetSecondAdjEnergyMinForProdToProjCase (G4double PrimAdjEnergy)

void	DefineCurrentMaterial (const G4MaterialCutsCouple *couple)

std::vector< std::vector< double > *>	ComputeAdjointCrossSectionVectorPerAtomForSecond (G4double kinEnergyProd, G4double Z, G4double A=0., G4int nbin_pro_decade=10)

std::vector< std::vector< double > *>	ComputeAdjointCrossSectionVectorPerAtomForScatProj (G4double kinEnergyProd, G4double Z, G4double A=0., G4int nbin_pro_decade=10)

std::vector< std::vector< double > *>	ComputeAdjointCrossSectionVectorPerVolumeForSecond (G4Material *aMaterial, G4double kinEnergyProd, G4int nbin_pro_decade=10)

std::vector< std::vector< double > *>	ComputeAdjointCrossSectionVectorPerVolumeForScatProj (G4Material *aMaterial, G4double kinEnergyProd, G4int nbin_pro_decade=10)

void	SetCSMatrices (std::vector< G4AdjointCSMatrix * > Vec1CSMatrix, std::vector< G4AdjointCSMatrix > *Vec2CSMatrix)

G4ParticleDefinition *	GetAdjointEquivalentOfDirectPrimaryParticleDefinition ()

G4ParticleDefinition *	GetAdjointEquivalentOfDirectSecondaryParticleDefinition ()

G4double	GetHighEnergyLimit ()

G4double	GetLowEnergyLimit ()

void	SetHighEnergyLimit (G4double aVal)

void	SetLowEnergyLimit (G4double aVal)

void	DefineDirectEMModel (G4VEmModel *aModel)

void	SetAdjointEquivalentOfDirectPrimaryParticleDefinition (G4ParticleDefinition *aPart)

void	SetAdjointEquivalentOfDirectSecondaryParticleDefinition (G4ParticleDefinition *aPart)

void	SetSecondPartOfSameType (G4bool aBool)

G4bool	GetSecondPartOfSameType ()

void	SetUseMatrix (G4bool aBool)

void	SetUseMatrixPerElement (G4bool aBool)

void	SetUseOnlyOneMatrixForAllElements (G4bool aBool)

void	SetApplyCutInRange (G4bool aBool)

G4bool	GetUseMatrix ()

G4bool	GetUseMatrixPerElement ()

G4bool	GetUseOnlyOneMatrixForAllElements ()

G4bool	GetApplyCutInRange ()

G4String	GetName ()

virtual void	SetCSBiasingFactor (G4double aVal)

Private Member Functions
void	DefineCurrentMaterialAndElectronEnergy (const G4MaterialCutsCouple *aCouple, G4double eEnergy)

Private Attributes
G4double	xsec [40]

G4double	totAdjointCS

G4double	totBiasedAdjointCS

G4double	factorCSBiasing

G4double	pre_step_AdjointCS

G4double	post_step_AdjointCS

G4double	shell_prob [40][40]

G4PEEffectFluoModel *	theDirectPEEffectModel

size_t	index_element

G4double	current_eEnergy

Additional Inherited Members
Protected Member Functions inherited from G4VEmAdjointModel
G4double	DiffCrossSectionFunction1 (G4double kinEnergyProj)

G4double	DiffCrossSectionFunction2 (G4double kinEnergyProj)

G4double	DiffCrossSectionPerVolumeFunctionForIntegrationOverEkinProj (G4double EkinProd)

G4double	SampleAdjSecEnergyFromCSMatrix (size_t MatrixIndex, G4double prim_energy, G4bool IsScatProjToProjCase)

G4double	SampleAdjSecEnergyFromCSMatrix (G4double prim_energy, G4bool IsScatProjToProjCase)

void	SelectCSMatrix (G4bool IsScatProjToProjCase)

virtual G4double	SampleAdjSecEnergyFromDiffCrossSectionPerAtom (G4double prim_energy, G4bool IsScatProjToProjCase)

Protected Attributes inherited from G4VEmAdjointModel
G4VEmModel *	theDirectEMModel

G4VParticleChange *	pParticleChange

const G4String	name

G4int	ASelectedNucleus

G4int	ZSelectedNucleus

G4Material *	SelectedMaterial

G4double	kinEnergyProdForIntegration

G4double	kinEnergyScatProjForIntegration

G4double	kinEnergyProjForIntegration

std::vector< G4AdjointCSMatrix >	pOnCSMatrixForProdToProjBackwardScattering

std::vector< G4AdjointCSMatrix >	pOnCSMatrixForScatProjToProjBackwardScattering

std::vector< G4double >	CS_Vs_ElementForScatProjToProjCase

std::vector< G4double >	CS_Vs_ElementForProdToProjCase

G4double	lastCS

G4double	lastAdjointCSForScatProjToProjCase

G4double	lastAdjointCSForProdToProjCase

G4ParticleDefinition *	theAdjEquivOfDirectPrimPartDef

G4ParticleDefinition *	theAdjEquivOfDirectSecondPartDef

G4ParticleDefinition *	theDirectPrimaryPartDef

G4bool	second_part_of_same_type

G4double	preStepEnergy

G4Material *	currentMaterial

G4MaterialCutsCouple *	currentCouple

size_t	currentMaterialIndex

size_t	currentCoupleIndex

G4double	currentTcutForDirectPrim

G4double	currentTcutForDirectSecond

G4bool	ApplyCutInRange

G4double	mass_ratio_product

G4double	mass_ratio_projectile

G4double	HighEnergyLimit

G4double	LowEnergyLimit

G4double	CS_biasing_factor

G4bool	UseMatrix

G4bool	UseMatrixPerElement

G4bool	UseOnlyOneMatrixForAllElements

size_t	indexOfUsedCrossSectionMatrix

size_t	model_index

Detailed Description

Definition at line 62 of file G4AdjointPhotoElectricModel.hh.

Constructor & Destructor Documentation

◆ G4AdjointPhotoElectricModel()

G4AdjointPhotoElectricModel::G4AdjointPhotoElectricModel ( )

Definition at line 42 of file G4AdjointPhotoElectricModel.cc.

                                                         :
  G4VEmAdjointModel("AdjointPEEffect")
 
 { SetUseMatrix(false);
   SetApplyCutInRange(false);
 
   //Initialization
   current_eEnergy =0.;
   totAdjointCS=0.;
   factorCSBiasing =1.;
   post_step_AdjointCS =0.;
   pre_step_AdjointCS =0.;
   totBiasedAdjointCS =0.;
 
   index_element=0;
 
   theAdjEquivOfDirectPrimPartDef =G4AdjointGamma::AdjointGamma();
   theAdjEquivOfDirectSecondPartDef=G4AdjointElectron::AdjointElectron();
   theDirectPrimaryPartDef=G4Gamma::Gamma();
   second_part_of_same_type=false;
   theDirectPEEffectModel = new G4PEEffectFluoModel();
 }

Here is the call graph for this function:

◆ ~G4AdjointPhotoElectricModel()

G4AdjointPhotoElectricModel::~G4AdjointPhotoElectricModel ( )

Definition at line 66 of file G4AdjointPhotoElectricModel.cc.

67 {;}

Member Function Documentation

◆ AdjointCrossSection()

G4double G4AdjointPhotoElectricModel::AdjointCrossSection	(	const G4MaterialCutsCouple *	aCouple,
		G4double	primEnergy,
		G4bool	IsScatProjToProjCase
	)

virtual

Reimplemented from G4VEmAdjointModel.

Definition at line 191 of file G4AdjointPhotoElectricModel.cc.

 { 
   
 
   if (IsScatProjToProjCase)  return 0.;
 
     
   if (aCouple !=currentCouple || current_eEnergy !=electronEnergy) {
     totAdjointCS = 0.;
     DefineCurrentMaterialAndElectronEnergy(aCouple, electronEnergy);
     const G4ElementVector* theElementVector = currentMaterial->GetElementVector();
     const double* theAtomNumDensityVector = currentMaterial->GetVecNbOfAtomsPerVolume();
     size_t nelm =  currentMaterial->GetNumberOfElements();
     for (index_element=0;index_element<nelm;index_element++){
         
         totAdjointCS +=AdjointCrossSectionPerAtom((*theElementVector)[index_element],electronEnergy)*theAtomNumDensityVector[index_element];
         xsec[index_element] = totAdjointCS;
     } 
 
     totBiasedAdjointCS=std::min(totAdjointCS,0.01);
 //  totBiasedAdjointCS=totAdjointCS;
     factorCSBiasing = totBiasedAdjointCS/totAdjointCS;
     lastCS=totBiasedAdjointCS;
     
     
   }
   return totBiasedAdjointCS;
 
   
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ AdjointCrossSectionPerAtom()

G4double G4AdjointPhotoElectricModel::AdjointCrossSectionPerAtom	(	const G4Element *	anElement,
		G4double	electronEnergy
	)

Definition at line 234 of file G4AdjointPhotoElectricModel.cc.

 { 
   G4int nShells = anElement->GetNbOfAtomicShells();
   G4double Z= anElement->GetZ();
   G4int i  = 0;  
   G4double B0=anElement->GetAtomicShell(0);
   G4double gammaEnergy = electronEnergy+B0;
   G4double CS= theDirectPEEffectModel->ComputeCrossSectionPerAtom(G4Gamma::Gamma(),gammaEnergy,Z,0.,0.,0.);
   G4double adjointCS =0.;
   if (CS >0) adjointCS += CS/gammaEnergy; 
   shell_prob[index_element][0] = adjointCS;                                          
   for (i=1;i<nShells;i++){
     //G4cout<<i<<G4endl;
     G4double Bi_= anElement->GetAtomicShell(i-1);
     G4double Bi = anElement->GetAtomicShell(i);
     //G4cout<<Bi_<<'\t'<<Bi<<G4endl;
     if (electronEnergy <Bi_-Bi) {
         gammaEnergy = electronEnergy+Bi;
         
         CS=theDirectPEEffectModel->ComputeCrossSectionPerAtom(G4Gamma::Gamma(),gammaEnergy,Z,0.,0.,0.);
         if (CS>0) adjointCS +=CS/gammaEnergy;
     }
     shell_prob[index_element][i] = adjointCS;   
   
   }
   adjointCS*=electronEnergy;
   return adjointCS;
   
 }               

Here is the call graph for this function:

Here is the caller graph for this function:

◆ CorrectPostStepWeight()

void G4AdjointPhotoElectricModel::CorrectPostStepWeight	(	G4ParticleChange *	fParticleChange,
		G4double	old_weight,
		G4double	adjointPrimKinEnergy,
		G4double	projectileKinEnergy,
		G4bool	IsScatProjToProjCase
	)

virtual

Reimplemented from G4VEmAdjointModel.

Definition at line 169 of file G4AdjointPhotoElectricModel.cc.

 {
  G4double new_weight=old_weight;
 
  G4double w_corr =G4AdjointCSManager::GetAdjointCSManager()->GetPostStepWeightCorrection()/factorCSBiasing;
  w_corr*=post_step_AdjointCS/pre_step_AdjointCS; 
 
 
  new_weight*=w_corr; 
  new_weight*=projectileKinEnergy/adjointPrimKinEnergy;
  fParticleChange->SetParentWeightByProcess(false);
  fParticleChange->SetSecondaryWeightByProcess(false);
  fParticleChange->ProposeParentWeight(new_weight);  
 }   

Here is the call graph for this function:

Here is the caller graph for this function:

◆ DefineCurrentMaterialAndElectronEnergy()

void G4AdjointPhotoElectricModel::DefineCurrentMaterialAndElectronEnergy	(	const G4MaterialCutsCouple *	aCouple,
		G4double	eEnergy
	)

private

Definition at line 266 of file G4AdjointPhotoElectricModel.cc.

 { currentCouple   = const_cast<G4MaterialCutsCouple*> (couple);
   currentMaterial = const_cast<G4Material*> (couple->GetMaterial());
   currentCoupleIndex = couple->GetIndex();
   currentMaterialIndex = currentMaterial->GetIndex();
   current_eEnergy = anEnergy;   
   theDirectPEEffectModel->SetCurrentCouple(couple);
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ GetAdjointCrossSection()

G4double G4AdjointPhotoElectricModel::GetAdjointCrossSection	(	const G4MaterialCutsCouple *	aCouple,
		G4double	primEnergy,
		G4bool	IsScatProjToProjCase
	)

virtual

Reimplemented from G4VEmAdjointModel.

Definition at line 226 of file G4AdjointPhotoElectricModel.cc.

229 { return AdjointCrossSection(aCouple,electronEnergy,IsScatProjToProjCase);

230 }

G4AdjointPhotoElectricModel::AdjointCrossSection

virtual G4double AdjointCrossSection(const G4MaterialCutsCouple *aCouple, G4double primEnergy, G4bool IsScatProjToProjCase)

Definition: G4AdjointPhotoElectricModel.cc:191

Here is the call graph for this function:

◆ SampleSecondaries()

void G4AdjointPhotoElectricModel::SampleSecondaries	(	const G4Track &	aTrack,
		G4bool	IsScatProjToProjCase,
		G4ParticleChange *	fParticleChange
	)

virtual

Implements G4VEmAdjointModel.

Definition at line 71 of file G4AdjointPhotoElectricModel.cc.

 { if (IsScatProjToProjCase) return ;
 
   //Compute the totAdjointCS vectors if not already done for the current couple and electron energy
   //-----------------------------------------------------------------------------------------------
   const G4MaterialCutsCouple* aCouple = aTrack.GetMaterialCutsCouple();
   const G4DynamicParticle* aDynPart =  aTrack.GetDynamicParticle() ;
   G4double electronEnergy = aDynPart->GetKineticEnergy();
   G4ThreeVector electronDirection= aDynPart->GetMomentumDirection() ;
   pre_step_AdjointCS = totAdjointCS; //The last computed CS was  at pre step point
   post_step_AdjointCS = AdjointCrossSection(aCouple, electronEnergy,IsScatProjToProjCase);
   post_step_AdjointCS = totAdjointCS; 
                 
 
 
 
   //Sample element
   //-------------
    const G4ElementVector* theElementVector = currentMaterial->GetElementVector();
    size_t nelm =  currentMaterial->GetNumberOfElements();
    G4double rand_CS= G4UniformRand()*xsec[nelm-1];
    for (index_element=0; index_element<nelm-1; index_element++){
     if (rand_CS<xsec[index_element]) break;
    }
     
    //Sample shell and binding energy
    //-------------
    G4int nShells = (*theElementVector)[index_element]->GetNbOfAtomicShells();
    rand_CS= shell_prob[index_element][nShells-1]*G4UniformRand();
    G4int i  = 0;  
    for (i=0; i<nShells-1; i++){
     if (rand_CS<shell_prob[index_element][i]) break;
    }
    G4double gammaEnergy= electronEnergy+(*theElementVector)[index_element]->GetAtomicShell(i);
     
   //Sample cos theta
   //Copy of the G4PEEfectFluoModel cos theta sampling method ElecCosThetaDistribution.   
   //This method cannot be used directly from G4PEEfectFluoModel because it is a friend method. I should ask Vladimir to change that  
   //------------------------------------------------------------------------------------------------    
   //G4double cos_theta = theDirectPEEffectModel->ElecCosThetaDistribution(electronEnergy);
     
    G4double  cos_theta = 1.;
    G4double gamma   = 1. + electronEnergy/electron_mass_c2;
    if (gamma <= 5.) {
     G4double beta  = std::sqrt(gamma*gamma-1.)/gamma;
     G4double b     = 0.5*gamma*(gamma-1.)*(gamma-2);
 
     G4double rndm,term,greject,grejsup;
     if (gamma < 2.) grejsup = gamma*gamma*(1.+b-beta*b);
     else            grejsup = gamma*gamma*(1.+b+beta*b);
 
     do {    rndm = 1.-2*G4UniformRand();
             cos_theta = (rndm+beta)/(rndm*beta+1.);
             term = 1.-beta*cos_theta;
             greject = (1.-cos_theta*cos_theta)*(1.+b*term)/(term*term);
           // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
     } while(greject < G4UniformRand()*grejsup);
   }
     
   // direction of the adjoint gamma electron
   //---------------------------------------
   
      
   G4double sin_theta = std::sqrt(1.-cos_theta*cos_theta);
   G4double Phi     = twopi * G4UniformRand();
   G4double dirx = sin_theta*std::cos(Phi),diry = sin_theta*std::sin(Phi),dirz = cos_theta;
   G4ThreeVector adjoint_gammaDirection(dirx,diry,dirz);
   adjoint_gammaDirection.rotateUz(electronDirection);
   
   
   
   //Weight correction
  //-----------------------                     
   CorrectPostStepWeight(fParticleChange, aTrack.GetWeight(), electronEnergy,gammaEnergy,IsScatProjToProjCase);  
  
   
   
   //Create secondary and modify fParticleChange 
   //--------------------------------------------
   G4DynamicParticle* anAdjointGamma = new G4DynamicParticle (
                        G4AdjointGamma::AdjointGamma(),adjoint_gammaDirection, gammaEnergy);
   
   
   
  
   
   fParticleChange->ProposeTrackStatus(fStopAndKill);
   fParticleChange->AddSecondary(anAdjointGamma);    
         
 
   
 
 }