49 G4int midBin = (je1 + je2)/2;
50 if (ek <
elab[midBin])
54 }
while ( (je2 - je1) > 1 && --iterationsLeft > 0 );
55 if ( iterationsLeft <= 0 ) {
57 "Problem with energy bin (elab) data");
71 G4double b = sig[je1][0] - rc*elab[je1];
75 if (verboseLevel > 1)
G4cout <<
"sample=" << sample <<
G4endl
76 << ek <<
" " << ke1 <<
" " << ke2 <<
" "
77 << sigint1 <<
" " << sigint2 <<
G4endl;
79 iterationsLeft= 2*
NANGLE +1;
81 G4int midBin = (ke1 + ke2)/2;
82 dsig = sig[je2][midBin] - sig[je1][midBin];
84 b = sig[je1][midBin] - rc*elab[je1];
86 if (sample < sigint) {
94 if (verboseLevel > 1)
G4cout << ke1 <<
" " << ke2 <<
" "
95 << sigint1 <<
" " << sigint2 <<
G4endl;
96 }
while ( (ke2 - ke1) > 1 && --iterationsLeft > 0);
97 if ( iterationsLeft <= 0 ) {
99 "Problem with angular distribution (sig) data");
114 dsig = sigint2 - sigint1;
116 b = ke1 - rc*sigint1;
124 if (verboseLevel > 1) {
125 G4cout <<
" energy bin " << je1 <<
" energy=" << elab[je1] <<
G4endl;
static const G4float sig[NENERGY][NANGLE]
static const G4float elab[NENERGY]
static constexpr double twopi
virtual G4double Phi() const
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
static constexpr double degree
G4GLOB_DLL std::ostream G4cout
virtual G4double CosTheta(G4double s, G4double m1, G4double m2) const
static constexpr double pi
static constexpr double GeV