GB01BOptrChangeCrossSection Class Reference

#include <GB01BOptrChangeCrossSection.hh>

Inheritance diagram for GB01BOptrChangeCrossSection:

Collaboration diagram for GB01BOptrChangeCrossSection:

Public Member Functions
	GB01BOptrChangeCrossSection (G4String particleToBias, G4String name="ChangeXS")
virtual	~GB01BOptrChangeCrossSection ()
virtual void	StartRun ()
Public Member Functions inherited from G4VBiasingOperator
	G4VBiasingOperator (G4String name)
virtual	~G4VBiasingOperator ()
virtual void	Configure ()
virtual void	ConfigureForWorker ()
virtual void	StartTracking (const G4Track *)
virtual void	EndTracking ()
const G4String	GetName () const
void	AttachTo (const G4LogicalVolume *)
G4BiasingAppliedCase	GetPreviousBiasingAppliedCase () const
G4VBiasingOperation *	GetProposedOccurenceBiasingOperation (const G4Track *track, const G4BiasingProcessInterface *callingProcess)
G4VBiasingOperation *	GetProposedFinalStateBiasingOperation (const G4Track *track, const G4BiasingProcessInterface *callingProcess)
G4VBiasingOperation *	GetProposedNonPhysicsBiasingOperation (const G4Track *track, const G4BiasingProcessInterface *callingProcess)
void	ExitingBiasing (const G4Track *track, const G4BiasingProcessInterface *callingProcess)
void	ReportOperationApplied (const G4BiasingProcessInterface *callingProcess, G4BiasingAppliedCase biasingCase, G4VBiasingOperation *operationApplied, const G4VParticleChange *particleChangeProduced)
void	ReportOperationApplied (const G4BiasingProcessInterface *callingProcess, G4BiasingAppliedCase biasingCase, G4VBiasingOperation *occurenceOperationApplied, G4double weightForOccurenceInteraction, G4VBiasingOperation *finalStateOperationApplied, const G4VParticleChange *particleChangeProduced)
const G4VBiasingOperation *	GetPreviousNonPhysicsAppliedOperation ()

Private Member Functions
virtual G4VBiasingOperation *	ProposeOccurenceBiasingOperation (const G4Track *track, const G4BiasingProcessInterface *callingProcess)
virtual G4VBiasingOperation *	ProposeFinalStateBiasingOperation (const G4Track *, const G4BiasingProcessInterface *)
virtual G4VBiasingOperation *	ProposeNonPhysicsBiasingOperation (const G4Track *, const G4BiasingProcessInterface *)
virtual void	OperationApplied (const G4BiasingProcessInterface *callingProcess, G4BiasingAppliedCase biasingCase, G4VBiasingOperation *occurenceOperationApplied, G4double weightForOccurenceInteraction, G4VBiasingOperation *finalStateOperationApplied, const G4VParticleChange *particleChangeProduced)

Private Attributes
std::map< const G4BiasingProcessInterface *, G4BOptnChangeCrossSection *>	fChangeCrossSectionOperations
G4bool	fSetup
const G4ParticleDefinition *	fParticleToBias

Additional Inherited Members
Static Public Member Functions inherited from G4VBiasingOperator
static const std::vector< G4VBiasingOperator *> &	GetBiasingOperators ()
static G4VBiasingOperator *	GetBiasingOperator (const G4LogicalVolume *)
Protected Member Functions inherited from G4VBiasingOperator
virtual void	OperationApplied (const G4BiasingProcessInterface *callingProcess, G4BiasingAppliedCase biasingCase, G4VBiasingOperation *operationApplied, const G4VParticleChange *particleChangeProduced)
virtual void	ExitBiasing (const G4Track *track, const G4BiasingProcessInterface *callingProcess)

Detailed Description

Definition at line 53 of file GB01BOptrChangeCrossSection.hh.

Constructor & Destructor Documentation

GB01BOptrChangeCrossSection()

GB01BOptrChangeCrossSection::GB01BOptrChangeCrossSection	(	G4String	particleToBias,
		G4String	name = "ChangeXS"
	)

