FiberSD.cc

Go to the documentation of this file.

#include "SimG4CMS/ShowerLibraryProducer/interface/FiberSD.h"
#include "SimDataFormats/CaloHit/interface/HFShowerPhoton.h"
#include "DataFormats/Math/interface/Point3D.h"
#include "Geometry/HcalCommonData/interface/HcalDDDSimConstants.h"
#include "Geometry/HcalCommonData/interface/HcalSimulationConstants.h"
#include "Geometry/Records/interface/HcalSimNumberingRecord.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"
 
//#define EDM_ML_DEBUG
 
FiberSD::FiberSD(const std::string& iname,
                                               const HcalSimulationConstants* hsps,
                                               const HcalDDDSimConstants* hdc,
                                               const SensitiveDetectorCatalog& clg,
                                               edm::ParameterSet const& p,
                                               const SimTrackManager* manager)
    : SensitiveCaloDetector(iname, clg), theShower(nullptr), theHCID(-1), theHC(nullptr) {
  edm::LogVerbatim("FiberSim") << "FiberSD : Instantiating for " << iname;
  // Get pointer to HcalDDDConstants and HcalSimulationConstants
  theShower = new HFShower(iname, hdc, hsps->hcalsimpar(), p, 1);
}
 
FiberSD::~FiberSD() {
  delete theShower;
  delete theHC;
}
 
void FiberSD::Initialize(G4HCofThisEvent* HCE) {
  edm::LogVerbatim("FiberSim") << "FiberSD : Initialize called for " << GetName() << " in collection "
                               << HCE;
  theHC = new FiberG4HitsCollection(GetName(), collectionName[0]);
  if (theHCID < 0)
    theHCID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
  HCE->AddHitsCollection(theHCID, theHC);
  edm::LogVerbatim("FiberSim") << "FiberSD : Add hit collectrion for " << collectionName[0] << ":"
                               << theHCID << ":" << theHC;
}
 
G4bool FiberSD::ProcessHits(G4Step* aStep, G4TouchableHistory*) {
  //std::vector<HFShower::Hit> hits = theShower->getHits(aStep);
  double zoffset = 1000;
  std::vector<HFShower::Hit> hits = theShower->getHits(aStep, true, zoffset);
 
  if (!hits.empty()) {
    std::vector<HFShowerPhoton> thePE;
    for (unsigned int i = 0; i < hits.size(); i++) {
      //edm::LogVerbatim("FiberSim") << "FiberSD :hit position z " << hits[i].position.z();
      HFShowerPhoton pe = HFShowerPhoton(
          hits[i].position.x(), hits[i].position.y(), hits[i].position.z(), hits[i].wavelength, hits[i].time);
      thePE.push_back(pe);
    }
    int trackID = aStep->GetTrack()->GetTrackID();
    G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
    const G4VTouchable* touch = preStepPoint->GetTouchable();
    G4LogicalVolume* lv = touch->GetVolume(0)->GetLogicalVolume();
    int depth = (touch->GetReplicaNumber(0)) % 10;
    int detID = setDetUnitId(aStep);
    math::XYZPoint theHitPos(
        preStepPoint->GetPosition().x(), preStepPoint->GetPosition().y(), preStepPoint->GetPosition().z());
    //edm::LogVerbatim("FiberSim") << "FiberSD :presteppoint position z " << preStepPoint->GetPosition().z();
 
    FiberG4Hit* aHit = new FiberG4Hit(lv, detID, depth, trackID);
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("FiberSim") << "FiberSD :hit size " << hits.size() << "  npe" << aHit->npe();
    edm::LogVerbatim("FiberSim") << "FiberSD :pre hit position " << aHit->hitPos();
#endif
    aHit->setNpe(hits.size());
    aHit->setPos(theHitPos);
    aHit->setTime(preStepPoint->GetGlobalTime());
    aHit->setPhoton(thePE);
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("FiberSim") << "FiberSD :ShowerPhoton position " << thePE[0].x() << " "
                                 << thePE[0].y() << " " << thePE[0].z();
 
    edm::LogVerbatim("FiberSim") << "FiberSD: Hit created at " << lv->GetName()
                                 << " DetID: " << aHit->towerId() << " Depth: " << aHit->depth()
                                 << " Track ID: " << aHit->trackId() << " Nb. of Cerenkov Photons: " << aHit->npe()
                                 << " Time: " << aHit->time() << " at " << aHit->hitPos();
    for (unsigned int i = 0; i < thePE.size(); i++)
      edm::LogVerbatim("FiberSim") << "FiberSD: PE[" << i << "] " << thePE[i];
#endif
    theHC->insert(aHit);
  }
  return true;
}
 
void FiberSD::EndOfEvent(G4HCofThisEvent* HCE) {
  edm::LogVerbatim("FiberSim") << "FiberSD: finds " << theHC->entries() << " hits";
  clear();
  edm::LogVerbatim("FiberSim") << "theHC entries = " << theHC->entries();
}
 
void FiberSD::clear() {}
 
void FiberSD::DrawAll() {}
 
void FiberSD::PrintAll() {}
 
void FiberSD::update(const BeginOfJob* job) {}
 
void FiberSD::update(const BeginOfRun*) {}
 
void FiberSD::update(const BeginOfEvent*) {}
 
void FiberSD::update(const ::EndOfEvent*) {}
 
void FiberSD::clearHits() {}
 
uint32_t FiberSD::setDetUnitId(const G4Step* aStep) {
  const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
  int fibre = (touch->GetReplicaNumber(1)) % 10;
  int cell = (touch->GetReplicaNumber(2));
  int tower = (touch->GetReplicaNumber(3));
  return ((tower * 1000 + cell) * 10 + fibre);
}
 
void FiberSD::fillHits(edm::PCaloHitContainer&, const std::string&) {}