G4StoppingPhysics Class Reference

#include <G4StoppingPhysics.hh>

Inheritance diagram for G4StoppingPhysics:

Public Member Functions
	G4StoppingPhysics (G4int ver=1)
	G4StoppingPhysics (const G4String &name, G4int ver=1, G4bool UseMuonMinusCapture=true)
virtual	~G4StoppingPhysics ()
virtual void	ConstructParticle ()
virtual void	ConstructProcess ()
Public Member Functions inherited from G4VPhysicsConstructor
	G4VPhysicsConstructor (const G4String &="")
	G4VPhysicsConstructor (const G4String &name, G4int physics_type)
virtual	~G4VPhysicsConstructor ()
void	SetPhysicsName (const G4String &="")
const G4String &	GetPhysicsName () const
void	SetPhysicsType (G4int)
G4int	GetPhysicsType () const
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const
G4int	GetInstanceID () const

Additional Inherited Members
Static Public Member Functions inherited from G4VPhysicsConstructor
static const G4VPCManager &	GetSubInstanceManager ()
Protected Member Functions inherited from G4VPhysicsConstructor
G4bool	RegisterProcess (G4VProcess *process, G4ParticleDefinition *particle)
Protected Attributes inherited from G4VPhysicsConstructor
G4int	verboseLevel
G4String	namePhysics
G4int	typePhysics
G4ParticleTable *	theParticleTable
G4int	g4vpcInstanceID
Static Protected Attributes inherited from G4VPhysicsConstructor
static G4RUN_DLL G4VPCManager	subInstanceManager

Detailed Description

Definition at line 63 of file G4StoppingPhysics.hh.

Constructor & Destructor Documentation

G4StoppingPhysics::G4StoppingPhysics

(

G4int

ver = 1

)

Definition at line 67 of file G4StoppingPhysics.cc.

References G4cout, and G4endl.

                               :  
  G4VPhysicsConstructor( "stopping" ),
//muProcess( 0 ), hBertiniProcess( 0 ), hFritiofProcess( 0 ),
  verbose( ver ),
  useMuonMinusCapture( true ) 
{
  if ( verbose > 1 ) G4cout << "### G4StoppingPhysics" << G4endl;
}

G4StoppingPhysics::G4StoppingPhysics	(	const G4String &	name,
		G4int	ver = 1,
		G4bool	UseMuonMinusCapture = true
	)

Definition at line 78 of file G4StoppingPhysics.cc.

References G4cout, and G4endl.

                                                                 :
  G4VPhysicsConstructor( name ),
  // muProcess( 0 ), hBertiniProcess( 0 ), hFritiofProcess( 0 ),
  verbose( ver ), 
  useMuonMinusCapture( UseMuonMinusCapture ) 
{
  if ( verbose > 1 ) G4cout << "### G4StoppingPhysics" << G4endl;
}

G4StoppingPhysics::~G4StoppingPhysics ( )

virtual

Definition at line 89 of file G4StoppingPhysics.cc.

{}

Member Function Documentation

void G4StoppingPhysics::ConstructParticle ( void  )

virtual

Implements G4VPhysicsConstructor.

Definition at line 92 of file G4StoppingPhysics.cc.

References G4LeptonConstructor::ConstructParticle(), G4BaryonConstructor::ConstructParticle(), and G4MesonConstructor::ConstructParticle().

                                          {
  // G4cout << "G4StoppingPhysics::ConstructParticle" << G4endl;
  G4LeptonConstructor pLeptonConstructor;
  pLeptonConstructor.ConstructParticle();

  G4MesonConstructor pMesonConstructor;
  pMesonConstructor.ConstructParticle();

  G4BaryonConstructor pBaryonConstructor;
  pBaryonConstructor.ConstructParticle();
}

void G4StoppingPhysics::ConstructProcess ( void  )

virtual

Implements G4VPhysicsConstructor.

Definition at line 105 of file G4StoppingPhysics.cc.

