47 G4RunManager::GetRunManager()->SetPrintProgress(1);
52 auto analysisManager = G4AnalysisManager::Instance();
57 analysisManager->SetVerboseLevel(1);
58 analysisManager->SetNtupleMerging(
true);
65 analysisManager->CreateH1(
"Eabs" ,
"Edep in absorber", 110, 0., 330*MeV);
66 analysisManager->CreateH1(
"Egap" ,
"Edep in gap", 100, 0., 30*MeV);
67 analysisManager->CreateH1(
"Labs" ,
"trackL in absorber", 100, 0., 50*cm);
68 analysisManager->CreateH1(
"Lgap" ,
"trackL in gap", 100, 0., 50*cm);
72 analysisManager->CreateNtuple(
"B4",
"Edep and TrackL");
73 analysisManager->CreateNtupleDColumn(
"Eabs");
74 analysisManager->CreateNtupleDColumn(
"Egap");
75 analysisManager->CreateNtupleDColumn(
"Labs");
76 analysisManager->CreateNtupleDColumn(
"Lgap");
77 analysisManager->FinishNtuple();
116 auto analysisManager = G4AnalysisManager::Instance();
117 if ( analysisManager->GetH1(1) ) {
118 G4cout << G4endl <<
" ----> print histograms statistic ";
120 G4cout <<
"for the entire run " << G4endl << G4endl;
123 G4cout <<
"for the local thread " << G4endl << G4endl;
126 G4cout <<
" EAbs : mean = "
127 <<
G4BestUnit(analysisManager->GetH1(0)->mean(),
"Energy")
129 <<
G4BestUnit(analysisManager->GetH1(0)->rms(),
"Energy") << G4endl;
131 G4cout <<
" EGap : mean = "
132 <<
G4BestUnit(analysisManager->GetH1(1)->mean(),
"Energy")
134 <<
G4BestUnit(analysisManager->GetH1(1)->rms(),
"Energy") << G4endl;
136 G4cout <<
" LAbs : mean = "
137 <<
G4BestUnit(analysisManager->GetH1(2)->mean(),
"Length")
139 <<
G4BestUnit(analysisManager->GetH1(2)->rms(),
"Length") << G4endl;
141 G4cout <<
" LGap : mean = "
142 <<
G4BestUnit(analysisManager->GetH1(3)->mean(),
"Length")
144 <<
G4BestUnit(analysisManager->GetH1(3)->rms(),
"Length") << G4endl;
149 analysisManager->Write();
150 analysisManager->CloseFile();