G4QMDNucleus Class Reference

#include <G4QMDNucleus.hh>

Inheritance diagram for G4QMDNucleus:

Public Member Functions
void	CalEnergyAndAngularMomentumInCM ()
void	Clear ()
void	DeleteParticipant (G4int i)
G4QMDParticipant *	EraseParticipant (G4int i)
	G4QMDNucleus ()
G4LorentzVector	Get4Momentum ()
G4int	GetAngularMomentum ()
G4int	GetAtomicNumber ()
G4double	GetExcitationEnergy ()
G4int	GetMassNumber ()
G4int	GetNOCollision ()
G4double	GetNuclearMass ()
G4QMDParticipant *	GetParticipant (G4int i)
G4int	GetTotalNumberOfParticipant ()
void	IncrementCollisionCounter ()
void	InsertParticipant (G4QMDParticipant *particle, G4int j)
void	SetParticipant (G4QMDParticipant *particle)
void	SetSystem (G4QMDSystem *, G4ThreeVector, G4ThreeVector)
void	SetTotalPotential (G4double x)
void	ShowParticipants ()
void	SubtractSystem (G4QMDSystem *)

Protected Attributes
std::vector< G4QMDParticipant * >	participants

Private Attributes
std::vector< G4double >	es
G4double	excitationEnergy
G4double	hbc
G4int	jj
G4int	numberOfCollision
std::vector< G4ThreeVector >	pcm
G4double	potentialEnergy
std::vector< G4ThreeVector >	rcm

Detailed Description

Definition at line 43 of file G4QMDNucleus.hh.

Constructor & Destructor Documentation

G4QMDNucleus()

G4QMDNucleus::G4QMDNucleus

(

)

Definition at line 38 of file G4QMDNucleus.cc.

{
   G4QMDParameters* parameters = G4QMDParameters::GetInstance();
   hbc = parameters->Get_hbc();
   
   jj = 0; // will be calcualted in CalEnergyAndAngularMomentumInCM;
   potentialEnergy = 0.0; // will be set through set method 
   excitationEnergy = 0.0;
}

References excitationEnergy, G4QMDParameters::Get_hbc(), G4QMDParameters::GetInstance(), hbc, jj, and potentialEnergy.

Member Function Documentation

CalEnergyAndAngularMomentumInCM()

void G4QMDNucleus::CalEnergyAndAngularMomentumInCM

(

)

Definition at line 138 of file G4QMDNucleus.cc.

{
 
   //G4cout << "CalEnergyAndAngularMomentumInCM " << this->GetAtomicNumber() << " " << GetMassNumber() << G4endl;
 
   G4double gamma = Get4Momentum().gamma();
   G4ThreeVector beta = Get4Momentum().v()/ Get4Momentum().e();
 
   G4ThreeVector pcm0( 0.0 ) ;
 
   G4int n = GetTotalNumberOfParticipant();
   pcm.resize( n );
 
   for ( G4int i= 0; i < n ; i++ ) 
   {
      G4ThreeVector p_i = GetParticipant( i )->GetMomentum();
 
      G4double trans = gamma / ( gamma + 1.0 ) * p_i * beta; 
      pcm[i] = p_i - trans*beta;
 
      pcm0 += pcm[i];
   }
 
   pcm0 = pcm0 / double ( n );
 
   //G4cout << "pcm0 " << pcm0 << G4endl;
 
   for ( G4int i= 0; i < n ; i++ ) 
   {
      pcm[i] += -pcm0;
      //G4cout << "pcm " << i << " " << pcm[i] << G4endl;
   }
 
 
   G4double tmass = 0;
   G4ThreeVector rcm0( 0.0 ) ;
   rcm.resize( n );
   es.resize( n );
 
   for ( G4int i= 0; i < n ; i++ ) 
   {
      G4ThreeVector ri = GetParticipant( i )->GetPosition();
      G4double trans = gamma / ( gamma + 1.0 ) * ri * beta; 
 
      es[i] = std::sqrt ( G4Pow::GetInstance()->powN ( GetParticipant( i )->GetMass() , 2 ) + pcm[i]*pcm[i] );
 
      rcm[i] = ri + trans*beta;
 
      rcm0 += rcm[i]*es[i];
 
      tmass += es[i];
   }
 
   rcm0 = rcm0/tmass;
 
   for ( G4int i= 0; i < n ; i++ ) 
   {
      rcm[i] += -rcm0;
      //G4cout << "rcm " << i << " " << rcm[i] << G4endl;
   }
 
// Angular momentum
 
   G4ThreeVector rl ( 0.0 ); 
   for ( G4int i= 0; i < n ; i++ ) 
   {
      rl += rcm[i].cross ( pcm[i] );
   }
 
// DHW: move hbc outside of sqrt to get correct units
//  jj = int ( std::sqrt ( rl*rl / hbc ) + 0.5 );
 
   jj = int (std::sqrt(rl*rl)/hbc + 0.5);
 
// kinetic energy per nucleon in CM
 
   G4double totalMass = 0.0;
   for ( G4int i= 0; i < n ; i++ ) 
   {
      // following two lines are equivalent
      //totalMass += GetParticipant( i )->GetDefinition()->GetPDGMass()/GeV;
      totalMass += GetParticipant( i )->GetMass();
   }
 
   //G4double kineticEnergyPerNucleon = ( std::accumulate ( es.begin() , es.end() , 0.0 ) - totalMass )/n;
 
// Total (not per nucleion ) Binding Energy 
   G4double bindingEnergy =  ( std::accumulate ( es.begin() , es.end() , 0.0 ) -totalMass ) + potentialEnergy;
 
   //G4cout << "KineticEnergyPerNucleon in GeV " << kineticEnergyPerNucleon << G4endl;
   //G4cout << "KineticEnergySum in GeV " << std::accumulate ( es.begin() , es.end() , 0.0 ) - totalMass << G4endl;
   //G4cout << "PotentialEnergy in GeV " << potentialEnergy << G4endl;
   //G4cout << "BindingEnergy in GeV " << bindingEnergy << G4endl;
   //G4cout << "G4BindingEnergy in GeV " << G4NucleiProperties::GetBindingEnergy( GetAtomicNumber() , GetMassNumber() )/GeV << G4endl;
 
   excitationEnergy = bindingEnergy + G4NucleiProperties::GetBindingEnergy( GetMassNumber() , GetAtomicNumber() )/GeV;
   //G4cout << "excitationEnergy in GeV " << excitationEnergy << G4endl;
   if ( excitationEnergy < 0 ) excitationEnergy = 0.0; 
 
}

References anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, G4InuclSpecialFunctions::bindingEnergy(), CLHEP::HepLorentzVector::e(), es, excitationEnergy, CLHEP::HepLorentzVector::gamma(), Get4Momentum(), GetAtomicNumber(), G4NucleiProperties::GetBindingEnergy(), G4Pow::GetInstance(), G4QMDParticipant::GetMass(), GetMassNumber(), G4QMDParticipant::GetMomentum(), G4QMDSystem::GetParticipant(), G4QMDParticipant::GetPosition(), G4QMDSystem::GetTotalNumberOfParticipant(), GeV, hbc, jj, CLHEP::detail::n, pcm, potentialEnergy, rcm, and CLHEP::HepLorentzVector::v().

Referenced by G4QMDReaction::ApplyYourself(), and G4QMDMeanField::SetNucleus().

Clear()

void G4QMDSystem::Clear ( )

inherited

Definition at line 68 of file G4QMDSystem.cc.

{
   for ( G4int i = 0 ; i < this->GetTotalNumberOfParticipant() ; i++ )
   {
      delete participants[i];
   }
   participants.clear();
}

References G4QMDSystem::GetTotalNumberOfParticipant(), and G4QMDSystem::participants.

Referenced by G4QMDReaction::ApplyYourself(), and G4QMDSystem::~G4QMDSystem().

DeleteParticipant()

void G4QMDSystem::DeleteParticipant ( G4int  i )

inlineinherited

Definition at line 57 of file G4QMDSystem.hh.

{ delete participants[ i ] ; participants.erase( std::find ( participants.begin() , participants.end() , participants[ i ] ) ); };

References G4QMDSystem::participants.

Referenced by G4QMDCollision::CalKinematicsOfBinaryCollisions().

EraseParticipant()

G4QMDParticipant * G4QMDSystem::EraseParticipant ( G4int  i )

inlineinherited

Definition at line 56 of file G4QMDSystem.hh.

{ G4QMDParticipant* particle =  participants[ i ]; participants.erase( std::find ( participants.begin() , participants.end() , participants[ i ] ) ) ; return particle; };

References G4QMDSystem::participants.

Referenced by G4QMDCollision::CalFinalStateOfTheBinaryCollision().

Get4Momentum()

G4LorentzVector G4QMDNucleus::Get4Momentum

(

)

Definition at line 56 of file G4QMDNucleus.cc.

{
   G4LorentzVector p( 0 );
   std::vector< G4QMDParticipant* >::iterator it; 
   for ( it = participants.begin() ; it != participants.end() ; it++ ) 
      p += (*it)->Get4Momentum();   
 
   return p;
}

