G4QMDNucleus Class Reference

#include <G4QMDNucleus.hh>

Inheritance diagram for G4QMDNucleus:

Collaboration 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

(

)

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

Here is the call graph for this function:

Member Function Documentation

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

// Angluar momentum

   G4ThreeVector rl ( 0.0 ); 
   for ( G4int i= 0; i < n ; i++ ) 
   {
      rl += rcm[i].cross ( pcm[i] );
   }

   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; 

}

Here is the call graph for this function:

Here is the caller graph for this function:

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

Here is the caller graph for this function:

G4int G4QMDNucleus::GetAngularMomentum ( )

inline

Definition at line 65 of file G4QMDNucleus.hh.

{ return jj; };

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

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4QMDNucleus::GetExcitationEnergy ( void  )

inline

Definition at line 63 of file G4QMDNucleus.hh.

{ return excitationEnergy; };

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

Here is the call graph for this function:

Here is the caller graph for this function:

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

Here is the call graph for this function:

void G4QMDNucleus::SetTotalPotential ( G4double  x )

inline

Definition at line 62 of file G4QMDNucleus.hh.

{ potentialEnergy = x; };

Here is the caller graph for this function:

---

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

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