G4VAdjointReverseReaction Class Reference

#include <G4VAdjointReverseReaction.hh>

Inheritance diagram for G4VAdjointReverseReaction:

Public Member Functions
	G4VAdjointReverseReaction (G4String process_name, G4bool whichScatCase)
virtual	~G4VAdjointReverseReaction ()
void	PreparePhysicsTable (const G4ParticleDefinition &)
void	BuildPhysicsTable (const G4ParticleDefinition &)
virtual G4VParticleChange *	PostStepDoIt (const G4Track &, const G4Step &)
void	SetIntegralMode (G4bool aBool)
Protected Member Functions
virtual G4double	GetMeanFreePath (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
Protected Attributes
G4VEmAdjointModel *	theAdjointEMModel
G4ParticleChange *	fParticleChange
G4AdjointCSManager *	theAdjointCSManager
G4bool	IsScatProjToProjCase

---

Detailed Description

Definition at line 74 of file G4VAdjointReverseReaction.hh.

---

Constructor & Destructor Documentation

G4VAdjointReverseReaction::G4VAdjointReverseReaction	(	G4String	process_name,
		G4bool	whichScatCase
	)

Definition at line 45 of file G4VAdjointReverseReaction.cc.

References fParticleChange, G4AdjointCSManager::GetAdjointCSManager(), IsScatProjToProjCase, and theAdjointCSManager.

                                                                              :
                        G4VDiscreteProcess(process_name)
{theAdjointCSManager = G4AdjointCSManager::GetAdjointCSManager();
 IsScatProjToProjCase=whichScatCase;
 fParticleChange=new G4ParticleChange();
 IsFwdCSUsed=false;
 IsIntegralModeUsed=false;
 lastCS=0.;
}

G4VAdjointReverseReaction::~G4VAdjointReverseReaction

(

)

[virtual]

Definition at line 57 of file G4VAdjointReverseReaction.cc.

References fParticleChange.

{ if (fParticleChange) delete fParticleChange;
}                       

---

Member Function Documentation

void G4VAdjointReverseReaction::BuildPhysicsTable

(

const G4ParticleDefinition &

)

[virtual]

Reimplemented from G4VProcess.

Definition at line 67 of file G4VAdjointReverseReaction.cc.

References G4AdjointCSManager::BuildCrossSectionMatrices(), G4AdjointCSManager::BuildTotalSigmaTables(), and theAdjointCSManager.

{

 theAdjointCSManager->BuildCrossSectionMatrices(); //do not worry it will be done just once
 theAdjointCSManager->BuildTotalSigmaTables();

}

G4double G4VAdjointReverseReaction::GetMeanFreePath	(	const G4Track &	track,
		G4double	previousStepSize,
		G4ForceCondition *	condition
	)			[protected, virtual]

Implements G4VDiscreteProcess.

Definition at line 107 of file G4VAdjointReverseReaction.cc.

References G4VEmAdjointModel::GetAdjointCrossSection(), G4AdjointCSManager::GetCrossSectionCorrection(), G4Track::GetDefinition(), G4Track::GetKineticEnergy(), G4Track::GetMaterialCutsCouple(), IsScatProjToProjCase, NotForced, theAdjointCSManager, and theAdjointEMModel.

{ *condition = NotForced;
  G4double preStepKinEnergy = track.GetKineticEnergy();
  
  /*G4double Sigma =
                theAdjointEMModel->AdjointCrossSection(track.GetMaterialCutsCouple(),preStepKinEnergy,IsScatProjToProjCase);*/
                
  G4double Sigma =
                theAdjointEMModel->GetAdjointCrossSection(track.GetMaterialCutsCouple(),preStepKinEnergy,IsScatProjToProjCase); 
  G4double fwd_TotCS;
  Sigma *=  theAdjointCSManager->GetCrossSectionCorrection(track.GetDefinition(),preStepKinEnergy,track.GetMaterialCutsCouple(),IsFwdCSUsed, fwd_TotCS);
  //G4cout<<fwd_TotCS<<G4endl;
  /*if (IsFwdCSUsed && IsIntegralModeUsed){ //take the maximum cross section only for charged particle          
        G4double e_sigma_max, sigma_max;
        theAdjointCSManager->GetMaxFwdTotalCS(track.GetDefinition(),
                                     track.GetMaterialCutsCouple(), e_sigma_max, sigma_max);
        if (e_sigma_max > preStepKinEnergy){
                Sigma*=sigma_max/fwd_TotCS;
        }                    
  }
  */            

  G4double mean_free_path = 1.e60 *mm; 
  if (Sigma>0) mean_free_path = 1./Sigma;
  lastCS=Sigma;
  
  /*G4cout<<"Sigma  "<<Sigma<<G4endl;
  G4cout<<"mean_free_path [mm] "<<mean_free_path/mm<<G4endl;
  */
  

  return mean_free_path;
}                                        

G4VParticleChange * G4VAdjointReverseReaction::PostStepDoIt	(	const G4Track &	,
		const G4Step &
	)			[virtual]

Reimplemented from G4VDiscreteProcess.

Definition at line 76 of file G4VAdjointReverseReaction.cc.

References G4VProcess::ClearNumberOfInteractionLengthLeft(), fParticleChange, G4ParticleChange::Initialize(), IsScatProjToProjCase, G4VEmAdjointModel::SampleSecondaries(), and theAdjointEMModel.

{ 
  
  
  
  fParticleChange->Initialize(track);
 
 /* if (IsFwdCSUsed && IsIntegralModeUsed){ //INtegral mode still unstable
         G4double Tkin = step.GetPostStepPoint()->GetKineticEnergy();
         G4double fwdCS = theAdjointCSManager->GetTotalForwardCS(track.GetDefinition(), Tkin, track.GetMaterialCutsCouple());
         //G4cout<<"lastCS "<<lastCS<<G4endl;
         if (fwdCS<lastCS*G4UniformRand()) { // the reaction does not take place, same integral method as the one used for forward ionisation in  G4 
                ClearNumberOfInteractionLengthLeft();
                return fParticleChange;
         } 
         
  }
 */
 
  theAdjointEMModel->SampleSecondaries(track,
                                       IsScatProjToProjCase,
                                        fParticleChange);
  
  ClearNumberOfInteractionLengthLeft();
  return fParticleChange;
                        
   
  
}

void G4VAdjointReverseReaction::PreparePhysicsTable

(

const G4ParticleDefinition &

)

[virtual]

Reimplemented from G4VProcess.

Definition at line 62 of file G4VAdjointReverseReaction.cc.

{;
}

void G4VAdjointReverseReaction::SetIntegralMode

(

G4bool

aBool

)

[inline]

Definition at line 89 of file G4VAdjointReverseReaction.hh.

Referenced by G4eInverseBremsstrahlung::G4eInverseBremsstrahlung(), G4eInverseCompton::G4eInverseCompton(), G4eInverseIonisation::G4eInverseIonisation(), G4hInverseIonisation::G4hInverseIonisation(), G4InversePEEffect::G4InversePEEffect(), and G4IonInverseIonisation::G4IonInverseIonisation().

{IsIntegralModeUsed = aBool;}

---

Field Documentation

G4ParticleChange* G4VAdjointReverseReaction::fParticleChange [protected]

Definition at line 99 of file G4VAdjointReverseReaction.hh.

Referenced by G4VAdjointReverseReaction(), PostStepDoIt(), and ~G4VAdjointReverseReaction().

G4bool G4VAdjointReverseReaction::IsScatProjToProjCase [protected]

Definition at line 101 of file G4VAdjointReverseReaction.hh.

Referenced by G4eInverseBremsstrahlung::G4eInverseBremsstrahlung(), G4VAdjointReverseReaction(), GetMeanFreePath(), and PostStepDoIt().

G4AdjointCSManager* G4VAdjointReverseReaction::theAdjointCSManager [protected]

Definition at line 100 of file G4VAdjointReverseReaction.hh.

Referenced by BuildPhysicsTable(), G4VAdjointReverseReaction(), and GetMeanFreePath().

G4VEmAdjointModel* G4VAdjointReverseReaction::theAdjointEMModel [protected]

Definition at line 98 of file G4VAdjointReverseReaction.hh.

Referenced by G4eInverseBremsstrahlung::G4eInverseBremsstrahlung(), G4eInverseCompton::G4eInverseCompton(), G4eInverseIonisation::G4eInverseIonisation(), G4hInverseIonisation::G4hInverseIonisation(), G4InversePEEffect::G4InversePEEffect(), G4IonInverseIonisation::G4IonInverseIonisation(), GetMeanFreePath(), and PostStepDoIt().

---

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

• G4VAdjointReverseReaction.hh
• G4VAdjointReverseReaction.cc

---