References G4QMDSystem::participants.

Referenced by CalEnergyAndAngularMomentumInCM().

GetAngularMomentum()

G4int G4QMDNucleus::GetAngularMomentum ( )

inline

Definition at line 65 of file G4QMDNucleus.hh.

{ return jj; };

References jj.

GetAtomicNumber()

G4int G4QMDNucleus::GetAtomicNumber

(

)

Definition at line 89 of file G4QMDNucleus.cc.

{
   G4int Z = 0; 
   std::vector< G4QMDParticipant* >::iterator it; 
   for ( it = participants.begin() ; it != participants.end() ; it++ ) 
   {
      if ( (*it)->GetDefinition() == G4Proton::Proton() ) 
         Z++; 
   }
   return Z;
}

References G4QMDSystem::participants, G4Proton::Proton(), and Z.

Referenced by CalEnergyAndAngularMomentumInCM(), GetNuclearMass(), and G4QMDGroundStateNucleus::packNucleons().

GetExcitationEnergy()

G4double G4QMDNucleus::GetExcitationEnergy ( void  )

inline

Definition at line 63 of file G4QMDNucleus.hh.

{ return excitationEnergy; };

References excitationEnergy.

GetMassNumber()

G4int G4QMDNucleus::GetMassNumber

(

)

Definition at line 68 of file G4QMDNucleus.cc.

{
 
   G4int A = 0; 
   std::vector< G4QMDParticipant* >::iterator it; 
   for ( it = participants.begin() ; it != participants.end() ; it++ ) 
   {
      if ( (*it)->GetDefinition() == G4Proton::Proton() 
        || (*it)->GetDefinition() == G4Neutron::Neutron() ) 
         A++; 
   }
 
   if ( A == 0 ) {
      throw G4HadronicException(__FILE__, __LINE__, "G4QMDNucleus has the mass number of 0!");
   }
 
   return A;
}

References A, G4Neutron::Neutron(), G4QMDSystem::participants, and G4Proton::Proton().

Referenced by CalEnergyAndAngularMomentumInCM(), G4QMDGroundStateNucleus::G4QMDGroundStateNucleus(), GetNuclearMass(), G4QMDGroundStateNucleus::killCMMotionAndAngularM(), G4QMDGroundStateNucleus::packNucleons(), and G4QMDGroundStateNucleus::samplingMomentum().

GetNOCollision()

G4int G4QMDSystem::GetNOCollision ( )

inlineinherited

Definition at line 65 of file G4QMDSystem.hh.

{ return numberOfCollision; };

References G4QMDSystem::numberOfCollision.

Referenced by G4QMDReaction::ApplyYourself().

GetNuclearMass()

G4double G4QMDNucleus::GetNuclearMass

(

void

)

Definition at line 103 of file G4QMDNucleus.cc.

{
 
   G4double mass = G4NucleiProperties::GetNuclearMass( GetMassNumber() , GetAtomicNumber() );
   
   if ( mass == 0.0 )
   {
 
      G4int Z = GetAtomicNumber();
      G4int A = GetMassNumber();
      G4int N = A - Z;
 
// Weizsacker-Bethe 
 
      G4double Av = 16*MeV; 
      G4double As = 17*MeV; 
      G4double Ac = 0.7*MeV; 
      G4double Asym = 23*MeV; 
 
      G4double BE = Av * A 
                  - As * G4Pow::GetInstance()->A23 ( G4double ( A ) ) 
                  - Ac * Z*Z/G4Pow::GetInstance()->A13 ( G4double ( A ) )
                  - Asym * ( N - Z )* ( N - Z ) / A; 
 
      mass = Z * G4Proton::Proton()->GetPDGMass() 
           + N * G4Neutron::Neutron()->GetPDGMass()
           - BE;
 
   }
 
   return mass;
}

References A, G4Pow::A13(), G4Pow::A23(), GetAtomicNumber(), G4Pow::GetInstance(), GetMassNumber(), G4NucleiProperties::GetNuclearMass(), G4ParticleDefinition::GetPDGMass(), MeV, G4Neutron::Neutron(), G4Proton::Proton(), and Z.

GetParticipant()

G4QMDParticipant * G4QMDSystem::GetParticipant ( G4int  i )

inlineinherited

Definition at line 62 of file G4QMDSystem.hh.

{ return participants[i]; };

References G4QMDSystem::participants.

