G4PhotonEvaporation Class Reference

#include <G4PhotonEvaporation.hh>

Inheritance diagram for G4PhotonEvaporation:

Public Member Functions
	G4PhotonEvaporation (const G4String &aName="Anonymous", G4EvaporationChannelType timeType=fDelayedEmission)
virtual	~G4PhotonEvaporation ()
virtual G4Fragment *	EmittedFragment (G4Fragment *theNucleus)
virtual G4FragmentVector *	BreakUpFragment (G4Fragment *theNucleus)
virtual G4FragmentVector *	BreakItUp (const G4Fragment &nucleus)
virtual G4FragmentVector *	BreakUp (const G4Fragment &nucleus)
virtual G4double	GetEmissionProbability (G4Fragment *theNucleus)
virtual void	SetEmissionStrategy (G4VEmissionProbability *probAlgorithm)
void	SetVerboseLevel (G4int verbose)
void	SetICM (G4bool)
void	RDMForced (G4bool)
void	SetMaxHalfLife (G4double)
void	SetTimeLimit (G4double value)
void	SetEOccupancy (G4ElectronOccupancy eOccupancy)
G4ElectronOccupancy	GetEOccupancy ()
G4int	GetVacantShellNumber ()
Public Member Functions inherited from G4VEvaporationChannel
	G4VEvaporationChannel (const G4String &aName="Anonymous", G4EvaporationChannelType timeType=fDelayedEmission)
virtual	~G4VEvaporationChannel ()
virtual G4double	GetLifeTime (G4Fragment *theNucleus)
G4String	GetName () const
void	SetName (const G4String &aName)
void	SetOPTxs (G4int opt)
void	UseSICB (G4bool use)

Additional Inherited Members
Protected Attributes inherited from G4VEvaporationChannel
G4EvaporationChannelType	sampleDecayTime
G4int	OPTxs
G4bool	useSICB

Detailed Description

Definition at line 74 of file G4PhotonEvaporation.hh.

Constructor & Destructor Documentation

G4PhotonEvaporation::G4PhotonEvaporation	(	const G4String &	aName = "Anonymous",
		G4EvaporationChannelType	timeType = fDelayedEmission
	)

Definition at line 72 of file G4PhotonEvaporation.cc.

References DBL_MAX, python.hepunit::second, G4DiscreteGammaDeexcitation::SetICM(), SetMaxHalfLife(), and G4VGammaDeexcitation::SetTimeLimit().

  : G4VEvaporationChannel(aName, timeType),
   verbose(0), myOwnProbAlgorithm(true),
   eOccupancy(0), vShellNumber(-1), gammaE(0.)
{ 
  probAlgorithm = new G4E1Probability;
  contDeexcitation = new G4ContinuumGammaDeexcitation;
  G4DiscreteGammaDeexcitation* p = new G4DiscreteGammaDeexcitation();
  discrDeexcitation = p;

  p->SetICM(false);

  // Time limits
  G4double timeLimit = DBL_MAX;
  char* env = getenv("G4AddTimeLimitToPhotonEvaporation"); 
  if(env) { timeLimit = 1.e-16*second; }

  // Time for short-cut simulation of photon evaporation
  p->SetTimeLimit(timeLimit);

  // Time limit for isomere production 
  SetMaxHalfLife(1.e-6*second);

  nucleus = 0;
}

G4PhotonEvaporation::~G4PhotonEvaporation ( )

virtual

Definition at line 99 of file G4PhotonEvaporation.cc.

{ 
  if(myOwnProbAlgorithm) delete probAlgorithm;
  delete discrDeexcitation;
  delete contDeexcitation;
}

Member Function Documentation

G4FragmentVector * G4PhotonEvaporation::BreakItUp ( const G4Fragment &  nucleus )

virtual

Definition at line 252 of file G4PhotonEvaporation.cc.

References G4VGammaDeexcitation::DoChain(), G4cout, G4endl, G4VGammaDeexcitation::GetEO(), G4VGammaDeexcitation::GetVacantSN(), G4VGammaDeexcitation::SetNucleus(), and G4VGammaDeexcitation::SetVaccantSN().

Referenced by G4NeutronHPCaptureFS::ApplyYourself().

{
  // The same pointer of primary nucleus
  nucleus = new G4Fragment(aNucleus);
  contDeexcitation->SetNucleus(nucleus);
  discrDeexcitation->SetNucleus(nucleus);

  //G4cout << "G4PhotonEvaporation::BreakItUp:  " << nucleus << G4endl;

  // Do the whole gamma chain 
  G4FragmentVector* products = contDeexcitation->DoChain();  
  if( !products ) { products = new G4FragmentVector; }

  // Products from continuum gamma transitions
  if (verbose > 0) {
    G4cout << " = BreakItUp = " << products->size()
       << " gammas from ContinuumDeexcitation " << G4endl;
  }

  // Products from discrete gamma transitions
  G4FragmentVector* discrProducts = discrDeexcitation->DoChain();
  if(discrProducts) {
    eOccupancy = discrDeexcitation->GetEO();
    vShellNumber = discrDeexcitation->GetVacantSN();

    // not sure if the following line is needed!
    discrDeexcitation->SetVaccantSN(-1);

    if (verbose > 0) {
      G4cout << " = BreakItUp = " << discrProducts->size() 
         << " gammas from DiscreteDeexcitation " << G4endl;
    }
    G4FragmentVector::iterator i;
    for (i = discrProducts->begin(); i != discrProducts->end(); ++i)
      {
    products->push_back(*i);
      }
    delete discrProducts;
  }
  // Add deexcited nucleus to products
  products->push_back(nucleus);

  if (verbose > 0) {
    G4cout << "*-*-* Photon evaporation: " << products->size() << G4endl;
  }
  return products;
}

G4FragmentVector * G4PhotonEvaporation::BreakUp ( const G4Fragment &  nucleus )

virtual

Implements G4VEvaporationChannel.

Definition at line 195 of file G4PhotonEvaporation.cc.

References G4VGammaDeexcitation::DoTransition(), G4cout, G4endl, G4VGammaDeexcitation::GetEO(), G4VGammaDeexcitation::GetVacantSN(), G4VGammaDeexcitation::SetNucleus(), and G4VGammaDeexcitation::SetVaccantSN().

Referenced by G4NuclearDecayChannel::DecayIt().

{
  //G4cout << "G4PhotonEvaporation::BreakUp" << G4endl;
  nucleus = new G4Fragment(aNucleus);

  contDeexcitation->SetNucleus(nucleus);
  discrDeexcitation->SetNucleus(nucleus);
  
  // Do one photon emission

  // Products from continuum gamma transitions

  G4FragmentVector* products = contDeexcitation->DoTransition();  
  if( !products ) { products = new G4FragmentVector(); }
  else if(verbose > 0) {
    G4cout << "G4PhotonEvaporation::BreakUp " << products->size() 
       << " gammas from ContinuesDeexcitation " << G4endl;
    G4cout << "   Residual: " << nucleus << G4endl;
  }

  if (0 == products->size())
    {
      // Products from discrete gamma transitions
      G4FragmentVector* discrProducts = discrDeexcitation->DoTransition();

      if (discrProducts) {
    eOccupancy = discrDeexcitation->GetEO();
    vShellNumber = discrDeexcitation->GetVacantSN();

    // not sure if the following line is needed!
    discrDeexcitation->SetVaccantSN(-1);
    //
    if (verbose > 0) {
      G4cout << " = BreakUp = " << discrProducts->size() 
         << " gammas from DiscreteDeexcitation " 
         << G4endl;
      G4cout << "   Residual: " << nucleus << G4endl;
    }
    G4FragmentVector::iterator i;
    for (i = discrProducts->begin(); i != discrProducts->end(); ++i)
      {
        products->push_back(*i);
      }
    delete discrProducts;
      }
    }
  
  // Add deexcited nucleus to products
  products->push_back(nucleus);

  if (verbose > 0) {
    G4cout << "*-*-*-* Photon evaporation: " << products->size() << G4endl;
  }

  return products;
}

