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 ()
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 52 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 60 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 237 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 496 of file G4Fancy3DNucleus.cc.

References python.hepunit::MeV.

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

void G4Fancy3DNucleus::DoLorentzBoost ( const G4LorentzVector &  theBoost )

virtual

Implements G4V3DNucleus.

Definition at line 199 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 206 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 227 of file G4Fancy3DNucleus.cc.

References CLHEP::HepLorentzVector::e(), and CLHEP::HepLorentzVector::vect().

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

void G4Fancy3DNucleus::DoLorentzContraction ( const G4ThreeVector &  theBeta )

virtual

Implements G4V3DNucleus.

Definition at line 213 of file G4Fancy3DNucleus.cc.

References CLHEP::Hep3Vector::mag2().

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

void G4Fancy3DNucleus::DoTranslation ( const G4ThreeVector &  theShift )

virtual

Implements G4V3DNucleus.

Definition at line 249 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 ( )

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

inlinevirtual

Implements G4V3DNucleus.

Definition at line 126 of file G4Fancy3DNucleus.hh.

{
    return myA;
}

G4Nucleon * G4Fancy3DNucleus::GetNextNucleon ( )

virtual

Implements G4V3DNucleus.

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

{
    return theDensity;
}

G4double G4Fancy3DNucleus::GetNuclearRadius ( )

virtual

Implements G4V3DNucleus.

Definition at line 166 of file G4Fancy3DNucleus.cc.

{
    return GetNuclearRadius(0.5);
}

G4double G4Fancy3DNucleus::GetNuclearRadius ( const G4double  maxRelativeDensity )

virtual

Implements G4V3DNucleus.

Definition at line 171 of file G4Fancy3DNucleus.cc.

References G4VNuclearDensity::GetRadius().

{   
    return theDensity->GetRadius(maxRelativeDensity);
}

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

virtual

Implements G4V3DNucleus.

Definition at line 131 of file G4Fancy3DNucleus.cc.

Referenced by G4MuMinusCapturePrecompound::ApplyYourself().

{
  return theNucleons;
}

G4double G4Fancy3DNucleus::GetOuterRadius ( )

virtual

Implements G4V3DNucleus.

Definition at line 176 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 G4LowEIonFragmentation::ApplyYourself(), and G4MuMinusCapturePrecompound::ApplyYourself().

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

void G4Fancy3DNucleus::SortNucleonsDecZ ( )

virtual

Implements G4V3DNucleus.

Definition at line 151 of file G4Fancy3DNucleus.cc.

References reverse, and 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 143 of file G4Fancy3DNucleus.cc.

References G4Fancy3DNucleusHelperForSortInZ(), and sort().

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 118 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