G4GEMChannel Class Reference

#include <G4GEMChannel.hh>

Inheritance diagram for G4GEMChannel:

Public Member Functions
	G4GEMChannel (const G4int theA, const G4int theZ, const G4String &aName, G4GEMProbability *aEmissionStrategy, G4VCoulombBarrier *aCoulombBarrier)
virtual	~G4GEMChannel ()
virtual G4double	GetEmissionProbability (G4Fragment *theNucleus)
virtual G4FragmentVector *	BreakUp (const G4Fragment &theNucleus)
void	SetLevelDensityParameter (G4VLevelDensityParameter *aLevelDensity)
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)

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

Detailed Description

Definition at line 49 of file G4GEMChannel.hh.

Constructor & Destructor Documentation

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

Definition at line 46 of file G4GEMChannel.cc.

References G4Pow::GetInstance(), and G4NucleiProperties::GetNuclearMass().

                                                                :
  G4VEvaporationChannel(aName),
  A(theA),
  Z(theZ),
  theEvaporationProbabilityPtr(aEmissionStrategy),
  theCoulombBarrierPtr(aCoulombBarrier),
  EmissionProbability(0.0),
  MaximalKineticEnergy(-CLHEP::GeV)
{ 
  theLevelDensityPtr = new G4EvaporationLevelDensityParameter;
  MyOwnLevelDensity = true;
  EvaporatedMass = G4NucleiProperties::GetNuclearMass(A, Z);
  ResidualMass = CoulombBarrier = 0.0;
  fG4pow = G4Pow::GetInstance(); 
  ResidualZ = ResidualA = 0;
}

G4GEMChannel::~G4GEMChannel ( )

virtual

Definition at line 65 of file G4GEMChannel.cc.

{
  if (MyOwnLevelDensity) { delete theLevelDensityPtr; }
}

Member Function Documentation

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

virtual

Implements G4VEvaporationChannel.

Definition at line 137 of file G4GEMChannel.cc.

References CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), G4Fragment::GetExcitationEnergy(), G4Fragment::GetGroundStateMass(), G4Fragment::GetMomentum(), CLHEP::Hep3Vector::rotateUz(), CLHEP::Hep3Vector::unit(), and CLHEP::HepLorentzVector::vect().

{
  G4double EvaporatedKineticEnergy = CalcKineticEnergy(theNucleus);
  G4double EvaporatedEnergy = EvaporatedKineticEnergy + EvaporatedMass;
  
  G4ThreeVector momentum(IsotropicVector(std::sqrt(EvaporatedKineticEnergy*
                           (EvaporatedKineticEnergy+2.0*EvaporatedMass))));
    
  momentum.rotateUz(theNucleus.GetMomentum().vect().unit());

  G4LorentzVector EvaporatedMomentum(momentum,EvaporatedEnergy);
  EvaporatedMomentum.boost(theNucleus.GetMomentum().boostVector());
  G4Fragment * EvaporatedFragment = new G4Fragment(A,Z,EvaporatedMomentum);
  // ** And now the residual nucleus ** 
  G4double theExEnergy = theNucleus.GetExcitationEnergy();
  G4double theMass = theNucleus.GetGroundStateMass();
  G4double ResidualEnergy = 
    theMass + (theExEnergy - EvaporatedKineticEnergy) - EvaporatedMass;
    
  G4LorentzVector ResidualMomentum(-momentum,ResidualEnergy);
  ResidualMomentum.boost(theNucleus.GetMomentum().boostVector());
    
  G4Fragment * ResidualFragment = new G4Fragment( ResidualA, ResidualZ, ResidualMomentum );
    
  G4FragmentVector * theResult = new G4FragmentVector;
    
  theResult->push_back(EvaporatedFragment);
  theResult->push_back(ResidualFragment);
  return theResult; 
} 

G4double G4GEMChannel::GetEmissionProbability ( G4Fragment *  theNucleus )

virtual

Implements G4VEvaporationChannel.

Definition at line 70 of file G4GEMChannel.cc.

References G4GEMProbability::EmissionProbability(), G4Fragment::GetA_asInt(), G4VCoulombBarrier::GetCoulombBarrier(), G4Fragment::GetExcitationEnergy(), G4Fragment::GetGroundStateMass(), G4PairingCorrection::GetInstance(), G4NucleiProperties::GetNuclearMass(), G4PairingCorrection::GetPairingCorrection(), G4Fragment::GetZ_asInt(), and python.hepunit::MeV.

{
  G4int anA = fragment->GetA_asInt();
  G4int aZ  = fragment->GetZ_asInt();
  ResidualA = anA - A;
  ResidualZ = aZ - Z;
  //G4cout << "G4GEMChannel::Initialize: Z= " << aZ << " A= " << anA
  //       << " Zres= " << ResidualZ << " Ares= " << ResidualA << G4endl; 

  // We only take into account channels which are physically allowed
  if (ResidualA <= 0 || ResidualZ <= 0 || ResidualA < ResidualZ ||
      (ResidualA == ResidualZ && ResidualA > 1)) 
    {
      CoulombBarrier = 0.0;
      MaximalKineticEnergy = -CLHEP::GeV;
      EmissionProbability = 0.0;
    } 
  else 
    {
      // Effective excitation energy
      // JMQ 071009: pairing in ExEnergy should be the one of parent compound nucleus 
      // FIXED the bug causing reported crash by VI (negative Probabilities 
      // due to inconsistency in Coulomb barrier calculation (CoulombBarrier and -Beta 
      // param for protons must be the same)   
      //    G4double ExEnergy = fragment.GetExcitationEnergy() -
      //    G4PairingCorrection::GetInstance()->GetPairingCorrection(ResidualA,ResidualZ);
      G4double ExEnergy = fragment->GetExcitationEnergy() -
    G4PairingCorrection::GetInstance()->GetPairingCorrection(anA,aZ);

      //G4cout << "Eexc(MeV)= " << ExEnergy/MeV << G4endl;

      if( ExEnergy <= 0.0) {
    CoulombBarrier = 0.0;
    MaximalKineticEnergy = -1000.0*MeV;
    EmissionProbability = 0.0;

      } else {

    ResidualMass = G4NucleiProperties::GetNuclearMass(ResidualA, ResidualZ);

        // Coulomb Barrier calculation
    CoulombBarrier = theCoulombBarrierPtr->GetCoulombBarrier(ResidualA,ResidualZ,ExEnergy);
    //G4cout << "CBarrier(MeV)= " << CoulombBarrier/MeV << G4endl;

    //Maximal kinetic energy (JMQ : at the Coulomb barrier)
    MaximalKineticEnergy = 
      CalcMaximalKineticEnergy(fragment->GetGroundStateMass()+ExEnergy);
    //G4cout << "MaxE(MeV)= " << MaximalKineticEnergy/MeV << G4endl;
        
    // Emission probability
    if (MaximalKineticEnergy <= 0.0) 
      {
        EmissionProbability = 0.0;
      }
    else 
      { 
        // Total emission probability for this channel
        EmissionProbability = 
          theEvaporationProbabilityPtr->EmissionProbability(*fragment,
                                MaximalKineticEnergy);
      }
      }
    }   
  //G4cout << "Prob= " << EmissionProbability << G4endl;
  return EmissionProbability;
}

G4double G4GEMChannel::GetMaximalKineticEnergy ( void  ) const

inline

Definition at line 71 of file G4GEMChannel.hh.

  { return MaximalKineticEnergy; }

void G4GEMChannel::SetLevelDensityParameter ( G4VLevelDensityParameter *  aLevelDensity )

inline

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