G4Fancy3DNucleus Class Reference

#include <G4Fancy3DNucleus.hh>

Inheritance 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 ()
virtual void	Init (G4int theA, G4int theZ)=0
virtual G4bool	StartLoop ()=0
virtual G4Nucleon *	GetNextNucleon ()=0
virtual const std::vector< G4Nucleon > &	GetNucleons ()=0
virtual G4int	GetMassNumber ()=0
virtual G4double	GetMass ()=0
virtual G4int	GetCharge ()=0
virtual G4double	GetNuclearRadius ()=0
virtual G4double	GetNuclearRadius (const G4double maxRelativeDensity)=0
virtual G4double	GetOuterRadius ()=0
virtual G4double	CoulombBarrier ()=0
virtual void	DoLorentzBoost (const G4LorentzVector &theBoost)=0
virtual void	DoLorentzBoost (const G4ThreeVector &theBeta)=0
virtual void	DoLorentzContraction (const G4LorentzVector &theBoost)=0
virtual void	DoLorentzContraction (const G4ThreeVector &theBeta)=0
virtual void	DoTranslation (const G4ThreeVector &theShift)=0
virtual const G4VNuclearDensity *	GetNuclearDensity () const =0
virtual void	SortNucleonsIncZ ()=0
virtual void	SortNucleonsDecZ ()=0
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::G4Fancy3DNucleus

(

)

Definition at line 59 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)
{
}

~G4Fancy3DNucleus()

G4Fancy3DNucleus::~G4Fancy3DNucleus

(

)

Definition at line 66 of file G4Fancy3DNucleus.cc.

{
  if(theDensity) delete theDensity;
}

Member Function Documentation

AddExcitationEnergy()

G4double G4Fancy3DNucleus::AddExcitationEnergy ( G4double  anE )

inline

Definition at line 131 of file G4Fancy3DNucleus.hh.

{
   excitationEnergy +=anE;
   return excitationEnergy;
}

CenterNucleons()

void G4Fancy3DNucleus::CenterNucleons

(

)

Definition at line 245 of file G4Fancy3DNucleus.cc.

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

Referenced by Init().

CoulombBarrier()

G4double G4Fancy3DNucleus::CoulombBarrier ( )

virtual

Implements G4V3DNucleus.

Definition at line 595 of file G4Fancy3DNucleus.cc.

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

DoLorentzBoost() [1/2]

void G4Fancy3DNucleus::DoLorentzBoost ( const G4LorentzVector &  theBoost )

virtual

Implements G4V3DNucleus.

Definition at line 207 of file G4Fancy3DNucleus.cc.

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

DoLorentzBoost() [2/2]

void G4Fancy3DNucleus::DoLorentzBoost ( const G4ThreeVector &  theBeta )

virtual

Implements G4V3DNucleus.

Definition at line 214 of file G4Fancy3DNucleus.cc.

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

DoLorentzContraction() [1/2]

void G4Fancy3DNucleus::DoLorentzContraction ( const G4LorentzVector &  theBoost )

virtual

Implements G4V3DNucleus.

Definition at line 235 of file G4Fancy3DNucleus.cc.

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

Referenced by DoLorentzContraction().

DoLorentzContraction() [2/2]

void G4Fancy3DNucleus::DoLorentzContraction ( const G4ThreeVector &  theBeta )

virtual

Implements G4V3DNucleus.

Definition at line 221 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);
       }
    }    
}

DoTranslation()

void G4Fancy3DNucleus::DoTranslation ( const G4ThreeVector &  theShift )

virtual

Implements G4V3DNucleus.

Definition at line 257 of file G4Fancy3DNucleus.cc.

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

Referenced by CenterNucleons().

GetCharge()

G4int G4Fancy3DNucleus::GetCharge ( )

inlinevirtual

Implements G4V3DNucleus.

Definition at line 122 of file G4Fancy3DNucleus.hh.

{
    return myZ;
}

GetExcitationEnergy()

G4double G4Fancy3DNucleus::GetExcitationEnergy ( )

inline

Definition at line 137 of file G4Fancy3DNucleus.hh.

{
   return excitationEnergy;
}

GetMass()

G4double G4Fancy3DNucleus::GetMass ( )

virtual

Implements G4V3DNucleus.

Definition at line 198 of file G4Fancy3DNucleus.cc.

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

GetMassNumber()

G4int G4Fancy3DNucleus::GetMassNumber ( )

inlinevirtual

Implements G4V3DNucleus.

