4 #include "G4VPhysicalVolume.hh" 5 #include "G4PVPlacement.hh" 6 #include "G4HCofThisEvent.hh" 7 #include "G4TouchableHistory.hh" 10 #include "G4VSolid.hh" 11 #include "G4DynamicParticle.hh" 12 #include "G4ParticleDefinition.hh" 13 #include "G4SDManager.hh" 16 #include "G4PhysicalConstants.hh" 17 #include "G4SystemOfUnits.hh" 31 edm::LogVerbatim(
"FiberSim") <<
"HFWedgeSD : Initialize called for " << GetName() <<
" in collection " << HCE;
42 G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
43 const G4VTouchable* touch = preStepPoint->GetTouchable();
45 trackID = aStep->GetTrack()->GetTrackID();
46 edep = aStep->GetTotalEnergyDeposit();
47 time = (preStepPoint->GetGlobalTime()) / ns;
51 const G4DynamicParticle* particle = aStep->GetTrack()->GetDynamicParticle();
52 momDir = particle->GetMomentumDirection();
78 std::map<int, HFShowerG4Hit*>::const_iterator it =
hitMap.find(
currentID);
90 <<
" Edep: " <<
edep / CLHEP::MeV <<
" MeV; Time: " <<
time <<
" ns; Position (local) " 124 const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
125 return (touch->GetReplicaNumber(0));
Log< level::Info, true > LogVerbatim
HFWedgeSD(const std::string &, const SensitiveDetectorCatalog &clg, const SimTrackManager *)
std::vector< PCaloHit > PCaloHitContainer
void updateHit(HFShowerG4Hit *)
void Initialize(G4HCofThisEvent *HCE) override
bool ProcessHits(G4Step *step, G4TouchableHistory *tHistory) override
void setLocalPos(const G4ThreeVector &xyz)
void fillHits(edm::PCaloHitContainer &, const std::string &) override
std::map< int, HFShowerG4Hit * > hitMap
void setTrackId(G4int trackId)
void clearHits() override
void setHitId(G4int hitId)
HFShowerG4HitsCollection * theHC
void EndOfEvent(G4HCofThisEvent *eventHC) override
void setPrimMomDir(const G4ThreeVector &xyz)
uint32_t setDetUnitId(const G4Step *) override
HFShowerG4Hit * currentHit
void setGlobalPos(const G4ThreeVector &xyz)
void updateEnergy(G4double edep)
HFShowerG4Hit * createNewHit()
G4THitsCollection< HFShowerG4Hit > HFShowerG4HitsCollection