Definition at line 38 of file GB01BOptrChangeCrossSection.cc.

  : G4VBiasingOperator(name),
    fSetup(true)
{
  fParticleToBias = G4ParticleTable::GetParticleTable()->FindParticle(particleName);
  
  if ( fParticleToBias == 0 )
    {
      G4ExceptionDescription ed;
      ed << "Particle `" << particleName << "' not found !" << G4endl;
      G4Exception("GB01BOptrChangeCrossSection(...)",
                  "exGB01.01",
                  JustWarning,
                  ed);
    }
}

Here is the call graph for this function:

~GB01BOptrChangeCrossSection()

GB01BOptrChangeCrossSection::~GB01BOptrChangeCrossSection ( )

virtual

Definition at line 56 of file GB01BOptrChangeCrossSection.cc.

{
  for ( std::map< const G4BiasingProcessInterface*, G4BOptnChangeCrossSection* >::iterator 
          it = fChangeCrossSectionOperations.begin() ;
        it != fChangeCrossSectionOperations.end() ;
        it++ ) delete (*it).second;
}

Member Function Documentation

OperationApplied()

void GB01BOptrChangeCrossSection::OperationApplied ( const G4BiasingProcessInterface *  callingProcess, G4BiasingAppliedCase  biasingCase, G4VBiasingOperation *  occurenceOperationApplied, G4double  weightForOccurenceInteraction, G4VBiasingOperation *  finalStateOperationApplied, const G4VParticleChange *  particleChangeProduced  )

privatevirtual

Reimplemented from G4VBiasingOperator.

Definition at line 183 of file GB01BOptrChangeCrossSection.cc.

{
  G4BOptnChangeCrossSection* operation = fChangeCrossSectionOperations[callingProcess];
  if ( operation ==  occurenceOperationApplied ) operation->SetInteractionOccured();
}

Here is the call graph for this function:

Here is the caller graph for this function:

ProposeFinalStateBiasingOperation()

virtual G4VBiasingOperation* GB01BOptrChangeCrossSection::ProposeFinalStateBiasingOperation ( const G4Track *  , const G4BiasingProcessInterface *    )

inlineprivatevirtual

Implements G4VBiasingOperator.

Definition at line 74 of file GB01BOptrChangeCrossSection.hh.

  {return 0;}

ProposeNonPhysicsBiasingOperation()

virtual G4VBiasingOperation* GB01BOptrChangeCrossSection::ProposeNonPhysicsBiasingOperation ( const G4Track *  , const G4BiasingProcessInterface *    )

inlineprivatevirtual

Implements G4VBiasingOperator.

Definition at line 77 of file GB01BOptrChangeCrossSection.hh.

  {return 0;}

Here is the call graph for this function:

ProposeOccurenceBiasingOperation()

G4VBiasingOperation * GB01BOptrChangeCrossSection::ProposeOccurenceBiasingOperation ( const G4Track *  track, const G4BiasingProcessInterface *  callingProcess  )

privatevirtual

Implements G4VBiasingOperator.

Definition at line 96 of file GB01BOptrChangeCrossSection.cc.

