61 {
62
64
65
66 G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
67
69 fGain = garfieldPhysics->
GetGain();
70
71
72 analysisManager->FillH1(1, fEnergyAbs);
73 analysisManager->FillH1(2, fTrackLAbs);
74 analysisManager->FillH1(3, fEnergyGas);
75 analysisManager->FillH1(4, fAvalancheSize);
76 analysisManager->FillH1(5, fGain);
77
78
79 analysisManager->FillNtupleDColumn(0, fEnergyAbs);
80 analysisManager->FillNtupleDColumn(1, fTrackLAbs);
81 analysisManager->FillNtupleDColumn(2, fEnergyGas);
82 analysisManager->FillNtupleDColumn(3, fAvalancheSize);
83 analysisManager->FillNtupleDColumn(4, fGain);
84
85
86
87 G4int eventID = event->GetEventID();
88 G4int printModulo = G4RunManager::GetRunManager()->GetPrintProgress();
89 if ((printModulo > 0) && (eventID % printModulo == 0)) {
90 G4cout << "---> End of event: " << eventID << G4endl;
91
92 G4cout << " Absorber: total energy: " << std::setw(7)
93 << G4BestUnit(fEnergyAbs, "Energy")
94 << " total track length: " << std::setw(7)
95 << G4BestUnit(fTrackLAbs, "Length") << G4endl;
96
97 G4cout << " Gas: total energy: " << std::setw(7)
98 << G4BestUnit(fEnergyGas, "Energy")
99 << " avalanche size: " << fAvalancheSize
100 << " gain: " << fGain << G4endl;
101 }
102}
double GetAvalancheSize() const