G4GEMChannel Class Reference

#include <G4GEMChannel.hh>

Inherits G4VEvaporationChannel.

Inherited by G4AlphaGEMChannel, G4B10GEMChannel, G4B11GEMChannel, G4B12GEMChannel, G4B13GEMChannel, G4B8GEMChannel, G4Be10GEMChannel, G4Be11GEMChannel, G4Be12GEMChannel, G4Be7GEMChannel, G4Be9GEMChannel, G4C10GEMChannel, G4C11GEMChannel, G4C12GEMChannel, G4C13GEMChannel, G4C14GEMChannel, G4C15GEMChannel, G4C16GEMChannel, G4DeuteronGEMChannel, G4F17GEMChannel, G4F18GEMChannel, G4F19GEMChannel, G4F20GEMChannel, G4F21GEMChannel, G4He3GEMChannel, G4He6GEMChannel, G4He8GEMChannel, G4Li6GEMChannel, G4Li7GEMChannel, G4Li8GEMChannel, G4Li9GEMChannel, G4Mg22GEMChannel, G4Mg23GEMChannel, G4Mg24GEMChannel, G4Mg25GEMChannel, G4Mg26GEMChannel, G4Mg27GEMChannel, G4Mg28GEMChannel, G4N12GEMChannel, G4N13GEMChannel, G4N14GEMChannel, G4N15GEMChannel, G4N16GEMChannel, G4N17GEMChannel, G4Na21GEMChannel, G4Na22GEMChannel, G4Na23GEMChannel, G4Na24GEMChannel, G4Na25GEMChannel, G4Ne18GEMChannel, G4Ne19GEMChannel, G4Ne20GEMChannel, G4Ne21GEMChannel, G4Ne22GEMChannel, G4Ne23GEMChannel, G4Ne24GEMChannel, G4NeutronGEMChannel, G4O14GEMChannel, G4O15GEMChannel, G4O16GEMChannel, G4O17GEMChannel, G4O18GEMChannel, G4O19GEMChannel, G4O20GEMChannel, G4ProtonGEMChannel, and G4TritonGEMChannel.

Collaboration 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 G4Fragment *	EmittedFragment (G4Fragment *theNucleus)
virtual G4FragmentVector *	BreakUp (const G4Fragment &theNucleus)
virtual void	Dump () const
void	SetLevelDensityParameter (G4VLevelDensityParameter *aLevelDensity)
Public Member Functions inherited from G4VEvaporationChannel
	G4VEvaporationChannel (const G4String &aName="")
virtual	~G4VEvaporationChannel ()
virtual void	Initialise ()
virtual G4double	GetLifeTime (G4Fragment *theNucleus)
virtual G4FragmentVector *	BreakUpFragment (G4Fragment *theNucleus)
virtual G4bool	BreakUpChain (G4FragmentVector *theResult, G4Fragment *theNucleus)
virtual void	SetICM (G4bool)
virtual void	RDMForced (G4bool)
virtual G4double	GetFinalLevelEnergy (G4int Z, G4int A, G4double energy)
virtual G4double	GetUpperLevelEnergy (G4int Z, G4int A)
G4double	GetMaxLevelEnergy (G4int Z, G4int A)
G4double	GetNearestLevelEnergy (G4int Z, G4int A, G4double energy)
void	SetPhotonEvaporation (G4VEvaporationChannel *p)
void	SetOPTxs (G4int opt)
void	UseSICB (G4bool use)

Private Member Functions
G4double	SampleKineticEnergy (const G4Fragment &fragment)
G4ThreeVector	IsotropicVector (G4double Magnitude=1.0)
	G4GEMChannel (const G4GEMChannel &right)
const G4GEMChannel &	operator= (const G4GEMChannel &right)
G4bool	operator== (const G4GEMChannel &right) const
G4bool	operator!= (const G4GEMChannel &right) const

Private Attributes
G4int	A
G4int	Z
G4double	EvaporatedMass
G4double	ResidualMass
G4Pow *	fG4pow
G4GEMProbability *	theEvaporationProbabilityPtr
G4bool	MyOwnLevelDensity
G4VLevelDensityParameter *	theLevelDensityPtr
G4VCoulombBarrier *	theCoulombBarrierPtr
G4double	CoulombBarrier
G4PairingCorrection *	pairingCorrection
G4int	ResidualA
G4int	ResidualZ
G4double	EmissionProbability
G4double	MaximalKineticEnergy

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() [1/2]

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

Definition at line 46 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;
  pairingCorrection = G4PairingCorrection::GetInstance();
}

Here is the call graph for this function:

Here is the caller graph for this function:

~G4GEMChannel()

G4GEMChannel::~G4GEMChannel ( )

virtual

Definition at line 66 of file G4GEMChannel.cc.

{
  if (MyOwnLevelDensity) { delete theLevelDensityPtr; }
}

G4GEMChannel() [2/2]

G4GEMChannel::G4GEMChannel ( const G4GEMChannel &  right )

private

Member Function Documentation

BreakUp()

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

virtual

Implements G4VEvaporationChannel.

Definition at line 144 of file G4GEMChannel.cc.

