G4GammaTransition Class Reference

#include <G4GammaTransition.hh>

Inheritance diagram for G4GammaTransition:

Collaboration diagram for G4GammaTransition:

Public Member Functions
	G4GammaTransition ()
virtual	~G4GammaTransition ()
virtual G4Fragment *	SampleTransition (G4Fragment *nucleus, G4double newExcEnergy, G4int deltaS, size_t shell, G4bool isGamma, G4bool isLongLived)

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

Detailed Description

Definition at line 50 of file G4GammaTransition.hh.

Constructor & Destructor Documentation

G4GammaTransition() [1/2]

G4GammaTransition::G4GammaTransition

(

)

Definition at line 48 of file G4GammaTransition.cc.

{}

~G4GammaTransition()

G4GammaTransition::~G4GammaTransition ( )

virtual

Definition at line 51 of file G4GammaTransition.cc.

{}

G4GammaTransition() [2/2]

G4GammaTransition::G4GammaTransition ( const G4GammaTransition &  right )

private

Member Function Documentation

operator!=()

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

private

operator=()

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

private

operator==()

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

private

SampleTransition()

G4Fragment * G4GammaTransition::SampleTransition ( G4Fragment *  nucleus, G4double  newExcEnergy, G4int  deltaS, size_t  shell, G4bool  isGamma, G4bool  isLongLived  )

virtual

Reimplemented in G4PolarizedGammaTransition.

Definition at line 55 of file G4GammaTransition.cc.

{
  G4Fragment* result = 0;
  G4double bond_energy = 0.;

  if (!isGamma) { 
    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;
  //  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;

  //G4double e0 = lv.e();

  // select secondary
  G4ParticleDefinition* part;

  if(isGamma) { part =  G4Gamma::Gamma(); }
  else {
    part = G4Electron::Electron();
    G4int ne = nucleus->GetNumberOfElectrons() - 1;
    if(ne < 0) { ne = 0; } 
    nucleus->SetNumberOfElectrons(ne);
    lv += G4LorentzVector(0.0,0.0,0.0,
                          CLHEP::electron_mass_c2 - bond_energy);
  }

  G4double cosTheta = 1. - 2. * G4UniformRand(); 
  G4double sinTheta = std::sqrt(1. - cosTheta * cosTheta);
  G4double phi = twopi * G4UniformRand();

  G4double emass = part->GetPDGMass();

  // 2-body decay in rest frame
  G4double ecm       = lv.mag();
  G4ThreeVector bst  = lv.boostVector();

  //G4cout << "Ecm= " << ecm << " mass= " << mass << " emass= " << emass 
  //     << "  isLongLived: " << isLongLived << G4endl;

  G4double energy = 0.5*((ecm - mass)*(ecm + mass) + emass*emass)/ecm;
  energy = std::max(energy, emass);

  G4double mom = std::sqrt((energy - emass)*(energy + emass));

  G4LorentzVector res4mom(mom * sinTheta * std::cos(phi),
              mom * sinTheta * std::sin(phi),
              mom * cosTheta, energy);

  // Lab system transform for short lived level
  if(!isLongLived) { 
    res4mom.boost(bst); 
    lv -= res4mom;
  } else {  
    // In exceptional case sample decay at rest at not correct position 
    // of stopping ion, 4-momentum balance is breaked but gamma energy
    // is correct
    lv -= res4mom;
    G4double E = lv.e();
    G4double P2= (E - mass)*(E + mass);
    G4ThreeVector v = lv.vect().unit();
    G4double p = 0.0;
    if(P2 > 0.0) { p = std::sqrt(P2); } 
    else { E = mass; }
    lv.set(v.x()*p, v.y()*p, v.z()*p, E);  
  }

  // modified primary fragment 
  nucleus->SetMomentum(lv);

  // gamma or e- are produced
  result = new G4Fragment(res4mom, part);

  //  G4cout << " DeltaE= " << e0 - lv.e() - res4mom.e() << G4endl;
  //G4cout << "G4GammaTransition::GenerateGamma : " << thePhoton << G4endl;
  //G4cout << "       Left nucleus: " << aNucleus << G4endl;
  return result;
}

Here is the call graph for this function:

Here is the caller graph for this function:

---

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

• G4GammaTransition.hh
• G4GammaTransition.cc