G4LowEIonFragmentation Class Reference

#include <G4LowEIonFragmentation.hh>

Inheritance diagram for G4LowEIonFragmentation:

Collaboration diagram for G4LowEIonFragmentation:

Public Member Functions
	G4LowEIonFragmentation (G4ExcitationHandler *const value)
	G4LowEIonFragmentation ()
virtual	~G4LowEIonFragmentation ()
virtual G4HadFinalState *	ApplyYourself (const G4HadProjectile &thePrimary, G4Nucleus &theNucleus)
G4double	GetCrossSection ()
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 &, G4Nucleus &)
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)
const G4HadronicInteraction *	GetMyPointer () const
virtual G4int	GetVerboseLevel () const
virtual 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
G4bool	operator== (const G4HadronicInteraction &right) const
G4bool	operator!= (const G4HadronicInteraction &right) 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 &)

Private Member Functions
	G4LowEIonFragmentation (const G4LowEIonFragmentation &)
const G4LowEIonFragmentation &	operator= (const G4LowEIonFragmentation &right)
G4bool	operator== (const G4LowEIonFragmentation &right) const
G4bool	operator!= (const G4LowEIonFragmentation &right) const

Private Attributes
G4HadFinalState	theResult
const G4ParticleDefinition *	proton
G4PreCompoundModel *	theModel
G4ExcitationHandler *	theHandler
G4int	hits
G4int	totalTries
G4double	area

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 51 of file G4LowEIonFragmentation.hh.

Constructor & Destructor Documentation

G4LowEIonFragmentation() [1/3]

G4LowEIonFragmentation::G4LowEIonFragmentation

(

G4ExcitationHandler *const

value

)

Definition at line 49 of file G4LowEIonFragmentation.cc.

{
  theHandler = value;
  theModel = new G4PreCompoundModel(theHandler);
  proton = G4Proton::Proton();
  hits = 0;
  totalTries = 1;
  area = 0.0;
}

Here is the call graph for this function:

G4LowEIonFragmentation() [2/3]

G4LowEIonFragmentation::G4LowEIonFragmentation

(

)

Definition at line 59 of file G4LowEIonFragmentation.cc.

{
  theHandler = new G4ExcitationHandler;
  theModel = new G4PreCompoundModel(theHandler);
  proton = G4Proton::Proton();
  hits = 0;
  totalTries = 1;
  area = 0.0;
}

Here is the call graph for this function:

Here is the caller graph for this function:

~G4LowEIonFragmentation()

G4LowEIonFragmentation::~G4LowEIonFragmentation ( )

virtual

Definition at line 69 of file G4LowEIonFragmentation.cc.

{
  delete theModel;
}

Here is the call graph for this function:

G4LowEIonFragmentation() [3/3]

G4LowEIonFragmentation::G4LowEIonFragmentation ( const G4LowEIonFragmentation &  )

private

Member Function Documentation

ApplyYourself()

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

virtual

Implements G4HadronicInteraction.

Definition at line 75 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.SetNumberOfParticles(particlesFromProjectile);
  anInitialState.SetNumberOfExcitedParticle(particlesFromProjectile
                        + particlesFromTarget,
                        chargedFromProjectile 
                        + chargedFromTarget);
  anInitialState.SetNumberOfHoles(particlesFromProjectile+particlesFromTarget,
                  chargedFromProjectile + chargedFromTarget);
  G4double time = thePrimary.GetGlobalTime();
  anInitialState.SetCreationTime(time);

  // Fragment the Fragment using Pre-compound
  G4ReactionProductVector* thePreCompoundResult = 
    theModel->DeExcite(anInitialState);
  
  // De-excite the projectile using ExcitationHandler
  G4ReactionProductVector * theExcitationResult = 0; 
  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);

    theExcitationResult = theHandler->BreakItUp(initialState2);
  }

  // Fill the particle change and clear intermediate vectors
  G4int nexc = 0;
  G4int npre = 0;
  if(theExcitationResult)  { nexc = theExcitationResult->size(); }
  if(thePreCompoundResult) { npre = thePreCompoundResult->size();}
  
  if(nexc > 0) {
    for(G4int k=0; k<nexc; ++k) {
      G4ReactionProduct* p = (*theExcitationResult)[k];
      theResult.AddSecondary(new G4DynamicParticle(p->GetDefinition(),
                           p->GetMomentum()));
      delete p;
    }
  }
  
  if(npre > 0) {
    for(G4int k=0; k<npre; ++k) {
      G4ReactionProduct* p = (*thePreCompoundResult)[k];
      theResult.AddSecondary(new G4DynamicParticle(p->GetDefinition(),
                           p->GetMomentum()));
      delete p;
    }
  }
  
  delete thePreCompoundResult;
  delete theExcitationResult;

  // return the particle change
  return &theResult;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetCrossSection()

G4double G4LowEIonFragmentation::GetCrossSection ( )

inline

Definition at line 64 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);
  }

Here is the call graph for this function:

operator!=()

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

private

Here is the caller graph for this function:

operator=()

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

private

Here is the caller graph for this function:

operator==()

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

private

Here is the caller graph for this function:

Member Data Documentation

area

G4double G4LowEIonFragmentation::area

private

Definition at line 91 of file G4LowEIonFragmentation.hh.

hits

G4int G4LowEIonFragmentation::hits

private

Definition at line 89 of file G4LowEIonFragmentation.hh.

proton

const G4ParticleDefinition* G4LowEIonFragmentation::proton

private

Definition at line 82 of file G4LowEIonFragmentation.hh.

theHandler

G4ExcitationHandler* G4LowEIonFragmentation::theHandler

private

Definition at line 85 of file G4LowEIonFragmentation.hh.

theModel

G4PreCompoundModel* G4LowEIonFragmentation::theModel

private

Definition at line 84 of file G4LowEIonFragmentation.hh.

theResult

G4HadFinalState G4LowEIonFragmentation::theResult

private

Definition at line 81 of file G4LowEIonFragmentation.hh.

totalTries

G4int G4LowEIonFragmentation::totalTries

private

Definition at line 90 of file G4LowEIonFragmentation.hh.

---

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

• G4LowEIonFragmentation.hh
• G4LowEIonFragmentation.cc