HFWedgeSD.cc

Go to the documentation of this file.

#include "SimG4CMS/ShowerLibraryProducer/interface/HFWedgeSD.h"
#include "DataFormats/Math/interface/Point3D.h"
 
#include "G4VPhysicalVolume.hh"
#include "G4PVPlacement.hh"
#include "G4HCofThisEvent.hh"
#include "G4TouchableHistory.hh"
#include "G4Track.hh"
#include "G4Step.hh"
#include "G4VSolid.hh"
#include "G4DynamicParticle.hh"
#include "G4ParticleDefinition.hh"
#include "G4SDManager.hh"
#include "G4ios.hh"
 
#include "G4PhysicalConstants.hh"
#include "G4SystemOfUnits.hh"
 
HFWedgeSD::HFWedgeSD(const std::string& iname,
                     const edm::EventSetup& es,
                     const SensitiveDetectorCatalog& clg,
                     edm::ParameterSet const& p,
                     const SimTrackManager* manager)
    : SensitiveCaloDetector(iname, es, clg, p),
      m_trackManager(manager),
      hcID(-1),
      theHC(nullptr),
      currentHit(nullptr) {}
 
HFWedgeSD::~HFWedgeSD() { delete theHC; }
 
void HFWedgeSD::Initialize(G4HCofThisEvent* HCE) {
  LogDebug("FiberSim") << "HFWedgeSD : Initialize called for " << GetName();
  theHC = new HFShowerG4HitsCollection(GetName(), collectionName[0]);
  if (hcID < 0)
    hcID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
  HCE->AddHitsCollection(hcID, theHC);
 
  clearHits();
}
 
G4bool HFWedgeSD::ProcessHits(G4Step* aStep, G4TouchableHistory*) {
  G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
  const G4VTouchable* touch = preStepPoint->GetTouchable();
  currentID = setDetUnitId(aStep);
  trackID = aStep->GetTrack()->GetTrackID();
  edep = aStep->GetTotalEnergyDeposit();
  time = (preStepPoint->GetGlobalTime()) / ns;
 
  globalPos = preStepPoint->GetPosition();
  localPos = touch->GetHistory()->GetTopTransform().TransformPoint(globalPos);
  const G4DynamicParticle* particle = aStep->GetTrack()->GetDynamicParticle();
  momDir = particle->GetMomentumDirection();
 
  if (hitExists() == false && edep > 0.)
    currentHit = createNewHit();
 
  return true;
}
 
void HFWedgeSD::EndOfEvent(G4HCofThisEvent* HCE) {
  LogDebug("FiberSim") << "HFWedgeSD: Sees" << theHC->entries() << " hits";
  clear();
}
 
void HFWedgeSD::clear() {}
 
void HFWedgeSD::DrawAll() {}
 
void HFWedgeSD::PrintAll() {}
 
G4bool HFWedgeSD::hitExists() {
  // Update if in the same detector, time-slice and for same track
  if (currentID == previousID) {
    updateHit(currentHit);
    return true;
  }
 
  std::map<int, HFShowerG4Hit*>::const_iterator it = hitMap.find(currentID);
  if (it != hitMap.end()) {
    updateHit(currentHit);
    return true;
  }
 
  return false;
}
 
HFShowerG4Hit* HFWedgeSD::createNewHit() {
  LogDebug("FiberSim") << "HFWedgeSD::CreateNewHit for ID " << currentID << " Track " << trackID
                       << " Edep: " << edep / MeV << " MeV; Time: " << time << " ns; Position (local) " << localPos
                       << " (global ) " << globalPos << " direction " << momDir;
  HFShowerG4Hit* aHit = new HFShowerG4Hit;
  aHit->setHitId(currentID);
  aHit->setTrackId(trackID);
  aHit->setTime(time);
  aHit->setLocalPos(localPos);
  aHit->setGlobalPos(globalPos);
  aHit->setPrimMomDir(momDir);
  updateHit(aHit);
 
  theHC->insert(aHit);
  hitMap.insert(std::pair<int, HFShowerG4Hit*>(previousID, aHit));
 
  return aHit;
}
 
void HFWedgeSD::updateHit(HFShowerG4Hit* aHit) {
  if (edep != 0) {
    aHit->updateEnergy(edep);
    LogDebug("FiberSim") << "HFWedgeSD: Add energy deposit in " << currentID << " edep " << edep / MeV << " MeV";
  }
  previousID = currentID;
}
 
void HFWedgeSD::clearHits() {
  hitMap.erase(hitMap.begin(), hitMap.end());
  previousID = -1;
}
 
uint32_t HFWedgeSD::setDetUnitId(const G4Step* aStep) {
  const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
  return (touch->GetReplicaNumber(0));
}
 
void HFWedgeSD::fillHits(edm::PCaloHitContainer&, const std::string&) {}