References G4ProcessManager::AddRestProcess(), G4AntiProton::AntiProton(), G4AntiSigmaPlus::AntiSigmaPlus(), aParticleIterator, G4cout, G4endl, G4ParticleDefinition::GetBaryonNumber(), G4ParticleDefinition::GetParticleName(), G4ParticleDefinition::GetPDGCharge(), G4ParticleDefinition::GetPDGMass(), G4ParticleDefinition::GetProcessManager(), G4HadronicAbsorptionBertini::IsApplicable(), G4HadronicAbsorptionFritiof::IsApplicable(), G4ParticleDefinition::IsShortLived(), G4KaonMinus::KaonMinus(), python.hepunit::MeV, G4MuonMinus::MuonMinus(), G4OmegaMinus::OmegaMinus(), G4PionMinus::PionMinus(), G4SigmaMinus::SigmaMinus(), and G4XiMinus::XiMinus().

                                         {
  if ( verbose > 1 ) G4cout << "### G4StoppingPhysics::ConstructProcess " 
                << wasActivated << G4endl;
  if ( wasActivated ) return;
  wasActivated = true;

  G4MuonMinusCapture* muProcess;
  G4HadronicAbsorptionBertini* hBertiniProcess;
  G4HadronicAbsorptionFritiof* hFritiofProcess;

  if ( useMuonMinusCapture ) {
    muProcess = new G4MuonMinusCapture();
  } else {
    muProcess = 0;
  }   

  hBertiniProcess = new G4HadronicAbsorptionBertini();
  hFritiofProcess = new G4HadronicAbsorptionFritiof();

  G4double mThreshold = 130.0*MeV;

  // Add Stopping Process
  G4ParticleDefinition* particle = 0;
  G4ProcessManager* pmanager = 0;

  aParticleIterator->reset();

  while ( (*aParticleIterator)() ) {

    particle = aParticleIterator->value();
    pmanager = particle->GetProcessManager();

    if ( particle == G4MuonMinus::MuonMinus() ) {
      if ( useMuonMinusCapture ) {
     pmanager->AddRestProcess( muProcess );
         if ( verbose > 1 ) {
           G4cout << "### G4StoppingPhysics added G4MuonMinusCapture for " 
              << particle->GetParticleName() << G4endl;
         }
      }
    }

    if ( particle->GetPDGCharge() < 0.0       && 
         particle->GetPDGMass() > mThreshold  &&
         ! particle->IsShortLived() ) {

      // Use Fritiof/Precompound for: anti-protons, anti-sigma+, and
      // anti-nuclei.
      if ( particle == G4AntiProton::AntiProton() ||
           particle == G4AntiSigmaPlus::AntiSigmaPlus() ||
           particle->GetBaryonNumber() < -1 ) {  // Anti-nuclei
        if ( hFritiofProcess->IsApplicable( *particle ) ) {
          pmanager->AddRestProcess( hFritiofProcess );
          if ( verbose > 1 ) {
        G4cout << "### G4HadronicAbsorptionFritiof added for "
                   << particle->GetParticleName() << G4endl;
          }
        }

      // Use Bertini/Precompound for pi-, K-, Sigma-, Xi-, and Omega-
      } else if ( particle == G4PionMinus::PionMinus() ||
                  particle == G4KaonMinus::KaonMinus() ||
                  particle == G4SigmaMinus::SigmaMinus() ||
                  particle == G4XiMinus::XiMinus() ||
                  particle == G4OmegaMinus::OmegaMinus() ) {
        if ( hBertiniProcess->IsApplicable( *particle ) ) {
          pmanager->AddRestProcess( hBertiniProcess );
          if ( verbose > 1 ) {
        G4cout << "### G4HadronicAbsorptionBertini added for "
                   << particle->GetParticleName() << G4endl;
          }
        }

      } else {
        if ( verbose > 1 ) {
          G4cout << "WARNING in G4StoppingPhysics::ConstructProcess: \
                     not able to deal with nuclear stopping of " 
                 << particle->GetParticleName() << G4endl;
        }
      }
    }

  } // end of while loop
}

---

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

• G4StoppingPhysics.hh
• G4StoppingPhysics.cc