G4QMDGroundStateNucleus Class Reference

#include <G4QMDGroundStateNucleus.hh>

Inheritance diagram for G4QMDGroundStateNucleus:

Public Member Functions
	G4QMDGroundStateNucleus ()
	G4QMDGroundStateNucleus (G4int z, G4int a)
	~G4QMDGroundStateNucleus ()
Public Member Functions inherited from G4QMDNucleus
	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 45 of file G4QMDGroundStateNucleus.hh.

Constructor & Destructor Documentation

G4QMDGroundStateNucleus() [1/2]

G4QMDGroundStateNucleus::G4QMDGroundStateNucleus

(

)

G4QMDGroundStateNucleus() [2/2]

G4QMDGroundStateNucleus::G4QMDGroundStateNucleus	(	G4int	z,
		G4int	a )

Definition at line 38 of file G4QMDGroundStateNucleus.cc.

: maxTrial ( 1000 )
, r00 ( 1.124 )  // radius parameter for Woods-Saxon [fm] 
, r01 ( 0.5 )    // radius parameter for Woods-Saxon 
, saa ( 0.2 )    // diffuse parameter for initial Woods-Saxon shape
, rada ( 0.9 )   // cutoff parameter
, radb ( 0.3 )   // cutoff parameter
, dsam ( 1.5 )   // minimum distance for same particle [fm]
, ddif ( 1.0 )   // minimum distance for different particle
, edepth ( 0.0 )
, epse ( 0.000001 )  // torelance for energy in [GeV]
, meanfield ( NULL ) 
{
 
   //std::cout << " G4QMDGroundStateNucleus( G4int z , G4int a ) Begin " << z << " " << a << std::endl;
 
   dsam2 = dsam*dsam;
   ddif2 = ddif*ddif;
 
   G4QMDParameters* parameters = G4QMDParameters::GetInstance();
 
   hbc = parameters->Get_hbc();
   gamm = parameters->Get_gamm();
   cpw = parameters->Get_cpw();
   cph = parameters->Get_cph();
   epsx = parameters->Get_epsx();
   cpc = parameters->Get_cpc();
 
   cdp = parameters->Get_cdp();
   c0p = parameters->Get_c0p();
   c3p = parameters->Get_c3p();
   csp = parameters->Get_csp();
   clp = parameters->Get_clp();
 
   ebini = 0.0;
   // Following 10 lines should be here, right before the line 90.
   // Otherwise, mass number cannot be conserved if the projectile or
   // the target are nucleons.
   //Nucleon primary or target case;
   if ( z == 1 && a == 1 ) {  // Hydrogen  Case or proton primary 
      SetParticipant( new G4QMDParticipant( G4Proton::Proton() , G4ThreeVector( 0.0 ) , G4ThreeVector( 0.0 ) ) );
//      ebini = 0.0; 
      return;
   }
   else if ( z == 0 && a == 1 ) { // Neutron primary 
      SetParticipant( new G4QMDParticipant( G4Neutron::Neutron() , G4ThreeVector( 0.0 ) , G4ThreeVector( 0.0 ) ) );
//      ebini = 0.0; 
      return;
   }
 
 
   //edepth = 0.0; 
 
   for ( int i = 0 ; i < a ; i++ )
   {
 
      G4ParticleDefinition* pd; 
 
      if ( i < z )
      { 
         pd = G4Proton::Proton();
      }
      else
      {
         pd = G4Neutron::Neutron();
      }
         
      G4ThreeVector p( 0.0 );
      G4ThreeVector r( 0.0 );
      G4QMDParticipant* aParticipant = new G4QMDParticipant( pd , p , r );
      SetParticipant( aParticipant );
 
   }
 
   G4double radious = r00 * G4Pow::GetInstance()->A13( double ( GetMassNumber() ) ); 
 
   rt00 = radious - r01; 
   radm = radious - rada * ( gamm - 1.0 ) + radb;
   rmax = 1.0 / ( 1.0 + G4Exp ( -rt00/saa ) );
 
   //maxTrial = 1000;
   
   
   meanfield = new G4QMDMeanField();
   meanfield->SetSystem( this );
 
   //std::cout << "G4QMDGroundStateNucleus( G4int z , G4int a ) packNucleons Begin ( z , a ) ( " << z << ", " << a << " )" << std::endl;
   packNucleons();
   //std::cout << "G4QMDGroundStateNucleus( G4int z , G4int a ) packNucleons End" << std::endl;
 
   delete meanfield;
 
}

~G4QMDGroundStateNucleus()

G4QMDGroundStateNucleus::~G4QMDGroundStateNucleus ( )

inline

Definition at line 50 of file G4QMDGroundStateNucleus.hh.

      {
         rho_l.clear();  
         d_pot.clear();
      };

---

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

• G4QMDGroundStateNucleus.hh
• G4QMDGroundStateNucleus.cc