G4EvaporationChannel Class Reference

#include <G4EvaporationChannel.hh>

Inheritance diagram for G4EvaporationChannel:

Public Member Functions
	G4EvaporationChannel (G4int theA, G4int theZ, const G4String &aName, G4EvaporationProbability *aEmissionStrategy, G4VCoulombBarrier *aCoulombBarrier)
virtual	~G4EvaporationChannel ()
void	SetEmissionStrategy (G4EvaporationProbability *aEmissionStrategy)
void	SetCoulombBarrierStrategy (G4VCoulombBarrier *aCoulombBarrier)
virtual G4double	GetEmissionProbability (G4Fragment *fragment)
G4FragmentVector *	BreakUp (const G4Fragment &theNucleus)
G4double	GetMaximalKineticEnergy (void) const
Public Member Functions inherited from G4VEvaporationChannel
	G4VEvaporationChannel (const G4String &aName="Anonymous", G4EvaporationChannelType timeType=fDelayedEmission)
virtual	~G4VEvaporationChannel ()
virtual G4double	GetLifeTime (G4Fragment *theNucleus)
virtual G4Fragment *	EmittedFragment (G4Fragment *theNucleus)
virtual G4FragmentVector *	BreakUpFragment (G4Fragment *theNucleus)
G4String	GetName () const
void	SetName (const G4String &aName)
void	SetOPTxs (G4int opt)
void	UseSICB (G4bool use)

Protected Member Functions
	G4EvaporationChannel ()

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

Detailed Description

Definition at line 47 of file G4EvaporationChannel.hh.

Constructor & Destructor Documentation

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

Definition at line 53 of file G4EvaporationChannel.cc.

References G4NucleiProperties::GetNuclearMass().

                                                                              :
    G4VEvaporationChannel(aName),
    theA(anA),
    theZ(aZ),
    theEvaporationProbabilityPtr(aEmissionStrategy),
    theCoulombBarrierPtr(aCoulombBarrier),
    EmissionProbability(0.0),
    MaximalKineticEnergy(-1000.0)
{ 
  ResidualA = 0;
  ResidualZ = 0;
  ResidualMass = CoulombBarrier=0.0;
  EvaporatedMass = G4NucleiProperties::GetNuclearMass(theA, theZ);
  theLevelDensityPtr = new G4EvaporationLevelDensityParameter;
}

G4EvaporationChannel::~G4EvaporationChannel ( )

virtual

Definition at line 87 of file G4EvaporationChannel.cc.

{
  delete theLevelDensityPtr;
}

G4EvaporationChannel::G4EvaporationChannel ( )

protected

Definition at line 72 of file G4EvaporationChannel.cc.

                                          :
    G4VEvaporationChannel(""),
    theA(0),
    theZ(0),
    theEvaporationProbabilityPtr(0),
    theCoulombBarrierPtr(0),
    EmissionProbability(0.0),
    MaximalKineticEnergy(-1000.0)
{ 
  ResidualA = 0;
  ResidualZ = 0;
  EvaporatedMass = ResidualMass = CoulombBarrier = 0.0;
  theLevelDensityPtr = new G4EvaporationLevelDensityParameter;
}

Member Function Documentation

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

virtual

Implements G4VEvaporationChannel.

Definition at line 154 of file G4EvaporationChannel.cc.

References CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), CLHEP::HepLorentzVector::e(), G4cout, G4endl, G4Fragment::GetMomentum(), python.hepunit::keV, python.hepunit::MeV, CLHEP::HepLorentzVector::vect(), CLHEP::Hep3Vector::x(), CLHEP::HepLorentzVector::x(), CLHEP::Hep3Vector::y(), CLHEP::HepLorentzVector::y(), CLHEP::Hep3Vector::z(), and CLHEP::HepLorentzVector::z().

{
  /*
  G4double Ecm = GetKineticEnergy(theNucleus) + ResidualMass + EvaporatedMass;
  
  G4double EvaporatedEnergy = 
    ((Ecm-ResidualMass)*(Ecm+ResidualMass) + EvaporatedMass*EvaporatedMass)/(2*Ecm);
  */
  G4double EvaporatedEnergy = GetKineticEnergy(theNucleus) + EvaporatedMass;

  G4ThreeVector momentum(IsotropicVector
                         (std::sqrt((EvaporatedEnergy - EvaporatedMass)*
                                    (EvaporatedEnergy + EvaporatedMass))));
  
  G4LorentzVector EvaporatedMomentum(momentum,EvaporatedEnergy);
  G4LorentzVector ResidualMomentum = theNucleus.GetMomentum();
  EvaporatedMomentum.boost(ResidualMomentum.boostVector());
  
  G4Fragment * EvaporatedFragment = new G4Fragment(theA,theZ,EvaporatedMomentum);
  ResidualMomentum -= EvaporatedMomentum;

  G4Fragment * ResidualFragment = new G4Fragment(ResidualA, ResidualZ, ResidualMomentum);
 
  G4FragmentVector * theResult = new G4FragmentVector;
  
#ifdef debug
  G4double Efinal = ResidualMomentum.e() + EvaporatedMomentum.e();
  G4ThreeVector Pfinal = ResidualMomentum.vect() + EvaporatedMomentum.vect();
  if (std::abs(Efinal-theNucleus.GetMomentum().e()) > 1.0*keV) {
    G4cout << "@@@@@@@@@@@@@@@@@@@@@ G4Evaporation Chanel: ENERGY @@@@@@@@@@@@@@@@" 
       << G4endl;
    G4cout << "Initial : " << theNucleus.GetMomentum().e()/MeV << " MeV    Final :" 
           <<Efinal/MeV << " MeV    Delta: " 
       <<  (Efinal-theNucleus.GetMomentum().e())/MeV 
           << " MeV" << G4endl;
  }
  if (std::abs(Pfinal.x()-theNucleus.GetMomentum().x()) > 1.0*keV ||
      std::abs(Pfinal.y()-theNucleus.GetMomentum().y()) > 1.0*keV ||
      std::abs(Pfinal.z()-theNucleus.GetMomentum().z()) > 1.0*keV ) {
    G4cout << "@@@@@@@@@@@@@@@@@@@@@ G4Evaporation Chanel: MOMENTUM @@@@@@@@@@@@@@@@" 
       << G4endl;
    G4cout << "Initial : " << theNucleus.GetMomentum().vect() << " MeV    Final :" 
           <<Pfinal/MeV << " MeV    Delta: " <<  Pfinal-theNucleus.GetMomentum().vect()
           << " MeV" << G4endl;   
    
  }
#endif
  theResult->push_back(EvaporatedFragment);
  theResult->push_back(ResidualFragment);
  return theResult; 
} 