Referenced by G4QMDReaction::ApplyYourself(), G4QMDMeanField::Cal2BodyQuantities(), CalEnergyAndAngularMomentumInCM(), G4QMDCollision::CalFinalStateOfTheBinaryCollision(), G4QMDCollision::CalFinalStateOfTheBinaryCollisionJQMD(), G4QMDMeanField::CalGraduate(), G4QMDCollision::CalKinematicsOfBinaryCollisions(), G4QMDMeanField::calPauliBlockingFactor(), G4QMDMeanField::DoClusterJudgment(), G4QMDMeanField::DoPropagation(), G4QMDMeanField::GetPotential(), G4QMDMeanField::GetTotalPotential(), G4QMDMeanField::IsPauliBlocked(), and G4QMDSystem::SubtractSystem().

GetTotalNumberOfParticipant()

G4int G4QMDSystem::GetTotalNumberOfParticipant ( )

inlineinherited

Definition at line 60 of file G4QMDSystem.hh.

{ return participants.size(); };

References G4QMDSystem::participants.

Referenced by G4QMDReaction::ApplyYourself(), G4QMDMeanField::Cal2BodyQuantities(), CalEnergyAndAngularMomentumInCM(), G4QMDMeanField::CalGraduate(), G4QMDCollision::CalKinematicsOfBinaryCollisions(), G4QMDMeanField::calPauliBlockingFactor(), G4QMDSystem::Clear(), G4QMDMeanField::DoClusterJudgment(), G4QMDMeanField::DoPropagation(), G4QMDMeanField::GetPotential(), G4QMDMeanField::GetTotalPotential(), G4QMDMeanField::SetSystem(), and G4QMDSystem::SubtractSystem().

IncrementCollisionCounter()

void G4QMDSystem::IncrementCollisionCounter ( )

inlineinherited

Definition at line 64 of file G4QMDSystem.hh.

{ numberOfCollision++; };

References G4QMDSystem::numberOfCollision.

Referenced by G4QMDCollision::CalKinematicsOfBinaryCollisions().

InsertParticipant()

void G4QMDSystem::InsertParticipant ( G4QMDParticipant *  particle, G4int  j  )

inherited

Definition at line 110 of file G4QMDSystem.cc.

{
 
   if ( (size_t) n > participants.size()+1 )
      G4cout << "G4QMDSystem::InsertParticipant size error" << G4endl;
 
   std::vector< G4QMDParticipant* >::iterator it; 
   it = participants.begin();
 
   for ( G4int i = 0; i < n ; i++ )
      it++;
 
   participants.insert( it, particle );
}

References G4cout, G4endl, CLHEP::detail::n, and G4QMDSystem::participants.

Referenced by G4QMDCollision::CalFinalStateOfTheBinaryCollision().

SetParticipant()

void G4QMDSystem::SetParticipant ( G4QMDParticipant *  particle )

inlineinherited

Definition at line 51 of file G4QMDSystem.hh.

{ participants.push_back ( particle ); };

References G4QMDSystem::participants.

Referenced by G4QMDReaction::ApplyYourself(), G4QMDCollision::CalKinematicsOfBinaryCollisions(), G4QMDMeanField::DoClusterJudgment(), G4QMDGroundStateNucleus::G4QMDGroundStateNucleus(), and G4QMDSystem::SetSystem().

SetSystem()

void G4QMDSystem::SetSystem ( G4QMDSystem *  nucleus, G4ThreeVector  dp, G4ThreeVector  dr  )

inherited

Definition at line 46 of file G4QMDSystem.cc.

{
   std::vector< G4QMDParticipant* >::iterator it; 
   for ( it = nucleus->participants.begin() ; it != nucleus->participants.end() ; it++ ) 
   {
      G4ThreeVector r = (*it)->GetPosition() + dr;
      (*it)->SetPosition ( r );
      G4ThreeVector p = (*it)->GetMomentum() + dp;
      (*it)->SetMomentum ( p );
      this->SetParticipant( *it );
   }
}

References G4QMDSystem::participants, and G4QMDSystem::SetParticipant().

SetTotalPotential()

void G4QMDNucleus::SetTotalPotential ( G4double  x )

inline

Definition at line 62 of file G4QMDNucleus.hh.

{ potentialEnergy = x; };

References potentialEnergy.

Referenced by G4QMDReaction::ApplyYourself(), and G4QMDMeanField::SetNucleus().

ShowParticipants()

void G4QMDSystem::ShowParticipants ( )

inherited

Definition at line 79 of file G4QMDSystem.cc.

