91 G4int nbins = del*den;
99 for (
G4int i=1; i<nbins; ++i) {
105 if(y < sum*0.01) {
break; }
119 }
else if(
OPTxs <= 2) {
136 G4double r0 = theParameters->GetR0()*theResA13;
139 return 1.e+25*
CLHEP::pi*r0*r0*theResA13*GetAlpha()*(1.0 + GetBeta()/ekin);
150 for(i=0; i<100; ++i) {
static G4double ComputeCrossSection(G4double K, G4double cb, G4double resA13, G4double amu1, G4int idx, G4int Z, G4int resA)
static G4double ComputePowerParameter(G4int resA, G4int idx)
static G4double ComputeCrossSection(G4double K, G4double cb, G4double resA13, G4double amu1, G4int idx, G4int Z, G4int A, G4int resA)
G4double CalcEmissionProbability(const G4Fragment &aFragment) override
G4double GetOpt0(G4double ekin) const
G4double IntegrateEmissionProbability(G4double low, G4double up, const G4Fragment &aFragment)
G4PreCompoundFragment(const G4ParticleDefinition *, G4VCoulombBarrier *aCoulombBarrier)
G4double CrossSection(G4double ekin) const
virtual ~G4PreCompoundFragment()
G4double SampleKineticEnergy(const G4Fragment &aFragment) override
virtual G4double ProbabilityDistributionFunction(G4double ekin, const G4Fragment &aFragment)=0
G4double theEmissionProbability
G4double theCoulombBarrier
static constexpr double MeV
static constexpr double pi
T max(const T t1, const T t2)
brief Return the largest of the two arguments