00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026 #include "G4LENDFission.hh"
00027 #include "G4SystemOfUnits.hh"
00028 #include "G4Nucleus.hh"
00029 #include "G4ParticleTable.hh"
00030
00031 G4HadFinalState * G4LENDFission::ApplyYourself(const G4HadProjectile& aTrack, G4Nucleus& aTarg )
00032 {
00033
00034 G4double temp = aTrack.GetMaterial()->GetTemperature();
00035
00036
00037 G4int iZ = aTarg.GetZ_asInt();
00038 G4int iA = aTarg.GetA_asInt();
00039
00040 G4double ke = aTrack.GetKineticEnergy();
00041
00042 G4HadFinalState* theResult = &theParticleChange;
00043 theResult->Clear();
00044
00045 G4GIDI_target* aTarget = usedTarget_map.find( lend_manager->GetNucleusEncoding( iZ , iA ) )->second->GetTarget();
00046 std::vector<G4GIDI_Product>* products = aTarget->getFissionFinalState( ke*MeV, temp, NULL, NULL );
00047 if ( products != NULL )
00048 {
00049 for ( G4int j = 0; j < int( products->size() ); j++ )
00050 {
00051 G4int jZ = (*products)[j].Z;
00052 G4int jA = (*products)[j].A;
00053
00054
00055
00056
00057
00058
00059
00060
00061
00062 G4DynamicParticle* theSec = new G4DynamicParticle;
00063
00064 if ( jZ > 0 )
00065 {
00066
00067 theSec->SetDefinition( G4ParticleTable::GetParticleTable()->FindIon( jZ, jA , 0, 0 ) );
00068 }
00069 else if ( jA == 1 && jZ == 0 )
00070 {
00071 theSec->SetDefinition( G4Neutron::Neutron() );
00072 }
00073 else
00074 {
00075 theSec->SetDefinition( G4Gamma::Gamma() );
00076 }
00077
00078 theSec->SetMomentum( G4ThreeVector( (*products)[j].px*MeV , (*products)[j].py*MeV , (*products)[j].pz*MeV ) );
00079
00080 theResult->AddSecondary( theSec );
00081 }
00082 }
00083 delete products;
00084
00085 theResult->SetStatusChange( stopAndKill );
00086
00087 return theResult;
00088
00089 }