G4Fancy3DNucleus Class Reference

#include <G4Fancy3DNucleus.hh>

Inheritance diagram for G4Fancy3DNucleus:

Public Member Functions
	G4Fancy3DNucleus ()
	~G4Fancy3DNucleus ()
void	Init (G4int theA, G4int theZ, G4int numberOfLambdas=0)
G4bool	StartLoop ()
G4Nucleon *	GetNextNucleon ()
const std::vector< G4Nucleon > &	GetNucleons ()
G4int	GetMassNumber ()
G4double	GetMass ()
G4int	GetCharge ()
G4int	GetNumberOfLambdas ()
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::G4Fancy3DNucleus

(

)

Definition at line 60 of file G4Fancy3DNucleus.cc.

  : myA(0), myZ(0), myL(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 67 of file G4Fancy3DNucleus.cc.

{
  if(theDensity) delete theDensity;
}

Member Function Documentation

AddExcitationEnergy()

G4double G4Fancy3DNucleus::AddExcitationEnergy ( G4double anE )

inline

Definition at line 138 of file G4Fancy3DNucleus.hh.

{
   excitationEnergy +=anE;
   return excitationEnergy;
}

CenterNucleons()

void G4Fancy3DNucleus::CenterNucleons

(

)

Definition at line 251 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 607 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 213 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 220 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 241 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 227 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 263 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 123 of file G4Fancy3DNucleus.hh.

{
    return myZ;
}

GetExcitationEnergy()

G4double G4Fancy3DNucleus::GetExcitationEnergy ( )

inline

Definition at line 144 of file G4Fancy3DNucleus.hh.

{
   return excitationEnergy;
}

GetMass()

G4double G4Fancy3DNucleus::GetMass ( )

virtual

Implements G4V3DNucleus.

Definition at line 203 of file G4Fancy3DNucleus.cc.

{
  if ( myL <= 0 ) return myZ*G4Proton::Proton()->GetPDGMass() + 
                         (myA-myZ)*G4Neutron::Neutron()->GetPDGMass() -
                         BindingEnergy();
  else            return G4HyperNucleiProperties::GetNuclearMass(myA, myZ, myL); 
}

GetMassNumber()

G4int G4Fancy3DNucleus::GetMassNumber ( )

inlinevirtual

Implements G4V3DNucleus.

Definition at line 128 of file G4Fancy3DNucleus.hh.

{
    return myA;
}

GetNextNucleon()

G4Nucleon * G4Fancy3DNucleus::GetNextNucleon ( )

virtual

Implements G4V3DNucleus.

Definition at line 138 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 273 of file G4Fancy3DNucleus.cc.

{
    return theDensity;
}

GetNuclearRadius() [1/2]

G4double G4Fancy3DNucleus::GetNuclearRadius ( )

virtual

Implements G4V3DNucleus.

Definition at line 179 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 184 of file G4Fancy3DNucleus.cc.

{   
    return theDensity->GetRadius(maxRelativeDensity);
}

GetNucleons()

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

virtual

Implements G4V3DNucleus.

Definition at line 144 of file G4Fancy3DNucleus.cc.

{
  return theNucleons;
}

Referenced by G4MuMinusCapturePrecompound::ApplyYourself().

GetNumberOfLambdas()

G4int G4Fancy3DNucleus::GetNumberOfLambdas ( )

inlinevirtual

Implements G4V3DNucleus.

Definition at line 133 of file G4Fancy3DNucleus.hh.

{
        return myL;
} 

GetOuterRadius()

G4double G4Fancy3DNucleus::GetOuterRadius ( )

virtual

Implements G4V3DNucleus.

Definition at line 189 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, G4int numberOfLambdas = 0 )

virtual

Implements G4V3DNucleus.

Definition at line 83 of file G4Fancy3DNucleus.cc.

{
  currentNucleon=-1;
  theNucleons.clear();
  nucleondistance = 0.8*fermi;
  places.clear();
  momentum.clear();
  fermiM.clear();
  testSums.clear();
 
  myZ = theZ;
  myA = theA;
  myL = std::max(numberOfLambdas, 0);  // Cannot be negative
  excitationEnergy=0;
 
  theNucleons.resize(myA);  // Pre-loads vector with empty elements
 
  // For simplicity, we neglect eventual Lambdas in the nucleus as far as the
  // density of nucler levels and the Fermi level are concerned.
  
  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 G4LowEIonFragmentation::ApplyYourself(), and G4MuMinusCapturePrecompound::ApplyYourself().

SortNucleonsDecZ()

void G4Fancy3DNucleus::SortNucleonsDecZ ( )

virtual

Implements G4V3DNucleus.

Definition at line 164 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 156 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 131 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