G4FissionLibrary Class Reference

#include <G4FissionLibrary.hh>

Inheritance diagram for G4FissionLibrary:

Public Member Functions
	G4FissionLibrary ()
	~G4FissionLibrary ()
void	Init (G4double A, G4double Z, G4int M, G4String &dirName, G4String &)
G4HadFinalState *	ApplyYourself (const G4HadProjectile &theTrack)
G4NeutronHPFinalState *	New ()
Public Member Functions inherited from G4NeutronHPFinalState
	G4NeutronHPFinalState ()
virtual	~G4NeutronHPFinalState ()
void	Init (G4double A, G4double Z, G4String &dirName, G4String &aFSType)
G4bool	HasXsec ()
G4bool	HasFSData ()
G4bool	HasAnyData ()
virtual G4double	GetXsec (G4double)
virtual G4NeutronHPVector *	GetXsec ()
void	SetA_Z (G4double anA, G4double aZ, G4int aM=0)
G4double	GetZ ()
G4double	GetN ()
G4int	GetM ()

Additional Inherited Members
Protected Member Functions inherited from G4NeutronHPFinalState
void	SetAZMs (G4double anA, G4double aZ, G4int aM, G4NeutronHPDataUsed used)
void	adjust_final_state (G4LorentzVector)
G4bool	DoNotAdjustFinalState ()
Protected Attributes inherited from G4NeutronHPFinalState
G4bool	hasXsec
G4bool	hasFSData
G4bool	hasAnyData
G4NeutronHPNames	theNames
G4HadFinalState	theResult
G4double	theBaseA
G4double	theBaseZ
G4int	theBaseM
G4int	theNDLDataZ
G4int	theNDLDataA
G4int	theNDLDataM

Detailed Description

Definition at line 76 of file G4FissionLibrary.hh.

Constructor & Destructor Documentation

G4FissionLibrary::G4FissionLibrary

(

)

Definition at line 66 of file G4FissionLibrary.cc.

References G4NeutronHPFinalState::hasXsec.

Referenced by New().

  : G4NeutronHPFinalState(), theIsotope(0), targetMass(0.0)
{
  hasXsec = false;
}

G4FissionLibrary::~G4FissionLibrary

(

)

Definition at line 72 of file G4FissionLibrary.cc.

{}

Member Function Documentation

G4HadFinalState * G4FissionLibrary::ApplyYourself ( const G4HadProjectile &  theTrack )

virtual

Reimplemented from G4NeutronHPFinalState.

Definition at line 138 of file G4FissionLibrary.cc.

References G4HadFinalState::AddSecondary(), G4HadFinalState::Clear(), G4Gamma::Gamma(), G4HadProjectile::Get4Momentum(), G4Nucleus::GetBiasedThermalNucleus(), G4HadProjectile::GetDefinition(), G4ReactionProduct::GetDefinition(), G4NeutronHPFissionERelease::GetFragmentKinetic(), G4HadProjectile::GetKineticEnergy(), G4ReactionProduct::GetKineticEnergy(), G4HadProjectile::GetMaterial(), G4ReactionProduct::GetMomentum(), getndircosu_(), getndircosv_(), getndircosw_(), getneng_(), G4ParticleDefinition::GetPDGMass(), getpdircosu_(), getpdircosv_(), getpdircosw_(), getpeng_(), G4Material::GetTemperature(), G4ReactionProduct::GetTotalMomentum(), G4DynamicParticle::GetTotalMomentum(), G4ReactionProduct::Lorentz(), python.hepunit::MeV, G4Neutron::Neutron(), G4ReactionProduct::SetDefinition(), G4DynamicParticle::SetDefinition(), G4ReactionProduct::SetKineticEnergy(), G4DynamicParticle::SetKineticEnergy(), G4HadFinalState::SetLocalEnergyDeposit(), G4ReactionProduct::SetMomentum(), G4DynamicParticle::SetMomentum(), G4NeutronHPAngular::SetNeutron(), G4HadFinalState::SetStatusChange(), G4NeutronHPAngular::SetTarget(), stopAndKill, G4NeutronHPFinalState::theResult, theTarget, and CLHEP::HepLorentzVector::vect().