G4FragmentVector * G4PhotonEvaporation::BreakUpFragment ( G4Fragment *  theNucleus )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 150 of file G4PhotonEvaporation.cc.

References G4VGammaDeexcitation::DoChain(), G4cout, G4endl, G4VGammaDeexcitation::GetEO(), G4VGammaDeexcitation::GetVacantSN(), G4VGammaDeexcitation::SetNucleus(), and G4VGammaDeexcitation::SetVaccantSN().

Referenced by G4NeutronRadCapture::ApplyYourself().

{
  //G4cout << "G4PhotonEvaporation::BreakUpFragment" << G4endl;
  // The same pointer of primary nucleus
  nucleus = aNucleus;
  contDeexcitation->SetNucleus(nucleus);
  discrDeexcitation->SetNucleus(nucleus);

  // Do the whole gamma chain 
  G4FragmentVector* products = contDeexcitation->DoChain();  
  if( !products ) { products = new G4FragmentVector(); }

  if (verbose > 0) {
    G4cout << "G4PhotonEvaporation::BreakUpFragment " << products->size() 
       << " gammas from ContinuumDeexcitation " << G4endl;
    G4cout << "   Residual: " << nucleus << G4endl;
  }
  // Products from discrete gamma transitions
  G4FragmentVector* discrProducts = discrDeexcitation->DoChain();
  if(discrProducts) {
    eOccupancy = discrDeexcitation->GetEO();
    vShellNumber = discrDeexcitation->GetVacantSN();

    // not sure if the following line is needed!
    discrDeexcitation->SetVaccantSN(-1);

    if (verbose > 0) {
      G4cout << "G4PhotonEvaporation::BreakUpFragment " << discrProducts->size() 
         << " gammas from DiscreteDeexcitation " << G4endl;
      G4cout << "   Residual: " << nucleus << G4endl;
    }
    G4FragmentVector::iterator i;
    for (i = discrProducts->begin(); i != discrProducts->end(); ++i)
      {
    products->push_back(*i);
      }
    delete discrProducts;
  }

  if (verbose > 0) {
    G4cout << "*-*-* Photon evaporation: " << products->size() << G4endl;
  }
  return products;
}

G4Fragment * G4PhotonEvaporation::EmittedFragment ( G4Fragment *  theNucleus )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 106 of file G4PhotonEvaporation.cc.

References G4VGammaDeexcitation::CanDoTransition(), G4cout, G4endl, G4VGammaDeexcitation::GenerateGamma(), G4VGammaDeexcitation::Initialize(), and G4VGammaDeexcitation::SetNucleus().

{
  //G4cout << "G4PhotonEvaporation::EmittedFragment" << G4endl;
  nucleus = aNucleus;
  
  // Do one photon emission by the continues deexcitation  
  contDeexcitation->SetNucleus(nucleus);
  contDeexcitation->Initialize();

  if(contDeexcitation->CanDoTransition()) {  
    G4Fragment* gamma = contDeexcitation->GenerateGamma();
    if(gamma) { 
      if (verbose > 0) {
    G4cout << "G4PhotonEvaporation::EmittedFragment continium deex: "   
           << gamma << G4endl;
        G4cout << "   Residual: " << nucleus << G4endl;
      }
      return gamma; 
    }
  }

  // Do one photon emission by the discrete deexcitation 
  discrDeexcitation->SetNucleus(nucleus);
  discrDeexcitation->Initialize();

  if(discrDeexcitation->CanDoTransition()) {  
    G4Fragment* gamma = discrDeexcitation->GenerateGamma();
    if(gamma) { 
      if (verbose > 0) {
    G4cout << "G4PhotonEvaporation::EmittedFragment discrete deex: "   
           << gamma << G4endl;
        G4cout << "   Residual: " << nucleus << G4endl;
      }
      return gamma; 
    }
  }

  if (verbose > 0) {
    G4cout << "G4PhotonEvaporation unable emit gamma: " 
       << nucleus << G4endl;
  }
  return 0;
}

