G4EvaporationChannel Class Reference

#include <G4EvaporationChannel.hh>

Inheritance diagram for G4EvaporationChannel:

Collaboration diagram for G4EvaporationChannel:

Public Member Functions
	G4EvaporationChannel (G4int theA, G4int theZ, const G4String &aName, G4EvaporationProbability *aEmissionStrategy, G4VCoulombBarrier *aCoulombBarrier)
virtual	~G4EvaporationChannel ()
void	Initialise ()
virtual G4double	GetEmissionProbability (G4Fragment *fragment)
virtual G4Fragment *	EmittedFragment (G4Fragment *theNucleus)
virtual G4FragmentVector *	BreakUp (const G4Fragment &theNucleus)
Public Member Functions inherited from G4VEvaporationChannel
	G4VEvaporationChannel (const G4String &aName="")
virtual	~G4VEvaporationChannel ()
virtual G4double	GetLifeTime (G4Fragment *theNucleus)
virtual G4FragmentVector *	BreakUpFragment (G4Fragment *theNucleus)
virtual G4bool	BreakUpChain (G4FragmentVector *theResult, G4Fragment *theNucleus)
virtual void	Dump () const
virtual void	SetICM (G4bool)
virtual void	RDMForced (G4bool)
virtual G4double	GetFinalLevelEnergy (G4int Z, G4int A, G4double energy)
virtual G4double	GetUpperLevelEnergy (G4int Z, G4int A)
G4double	GetMaxLevelEnergy (G4int Z, G4int A)
G4double	GetNearestLevelEnergy (G4int Z, G4int A, G4double energy)
void	SetPhotonEvaporation (G4VEvaporationChannel *p)
void	SetOPTxs (G4int opt)
void	UseSICB (G4bool use)

Private Member Functions
G4ThreeVector	IsotropicVector (G4double Magnitude=1.0)
	G4EvaporationChannel (const G4EvaporationChannel &right)
const G4EvaporationChannel &	operator= (const G4EvaporationChannel &right)
G4bool	operator== (const G4EvaporationChannel &right) const
G4bool	operator!= (const G4EvaporationChannel &right) const

Private Attributes
G4int	theA
G4int	theZ
G4double	EvaporatedMass
G4double	ResidualMass
G4EvaporationProbability *	theProbability
G4VLevelDensityParameter *	theLevelDensityPtr
G4VCoulombBarrier *	theCoulombBarrier
G4double	CoulombBarrier
G4PairingCorrection *	pairingCorrection
G4int	ResidualA
G4int	ResidualZ
G4double	EmissionProbability
G4double	MaximalKineticEnergy

Additional Inherited Members
Protected Attributes inherited from G4VEvaporationChannel
G4int	OPTxs
G4bool	useSICB

Detailed Description

Definition at line 48 of file G4EvaporationChannel.hh.

Constructor & Destructor Documentation

G4EvaporationChannel() [1/2]

G4EvaporationChannel::G4EvaporationChannel	(	G4int	theA,
		G4int	theZ,
		const G4String &	aName,
		G4EvaporationProbability *	aEmissionStrategy,
		G4VCoulombBarrier *	aCoulombBarrier
	)

Definition at line 53 of file G4EvaporationChannel.cc.

                                                                      :
    G4VEvaporationChannel(aName),
    theA(anA),
    theZ(aZ),
    theProbability(aprob),
    theCoulombBarrier(barrier),
    EmissionProbability(0.0),
    MaximalKineticEnergy(-1000.0)
{ 
  ResidualA = 0;
  ResidualZ = 0;
  ResidualMass = CoulombBarrier = 0.0;
  EvaporatedMass = G4NucleiProperties::GetNuclearMass(theA, theZ);
  theLevelDensityPtr = new G4EvaporationLevelDensityParameter;
  pairingCorrection = G4PairingCorrection::GetInstance();
}

Here is the call graph for this function:

~G4EvaporationChannel()

G4EvaporationChannel::~G4EvaporationChannel ( )

virtual

