G4BiasingHelper Class Reference

#include <G4BiasingHelper.hh>

Collaboration diagram for G4BiasingHelper:

Static Public Member Functions
static G4bool	ActivatePhysicsBiasing (G4ProcessManager *pmanager, G4String physicsProcessToBias, G4String wrappedName="")
static void	ActivateNonPhysicsBiasing (G4ProcessManager *pmanager, G4String nonPhysicsProcessName="")

Detailed Description

Definition at line 52 of file G4BiasingHelper.hh.

Member Function Documentation

ActivateNonPhysicsBiasing()

void G4BiasingHelper::ActivateNonPhysicsBiasing ( G4ProcessManager *  pmanager, G4String  nonPhysicsProcessName = ""  )

static

Definition at line 94 of file G4BiasingHelper.cc.

{
  G4BiasingProcessInterface* biasingNonPhys(0);
  if ( nonPhysicsProcessName == "" ) biasingNonPhys = new G4BiasingProcessInterface();
  else                               biasingNonPhys = new G4BiasingProcessInterface(nonPhysicsProcessName );
  pmanager->AddProcess( biasingNonPhys,
            ordInActive,
            ordInActive,
            ordDefault);
}

Here is the call graph for this function:

Here is the caller graph for this function:

ActivatePhysicsBiasing()

G4bool G4BiasingHelper::ActivatePhysicsBiasing ( G4ProcessManager *  pmanager, G4String  physicsProcessToBias, G4String  wrappedName = ""  )

static

Definition at line 31 of file G4BiasingHelper.cc.

{
  G4VProcess* physicsProcess(0);
  
  G4ProcessVector* vprocess = pmanager->GetProcessList();
  for (G4int ip = 0 ; ip < vprocess->size() ; ip++)
    {
      if ( (*vprocess)[ip]->GetProcessName() == physicsProcessToBias )
    {
      physicsProcess = (*vprocess)[ip];
      break;
    }
    }
  
  // -- process not found, return "false" to tell about failure
  if ( physicsProcess == 0 ) return false;
  
  // -- process is not a physics one, return "false" to tell about failure
  G4int processType = physicsProcess->GetProcessType();
  if ( ( processType != 2 ) &&  // EM
       ( processType != 3 ) &&  // Optical
       ( processType != 4 ) &&  // Hadronic
       ( processType != 6 ) )   // Decay
    return false;

  // -- prevent wrapper of wrapper...
  if ( dynamic_cast< G4BiasingProcessInterface* >( physicsProcess ) ) return false;

  // -- remember process indeces:
  G4int    atRestIndex = pmanager->GetProcessOrdering(physicsProcess, idxAtRest   );
  G4int alongStepIndex = pmanager->GetProcessOrdering(physicsProcess, idxAlongStep);
  G4int  postStepIndex = pmanager->GetProcessOrdering(physicsProcess, idxPostStep );

  // -- now remove the physic process, that will be replaced by a wrapped version:
  G4VProcess* removed = pmanager->RemoveProcess(physicsProcess);
  if ( removed != physicsProcess )
    {
      G4ExceptionDescription ed;
      ed << "Internal inconsistency in processes handling. Please report !" << G4endl;
      G4Exception("G4BiasingHelper::ActivatePhysicsBiasing(...)",
          "BIAS.GEN.01",
          FatalException,
          ed);
    }
  
  G4BiasingProcessInterface* biasingWrapper = new G4BiasingProcessInterface( physicsProcess,
                                         atRestIndex    != ordInActive,
                                         alongStepIndex != ordInActive,
                                         postStepIndex  != ordInActive,
                                         wrappedName);

  if ( alongStepIndex == -1 ) alongStepIndex = ordDefault;
  
  pmanager->AddProcess( biasingWrapper,
            atRestIndex,
            alongStepIndex,
            postStepIndex);

  return true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

---

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

• G4BiasingHelper.hh
• G4BiasingHelper.cc