G4PreCompoundTransitions Class Reference

#include <G4PreCompoundTransitions.hh>

Inheritance diagram for G4PreCompoundTransitions:

Collaboration diagram for G4PreCompoundTransitions:

Public Member Functions
	G4PreCompoundTransitions ()
virtual	~G4PreCompoundTransitions ()
virtual G4double	CalculateProbability (const G4Fragment &aFragment)
virtual void	PerformTransition (G4Fragment &aFragment)
Public Member Functions inherited from G4VPreCompoundTransitions
	G4VPreCompoundTransitions ()
virtual	~G4VPreCompoundTransitions ()
G4double	GetTransitionProb1 () const
G4double	GetTransitionProb2 () const
G4double	GetTransitionProb3 () const
void	UseNGB (G4bool use)
void	UseCEMtr (G4bool use)

Private Member Functions
	G4PreCompoundTransitions (const G4PreCompoundTransitions &)
const G4PreCompoundTransitions &	operator= (const G4PreCompoundTransitions &right)
G4bool	operator== (const G4PreCompoundTransitions &right) const
G4bool	operator!= (const G4PreCompoundTransitions &right) const

Private Attributes
const G4ParticleDefinition *	proton
G4double	FermiEnergy
G4double	r0
G4double	aLDP

Additional Inherited Members
Protected Attributes inherited from G4VPreCompoundTransitions
G4bool	useNGB
G4bool	useCEMtr
G4double	TransitionProb1
G4double	TransitionProb2
G4double	TransitionProb3

Detailed Description

Definition at line 50 of file G4PreCompoundTransitions.hh.

Constructor & Destructor Documentation

G4PreCompoundTransitions() [1/2]

G4PreCompoundTransitions::G4PreCompoundTransitions

(

)

Definition at line 60 of file G4PreCompoundTransitions.cc.

{
  proton = G4Proton::Proton();
  G4PreCompoundParameters param;
  FermiEnergy = param.GetFermiEnergy();
  r0 = param.GetTransitionsr0();
  aLDP = param.GetLevelDensity();
}

Here is the call graph for this function:

~G4PreCompoundTransitions()

G4PreCompoundTransitions::~G4PreCompoundTransitions ( )

virtual

Definition at line 69 of file G4PreCompoundTransitions.cc.

{}

Here is the call graph for this function:

G4PreCompoundTransitions() [2/2]

G4PreCompoundTransitions::G4PreCompoundTransitions ( const G4PreCompoundTransitions &  )

private

Member Function Documentation

CalculateProbability()

G4double G4PreCompoundTransitions::CalculateProbability ( const G4Fragment &  aFragment )

virtual

Implements G4VPreCompoundTransitions.

Definition at line 75 of file G4PreCompoundTransitions.cc.

