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	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 ()

---

Detailed Description

Definition at line 54 of file G4Fancy3DNucleus.hh.

---

Constructor & Destructor Documentation

G4Fancy3DNucleus::G4Fancy3DNucleus

(

)

Definition at line 52 of file G4Fancy3DNucleus.cc.

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

G4Fancy3DNucleus::~G4Fancy3DNucleus

(

)

Definition at line 60 of file G4Fancy3DNucleus.cc.

{
  if(theDensity) delete theDensity;
}

---

Member Function Documentation

void G4Fancy3DNucleus::CenterNucleons

(

)

Definition at line 232 of file G4Fancy3DNucleus.cc.

References DoTranslation().

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

G4double G4Fancy3DNucleus::CoulombBarrier

(

)

[virtual]

Implements G4V3DNucleus.

Definition at line 490 of file G4Fancy3DNucleus.cc.

{
  G4double coulombBarrier = (1.44/1.14) * MeV * myZ / (1.0 + std::pow(G4double(myA),1./3.));
  return coulombBarrier;
}

void G4Fancy3DNucleus::DoLorentzBoost

(

const G4ThreeVector &

theBeta

)

[virtual]

Implements G4V3DNucleus.

Definition at line 205 of file G4Fancy3DNucleus.cc.

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

void G4Fancy3DNucleus::DoLorentzBoost

(

const G4LorentzVector &

theBoost

)

[virtual]

Implements G4V3DNucleus.

Definition at line 198 of file G4Fancy3DNucleus.cc.

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

void G4Fancy3DNucleus::DoLorentzContraction

(

const G4ThreeVector &

theBeta

)

[virtual]

Implements G4V3DNucleus.

Definition at line 212 of file G4Fancy3DNucleus.cc.

{
        G4double factor=(1-std::sqrt(1-theBeta.mag2()))/theBeta.mag2(); // (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);
        }    
}

void G4Fancy3DNucleus::DoLorentzContraction

(

const G4LorentzVector &

theBoost

)

[virtual]

Implements G4V3DNucleus.

Definition at line 223 of file G4Fancy3DNucleus.cc.

{
        G4ThreeVector beta = theBoost.vect()/theBoost.e();
        // DoLorentzBoost(beta); 
        DoLorentzContraction(beta);
}

void G4Fancy3DNucleus::DoTranslation

(

const G4ThreeVector &

theShift

)

[virtual]

Implements G4V3DNucleus.

Definition at line 244 of file G4Fancy3DNucleus.cc.

Referenced by CenterNucleons().

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

G4int G4Fancy3DNucleus::GetCharge

(

)

[inline, virtual]

Implements G4V3DNucleus.

Definition at line 118 of file G4Fancy3DNucleus.hh.

{
        return myZ;
}

G4double G4Fancy3DNucleus::GetMass

(

)

[virtual]

Implements G4V3DNucleus.

Definition at line 189 of file G4Fancy3DNucleus.cc.

References G4ParticleDefinition::GetPDGMass(), G4Neutron::Neutron(), and G4Proton::Proton().

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

G4int G4Fancy3DNucleus::GetMassNumber

(

)

[inline, virtual]

Implements G4V3DNucleus.

Definition at line 123 of file G4Fancy3DNucleus.hh.

{
        return myA;
}

G4Nucleon * G4Fancy3DNucleus::GetNextNucleon

(

)

[virtual]

Implements G4V3DNucleus.

Definition at line 124 of file G4Fancy3DNucleus.cc.

Referenced by G4LowEIonFragmentation::ApplyYourself().

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

const G4VNuclearDensity * G4Fancy3DNucleus::GetNuclearDensity

(

)

const [virtual]

Implements G4V3DNucleus.

Definition at line 254 of file G4Fancy3DNucleus.cc.

{
        return theDensity;
}

G4double G4Fancy3DNucleus::GetNuclearRadius

(

const G4double

maxRelativeDensity

)

[virtual]

Implements G4V3DNucleus.

Definition at line 170 of file G4Fancy3DNucleus.cc.

References G4VNuclearDensity::GetRadius().

{       
        return theDensity->GetRadius(maxRelativeDensity);
}

G4double G4Fancy3DNucleus::GetNuclearRadius

(

)

[virtual]

Implements G4V3DNucleus.

Definition at line 165 of file G4Fancy3DNucleus.cc.

{
        return GetNuclearRadius(0.5);
}

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

(

)

[virtual]

Implements G4V3DNucleus.

Definition at line 130 of file G4Fancy3DNucleus.cc.

Referenced by G4MuMinusCapturePrecompound::ApplyYourself().

{
  return theNucleons;
}

G4double G4Fancy3DNucleus::GetOuterRadius

(

)

[virtual]

Implements G4V3DNucleus.

Definition at line 175 of file G4Fancy3DNucleus.cc.

Referenced by G4LowEIonFragmentation::ApplyYourself().

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

void G4Fancy3DNucleus::Init	(	G4int	theA,
		G4int	theZ
	)			[virtual]

Implements G4V3DNucleus.

Definition at line 76 of file G4Fancy3DNucleus.cc.

References G4FermiMomentum::Init().

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

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

  myZ = theZ;
  myA= theA;

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

void G4Fancy3DNucleus::SortNucleonsDecZ

(

)

[virtual]

Implements G4V3DNucleus.

Definition at line 150 of file G4Fancy3DNucleus.cc.

References SortNucleonsIncZ().

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

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

void G4Fancy3DNucleus::SortNucleonsIncZ

(

)

[virtual]

Implements G4V3DNucleus.

Definition at line 142 of file G4Fancy3DNucleus.cc.

References G4Fancy3DNucleusHelperForSortInZ().

Referenced by SortNucleonsDecZ().

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

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

G4bool G4Fancy3DNucleus::StartLoop

(

)

[virtual]

Implements G4V3DNucleus.

Definition at line 117 of file G4Fancy3DNucleus.cc.

Referenced by G4LowEIonFragmentation::ApplyYourself().

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

---

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

• G4Fancy3DNucleus.hh
• G4Fancy3DNucleus.cc

---