22 #include <Geant4/G4IonisParamMat.hh>
23 #include <Geant4/G4Material.hh>
24 #include <Geant4/G4MaterialCutsCouple.hh>
25 #include <Geant4/G4ParticleDefinition.hh>
26 #include <Geant4/G4Step.hh>
27 #include <Geant4/G4StepPoint.hh>
28 #include <Geant4/G4StepStatus.hh>
29 #include <Geant4/G4String.hh>
30 #include <Geant4/G4SystemOfUnits.hh>
31 #include <Geant4/G4ThreeVector.hh>
32 #include <Geant4/G4TouchableHandle.hh>
33 #include <Geant4/G4Track.hh>
34 #include <Geant4/G4TrackStatus.hh>
35 #include <Geant4/G4Types.hh>
36 #include <Geant4/G4VTouchable.hh>
37 #include <Geant4/G4VUserTrackInformation.hh>
52 absorberhits_(nullptr),
54 savehitcontainer(nullptr),
57 savepoststepstatus(-1)
90 bool geantino =
false;
212 cout <<
"PHG4SpacalSteppingAction: hit outside acceptance, layer: "
224 cout <<
GetName() <<
": hit was not created" << endl;
233 cout <<
GetName() <<
": hits do not belong to the same track" << endl;
235 <<
", current trackid: " << aTrack->
GetTrackID()
259 static bool once =
true;
260 if (once and edep > 0)
266 cout <<
"PHG4SpacalSteppingAction::UserSteppingAction::"
269 <<
" use scintillating light model at each Geant4 steps. "
274 <<
"Birk Constant = "
277 <<
"edep = " << edep <<
", "
280 <<
"light_yield = " << light_yield << endl;
289 cout <<
"PHG4SpacalSteppingAction: hit outside acceptance get_zmin() "
353 string absorbernodename;
366 hits_ = findNode::getClass<PHG4HitContainer>(topNode, hitnodename.c_str());
367 absorberhits_ = findNode::getClass<PHG4HitContainer>(topNode,
368 absorbernodename.c_str());
372 std::cout <<
"PHG4SpacalSteppingAction::SetTopNode - unable to find "
373 << hitnodename << std::endl;
379 std::cout <<
"PHG4SpacalSteppingAction::SetTopNode - unable to find "
380 << absorbernodename << std::endl;