82{
83
85
86
89 if (SDman->FindSensitiveDetector(SensitiveDetectorName="det",0)) {
91 SDman->GetCollectionID(SensitiveDetectorName="det/collection");
92 }
93 }
94
97
98 if (HCE) {
102 }
103 }
104
105
106 G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
107
108
109 G4int eventNumber =
110 G4RunManager::GetRunManager()->GetCurrentEvent()->GetEventID();
111
112
113 if (fSensitiveDetectorHC) {
114 size_t vNumberOfHit = fSensitiveDetectorHC->entries();
115 for (size_t i=0; i<vNumberOfHit; i++) {
117 analysisManager->FillNtupleDColumn(0,0,aHit->
GetEnergy()/CLHEP::keV);
118 analysisManager->FillNtupleDColumn(0,1,aHit->
GetPos().x()/CLHEP::mm);
119 analysisManager->FillNtupleDColumn(0,2,aHit->
GetPos().y()/CLHEP::mm);
120 analysisManager->FillNtupleDColumn(0,3,aHit->
GetPos().z()/CLHEP::mm);
121 analysisManager->FillNtupleDColumn(0,4,aHit->
GetMom().x());
122 analysisManager->FillNtupleDColumn(0,5,aHit->
GetMom().y());
123 analysisManager->FillNtupleDColumn(0,6,aHit->
GetMom().z());
124 analysisManager->FillNtupleDColumn(0,7,aHit->
GetTime()/CLHEP::ns);
125 analysisManager->FillNtupleIColumn(0,8,aHit->
GetType());
126 analysisManager->FillNtupleIColumn(0,9,aHit->
GetTrackID());
127 analysisManager->FillNtupleIColumn(0,10,
fNRi);
128 analysisManager->FillNtupleIColumn(0,11,
fNCi);
129 analysisManager->FillNtupleIColumn(0,12,
fNDi);
130 analysisManager->FillNtupleIColumn(0,13,eventNumber);
131 analysisManager->FillNtupleDColumn(0,14,aHit->
GetWeight());
132 analysisManager->AddNtupleRow(0);
133 }
134 }
135
136}
G4THitsCollection< SAXSSensitiveDetectorHit > SensitiveDetectorHitsCollection
G4ThreeVector GetMom() const
G4ThreeVector GetPos() const
G4double GetEnergy() const
G4double GetWeight() const