G4PiMinusStopAbsorption Class Reference

#include <G4PiMinusStopAbsorption.hh>

Public Member Functions
	G4PiMinusStopAbsorption (G4PiMinusStopMaterial *materialAlgo, const G4double Z, const G4double A)
	~G4PiMinusStopAbsorption ()
G4DynamicParticleVector *	DoAbsorption ()
G4double	Energy ()
G4ThreeVector	RecoilMomentum ()
G4int	NProtons ()
G4int	NNeutrons ()
void	SetVerboseLevel (G4int level)

Detailed Description

Definition at line 46 of file G4PiMinusStopAbsorption.hh.

Constructor & Destructor Documentation

G4PiMinusStopAbsorption()

G4PiMinusStopAbsorption::G4PiMinusStopAbsorption	(	G4PiMinusStopMaterial *	materialAlgo,
		const G4double	Z,
		const G4double	A
	)

Definition at line 53 of file G4PiMinusStopAbsorption.cc.

{
  G4HadronicDeprecate("G4PiMinusStopAbsorption");
  _materialAlgo = materialAlgo;
  _nucleusZ = Z;
  _nucleusA = A;
  _level = 0;
  _absorptionProducts = new G4DynamicParticleVector();
}

~G4PiMinusStopAbsorption()

G4PiMinusStopAbsorption::~G4PiMinusStopAbsorption

(

)

Definition at line 67 of file G4PiMinusStopAbsorption.cc.

{
  // Memory management of materialAlgo needs better thought (MGP)
  delete _materialAlgo;
  // Who owns it? Memory management is not clear... (MGP)
  //  _absorptionProducts->clearAndDestroy();
  delete _absorptionProducts;
}

Member Function Documentation

DoAbsorption()

G4DynamicParticleVector * G4PiMinusStopAbsorption::DoAbsorption

(

)

Definition at line 76 of file G4PiMinusStopAbsorption.cc.

{
  std::vector<G4ParticleDefinition*>* defNucleons = _materialAlgo->DefinitionVector();
 
  G4double newA = _nucleusA;
  G4double newZ = _nucleusZ;
 
  if (defNucleons != 0)
    {
      for (unsigned int i=0; i<defNucleons->size(); i++)
    {
      if ( (*defNucleons)[i] == G4Proton::Proton())
        {
          newA = newA - 1;
          newZ = newZ - 1;
        }
      if ((*defNucleons)[i] == G4Neutron::Neutron()) 
        { newA = newA - 1; }
    }
    }
 
  G4double binding = G4NucleiProperties::GetBindingEnergy(static_cast<G4int>(_nucleusA) ,static_cast<G4int>(_nucleusZ)) / _nucleusA;
  G4double mass = G4NucleiProperties::GetNuclearMass(static_cast<G4int>(newA),static_cast<G4int>(newZ));
 
 
  std::vector<G4LorentzVector*>* p4Nucleons = _materialAlgo->P4Vector(binding,mass);
 
  if (defNucleons != 0 && p4Nucleons != 0)
    {
      unsigned int nNucleons = p4Nucleons->size();
      
      G4double seen = _materialAlgo->FinalNucleons() / 2.;
      G4int maxN = nNucleons;
      if (defNucleons->size() < nNucleons) { maxN = defNucleons->size(); }
 
      for (G4int i=0; i<maxN; i++)
    {
      G4DynamicParticle* product;
      if ((*defNucleons)[i] == G4Proton::Proton()) 
        { product = new G4DynamicParticle(G4Proton::Proton(),*((*p4Nucleons)[i])); }
      else
        { product = new G4DynamicParticle(G4Neutron::Neutron(),*((*p4Nucleons)[i])); }
      G4double ranflat = G4UniformRand();
      
      if (ranflat < seen) 
        { _absorptionProducts->push_back(product); }
          else
        { delete product; }
    }
    }
 
  return _absorptionProducts;
 
}

Referenced by G4PiMinusAbsorptionAtRest::AtRestDoIt().

Energy()

G4double G4PiMinusStopAbsorption::Energy

(

)

