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")
virtual	~G4VEvaporationChannel ()
virtual G4Fragment *	EmittedFragment (G4Fragment *theNucleus)
virtual G4FragmentVector *	BreakUpFragment (G4Fragment *theNucleus)
virtual G4FragmentVector *	BreakUp (const G4Fragment &theNucleus)=0
virtual G4double	GetEmissionProbability (G4Fragment *theNucleus)=0
G4String	GetName () const
void	SetName (const G4String &aName)
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 49 of file G4GEMChannel.hh.

Constructor & Destructor Documentation

G4GEMChannel()

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

Definition at line 44 of file G4GEMChannel.cc.

                                                                :
  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::~G4GEMChannel ( )

virtual

Definition at line 63 of file G4GEMChannel.cc.

{
  if (MyOwnLevelDensity) { delete theLevelDensityPtr; }
}

Member Function Documentation

BreakUp()

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

virtual

Implements G4VEvaporationChannel.

Definition at line 135 of file G4GEMChannel.cc.

{
  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; 
} 

GetEmissionProbability()

G4double G4GEMChannel::GetEmissionProbability ( G4Fragment *  theNucleus )

virtual

Implements G4VEvaporationChannel.

Definition at line 68 of file G4GEMChannel.cc.

{
  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;
}

GetMaximalKineticEnergy()

G4double G4GEMChannel::GetMaximalKineticEnergy ( void  ) const

inline

Definition at line 71 of file G4GEMChannel.hh.

  { return MaximalKineticEnergy; }

SetLevelDensityParameter()

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