G4VGammaDeexcitation Class Reference

#include <G4VGammaDeexcitation.hh>

Inheritance diagram for G4VGammaDeexcitation:

Public Member Functions
	G4VGammaDeexcitation ()
virtual	~G4VGammaDeexcitation ()
virtual G4VGammaTransition *	CreateTransition ()=0
virtual G4bool	CanDoTransition ()=0
G4FragmentVector *	DoTransition ()
G4FragmentVector *	DoChain ()
G4Fragment *	GenerateGamma ()
G4Fragment *	GetNucleus ()
void	SetNucleus (G4Fragment *nucleus)
void	SetVerboseLevel (G4int verbose)
void	Initialize ()
void	SetEO (G4ElectronOccupancy eo)
void	SetVaccantSN (G4int val)
G4ElectronOccupancy	GetEO ()
G4int	GetVacantSN ()
void	SetTimeLimit (G4double value)

Protected Member Functions
void	Update ()

Protected Attributes
G4VGammaTransition *	_transition
G4int	_verbose

Detailed Description

Definition at line 72 of file G4VGammaDeexcitation.hh.

Constructor & Destructor Documentation

G4VGammaDeexcitation()

G4VGammaDeexcitation::G4VGammaDeexcitation

(

)

Definition at line 76 of file G4VGammaDeexcitation.cc.

                                          : _transition(0), _verbose(0),
                          _electronO (0), _vSN(-1)
{ 
  _nucleus = 0;
  fTimeLimit = DBL_MAX;
}

~G4VGammaDeexcitation()

G4VGammaDeexcitation::~G4VGammaDeexcitation ( )

virtual

Definition at line 83 of file G4VGammaDeexcitation.cc.

{ 
  if (_transition != 0) { delete _transition; }
}

Member Function Documentation

CanDoTransition()

virtual G4bool G4VGammaDeexcitation::CanDoTransition ( )

pure virtual

Implemented in G4ContinuumGammaDeexcitation, and G4DiscreteGammaDeexcitation.

Referenced by DoChain(), DoTransition(), and G4PhotonEvaporation::EmittedFragment().

CreateTransition()

virtual G4VGammaTransition * G4VGammaDeexcitation::CreateTransition ( )

pure virtual

Implemented in G4ContinuumGammaDeexcitation, and G4DiscreteGammaDeexcitation.

Referenced by Initialize(), and Update().

DoChain()

G4FragmentVector * G4VGammaDeexcitation::DoChain

(

)

Definition at line 106 of file G4VGammaDeexcitation.cc.

{
  if (_verbose > 1) { G4cout << "G4VGammaDeexcitation::DoChain" << G4endl; }
  const G4double tolerance = CLHEP::keV;
 
  Initialize();
  G4FragmentVector* products = new G4FragmentVector();
  
  while (CanDoTransition())
    {      
      _transition->SetEnergyFrom(_nucleus->GetExcitationEnergy());
      G4Fragment* gamma = GenerateGamma();
      if (gamma != 0) 
    {
      products->push_back(gamma);
      //G4cout << "Eex(keV)= " << _nucleus->GetExcitationEnergy()/keV << G4endl;
      if(_nucleus->GetExcitationEnergy() <= tolerance) { break; }
      Update();
    }
    } 
  
  if (_verbose > 1) {
    G4cout << "G4VGammaDeexcitation::DoChain - Transition deleted, end of chain " << G4endl;
  }
  
  return products;
}

Referenced by G4PhotonEvaporation::BreakItUp(), and G4PhotonEvaporation::BreakUpFragment().

DoTransition()

G4FragmentVector * G4VGammaDeexcitation::DoTransition

(

)

Definition at line 88 of file G4VGammaDeexcitation.cc.

{
  Initialize();
  G4FragmentVector* products = new G4FragmentVector();
 
  if (CanDoTransition())
    {
      G4Fragment* gamma = GenerateGamma();
      if (gamma != 0) { products->push_back(gamma); }
    }
 
  if (_verbose > 1) {
    G4cout << "G4VGammaDeexcitation::DoTransition - Transition deleted " << G4endl;
  }
 
  return products;
}

Referenced by G4PhotonEvaporation::BreakUp().

GenerateGamma()

G4Fragment * G4VGammaDeexcitation::GenerateGamma

(

)

Definition at line 134 of file G4VGammaDeexcitation.cc.

