G4NeutronHPFissionBaseFS Class Reference

#include <G4NeutronHPFissionBaseFS.hh>

Inheritance diagram for G4NeutronHPFissionBaseFS:

Public Member Functions
	G4NeutronHPFissionBaseFS ()
virtual	~G4NeutronHPFissionBaseFS ()
void	Init (G4double A, G4double Z, G4int M, G4String &dirName, G4String &bit)
G4DynamicParticleVector *	ApplyYourself (G4int Prompt)
virtual G4double	GetXsec (G4double anEnergy)
virtual G4NeutronHPVector *	GetXsec ()
void	SetNeutron (const G4ReactionProduct &aNeutron)
void	SetTarget (const G4ReactionProduct &aTarget)
Public Member Functions inherited from G4NeutronHPFinalState
	G4NeutronHPFinalState ()
virtual	~G4NeutronHPFinalState ()
void	Init (G4double A, G4double Z, G4String &dirName, G4String &aFSType)
virtual G4NeutronHPFinalState *	New ()=0
G4bool	HasXsec ()
G4bool	HasFSData ()
G4bool	HasAnyData ()
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 40 of file G4NeutronHPFissionBaseFS.hh.

Constructor & Destructor Documentation

G4NeutronHPFissionBaseFS::G4NeutronHPFissionBaseFS ( )

inline

Definition at line 44 of file G4NeutronHPFissionBaseFS.hh.

References G4NeutronHPFinalState::hasXsec.

  { 
    hasXsec = true; 
    theXsection = new G4NeutronHPVector;
  }

virtual G4NeutronHPFissionBaseFS::~G4NeutronHPFissionBaseFS ( )

inlinevirtual

Definition at line 49 of file G4NeutronHPFissionBaseFS.hh.

  {
    delete theXsection;
  }

Member Function Documentation

G4DynamicParticleVector * G4NeutronHPFissionBaseFS::ApplyYourself

(

G4int

Prompt

)

Definition at line 94 of file G4NeutronHPFissionBaseFS.cc.

References G4ReactionProduct::GetKineticEnergy(), G4NeutronHPFinalState::HasFSData(), G4ReactionProduct::Lorentz(), G4Neutron::Neutron(), G4NeutronHPEnergyDistribution::Sample(), G4NeutronHPAngular::SampleAndUpdate(), G4ReactionProduct::SetDefinition(), G4DynamicParticle::SetDefinition(), G4ReactionProduct::SetKineticEnergy(), and G4DynamicParticle::SetMomentum().

Referenced by G4NeutronHPFCFissionFS::ApplyYourself(), G4NeutronHPTCFissionFS::ApplyYourself(), G4NeutronHPSCFissionFS::ApplyYourself(), G4NeutronHPLCFissionFS::ApplyYourself(), and G4NeutronHPFFFissionFS::ApplyYourself().

  {  
// if therere were no data for this isotope, break out.    
    if(!HasFSData()) { return 0; }
    
    G4int i;
    G4DynamicParticleVector * aResult = new G4DynamicParticleVector;
    G4ReactionProduct boosted;
    boosted.Lorentz(theNeutron, theTarget);
    G4double eKinetic = boosted.GetKineticEnergy();
    
// Build neutrons
    G4ReactionProduct * theNeutrons = new G4ReactionProduct[nPrompt];
    for(i=0; i<nPrompt; i++)
    {
      theNeutrons[i].SetDefinition(G4Neutron::Neutron());
    }
    
// sample energies
   G4int dummy;
   for(i=0; i<nPrompt; i++)
   {
     // always in the lab system (if file-5)
     theNeutrons[i].SetKineticEnergy(theEnergyDistribution.Sample(eKinetic, dummy));
   }

// sample neutron angular distribution
   for(i=0; i<nPrompt; i++)
   {
     theAngularDistribution.SampleAndUpdate(theNeutrons[i]);
   }
   
// already in lab. Add neutrons to dynamic particle vector
   for(i=0; i<nPrompt; i++)
   {
      G4DynamicParticle * it = new G4DynamicParticle;
      it->SetDefinition(theNeutrons[i].GetDefinition());
      it->SetMomentum(theNeutrons[i].GetMomentum());
      aResult->push_back(it);
   }
   delete [] theNeutrons;

// return the result
    return aResult;
  }

virtual G4double G4NeutronHPFissionBaseFS::GetXsec ( G4double  anEnergy )

inlinevirtual

Reimplemented from G4NeutronHPFinalState.

