Geant4_10
G4FissionProductYieldDist.hh
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 /*
27  * File: G4FissionProductYieldDist.hh
28  * Author: B. Wendt (wendbryc@isu.edu)
29  *
30  * Created on May 11, 2011, 12:04 PM
31  */
32 
33 #ifndef G4FISSIONPRODUCTYIELDDIST_HH
34 #define G4FISSIONPRODUCTYIELDDIST_HH
35 
36 #include "G4Ions.hh"
37 #include "G4IonTable.hh"
38 #include "G4NeutronHPNames.hh"
39 #include "globals.hh"
40 #include "G4DynamicParticle.hh"
42 #include "G4ReactionProduct.hh"
43 
44 #include "G4ENDFTapeRead.hh"
46 #include "G4FFGEnumerations.hh"
47 #include "G4FPYNubarValues.hh"
48 #include "G4FPYSamplingOps.hh"
49 #include "G4FPYTreeStructures.hh"
50 
54 {
55 public:
56 // Constructor definition
66  G4FissionProductYieldDist( G4int WhichIsotope,
67  G4FFGEnumerations::MetaState WhichMetaState,
69  G4FFGEnumerations::YieldType WhichYieldType,
70  std::istringstream& dataStream);
81  G4FissionProductYieldDist( G4int WhichIsotope,
82  G4FFGEnumerations::MetaState WhichMetaState,
84  G4FFGEnumerations::YieldType WhichYieldType,
86  std::istringstream& dataStream);
87 private:
89  void Initialize( std::istringstream& dataStream );
90 
91 public:
92 // Functions
114  G4Ions* G4GetFissionProduct( void );
126  void G4SetAlphaProduction( G4double WhatAlphaProduction );
134  void G4SetEnergy( G4double WhatIncidentEnergy );
142  void G4SetTernaryProbability( G4double TernaryProbability );
159  void G4SetVerbosity( G4int WhatVerbosity );
160 
161 protected:
162 // Variables
163  // Class descriptor variables
177 
178  // Datafile variables
183 
184  // Fission generation variables
219 
220  // Pointers to the field of trees and relevant normalization data
243 
244  // Pointers to runtime classes that G4FissionProductYieldDist needs to
245  // function properly
257 
258 //Functions
262  void CheckAlphaSanity( void );
266  G4Ions* FindParticle( G4double RandomParticle );
273  G4Ions* FindParticleExtrapolation( G4double RandomParticle,
274  G4bool LowerEnergyGroupExists );
281  G4Ions* FindParticleInterpolation( G4double RandomParticle,
282  G4int LowerEnergyGroup );
289  G4double RandomParticle,
290  G4int EnergyGroup1,
291  G4int EnergyGroup2 );
294  virtual void GenerateAlphas( std::vector< G4ReactionProduct* >* Alphas );
298  virtual void GenerateNeutrons( std::vector< G4ReactionProduct* >* Neutrons );
302  virtual G4Ions* GetFissionProduct( void ) = 0;
316  G4String MakeDirectoryName( void );
318  G4String MakeFileName( G4int Isotope,
325  G4String MakeIsotopeName( G4int Isotope,
330  virtual void MakeTrees( void );
332  virtual void ReadProbabilities( void );
338  void Renormalize( ProbabilityBranch* Branch );
342  void SampleAlphaEnergies( std::vector< G4ReactionProduct* >* Alphas );
344  void SampleGammaEnergies( std::vector< G4ReactionProduct* >* Gammas );
348  void SampleNeutronEnergies( std::vector< G4ReactionProduct* >* Neutrons );
353  void SetNubar( void );
355  virtual void SortProbability( G4ENDFYieldDataContainer* YieldData );
356 
357 // Destructor function(s)
358 public:
362  virtual ~G4FissionProductYieldDist( void );
363 protected:
365  void BurnTree( ProbabilityBranch* Branch );
366 };
367 
368 #endif /* G4FISSIONPRODUCTYIELDDIST_HH */
369 
const G4FFGEnumerations::YieldType YieldType_
G4FissionProductYieldDist(G4int WhichIsotope, G4FFGEnumerations::MetaState WhichMetaState, G4FFGEnumerations::FissionCause WhichCause, G4FFGEnumerations::YieldType WhichYieldType, std::istringstream &dataStream)
G4Ions * FindParticle(G4double RandomParticle)
void G4SetTernaryProbability(G4double TernaryProbability)
G4DynamicParticleVector * G4GetFission(void)
void G4SetAlphaProduction(G4double WhatAlphaProduction)
G4String MakeFileName(G4int Isotope, G4FFGEnumerations::MetaState MetaState)
virtual void SortProbability(G4ENDFYieldDataContainer *YieldData)
int G4int
Definition: G4Types.hh:78
void G4SetVerbosity(G4int WhatVerbosity)
virtual void GenerateNeutrons(std::vector< G4ReactionProduct * > *Neutrons)
virtual void GenerateAlphas(std::vector< G4ReactionProduct * > *Alphas)
G4DynamicParticle * MakeG4DynamicParticle(G4ReactionProduct *)
void SampleAlphaEnergies(std::vector< G4ReactionProduct * > *Alphas)
Definition: G4Ions.hh:51
bool G4bool
Definition: G4Types.hh:79
std::vector< G4DynamicParticle * > G4DynamicParticleVector
G4Ions * GetParticleDefinition(G4int Product, G4FFGEnumerations::MetaState MetaState)
G4Ions * FindParticleBranchSearch(ProbabilityBranch *Branch, G4double RandomParticle, G4int EnergyGroup1, G4int EnergyGroup2)
void G4SetEnergy(G4double WhatIncidentEnergy)
void Renormalize(ProbabilityBranch *Branch)
void SampleGammaEnergies(std::vector< G4ReactionProduct * > *Gammas)
G4Ions * FindParticleExtrapolation(G4double RandomParticle, G4bool LowerEnergyGroupExists)
G4Ions * FindParticleInterpolation(G4double RandomParticle, G4int LowerEnergyGroup)
const G4FFGEnumerations::FissionCause Cause_
void SampleNeutronEnergies(std::vector< G4ReactionProduct * > *Neutrons)
void BurnTree(ProbabilityBranch *Branch)
double G4double
Definition: G4Types.hh:76
virtual G4Ions * GetFissionProduct(void)=0
G4String MakeIsotopeName(G4int Isotope, G4FFGEnumerations::MetaState MetaState)
const G4FFGEnumerations::MetaState MetaState_