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

#include <G4LENDFission.hh>

Inheritance diagram for G4LENDFission:
Collaboration diagram for G4LENDFission:

Public Member Functions

 G4LENDFission (G4ParticleDefinition *pd)
 
 ~G4LENDFission ()
 
virtual const std::pair
< G4double, G4double
GetFatalEnergyCheckLevels () const
 
G4HadFinalStateApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &aTargetNucleus)
 
- Public Member Functions inherited from G4LENDModel
 G4LENDModel (G4String name="LENDModel")
 
 ~G4LENDModel ()
 
void ChangeDefaultEvaluation (G4String name)
 
void AllowNaturalAbundanceTarget ()
 
void AllowAnyCandidateTarget ()
 
void BuildPhysicsTable (const G4ParticleDefinition &)
 
- 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 InitialiseModel ()
 

Additional Inherited Members

- Protected Member Functions inherited from G4LENDModel
void create_used_target_map ()
 
void recreate_used_target_map ()
 
- Protected Member Functions inherited from G4HadronicInteraction
void SetModelName (const G4String &nam)
 
G4bool IsBlocked () const
 
void Block ()
 
- Protected Attributes inherited from G4LENDModel
G4ParticleDefinitionproj
 
G4LENDManagerlend_manager
 
std::map< G4int,
G4LENDUsedTarget * > 
usedTarget_map
 
- Protected Attributes inherited from G4HadronicInteraction
G4HadFinalState theParticleChange
 
G4int verboseLevel
 
G4double theMinEnergy
 
G4double theMaxEnergy
 
G4bool isBlocked
 

Detailed Description

Definition at line 46 of file G4LENDFission.hh.

Constructor & Destructor Documentation

G4LENDFission::G4LENDFission ( G4ParticleDefinition pd)
inline

Definition at line 51 of file G4LENDFission.hh.

52  :G4LENDModel( "LENDFission" )
53  {
54  proj = pd;
55 
56 // theModelName = "LENDFission for ";
57 // theModelName += proj->GetParticleName();
59  };
G4ParticleDefinition * proj
Definition: G4LENDModel.hh:79
G4LENDModel(G4String name="LENDModel")
Definition: G4LENDModel.cc:47
void create_used_target_map()
Definition: G4LENDModel.cc:93

Here is the call graph for this function:

G4LENDFission::~G4LENDFission ( )
inline

Definition at line 61 of file G4LENDFission.hh.

61 {;};

Member Function Documentation

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

Reimplemented from G4LENDModel.

Definition at line 31 of file G4LENDFission.cc.

32 {
33 
34  G4double temp = aTrack.GetMaterial()->GetTemperature();
35 
36  //migrate to integer A and Z (GetN_asInt returns number of neutrons in the nucleus since this)
37  G4int iZ = aTarg.GetZ_asInt();
38  G4int iA = aTarg.GetA_asInt();
39  //G4int iM = aTarg.GetM_asInt();
40  G4int iM = 0;
41  if ( aTarg.GetIsotope() != NULL ) {
42  iM = aTarg.GetIsotope()->Getm();
43  }
44 
45  G4double ke = aTrack.GetKineticEnergy();
46 
47  G4HadFinalState* theResult = &theParticleChange;
48  theResult->Clear();
49 
50  G4GIDI_target* aTarget = usedTarget_map.find( lend_manager->GetNucleusEncoding( iZ , iA , iM ) )->second->GetTarget();
51  std::vector<G4GIDI_Product>* products = aTarget->getFissionFinalState( ke*MeV, temp, MyRNG, NULL );
52  if ( products != NULL )
53  {
54  for ( G4int j = 0; j < int( products->size() ); j++ )
55  {
56  G4int jZ = (*products)[j].Z;
57  G4int jA = (*products)[j].A;
58  G4int jM = (*products)[j].m;
59 
60  //G4cout << "Z = " << (*products)[j].Z
61  // << ", A = " << (*products)[j].A
62  // << ", EK = " << (*products)[j].kineticEnergy << " [MeV]"
63  // << ", px = " << (*products)[j].px
64  // << ", py = " << (*products)[j].py
65  // << ", pz = " << (*products)[j].pz
66  // << ", birthTimeSec = " << (*products)[j].birthTimeSec << " [second]"
67  // << G4endl;
68 
70 
71  if ( jZ > 0 )
72  {
73  theSec->SetDefinition( G4IonTable::GetIonTable()->GetIon( jZ, jA , jM ) );
74  }
75  else if ( jA == 1 && jZ == 0 )
76  {
77  theSec->SetDefinition( G4Neutron::Neutron() );
78  }
79  else
80  {
81  theSec->SetDefinition( G4Gamma::Gamma() );
82  }
83 
84  theSec->SetMomentum( G4ThreeVector( (*products)[j].px*MeV , (*products)[j].py*MeV , (*products)[j].pz*MeV ) );
85  //G4cout << theSec->GetDefinition()->GetParticleName() << G4endl;
86  theResult->AddSecondary( theSec );
87  //Set time for delayed neutrons
88  //Current implementation is a little tricky,
89  if ( (*products)[j].birthTimeSec != 0 ) {
90  G4double time = (*products)[j].birthTimeSec*second + aTrack.GetGlobalTime();
91  theResult->GetSecondary(theResult->GetNumberOfSecondaries()-1)->SetTime(time);
92  }
93  }
94  }
95  delete products;
96 
97  theResult->SetStatusChange( stopAndKill );
98 
99  return theResult;
100 
101 }
void SetMomentum(const G4ThreeVector &momentum)
G4HadSecondary * GetSecondary(size_t i)
CLHEP::Hep3Vector G4ThreeVector
double MyRNG(void *)
Definition: G4LENDModel.cc:45
static constexpr double second
Definition: G4SIunits.hh:157
int G4int
Definition: G4Types.hh:78
void SetStatusChange(G4HadFinalStateStatus aS)
std::vector< G4GIDI_Product > * getFissionFinalState(double e_in, double temperature, double(*rng)(void *), void *rngState)
G4int GetNucleusEncoding(G4int iZ, G4int iA, G4int iM)
std::map< G4int, G4LENDUsedTarget * > usedTarget_map
Definition: G4LENDModel.hh:81
G4double GetKineticEnergy() const
G4double GetGlobalTime() const
typedef int(XMLCALL *XML_NotStandaloneHandler)(void *userData)
static G4Neutron * Neutron()
Definition: G4Neutron.cc:104
static G4Gamma * Gamma()
Definition: G4Gamma.cc:86
static G4IonTable * GetIonTable()
Definition: G4IonTable.hh:78
G4double GetTemperature() const
Definition: G4Material.hh:182
static constexpr double MeV
Definition: G4SIunits.hh:214
const G4Material * GetMaterial() const
void AddSecondary(G4DynamicParticle *aP, G4int mod=-1)
double G4double
Definition: G4Types.hh:76
void SetDefinition(const G4ParticleDefinition *aParticleDefinition)
G4LENDManager * lend_manager
Definition: G4LENDModel.hh:80
G4int GetNumberOfSecondaries() const

Here is the call graph for this function:

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

Reimplemented from G4HadronicInteraction.

Definition at line 102 of file G4LENDFission.cc.

103 {
104  // max energy non-conservation is mass of heavy nucleus
105  //return std::pair<G4double, G4double>(5*perCent,250*GeV);
106  return std::pair<G4double, G4double>(5*perCent,DBL_MAX);
107 }
static constexpr double perCent
Definition: G4SIunits.hh:332
#define DBL_MAX
Definition: templates.hh:83

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