G4ParticleHPFissionBaseFS Class Reference

#include <G4ParticleHPFissionBaseFS.hh>

Inheritance diagram for G4ParticleHPFissionBaseFS:

Public Member Functions
	G4ParticleHPFissionBaseFS ()
	~G4ParticleHPFissionBaseFS () override
void	Init (G4double A, G4double Z, G4int M, G4String &dirName, G4String &bit, G4ParticleDefinition *) override
G4DynamicParticleVector *	ApplyYourself (G4int Prompt)
G4double	GetXsec (G4double anEnergy) override
G4ParticleHPVector *	GetXsec () override
void	SetNeutronRP (const G4ReactionProduct &aNeutron)
void	SetTarget (const G4ReactionProduct &aTarget)
Public Member Functions inherited from G4ParticleHPFinalState
	G4ParticleHPFinalState ()
virtual	~G4ParticleHPFinalState ()
void	Init (G4double A, G4double Z, G4String &dirName, G4String &aFSType, G4ParticleDefinition *p)
virtual G4ParticleHPFinalState *	New ()=0
G4bool	HasXsec ()
G4bool	HasFSData ()
G4bool	HasAnyData ()
void	SetA_Z (G4double anA, G4double aZ, G4int aM=0)
G4double	GetZ ()
G4double	GetN ()
G4double	GetA ()
G4int	GetM ()
void	SetAZMs (G4ParticleHPDataUsed used)
void	SetAZMs (G4double anA, G4double aZ, G4int aM, G4ParticleHPDataUsed used)
void	SetProjectile (G4ParticleDefinition *projectile)
G4ParticleHPFinalState &	operator= (const G4ParticleHPFinalState &right)=delete
	G4ParticleHPFinalState (const G4ParticleHPFinalState &)=delete

Additional Inherited Members
Protected Member Functions inherited from G4ParticleHPFinalState
void	adjust_final_state (G4LorentzVector)
Protected Attributes inherited from G4ParticleHPFinalState
G4ParticleDefinition *	theProjectile {nullptr}
G4ParticleHPManager *	fManager
G4IonTable *	ionTable
G4int	theBaseA {0}
G4int	theBaseZ {0}
G4int	theBaseM {0}
G4int	theNDLDataZ {0}
G4int	theNDLDataA {0}
G4int	theNDLDataM {0}
G4int	secID {-1}
G4bool	hasXsec {true}
G4bool	hasFSData {true}
G4bool	hasAnyData {true}
G4ParticleHPNames	theNames
G4Cache< G4HadFinalState * >	theResult

Detailed Description

Definition at line 42 of file G4ParticleHPFissionBaseFS.hh.

Constructor & Destructor Documentation

G4ParticleHPFissionBaseFS()

G4ParticleHPFissionBaseFS::G4ParticleHPFissionBaseFS ( )

inline

Definition at line 52 of file G4ParticleHPFissionBaseFS.hh.

    {
      hasXsec = true;
      theXsection = new G4ParticleHPVector;
    }

~G4ParticleHPFissionBaseFS()

G4ParticleHPFissionBaseFS::~G4ParticleHPFissionBaseFS ( )

inlineoverride

Definition at line 58 of file G4ParticleHPFissionBaseFS.hh.

{ delete theXsection; }

Member Function Documentation

ApplyYourself()

G4DynamicParticleVector * G4ParticleHPFissionBaseFS::ApplyYourself

(

G4int

Prompt

)

Definition at line 97 of file G4ParticleHPFissionBaseFS.cc.

{
  // if therere were no data for this isotope, break out.
  if (!HasFSData()) {
    return nullptr;
  }
 
  G4int i;
  auto aResult = new G4DynamicParticleVector;
  G4ReactionProduct boosted;
  boosted.Lorentz(*(fCache.Get().theNeutronRP), *(fCache.Get().theTarget));
  G4double eKinetic = boosted.GetKineticEnergy();
 
  // Build neutrons
  std::vector<G4ReactionProduct> theNeutrons;
  for (i = 0; i < nPrompt; i++) {
    theNeutrons.emplace_back();
    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++) {
    auto it = new G4DynamicParticle;
    it->SetDefinition(theNeutrons[i].GetDefinition());
    it->SetMomentum(theNeutrons[i].GetMomentum());
    aResult->push_back(it);
  }
 
  // return the result
  return aResult;
}

Referenced by G4ParticleHPFCFissionFS::ApplyYourself(), G4ParticleHPFFFissionFS::ApplyYourself(), G4ParticleHPLCFissionFS::ApplyYourself(), G4ParticleHPSCFissionFS::ApplyYourself(), and G4ParticleHPTCFissionFS::ApplyYourself().

GetXsec() [1/2]

G4ParticleHPVector * G4ParticleHPFissionBaseFS::GetXsec ( )

inlineoverridevirtual

Reimplemented from G4ParticleHPFinalState.

Definition at line 70 of file G4ParticleHPFissionBaseFS.hh.

{ return theXsection; }

GetXsec() [2/2]

G4double G4ParticleHPFissionBaseFS::GetXsec ( G4double anEnergy )

inlineoverridevirtual

Reimplemented from G4ParticleHPFinalState.

Definition at line 65 of file G4ParticleHPFissionBaseFS.hh.

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

Referenced by G4ParticleHPFissionFS::ApplyYourself().

Init()

void G4ParticleHPFissionBaseFS::Init ( G4double A, G4double Z, G4int M, G4String & dirName, G4String & bit, G4ParticleDefinition *  )

overridevirtual

Implements G4ParticleHPFinalState.

Reimplemented in G4ParticleHPLCFissionFS, G4ParticleHPSCFissionFS, and G4ParticleHPTCFissionFS.

Definition at line 46 of file G4ParticleHPFissionBaseFS.cc.

{
  G4String tString = dirName;
  G4bool dbool;
  G4ParticleHPDataUsed 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);
  G4ParticleHPManager::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();
}

Referenced by G4ParticleHPFCFissionFS::Init(), G4ParticleHPLCFissionFS::Init(), G4ParticleHPSCFissionFS::Init(), and G4ParticleHPTCFissionFS::Init().

SetNeutronRP()

void G4ParticleHPFissionBaseFS::SetNeutronRP ( const G4ReactionProduct & aNeutron )

inline

Definition at line 72 of file G4ParticleHPFissionBaseFS.hh.

    {
      fCache.Get().theNeutronRP = &aNeutron;
      theAngularDistribution.SetProjectileRP(aNeutron);
    }

Referenced by G4ParticleHPFissionFS::ApplyYourself().

SetTarget()

void G4ParticleHPFissionBaseFS::SetTarget ( const G4ReactionProduct & aTarget )

inline

Definition at line 78 of file G4ParticleHPFissionBaseFS.hh.

    {
      fCache.Get().theTarget = &aTarget;
      theAngularDistribution.SetTarget(aTarget);
    }

Referenced by G4ParticleHPFissionFS::ApplyYourself().

---

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

• G4ParticleHPFissionBaseFS.hh
• G4ParticleHPFissionBaseFS.cc