Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Public Member Functions | List of all members
G4LENDFission Class Reference

#include <G4LENDFission.hh>

Inheritance diagram for G4LENDFission:
Inheritance graph
[legend]
Collaboration diagram for G4LENDFission:
Collaboration graph
[legend]

Public Member Functions

 G4LENDFission (G4ParticleDefinition *pd)
 
 ~G4LENDFission ()
 
virtual const std::pair
< G4double, G4double > 		GetFatalEnergyCheckLevels () const
 
G4HadFinalState * ApplyYourself (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 G4String & GetModelName () 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
G4ParticleDefinition * proj
 
G4LENDManager * lend_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();
58  create_used_target_map();
59  };
G4LENDModel::proj
G4ParticleDefinition * proj
Definition: G4LENDModel.hh:79
G4LENDModel::G4LENDModel
G4LENDModel(G4String name="LENDModel")
Definition: G4LENDModel.cc:47
G4LENDModel::create_used_target_map
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 
69  G4DynamicParticle* theSec = new G4DynamicParticle;
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 }
stopAndKill
Definition: G4HadFinalState.hh:42
G4DynamicParticle::SetMomentum
void SetMomentum(const G4ThreeVector &momentum)
Definition: G4DynamicParticle.cc:371
G4HadFinalState::GetSecondary
G4HadSecondary * GetSecondary(size_t i)
Definition: G4HadFinalState.cc:79
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:42
MyRNG
double MyRNG(void *)
Definition: G4LENDModel.cc:45
second
static constexpr double second
Definition: G4SIunits.hh:157
G4int
int G4int
Definition: G4Types.hh:78
G4HadFinalState::SetStatusChange
void SetStatusChange(G4HadFinalStateStatus aS)
Definition: G4HadFinalState.hh:67
G4GIDI_target::getFissionFinalState
std::vector< G4GIDI_Product > * getFissionFinalState(double e_in, double temperature, double(*rng)(void *), void *rngState)
Definition: G4GIDI_target.cc:378
G4LENDManager::GetNucleusEncoding
G4int GetNucleusEncoding(G4int iZ, G4int iA, G4int iM)
Definition: G4LENDManager.cc:343
G4LENDModel::usedTarget_map
std::map< G4int, G4LENDUsedTarget * > usedTarget_map
Definition: G4LENDModel.hh:81
G4GIDI_target
Definition: G4GIDI_target.hh:58
G4HadProjectile::GetKineticEnergy
G4double GetKineticEnergy() const
Definition: G4HadProjectile.hh:106
G4HadProjectile::GetGlobalTime
G4double GetGlobalTime() const
Definition: G4HadProjectile.hh:113
int
typedef int(XMLCALL *XML_NotStandaloneHandler)(void *userData)
G4Neutron::Neutron
static G4Neutron * Neutron()
Definition: G4Neutron.cc:104
G4Gamma::Gamma
static G4Gamma * Gamma()
Definition: G4Gamma.cc:86
G4IonTable::GetIonTable
static G4IonTable * GetIonTable()
Definition: G4IonTable.hh:78
G4Material::GetTemperature
G4double GetTemperature() const
Definition: G4Material.hh:182
MeV
static constexpr double MeV
Definition: G4SIunits.hh:214
G4HadProjectile::GetMaterial
const G4Material * GetMaterial() const
Definition: G4HadProjectile.hh:76
G4HadFinalState::AddSecondary
void AddSecondary(G4DynamicParticle *aP, G4int mod=-1)
Definition: G4HadFinalState.hh:61
G4double
double G4double
Definition: G4Types.hh:76
G4DynamicParticle::SetDefinition
void SetDefinition(const G4ParticleDefinition *aParticleDefinition)
Definition: G4DynamicParticle.cc:305
G4LENDModel::lend_manager
G4LENDManager * lend_manager
Definition: G4LENDModel.hh:80
G4HadFinalState::GetNumberOfSecondaries
G4int GetNumberOfSecondaries() const
Definition: G4HadFinalState.hh:51

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 }
perCent
static constexpr double perCent
Definition: G4SIunits.hh:332
DBL_MAX
#define DBL_MAX
Definition: templates.hh:83
The documentation for this class was generated from the following files: