G4BOptnForceCommonTruncatedExp Class Reference

#include <G4BOptnForceCommonTruncatedExp.hh>

Inheritance diagram for G4BOptnForceCommonTruncatedExp:

Collaboration diagram for G4BOptnForceCommonTruncatedExp:

Public Member Functions
	G4BOptnForceCommonTruncatedExp (G4String name)
virtual	~G4BOptnForceCommonTruncatedExp ()
virtual const G4VBiasingInteractionLaw *	ProvideOccurenceBiasingInteractionLaw (const G4BiasingProcessInterface *, G4ForceCondition &)
virtual G4double	ProposeAlongStepLimit (const G4BiasingProcessInterface *)
virtual G4GPILSelection	ProposeGPILSelection (const G4GPILSelection processSelection)
virtual G4VParticleChange *	ApplyFinalStateBiasing (const G4BiasingProcessInterface *, const G4Track *, const G4Step *, G4bool &)
virtual G4double	DistanceToApplyOperation (const G4Track *, G4double, G4ForceCondition *)
virtual G4VParticleChange *	GenerateBiasingFinalState (const G4Track *, const G4Step *)
G4ILawCommonTruncatedExp *	GetCommonTruncatedExpLaw ()
G4ILawForceFreeFlight *	GetForceFreeFlightLaw ()
void	Initialize (const G4Track *)
void	UpdateForStep (const G4Step *)
void	Sample ()
const G4ThreeVector &	GetInitialMomentum () const
G4double	GetMaximumDistance () const
void	ChooseProcessToApply ()
const G4VProcess *	GetProcessToApply () const
void	AddCrossSection (const G4VProcess *, G4double)
size_t	GetNumberOfSharing () const
void	SetInteractionOccured (G4bool b)
G4bool	GetInteractionOccured () const
Public Member Functions inherited from G4VBiasingOperation
	G4VBiasingOperation (G4String name)
virtual	~G4VBiasingOperation ()
virtual void	AlongMoveBy (const G4BiasingProcessInterface *, const G4Step *, G4double)
const G4String &	GetName () const
std::size_t	GetUniqueID () const

Detailed Description

Definition at line 60 of file G4BOptnForceCommonTruncatedExp.hh.

Constructor & Destructor Documentation

G4BOptnForceCommonTruncatedExp::G4BOptnForceCommonTruncatedExp

(

G4String

name

)

Definition at line 34 of file G4BOptnForceCommonTruncatedExp.cc.

  : G4VBiasingOperation(name),
    fNumberOfSharing(0),
    fProcessToApply(nullptr),
    fInteractionOccured(false),
    fMaximumDistance(-1.0)
{
  fCommonTruncatedExpLaw = new G4ILawCommonTruncatedExp("ExpLawForOperation"+name);
  fForceFreeFlightLaw    = new G4ILawForceFreeFlight   ("FFFLawForOperation"+name);
  
  fTotalCrossSection = 0.0;
}

G4BOptnForceCommonTruncatedExp::~G4BOptnForceCommonTruncatedExp ( )

virtual

Definition at line 47 of file G4BOptnForceCommonTruncatedExp.cc.

{
  if ( fCommonTruncatedExpLaw ) delete fCommonTruncatedExpLaw;
  if ( fForceFreeFlightLaw )    delete fForceFreeFlightLaw;
}

Member Function Documentation

void G4BOptnForceCommonTruncatedExp::AddCrossSection	(	const G4VProcess *	process,
		G4double	crossSection
	)

Definition at line 116 of file G4BOptnForceCommonTruncatedExp.cc.

{
  fTotalCrossSection      += crossSection;
  fCrossSections[process]  = crossSection;
  fNumberOfSharing         = fCrossSections.size();
}

G4VParticleChange * G4BOptnForceCommonTruncatedExp::ApplyFinalStateBiasing ( const G4BiasingProcessInterface *  callingProcess, const G4Track *  track, const G4Step *  step, G4bool &  forceFinalState  )

virtual

Implements G4VBiasingOperation.

Definition at line 76 of file G4BOptnForceCommonTruncatedExp.cc.