{  
  theResult.Clear();

// prepare neutron
  G4double eKinetic = theTrack.GetKineticEnergy();
  const G4HadProjectile *incidentParticle = &theTrack;
  G4ReactionProduct theNeutron( const_cast<G4ParticleDefinition *>(incidentParticle->GetDefinition()) );
  theNeutron.SetMomentum( incidentParticle->Get4Momentum().vect() );
  theNeutron.SetKineticEnergy( eKinetic );

// prepare target
  G4Nucleus aNucleus;
  G4ReactionProduct theTarget; 
  G4ThreeVector neuVelo = (1./incidentParticle->GetDefinition()->GetPDGMass())*theNeutron.GetMomentum();
  theTarget = aNucleus.GetBiasedThermalNucleus( targetMass, neuVelo, theTrack.GetMaterial()->GetTemperature());

// set neutron and target in the FS classes 
  theNeutronAngularDis.SetNeutron(theNeutron);
  theNeutronAngularDis.SetTarget(theTarget);

// boost to target rest system
  theNeutron.Lorentz(theNeutron, -1*theTarget);

  eKinetic = theNeutron.GetKineticEnergy();    

// dice neutron and gamma multiplicities, energies and momenta in Lab. @@
// no energy conservation on an event-to-event basis. we rely on the data to be ok. @@
// also for mean, we rely on the consistency of the data. @@

  G4int nPrompt=0, gPrompt=0;
  SampleMult(theTrack, &nPrompt, &gPrompt, eKinetic);

// Build neutrons and add them to dynamic particle vector
  G4double momentum;
  for(G4int i=0; i<nPrompt; i++)
  {
    G4DynamicParticle * it = new G4DynamicParticle;
    it->SetDefinition(G4Neutron::Neutron());
    it->SetKineticEnergy(getneng_(&i)*MeV);
    momentum = it->GetTotalMomentum();
    G4ThreeVector temp(momentum*getndircosu_(&i), 
                       momentum*getndircosv_(&i), 
                       momentum*getndircosw_(&i));
    it->SetMomentum( temp );
//    it->SetGlobalTime(getnage_(&i)*second);
    theResult.AddSecondary(it);
//    G4cout <<"G4FissionLibrary::ApplyYourself: energy of prompt neutron " << i << " = " << it->GetKineticEnergy()<<G4endl;
  }

// Build gammas, lorentz transform them, and add them to dynamic particle vector
  for(G4int i=0; i<gPrompt; i++)
  {
    G4ReactionProduct * thePhoton = new G4ReactionProduct;
    thePhoton->SetDefinition(G4Gamma::Gamma());
    thePhoton->SetKineticEnergy(getpeng_(&i)*MeV);
    momentum = thePhoton->GetTotalMomentum();
    G4ThreeVector temp(momentum*getpdircosu_(&i), 
                       momentum*getpdircosv_(&i), 
                       momentum*getpdircosw_(&i));
    thePhoton->SetMomentum( temp );
    thePhoton->Lorentz(*thePhoton, -1.*theTarget);
    
    G4DynamicParticle * it = new G4DynamicParticle;
    it->SetDefinition(thePhoton->GetDefinition());
    it->SetMomentum(thePhoton->GetMomentum());
//    it->SetGlobalTime(getpage_(&i)*second);
//    G4cout <<"G4FissionLibrary::ApplyYourself: energy of prompt photon " << i << " = " << it->GetKineticEnergy()<<G4endl;
    theResult.AddSecondary(it);
    delete thePhoton;  
  }
//  G4cout <<"G4FissionLibrary::ApplyYourself: Number of secondaries = "<<theResult.GetNumberOfSecondaries()<< G4endl;
//  G4cout <<"G4FissionLibrary::ApplyYourself: Number of induced prompt neutron = "<<nPrompt<<G4endl;
//  G4cout <<"G4FissionLibrary::ApplyYourself: Number of induced prompt photons = "<<gPrompt<<G4endl;

// finally deal with local energy depositions.
  G4double eDepByFragments = theEnergyRelease.GetFragmentKinetic();
  theResult.SetLocalEnergyDeposit(eDepByFragments);
//   G4cout << "G4FissionLibrary::local energy deposit" << eDepByFragments<<G4endl;
// clean up the primary neutron
  theResult.SetStatusChange(stopAndKill);
  return &theResult;
}