G4double G4PhotonEvaporation::GetEmissionProbability ( G4Fragment *  theNucleus )

virtual

Implements G4VEvaporationChannel.

Definition at line 301 of file G4PhotonEvaporation.cc.

References G4VEmissionProbability::EmissionProbability(), and G4Fragment::GetExcitationEnergy().

{
  nucleus = theNucleus;
  G4double prob = 
    probAlgorithm->EmissionProbability(*nucleus, nucleus->GetExcitationEnergy());
  return prob;
}

G4ElectronOccupancy G4PhotonEvaporation::GetEOccupancy ( )

inline

Definition at line 107 of file G4PhotonEvaporation.hh.

{ return eOccupancy;} ;

G4int G4PhotonEvaporation::GetVacantShellNumber ( )

inline

Definition at line 109 of file G4PhotonEvaporation.hh.

Referenced by G4NuclearDecayChannel::DecayIt().

{ return vShellNumber;};

void G4PhotonEvaporation::RDMForced

(

G4bool

fromRDM

)

Definition at line 341 of file G4PhotonEvaporation.cc.

Referenced by G4NuclearDecayChannel::DecayIt().

{
  static_cast<G4DiscreteGammaDeexcitation*>(discrDeexcitation)->SetRDM(fromRDM);
}

void G4PhotonEvaporation::SetEmissionStrategy ( G4VEmissionProbability *  probAlgorithm )

virtual

Definition at line 310 of file G4PhotonEvaporation.cc.

{
  // CD - not sure about always wanting to delete this pointer....
  if(myOwnProbAlgorithm) { delete probAlgorithm; }
  probAlgorithm = alg;
  myOwnProbAlgorithm = false;
}

void G4PhotonEvaporation::SetEOccupancy

(

G4ElectronOccupancy

eOccupancy

)

Definition at line 346 of file G4PhotonEvaporation.cc.

References G4VGammaDeexcitation::SetEO().

{
  discrDeexcitation->SetEO(eo);
}

void G4PhotonEvaporation::SetICM

(

G4bool

ic

)

Definition at line 326 of file G4PhotonEvaporation.cc.

Referenced by G4NeutronHPCaptureFS::ApplyYourself(), and G4NuclearDecayChannel::DecayIt().

{
  static_cast<G4DiscreteGammaDeexcitation*>(discrDeexcitation)->SetICM(ic);
}

void G4PhotonEvaporation::SetMaxHalfLife

(

G4double

hl

)

Definition at line 331 of file G4PhotonEvaporation.cc.

Referenced by G4PhotonEvaporation().

{
  static_cast<G4DiscreteGammaDeexcitation*>(discrDeexcitation)->SetHL(hl);
}

void G4PhotonEvaporation::SetTimeLimit

(

G4double

value

)

Definition at line 336 of file G4PhotonEvaporation.cc.

References G4VGammaDeexcitation::SetTimeLimit().

{
  discrDeexcitation->SetTimeLimit(val);
}

void G4PhotonEvaporation::SetVerboseLevel

(

G4int

verbose

)

Definition at line 319 of file G4PhotonEvaporation.cc.

References G4VGammaDeexcitation::SetVerboseLevel().

Referenced by G4NuclearDecayChannel::DecayIt().

{
  verbose = verb;
  contDeexcitation->SetVerboseLevel(verbose);
  discrDeexcitation->SetVerboseLevel(verbose);
}

---

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

• G4PhotonEvaporation.hh
• G4PhotonEvaporation.cc