GB05BOptnSplitAndKillByCrossSection Class Reference

#include <GB05BOptnSplitAndKillByCrossSection.hh>

Inheritance diagram for GB05BOptnSplitAndKillByCrossSection:

Collaboration 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)
Public Member Functions inherited from G4VBiasingOperation
	G4VBiasingOperation (G4String name)
virtual	~G4VBiasingOperation ()
virtual G4double	ProposeAlongStepLimit (const G4BiasingProcessInterface *)
virtual G4GPILSelection	ProposeGPILSelection (const G4GPILSelection wrappedProcessSelection)
virtual void	AlongMoveBy (const G4BiasingProcessInterface *, const G4Step *, G4double)
const G4String &	GetName () const
std::size_t	GetUniqueID () const

Detailed Description

Definition at line 39 of file GB05BOptnSplitAndKillByCrossSection.hh.

Constructor & Destructor Documentation

GB05BOptnSplitAndKillByCrossSection::GB05BOptnSplitAndKillByCrossSection

(

G4String

name

)

Definition at line 36 of file GB05BOptnSplitAndKillByCrossSection.cc.

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

GB05BOptnSplitAndKillByCrossSection::~GB05BOptnSplitAndKillByCrossSection ( )

virtual

Definition at line 44 of file GB05BOptnSplitAndKillByCrossSection.cc.

{}

Member Function Documentation

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

inlinefinalvirtual

Implements G4VBiasingOperation.

Definition at line 56 of file GB05BOptnSplitAndKillByCrossSection.hh.

  {return 0;}

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

finalvirtual

Implements G4VBiasingOperation.

Definition at line 50 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;
}

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

finalvirtual

Implements G4VBiasingOperation.

Definition at line 66 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;
  
}

Here is the call graph for this function:

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

inlinefinalvirtual

Implements G4VBiasingOperation.

Definition at line 52 of file GB05BOptnSplitAndKillByCrossSection.hh.

  {return 0;}

void GB05BOptnSplitAndKillByCrossSection::SetInteractionLength ( G4double  interactionLength )

inline

Definition at line 81 of file GB05BOptnSplitAndKillByCrossSection.hh.

  {
    fInteractionLength = interactionLength;
  }

---

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

• source/geant4.10.03.p03/examples/extended/biasing/GB05/include/GB05BOptnSplitAndKillByCrossSection.hh
• source/geant4.10.03.p03/examples/extended/biasing/GB05/src/GB05BOptnSplitAndKillByCrossSection.cc