Definition at line 127 of file G4Fancy3DNucleus.hh.

{
    return myA;
}

GetNextNucleon()

G4Nucleon * G4Fancy3DNucleus::GetNextNucleon ( )

virtual

Implements G4V3DNucleus.

Definition at line 133 of file G4Fancy3DNucleus.cc.

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

Referenced by G4LowEIonFragmentation::ApplyYourself().

GetNuclearDensity()

const G4VNuclearDensity * G4Fancy3DNucleus::GetNuclearDensity ( ) const

virtual

Implements G4V3DNucleus.

Definition at line 267 of file G4Fancy3DNucleus.cc.

{
    return theDensity;
}

GetNuclearRadius() [1/2]

G4double G4Fancy3DNucleus::GetNuclearRadius ( )

virtual

Implements G4V3DNucleus.

Definition at line 174 of file G4Fancy3DNucleus.cc.

{
    return GetNuclearRadius(0.5);
}

Referenced by GetNuclearRadius().

GetNuclearRadius() [2/2]

G4double G4Fancy3DNucleus::GetNuclearRadius ( const G4double  maxRelativeDensity )

virtual

Implements G4V3DNucleus.

Definition at line 179 of file G4Fancy3DNucleus.cc.

{   
    return theDensity->GetRadius(maxRelativeDensity);
}

GetNucleons()

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

virtual

Implements G4V3DNucleus.

Definition at line 139 of file G4Fancy3DNucleus.cc.

{
  return theNucleons;
}

Referenced by G4MuMinusCapturePrecompound::ApplyYourself().

GetOuterRadius()

G4double G4Fancy3DNucleus::GetOuterRadius ( )

virtual

Implements G4V3DNucleus.

Definition at line 184 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;
}

Referenced by G4LowEIonFragmentation::ApplyYourself().

Init()

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

virtual

Implements G4V3DNucleus.

Definition at line 82 of file G4Fancy3DNucleus.cc.

{
  currentNucleon=-1;
  theNucleons.clear();
  nucleondistance = 0.8*fermi;
  places.clear();
  momentum.clear();
  fermiM.clear();
  testSums.clear();
 
  myZ = theZ;
  myA= theA;
  excitationEnergy=0;
 
  theNucleons.resize(myA);  // Pre-loads vector with empty elements
  
  if(theDensity) delete theDensity;
  if ( myA < 17 ) {
     theDensity = new G4NuclearShellModelDensity(myA, myZ);
     if( myA == 12 ) nucleondistance=0.9*fermi; 
  } else {
     theDensity = new G4NuclearFermiDensity(myA, myZ);
  }
 
  theFermi.Init(myA, myZ);
  
  ChooseNucleons();
  
  ChoosePositions();
  
  if( myA == 12 ) CenterNucleons();   // This would introduce a bias
 
  ChooseFermiMomenta();
  
  G4double Ebinding= BindingEnergy()/myA;
  
  for (G4int aNucleon=0; aNucleon < myA; aNucleon++)
  {
    theNucleons[aNucleon].SetBindingEnergy(Ebinding);
  }
  
  return;
}

Referenced by G4MuMinusCapturePrecompound::ApplyYourself(), and G4LowEIonFragmentation::ApplyYourself().

SortNucleonsDecZ()

void G4Fancy3DNucleus::SortNucleonsDecZ ( )

virtual

Implements G4V3DNucleus.

Definition at line 159 of file G4Fancy3DNucleus.cc.

{
  if (theNucleons.size() < 2 ) return;      // Avoid unnecessary work
  SortNucleonsIncZ();
 
  std::reverse(theNucleons.begin(), theNucleons.end());
}

SortNucleonsIncZ()

void G4Fancy3DNucleus::SortNucleonsIncZ ( )

virtual

Implements G4V3DNucleus.

Definition at line 151 of file G4Fancy3DNucleus.cc.

{
  if (theNucleons.size() < 2 ) return;      // Avoid unnecesary work
 
  std::sort(theNucleons.begin(), theNucleons.end(),
        G4Fancy3DNucleusHelperForSortInZ); 
}

Referenced by SortNucleonsDecZ().

StartLoop()

G4bool G4Fancy3DNucleus::StartLoop ( )

virtual

Implements G4V3DNucleus.

Definition at line 126 of file G4Fancy3DNucleus.cc.

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

Referenced by G4LowEIonFragmentation::ApplyYourself().

---

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

• G4Fancy3DNucleus.hh
• G4Fancy3DNucleus.cc