G4Fancy3DNucleus Class Reference

#include <G4Fancy3DNucleus.hh>

Inheritance diagram for G4Fancy3DNucleus:

Collaboration diagram for G4Fancy3DNucleus:

Public Member Functions
	G4Fancy3DNucleus ()
	~G4Fancy3DNucleus ()
void	Init (G4int theA, G4int theZ)
G4bool	StartLoop ()
G4Nucleon *	GetNextNucleon ()
const std::vector< G4Nucleon > &	GetNucleons ()
G4int	GetMassNumber ()
G4double	GetMass ()
G4int	GetCharge ()
G4double	GetNuclearRadius ()
G4double	GetNuclearRadius (const G4double maxRelativeDensity)
G4double	GetOuterRadius ()
G4double	AddExcitationEnergy (G4double)
G4double	GetExcitationEnergy ()
G4double	CoulombBarrier ()
void	DoLorentzBoost (const G4LorentzVector &theBoost)
void	DoLorentzBoost (const G4ThreeVector &theBeta)
void	DoLorentzContraction (const G4LorentzVector &theBoost)
void	DoLorentzContraction (const G4ThreeVector &theBeta)
void	CenterNucleons ()
void	DoTranslation (const G4ThreeVector &theShift)
const G4VNuclearDensity *	GetNuclearDensity () const
void	SortNucleonsIncZ ()
void	SortNucleonsDecZ ()
Public Member Functions inherited from G4V3DNucleus
	G4V3DNucleus ()
virtual	~G4V3DNucleus ()
std::pair< G4double, G4double >	ChooseImpactXandY (G4double maxImpact)
std::pair< G4double, G4double >	RefetchImpactXandY ()

Detailed Description

Definition at line 54 of file G4Fancy3DNucleus.hh.

Constructor & Destructor Documentation

G4Fancy3DNucleus::G4Fancy3DNucleus

(

)

Definition at line 53 of file G4Fancy3DNucleus.cc.

  : myA(0), myZ(0), theNucleons(250), currentNucleon(-1), theDensity(0), 
    nucleondistance(0.8*fermi),excitationEnergy(0.),
    places(250), momentum(250), fermiM(250), testSums(250)
{
//G4cout <<"G4Fancy3DNucleus::G4Fancy3DNucleus()"<<G4endl;
}

G4Fancy3DNucleus::~G4Fancy3DNucleus

(

)

Definition at line 61 of file G4Fancy3DNucleus.cc.

{
  if(theDensity) delete theDensity;
}

Member Function Documentation

G4double G4Fancy3DNucleus::AddExcitationEnergy ( G4double  anE )

inline

Definition at line 130 of file G4Fancy3DNucleus.hh.

{
   excitationEnergy +=anE;
   return excitationEnergy;
}

void G4Fancy3DNucleus::CenterNucleons

(

)

Definition at line 238 of file G4Fancy3DNucleus.cc.

{
    G4ThreeVector   center;
    
    for (G4int i=0; i<myA; i++ )
    {
       center+=theNucleons[i].GetPosition();
    }   
    center /= -myA;
    DoTranslation(center);
}

Here is the call graph for this function:

G4double G4Fancy3DNucleus::CoulombBarrier ( )

virtual

Implements G4V3DNucleus.

Definition at line 501 of file G4Fancy3DNucleus.cc.

{
  static const G4double cfactor = (1.44/1.14) * MeV;
  return cfactor*myZ/(1.0 + G4Pow::GetInstance()->Z13(myA));
}

Here is the call graph for this function:

void G4Fancy3DNucleus::DoLorentzBoost ( const G4LorentzVector &  theBoost )

virtual

Implements G4V3DNucleus.

Definition at line 200 of file G4Fancy3DNucleus.cc.

{
    for (G4int i=0; i<myA; i++){
        theNucleons[i].Boost(theBoost);
    }
}

void G4Fancy3DNucleus::DoLorentzBoost ( const G4ThreeVector &  theBeta )

virtual

Implements G4V3DNucleus.

Definition at line 207 of file G4Fancy3DNucleus.cc.

{
    for (G4int i=0; i<myA; i++){
        theNucleons[i].Boost(theBeta);
    }
}

void G4Fancy3DNucleus::DoLorentzContraction ( const G4LorentzVector &  theBoost )

virtual

Implements G4V3DNucleus.

Definition at line 228 of file G4Fancy3DNucleus.cc.

{
    if (theBoost.e() !=0 ) {
       G4ThreeVector beta = theBoost.vect()/theBoost.e();
       DoLorentzContraction(beta);
    }
}

Here is the call graph for this function:

void G4Fancy3DNucleus::DoLorentzContraction ( const G4ThreeVector &  theBeta )

virtual

Implements G4V3DNucleus.

Definition at line 214 of file G4Fancy3DNucleus.cc.

{
    G4double beta2=theBeta.mag2();
    if (beta2 > 0) {
       G4double factor=(1-std::sqrt(1-beta2))/beta2; // (gamma-1)/gamma/beta**2
       G4ThreeVector rprime;
       for (G4int i=0; i< myA; i++) {
         rprime = theNucleons[i].GetPosition() - 
           factor * (theBeta*theNucleons[i].GetPosition()) * theBeta;  
         theNucleons[i].SetPosition(rprime);
       }
    }    
}

Here is the call graph for this function:

void G4Fancy3DNucleus::DoTranslation ( const G4ThreeVector &  theShift )

virtual

Implements G4V3DNucleus.

Definition at line 250 of file G4Fancy3DNucleus.cc.

{
  G4ThreeVector tempV;
  for (G4int i=0; i<myA; i++ )
    {
      tempV = theNucleons[i].GetPosition() + theShift;
      theNucleons[i].SetPosition(tempV);
    }   
}