G4double G4EvaporationChannel::GetEmissionProbability ( G4Fragment *  fragment )

virtual

Implements G4VEvaporationChannel.

Definition at line 95 of file G4EvaporationChannel.cc.

References G4Fragment::GetA_asInt(), G4VCoulombBarrier::GetCoulombBarrier(), G4Fragment::GetExcitationEnergy(), G4Fragment::GetGroundStateMass(), G4PairingCorrection::GetInstance(), G4NucleiProperties::GetNuclearMass(), G4PairingCorrection::GetPairingCorrection(), G4Fragment::GetZ_asInt(), G4INCL::Math::max(), python.hepunit::MeV, G4VEvaporationChannel::OPTxs, G4VEmissionProbability::SetOPTxs(), G4VEmissionProbability::UseSICB(), and G4VEvaporationChannel::useSICB.

{
  //for inverse cross section choice
  theEvaporationProbabilityPtr->SetOPTxs(OPTxs);
  // for superimposed Coulomb Barrier for inverse cross sections
  theEvaporationProbabilityPtr->UseSICB(useSICB);
  
  G4int FragmentA = fragment->GetA_asInt();
  G4int FragmentZ = fragment->GetZ_asInt();
  ResidualA = FragmentA - theA;
  ResidualZ = FragmentZ - theZ;
  //G4cout << "G4EvaporationChannel::Initialize Z= " << theZ << " A= " << theA 
  //     << " FragZ= " << FragmentZ << " FragA= " << FragmentA << G4endl;
  
  //Effective excitation energy
  G4double ExEnergy = fragment->GetExcitationEnergy() - 
    G4PairingCorrection::GetInstance()->GetPairingCorrection(FragmentA,FragmentZ);
  
  // Only channels which are physically allowed are taken into account 
  if (ResidualA <= 0 || ResidualZ <= 0 || ResidualA < ResidualZ ||
      (ResidualA == ResidualZ && ResidualA > 1) || ExEnergy <= 0.0) {
    CoulombBarrier = ResidualMass = 0.0;
    MaximalKineticEnergy = -1000.0*MeV;
    EmissionProbability = 0.0;
  } else {
    ResidualMass = G4NucleiProperties::GetNuclearMass(ResidualA, ResidualZ);
    G4double FragmentMass = fragment->GetGroundStateMass();
    CoulombBarrier = 
      theCoulombBarrierPtr->GetCoulombBarrier(ResidualA,ResidualZ,ExEnergy);
    // Maximal Kinetic Energy
    MaximalKineticEnergy = CalcMaximalKineticEnergy(FragmentMass + ExEnergy);
    //MaximalKineticEnergy = ExEnergy + fragment.GetGroundStateMass() 
    //  - EvaporatedMass - ResidualMass;
    
    // Emission probability
    // Protection for the case Tmax<V. If not set in this way we could end up in an 
    // infinite loop in  the method GetKineticEnergy if OPTxs!=0 && useSICB=true. 
    // Of course for OPTxs=0 we have the Coulomb barrier 
    
    G4double limit;
    if (OPTxs==0 || (OPTxs!=0 && useSICB)) 
      limit= CoulombBarrier;
    else limit=0.;
    limit = std::max(limit, FragmentMass - ResidualMass - EvaporatedMass);
  
    // The threshold for charged particle emission must be  set to 0 if Coulomb 
    //cutoff  is included in the cross sections
    if (MaximalKineticEnergy <= limit) EmissionProbability = 0.0;
    else { 
      // Total emission probability for this channel
      EmissionProbability = theEvaporationProbabilityPtr->
        EmissionProbability(*fragment, MaximalKineticEnergy);
    }
  }
  //G4cout << "G4EvaporationChannel:: probability= " << EmissionProbability << G4endl; 
  
  return EmissionProbability;
}

G4double G4EvaporationChannel::GetMaximalKineticEnergy ( void  ) const

inline

Definition at line 76 of file G4EvaporationChannel.hh.

  { return MaximalKineticEnergy; }

void G4EvaporationChannel::SetCoulombBarrierStrategy ( G4VCoulombBarrier *  aCoulombBarrier )

inline

Definition at line 61 of file G4EvaporationChannel.hh.

  {theCoulombBarrierPtr = aCoulombBarrier;} 

void G4EvaporationChannel::SetEmissionStrategy ( G4EvaporationProbability *  aEmissionStrategy )

inline

Definition at line 58 of file G4EvaporationChannel.hh.

  {theEvaporationProbabilityPtr = aEmissionStrategy;}

---

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

• G4EvaporationChannel.hh
• G4EvaporationChannel.cc