G4GEMChannel Class Reference

#include <G4GEMChannel.hh>

Inheritance diagram for G4GEMChannel:

Public Member Functions
	G4GEMChannel (G4int theA, G4int theZ, const G4String &aName, G4GEMProbability *aEmissionStrategy)
virtual	~G4GEMChannel ()
virtual G4double	GetEmissionProbability (G4Fragment *theNucleus)
virtual G4Fragment *	EmittedFragment (G4Fragment *theNucleus)
virtual void	Dump () const
void	SetLevelDensityParameter (G4VLevelDensityParameter *aLevelDensity)
Public Member Functions inherited from G4VEvaporationChannel
	G4VEvaporationChannel (const G4String &aName="")
virtual	~G4VEvaporationChannel ()=default
virtual G4double	GetEmissionProbability (G4Fragment *theNucleus)=0
virtual void	Initialise ()
virtual G4double	GetLifeTime (G4Fragment *theNucleus)
virtual G4Fragment *	EmittedFragment (G4Fragment *theNucleus)
virtual G4bool	BreakUpChain (G4FragmentVector *theResult, G4Fragment *theNucleus)
G4FragmentVector *	BreakUpFragment (G4Fragment *theNucleus)
virtual G4double	ComputeInverseXSection (G4Fragment *theNucleus, G4double kinEnergy)
virtual G4double	ComputeProbability (G4Fragment *theNucleus, G4double kinEnergy)
virtual void	Dump () const
virtual void	SetICM (G4bool)
virtual void	RDMForced (G4bool)
void	SetOPTxs (G4int opt)
void	UseSICB (G4bool use)
	G4VEvaporationChannel (const G4VEvaporationChannel &right)=delete
const G4VEvaporationChannel &	operator= (const G4VEvaporationChannel &right)=delete
G4bool	operator== (const G4VEvaporationChannel &right) const =delete
G4bool	operator!= (const G4VEvaporationChannel &right) const =delete

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

Detailed Description

Definition at line 47 of file G4GEMChannel.hh.

Constructor & Destructor Documentation

G4GEMChannel()

G4GEMChannel::G4GEMChannel ( G4int  theA, G4int  theZ, const G4String &  aName, G4GEMProbability *  aEmissionStrategy  )

explicit

Definition at line 49 of file G4GEMChannel.cc.

                                                                 :
  G4VEvaporationChannel(aName),
  A(theA),
  Z(theZ),
  EmissionProbability(0.0),
  MaximalKineticEnergy(-CLHEP::GeV),
  theEvaporationProbabilityPtr(aEmissionStrategy),
  secID(-1)
{ 
  theCoulombBarrierPtr = new G4GEMCoulombBarrier(theA, theZ);
  theEvaporationProbabilityPtr->SetCoulomBarrier(theCoulombBarrierPtr);
  theLevelDensityPtr = new G4EvaporationLevelDensityParameter;
  MyOwnLevelDensity = true;
  EvaporatedMass = G4NucleiProperties::GetNuclearMass(A, Z);
  ResidualMass = CoulombBarrier = 0.0;
  fG4pow = G4Pow::GetInstance(); 
  ResidualZ = ResidualA = 0;
  fNucData = G4NuclearLevelData::GetInstance();
  secID = G4PhysicsModelCatalog::GetModelID("model_G4GEMChannel");
}

~G4GEMChannel()

G4GEMChannel::~G4GEMChannel ( )

virtual

Definition at line 71 of file G4GEMChannel.cc.

{
  if (MyOwnLevelDensity) { delete theLevelDensityPtr; }
  delete theCoulombBarrierPtr;
}

Member Function Documentation

Dump()

void G4GEMChannel::Dump ( ) const

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 254 of file G4GEMChannel.cc.

{
  theEvaporationProbabilityPtr->Dump();
}

EmittedFragment()

G4Fragment * G4GEMChannel::EmittedFragment ( G4Fragment *  theNucleus )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 130 of file G4GEMChannel.cc.