Definition at line 58 of file G4NeutronHPFissionBaseFS.hh.

References G4NeutronHPVector::GetY(), and G4INCL::Math::max().

Referenced by G4NeutronHPFissionFS::ApplyYourself().

  {
    return std::max(0., theXsection->GetY(anEnergy));
  }

virtual G4NeutronHPVector* G4NeutronHPFissionBaseFS::GetXsec ( )

inlinevirtual

Reimplemented from G4NeutronHPFinalState.

Definition at line 62 of file G4NeutronHPFissionBaseFS.hh.

{ return theXsection; }

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

virtual

Implements G4NeutronHPFinalState.

Reimplemented in G4NeutronHPLCFissionFS, G4NeutronHPSCFissionFS, and G4NeutronHPTCFissionFS.

Definition at line 43 of file G4NeutronHPFissionBaseFS.cc.

References python.hepunit::eV, G4NeutronHPManager::GetDataStream(), G4NeutronHPManager::GetInstance(), G4NeutronHPNames::GetName(), G4NeutronHPDataUsed::GetName(), G4NeutronHPFinalState::hasAnyData, G4NeutronHPFinalState::hasFSData, G4NeutronHPFinalState::hasXsec, G4NeutronHPAngular::Init(), G4NeutronHPEnergyDistribution::Init(), G4NeutronHPVector::Init(), G4NeutronHPFinalState::SetAZMs(), G4NeutronHPFinalState::theNames, G4NeutronHPFinalState::theNDLDataA, G4NeutronHPFinalState::theNDLDataZ, and G4INCL::CrossSections::total().

Referenced by G4NeutronHPFCFissionFS::Init(), G4NeutronHPTCFissionFS::Init(), G4NeutronHPSCFissionFS::Init(), and G4NeutronHPLCFissionFS::Init().

  {
    G4String tString = dirName;
    G4bool dbool;
    G4NeutronHPDataUsed aFile = theNames.GetName(static_cast<G4int>(A), static_cast<G4int>(Z), M, tString, bit, dbool);
    G4String filename = aFile.GetName();
    SetAZMs( A, Z, M, aFile ); 
    //theBaseA = aFile.GetA();
    //theBaseZ = aFile.GetZ();
    //if(!dbool  || ( Z<2.5 && ( std::abs(theBaseZ - Z)>0.0001 || std::abs(theBaseA - A)>0.0001) ) )
    if ( !dbool || ( Z<2.5 && ( std::abs(theNDLDataZ - Z)>0.0001 || std::abs(theNDLDataA - A)>0.0001)) )
    {
      hasAnyData = false;
      hasFSData = false; 
      hasXsec = false;
      return; // no data for exactly this isotope.
    }

    //std::ifstream theData(filename, std::ios::in);
   std::istringstream theData(std::ios::in);
   G4NeutronHPManager::GetInstance()->GetDataStream(filename,theData);
    G4int dummy;
    if(!(theData))
    {
      //theData.close();
      hasFSData = false;
      hasXsec = false;
      hasAnyData = false;
      return; // no data for this FS for this isotope
    }
    theData >> dummy>>dummy;
    G4int total;
    theData >> total;
    theXsection->Init(theData, total, eV);
    if (!(theData >> dummy))
    {
      hasFSData = false;
      //theData.close();
      return;
    }
    theData >> dummy;

    theAngularDistribution.Init(theData);

    theData >> dummy >> dummy;

    theEnergyDistribution.Init(theData);
    //theData.close();
    
  }

void G4NeutronHPFissionBaseFS::SetNeutron ( const G4ReactionProduct &  aNeutron )

inline

Definition at line 64 of file G4NeutronHPFissionBaseFS.hh.

References G4NeutronHPAngular::SetNeutron().

Referenced by G4NeutronHPFissionFS::ApplyYourself().

                        { 
                          theNeutron = aNeutron;
                          theAngularDistribution.SetNeutron(aNeutron);
                        }

void G4NeutronHPFissionBaseFS::SetTarget ( const G4ReactionProduct &  aTarget )

inline

Definition at line 70 of file G4NeutronHPFissionBaseFS.hh.

References G4NeutronHPAngular::SetTarget().

Referenced by G4NeutronHPFissionFS::ApplyYourself().

                        { 
                          theTarget = aTarget; 
                          theAngularDistribution.SetTarget(aTarget);
                        }

---

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

• G4NeutronHPFissionBaseFS.hh
• G4NeutronHPFissionBaseFS.cc