34 #include "RunAction.hh"
53 analysisManager->OpenFile(fileName);
54 analysisManager->SetVerboseLevel(1);
55 G4String extension = analysisManager->GetFileType();
56 fileName = fileName +
"." + extension;
60 analysisManager->SetFirstHistoId(1);
61 analysisManager->CreateH1(
"1",
"1/(1+(theta+[g]+)**2)",100, 0 ,1.);
62 analysisManager->CreateH1(
"2",
"log10(theta+ [g]+)", 100,-3.,1.);
63 analysisManager->CreateH1(
"3",
"log10(theta- [g]-)", 100,-3.,1.);
64 analysisManager->CreateH1(
"4",
"log10(theta+ [g]+ -theta- [g]-)", 100,-3.,1.);
65 analysisManager->CreateH1(
"5",
"xPlus" ,100,0.,1.);
66 analysisManager->CreateH1(
"6",
"xMinus",100,0.,1.);
68 G4cout <<
"\n----> Histogram file is opened in " << fileName <<
G4endl;
77 analysisManager->Write();
78 analysisManager->CloseFile();
80 delete G4AnalysisManager::Instance();
100 size_t nbProc = fProcCounter->size();
102 while ((i<nbProc)&&((*fProcCounter)[i]->GetName()!=procName)) i++;
103 if (i == nbProc) fProcCounter->push_back(
new OneProcessCount(procName));
105 (*fProcCounter)[i]->Count();
116 G4cout <<
"\n Number of process calls --->";
117 for (
size_t i=0; i< fProcCounter->size();i++) {
118 G4String procName = (*fProcCounter)[i]->GetName();
119 if (procName !=
"Transportation") {
120 G4int count = (*fProcCounter)[i]->GetCounter();
121 G4cout <<
"\t" << procName <<
" : " << count;