G4LowEIonFragmentation Class Reference

#include <G4LowEIonFragmentation.hh>

Inheritance diagram for G4LowEIonFragmentation:

Public Member Functions
	G4LowEIonFragmentation (G4ExcitationHandler *const value)
	G4LowEIonFragmentation ()
virtual	~G4LowEIonFragmentation ()
virtual G4HadFinalState *	ApplyYourself (const G4HadProjectile &thePrimary, G4Nucleus &theNucleus)
Static Public Member Functions
static G4double	GetCrossSection ()

---

Detailed Description

Definition at line 50 of file G4LowEIonFragmentation.hh.

---

Constructor & Destructor Documentation

G4LowEIonFragmentation::G4LowEIonFragmentation

(

G4ExcitationHandler *const

value

)

Definition at line 53 of file G4LowEIonFragmentation.cc.

References G4Proton::Proton().

{
  theHandler = value;
  theModel = new G4PreCompoundModel(theHandler);
  proton = G4Proton::Proton();
}

G4LowEIonFragmentation::G4LowEIonFragmentation

(

)

Definition at line 60 of file G4LowEIonFragmentation.cc.

References G4Proton::Proton().

{
  theHandler = new G4ExcitationHandler;
  theModel = new G4PreCompoundModel(theHandler);
  proton = G4Proton::Proton();
}

G4LowEIonFragmentation::~G4LowEIonFragmentation

(

)

[virtual]

Definition at line 67 of file G4LowEIonFragmentation.cc.

{
  delete theModel;
}

---

Member Function Documentation

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

Implements G4HadronicInteraction.

Definition at line 73 of file G4LowEIonFragmentation.cc.

References G4HadFinalState::AddSecondary(), G4ExcitationHandler::BreakItUp(), G4HadFinalState::Clear(), G4PreCompoundModel::DeExcite(), G4lrint(), G4UniformRand, G4HadProjectile::Get4Momentum(), G4Nucleus::GetA_asInt(), G4ParticleDefinition::GetBaryonNumber(), G4HadProjectile::GetDefinition(), G4HadProjectile::GetGlobalTime(), G4Fancy3DNucleus::GetNextNucleon(), G4NucleiProperties::GetNuclearMass(), G4Fancy3DNucleus::GetOuterRadius(), G4Nucleon::GetParticleType(), G4ParticleDefinition::GetPDGCharge(), G4Nucleon::GetPosition(), G4HadProjectile::GetTotalEnergy(), G4Nucleus::GetZ_asInt(), G4Fancy3DNucleus::Init(), G4INCL::Math::pi, G4Fragment::SetCreationTime(), G4HadFinalState::SetEnergyChange(), G4Fragment::SetNumberOfExcitedParticle(), G4Fragment::SetNumberOfHoles(), G4HadFinalState::SetStatusChange(), G4Fancy3DNucleus::StartLoop(), and stopAndKill.

{
  area = 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.
  while(0==particlesFromProjectile)
  {
    do
    {
      x = 2*G4UniformRand() - 1;
      y = 2*G4UniformRand() - 1;
    }
    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();
    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();  
  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;
  
}

static G4double G4LowEIonFragmentation::GetCrossSection

(

)

[inline, static]

Definition at line 63 of file G4LowEIonFragmentation.hh.

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

---

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

• G4LowEIonFragmentation.hh
• G4LowEIonFragmentation.cc

---