52{
55 }
58 }
59
61
62
63
64 const G4StepPoint *pPreStepPoint(pStep->GetPreStepPoint());
65 const G4StepPoint *pPostStepPoint(pStep->GetPostStepPoint());
66
67 const double kineticEnergy(pPreStepPoint->GetKineticEnergy());
68
69 const G4ThreeVector preStepPoint(pPreStepPoint->GetPosition());
70 const G4ThreeVector postStepPoint(pPostStepPoint->GetPosition());
71
72
73 const G4String vol_name(pPreStepPoint->GetTouchableHandle()->GetVolume()->GetName());
74
75 const double preStepX(preStepPoint.x());
76 const double preStepY(preStepPoint.y());
77 const double preStepZ(preStepPoint.z());
78
79 const double postStepX(postStepPoint.x());
80 const double postStepY(postStepPoint.y());
81 const double postStepZ(postStepPoint.z());
82
83 int n_plate(0);
84 int n_fiber(0);
85 int n_station(0);
86
88
89
90 const G4Track *pTrack(pStep->GetTrack());
91 const int particleEncoding(pTrack->GetDefinition()->GetPDGEncoding());
92
93 const double globalTime(pTrack->GetGlobalTime());
94 const int trackID(pTrack->GetTrackID());
95
96 std::string vol_test_str = vol_name.substr(0,7);
97
98
99 if(vol_name.find("GVS") != std::string::npos)
100 {
101
102 if(vol_name.find("B7L1") != std::string::npos) n_station=0;
103 else if(vol_name.find("A7L1") != std::string::npos) n_station=1;
104 else if(vol_name.find("A7R1") != std::string::npos) n_station=2;
105 else if(vol_name.find("B7R1") != std::string::npos) n_station=3;
106 else n_station=-1;
108 {
110 (float) energyDeposit,(float) preStepX, (float) preStepY, (float) preStepZ,
111 (float) postStepX, (float) postStepY, (float) postStepZ,(float) globalTime,
112 -1, 100, -1, (int) n_station);
113 }
114 else
115 {
118 G4Exception(
"ALFA_SensitiveDetector",
"InvalidHitColl1", FatalException,
description);
119 return false;
120 }
121 }
122
123 if (vol_test_str.compare("ALFA_Fi") == 0)
124 {
125 if (std::abs(energyDeposit)<std::numeric_limits<double>::epsilon()) { return true; }
127 std::string substring (vol_name);
128 std::string num_string (vol_name);
129
130
131 std::string test_str ("A");
132 test_str = substring.substr(
m_pos2,1);
133 int sign_fiber(0);
134 if (test_str.compare("U") == 0)
135 {
136 sign_fiber = 1;
137 }
138
139 if (test_str.compare("V") == 0)
140 {
141 sign_fiber = -1;
142 }
143
144
145 for (
int k = 0;
k < 3;
k++ )
146 {
147 substring = substring.substr(
m_pos2+1);
151 std::istringstream is(num_string);
153 }
154
155 n_station =
m_num[0];
158
159
161 {
163 trackID,
164 particleEncoding,
165 (float) kineticEnergy,
166 (float) energyDeposit,
167 (float) preStepX, (float) preStepY, (float) preStepZ,
168 (float) postStepX, (float) postStepY, (float) postStepZ,
169 (float) globalTime,
170 (int) sign_fiber, (int) n_plate, (int) n_fiber, (int) n_station
171 );
172 }
173 else
174 {
177 G4Exception(
"ALFA_SensitiveDetector",
"InvalidHitColl2", FatalException,
description);
178 return false;
179 }
180 }
181
182 if (vol_test_str.compare("ODFiber") == 0)
183 {
184 if (std::abs(energyDeposit)<std::numeric_limits<double>::epsilon()) { return true; }
186 std::string substring (vol_name);
187 std::string num_string (std::move(vol_name));
188
189 std::string test_str = substring.substr(
m_pos2,1);
190 int sign_fiber(0);
191 if (test_str.compare("U") == 0)
192 {
193 sign_fiber = 1;
194 }
195
196 if (test_str.compare("V") == 0)
197 {
198 sign_fiber = -1;
199 }
200
201 std::string test_str_side ("A");
202 test_str_side = substring.substr(
m_pos2+1,1);
203
204 int OD_side(0);
205 if (test_str_side.compare("0") == 0)
206 {
207 OD_side = 0;
208 }
209
210 if (test_str_side.compare("1") == 0)
211 {
212 OD_side = 1;
213 }
214
215 for (
int k = 0;
k < 3;
k++ )
216 {
217 substring = substring.substr(
m_pos2+1);
221 std::istringstream is(num_string);
223 }
224
225 n_station =
m_num[0];
228
230 {
232 trackID,
233 particleEncoding,
234 (float) kineticEnergy,
235 (float) energyDeposit,
236 (float) preStepX, (float) preStepY, (float) preStepZ,
237 (float) postStepX, (float) postStepY, (float) postStepZ,
238 (float) globalTime,
239 (int) sign_fiber, (int) OD_side, (int) n_plate, (int) n_fiber, (int) n_station
240 );
241 }
242 else
243 {
246 G4Exception(
"ALFA_SensitiveDetector",
"InvalidHitColl3", FatalException,
description);
247 return false;
248 }
249
250 }
251
252
253 return true;
254}
std::string description
glabal timer - how long have I taken so far?
1.16.1