{

  // -----------------------------------------------------
  // -- Check if current particle type is the one to bias:
  // -----------------------------------------------------
  if ( track->GetDefinition() != fParticleToBias ) return 0;
  
  
  // ---------------------------------------------------------------------
  // -- select and setup the biasing operation for current callingProcess:
  // ---------------------------------------------------------------------
  // -- Check if the analog cross-section well defined : for example, the conversion
  // -- process for a gamma below e+e- creation threshold has an DBL_MAX interaction
  // -- length. Nothing is done in this case (ie, let analog process to deal with the case)
  G4double analogInteractionLength =  
    callingProcess->GetWrappedProcess()->GetCurrentInteractionLength();
  if ( analogInteractionLength > DBL_MAX/10. ) return 0;

  // -- Analog cross-section is well-defined:
  G4double analogXS = 1./analogInteractionLength;

  // -- Choose a constant cross-section bias. But at this level, this factor can be made
  // -- direction dependent, like in the exponential transform MCNP case, or it
  // -- can be chosen differently, depending on the process, etc.
  G4double XStransformation = 2.0 ;
  
  // -- fetch the operation associated to this callingProcess:
  G4BOptnChangeCrossSection*   operation = fChangeCrossSectionOperations[callingProcess];
  // -- get the operation that was proposed to the process in the previous step:
  G4VBiasingOperation* previousOperation = callingProcess->GetPreviousOccurenceBiasingOperation();
  
  // -- now setup the operation to be returned to the process: this
  // -- consists in setting the biased cross-section, and in asking
  // -- the operation to sample its exponential interaction law.
  // -- To do this, to first order, the two lines:
  //        operation->SetBiasedCrossSection( XStransformation * analogXS );
  //        operation->Sample();
  // -- are correct and sufficient.
  // -- But, to avoid having to shoot a random number each time, we sample
  // -- only on the first time the operation is proposed, or if the interaction
  // -- occured. If the interaction did not occur for the process in the previous,
  // -- we update the number of interaction length instead of resampling.
  if ( previousOperation == 0 )
    {
      operation->SetBiasedCrossSection( XStransformation * analogXS );
      operation->Sample();
    }
  else
    {
      if (  previousOperation != operation )
        {
          // -- should not happen !
          G4ExceptionDescription ed;
          ed << " Logic problem in operation handling !" << G4endl;
          G4Exception("GB01BOptrChangeCrossSection::ProposeOccurenceBiasingOperation(...)",
                      "exGB01.02",
                      JustWarning,
                      ed);
          return 0;
        }
      if ( operation->GetInteractionOccured() )
        {
          operation->SetBiasedCrossSection( XStransformation * analogXS );
          operation->Sample();
        }
      else
        {
          // -- update the 'interaction length' and underneath 'number of interaction lengths'
          // -- for past step  (this takes into accout the previous step cross-section value)
          operation->UpdateForStep( callingProcess->GetPreviousStepSize() );
          // -- update the cross-section value:
          operation->SetBiasedCrossSection( XStransformation * analogXS );
          // -- forces recomputation of the 'interaction length' taking into account above
          // -- new cross-section value [tricky : to be improved]
          operation->UpdateForStep( 0.0 );
        }
    }
  
  return operation;
  
}

Here is the call graph for this function:

StartRun()

void GB01BOptrChangeCrossSection::StartRun ( )

virtual

Reimplemented from G4VBiasingOperator.

Definition at line 65 of file GB01BOptrChangeCrossSection.cc.

{
  // --------------
  // -- Setup stage:
  // ---------------
  // -- Start by collecting processes under biasing, create needed biasing
  // -- operations and associate these operations to the processes:
  if ( fSetup )
    {
      const G4ProcessManager* processManager = fParticleToBias->GetProcessManager();
      const G4BiasingProcessSharedData* sharedData =
        G4BiasingProcessInterface::GetSharedData( processManager );
      if ( sharedData ) // -- sharedData tested, as is can happen a user attaches an operator to a
                        // -- volume but without defined BiasingProcessInterface processes.
        {
          for ( size_t i = 0 ; i < (sharedData->GetPhysicsBiasingProcessInterfaces()).size(); i++ )
            {
              const G4BiasingProcessInterface* wrapperProcess =
                (sharedData->GetPhysicsBiasingProcessInterfaces())[i];
              G4String operationName = "XSchange-" + 
                wrapperProcess->GetWrappedProcess()->GetProcessName();
              fChangeCrossSectionOperations[wrapperProcess] = 
                new G4BOptnChangeCrossSection(operationName);
            }
        }
      fSetup = false;
    }
}

Here is the call graph for this function:

Member Data Documentation

fChangeCrossSectionOperations

std::map< const G4BiasingProcessInterface*, G4BOptnChangeCrossSection* > GB01BOptrChangeCrossSection::fChangeCrossSectionOperations

private

Definition at line 98 of file GB01BOptrChangeCrossSection.hh.

fParticleToBias

const G4ParticleDefinition* GB01BOptrChangeCrossSection::fParticleToBias

private

Definition at line 100 of file GB01BOptrChangeCrossSection.hh.

fSetup

G4bool GB01BOptrChangeCrossSection::fSetup

private

Definition at line 99 of file GB01BOptrChangeCrossSection.hh.

---

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

• GB01BOptrChangeCrossSection.hh
• GB01BOptrChangeCrossSection.cc