{
  if ( callingProcess->GetWrappedProcess() != fProcessToApply )
    {
      forceFinalState = true;
      fDummyParticleChange.Initialize( *track );
      return &fDummyParticleChange; 
    }
  if ( fInteractionOccured )
    {
      forceFinalState = true;
      fDummyParticleChange.Initialize( *track );
      return &fDummyParticleChange;
    }
  
  // -- checks if process won the GPIL race:
  G4double processGPIL = callingProcess->GetPostStepGPIL() < callingProcess->GetAlongStepGPIL() ?
    callingProcess->GetPostStepGPIL() : callingProcess->GetAlongStepGPIL() ;
  if ( processGPIL <= step->GetStepLength() )
    {
      // -- if process won, wrapped process produces the final state.
      // -- In this case, the weight for occurence biasing is applied
      // -- by the callingProcess, at exit of present method. This is
      // -- selected by "forceFinalState = false":
      forceFinalState     = false;
      fInteractionOccured = true;
      return callingProcess->GetWrappedProcess()->PostStepDoIt( *track, *step );
    }
  else
    {
      forceFinalState = true;
      fDummyParticleChange.Initialize( *track );
      return &fDummyParticleChange; 
    }
}

Here is the call graph for this function:

void G4BOptnForceCommonTruncatedExp::ChooseProcessToApply

(

)

Definition at line 167 of file G4BOptnForceCommonTruncatedExp.cc.

{
  G4double sigmaRand   = G4UniformRand() * fTotalCrossSection;
  G4double sigmaSelect = 0.0;
  for ( std::map< const G4VProcess*, G4double>::const_iterator it = fCrossSections.begin();
    it != fCrossSections.end();
    it++)
    {
      sigmaSelect += (*it).second;
      if ( sigmaRand <= sigmaSelect )
    {
      fProcessToApply = (*it).first;
      break;
    }
    }
}

Here is the caller graph for this function:

virtual G4double G4BOptnForceCommonTruncatedExp::DistanceToApplyOperation ( const G4Track *  , G4double  , G4ForceCondition *    )

inlinevirtual

Implements G4VBiasingOperation.

Definition at line 79 of file G4BOptnForceCommonTruncatedExp.hh.

                                                           {return DBL_MAX;}

virtual G4VParticleChange* G4BOptnForceCommonTruncatedExp::GenerateBiasingFinalState ( const G4Track *  , const G4Step *    )

inlinevirtual

Implements G4VBiasingOperation.

Definition at line 82 of file G4BOptnForceCommonTruncatedExp.hh.

                                                           {return 0;}

G4ILawCommonTruncatedExp* G4BOptnForceCommonTruncatedExp::GetCommonTruncatedExpLaw ( )

inline

Definition at line 89 of file G4BOptnForceCommonTruncatedExp.hh.

  {
    return fCommonTruncatedExpLaw;
  }

G4ILawForceFreeFlight* G4BOptnForceCommonTruncatedExp::GetForceFreeFlightLaw ( )

inline

Definition at line 93 of file G4BOptnForceCommonTruncatedExp.hh.

  {
    return fForceFreeFlightLaw;
  }

const G4ThreeVector& G4BOptnForceCommonTruncatedExp::GetInitialMomentum ( ) const

inline

Definition at line 101 of file G4BOptnForceCommonTruncatedExp.hh.

{ return fInitialMomentum; }

G4bool G4BOptnForceCommonTruncatedExp::GetInteractionOccured ( ) const

inline

Definition at line 108 of file G4BOptnForceCommonTruncatedExp.hh.

{ return fInteractionOccured; }

Here is the caller graph for this function:

G4double G4BOptnForceCommonTruncatedExp::GetMaximumDistance ( ) const

inline

Definition at line 102 of file G4BOptnForceCommonTruncatedExp.hh.

{ return fMaximumDistance; }

size_t G4BOptnForceCommonTruncatedExp::GetNumberOfSharing ( ) const

inline

Definition at line 106 of file G4BOptnForceCommonTruncatedExp.hh.

{ return fNumberOfSharing;}

const G4VProcess* G4BOptnForceCommonTruncatedExp::GetProcessToApply ( ) const

inline

Definition at line 104 of file G4BOptnForceCommonTruncatedExp.hh.

{ return fProcessToApply; }