{
  G4FragmentVector * theResult = new G4FragmentVector();
  G4Fragment* frag0 = new G4Fragment(theNucleus);
  G4Fragment* frag1 = EmittedFragment(frag0);
  if(frag1) { theResult->push_back(frag1); }
  theResult->push_back(frag0);
  return theResult;
} 

Here is the call graph for this function:

Dump()

void G4GEMChannel::Dump ( ) const

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 270 of file G4GEMChannel.cc.

{
  theEvaporationProbabilityPtr->Dump();
}

Here is the call graph for this function:

EmittedFragment()

G4Fragment * G4GEMChannel::EmittedFragment ( G4Fragment *  theNucleus )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 124 of file G4GEMChannel.cc.

{
  G4Fragment* evFragment = 0;
  G4double evEnergy = SampleKineticEnergy(*theNucleus) + EvaporatedMass;

  G4ThreeVector momentum(IsotropicVector
    (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; 
} 

Here is the call graph for this function:

Here is the caller graph for this function:

GetEmissionProbability()

G4double G4GEMChannel::GetEmissionProbability ( G4Fragment *  theNucleus )

virtual

Implements G4VEvaporationChannel.

Definition at line 71 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()
      - pairingCorrection->GetPairingCorrection(anA, aZ);
    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;
}

Here is the call graph for this function:

IsotropicVector()

G4ThreeVector G4GEMChannel::IsotropicVector ( G4double  Magnitude = 1.0 )

private

Definition at line 257 of file G4GEMChannel.cc.

{
  G4double CosTheta = 1.0 - 2.0*G4UniformRand();
  G4double SinTheta = std::sqrt(1.0 - CosTheta*CosTheta);
  G4double Phi = twopi*G4UniformRand();
  G4ThreeVector Vector(Magnitude*std::cos(Phi)*SinTheta,
               Magnitude*std::sin(Phi)*SinTheta,
               Magnitude*CosTheta);
  return Vector;
}

Here is the caller graph for this function:

operator!=()

G4bool G4GEMChannel::operator!= ( const G4GEMChannel &  right ) const

private

Here is the caller graph for this function:

operator=()

const G4GEMChannel& G4GEMChannel::operator= ( const G4GEMChannel &  right )

private

Here is the caller graph for this function:

operator==()

G4bool G4GEMChannel::operator== ( const G4GEMChannel &  right ) const

private

Here is the caller graph for this function:

SampleKineticEnergy()

G4double G4GEMChannel::SampleKineticEnergy ( const G4Fragment &  fragment )

private

Definition at line 154 of file G4GEMChannel.cc.

{
  G4double U = fragment.GetExcitationEnergy();
  
  G4double Alpha = theEvaporationProbabilityPtr->CalcAlphaParam(fragment);
  G4double Beta = theEvaporationProbabilityPtr->CalcBetaParam(fragment);

  //                       ***RESIDUAL***
  //JMQ (September 2009) the following quantities  refer to the RESIDUAL:
  G4double delta0 = pairingCorrection->GetPairingCorrection(ResidualA,ResidualZ);
  G4double Ux = (2.5 + 150.0/ResidualA)*MeV;
  G4double Ex = Ux + delta0;
  G4double InitialLevelDensity;
  //                    ***end RESIDUAL ***
  
  //                       ***PARENT***
  //JMQ (September 2009) the following quantities   refer to the PARENT:
  
  G4double deltaCN = pairingCorrection->GetPairingCorrection(fragment.GetA_asInt(),
                                 fragment.GetZ_asInt());
  G4double aCN = theLevelDensityPtr->LevelDensityParameter(fragment.GetA_asInt(),
                               fragment.GetZ_asInt(),
                               U-deltaCN);   
  G4double UxCN = (2.5 + 150.0/G4double(fragment.GetA_asInt()))*MeV;
  G4double ExCN = UxCN + deltaCN;
  G4double TCN = 1.0/(std::sqrt(aCN/UxCN) - 1.5/UxCN);
  //                       ***end PARENT***
  
  //JMQ quantities calculated for  CN in InitialLevelDensity
  if ( U < ExCN ) 
    {
      G4double E0CN = ExCN - TCN*(G4Log(TCN/MeV) - G4Log(aCN*MeV)/4.0 
                  - 1.25*G4Log(UxCN/MeV) + 2.0*std::sqrt(aCN*UxCN));
      InitialLevelDensity = (pi/12.0)*G4Exp((U-E0CN)/TCN)/TCN;
    } 
  else 
    {
      G4double x  = U-deltaCN;
      G4double x1 = std::sqrt(aCN*x);
      InitialLevelDensity = (pi/12.0)*G4Exp(2*x1)/(x*std::sqrt(x1));
    }
  
  G4double Spin = theEvaporationProbabilityPtr->GetSpin();
  //JMQ  BIG BUG fixed: hbarc instead of hbar_Planck !!!!
  //     it was fixed in total probability (for this channel) but remained still here!!
  //    G4double g = (2.0*Spin+1.0)*NuclearMass/(pi2* hbar_Planck*hbar_Planck);
  G4double gg = (2.0*Spin+1.0)*EvaporatedMass/(pi2* hbarc*hbarc);
  //
  //JMQ  fix on Rb and  geometrical cross sections according to Furihata's paper 
  //                      (JAERI-Data/Code 2001-105, p6)
  G4double Rb = 0.0; 
  if (A > 4) 
    {
      G4double Ad = fG4pow->Z13(ResidualA);
      G4double Aj = fG4pow->Z13(A);
      Rb = (1.12*(Aj + Ad) - 0.86*((Aj+Ad)/(Aj*Ad))+2.85)*fermi;
    }
  else if (A>1)
    {
      G4double Ad = fG4pow->Z13(ResidualA);
      G4double Aj = fG4pow->Z13(A);
      Rb=1.5*(Aj+Ad)*fermi;
    }
  else 
    {
      G4double Ad = fG4pow->Z13(ResidualA);
      Rb = 1.5*Ad*fermi;
    }
  G4double GeometricalXS = pi*Rb*Rb; 
    
  G4double ConstantFactor = gg*GeometricalXS*Alpha*pi/(InitialLevelDensity*12);
  //JMQ : this is the assimptotic maximal kinetic energy of the ejectile 
  //      (obtained by energy-momentum conservation). 
  //      In general smaller than U-theSeparationEnergy 
  G4double theEnergy = MaximalKineticEnergy + CoulombBarrier;
  G4double KineticEnergy;
  G4double Probability;

  for(G4int i=0; i<100; ++i) {
    KineticEnergy =  CoulombBarrier + G4UniformRand()*(MaximalKineticEnergy);
    G4double edelta = theEnergy-KineticEnergy-delta0;
    Probability = ConstantFactor*(KineticEnergy + Beta);
    G4double a = 
      theLevelDensityPtr->LevelDensityParameter(ResidualA,ResidualZ,edelta);
    G4double T = 1.0/(std::sqrt(a/Ux) - 1.5/Ux);
    //JMQ fix in units
    
    if (theEnergy - KineticEnergy < Ex) {
      G4double E0 = Ex - T*(G4Log(T) - G4Log(a)*0.25
                - 1.25*G4Log(Ux) + 2.0*std::sqrt(a*Ux));
      Probability *= G4Exp((theEnergy-KineticEnergy-E0)/T)/T;
    } else {
      G4double e2 = edelta*edelta;
      Probability *= 
    G4Exp(2*std::sqrt(a*edelta) - 0.25*G4Log(a*edelta*e2*e2));
    }
    if(EmissionProbability*G4UniformRand() <= Probability) { break; }
  }
    
  return KineticEnergy;
} 

Here is the call graph for this function:

Here is the caller graph for this function:

SetLevelDensityParameter()

void G4GEMChannel::SetLevelDensityParameter ( G4VLevelDensityParameter *  aLevelDensity )

inline

Definition at line 68 of file G4GEMChannel.hh.

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

Here is the call graph for this function:

Member Data Documentation

A

G4int G4GEMChannel::A

private

Definition at line 93 of file G4GEMChannel.hh.

CoulombBarrier

G4double G4GEMChannel::CoulombBarrier

private

Definition at line 112 of file G4GEMChannel.hh.

EmissionProbability

G4double G4GEMChannel::EmissionProbability

private

Definition at line 129 of file G4GEMChannel.hh.

EvaporatedMass

G4double G4GEMChannel::EvaporatedMass

private

Definition at line 98 of file G4GEMChannel.hh.

fG4pow

G4Pow* G4GEMChannel::fG4pow

private

Definition at line 101 of file G4GEMChannel.hh.

MaximalKineticEnergy

G4double G4GEMChannel::MaximalKineticEnergy

private

Definition at line 132 of file G4GEMChannel.hh.

MyOwnLevelDensity

G4bool G4GEMChannel::MyOwnLevelDensity

private

Definition at line 107 of file G4GEMChannel.hh.

pairingCorrection

G4PairingCorrection* G4GEMChannel::pairingCorrection

private

Definition at line 114 of file G4GEMChannel.hh.

ResidualA

G4int G4GEMChannel::ResidualA

private

Definition at line 123 of file G4GEMChannel.hh.

ResidualMass

G4double G4GEMChannel::ResidualMass

private

Definition at line 99 of file G4GEMChannel.hh.

ResidualZ

G4int G4GEMChannel::ResidualZ

private

Definition at line 126 of file G4GEMChannel.hh.

theCoulombBarrierPtr

G4VCoulombBarrier* G4GEMChannel::theCoulombBarrierPtr

private

Definition at line 111 of file G4GEMChannel.hh.

theEvaporationProbabilityPtr

G4GEMProbability* G4GEMChannel::theEvaporationProbabilityPtr

private

Definition at line 104 of file G4GEMChannel.hh.

theLevelDensityPtr

G4VLevelDensityParameter* G4GEMChannel::theLevelDensityPtr

private

Definition at line 108 of file G4GEMChannel.hh.

Z

G4int G4GEMChannel::Z

private

Definition at line 96 of file G4GEMChannel.hh.

---

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

• G4GEMChannel.hh
• G4GEMChannel.cc