} for(int i=0;iE2) { G4cout << \above \ return; } G4int index = static_cast ((E2-Energy)/dE); index = NCHN-1-index; if(index<0) index = 0; TData[index] += Data; } //Record tally data(single-channel). void MyDataRecorder::AddData(G4double Data) { TData[0] += Data; } //Calculating relative errors. void MyDataRecorder::CalErr() { for(G4int i=0;i1 && MyData[i]>0) { rerr[i]=sqrt( sumD2[i]- MyData[i]*MyData[i]/NPS )/MyData[i]; } } } //Updating data after each event, it should be called at the end of each event. void MyDataRecorder::UpdateThisEvent() { for(int i=0; i 1) { if(total>0) totalerr = sqrt( totalSD2 - total*total/NPS )/total; dataFile << std::scientific << std::setprecision(4) << \ << std::setprecision(5) << total/NPS/factor << \ << std::fixed << std::setprecision(4) << totalerr << G4endl; } //cnscott的博客 cnscott.blog.163.com dataFile << \ dataFile << \ dataFile << \ if(!message.isNull()) dataFile << message << G4endl; dataFile.close(); G4cout<< \\\\<< FileName << \is created instead.\<
