G4ParticleHPFissionBaseFS Class Reference

#include <G4ParticleHPFissionBaseFS.hh>

Inheritance diagram for G4ParticleHPFissionBaseFS:

Public Member Functions
	G4ParticleHPFissionBaseFS ()
virtual	~G4ParticleHPFissionBaseFS ()
void	Init (G4double A, G4double Z, G4int M, G4String &dirName, G4String &bit, G4ParticleDefinition *)
G4DynamicParticleVector *	ApplyYourself (G4int Prompt)
virtual G4double	GetXsec (G4double anEnergy)
virtual G4ParticleHPVector *	GetXsec ()
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 *projectile)
virtual void	Init (G4double A, G4double Z, G4int M, G4String &dirName, G4String &aFSType, G4ParticleDefinition *)=0
virtual G4HadFinalState *	ApplyYourself (const G4HadProjectile &)
virtual G4ParticleHPFinalState *	New ()=0
G4bool	HasXsec ()
G4bool	HasFSData ()
G4bool	HasAnyData ()
virtual G4double	GetXsec (G4double)
virtual G4ParticleHPVector *	GetXsec ()
void	SetA_Z (G4double anA, G4double aZ, G4int aM=0)
G4double	GetZ ()
G4double	GetN ()
G4double	GetA ()
G4int	GetM ()
void	SetAZMs (G4double anA, G4double aZ, G4int aM, G4ParticleHPDataUsed used)
void	SetProjectile (G4ParticleDefinition *projectile)

Additional Inherited Members
Protected Member Functions inherited from G4ParticleHPFinalState
void	adjust_final_state (G4LorentzVector)
G4bool	DoNotAdjustFinalState ()
Protected Attributes inherited from G4ParticleHPFinalState
G4bool	hasXsec
G4bool	hasFSData
G4bool	hasAnyData
G4ParticleHPNames	theNames
G4Cache< G4HadFinalState * >	theResult
G4ParticleDefinition *	theProjectile
G4double	theBaseA
G4double	theBaseZ
G4int	theBaseM
G4int	theNDLDataZ
G4int	theNDLDataA
G4int	theNDLDataM

Detailed Description

Definition at line 42 of file G4ParticleHPFissionBaseFS.hh.

Constructor & Destructor Documentation

G4ParticleHPFissionBaseFS()

G4ParticleHPFissionBaseFS::G4ParticleHPFissionBaseFS ( )

inline

Definition at line 53 of file G4ParticleHPFissionBaseFS.hh.

  { 
    hasXsec = true; 
    theXsection = new G4ParticleHPVector;
  }

~G4ParticleHPFissionBaseFS()

virtual G4ParticleHPFissionBaseFS::~G4ParticleHPFissionBaseFS ( )

inlinevirtual

Definition at line 59 of file G4ParticleHPFissionBaseFS.hh.

  {
    delete theXsection;
  }

Member Function Documentation

ApplyYourself()

G4DynamicParticleVector * G4ParticleHPFissionBaseFS::ApplyYourself

(

G4int

Prompt

)

Definition at line 96 of file G4ParticleHPFissionBaseFS.cc.

  {  
// if therere were no data for this isotope, break out.    
    if(!HasFSData()) { return 0; }
    
    G4int i;
    G4DynamicParticleVector * aResult = new G4DynamicParticleVector;
    G4ReactionProduct boosted;
    boosted.Lorentz( *(fCache.Get().theNeutronRP) , *(fCache.Get().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;
  }

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

GetXsec() [1/2]

virtual G4ParticleHPVector * G4ParticleHPFissionBaseFS::GetXsec ( )

inlinevirtual

Reimplemented from G4ParticleHPFinalState.

Definition at line 74 of file G4ParticleHPFissionBaseFS.hh.

  {
    return theXsection;
  }

GetXsec() [2/2]

virtual G4double G4ParticleHPFissionBaseFS::GetXsec ( G4double  anEnergy )

inlinevirtual

Reimplemented from G4ParticleHPFinalState.

Definition at line 69 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 *    )

virtual

Implements G4ParticleHPFinalState.

Reimplemented in G4ParticleHPLCFissionFS, G4ParticleHPSCFissionFS, and G4ParticleHPTCFissionFS.

Definition at line 45 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 79 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 85 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