{
  G4Fragment* evFragment = 0;
  G4double evEnergy = SampleKineticEnergy(*theNucleus) + EvaporatedMass;
 
  G4ThreeVector momentum = G4RandomDirection()*
    std::sqrt((evEnergy - EvaporatedMass)*(evEnergy + EvaporatedMass));
  
  G4LorentzVector EvaporatedMomentum(momentum, evEnergy);
  G4LorentzVector ResidualMomentum = theNucleus->GetMomentum();
  EvaporatedMomentum.boost(ResidualMomentum.boostVector());
  
  evFragment = new G4Fragment(A, Z, EvaporatedMomentum);
  if ( evFragment != nullptr ) { evFragment->SetCreatorModelID(secID); }
  ResidualMomentum -= EvaporatedMomentum;
  theNucleus->SetZandA_asInt(ResidualZ, ResidualA);
  theNucleus->SetMomentum(ResidualMomentum);
  theNucleus->SetCreatorModelID(secID);
  
  return evFragment; 
} 

GetEmissionProbability()

G4double G4GEMChannel::GetEmissionProbability ( G4Fragment *  theNucleus )

virtual

Implements G4VEvaporationChannel.

Definition at line 77 of file G4GEMChannel.cc.

{
  G4int anA = fragment->GetA_asInt();
  G4int aZ  = fragment->GetZ_asInt();
  ResidualA = anA - A;
  ResidualZ = aZ - Z;
  /*
  G4cout << "G4GEMChannel: Z= " << Z << "  A= " << A 
     << " FragmentZ= " << aZ << " FragmentA= " << anA
     << " Zres= " << ResidualZ << " Ares= " << ResidualA 
     << G4endl; 
  */
  // We only take into account channels which are physically allowed
  EmissionProbability = 0.0;
 
  // Only channels which are physically allowed are taken into account 
  if (ResidualA >= ResidualZ && ResidualZ >= 0 && ResidualA >= A) {
  
    //Effective excitation energy
    G4double ExEnergy = fragment->GetExcitationEnergy()
      - fNucData->GetPairingCorrection(aZ, anA);
    if(ExEnergy > 0.0) { 
      ResidualMass = G4NucleiProperties::GetNuclearMass(ResidualA, ResidualZ);
      G4double FragmentMass = fragment->GetGroundStateMass();
      G4double Etot = FragmentMass + ExEnergy;
      // Coulomb Barrier calculation
      CoulombBarrier = 
    theCoulombBarrierPtr->GetCoulombBarrier(ResidualA,ResidualZ,ExEnergy);
      /*  
      G4cout << "Eexc(MeV)= " << ExEnergy/MeV 
         << " CoulBarrier(MeV)= " << CoulombBarrier/MeV << G4endl;
      */
      if(Etot > ResidualMass + EvaporatedMass + CoulombBarrier) {
  
    // Maximal Kinetic Energy
    MaximalKineticEnergy = ((Etot-ResidualMass)*(Etot+ResidualMass) 
      + EvaporatedMass*EvaporatedMass)/(2.0*Etot) 
      - EvaporatedMass - CoulombBarrier;
 
    //G4cout << "CBarrier(MeV)= " << CoulombBarrier/MeV << G4endl;
 
    if (MaximalKineticEnergy > 0.0) { 
      // Total emission probability for this channel
      EmissionProbability = theEvaporationProbabilityPtr->
        EmissionProbability(*fragment, MaximalKineticEnergy);
    }
      }
    }
  }   
  //G4cout << "Prob= " << EmissionProbability << G4endl;
  return EmissionProbability;
}

SetLevelDensityParameter()

void G4GEMChannel::SetLevelDensityParameter ( G4VLevelDensityParameter *  aLevelDensity )

inline

Definition at line 63 of file G4GEMChannel.hh.

  {
    if (MyOwnLevelDensity) { delete theLevelDensityPtr; }
    theLevelDensityPtr = aLevelDensity;
    MyOwnLevelDensity = false;
  }

---

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

• G4GEMChannel.hh
• G4GEMChannel.cc