G4AdjointAlongStepWeightCorrection Class Reference

#include <G4AdjointAlongStepWeightCorrection.hh>

Inheritance diagram for G4AdjointAlongStepWeightCorrection:

Public Member Functions
	G4AdjointAlongStepWeightCorrection (const G4String &name="ContinuousWeightCorrection", G4ProcessType type=fElectromagnetic)
virtual	~G4AdjointAlongStepWeightCorrection ()
void	PreparePhysicsTable (const G4ParticleDefinition &)
void	BuildPhysicsTable (const G4ParticleDefinition &)
G4VParticleChange *	AlongStepDoIt (const G4Track &, const G4Step &)
Protected Member Functions
virtual G4double	GetContinuousStepLimit (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &currentSafety)
Protected Attributes
G4ParticleChange *	fParticleChange

---

Detailed Description

Definition at line 70 of file G4AdjointAlongStepWeightCorrection.hh.

---

Constructor & Destructor Documentation

G4AdjointAlongStepWeightCorrection::G4AdjointAlongStepWeightCorrection	(	const G4String &	name = "ContinuousWeightCorrection",
		G4ProcessType	type = fElectromagnetic
	)

Definition at line 37 of file G4AdjointAlongStepWeightCorrection.cc.

References fParticleChange.

                     : G4VContinuousProcess(name, type)
{fParticleChange = new G4ParticleChange();
 currentMaterialIndex=0;
 preStepKinEnergy=1.;
 currentCouple=0;
}

G4AdjointAlongStepWeightCorrection::~G4AdjointAlongStepWeightCorrection

(

)

[virtual]

Definition at line 48 of file G4AdjointAlongStepWeightCorrection.cc.

References fParticleChange.

{delete fParticleChange;
}

---

Member Function Documentation

G4VParticleChange * G4AdjointAlongStepWeightCorrection::AlongStepDoIt	(	const G4Track &	,
		const G4Step &
	)			[virtual]

Reimplemented from G4VContinuousProcess.

Definition at line 66 of file G4AdjointAlongStepWeightCorrection.cc.

References fParticleChange, G4AdjointCSManager::GetAdjointCSManager(), G4AdjointCSManager::GetContinuousWeightCorrection(), G4DynamicParticle::GetDefinition(), G4Track::GetDynamicParticle(), G4StepPoint::GetKineticEnergy(), G4Step::GetPostStepPoint(), G4Step::GetStepLength(), G4StepPoint::GetWeight(), G4ParticleChange::Initialize(), G4VParticleChange::ProposeParentWeight(), G4VParticleChange::SetParentWeightByProcess(), and G4VParticleChange::SetSecondaryWeightByProcess().

{
   
  fParticleChange->Initialize(track);
  
  // Get the actual (true) Step length
  //----------------------------------
  G4double length = step.GetStepLength();
 

  G4double Tkin = step.GetPostStepPoint()->GetKineticEnergy();
  G4ParticleDefinition* thePartDef= const_cast<G4ParticleDefinition*>  (track.GetDynamicParticle()->GetDefinition());
  G4double weight_correction=G4AdjointCSManager::GetAdjointCSManager()->GetContinuousWeightCorrection(thePartDef,
                                                                        preStepKinEnergy,Tkin, currentCouple,length);
        
  
  

  //Caution!!!
  // It is important  to select the weight of the post_step_point
  // as the current weight and not the weight of the track, as t
  // the  weight of the track is changed after having applied all
  // the along_step_do_it.

  // G4double new_weight=weight_correction*track.GetWeight(); //old
  G4double new_weight=weight_correction*step.GetPostStepPoint()->GetWeight();



  //if (weight_correction >2.) new_weight=1.e-300;
  
  
  //The following test check for zero weight.
  //This happens after weight correction of gamma for photo electric effect.
  //When the new weight is 0 it will be later on consider as nan by G4.
  //Therefore we do put a lower limit of 1.e-300. for new_weight 
  //Correction by L.Desorgher on 15 July 2009 
  if (new_weight==0 || (new_weight<=0 && new_weight>0)){
                //G4cout<<new_weight<<'\t'<<weight_correction<<'\t'<<track.GetWeight()<<G4endl;
                new_weight=1.e-300;
  }
  
  //G4cout<<new_weight<<'\t'<<weight_correction<<'\t'<<track.GetWeight()<<G4endl;
  fParticleChange->SetParentWeightByProcess(false);
  fParticleChange->SetSecondaryWeightByProcess(false);
  fParticleChange->ProposeParentWeight(new_weight);
  

  return fParticleChange;

}

void G4AdjointAlongStepWeightCorrection::BuildPhysicsTable

(

const G4ParticleDefinition &

)

[virtual]

Reimplemented from G4VProcess.

Definition at line 61 of file G4AdjointAlongStepWeightCorrection.cc.

{;
}

G4double G4AdjointAlongStepWeightCorrection::GetContinuousStepLimit	(	const G4Track &	track,
		G4double	previousStepSize,
		G4double	currentMinimumStep,
		G4double &	currentSafety
	)			[protected, virtual]

Implements G4VContinuousProcess.

Definition at line 120 of file G4AdjointAlongStepWeightCorrection.cc.

References DBL_MAX, G4Track::GetKineticEnergy(), and G4Track::GetMaterialCutsCouple().

{ 
  G4double x = DBL_MAX;
  DefineMaterial(track.GetMaterialCutsCouple());
  preStepKinEnergy = track.GetKineticEnergy();
  return x;
}

void G4AdjointAlongStepWeightCorrection::PreparePhysicsTable

(

const G4ParticleDefinition &

)

[virtual]

Reimplemented from G4VProcess.

Definition at line 53 of file G4AdjointAlongStepWeightCorrection.cc.

{
; 
}

---

Field Documentation

G4ParticleChange* G4AdjointAlongStepWeightCorrection::fParticleChange [protected]

Definition at line 112 of file G4AdjointAlongStepWeightCorrection.hh.

Referenced by AlongStepDoIt(), G4AdjointAlongStepWeightCorrection(), and ~G4AdjointAlongStepWeightCorrection().

---

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

• G4AdjointAlongStepWeightCorrection.hh
• G4AdjointAlongStepWeightCorrection.cc

---