Definition at line 73 of file G4EvaporationChannel.cc.

{
  delete theLevelDensityPtr;
}

G4EvaporationChannel() [2/2]

G4EvaporationChannel::G4EvaporationChannel ( const G4EvaporationChannel &  right )

private

Member Function Documentation

BreakUp()

G4FragmentVector * G4EvaporationChannel::BreakUp ( const G4Fragment &  theNucleus )

virtual

Implements G4VEvaporationChannel.

Definition at line 160 of file G4EvaporationChannel.cc.

{
  G4FragmentVector * theResult = new G4FragmentVector();
  G4Fragment* frag0 = new G4Fragment(theNucleus);
  G4Fragment* frag1 = EmittedFragment(frag0);
  if(frag1) { theResult->push_back(frag1); }
  theResult->push_back(frag0);
  return theResult;
} 

Here is the call graph for this function:

EmittedFragment()

G4Fragment * G4EvaporationChannel::EmittedFragment ( G4Fragment *  theNucleus )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 138 of file G4EvaporationChannel.cc.

{
  G4Fragment* evFragment = 0;
  G4double evEnergy = EvaporatedMass +
    theProbability->SampleKineticEnergy(CoulombBarrier,
                    MaximalKineticEnergy);

  G4ThreeVector momentum(IsotropicVector
    (std::sqrt((evEnergy - EvaporatedMass)*(evEnergy + EvaporatedMass))));
  
  G4LorentzVector EvaporatedMomentum(momentum, evEnergy);
  G4LorentzVector ResidualMomentum = theNucleus->GetMomentum();
  EvaporatedMomentum.boost(ResidualMomentum.boostVector());
  
  evFragment = new G4Fragment(theA,theZ,EvaporatedMomentum);
  ResidualMomentum -= EvaporatedMomentum;
  theNucleus->SetZandA_asInt(ResidualZ, ResidualA);
  theNucleus->SetMomentum(ResidualMomentum);

  return evFragment; 
} 

Here is the call graph for this function:

Here is the caller graph for this function:

GetEmissionProbability()

G4double G4EvaporationChannel::GetEmissionProbability ( G4Fragment *  fragment )

virtual

Implements G4VEvaporationChannel.

Definition at line 88 of file G4EvaporationChannel.cc.

{
  G4int FragA = fragment->GetA_asInt();
  G4int FragZ = fragment->GetZ_asInt();
  ResidualA = FragA - theA;
  ResidualZ = FragZ - theZ;
  //G4cout << "G4EvaporationChannel::Initialize Z= " << theZ << " A= " << theA 
  //     << " FragZ= " << FragZ << " FragA= " << FragA << G4endl;
  EmissionProbability = 0.0;

  // Only channels which are physically allowed are taken into account 
  if (ResidualA >= ResidualZ && ResidualZ > 0 && ResidualA >= theA) {
  
    //Effective excitation energy
    G4double ExEnergy = fragment->GetExcitationEnergy();
    G4double delta0 = 
      std::max(0.0,pairingCorrection->GetPairingCorrection(FragA,FragZ));
    G4double delta1 = 
      std::max(0.0,pairingCorrection->GetPairingCorrection(ResidualA,ResidualZ));
    ResidualMass = G4NucleiProperties::GetNuclearMass(ResidualA, ResidualZ);
    G4double FragmentMass = fragment->GetGroundStateMass();
    G4double Etot = FragmentMass + ExEnergy;
    G4double ResMass = ResidualMass + delta1;  

    if(ExEnergy >= delta0 && Etot > ResMass + EvaporatedMass) {
  
      // Maximal Kinetic Energy
      MaximalKineticEnergy = ((Etot-ResMass)*(Etot+ResMass) 
        + EvaporatedMass*EvaporatedMass)/(2.0*Etot) - EvaporatedMass;

      // The threshold for charged particle emission must be set to 
      // 0 if Coulomb cutoff  is included in the cross sections
      // Of course for OPTxs=0 we have the Coulomb barrier 

      CoulombBarrier = 0.0;
      if (OPTxs==0 || useSICB) {
    CoulombBarrier = 
      theCoulombBarrier->GetCoulombBarrier(ResidualA,ResidualZ,ExEnergy);
      }
      if (MaximalKineticEnergy > CoulombBarrier) {
    EmissionProbability = theProbability->
      TotalProbability(*fragment, CoulombBarrier, MaximalKineticEnergy);
      }
    }
  }
  //G4cout << "G4EvaporationChannel:: probability= " 
  // << EmissionProbability << G4endl;   
  return EmissionProbability;
}