{
  // Number of holes
  G4int H = aFragment.GetNumberOfHoles();
  // Number of Particles 
  G4int P = aFragment.GetNumberOfParticles();
  // Number of Excitons 
  G4int N = P+H;
  // Nucleus 
  G4int A = aFragment.GetA_asInt();
  G4int Z = aFragment.GetZ_asInt();
  G4double U = aFragment.GetExcitationEnergy();
  TransitionProb2 = 0.0;
  TransitionProb3 = 0.0;

  //G4cout << aFragment << G4endl;
  
  if(U < 10*eV || 0==N) { return 0.0; }
  
  //J. M. Quesada (Feb. 08) new physics
  // OPT=1 Transitions are calculated according to Gudima's paper 
  //       (original in G4PreCompound from VL) 
  // OPT=2 Transitions are calculated according to Gupta's formulae
  //
  if (useCEMtr) {
    // Relative Energy (T_{rel})
    G4double RelativeEnergy = 1.6*FermiEnergy + U/G4double(N);
    
    // Sample kind of nucleon-projectile 
    G4bool ChargedNucleon(false);
    if(G4int(P*G4UniformRand()) <= aFragment.GetNumberOfCharged()) {
      ChargedNucleon = true; 
    }
    
    // Relative Velocity: 
    // <V_{rel}>^2
    G4double RelativeVelocitySqr;
    if (ChargedNucleon) { 
      RelativeVelocitySqr = 2*RelativeEnergy/CLHEP::proton_mass_c2; 
    } else { 
      RelativeVelocitySqr = 2*RelativeEnergy/CLHEP::neutron_mass_c2; 
    }
    // <V_{rel}>
    G4double RelativeVelocity = std::sqrt(RelativeVelocitySqr);
    
    // Proton-Proton Cross Section
    G4double ppXSection = 
      (10.63/RelativeVelocitySqr - 29.92/RelativeVelocity + 42.9)
      * CLHEP::millibarn;
    // Proton-Neutron Cross Section
    G4double npXSection = 
      (34.10/RelativeVelocitySqr - 82.20/RelativeVelocity + 82.2)
      * CLHEP::millibarn;
    
    // Averaged Cross Section: \sigma(V_{rel})
    G4double AveragedXSection;
    if (ChargedNucleon)
      {
        //JMQ: small bug fixed
        AveragedXSection = ((Z-1)*ppXSection + (A-Z)*npXSection)/G4double(A-1);
      }
    else 
      {
        AveragedXSection = ((A-Z-1)*ppXSection + Z*npXSection)/G4double(A-1);
      }
    
    // Fermi relative energy ratio
    G4double FermiRelRatio = FermiEnergy/RelativeEnergy;
    
    // This factor is introduced to take into account the Pauli principle
    G4double PauliFactor = 1.0 - 1.4*FermiRelRatio;
    if (FermiRelRatio > 0.5) {
      G4double x = 2.0 - 1.0/FermiRelRatio;
      PauliFactor += 0.4*FermiRelRatio*x*x*std::sqrt(x);
    }
    // Interaction volume 
    G4double xx = 2*r0 + CLHEP::hbarc/(CLHEP::proton_mass_c2*RelativeVelocity);
    G4double Vint = CLHEP::pi*xx*xx*xx/0.75;
    
    // Transition probability for \Delta n = +2
    
    TransitionProb1 = std::max(0.0, AveragedXSection*PauliFactor
      *std::sqrt(2.0*RelativeEnergy/CLHEP::proton_mass_c2)/Vint);

    //JMQ 281009  phenomenological factor in order to increase 
    //   equilibrium contribution
    //   G4double factor=5.0;
    //   TransitionProb1 *= factor;
    
    // GE = g*E where E is Excitation Energy
    G4double GE = sixdpi2*aLDP*A*U;
    G4double Fph = G4double(P*P+H*H+P-3*H)*0.25;
    
    if(!useNGB) { 
        
      // F(p+1,h+1)
      G4double Fph1 = Fph + N*0.5;

      static const G4double plimit = 100;

      //JMQ/AH  bug fixed: if (U-Fph < 0.0) 
      if (GE-Fph1 > 0.0) { 
        G4double x0 = GE-Fph;
    G4double x1 = (N+1)*G4Log(x0/(GE-Fph1));
    if(x1 < plimit) {
      x1 = G4Exp(x1)*TransitionProb1/x0;
    
      // Transition probability for \Delta n = -2 (at F(p,h) = 0)
      TransitionProb2 = std::max(0.0, (P*H*(N+1)*(N-2))*x1/x0);
        
      // Transition probability for \Delta n = 0 (at F(p,h) = 0)
      TransitionProb3 = std::max(0.0,((N+1)*(P*(P-1) + 4*P*H + H*(H-1)))*x1
                     /G4double(N));
    }
      }
    }

  } else {
    //JMQ: Transition probabilities from Gupta's work    
    // GE = g*E where E is Excitation Energy
    TransitionProb1 = std::max(0.0, U*(4.2e+12 - 3.6e+10*U/G4double(N+1)))
      /(16*CLHEP::c_light); 

    if (!useNGB && N > 1) {
      G4double GE = sixdpi2*aLDP*A*U; 
      TransitionProb2 = ((N-1)*(N-2)*P*H)*TransitionProb1/(GE*GE);  
    }
  }
  //  G4cout<<"U = "<<U<<G4endl;
  //  G4cout<<"N="<<N<<"  P="<<P<<"  H="<<H<<G4endl;
  //  G4cout<<"l+ ="<<TransitionProb1<<"  l- ="<< TransitionProb2
  //   <<"  l0 ="<< TransitionProb3<<G4endl; 
  return TransitionProb1 + TransitionProb2 + TransitionProb3;
}

