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 ()
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 void	Dump () const
virtual void	SetICM (G4bool)
virtual void	RDMForced (G4bool)
void	SetOPTxs (G4int opt)
void	UseSICB (G4bool use)

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 48 of file G4GEMChannel.cc.

                                                                 :
  G4VEvaporationChannel(aName),
  A(theA),
  Z(theZ),
  EmissionProbability(0.0),
  MaximalKineticEnergy(-CLHEP::GeV),
  theEvaporationProbabilityPtr(aEmissionStrategy)
{ 
  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();
}

~G4GEMChannel()

G4GEMChannel::~G4GEMChannel ( )

virtual

Definition at line 68 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 249 of file G4GEMChannel.cc.

{
  theEvaporationProbabilityPtr->Dump();
}

EmittedFragment()

G4Fragment * G4GEMChannel::EmittedFragment ( G4Fragment *  theNucleus )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 127 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);
  ResidualMomentum -= EvaporatedMomentum;
  theNucleus->SetZandA_asInt(ResidualZ, ResidualA);
  theNucleus->SetMomentum(ResidualMomentum);
 
  return evFragment; 
} 

GetEmissionProbability()

G4double G4GEMChannel::GetEmissionProbability ( G4Fragment *  theNucleus )

virtual

Implements G4VEvaporationChannel.

Definition at line 74 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