#include <G4PromptPhotonEvaporation.hh>

Inheritance diagram for G4PromptPhotonEvaporation:

Public Member Functions
	G4PromptPhotonEvaporation ()

virtual	~G4PromptPhotonEvaporation ()

virtual G4double	GetEmissionProbability (G4Fragment *theNucleus)

virtual G4Fragment *	EmittedFragment (G4Fragment *theNucleus)

virtual G4FragmentVector *	BreakUpFragment (G4Fragment *theNucleus)

virtual G4FragmentVector *	BreakUp (const G4Fragment &theNucleus)

void	SetVerboseLevel (G4int verbose)

void	SetICM (G4bool)

void	RDMForced (G4bool)

void	SetMaxHalfLife (G4double)

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 55 of file G4PromptPhotonEvaporation.hh.

Constructor & Destructor Documentation

G4PromptPhotonEvaporation::G4PromptPhotonEvaporation ( )

Definition at line 64 of file G4PromptPhotonEvaporation.cc.

References G4NuclearLevelStore::GetInstance().

   :fVerbose(0), fICM(true), fRDM(false), fMaxHalfTime(DBL_MAX),
    fEmissionProbability(0.0),levelManager(0),nucleus(0)
 {
   fNuclearLevelStore = G4NuclearLevelStore::GetInstance(); 
   theA = theZ = 0;
   fEnergyFermi = fExcEnergyMax = gammaE = 0.0;
 }

G4PromptPhotonEvaporation::~G4PromptPhotonEvaporation ( )

virtual

Definition at line 73 of file G4PromptPhotonEvaporation.cc.

74 {

75 }

Member Function Documentation

G4FragmentVector * G4PromptPhotonEvaporation::BreakUp ( const G4Fragment & theNucleus )

virtual

Implements G4VEvaporationChannel.

Definition at line 171 of file G4PromptPhotonEvaporation.cc.

References BreakUpFragment(), and test::v.

 {
   //G4cout << "G4PromptPhotonEvaporation::BreakUp" << G4endl;
   G4Fragment* initialState = new G4Fragment(theNucleus);
   G4FragmentVector* v = BreakUpFragment(initialState);
   v->push_back(initialState);
   return v;
 }

G4FragmentVector * G4PromptPhotonEvaporation::BreakUpFragment ( G4Fragment * theNucleus )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 157 of file G4PromptPhotonEvaporation.cc.

References EmittedFragment(), and test::v.

Referenced by BreakUp().

 {
   //G4cout << "G4PromptPhotonEvaporation::BreakUpFragment" << G4endl;
   G4FragmentVector* v = new G4FragmentVector();
   G4Fragment* gamma = 0;
   do {
     gamma = EmittedFragment(theNucleus);
     if(gamma) { v->push_back(gamma); }
   } while(gamma);
 
   return v;
 }

G4Fragment * G4PromptPhotonEvaporation::EmittedFragment ( G4Fragment * theNucleus )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 111 of file G4PromptPhotonEvaporation.cc.

References GetEmissionProbability(), G4Fragment::GetExcitationEnergy(), and python.hepunit::keV.

Referenced by BreakUpFragment().

 {
   //G4cout << "G4PromptPhotonEvaporation::EmittedFragment" << G4endl;
   
   G4Fragment* gamma = 0;
   if(theNucleus->GetExcitationEnergy() <= keV) { return gamma; }
   if(GetEmissionProbability(theNucleus) <= 0.0){ return gamma; }
 
   //G4cout << "G4PromptPhotonEvaporation::EmittedFragment" << G4endl;
   /*
     G4Fragment* gamma = _contDeexcitation->GenerateGamma();
     if(gamma) { 
       if (_verbose > 0) {
     G4cout << "G4PromptPhotonEvaporation::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 << "G4PromptPhotonEvaporation::EmittedFragment discrete deex: "   
            << gamma << G4endl;
         G4cout << "   Residual: " << nucleus << G4endl;
       }
       return gamma; 
     }
   }
  
   if (_verbose > 0) {
     G4cout << "G4PromptPhotonEvaporation unable emit gamma: " 
        << nucleus << G4endl;
   }
   */
   return gamma;
 }

G4double G4PromptPhotonEvaporation::GetEmissionProbability ( G4Fragment * theNucleus )

virtual

Implements G4VEvaporationChannel.

Definition at line 78 of file G4PromptPhotonEvaporation.cc.

References G4Fragment::GetA_asInt(), G4Fragment::GetExcitationEnergy(), G4Fragment::GetGroundStateMass(), G4NuclearLevelStore::GetManager(), G4NucleiProperties::GetNuclearMass(), G4Fragment::GetZ_asInt(), python.hepunit::keV, python.hepunit::MeV, and python.hepunit::neutron_mass_c2.

Referenced by EmittedFragment().

 {
   fEmissionProbability = 0.0;
   nucleus = theNucleus;
   G4double ex = nucleus->GetExcitationEnergy();
 
   if(nucleus->GetZ_asInt() != theZ || nucleus->GetA_asInt() != theA) {
     G4int Z = nucleus->GetZ_asInt();
     G4int A = nucleus->GetA_asInt();
     fExcEnergyMax = -1.0;
     if(1 < A && ex > keV) {
       fEnergyFermi = G4NucleiProperties::GetNuclearMass(A-1, Z) 
     + neutron_mass_c2 - nucleus->GetGroundStateMass();
       fExcEnergyMax = fEnergyFermi + 15*MeV;
     }
     if(ex < fExcEnergyMax) {
 
       theZ = Z;
       theA = A;
       levelManager = fNuclearLevelStore->GetManager(Z,A);
 
       // continium transition
       if(ex >= fEnergyFermi) {
 
     // discrete transition
       } else {
       }
     }
   }
   return fEmissionProbability;
 }