Here is the caller graph for this function:

G4int G4Fancy3DNucleus::GetCharge ( )

inlinevirtual

Implements G4V3DNucleus.

Definition at line 121 of file G4Fancy3DNucleus.hh.

{
    return myZ;
}

G4double G4Fancy3DNucleus::GetExcitationEnergy ( void  )

inline

Definition at line 136 of file G4Fancy3DNucleus.hh.

{
   return excitationEnergy;
}

G4double G4Fancy3DNucleus::GetMass ( )

virtual

Implements G4V3DNucleus.

Definition at line 191 of file G4Fancy3DNucleus.cc.

{
        return   myZ*G4Proton::Proton()->GetPDGMass() + 
                 (myA-myZ)*G4Neutron::Neutron()->GetPDGMass() -
                 BindingEnergy();
}

Here is the call graph for this function:

G4int G4Fancy3DNucleus::GetMassNumber ( )

inlinevirtual

Implements G4V3DNucleus.

Definition at line 126 of file G4Fancy3DNucleus.hh.

{
    return myA;
}

G4Nucleon * G4Fancy3DNucleus::GetNextNucleon ( )

virtual

Implements G4V3DNucleus.

Definition at line 126 of file G4Fancy3DNucleus.cc.

{
  return ( (currentNucleon>=0 && currentNucleon<myA) ? 
       &theNucleons[currentNucleon++] : 0 );
}

Here is the caller graph for this function:

const G4VNuclearDensity * G4Fancy3DNucleus::GetNuclearDensity ( ) const

virtual

Implements G4V3DNucleus.

Definition at line 260 of file G4Fancy3DNucleus.cc.

{
    return theDensity;
}

G4double G4Fancy3DNucleus::GetNuclearRadius ( )

virtual

Implements G4V3DNucleus.

Definition at line 167 of file G4Fancy3DNucleus.cc.

{
    return GetNuclearRadius(0.5);
}

G4double G4Fancy3DNucleus::GetNuclearRadius ( const G4double  maxRelativeDensity )

virtual

Implements G4V3DNucleus.

Definition at line 172 of file G4Fancy3DNucleus.cc.

{   
    return theDensity->GetRadius(maxRelativeDensity);
}

Here is the call graph for this function:

const std::vector< G4Nucleon > & G4Fancy3DNucleus::GetNucleons ( )

virtual

Implements G4V3DNucleus.

Definition at line 132 of file G4Fancy3DNucleus.cc.

{
  return theNucleons;
}

Here is the caller graph for this function:

G4double G4Fancy3DNucleus::GetOuterRadius ( )

virtual

Implements G4V3DNucleus.

Definition at line 177 of file G4Fancy3DNucleus.cc.

{
    G4double maxradius2=0;
    
    for (int i=0; i<myA; i++)
    {
       if ( theNucleons[i].GetPosition().mag2() > maxradius2 )
       {
        maxradius2=theNucleons[i].GetPosition().mag2();
       }
    }
    return std::sqrt(maxradius2)+nucleondistance;
}

Here is the caller graph for this function:

void G4Fancy3DNucleus::Init ( G4int  theA, G4int  theZ  )

virtual

Implements G4V3DNucleus.

Definition at line 77 of file G4Fancy3DNucleus.cc.

{
//  G4cout << "G4Fancy3DNucleus::Init(theA, theZ) called"<<G4endl;
  currentNucleon=-1;
  theNucleons.clear();

  myZ = theZ;
  myA= theA;
  excitationEnergy=0;

  theNucleons.resize(myA);  // Pre-loads vector with empty elements
  
//  G4cout << "myA, myZ" << myA << ", " << myZ << G4endl;

  if(theDensity) delete theDensity;
  if ( myA < 17 ) {
     theDensity = new G4NuclearShellModelDensity(myA, myZ);
  } else {
     theDensity = new G4NuclearFermiDensity(myA, myZ);
  }

  theFermi.Init(myA, myZ);
  
  ChooseNucleons();
  
  ChoosePositions();
  
//  CenterNucleons();   // This would introduce a bias 

  ChooseFermiMomenta();
  
  G4double Ebinding= BindingEnergy()/myA;
  
  for (G4int aNucleon=0; aNucleon < myA; aNucleon++)
  {
    theNucleons[aNucleon].SetBindingEnergy(Ebinding);
  }
  
  
  return;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void G4Fancy3DNucleus::SortNucleonsDecZ ( )

virtual

Implements G4V3DNucleus.

Definition at line 152 of file G4Fancy3DNucleus.cc.

{
  if (theNucleons.size() < 2 ) return;      // Avoid unnecessary work
  SortNucleonsIncZ();

  std::reverse(theNucleons.begin(), theNucleons.end());
}

Here is the call graph for this function:

void G4Fancy3DNucleus::SortNucleonsIncZ ( )

virtual

Implements G4V3DNucleus.

Definition at line 144 of file G4Fancy3DNucleus.cc.

{
  if (theNucleons.size() < 2 ) return;      // Avoid unnecesary work

  std::sort(theNucleons.begin(), theNucleons.end(),
        G4Fancy3DNucleusHelperForSortInZ); 
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4bool G4Fancy3DNucleus::StartLoop ( )

virtual

Implements G4V3DNucleus.

Definition at line 119 of file G4Fancy3DNucleus.cc.

{
    currentNucleon=0;
    return (theNucleons.size()>0);
}   

Here is the caller graph for this function:

---

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

• source/geant4.10.03.p03/source/processes/hadronic/models/util/include/G4Fancy3DNucleus.hh
• source/geant4.10.03.p03/source/processes/hadronic/models/util/src/G4Fancy3DNucleus.cc