| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 /// \file XraySPOSteppingAction.cc
27 /// \brief Implementation of the SteppingActio
28 //
29 // Authors: P.Dondero (paolo.dondero@cern.ch),
30 //
31 //....oooOO0OOooo........oooOO0OOooo........oo
32 //....oooOO0OOooo........oooOO0OOooo........oo
33
34
35 #include "XraySPOSteppingAction.hh"
36 #include "XraySPODetectorConstruction.hh"
37 #include "XraySPOHistoManager.hh"
38 #include "G4SteppingManager.hh"
39 #include "G4RunManager.hh"
40 #include "Randomize.hh"
41
42 //....oooOO0OOooo........oooOO0OOooo........oo
43
44 void XraySPOSteppingAction::UserSteppingAction
45 {
46 G4int parentID = step->GetTrack()->GetParent
47
48 // Only primaries
49 if (parentID == 0)
50 {
51 G4RunManager *rm = G4RunManager::GetRunMan
52 G4int eventID = rm->GetCurrentEvent()->Get
53
54 G4double len_fraction = 0.1; // approx 1/2
55 G4bool proceed = false;
56
57 if (fPrevEventID != eventID)
58 {
59 fNumReflections = -2;
60 }
61
62 G4ThreeVector direction = step->GetTrack()
63 auto theta = (G4double)direction.theta();
64 auto phi = (G4double)direction.phi();
65
66 G4double x = step->GetPreStepPoint()->GetP
67 G4double y = step->GetPreStepPoint()->GetP
68 G4double z = step->GetPreStepPoint()->GetP
69
70 const G4String& vol_name = step->GetTrack(
71 G4int trackid = step->GetTrack()->GetTrack
72
73 G4StepPoint * postStep = step->GetPostStep
74 G4StepPoint* post_point = step->GetPostSte
75
76 const G4String& proc_name = postStep->GetP
77
78 // Count the number of reflections
79 G4String s = "";
80 G4String& post_phys_name = s;
81 if (fPrevEventID == eventID && step->GetSt
82 {
83 if (post_point != nullptr)
84 {
85 G4VPhysicalVolume* post_phys = post_po
86 if (post_phys != nullptr)
87 {
88 post_phys_name = post_phys->GetName(
89 }
90
91 if (post_phys_name == "cone_1" || post
92 {
93 if (fNumReflections < 0)
94 {
95 fNumReflections = 1;
96 }
97 else
98 {
99 fNumReflections += 1;
100 }
101 proceed = true;
102 }
103 }
104 }
105
106 if (vol_name == "pDummyEntrance" || vol_na
107 {
108 proceed = true;
109 }
110
111 if (proceed)
112 {
113 // G4cout << "Save with " << fNumReflect
114 G4AnalysisManager* analysisManager = G4A
115 analysisManager->FillNtupleIColumn(0, ev
116 analysisManager->FillNtupleSColumn(1, vo
117 analysisManager->FillNtupleIColumn(2, tr
118 analysisManager->FillNtupleDColumn(3, x)
119 analysisManager->FillNtupleDColumn(4, y)
120 analysisManager->FillNtupleDColumn(5, z)
121 analysisManager->FillNtupleDColumn(6, th
122 analysisManager->FillNtupleDColumn(7, ph
123 analysisManager->FillNtupleSColumn(8, pr
124 analysisManager->FillNtupleIColumn(9, pa
125 analysisManager->FillNtupleIColumn(10, f
126 analysisManager->AddNtupleRow();
127 }
128 fPrevx = x;
129 fPrevy = y;
130 fPrevz = z;
131 fPrevEventID = eventID;
132 }
133 }
134
135 //....oooOO0OOooo........oooOO0OOooo........oo
136