Definition at line 169 of file G4PiMinusStopAbsorption.cc.

{
  G4double energy = 0.;
  G4double productEnergy = 0.;
  G4ThreeVector pProducts(0.,0.,0.);
  G4int nN = 0;
  G4int nP = 0;
 
 
  G4int nAbsorptionProducts = _absorptionProducts->size();
 
  for (int i = 0; i<nAbsorptionProducts; i++)
    {
      productEnergy += (*_absorptionProducts)[i]->GetKineticEnergy();
      pProducts = pProducts + (*_absorptionProducts)[i]->GetMomentum();
      if ((*_absorptionProducts)[i]->GetDefinition() == G4Neutron::Neutron()) nN++;
      if ((*_absorptionProducts)[i]->GetDefinition() == G4Proton::Proton()) nP++;
    }
 
  G4double productBinding = (G4NucleiProperties::GetBindingEnergy(static_cast<G4int>(_nucleusA),static_cast<G4int>(_nucleusZ)) / _nucleusA) * nAbsorptionProducts;
  G4double mass = G4NucleiProperties::GetNuclearMass(_nucleusA - (nP + nN),_nucleusZ - nP);
  G4double pNucleus = pProducts.mag();
  G4double eNucleus = std::sqrt(pNucleus*pNucleus + mass*mass);
  G4double tNucleus = eNucleus - mass;
  G4double temp = 
    G4NucleiProperties::GetBindingEnergy(static_cast<G4int>(_nucleusA - (nP + nN)),static_cast<G4int>(_nucleusZ - nP)) - 
    G4NucleiProperties::GetBindingEnergy(static_cast<G4int>(_nucleusA),static_cast<G4int>(_nucleusZ));
  energy = productEnergy + productBinding + tNucleus;
  
  if (_level > 0)
    {
      std::cout << "E products " <<  productEnergy  
       << " Binding " << productBinding << " " << temp << " "
       << " Tnucleus " << tNucleus 
       << " energy = " << energy << G4endl;
    }
 
  return energy;
}

Referenced by G4PiMinusAbsorptionAtRest::AtRestDoIt().

NNeutrons()

G4int G4PiMinusStopAbsorption::NNeutrons

(

)

Definition at line 156 of file G4PiMinusStopAbsorption.cc.

{
  G4int n = 0;
  G4int entries = _absorptionProducts->size();
  for (int i = 0; i<entries; i++)
    {
      if ((*_absorptionProducts)[i]->GetDefinition() == G4Neutron::Neutron())
    { n = n + 1; }
    }
  return n;
}

Referenced by G4PiMinusAbsorptionAtRest::AtRestDoIt().

NProtons()

G4int G4PiMinusStopAbsorption::NProtons

(

)

Definition at line 143 of file G4PiMinusStopAbsorption.cc.

{
  G4int n = 0;
  G4int entries = _absorptionProducts->size();
  for (int i = 0; i<entries; i++)
    {
      if ((*_absorptionProducts)[i]->GetDefinition() == G4Proton::Proton())
    { n = n + 1; }
    }
  return n;
}

Referenced by G4PiMinusAbsorptionAtRest::AtRestDoIt().

RecoilMomentum()

G4ThreeVector G4PiMinusStopAbsorption::RecoilMomentum

(

)

Definition at line 131 of file G4PiMinusStopAbsorption.cc.

{
  G4ThreeVector pProducts(0.,0.,0.);
  
  for (unsigned int i = 0; i< _absorptionProducts->size(); i++)
    {
      pProducts = pProducts + (*_absorptionProducts)[i]->GetMomentum();
    }
  return pProducts;
}

Referenced by G4PiMinusAbsorptionAtRest::AtRestDoIt().

SetVerboseLevel()

void G4PiMinusStopAbsorption::SetVerboseLevel

(

G4int

level

)

Definition at line 209 of file G4PiMinusStopAbsorption.cc.

{
  _level = level;
  return;
}

Referenced by G4PiMinusAbsorptionAtRest::AtRestDoIt().

---

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

• G4PiMinusStopAbsorption.hh
• G4PiMinusStopAbsorption.cc