46 theFS.
Init(A, Z, M, dirName, aFSType, projectile);
47 theFC.
Init(A, Z, M, dirName, aFSType, projectile);
48 theSC.
Init(A, Z, M, dirName, aFSType, projectile);
49 theTC.
Init(A, Z, M, dirName, aFSType, projectile);
50 theLC.
Init(A, Z, M, dirName, aFSType, projectile);
52 theFF.
Init(A, Z, M, dirName, aFSType, projectile);
55 G4cout <<
"Fission fragment production is now activated in HP package for "
57 <<
", A = " << (
G4int)A
60 G4cout <<
"As currently modeled this option precludes production of delayed neutrons from fission fragments." <<
G4endl;
111 theNeutron.
Lorentz(theNeutron, -1*theTarget);
138 if(random<xSec[i]/xSec[3])
break;
146 G4int Prompt=0, delayed=0, all=0;
152 if(Prompt==0&&delayed==0) Prompt=all;
158 if(Prompt==0&&delayed==0) Prompt=all;
163 if(Prompt==0&&delayed==0) Prompt=all;
168 if(Prompt==0&&delayed==0) Prompt=all;
183 if( theNeutrons != 0)
185 theDecayConstants =
new G4double[delayed];
190 for(i=0; i<theNeutrons->size(); i++)
199 for(i=0; i<theDelayed->size(); i++)
212 theDecayConstants =
new G4double[delayed];
213 if(Prompt==0&&delayed==0) Prompt=all;
219 for(i0=0; i0<Prompt; i0++)
225 for(i0=Prompt; i0<Prompt+delayed; i0++)
229 if ( theDecayConstants[i0-Prompt] > 1.0
e-30 ) {
233 ed <<
" theDecayConstants[i0-Prompt]=" << theDecayConstants[i0-Prompt]
234 <<
" -> cannot sample the time : set it to 0.0 !" <<
G4endl;
244 delete [] theDecayConstants;
246 unsigned int nPhotons = 0;
249 nPhotons = thePhotons->size();
250 for(i=0; i<nPhotons; i++)
297 G4ThreeVector direction (sinth*std::cos(phi) , sinth*std::sin(phi), std::cos(theta) );
void Init(G4double A, G4double Z, G4int M, G4String &dirName, G4String &aFSType, G4ParticleDefinition *projectile)
static G4ParticleHPManager * GetInstance()
G4ParticleHPFCFissionFS theFC
void Init(G4double A, G4double Z, G4int M, G4String &dirName, G4String &aFSType, G4ParticleDefinition *projectile)
std::ostringstream G4ExceptionDescription
G4ParticleHPLCFissionFS theLC
void Init(G4double A, G4double Z, G4int M, G4String &dirName, G4String &aFSType, G4ParticleDefinition *projectile)
static G4IonTable * GetIonTable()
void GetAFissionFragment(G4double, G4int &, G4int &, G4int &)
void SetTarget(const G4ReactionProduct &aTarget)
G4HadSecondary * GetSecondary(size_t i)
void SetDefinition(const G4ParticleDefinition *aParticleDefinition)
void SetMomentum(const G4double x, const G4double y, const G4double z)
void SetNeutronRP(const G4ReactionProduct &aNeutron)
std::vector< G4DynamicParticle * > G4DynamicParticleVector
void Init(G4double A, G4double Z, G4int M, G4String &dirName, G4String &aFSType, G4ParticleDefinition *projectile)
void AddSecondary(G4DynamicParticle *aP, G4int mod=-1)
void Init(G4double A, G4double Z, G4int M, G4String &dirName, G4String &aFSType, G4ParticleDefinition *)
void Init(G4double A, G4double Z, G4int M, G4String &dirName, G4String &aFSType, G4ParticleDefinition *)
G4DynamicParticleVector * ApplyYourself(G4int NNeutrons)
G4double GetPDGMass() const
G4double G4Log(G4double x)
G4DynamicParticleVector * GetPhotons()
void SetKineticEnergy(const G4double en)
G4DynamicParticleVector * ApplyYourself(G4int Prompt, G4int delayed, G4double *decayconst)
G4double GetKineticEnergy() const
G4double GetGlobalTime() const
G4DynamicParticleVector * ApplyYourself(G4int NNeutrons)
const G4Material * GetMaterial() const
G4ParticleHPTCFissionFS theTC
G4ParticleHPFFFissionFS theFF
G4Cache< G4HadFinalState * > theResult
static constexpr double twopi
G4ErrorTarget * theTarget
G4ParticleHPSCFissionFS theSC
virtual G4double GetXsec(G4double anEnergy)
void SetNeutronRP(const G4ReactionProduct &aNeutron)
double A(double temperature)
G4ThreeVector GetMomentum() const
void Put(const value_type &val) const
const G4LorentzVector & Get4Momentum() const
G4int GetNumberOfSecondaries() const
void SetTarget(const G4ReactionProduct &aTarget)
G4double GetKineticEnergy() const
G4double GetFragmentKinetic()
void Init(G4double A, G4double Z, G4int M, G4String &dirName, G4String &aFSType, G4ParticleDefinition *)
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
const G4ParticleDefinition * GetDefinition() const
G4bool produceFissionFragments
G4HadFinalState * ApplyYourself(const G4HadProjectile &theTrack)
G4GLOB_DLL std::ostream G4cout
G4DynamicParticleVector * ApplyYourself(G4int NNeutrons)
static constexpr double pi
void SampleNeutronMult(G4int &all, G4int &Prompt, G4int &delayed, G4double energy, G4int off)
G4bool GetProduceFissionFragments()
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int lvl=0)
void SetLocalEnergyDeposit(G4double aE)
G4ParticleHPFissionERelease * GetEnergyRelease()
void Lorentz(const G4ReactionProduct &p1, const G4ReactionProduct &p2)
G4ReactionProduct GetBiasedThermalNucleus(G4double aMass, G4ThreeVector aVelocity, G4double temp=-1) const
G4DynamicParticleVector * ApplyYourself(G4int nNeutrons)
void SetStatusChange(G4HadFinalStateStatus aS)
G4double GetTemperature() const
G4ParticleHPFSFissionFS theFS