G4GammaTransition Class Reference

#include <G4GammaTransition.hh>

Public Member Functions
	G4GammaTransition ()
virtual void	SampleDirection (G4Fragment *nuc, G4double ratio, G4int twoJ1, G4int twoJ2, G4int mp)
virtual G4Fragment *	SampleTransition (G4Fragment *nucleus, G4double newExcEnergy, G4double mpRatio, G4int JP1, G4int JP2, G4int MP, G4int shell, G4bool isDiscrete, G4bool isGamma)
void	SetPolarizationFlag (G4bool val)
void	SetTwoJMAX (G4int val)
void	SetVerbose (G4int val)
virtual	~G4GammaTransition ()

Protected Attributes
G4ThreeVector	fDirection
G4PolarizationTransition	fPolTrans
G4int	fTwoJMAX
G4int	fVerbose

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

Private Attributes
G4bool	polarFlag

Detailed Description

Definition at line 52 of file G4GammaTransition.hh.

Constructor & Destructor Documentation

G4GammaTransition() [1/2]

G4GammaTransition::G4GammaTransition ( )

explicit

Definition at line 48 of file G4GammaTransition.cc.

  : polarFlag(false), fDirection(0.,0.,0.), fTwoJMAX(10), fVerbose(0)
{}

~G4GammaTransition()

G4GammaTransition::~G4GammaTransition ( )

virtual

Definition at line 52 of file G4GammaTransition.cc.

{}

G4GammaTransition() [2/2]

G4GammaTransition::G4GammaTransition ( const G4GammaTransition &  right )

privatedelete

Member Function Documentation

operator!=()

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

privatedelete

operator=()

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

privatedelete

operator==()

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

privatedelete

SampleDirection()

void G4GammaTransition::SampleDirection ( G4Fragment *  nuc, G4double  ratio, G4int  twoJ1, G4int  twoJ2, G4int  mp  )

virtual

Definition at line 152 of file G4GammaTransition.cc.

{
  G4double cosTheta, phi;
  G4NuclearPolarization* np = nuc->GetNuclearPolarization(); 
  if(fVerbose > 2) {
    G4cout << "G4GammaTransition::SampleDirection : 2J1= " << twoJ1 
       << " 2J2= " << twoJ2 << " ratio= " << ratio 
       << " mp= " << mp << G4endl;
    G4cout << "  Nucleus: " << *nuc << G4endl;
  }
  if(nullptr == np) {
    cosTheta = 2*G4UniformRand() - 1.0;
    phi = CLHEP::twopi*G4UniformRand();
 
  } else {
    // PhotonEvaporation dataset:
    // The multipolarity number with 1,2,3,4,5,6,7 representing E0,E1,M1,E2,M2,E3,M3
    // monopole transition and 100*Nx+Ny representing multipolarity transition with
    // Ny and Ny taking the value 1,2,3,4,5,6,7 referring to E0,E1,M1,E2,M2,E3,M3,..
    // For example a M1+E2 transition would be written 304.
    // M1 is the primary transition (L) and E2 is the secondary (L')
 
    G4double mpRatio = ratio;
 
    G4int L0 = 0, Lp = 0;
    if (mp > 99) {
      L0 = mp/200;
      Lp = (mp%100)/2;
    } else {
      L0 = mp/2;
      Lp = 0;
      mpRatio = 0.;
    } 
    fPolTrans.SampleGammaTransition(np, twoJ1, twoJ2, L0, Lp, mpRatio, cosTheta, phi);
  }
 
  G4double sinTheta = std::sqrt((1.-cosTheta)*(1.+cosTheta));
  fDirection.set(sinTheta*std::cos(phi),sinTheta*std::sin(phi),cosTheta);
  if(fVerbose > 3) {
    G4cout << "G4GammaTransition::SampleDirection done: " << fDirection << G4endl;
    if(np) { G4cout << *np << G4endl; }
  }
}

References fDirection, fPolTrans, fVerbose, G4cout, G4endl, G4UniformRand, G4InuclParticleNames::nuc, G4PolarizationTransition::SampleGammaTransition(), CLHEP::Hep3Vector::set(), and CLHEP::twopi.

Referenced by SampleTransition().

SampleTransition()

G4Fragment * G4GammaTransition::SampleTransition ( G4Fragment *  nucleus, G4double  newExcEnergy, G4double  mpRatio, G4int  JP1, G4int  JP2, G4int  MP, G4int  shell, G4bool  isDiscrete, G4bool  isGamma  )

virtual

Definition at line 56 of file G4GammaTransition.cc.

