Hi all,
I set water with RINDEX
const int n = 4;
double e[n] = {2.034*eV,2.82*eV,4.136*eV, 4.58*eV};
double rind[n] = {1.3435, 1.353, 1.3608, 1.3610};
fMPT_water->AddProperty("RINDEX",e, rind,n);
and use G4Cerenkov process
Found that it will cause Infinite loop@
G4Cerenkov.cc
at special situation
do
{
rand = G4UniformRand();
NumberOfPhotons = MeanNumberOfPhotons1 -
rand * (MeanNumberOfPhotons1 - MeanNumberOfPhotons2);
N =
G4UniformRand() * std::max(MeanNumberOfPhotons1, MeanNumberOfPhotons2);
// Loop checking, 07-Aug-2015, Vladimir Ivanchenko
} while(N > NumberOfPhotons);
which the MeanNumberOfPhotons1 <0
that is wrong
Geant4 version 11.1.1, 11.1.2, on debian11,
detail value are
dp: 5.77021e-07
nMin: 1.3435
nMax: 1.361
Pmin: 4.00298e-06
PMax: 4.58e-06
CaiMin: 1.07494e-06
CaiMax:1.38692e-06
ge: 3.11984e-07
NumOfPhotons: -0.00127874