{
  // 23/04/10 V.Ivanchenko rewrite complitely
  G4double eGamma = 0.;
  
  if (_transition) {
    _transition->SelectGamma();  // it can be conversion electron too
    eGamma = _transition->GetGammaEnergy(); 
    //G4cout << "G4VGammaDeexcitation::GenerateGamma - Egam(MeV)= " 
    //     << eGamma << G4endl; 
    if(eGamma <= 0.0) { return 0; }
  } else { return 0; }
 
  G4double excitation = _nucleus->GetExcitationEnergy() - eGamma;
  if(excitation < 0.0) { excitation = 0.0; } 
  if (_verbose > 1) 
    {
      G4cout << "G4VGammaDeexcitation::GenerateGamma - Edeexc(MeV)= " << eGamma 
         << " ** left Eexc(MeV)= " << excitation
         << G4endl;
    }
 
  G4double gammaTime = _transition->GetGammaCreationTime();
  
  // Do complete Lorentz computation 
  G4LorentzVector lv = _nucleus->GetMomentum();
  G4double Mass = _nucleus->GetGroundStateMass() + excitation;
 
  // select secondary
  G4ParticleDefinition* gamma = G4Gamma::Gamma();
 
  G4DiscreteGammaTransition* dtransition = 
    dynamic_cast <G4DiscreteGammaTransition*> (_transition);
 
  G4bool eTransition = false;
  if (dtransition && !( dtransition->IsAGamma()) ) {
    eTransition = true; 
    gamma = G4Electron::Electron(); 
    _vSN = dtransition->GetOrbitNumber();   
    _electronO.RemoveElectron(_vSN);
    lv += G4LorentzVector(0.0,0.0,0.0,CLHEP::electron_mass_c2 - dtransition->GetBondEnergy());
  }
 
  G4double cosTheta = 1. - 2. * G4UniformRand(); 
  G4double sinTheta = std::sqrt(1. - cosTheta * cosTheta);
  G4double phi = twopi * G4UniformRand();
 
  G4double eMass = gamma->GetPDGMass();
  G4LorentzVector Gamma4P;
 
  //G4cout << "Egamma= " << eGamma << " Mass= " << eMass << " t= " << gammaTime
  //     << " tlim= " << fTimeLimit << G4endl;
 
  if(gammaTime > fTimeLimit) {
    // shortcut for long lived levels
    // not correct position of stopping ion gamma emission
    // 4-momentum balance is breaked
    G4double e = eGamma + eMass;
    G4double mom = std::sqrt(eGamma*(eGamma + 2*eMass));
    Gamma4P.set(mom * sinTheta * std::cos(phi),
        mom * sinTheta * std::sin(phi),
        mom * cosTheta, e); 
    lv -= Gamma4P;
    e = lv.e();
    if(e < Mass) { e = Mass; }
    mom = std::sqrt((e - Mass)*(e + Mass));
    G4ThreeVector v = lv.vect().unit();
    lv.set(mom*v.x(), mom*v.y(), mom*v.z(), e); 
 
  } else {
    // 2-body decay in rest frame
    G4double Ecm       = lv.mag();
    G4ThreeVector bst  = lv.boostVector();
 
    G4double GammaEnergy = 0.5*((Ecm - Mass)*(Ecm + Mass) + eMass*eMass)/Ecm;
    if(GammaEnergy <= eMass) { return 0; }
 
    G4double mom = std::sqrt((GammaEnergy - eMass)*(GammaEnergy + eMass));
    Gamma4P.set(mom * sinTheta * std::cos(phi),
        mom * sinTheta * std::sin(phi),
        mom * cosTheta,
        GammaEnergy);
 
    Gamma4P.boost(bst);  
    lv -= Gamma4P;
  }
 
  // modified primary fragment 
  gammaTime += _nucleus->GetCreationTime();
 
  _nucleus->SetMomentum(lv);
  _nucleus->SetCreationTime(gammaTime);
 
  // e- is not produced
  if(eTransition && !dtransition->GetICM()) { return 0; }
 
  // gamma or e- are produced
  G4Fragment * thePhoton = new G4Fragment(Gamma4P,gamma);
  thePhoton->SetCreationTime(gammaTime);
 
  //G4cout << "G4VGammaDeexcitation::GenerateGamma : " << thePhoton << G4endl;
  //G4cout << "       Left nucleus: " << _nucleus << G4endl;
  return thePhoton;
}

Referenced by DoChain(), DoTransition(), and G4PhotonEvaporation::EmittedFragment().

GetEO()

G4ElectronOccupancy G4VGammaDeexcitation::GetEO ( )

inline

Definition at line 102 of file G4VGammaDeexcitation.hh.

{ return _electronO; };    

Referenced by G4PhotonEvaporation::BreakItUp(), G4PhotonEvaporation::BreakUp(), and G4PhotonEvaporation::BreakUpFragment().

GetNucleus()

G4Fragment * G4VGammaDeexcitation::GetNucleus ( )

inline

Definition at line 128 of file G4VGammaDeexcitation.hh.

{
  return _nucleus; 
}

