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 ()
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 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::~G4Fancy3DNucleus

(

)

Definition at line 60 of file G4Fancy3DNucleus.cc.

{
  if(theDensity) delete theDensity;
}

Member Function Documentation

CenterNucleons()

void G4Fancy3DNucleus::CenterNucleons

(

)

Definition at line 232 of file G4Fancy3DNucleus.cc.

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

CoulombBarrier()

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

DoLorentzBoost() [1/2]

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

DoLorentzBoost() [2/2]

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

DoLorentzContraction() [1/2]

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

Referenced by DoLorentzContraction().

DoLorentzContraction() [2/2]

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

DoTranslation()

void G4Fancy3DNucleus::DoTranslation ( const G4ThreeVector &  theShift )

virtual

Implements G4V3DNucleus.

Definition at line 244 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 118 of file G4Fancy3DNucleus.hh.

{
    return myZ;
}

GetMass()

G4double G4Fancy3DNucleus::GetMass ( )

virtual

Implements G4V3DNucleus.

Definition at line 189 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 123 of file G4Fancy3DNucleus.hh.

{
    return myA;
}

GetNextNucleon()

G4Nucleon * G4Fancy3DNucleus::GetNextNucleon ( )

virtual

Implements G4V3DNucleus.

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

{
    return theDensity;
}

GetNuclearRadius() [1/2]

G4double G4Fancy3DNucleus::GetNuclearRadius ( )

virtual

Implements G4V3DNucleus.

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

{   
    return theDensity->GetRadius(maxRelativeDensity);
}

GetNucleons()

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

virtual

Implements G4V3DNucleus.

Definition at line 130 of file G4Fancy3DNucleus.cc.

{
  return theNucleons;
}

Referenced by G4MuMinusCapturePrecompound::ApplyYourself().

GetOuterRadius()

G4double G4Fancy3DNucleus::GetOuterRadius ( )

virtual

Implements G4V3DNucleus.

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

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

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

SortNucleonsDecZ()

void G4Fancy3DNucleus::SortNucleonsDecZ ( )

virtual

Implements G4V3DNucleus.

Definition at line 150 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 142 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 117 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