5 #include <phparameter/PHParameters.h>
20 #include <Geant4/G4ParticleDefinition.hh>
21 #include <Geant4/G4ReferenceCountedHandle.hh>
22 #include <Geant4/G4Step.hh>
23 #include <Geant4/G4StepPoint.hh>
24 #include <Geant4/G4StepStatus.hh>
25 #include <Geant4/G4String.hh>
26 #include <Geant4/G4SystemOfUnits.hh>
27 #include <Geant4/G4ThreeVector.hh>
28 #include <Geant4/G4TouchableHandle.hh>
29 #include <Geant4/G4Track.hh>
30 #include <Geant4/G4TrackStatus.hh>
31 #include <Geant4/G4Types.hh>
32 #include <Geant4/G4VPhysicalVolume.hh>
33 #include <Geant4/G4VTouchable.hh>
34 #include <Geant4/G4VUserTrackInformation.hh>
47 , m_Detector(detector)
48 , m_Params(parameters)
49 , m_HitContainer(nullptr)
51 , m_SaveHitContainer(nullptr)
52 , m_SaveVolPre(nullptr)
53 , m_SaveVolPost(nullptr)
55 , m_SavePreStepStatus(-1)
56 , m_SavePostStepStatus(-1)
57 , m_ActiveFlag(m_Params->get_int_param(
"active"))
58 , m_BlackHoleFlag(m_Params->get_int_param(
"blackhole"))
104 int detector_id = -1;
105 if (volume->
GetName().find(
"Coil") != string::npos)
107 for (
int i = 0; i <= 4; i++)
116 else if (volume->
GetName().find(
"CryostatHe") != string::npos)
120 else if (volume->
GetName().find(
"CryostatAl") != string::npos)
124 else if (volume->
GetName().find(
"Yoke") != string::npos)
130 cout <<
"cannot extract detector id from " << volume->
GetName() << endl;
132 bool geantino =
false;
170 cout <<
GetName() <<
": New Hit for " << endl;
171 cout <<
"prestep status: "
173 <<
", poststep status: "
175 <<
", last pre step status: "
177 <<
", last post step status: "
180 <<
", current trackid: " << aTrack->
GetTrackID() << endl;
181 cout <<
"phys pre vol: " << volume->
GetName()
216 cout <<
"implement stuff for whichactive < 0 (inactive volumes)" << endl;
238 cout <<
GetName() <<
": hit was not created" << endl;
239 cout <<
"prestep status: "
241 <<
", poststep status: "
243 <<
", last pre step status: "
245 <<
", last post step status: "
248 <<
", current trackid: " << aTrack->
GetTrackID() << endl;
249 cout <<
"phys pre vol: " << volume->
GetName()
259 cout <<
GetName() <<
": hits do not belong to the same track" << endl;
261 <<
", current trackid: " << aTrack->
GetTrackID()
349 m_HitContainer = findNode::getClass<PHG4HitContainer>(topNode, hitnodename);
353 std::cout <<
"BeastMagnetSteppingAction::SetTopNode - unable to find "
354 << hitnodename << std::endl;