G4QMDNucleus Class Reference

#include <G4QMDNucleus.hh>

Inheritance diagram for G4QMDNucleus:

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

Additional Inherited Members
Protected Attributes inherited from G4QMDSystem
std::vector< G4QMDParticipant * >	participants

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;
}

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];
   }
 
   if (tmass != 0.)
   {
     rcm0 = rcm0/tmass;
   }
   else
   {
     rcm0 = G4ThreeVector(0.);
   }
 
   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; 
 
}

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

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;
}

Referenced by CalEnergyAndAngularMomentumInCM().

GetAngularMomentum()

G4int G4QMDNucleus::GetAngularMomentum ( )

inline

Definition at line 65 of file G4QMDNucleus.hh.

{ return 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;
}

Referenced by CalEnergyAndAngularMomentumInCM(), and GetNuclearMass().

GetExcitationEnergy()

G4double G4QMDNucleus::GetExcitationEnergy ( )

inline

Definition at line 63 of file G4QMDNucleus.hh.

{ return 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;
}

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

GetNuclearMass()

G4double G4QMDNucleus::GetNuclearMass

(

)

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;
}

SetTotalPotential()

void G4QMDNucleus::SetTotalPotential ( G4double x )

inline

Definition at line 62 of file G4QMDNucleus.hh.

{ potentialEnergy = x; };

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

---

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

• G4QMDNucleus.hh
• G4QMDNucleus.cc