G4LowEIonFragmentation Class Reference

#include <G4LowEIonFragmentation.hh>

Inheritance diagram for G4LowEIonFragmentation:

Public Member Functions
	G4LowEIonFragmentation (G4ExcitationHandler *const value=nullptr)
virtual	~G4LowEIonFragmentation ()
G4HadFinalState *	ApplyYourself (const G4HadProjectile &thePrimary, G4Nucleus &theNucleus) override
G4double	GetCrossSection ()
	G4LowEIonFragmentation (const G4LowEIonFragmentation &)=delete
const G4LowEIonFragmentation &	operator= (const G4LowEIonFragmentation &right)=delete
G4bool	operator== (const G4LowEIonFragmentation &right) const =delete
G4bool	operator!= (const G4LowEIonFragmentation &right) const =delete
Public Member Functions inherited from G4HadronicInteraction
	G4HadronicInteraction (const G4String &modelName="HadronicModel")
virtual	~G4HadronicInteraction ()
virtual G4double	SampleInvariantT (const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A)
virtual G4bool	IsApplicable (const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
G4double	GetMinEnergy () const
G4double	GetMinEnergy (const G4Material *aMaterial, const G4Element *anElement) const
void	SetMinEnergy (G4double anEnergy)
void	SetMinEnergy (G4double anEnergy, const G4Element *anElement)
void	SetMinEnergy (G4double anEnergy, const G4Material *aMaterial)
G4double	GetMaxEnergy () const
G4double	GetMaxEnergy (const G4Material *aMaterial, const G4Element *anElement) const
void	SetMaxEnergy (const G4double anEnergy)
void	SetMaxEnergy (G4double anEnergy, const G4Element *anElement)
void	SetMaxEnergy (G4double anEnergy, const G4Material *aMaterial)
G4int	GetVerboseLevel () const
void	SetVerboseLevel (G4int value)
const G4String &	GetModelName () const
void	DeActivateFor (const G4Material *aMaterial)
void	ActivateFor (const G4Material *aMaterial)
void	DeActivateFor (const G4Element *anElement)
void	ActivateFor (const G4Element *anElement)
G4bool	IsBlocked (const G4Material *aMaterial) const
G4bool	IsBlocked (const G4Element *anElement) const
void	SetRecoilEnergyThreshold (G4double val)
G4double	GetRecoilEnergyThreshold () const
virtual const std::pair< G4double, G4double >	GetFatalEnergyCheckLevels () const
virtual std::pair< G4double, G4double >	GetEnergyMomentumCheckLevels () const
void	SetEnergyMomentumCheckLevels (G4double relativeLevel, G4double absoluteLevel)
virtual void	ModelDescription (std::ostream &outFile) const
virtual void	BuildPhysicsTable (const G4ParticleDefinition &)
virtual void	InitialiseModel ()
	G4HadronicInteraction (const G4HadronicInteraction &right)=delete
const G4HadronicInteraction &	operator= (const G4HadronicInteraction &right)=delete
G4bool	operator== (const G4HadronicInteraction &right) const =delete
G4bool	operator!= (const G4HadronicInteraction &right) const =delete

Additional Inherited Members
Protected Member Functions inherited from G4HadronicInteraction
void	SetModelName (const G4String &nam)
G4bool	IsBlocked () const
void	Block ()
Protected Attributes inherited from G4HadronicInteraction
G4HadFinalState	theParticleChange
G4int	verboseLevel
G4double	theMinEnergy
G4double	theMaxEnergy
G4bool	isBlocked

Detailed Description

Definition at line 49 of file G4LowEIonFragmentation.hh.

Constructor & Destructor Documentation

G4LowEIonFragmentation() [1/2]

G4LowEIonFragmentation::G4LowEIonFragmentation

(

G4ExcitationHandler *const

value = nullptr

)

Definition at line 49 of file G4LowEIonFragmentation.cc.

  : G4HadronicInteraction("LowEIonPreco")
{
  theHandler = value;
  theModel = new G4PreCompoundModel(theHandler);
  proton = G4Proton::Proton();
  secID = G4PhysicsModelCatalog::GetModelID("model_" + GetModelName());
}

~G4LowEIonFragmentation()

G4LowEIonFragmentation::~G4LowEIonFragmentation ( )

virtual

Definition at line 58 of file G4LowEIonFragmentation.cc.

{
  theResult.Clear();
}

G4LowEIonFragmentation() [2/2]

G4LowEIonFragmentation::G4LowEIonFragmentation ( const G4LowEIonFragmentation & )

delete

Member Function Documentation

ApplyYourself()

G4HadFinalState * G4LowEIonFragmentation::ApplyYourself ( const G4HadProjectile & thePrimary, G4Nucleus & theNucleus )

overridevirtual

Reimplemented from G4HadronicInteraction.

Definition at line 63 of file G4LowEIonFragmentation.cc.

{
  area = 0.0;
  // initialize the particle change
  theResult.Clear();
  theResult.SetStatusChange( stopAndKill );
  theResult.SetEnergyChange( 0.0 );
 
  // Get Target A, Z
  G4int aTargetA = theNucleus.GetA_asInt();
  G4int aTargetZ = theNucleus.GetZ_asInt();
 
  // Get Projectile A, Z
  G4int aProjectileA = thePrimary.GetDefinition()->GetBaryonNumber();
  G4int aProjectileZ = 
    G4lrint(thePrimary.GetDefinition()->GetPDGCharge()/eplus);
 
  // Get Maximum radius of both
  
  G4Fancy3DNucleus aPrim;
  aPrim.Init(aProjectileA, aProjectileZ);
  G4double projectileOuterRadius = aPrim.GetOuterRadius();
  
  G4Fancy3DNucleus aTarg;
  aTarg.Init(aTargetA, aTargetZ);
  G4double targetOuterRadius = aTarg.GetOuterRadius();
 
  // Get the Impact parameter
  G4int particlesFromProjectile = 0;
  G4int chargedFromProjectile = 0;
  G4double impactParameter = 0;
  G4double x,y;
  G4Nucleon * pNucleon;
  // need at lease one particle from the projectile model beyond the 
  // projectileHorizon.
 
  // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
  while(0==particlesFromProjectile)
  {
    do
    {
      x = 2*G4UniformRand() - 1;
      y = 2*G4UniformRand() - 1;
    }
    // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
    while(x*x + y*y > 1);
    impactParameter = std::sqrt(x*x+y*y)*
      (targetOuterRadius+projectileOuterRadius);
    ++totalTries;
    area = pi*(targetOuterRadius+projectileOuterRadius)*
              (targetOuterRadius+projectileOuterRadius);
    G4double projectileHorizon = impactParameter-targetOuterRadius; 
    
    // Empirical boundary transparency.
    G4double empirical = G4UniformRand();
    if(projectileHorizon > empirical*projectileOuterRadius) { continue; }
    
    // Calculate the number of nucleons involved in collision
    // From projectile
    aPrim.StartLoop();
 
    // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
    while((pNucleon = aPrim.GetNextNucleon()))
    {
      if(pNucleon->GetPosition().y()>projectileHorizon)
      {
        // We have one
        ++particlesFromProjectile;
        if(pNucleon->GetParticleType() == proton) 
        {
          ++chargedFromProjectile;
        } 
      }
    }
  }
  ++hits;
 
  // From target:
  G4double targetHorizon = impactParameter-projectileOuterRadius;
  G4int chargedFromTarget = 0;
  G4int particlesFromTarget = 0;
  aTarg.StartLoop();  
  // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
  while((pNucleon = aTarg.GetNextNucleon()))
  {
    if(pNucleon->GetPosition().y()>targetHorizon)
    {
      // We have one
      ++particlesFromTarget;
      if(pNucleon->GetParticleType() == proton) 
      {
        ++chargedFromTarget;
      }
    }
  }
  
  // Energy sharing between projectile and target. 
  // Note that this is a quite simplistic kinetically.
  G4ThreeVector momentum = thePrimary.Get4Momentum().vect();
  G4double w = (G4double)particlesFromProjectile/(G4double)aProjectileA;
  
  G4double projTotEnergy = thePrimary.GetTotalEnergy();  
  G4double targetMass = G4NucleiProperties::GetNuclearMass(aTargetA, aTargetZ);
  G4LorentzVector fragment4Momentum(momentum*w, projTotEnergy*w + targetMass);
 
  // take the nucleons and fill the Fragments
  G4Fragment anInitialState(aTargetA+particlesFromProjectile,
                aTargetZ+chargedFromProjectile,
                fragment4Momentum);
  // M.A. Cortes fix
  anInitialState.SetNumberOfExcitedParticle(particlesFromProjectile
                        + particlesFromTarget,
                        chargedFromProjectile 
                        + chargedFromTarget);
  anInitialState.SetNumberOfHoles(particlesFromProjectile+particlesFromTarget,
                  chargedFromProjectile + chargedFromTarget);
  G4double time = thePrimary.GetGlobalTime();
  anInitialState.SetCreationTime(time);
  anInitialState.SetCreatorModelID(secID);
 
  // Fragment the Fragment using Pre-compound
  G4ReactionProductVector* thePreCompoundResult = 
    theModel->DeExcite(anInitialState);
  
  // De-excite the projectile using ExcitationHandler
  G4ReactionProductVector * theExcitationResult = nullptr;
  if(particlesFromProjectile < aProjectileA)
  {
    G4LorentzVector residual4Momentum(momentum*(1.0-w), projTotEnergy*(1.0-w));
 
    G4Fragment initialState2(aProjectileA-particlesFromProjectile,
                 aProjectileZ-chargedFromProjectile,
                 residual4Momentum );
 
    // half of particles are excited (?!)
    G4int pinit = (aProjectileA-particlesFromProjectile)/2;
    G4int cinit = (aProjectileZ-chargedFromProjectile)/2;
 
    initialState2.SetNumberOfExcitedParticle(pinit,cinit);
    initialState2.SetNumberOfHoles(pinit,cinit);
    initialState2.SetCreationTime(time);
    initialState2.SetCreatorModelID(secID);
 
    theExcitationResult = theHandler->BreakItUp(initialState2);
  }
 
  // Fill the particle change and clear intermediate vectors
  std::size_t nexc = (nullptr != theExcitationResult) ? 
    theExcitationResult->size() : 0;
  std::size_t npre = (nullptr != thePreCompoundResult) ?
    thePreCompoundResult->size() : 0;
  
  for(std::size_t k=0; k<nexc; ++k) {
    G4ReactionProduct* p = (*theExcitationResult)[k];
    G4HadSecondary secondary(new G4DynamicParticle(p->GetDefinition(), p->GetMomentum()));
    secondary.SetTime(p->GetTOF());
    secondary.SetCreatorModelID(secID);
    theResult.AddSecondary(secondary);
    delete p;
  }
  for(std::size_t k=0; k<npre; ++k) {
    G4ReactionProduct* p = (*thePreCompoundResult)[k];
    G4HadSecondary secondary(new G4DynamicParticle(p->GetDefinition(), p->GetMomentum()));
    secondary.SetTime(p->GetTOF());
    secondary.SetCreatorModelID(secID);
    theResult.AddSecondary(secondary);
    delete p;
  }
  
  delete thePreCompoundResult;
  delete theExcitationResult;
 
  // return the particle change
  return &theResult;
}

GetCrossSection()

G4double G4LowEIonFragmentation::GetCrossSection ( )

inline

Definition at line 60 of file G4LowEIonFragmentation.hh.

  {
    //    G4cout << "area/millibarn = "<<area/millibarn<<G4endl;
    //    G4cout << "hits = "<<hits<<G4endl;
    //    G4cout << "totalTries = "<<totalTries<<G4endl;
    return area*hits/(static_cast<G4double>(totalTries)*CLHEP::millibarn);
  }

operator!=()

G4bool G4LowEIonFragmentation::operator!= ( const G4LowEIonFragmentation & right ) const

delete

operator=()

const G4LowEIonFragmentation & G4LowEIonFragmentation::operator= ( const G4LowEIonFragmentation & right )

delete

operator==()

G4bool G4LowEIonFragmentation::operator== ( const G4LowEIonFragmentation & right ) const

delete

---

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

• G4LowEIonFragmentation.hh
• G4LowEIonFragmentation.cc