98 for (
G4int i=1; i<nbins; ++i) {
103 if(y < sum*0.01) {
break; }
117 }
else if(
OPTxs <= 2) {
147 for(i=0; i<100; ++i) {
148 T = theMinKinEnergy + delta*rndm->
flat();
G4double theEmissionProbability
G4double SampleKineticEnergy(const G4Fragment &aFragment)
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4DeexPrecoParameters * theParameters
virtual ~G4PreCompoundFragment()
static G4double ComputeCrossSection(G4double K, G4double cb, G4double resA13, G4double amu1, G4int idx, G4int Z, G4int A, G4int resA)
virtual G4double GetAlpha() const =0
virtual G4double ProbabilityDistributionFunction(G4double K, const G4Fragment &aFragment)=0
static G4double ComputeCrossSection(G4double K, G4double cb, G4double resA13, G4double amu1, G4int idx, G4int Z, G4int resA)
G4double CrossSection(G4double ekin) const
G4double theCoulombBarrier
static constexpr double MeV
G4double IntegrateEmissionProbability(G4double Low, G4double Up, const G4Fragment &aFragment)
G4PreCompoundFragment(const G4ParticleDefinition *, G4VCoulombBarrier *aCoulombBarrier)
virtual G4double GetBeta() const =0
G4double CalcEmissionProbability(const G4Fragment &aFragment)
static G4double ComputePowerParameter(G4int resA, G4int idx)
G4double GetOpt0(G4double ekin) const
static constexpr double pi