Hi~
When I run my case, I got the following messages:
G4WT5 > 080808 Something unexpected is happen in G4ParticleHPLabAngularEnergy
G4WT5 > 080808 Something unexpected is happen in G4ParticleHPLabAngularEnergy
G4WT7 > 080808 Something unexpected is happen in G4ParticleHPLabAngularEnergy
G4WT7 > 080808 Something unexpected is happen in G4ParticleHPLabAngularEnergy
G4WT6 > 080808 Something unexpected is happen in G4ParticleHPLabAngularEnergy
G4WT6 > 080808 Something unexpected is happen in G4ParticleHPLabAngularEnergy
G4WT12 > 080808 Something unexpected is happen in G4ParticleHPLabAngularEnergy
G4WT12 > 080808 Something unexpected is happen in G4ParticleHPLabAngularEnergy
G4WT14 > 080808 Something unexpected is happen in G4ParticleHPLabAngularEnergy
G4WT14 > 080808 Something unexpected is happen in G4ParticleHPLabAngularEnergy
G4WT3 > 080808 Something unexpected is happen in G4ParticleHPLabAngularEnergy
…
Then I found out this:
112 // get theta, E
113 G4double cosTh, secEnergy;
114 G4int i, it(0);
115 // find the energy bin
116 for(i=0; i<nEnergies; i++)
117 {
118 it = i;
119 if ( anEnergy < theEnergies[i] ) break;
120 }
121 //080808
122 //if ( it == 0 || it == nEnergies-1 ) // it marks the energy bin
123 if ( it == 0 ) // it marks the energy bin
124 {
125 G4cout << "080808 Something unexpected is happen in G4ParticleHPLabAngularEnergy " << G4endl;
126 // integrate the prob for each costh, and select theta.
127 G4double * running = new G4double [nCosTh[it]];
128 running[0]=0;
129 for(i=0;i<nCosTh[it]; i++)
130 {
131 if(i!=0) running[i] = running[i-1];
132 running[i]+=theData[it][i].GetIntegral(); // Does interpolated integral.
133 }
134 G4double random = running[nCosTh[it]-1]*G4UniformRand();
135 G4int ith(0);
136 for(i=0;i<nCosTh[it]; i++)
137 {
138 ith = i;
139 if(random<running[i]) break;
140 }
from this website: http://hurel.hanyang.ac.kr/Geant4/Doxygen/10.03.p01/html/class_g4_particle_h_p_lab_angular_energy.html
But I still cannot fingure out the cause of those messages and what effect they would make to my results.
Can anyone help me solve this issue? Thank you.
Best regards,
jcren.