void G4BOptnForceCommonTruncatedExp::Initialize

(

const G4Track *

track

)

Definition at line 124 of file G4BOptnForceCommonTruncatedExp.cc.

{
  fCrossSections.clear();
  fTotalCrossSection  = 0.0;
  fNumberOfSharing    = 0;
  fProcessToApply     = 0;
  fInteractionOccured = false;
  fInitialMomentum    = track->GetMomentum();

  G4VSolid* currentSolid = track->GetVolume()->GetLogicalVolume()->GetSolid();
  G4ThreeVector  localPosition = (G4TransportationManager::GetTransportationManager()->
                  GetNavigatorForTracking()->
                  GetGlobalToLocalTransform()).TransformPoint(track->GetPosition());
  G4ThreeVector localDirection = (G4TransportationManager::GetTransportationManager()->
                  GetNavigatorForTracking()->
                  GetGlobalToLocalTransform()).TransformAxis(track->GetMomentumDirection());
  fMaximumDistance = currentSolid->DistanceToOut(localPosition, localDirection);
  if ( fMaximumDistance <= DBL_MIN )  fMaximumDistance = 0.0;
  fCommonTruncatedExpLaw->SetMaximumDistance( fMaximumDistance );
}

Here is the call graph for this function:

virtual G4double G4BOptnForceCommonTruncatedExp::ProposeAlongStepLimit ( const G4BiasingProcessInterface *  )

inlinevirtual

Reimplemented from G4VBiasingOperation.

Definition at line 72 of file G4BOptnForceCommonTruncatedExp.hh.

{ return DBL_MAX; }

G4GPILSelection G4BOptnForceCommonTruncatedExp::ProposeGPILSelection ( const G4GPILSelection  processSelection )

virtual

Reimplemented from G4VBiasingOperation.

Definition at line 70 of file G4BOptnForceCommonTruncatedExp.cc.

{
  return NotCandidateForSelection;
}

const G4VBiasingInteractionLaw * G4BOptnForceCommonTruncatedExp::ProvideOccurenceBiasingInteractionLaw ( const G4BiasingProcessInterface *  callingProcess, G4ForceCondition &  proposeForceCondition  )

virtual

Implements G4VBiasingOperation.

Definition at line 54 of file G4BOptnForceCommonTruncatedExp.cc.

{
  if ( callingProcess->GetWrappedProcess() == fProcessToApply )
    {
      proposeForceCondition =                 Forced;
      return                  fCommonTruncatedExpLaw;
    }
  else
    {
      proposeForceCondition =                 Forced;
      return                     fForceFreeFlightLaw;
    }
}

Here is the call graph for this function:

void G4BOptnForceCommonTruncatedExp::Sample

(

)

Definition at line 158 of file G4BOptnForceCommonTruncatedExp.cc.

{
  fCommonTruncatedExpLaw->SetForceCrossSection( fTotalCrossSection );
  fCommonTruncatedExpLaw->Sample();
  ChooseProcessToApply();
  fCommonTruncatedExpLaw->SetSelectedProcessXSfraction(fCrossSections[fProcessToApply] / fTotalCrossSection);
}

Here is the call graph for this function:

void G4BOptnForceCommonTruncatedExp::SetInteractionOccured ( G4bool  b )

inline

Definition at line 107 of file G4BOptnForceCommonTruncatedExp.hh.

{ fInteractionOccured = b;    }

void G4BOptnForceCommonTruncatedExp::UpdateForStep

(

const G4Step *

step

)

Definition at line 146 of file G4BOptnForceCommonTruncatedExp.cc.

{
  fCrossSections.clear();
  fTotalCrossSection  = 0.0;
  fNumberOfSharing    = 0;
  fProcessToApply     = 0;
  
  fCommonTruncatedExpLaw->UpdateForStep( step->GetStepLength() );
  fMaximumDistance = fCommonTruncatedExpLaw->GetMaximumDistance();
}

Here is the call graph for this function:

---

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

• geant4.10.03.p01/source/processes/biasing/generic/include/G4BOptnForceCommonTruncatedExp.hh
• geant4.10.03.p01/source/processes/biasing/generic/src/G4BOptnForceCommonTruncatedExp.cc