5 #include <phparameter/PHParameters.h>
16 #include <Geant4/G4ParticleDefinition.hh>
17 #include <Geant4/G4ReferenceCountedHandle.hh>
18 #include <Geant4/G4Step.hh>
19 #include <Geant4/G4StepPoint.hh>
20 #include <Geant4/G4StepStatus.hh>
21 #include <Geant4/G4SystemOfUnits.hh>
22 #include <Geant4/G4ThreeVector.hh>
23 #include <Geant4/G4TouchableHandle.hh>
24 #include <Geant4/G4Track.hh>
25 #include <Geant4/G4TrackStatus.hh>
26 #include <Geant4/G4Types.hh>
27 #include <Geant4/G4VPhysicalVolume.hh>
28 #include <Geant4/G4VTouchable.hh>
29 #include <Geant4/G4VUserTrackInformation.hh>
31 #include <boost/io/ios_state.hpp>
43 , m_Detector(detector)
44 , m_Params(parameters)
45 , m_HitContainer(nullptr)
47 , m_ActiveGasHits(nullptr)
49 , m_SaveHitContainer(nullptr)
50 , m_SaveShower(nullptr)
51 , m_SaveVolPre(nullptr)
52 , m_SaveVolPost(nullptr)
54 , m_SavePreStepStatus(-1)
55 , m_SavePostStepStatus(-1)
56 , m_BlackHoleFlag(m_Params->get_int_param(
"blackhole"))
57 , m_ActiveFlag(m_Params->get_int_param(
"active"))
58 , m_UseG4StepsFlag(m_Params->get_int_param(
"use_g4steps"))
59 , m_Zmin(m_Params->get_double_param(
"z_begin") *
cm)
60 , m_Zmax(m_Params->get_double_param(
"z_begin") *
cm + 3.0)
102 bool geantino =
false;
178 std::cout <<
"hits do not belong to the same track" << std::endl;
180 <<
", current trackid: " << aTrack->
GetTrackID()
266 std::string hitnodename;
268 std::string gasnodename;
282 m_HitContainer = findNode::getClass<PHG4HitContainer>(topNode, hitnodename.c_str());
283 m_ActiveGasHits = findNode::getClass<PHG4HitContainer>(topNode, gasnodename.c_str());
288 std::cout <<
"PHG4ECAPToFSteppingAction::SetTopNode - unable to find " << hitnodename << std::endl;
291 if (!m_ActiveGasHits)
295 std::cout <<
"PHG4ECAPToFSteppingAction::SetTopNode - unable to find " << gasnodename << std::endl;