24 #include <phparameter/PHParameters.h>
40 #include <Geant4/G4NavigationHistory.hh>
41 #include <Geant4/G4ParticleDefinition.hh>
42 #include <Geant4/G4ReferenceCountedHandle.hh>
43 #include <Geant4/G4Step.hh>
44 #include <Geant4/G4StepPoint.hh>
45 #include <Geant4/G4StepStatus.hh>
46 #include <Geant4/G4String.hh>
47 #include <Geant4/G4SystemOfUnits.hh>
48 #include <Geant4/G4ThreeVector.hh>
49 #include <Geant4/G4TouchableHandle.hh>
50 #include <Geant4/G4Track.hh>
51 #include <Geant4/G4TrackStatus.hh>
52 #include <Geant4/G4Types.hh>
53 #include <Geant4/G4VPhysicalVolume.hh>
54 #include <Geant4/G4VTouchable.hh>
55 #include <Geant4/G4VUserTrackInformation.hh>
69 , m_Detector(detector)
70 , m_Params(parameters)
71 , m_HitContainer(nullptr)
73 , m_SaveShower(nullptr)
74 , m_SaveVolPre(nullptr)
75 , m_SaveVolPost(nullptr)
76 , m_SaveLightYieldFlag(m_Params->get_int_param(
"lightyield"))
78 , m_SavePreStepStatus(-1)
79 , m_SavePostStepStatus(-1)
80 , m_ActiveFlag(m_Params->get_int_param(
"active"))
81 , m_BlackHoleFlag(m_Params->get_int_param(
"blackhole"))
82 , m_UseG4StepsFlag(m_Params->get_int_param(
"use_g4steps"))
83 , m_Zmin(m_Params->get_double_param(
"place_z") *
cm - m_Params->get_double_param(
"length") *
cm / 2.)
84 , m_Zmax(m_Params->get_double_param(
"place_z") *
cm + m_Params->get_double_param(
"length") *
cm / 2.)
85 , m_Tmin(m_Params->get_double_param(
"tmin") *
ns)
86 , m_Tmax(m_Params->get_double_param(
"tmax") *
ns)
154 bool geantino =
false;
192 std::cout <<
GetName() <<
": New Hit for " << std::endl;
193 std::cout <<
"prestep status: "
195 <<
", poststep status: "
197 <<
", last pre step status: "
199 <<
", last post step status: "
202 <<
", current trackid: " << aTrack->
GetTrackID() << std::endl;
203 std::cout <<
"phys pre vol: " << volume->
GetName()
266 <<
" PHG4CylinderSteppingAction: Entry hit z " <<
m_Hit->
get_z(0) *
cm
267 <<
" outside acceptance, zmin " <<
m_Zmin
268 <<
", zmax " <<
m_Zmax <<
", layer: " << layer_id << std::endl;
280 std::cout <<
GetName() <<
": hit was not created" << std::endl;
281 std::cout <<
"prestep status: "
283 <<
", poststep status: "
285 <<
", last pre step status: "
287 <<
", last post step status: "
290 <<
", current trackid: " << aTrack->
GetTrackID() << std::endl;
291 std::cout <<
"phys pre vol: " << volume->
GetName()
302 std::cout <<
GetName() <<
": hits do not belong to the same track" << std::endl;
304 <<
", current trackid: " << aTrack->
GetTrackID()
332 <<
" PHG4CylinderSteppingAction: Exit hit z " <<
m_Hit->
get_z(1) *
cm
333 <<
" outside acceptance zmin " <<
m_Zmin
334 <<
", zmax " <<
m_Zmax <<
", layer: " << layer_id << std::endl;
422 j_0 = (
int) ((radius - det_x_pos + towersize + (prePoint->
GetPosition().
x() /
cm)) / towersize);
423 k_0 = (
int) ((radius - det_y_pos + towersize + (prePoint->
GetPosition().
y() /
cm)) / towersize);
435 std::cout <<
"EICG4B0ECALSteppingAction::SetTopNode - unable to find "