{
   //store orginal precision
   std::ios::fmtflags oldform = G4cout.flags();
 
   G4ThreeVector p_sum( 0.0 );
   std::vector< G4QMDParticipant* >::iterator it; 
   G4cout << "Momentum and Position of each participant " << G4endl; 
   G4int i = 0; 
   for ( it = participants.begin() ; it != participants.end() ; it++ ) 
   { 
      G4cout << i
             << " " 
             << (*it)->GetDefinition()->GetParticleName() 
             << " " 
             << std::setprecision( 8 )
             << (*it)->GetMomentum()
             << " " 
             << (*it)->GetPosition() 
             << G4endl;
      p_sum += (*it)->GetMomentum();
      i++;
   }
   G4cout << "Sum upped Momentum and its mag " << p_sum << " " << p_sum.mag() << G4endl;
 
   //restore orginal precision
   G4cout.flags( oldform );
}

References G4cout, G4endl, CLHEP::Hep3Vector::mag(), and G4QMDSystem::participants.

SubtractSystem()

void G4QMDSystem::SubtractSystem ( G4QMDSystem *  nucleus )

inherited

Definition at line 59 of file G4QMDSystem.cc.

{
   
   for ( G4int i = 0 ; i < nucleus->GetTotalNumberOfParticipant() ; i++ )
   {
      participants.erase ( std::find ( participants.begin() , participants.end() , nucleus->GetParticipant( i ) ) );
   }
}

References G4QMDSystem::GetParticipant(), G4QMDSystem::GetTotalNumberOfParticipant(), and G4QMDSystem::participants.

Referenced by G4QMDMeanField::DoClusterJudgment().

Field Documentation

es

std::vector< G4double > G4QMDNucleus::es

private

Definition at line 72 of file G4QMDNucleus.hh.

Referenced by CalEnergyAndAngularMomentumInCM().

excitationEnergy

G4double G4QMDNucleus::excitationEnergy

private

Definition at line 76 of file G4QMDNucleus.hh.

Referenced by CalEnergyAndAngularMomentumInCM(), G4QMDNucleus(), and GetExcitationEnergy().

hbc

G4double G4QMDNucleus::hbc

private

Definition at line 69 of file G4QMDNucleus.hh.

Referenced by CalEnergyAndAngularMomentumInCM(), and G4QMDNucleus().

jj

G4int G4QMDNucleus::jj

private

Definition at line 73 of file G4QMDNucleus.hh.

Referenced by CalEnergyAndAngularMomentumInCM(), G4QMDNucleus(), and GetAngularMomentum().

numberOfCollision

G4int G4QMDSystem::numberOfCollision

privateinherited

Definition at line 75 of file G4QMDSystem.hh.

Referenced by G4QMDSystem::G4QMDSystem(), G4QMDSystem::GetNOCollision(), and G4QMDSystem::IncrementCollisionCounter().

participants

std::vector< G4QMDParticipant* > G4QMDSystem::participants

protectedinherited

Definition at line 72 of file G4QMDSystem.hh.

Referenced by G4QMDSystem::Clear(), G4QMDSystem::DeleteParticipant(), G4QMDSystem::EraseParticipant(), G4QMDSystem::G4QMDSystem(), Get4Momentum(), GetAtomicNumber(), GetMassNumber(), G4QMDSystem::GetParticipant(), G4QMDSystem::GetTotalNumberOfParticipant(), G4QMDSystem::InsertParticipant(), G4QMDGroundStateNucleus::killCMMotionAndAngularM(), G4QMDGroundStateNucleus::packNucleons(), G4QMDGroundStateNucleus::samplingMomentum(), G4QMDGroundStateNucleus::samplingPosition(), G4QMDSystem::SetParticipant(), G4QMDSystem::SetSystem(), G4QMDSystem::ShowParticipants(), and G4QMDSystem::SubtractSystem().

pcm

std::vector< G4ThreeVector > G4QMDNucleus::pcm

private

Definition at line 71 of file G4QMDNucleus.hh.

Referenced by CalEnergyAndAngularMomentumInCM().

potentialEnergy

G4double G4QMDNucleus::potentialEnergy

private

Definition at line 75 of file G4QMDNucleus.hh.

Referenced by CalEnergyAndAngularMomentumInCM(), G4QMDNucleus(), and SetTotalPotential().

rcm

std::vector< G4ThreeVector > G4QMDNucleus::rcm

private

Definition at line 71 of file G4QMDNucleus.hh.

Referenced by CalEnergyAndAngularMomentumInCM().

---

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

• source/processes/hadronic/models/qmd/include/G4QMDNucleus.hh
• source/processes/hadronic/models/qmd/src/G4QMDNucleus.cc