Here is the call graph for this function:

Here is the caller graph for this function:

operator!=()

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

private

operator=()

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

private

operator==()

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

private

PerformTransition()

void G4PreCompoundTransitions::PerformTransition ( G4Fragment &  aFragment )

virtual

Implements G4VPreCompoundTransitions.

Definition at line 210 of file G4PreCompoundTransitions.cc.

{
  G4double ChosenTransition = 
    G4UniformRand()*(TransitionProb1 + TransitionProb2 + TransitionProb3);
  G4int deltaN = 0;
  G4int Npart     = result.GetNumberOfParticles();
  G4int Ncharged  = result.GetNumberOfCharged();
  G4int Nholes    = result.GetNumberOfHoles();
  if (ChosenTransition <= TransitionProb1) 
    {
      // Number of excitons is increased on \Delta n = +2
      deltaN = 2;
    } 
  else if (ChosenTransition <= TransitionProb1+TransitionProb2) 
    {
      // Number of excitons is increased on \Delta n = -2
      deltaN = -2;
    }

  // AH/JMQ: Randomly decrease the number of charges if deltaN is -2 and  
  // in proportion to the number charges w.r.t. number of particles,  
  // PROVIDED that there are charged particles
  deltaN /= 2;

  //G4cout << "deltaN= " << deltaN << G4endl;

  // JMQ the following lines have to be before SetNumberOfCharged, 
  //     otherwise the check on number of charged vs. number of particles fails
  result.SetNumberOfParticles(Npart+deltaN);
  result.SetNumberOfHoles(Nholes+deltaN); 

  if(deltaN < 0) {
    if( (Ncharged == Npart) ||
    (Ncharged >= 1 && G4int(Npart*G4UniformRand()) <= Ncharged)) 
      { 
    result.SetNumberOfCharged(Ncharged+deltaN); // deltaN is negative!
      }

  } else if ( deltaN > 0 ) {
    // With weight Z/A, number of charged particles is increased with +1
    G4int A = result.GetA_asInt() - Npart;
    G4int Z = result.GetZ_asInt() - Ncharged;
    if((Z == A) ||  (Z > 0 && G4int(A*G4UniformRand()) <= Z)) 
      {
    result.SetNumberOfCharged(Ncharged+deltaN);
      }
  }
  
  // Number of charged can not be greater that number of particles
  if ( Npart < Ncharged ) 
    {
      result.SetNumberOfCharged(Npart);
    }
  //G4cout << "### After transition" << G4endl;
  //G4cout << result << G4endl;
}

Here is the call graph for this function:

Member Data Documentation

aLDP

G4double G4PreCompoundTransitions::aLDP

private

Definition at line 73 of file G4PreCompoundTransitions.hh.

FermiEnergy

G4double G4PreCompoundTransitions::FermiEnergy

private

Definition at line 71 of file G4PreCompoundTransitions.hh.

proton

const G4ParticleDefinition* G4PreCompoundTransitions::proton

private

Definition at line 69 of file G4PreCompoundTransitions.hh.

r0

G4double G4PreCompoundTransitions::r0

private

Definition at line 72 of file G4PreCompoundTransitions.hh.

---

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

• G4PreCompoundTransitions.hh
• G4PreCompoundTransitions.cc