#include <Doxymodules_biasing.h>

Inheritance diagram for GB07BOptrLeadingParticle:

Public Member Functions
	GB07BOptrLeadingParticle (G4String operatorName="LeadingParticleBiasingOperator")

virtual	~GB07BOptrLeadingParticle ()

virtual void	StartRun () final

virtual void	StartTracking (const G4Track *track) final

Private Member Functions
virtual G4VBiasingOperation *	ProposeNonPhysicsBiasingOperation (const G4Track , const G4BiasingProcessInterface ) final

virtual G4VBiasingOperation *	ProposeOccurenceBiasingOperation (const G4Track , const G4BiasingProcessInterface ) final

virtual G4VBiasingOperation *	ProposeFinalStateBiasingOperation (const G4Track track, const G4BiasingProcessInterface callingProcess) final

Private Attributes
G4BOptnLeadingParticle *	fLeadingParticleBiasingOperation

const G4VProcess *	fAnnihilation

const G4VProcess *	fConversion

const G4VProcess *	fDecay

const G4VProcess *	fTwoParticleProcess

Detailed Description

Definition at line 154 of file Doxymodules_biasing.h.

Constructor & Destructor Documentation

◆ GB07BOptrLeadingParticle()

GB07BOptrLeadingParticle::GB07BOptrLeadingParticle ( G4String operatorName = "LeadingParticleBiasingOperator" )

Definition at line 43 of file GB07BOptrLeadingParticle.cc.

  : G4VBiasingOperator        ( operatorName ),
    fAnnihilation             ( nullptr ),
    fConversion               ( nullptr ),
    fDecay                    ( nullptr ),
    fTwoParticleProcess       ( nullptr )
{
  fLeadingParticleBiasingOperation =
    new G4BOptnLeadingParticle("LeadingParticleBiasingOperation");
}

◆ ~GB07BOptrLeadingParticle()

GB07BOptrLeadingParticle::~GB07BOptrLeadingParticle ( )

virtual

Definition at line 56 of file GB07BOptrLeadingParticle.cc.

{
  delete fLeadingParticleBiasingOperation;
}

Member Function Documentation

◆ ProposeNonPhysicsBiasingOperation()

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

inlinefinalprivatevirtual

Definition at line 49 of file GB07BOptrLeadingParticle.hh.

51 { return nullptr; }

◆ ProposeOccurenceBiasingOperation()

virtual G4VBiasingOperation * GB07BOptrLeadingParticle::ProposeOccurenceBiasingOperation	(	const G4Track *	,
		const G4BiasingProcessInterface *
	)

inlinefinalprivatevirtual

Definition at line 54 of file GB07BOptrLeadingParticle.hh.

56 { return nullptr; }

◆ ProposeFinalStateBiasingOperation()

G4VBiasingOperation * GB07BOptrLeadingParticle::ProposeFinalStateBiasingOperation	(	const G4Track *	track,
		const G4BiasingProcessInterface *	callingProcess
	)

finalprivatevirtual

Definition at line 64 of file GB07BOptrLeadingParticle.cc.

{
  // --   When the present method is called, we are at the process final state
  // -- generation level. The process is given by the callingProcess argument,
  // -- which, in our case,  wrappes a physics process, to control it.
  // --   To bias the final state generation, we return a biasing operation
  // -- which is fLeadingParticleBiasingOperation here. Before returning it, we
  // -- configure it depending on if the process is a two-particle final state
  // -- or if it is a many-particle final state process. For the two-particle
  // -- final state, one track is the leading and the other is alone in its category,
  // -- so always surviving by default. We play a Russian roulette on it to
  // -- trim also these two-particles final states.
  
  if ( callingProcess == fTwoParticleProcess )
    {
      // -- secondary particle accompagnying the leading one will be
      // -- killed with 2./3. probability (Russian roulette):
      fLeadingParticleBiasingOperation->SetFurtherKillingProbability( 2./3.);
    }
  else
    {
      // -- -1.0 means no effect : no further killing is applied to secondary
      // -- particles accompanying the leading one.
      fLeadingParticleBiasingOperation->SetFurtherKillingProbability( -1.0 );
    }
  
  return fLeadingParticleBiasingOperation;
}

◆ StartRun()

void GB07BOptrLeadingParticle::StartRun ( )

finalvirtual

Definition at line 98 of file GB07BOptrLeadingParticle.cc.

{
  // -- collect the two-particle final state processes:
  fAnnihilation = nullptr;
  fConversion   = nullptr;
  fDecay        = nullptr;
  
  // ---- collect e+ annihilation process:
  auto positronProcesses = G4Positron::Definition()->GetProcessManager()->GetProcessList();
  for ( size_t i = 0; i < positronProcesses->size(); ++i )
    {
      if ( (*positronProcesses)[i]->GetProcessName() == "biasWrapper(annihil)")
        {
          fAnnihilation = (*positronProcesses)[i];
          break;
        }
    }
  
  // ---- collect gamma conversion process:
  auto gammaProcesses = G4Gamma::Definition()->GetProcessManager()->GetProcessList();
  for ( size_t i = 0; i < gammaProcesses->size(); ++i )
    {
      if ( (*gammaProcesses)[i]->GetProcessName() == "biasWrapper(conv)")
        {
          fConversion = (*gammaProcesses)[i];
          break;
        }
    }
  
  // ---- collect pi0 decay process:
  auto pi0Processes = G4PionZero::Definition()->GetProcessManager()->GetProcessList();
  for ( size_t i = 0; i < pi0Processes->size(); ++i )
    {
      if ( (*pi0Processes)[i]->GetProcessName() == "biasWrapper(Decay)")
        {
          fDecay = (*pi0Processes)[i];
          break;
        }
    }
}

◆ StartTracking()

void GB07BOptrLeadingParticle::StartTracking ( const G4Track * track )

finalvirtual

Definition at line 143 of file GB07BOptrLeadingParticle.cc.

{
  // -- remember what is the two-particle final state process -if any- for this starting
  // -- track:
  fTwoParticleProcess = nullptr;
  if ( track->GetDefinition() == G4Gamma   ::Definition() ) fTwoParticleProcess =   fConversion;
  if ( track->GetDefinition() == G4Positron::Definition() ) fTwoParticleProcess = fAnnihilation;
  if ( track->GetDefinition() == G4PionZero::Definition() ) fTwoParticleProcess =        fDecay;
}