Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4FissLib Class Reference

#include <G4FissLib.hh>

Inheritance diagram for G4FissLib:
Collaboration diagram for G4FissLib:

Public Member Functions

 G4FissLib ()
 
 ~G4FissLib ()
 
G4HadFinalStateApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &aTargetNucleus)
 
const std::pair< G4double,
G4double
GetFatalEnergyCheckLevels () const
 
- Public Member Functions inherited from G4HadronicInteraction
 G4HadronicInteraction (const G4String &modelName="HadronicModel")
 
virtual ~G4HadronicInteraction ()
 
virtual G4double SampleInvariantT (const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A)
 
virtual G4bool IsApplicable (const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
 
G4double GetMinEnergy () const
 
G4double GetMinEnergy (const G4Material *aMaterial, const G4Element *anElement) const
 
void SetMinEnergy (G4double anEnergy)
 
void SetMinEnergy (G4double anEnergy, const G4Element *anElement)
 
void SetMinEnergy (G4double anEnergy, const G4Material *aMaterial)
 
G4double GetMaxEnergy () const
 
G4double GetMaxEnergy (const G4Material *aMaterial, const G4Element *anElement) const
 
void SetMaxEnergy (const G4double anEnergy)
 
void SetMaxEnergy (G4double anEnergy, const G4Element *anElement)
 
void SetMaxEnergy (G4double anEnergy, const G4Material *aMaterial)
 
G4int GetVerboseLevel () const
 
void SetVerboseLevel (G4int value)
 
const G4StringGetModelName () const
 
void DeActivateFor (const G4Material *aMaterial)
 
void ActivateFor (const G4Material *aMaterial)
 
void DeActivateFor (const G4Element *anElement)
 
void ActivateFor (const G4Element *anElement)
 
G4bool IsBlocked (const G4Material *aMaterial) const
 
G4bool IsBlocked (const G4Element *anElement) const
 
void SetRecoilEnergyThreshold (G4double val)
 
G4double GetRecoilEnergyThreshold () const
 
virtual std::pair< G4double,
G4double
GetEnergyMomentumCheckLevels () const
 
void SetEnergyMomentumCheckLevels (G4double relativeLevel, G4double absoluteLevel)
 
virtual void ModelDescription (std::ostream &outFile) const
 
virtual void BuildPhysicsTable (const G4ParticleDefinition &)
 
virtual void InitialiseModel ()
 

Additional Inherited Members

- Protected Member Functions inherited from G4HadronicInteraction
void SetModelName (const G4String &nam)
 
G4bool IsBlocked () const
 
void Block ()
 
- Protected Attributes inherited from G4HadronicInteraction
G4HadFinalState theParticleChange
 
G4int verboseLevel
 
G4double theMinEnergy
 
G4double theMaxEnergy
 
G4bool isBlocked
 

Detailed Description

Definition at line 133 of file G4FissLib.hh.

Constructor & Destructor Documentation

G4FissLib::G4FissLib ( )

Definition at line 65 of file G4FissLib.cc.

66  :xSec(0)
67 {
68  SetMinEnergy(0.0);
69  SetMaxEnergy(20.*MeV);
70  if(!getenv("G4NEUTRONHPDATA")) {
71  G4cout << "Please setenv G4NEUTRONHPDATA to point to the neutron cross-section files." << G4endl;
72  throw G4HadronicException(__FILE__, __LINE__, "Please setenv G4NEUTRONHPDATA to point to the neutron cross-section files.");
73  }
74  dirName = getenv("G4NEUTRONHPDATA");
75  G4String tString = "/Fission/";
76  dirName = dirName + tString;
78  theFission = new G4ParticleHPChannel[numEle];
79 
80  for (G4int i=0; i<numEle; i++)
81  {
82 // G4cout << "G4FissLib::G4FissLib(): element "<< i << " : " << (*(G4Element::GetElementTable()))[i]->GetZ()<< G4endl;
83  if((*(G4Element::GetElementTable()))[i]->GetZ()>89)
84  {
85  theFission[i].Init((*(G4Element::GetElementTable()))[i], dirName);
86  theFission[i].Register(&theLibrary);
87  }
88  }
89 }
int G4int
Definition: G4Types.hh:78
G4bool Register(G4ParticleHPFinalState *theFS)
void SetMinEnergy(G4double anEnergy)
G4GLOB_DLL std::ostream G4cout
static size_t GetNumberOfElements()
Definition: G4Element.cc:405
void Init(G4Element *theElement, const G4String dirName)
void SetMaxEnergy(const G4double anEnergy)
#define G4endl
Definition: G4ios.hh:61
static constexpr double MeV
Definition: G4SIunits.hh:214
static G4ElementTable * GetElementTable()
Definition: G4Element.cc:398

Here is the call graph for this function:

G4FissLib::~G4FissLib ( )

Definition at line 91 of file G4FissLib.cc.

92 {
93  delete [] theFission;
94 }

Member Function Documentation

G4HadFinalState * G4FissLib::ApplyYourself ( const G4HadProjectile aTrack,
G4Nucleus aTargetNucleus 
)
virtual

Implements G4HadronicInteraction.

Definition at line 97 of file G4FissLib.cc.

98 {
100 
101  const G4Material* theMaterial = aTrack.GetMaterial();
102  G4int n = theMaterial->GetNumberOfElements();
103  G4int index = theMaterial->GetElement(0)->GetIndex();
104 
105  if (n != 1) {
106  xSec = new G4double[n];
107  G4double sum = 0;
108  G4int i;
109  G4int imat;
110  const G4double * NumAtomsPerVolume = theMaterial->GetVecNbOfAtomsPerVolume();
111  G4double rWeight;
112  G4ParticleHPThermalBoost aThermalE;
113  for (i = 0; i < n; i++) {
114  imat = theMaterial->GetElement(i)->GetIndex();
115  rWeight = NumAtomsPerVolume[i];
116  xSec[i] = theFission[imat].GetXsec(aThermalE.GetThermalEnergy(aTrack,
117  theMaterial->GetElement(i),
118  theMaterial->GetTemperature()));
119  xSec[i] *= rWeight;
120  sum+=xSec[i];
121  }
122 
123  G4double random = G4UniformRand();
124  G4double running = 0;
125  for (i = 0; i < n; i++) {
126  running += xSec[i];
127  index = theMaterial->GetElement(i)->GetIndex();
128  if(random<=running/sum) break;
129  }
130  delete [] xSec;
131  }
132 
133  //return theFission[index].ApplyYourself(aTrack);
134  //Overwrite target parameters
135  G4HadFinalState* result = theFission[index].ApplyYourself(aTrack);
136  aNucleus.SetParameters(G4ParticleHPManager::GetInstance()->GetReactionWhiteBoard()->GetTargA(),G4ParticleHPManager::GetInstance()->GetReactionWhiteBoard()->GetTargZ());
137  const G4Element* target_element = (*G4Element::GetElementTable())[index];
138  const G4Isotope* target_isotope=NULL;
139  G4int iele = target_element->GetNumberOfIsotopes();
140  for ( G4int j = 0 ; j != iele ; j++ ) {
141  target_isotope=target_element->GetIsotope( j );
142  if ( target_isotope->GetN() == G4ParticleHPManager::GetInstance()->GetReactionWhiteBoard()->GetTargA() ) break;
143  }
144  aNucleus.SetIsotope( target_isotope );
146  return result;
147 }
G4double G4ParticleHPJENDLHEData::G4double result
static G4ParticleHPManager * GetInstance()
size_t GetNumberOfIsotopes() const
Definition: G4Element.hh:159
const G4Element * GetElement(G4int iel) const
Definition: G4Material.hh:202
int G4int
Definition: G4Types.hh:78
const G4double * GetVecNbOfAtomsPerVolume() const
Definition: G4Material.hh:206
G4int GetN() const
Definition: G4Isotope.hh:94
#define G4UniformRand()
Definition: Randomize.hh:97
size_t GetIndex() const
Definition: G4Element.hh:182
const G4int n
G4double GetThermalEnergy(const G4HadProjectile &aP, const G4Element *anE, G4double aT)
const G4Isotope * GetIsotope(G4int iso) const
Definition: G4Element.hh:170
G4double GetTemperature() const
Definition: G4Material.hh:182
const G4Material * GetMaterial() const
size_t GetNumberOfElements() const
Definition: G4Material.hh:186
double G4double
Definition: G4Types.hh:76
static G4ElementTable * GetElementTable()
Definition: G4Element.cc:398
G4ParticleHPReactionWhiteBoard * GetReactionWhiteBoard()
G4double GetXsec(G4double energy)
G4HadFinalState * ApplyYourself(const G4HadProjectile &theTrack, G4int isoNumber=-1)

Here is the call graph for this function:

const std::pair< G4double, G4double > G4FissLib::GetFatalEnergyCheckLevels ( ) const
virtual

Reimplemented from G4HadronicInteraction.

Definition at line 149 of file G4FissLib.cc.

150 {
151  // max energy non-conservation is mass of heavy nucleus (taken from G4LFission)
152  return std::pair<G4double, G4double>(5*perCent,250*GeV);
153 }
static constexpr double perCent
Definition: G4SIunits.hh:332
static constexpr double GeV
Definition: G4SIunits.hh:217

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