35 #include "HistoManager.hh"
70 analysisManager->SetVerboseLevel(2);
71 G4String extension = analysisManager->GetFileType();
75 analysisManager->SetHistoDirectoryName(
"histo");
76 analysisManager->SetNtupleDirectoryName(
"ntuple");
89 analysisManager->SetFirstHistoId(1);
90 analysisManager->SetFirstNtupleId(1);
92 fHistId[1] = analysisManager->CreateH1(
"1",
"Edep in absorber (MeV)",
96 fHistId[2] = analysisManager->CreateH1(
"2",
"Edep in gap (MeV)",
100 fHistId[3] = analysisManager->CreateH1(
"3",
"trackL in absorber (mm)",
104 fHistId[4] = analysisManager->CreateH1(
"4",
"trackL in gap (mm)",
110 analysisManager->CreateNtuple(
"101",
"Edep");
111 fNtColId[0] = analysisManager->CreateNtupleDColumn(
"Eabs");
112 fNtColId[1] = analysisManager->CreateNtupleDColumn(
"Egap");
113 analysisManager->FinishNtuple();
117 analysisManager->CreateNtuple(
"102",
"TrackL");
118 fNtColId[2] = analysisManager->CreateNtupleDColumn(
"Labs");
119 fNtColId[3] = analysisManager->CreateNtupleDColumn(
"Lgap");
120 analysisManager->FinishNtuple();
132 analysisManager->Write();
133 analysisManager->CloseFile();
136 delete G4AnalysisManager::Instance();
146 G4cout <<
"---> warning from HistoManager::FillHisto() : histo " << ih
147 <<
"does note xist; xbin= " << xbin <<
" w= " << weight <<
G4endl;
158 if (ih >= MaxHisto) {
159 G4cout <<
"---> warning from HistoManager::Normalize() : histo " << ih
160 <<
" fac= " << fac <<
G4endl;
174 analysisManager->FillNtupleDColumn(1,
fNtColId[0], energyAbs);
175 analysisManager->FillNtupleDColumn(1,
fNtColId[1], energyGap);
176 analysisManager->AddNtupleRow(1);
178 analysisManager->FillNtupleDColumn(2,
fNtColId[2], trackLAbs);
179 analysisManager->FillNtupleDColumn(2,
fNtColId[3], trackLGap);
180 analysisManager->AddNtupleRow(2);
188 G4cout <<
"\n ----> print histograms statistic \n" <<
G4endl;
void FillHisto(G4int id, G4double e, G4double weight=1.0)
#define G4BestUnit(a, b)
#define G4_USE_G4BESTUNIT_FOR_VERBOSE 1
void Normalize(G4int id, G4double fac)
G4AnaH1 * fHistPt[MaxHisto]
G4GLOB_DLL std::ostream G4cout
ExG4HbookAnalysisManager G4AnalysisManager
static const G4double fac
void FillNtuple(G4double EnergyAbs, G4double EnergyGap, G4double TrackLAbs, G4double TrackLGap)