G4LightIonQMDNucleus Class Reference

#include <G4LightIonQMDNucleus.hh>

Inheritance diagram for G4LightIonQMDNucleus:

Public Member Functions
	G4LightIonQMDNucleus ()
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 44 of file G4LightIonQMDNucleus.hh.

Constructor & Destructor Documentation

G4LightIonQMDNucleus()

G4LightIonQMDNucleus::G4LightIonQMDNucleus

(

)

Definition at line 46 of file G4LightIonQMDNucleus.cc.

{
   G4LightIonQMDParameters* parameters = G4LightIonQMDParameters::GetInstance();
   hbc = parameters->Get_hbc();
   
   jj = 0; // will be calcualted in CalEnergyAndAngularMomentumInCM;
   potentialEnergy = 0.0; // will be set through set method 
   excitationEnergy = 0.0;
 
   // Following Parameters are added (20230309)
   wl = parameters->Get_wl();
   cl = parameters->Get_cl();
   rho0 = parameters->Get_rho0();
   gamm = parameters->Get_gamm();
   eta = parameters->Get_eta(); // Skyrme-QMD
   kappas = parameters->Get_kappas(); // Skyrme-QMD
    
   cpw = parameters->Get_cpw();
   cph = parameters->Get_cph();
   cpc = parameters->Get_cpc();
    
   c0 = parameters->Get_c0();
   c3 = parameters->Get_c3();
   cs = parameters->Get_cs();
   g0 = parameters->Get_g0(); // Skyrme-QMD
   g0iso = parameters->Get_g0iso(); // Skyrme-QMD
   gtau0 = parameters->Get_gtau0(); // Skyrme-QMD
    
   // distance
   c0w = 1.0/4.0/wl;
   //c3w = 1.0/4.0/wl; //no need
   c0sw = std::sqrt( c0w );
   clw = 2.0 / std::sqrt ( 4.0 * pi * wl );
    
   // graduate
   c0g = - c0 / ( 2.0 * wl );
   c3g = - c3 / ( 4.0 * wl ) * gamm;
   csg = - cs / ( 2.0 * wl );
   pag = gamm - 1;
   pag_tau = eta - 1; // Skyrme-QMD
   cg0 = - g0 / ( 2.0 * wl );  // Skyrme-QMD
   cgtau0 = - gtau0 / ( 4.0 * wl ) * eta;  // Skyrme-QMD
}

Member Function Documentation

CalEnergyAndAngularMomentumInCM()

void G4LightIonQMDNucleus::CalEnergyAndAngularMomentumInCM

(

)

Definition at line 180 of file G4LightIonQMDNucleus.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 );
 
   // binding energy should be evaluated with a relativistic version: 20230308 by Y-H. Sato and A. Haga
   for ( G4int i= 0; i < n ; i++ )
   {
      G4ThreeVector ri = GetParticipant( i )->GetPosition();
      G4double trans = gamma / ( gamma + 1.0 ) * ri * beta;
      G4double nucpote = GetNuclPotential( i );
 
      es[i] = std::sqrt ( G4Pow::GetInstance()->powN ( GetParticipant( i )->GetMass() , 2 ) + pcm[i]*pcm[i] + 2.0*GetParticipant( i )->GetMass()*nucpote) - GetParticipant( i )->GetMass(); //R-JQMD
 
      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
   // relativistic version Y-H. Sato and A. Haga 20230309
   G4double bindingEnergy =  ( std::accumulate ( es.begin() , es.end() , 0.0 ) );
 
   //G4cout << "n " << n << "totalpote " << totalpote << " " << potentialEnergy << " " << bindingEnergy << G4endl;
   //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;
   if ( excitationEnergy < 0 ) excitationEnergy = 0.0;
 
 }

Referenced by G4LightIonQMDReaction::ApplyYourself(), and G4LightIonQMDMeanField::SetNucleus().

Get4Momentum()

G4LorentzVector G4LightIonQMDNucleus::Get4Momentum

(

)

Definition at line 98 of file G4LightIonQMDNucleus.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 G4LightIonQMDNucleus::GetAngularMomentum ( )

inline

Definition at line 66 of file G4LightIonQMDNucleus.hh.

{ return jj; };

GetAtomicNumber()

G4int G4LightIonQMDNucleus::GetAtomicNumber

(

)

Definition at line 131 of file G4LightIonQMDNucleus.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 G4LightIonQMDNucleus::GetExcitationEnergy ( )

inline

Definition at line 64 of file G4LightIonQMDNucleus.hh.

{ return excitationEnergy; };

GetMassNumber()

G4int G4LightIonQMDNucleus::GetMassNumber

(

)

Definition at line 110 of file G4LightIonQMDNucleus.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__, "G4LightIonQMDNucleus has the mass number of 0!");
   }
 
   return A;
}

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

GetNuclearMass()

G4double G4LightIonQMDNucleus::GetNuclearMass

(

)

Definition at line 145 of file G4LightIonQMDNucleus.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 G4LightIonQMDNucleus::SetTotalPotential ( G4double x )

inline

Definition at line 63 of file G4LightIonQMDNucleus.hh.

{ potentialEnergy = x; };

Referenced by G4LightIonQMDReaction::ApplyYourself(), and G4LightIonQMDMeanField::SetNucleus().

---

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

• G4LightIonQMDNucleus.hh
• G4LightIonQMDNucleus.cc