11 #include <phparameter/PHParameters.h>
26 #include <Geant4/G4NavigationHistory.hh>
27 #include <Geant4/G4ParticleDefinition.hh>
28 #include <Geant4/G4ReferenceCountedHandle.hh>
29 #include <Geant4/G4Step.hh>
30 #include <Geant4/G4StepPoint.hh>
31 #include <Geant4/G4StepStatus.hh>
32 #include <Geant4/G4String.hh>
33 #include <Geant4/G4SystemOfUnits.hh>
34 #include <Geant4/G4ThreeVector.hh>
35 #include <Geant4/G4TouchableHandle.hh>
36 #include <Geant4/G4Track.hh>
37 #include <Geant4/G4TrackStatus.hh>
38 #include <Geant4/G4Types.hh>
39 #include <Geant4/G4VPhysicalVolume.hh>
40 #include <Geant4/G4VTouchable.hh>
41 #include <Geant4/G4VUserTrackInformation.hh>
55 , m_Detector(detector)
56 , m_Params(parameters)
57 , m_HitContainer(nullptr)
59 , m_SaveShower(nullptr)
60 , m_SaveVolPre(nullptr)
61 , m_SaveVolPost(nullptr)
62 , m_SaveLightYieldFlag(m_Params->get_int_param(
"lightyield"))
64 , m_SavePreStepStatus(-1)
65 , m_SavePostStepStatus(-1)
66 , m_ActiveFlag(m_Params->get_int_param(
"active"))
67 , m_BlackHoleFlag(m_Params->get_int_param(
"blackhole"))
68 , m_UseG4StepsFlag(m_Params->get_int_param(
"use_g4steps"))
69 , m_Tmin(m_Params->get_double_param(
"tmin") *
ns)
70 , m_Tmax(m_Params->get_double_param(
"tmax") *
ns)
102 if ( ! activeMaterial && ! virtualMaterial )
133 bool geantino =
false;
171 std::cout <<
GetName() <<
": New Hit for " << std::endl;
172 std::cout <<
"prestep status: "
174 <<
", poststep status: "
176 <<
", last pre step status: "
178 <<
", last post step status: "
181 <<
", current trackid: " << aTrack->
GetTrackID() << std::endl;
182 std::cout <<
"phys pre vol: " << volume->
GetName()
262 std::cout <<
GetName() <<
": hit was not created" << std::endl;
263 std::cout <<
"prestep status: "
265 <<
", poststep status: "
267 <<
", last pre step status: "
269 <<
", last post step status: "
272 <<
", current trackid: " << aTrack->
GetTrackID() << std::endl;
273 std::cout <<
"phys pre vol: " << volume->
GetName()
284 std::cout <<
GetName() <<
": hits do not belong to the same track" << std::endl;
286 <<
", current trackid: " << aTrack->
GetTrackID()
376 if( virtualMaterial )
410 std::cout <<
"EICG4RPSteppingAction::SetTopNode - unable to find "