5 #include <phparameter/PHParameters.h>
19 #include <Geant4/G4ParticleDefinition.hh>
20 #include <Geant4/G4ReferenceCountedHandle.hh>
21 #include <Geant4/G4Step.hh>
22 #include <Geant4/G4StepPoint.hh>
23 #include <Geant4/G4StepStatus.hh>
24 #include <Geant4/G4String.hh>
25 #include <Geant4/G4SystemOfUnits.hh>
26 #include <Geant4/G4ThreeVector.hh>
27 #include <Geant4/G4TouchableHandle.hh>
28 #include <Geant4/G4Track.hh>
29 #include <Geant4/G4TrackStatus.hh>
30 #include <Geant4/G4Types.hh>
31 #include <Geant4/G4VPhysicalVolume.hh>
32 #include <Geant4/G4VTouchable.hh>
33 #include <Geant4/G4VUserTrackInformation.hh>
35 #include <boost/io/ios_state.hpp>
49 , m_Detector(detector)
50 , m_Params(parameters)
51 , m_HitContainer(nullptr)
53 , m_SaveShower(nullptr)
54 , m_SaveVolPre(nullptr)
55 , m_SaveVolPost(nullptr)
56 , m_SaveLightYieldFlag(m_Params->get_int_param(
"lightyield"))
58 , m_SavePreStepStatus(-1)
59 , m_SavePostStepStatus(-1)
60 , m_ActiveFlag(m_Params->get_int_param(
"active"))
61 , m_BlackHoleFlag(m_Params->get_int_param(
"blackhole"))
62 , m_UseG4StepsFlag(m_Params->get_int_param(
"use_g4steps"))
63 , m_Zmin(m_Params->get_double_param(
"place_z") *
cm - m_Params->get_double_param(
"length") *
cm / 2.)
64 , m_Zmax(m_Params->get_double_param(
"place_z") *
cm + m_Params->get_double_param(
"length") *
cm / 2.)
65 , m_Tmin(m_Params->get_double_param(
"tmin") *
ns)
66 , m_Tmax(m_Params->get_double_param(
"tmax") *
ns)
119 bool geantino =
false;
143 cout <<
GetName() <<
": New Hit for " << endl;
149 <<
", current trackid: " << aTrack->
GetTrackID() << endl;
150 cout <<
"phys pre vol: " << volume->
GetName()
201 boost::io::ios_precision_saver ips(cout);
203 <<
"PHG4CylinderStripSteppingAction: Entry hit z " <<
m_Hit->
get_z(0) *
cm
204 <<
" outside acceptance, zmin " <<
m_Zmin
205 <<
", zmax " <<
m_Zmax <<
", layer: " << layer_id << endl;
215 cout <<
GetName() <<
": hit was not created" << endl;
221 <<
", current trackid: " << aTrack->
GetTrackID() << endl;
222 cout <<
"phys pre vol: " << volume->
GetName()
232 cout <<
"hits do not belong to the same track" << endl;
234 <<
", current trackid: " << aTrack->
GetTrackID()
257 <<
" PHG4CylinderStripSteppingAction: Exit hit z " <<
m_Hit->
get_z(1) *
cm
258 <<
" outside acceptance zmin " <<
m_Zmin
259 <<
", zmax " <<
m_Zmax <<
", layer: " << layer_id << endl;
345 m_HitContainer = findNode::getClass<PHG4HitContainer>(topNode, hitnodename.c_str());
350 std::cout <<
"PHG4CylinderStripSteppingAction::SetTopNode - unable to find " << hitnodename << std::endl;