Here is the call graph for this function:

Initialise()

void G4EvaporationChannel::Initialise ( )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 78 of file G4EvaporationChannel.cc.

{
  //for inverse cross section choice
  theProbability->SetOPTxs(OPTxs);
  // for superimposed Coulomb Barrier for inverse cross sections
  theProbability->UseSICB(useSICB);

  G4VEvaporationChannel::Initialise();  
}

Here is the call graph for this function:

IsotropicVector()

G4ThreeVector G4EvaporationChannel::IsotropicVector ( G4double  Magnitude = 1.0 )

private

Definition at line 170 of file G4EvaporationChannel.cc.

{
  G4double CosTheta = 1.0 - 2.0*G4UniformRand();
  G4double SinTheta = std::sqrt(1.0 - CosTheta*CosTheta);
  G4double Phi = twopi*G4UniformRand();
  G4ThreeVector Vector(Magnitude*std::cos(Phi)*SinTheta,
               Magnitude*std::sin(Phi)*SinTheta,
               Magnitude*CosTheta);
  return Vector;
}

Here is the caller graph for this function:

operator!=()

G4bool G4EvaporationChannel::operator!= ( const G4EvaporationChannel &  right ) const

private

operator=()

const G4EvaporationChannel& G4EvaporationChannel::operator= ( const G4EvaporationChannel &  right )

private

operator==()

G4bool G4EvaporationChannel::operator== ( const G4EvaporationChannel &  right ) const

private

Member Data Documentation

CoulombBarrier

G4double G4EvaporationChannel::CoulombBarrier

private

Definition at line 98 of file G4EvaporationChannel.hh.

EmissionProbability

G4double G4EvaporationChannel::EmissionProbability

private

Definition at line 116 of file G4EvaporationChannel.hh.

EvaporatedMass

G4double G4EvaporationChannel::EvaporatedMass

private

Definition at line 87 of file G4EvaporationChannel.hh.

MaximalKineticEnergy

G4double G4EvaporationChannel::MaximalKineticEnergy

private

Definition at line 119 of file G4EvaporationChannel.hh.

pairingCorrection

G4PairingCorrection* G4EvaporationChannel::pairingCorrection

private

Definition at line 100 of file G4EvaporationChannel.hh.

ResidualA

G4int G4EvaporationChannel::ResidualA

private

Definition at line 110 of file G4EvaporationChannel.hh.

ResidualMass

G4double G4EvaporationChannel::ResidualMass

private

Definition at line 88 of file G4EvaporationChannel.hh.

ResidualZ

G4int G4EvaporationChannel::ResidualZ

private

Definition at line 113 of file G4EvaporationChannel.hh.

theA

G4int G4EvaporationChannel::theA

private

Definition at line 82 of file G4EvaporationChannel.hh.

theCoulombBarrier

G4VCoulombBarrier* G4EvaporationChannel::theCoulombBarrier

private

Definition at line 97 of file G4EvaporationChannel.hh.

theLevelDensityPtr

G4VLevelDensityParameter* G4EvaporationChannel::theLevelDensityPtr

private

Definition at line 94 of file G4EvaporationChannel.hh.

theProbability

G4EvaporationProbability* G4EvaporationChannel::theProbability

private

Definition at line 91 of file G4EvaporationChannel.hh.

theZ

G4int G4EvaporationChannel::theZ

private

Definition at line 85 of file G4EvaporationChannel.hh.

---

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

• G4EvaporationChannel.hh
• G4EvaporationChannel.cc