#include <Doxymodules_biasing.h>

Inheritance diagram for GB05BOptnSplitAndKillByCrossSection:

Public Member Functions
	GB05BOptnSplitAndKillByCrossSection (G4String name)

virtual	~GB05BOptnSplitAndKillByCrossSection ()

virtual const G4VBiasingInteractionLaw *	ProvideOccurenceBiasingInteractionLaw (const G4BiasingProcessInterface *, G4ForceCondition &) final

virtual G4VParticleChange *	ApplyFinalStateBiasing (const G4BiasingProcessInterface , const G4Track , const G4Step *, G4bool &) final

virtual G4double	DistanceToApplyOperation (const G4Track , G4double, G4ForceCondition condition) final

virtual G4VParticleChange *	GenerateBiasingFinalState (const G4Track , const G4Step ) final

void	SetInteractionLength (G4double interactionLength)

Private Attributes
G4ParticleChange	fParticleChange

G4double	fInteractionLength

Detailed Description

Definition at line 123 of file Doxymodules_biasing.h.

Constructor & Destructor Documentation

◆ GB05BOptnSplitAndKillByCrossSection()

GB05BOptnSplitAndKillByCrossSection::GB05BOptnSplitAndKillByCrossSection ( G4String name )

Definition at line 35 of file GB05BOptnSplitAndKillByCrossSection.cc.

: G4VBiasingOperation(name),
  fParticleChange(),
  fInteractionLength(-1.0)
{}

◆ ~GB05BOptnSplitAndKillByCrossSection()

GB05BOptnSplitAndKillByCrossSection::~GB05BOptnSplitAndKillByCrossSection ( )

virtual

Definition at line 43 of file GB05BOptnSplitAndKillByCrossSection.cc.

44{}

Member Function Documentation

◆ ProvideOccurenceBiasingInteractionLaw()

virtual const G4VBiasingInteractionLaw * GB05BOptnSplitAndKillByCrossSection::ProvideOccurenceBiasingInteractionLaw	(	const G4BiasingProcessInterface *	,
		G4ForceCondition &
	)

inlinefinalvirtual

Definition at line 51 of file GB05BOptnSplitAndKillByCrossSection.hh.

53 {return 0;}

◆ ApplyFinalStateBiasing()

virtual G4VParticleChange * GB05BOptnSplitAndKillByCrossSection::ApplyFinalStateBiasing	(	const G4BiasingProcessInterface *	,
		const G4Track *	,
		const G4Step *	,
		G4bool &
	)

inlinefinalvirtual

Definition at line 55 of file GB05BOptnSplitAndKillByCrossSection.hh.

59 {return 0;}

◆ DistanceToApplyOperation()

G4double GB05BOptnSplitAndKillByCrossSection::DistanceToApplyOperation	(	const G4Track *	,
		G4double	,
		G4ForceCondition *	condition
	)

finalvirtual

Definition at line 48 of file GB05BOptnSplitAndKillByCrossSection.cc.

{
  *condition = NotForced;
 
  // -- Sample the exponential law using the total interaction length of processes
  // -- to couterbalance for:
  G4double proposedStepLength =  -std::log( G4UniformRand() ) * fInteractionLength;
  
  return proposedStepLength;
}

◆ GenerateBiasingFinalState()

G4VParticleChange * GB05BOptnSplitAndKillByCrossSection::GenerateBiasingFinalState	(	const G4Track *	track,
		const G4Step *
	)

finalvirtual

Definition at line 65 of file GB05BOptnSplitAndKillByCrossSection.cc.

{
  
  // -- This method is called if we have limited the step.
  // -- We hence make the splitting or killing.
 
  // Get track weight:
  G4double initialWeight = track->GetWeight();
  
  // The "particle change" is the object to be used to communicate to
  // the tracking the update of the primary state and/or creation
  // secondary tracks.
  fParticleChange.Initialize(*track);
 
  // -- Splitting and killing factors.
  // -- They are taken the same, but the killing factor can be make bigger.
  G4double splittingFactor =  2.0;
  G4double   killingFactor =  2.0;
 
 
  if ( track->GetMomentumDirection().z() > 0 )
    {
      // -- We split if the track is moving forward:
      
      // Define the tracks weight:
      G4double weightOfTrack = initialWeight/splittingFactor;
      
      // Ask currect track weight to be changed to new value:
      fParticleChange.ProposeParentWeight( weightOfTrack );
      // Now make clones of this track (this is the actual splitting):
      // we will then have the primary and clone of it, hence the
      // splitting by a factor 2:
      G4Track* clone = new G4Track( *track );
      clone->SetWeight( weightOfTrack );
      fParticleChange.AddSecondary( clone );
      // -- Below's call added for safety & illustration : inform particle change to not
      // -- modify the clone (ie : daughter) weight to male it that of the
      // -- primary. Here call not mandatory and both tracks have same weights.
      fParticleChange.SetSecondaryWeightByProcess(true);
    }
  else
    {
      // -- We apply Russian roulette if the track is moving backward:
      
      // Shoot a random number (in ]0,1[ segment):
      G4double random = G4UniformRand();
      G4double killingProbability = 1.0 - 1.0/killingFactor;
      if ( random < killingProbability )
        {
          // We ask for the the track to be killed:
          fParticleChange.ProposeTrackStatus(fStopAndKill);
        }
      else
        {
          // In this case, the track survives. We change its weight
          // to conserve weight among killed and survival tracks:
          fParticleChange.ProposeParentWeight( initialWeight*killingFactor );
        }
    }
  
  return &fParticleChange;
  
}

◆ SetInteractionLength()

void GB05BOptnSplitAndKillByCrossSection::SetInteractionLength ( G4double interactionLength )

inline

Definition at line 80 of file GB05BOptnSplitAndKillByCrossSection.hh.

  {
    fInteractionLength = interactionLength;
  }

Member Data Documentation

◆ fParticleChange

G4ParticleChange GB05BOptnSplitAndKillByCrossSection::fParticleChange

private

Definition at line 87 of file GB05BOptnSplitAndKillByCrossSection.hh.

◆ fInteractionLength

G4double GB05BOptnSplitAndKillByCrossSection::fInteractionLength