46 G4RunManager::GetRunManager()->SetPrintProgress(1);
51 auto analysisManager = G4AnalysisManager::Instance();
56 analysisManager->SetVerboseLevel(1);
57 analysisManager->SetNtupleMerging(
true);
64 analysisManager->CreateH1(
"Eabs" ,
"Edep in absorber", 110, 0., 330*MeV);
65 analysisManager->CreateH1(
"Egap" ,
"Edep in gap", 100, 0., 30*MeV);
66 analysisManager->CreateH1(
"Labs" ,
"trackL in absorber", 100, 0., 50*cm);
67 analysisManager->CreateH1(
"Lgap" ,
"trackL in gap", 100, 0., 50*cm);
71 analysisManager->CreateNtuple(
"B4",
"Edep and TrackL");
72 analysisManager->CreateNtupleDColumn(
"Eabs");
73 analysisManager->CreateNtupleDColumn(
"Egap");
74 analysisManager->CreateNtupleDColumn(
"Labs");
75 analysisManager->CreateNtupleDColumn(
"Lgap");
76 analysisManager->FinishNtuple();
106 auto analysisManager = G4AnalysisManager::Instance();
107 if ( analysisManager->GetH1(1) ) {
108 G4cout << G4endl <<
" ----> print histograms statistic ";
110 G4cout <<
"for the entire run " << G4endl << G4endl;
113 G4cout <<
"for the local thread " << G4endl << G4endl;
116 G4cout <<
" EAbs : mean = "
117 <<
G4BestUnit(analysisManager->GetH1(0)->mean(),
"Energy")
119 <<
G4BestUnit(analysisManager->GetH1(0)->rms(),
"Energy") << G4endl;
121 G4cout <<
" EGap : mean = "
122 <<
G4BestUnit(analysisManager->GetH1(1)->mean(),
"Energy")
124 <<
G4BestUnit(analysisManager->GetH1(1)->rms(),
"Energy") << G4endl;
126 G4cout <<
" LAbs : mean = "
127 <<
G4BestUnit(analysisManager->GetH1(2)->mean(),
"Length")
129 <<
G4BestUnit(analysisManager->GetH1(2)->rms(),
"Length") << G4endl;
131 G4cout <<
" LGap : mean = "
132 <<
G4BestUnit(analysisManager->GetH1(3)->mean(),
"Length")
134 <<
G4BestUnit(analysisManager->GetH1(3)->rms(),
"Length") << G4endl;
139 analysisManager->Write();
140 analysisManager->CloseFile();