void G4FissionLibrary::Init ( G4double  A, G4double  Z, G4int  M, G4String &  dirName, G4String &    )

virtual

Implements G4NeutronHPFinalState.

Definition at line 82 of file G4FissionLibrary.cc.

References G4cout, G4endl, G4NeutronHPNames::GetName(), G4NeutronHPDataUsed::GetName(), G4NeutronHPNeutronYield::GetTargetMass(), G4NeutronHPFinalState::hasAnyData, G4NeutronHPFinalState::hasFSData, G4NeutronHPFinalState::hasXsec, G4NeutronHPFissionERelease::Init(), G4NeutronHPAngular::Init(), G4NeutronHPEnergyDistribution::Init(), G4NeutronHPPhotonDist::InitAngular(), G4NeutronHPNeutronYield::InitDelayed(), G4NeutronHPPhotonDist::InitEnergies(), G4NeutronHPNeutronYield::InitMean(), G4NeutronHPPhotonDist::InitMean(), G4NeutronHPNeutronYield::InitPrompt(), and G4NeutronHPFinalState::theNames.

{
  G4String tString = "/FS/";
  G4bool dbool;
  theIsotope = static_cast<G4int>(1000*Z+A);
  G4NeutronHPDataUsed aFile = theNames.GetName(static_cast<G4int>(A), static_cast<G4int>(Z), M, dirName, tString, dbool);
  G4String filename = aFile.GetName();

  if(!dbool)
  {
    hasAnyData = false;
    hasFSData = false;
    hasXsec = false;
    return;
  }
  std::ifstream theData(filename, std::ios::in);

  // here it comes
  G4int infoType, dataType;
  hasFSData = false;
  while (theData >> infoType)
  {
    hasFSData = true;
    theData >> dataType;
    switch(infoType)
    {
      case 1:
        if(dataType==4) theNeutronAngularDis.Init(theData);
        if(dataType==5) thePromptNeutronEnDis.Init(theData);
        if(dataType==12) theFinalStatePhotons.InitMean(theData);
        if(dataType==14) theFinalStatePhotons.InitAngular(theData);
        if(dataType==15) theFinalStatePhotons.InitEnergies(theData);
        break;
      case 2:
        if(dataType==1) theFinalStateNeutrons.InitMean(theData);
        break;
      case 3:
        if(dataType==1) theFinalStateNeutrons.InitDelayed(theData);
        if(dataType==5) theDelayedNeutronEnDis.Init(theData);
        break;
      case 4:
        if(dataType==1) theFinalStateNeutrons.InitPrompt(theData);
        break;
      case 5:
        if(dataType==1) theEnergyRelease.Init(theData);
        break;
      default:
        G4cout << "G4FissionLibrary::Init: unknown data type"<<dataType<<G4endl;
        throw G4HadronicException(__FILE__, __LINE__, "G4FissionLibrary::Init: unknown data type");
        break;
    }
  }
  targetMass = theFinalStateNeutrons.GetTargetMass();
  theData.close();
}

G4NeutronHPFinalState * G4FissionLibrary::New ( )

virtual

Implements G4NeutronHPFinalState.

Definition at line 75 of file G4FissionLibrary.cc.

References G4FissionLibrary().

{
  G4FissionLibrary * theNew = new G4FissionLibrary;
  return theNew;
}

---

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

• G4FissionLibrary.hh
• G4FissionLibrary.cc