{
  G4Fragment* result = nullptr;
  G4double bond_energy = 0.0;
 
  if (!isGamma) { 
    if(0 <= shell) {
      G4int Z = nucleus->GetZ_asInt();
      if(Z <= 100) {
    G4int idx = (G4int)shell;
    idx = std::min(idx, G4AtomicShells::GetNumberOfShells(Z)-1);
    bond_energy = G4AtomicShells::GetBindingEnergy(Z, idx);
      }
    }
  }
  G4double etrans = nucleus->GetExcitationEnergy() - newExcEnergy 
    - bond_energy;
  if(fVerbose > 2) {
    G4cout << "G4GammaTransition::GenerateGamma - Etrans(MeV)= " 
       << etrans << "  Eexnew= " << newExcEnergy 
       << " Ebond= " << bond_energy << G4endl;
  } 
  if(etrans <= 0.0) { 
    etrans += bond_energy;
    bond_energy = 0.0;
  }
  
  // Do complete Lorentz computation 
  G4LorentzVector lv = nucleus->GetMomentum();
  G4double mass = nucleus->GetGroundStateMass() + newExcEnergy;
 
  // select secondary
  G4ParticleDefinition* part;
 
  if(isGamma) { part =  G4Gamma::Gamma(); }
  else {
    part = G4Electron::Electron();
    G4int ne = std::max(nucleus->GetNumberOfElectrons() - 1, 0);
    nucleus->SetNumberOfElectrons(ne);
  }
 
  if(polarFlag && isDiscrete && JP1 <= fTwoJMAX) {
    SampleDirection(nucleus, mpRatio, JP1, JP2, MP);
  } else {
    fDirection = G4RandomDirection();
  }
 
  G4double emass = part->GetPDGMass();
 
  // 2-body decay in rest frame
  G4double ecm       = lv.mag();
  G4ThreeVector bst  = lv.boostVector();
  if(!isGamma) { ecm += (CLHEP::electron_mass_c2 - bond_energy); }
 
  //G4cout << "Ecm= " << ecm << " mass= " << mass << " emass= " << emass << G4endl;
 
  ecm = std::max(ecm, mass + emass);
  G4double energy = 0.5*((ecm - mass)*(ecm + mass) + emass*emass)/ecm;
  G4double mom = (emass > 0.0) ? std::sqrt((energy - emass)*(energy + emass))
    : energy;
 
  // emitted gamma or e-
  G4LorentzVector res4mom(mom * fDirection.x(),
              mom * fDirection.y(),
              mom * fDirection.z(), energy);
  // residual
  energy = std::max(ecm - energy, mass);
  lv.set(-mom*fDirection.x(), -mom*fDirection.y(), -mom*fDirection.z(), energy);
 
  // Lab system transform for short lived level
  lv.boost(bst);
 
  // modified primary fragment 
  nucleus->SetExcEnergyAndMomentum(newExcEnergy, lv);
 
  // gamma or e- are produced
  res4mom.boost(bst); 
  result = new G4Fragment(res4mom, part);
 
  //G4cout << " DeltaE= " << e0 - lv.e() - res4mom.e() + emass
  //     << "   Emass= " << emass << G4endl;
  if(fVerbose > 2) {
    G4cout << "G4GammaTransition::SampleTransition : " << *result << G4endl;
    G4cout << "       Left nucleus: " << *nucleus << G4endl;
  }
  return result;
}

References CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), G4Electron::Electron(), CLHEP::electron_mass_c2, G4INCL::KinematicsUtils::energy(), fDirection, fTwoJMAX, fVerbose, G4cout, G4endl, G4RandomDirection(), G4Gamma::Gamma(), G4AtomicShells::GetBindingEnergy(), G4Fragment::GetExcitationEnergy(), G4Fragment::GetGroundStateMass(), G4Fragment::GetMomentum(), G4Fragment::GetNumberOfElectrons(), G4AtomicShells::GetNumberOfShells(), G4ParticleDefinition::GetPDGMass(), G4Fragment::GetZ_asInt(), CLHEP::HepLorentzVector::mag(), G4INCL::Math::max(), G4INCL::Math::min(), polarFlag, SampleDirection(), CLHEP::HepLorentzVector::set(), G4Fragment::SetExcEnergyAndMomentum(), G4Fragment::SetNumberOfElectrons(), CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), CLHEP::Hep3Vector::z(), and Z.

Referenced by G4PhotonEvaporation::GenerateGamma().

SetPolarizationFlag()

void G4GammaTransition::SetPolarizationFlag ( G4bool  val )

inline

Definition at line 73 of file G4GammaTransition.hh.

{ polarFlag = val; };

References polarFlag.

Referenced by G4PhotonEvaporation::Initialise().

SetTwoJMAX()

void G4GammaTransition::SetTwoJMAX ( G4int  val )

inline

Definition at line 75 of file G4GammaTransition.hh.

{ fTwoJMAX = val; };

References fTwoJMAX.

Referenced by G4PhotonEvaporation::Initialise().

SetVerbose()

void G4GammaTransition::SetVerbose ( G4int  val )

inline

Definition at line 77 of file G4GammaTransition.hh.

{ fVerbose = val; fPolTrans.SetVerbose(val); };

References fPolTrans, fVerbose, and G4PolarizationTransition::SetVerbose().

Referenced by G4PhotonEvaporation::Initialise().

Field Documentation

fDirection

G4ThreeVector G4GammaTransition::fDirection

protected

Definition at line 90 of file G4GammaTransition.hh.

Referenced by SampleDirection(), and SampleTransition().

fPolTrans

G4PolarizationTransition G4GammaTransition::fPolTrans

protected

Definition at line 91 of file G4GammaTransition.hh.

Referenced by SampleDirection(), and SetVerbose().

fTwoJMAX

G4int G4GammaTransition::fTwoJMAX

protected

Definition at line 92 of file G4GammaTransition.hh.

Referenced by SampleTransition(), and SetTwoJMAX().

fVerbose

G4int G4GammaTransition::fVerbose

protected

Definition at line 93 of file G4GammaTransition.hh.

Referenced by SampleDirection(), SampleTransition(), and SetVerbose().

polarFlag

G4bool G4GammaTransition::polarFlag

private

Definition at line 86 of file G4GammaTransition.hh.

Referenced by SampleTransition(), and SetPolarizationFlag().

---

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

• source/processes/hadronic/models/de_excitation/photon_evaporation/include/G4GammaTransition.hh
• source/processes/hadronic/models/de_excitation/photon_evaporation/src/G4GammaTransition.cc