Referenced by G4ContinuumGammaDeexcitation::CanDoTransition(), G4DiscreteGammaDeexcitation::CanDoTransition(), G4ContinuumGammaDeexcitation::CreateTransition(), and G4DiscreteGammaDeexcitation::CreateTransition().

GetVacantSN()

G4int G4VGammaDeexcitation::GetVacantSN ( )

inline

Definition at line 103 of file G4VGammaDeexcitation.hh.

{return _vSN;};

Referenced by G4PhotonEvaporation::BreakItUp(), G4PhotonEvaporation::BreakUp(), and G4PhotonEvaporation::BreakUpFragment().

Initialize()

void G4VGammaDeexcitation::Initialize ( )

inline

Definition at line 143 of file G4VGammaDeexcitation.hh.

{
  if (_transition != 0) { delete _transition; }
  _transition = CreateTransition();
  if (_transition != 0) {
    _transition->SetEnergyFrom(_nucleus->GetExcitationEnergy());
  }
}

Referenced by DoChain(), DoTransition(), and G4PhotonEvaporation::EmittedFragment().

SetEO()

void G4VGammaDeexcitation::SetEO ( G4ElectronOccupancy  eo )

inline

Definition at line 99 of file G4VGammaDeexcitation.hh.

{ _electronO = eo; };

Referenced by G4PhotonEvaporation::SetEOccupancy().

SetNucleus()

void G4VGammaDeexcitation::SetNucleus ( G4Fragment *  nucleus )

inline

Definition at line 133 of file G4VGammaDeexcitation.hh.

{
  _nucleus = nucleus;
}

Referenced by G4PhotonEvaporation::BreakItUp(), G4PhotonEvaporation::BreakUp(), G4PhotonEvaporation::BreakUpFragment(), and G4PhotonEvaporation::EmittedFragment().

SetTimeLimit()

void G4VGammaDeexcitation::SetTimeLimit ( G4double  value )

inline

Definition at line 105 of file G4VGammaDeexcitation.hh.

{ fTimeLimit = value; }

Referenced by G4PhotonEvaporation::G4PhotonEvaporation(), and G4PhotonEvaporation::SetTimeLimit().

SetVaccantSN()

void G4VGammaDeexcitation::SetVaccantSN ( G4int  val )

inline

Definition at line 100 of file G4VGammaDeexcitation.hh.

{ _vSN = val;};

Referenced by G4PhotonEvaporation::BreakItUp(), G4PhotonEvaporation::BreakUp(), and G4PhotonEvaporation::BreakUpFragment().

SetVerboseLevel()

void G4VGammaDeexcitation::SetVerboseLevel ( G4int  verbose )

inline

Definition at line 138 of file G4VGammaDeexcitation.hh.

{
  _verbose = verbose;
}

Referenced by G4PhotonEvaporation::SetVerboseLevel().

Update()

void G4VGammaDeexcitation::Update ( )

protected

Definition at line 239 of file G4VGammaDeexcitation.cc.

{
  if (_transition !=  0) 
    { 
      delete _transition;
      _transition = 0;
      if (_verbose > 1) {
    G4cout << "G4VGammaDeexcitation::Update - Transition deleted " << G4endl;
      }
    }
  
  _transition = CreateTransition();
  if (_transition != 0) 
    {
      _transition->SetEnergyFrom(_nucleus->GetExcitationEnergy());
      // if ( _vSN != -1) (dynamic_cast <G4DiscreteGammaTransition*> (_transition))->SetICM(false);
      // the above line is commented out for bug fix #952. It was intruduced for reason that
      // the k-shell electron is most likely one to be kicked out and there is no time for 
      // the atom to deexcite before the next IC. But this limitation is causing other problems as 
      // reported in #952
    }
  
  return;
}

Referenced by DoChain().

Member Data Documentation

_transition

G4VGammaTransition* G4VGammaDeexcitation::_transition

protected

Definition at line 111 of file G4VGammaDeexcitation.hh.

Referenced by G4ContinuumGammaDeexcitation::CanDoTransition(), G4DiscreteGammaDeexcitation::CanDoTransition(), DoChain(), GenerateGamma(), Initialize(), Update(), and ~G4VGammaDeexcitation().

_verbose

G4int G4VGammaDeexcitation::_verbose

protected

Definition at line 112 of file G4VGammaDeexcitation.hh.

Referenced by G4ContinuumGammaDeexcitation::CanDoTransition(), G4DiscreteGammaDeexcitation::CanDoTransition(), G4ContinuumGammaDeexcitation::CreateTransition(), G4DiscreteGammaDeexcitation::CreateTransition(), DoChain(), DoTransition(), GenerateGamma(), SetVerboseLevel(), and Update().

---

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

• G4VGammaDeexcitation.hh
• G4VGammaDeexcitation.cc