#include <FiberSD.h>

Inheritance diagram for FiberSD:

Public Member Functions
void	clear () override

void	clearHits () override

void	DrawAll () override

void	EndOfEvent (G4HCofThisEvent *HCE) override

	FiberSD (const std::string &, const HcalSimulationConstants , const HcalDDDSimConstants , const SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager *)

void	fillHits (edm::PCaloHitContainer &, const std::string &) override

void	Initialize (G4HCofThisEvent *HCE) override

void	PrintAll () override

G4bool	ProcessHits (G4Step aStep, G4TouchableHistory ROhist) override

uint32_t	setDetUnitId (const G4Step *) override

	~FiberSD () override

Public Member Functions inherited from SensitiveCaloDetector
virtual void	reset ()

	SensitiveCaloDetector (const std::string &iname, const SensitiveDetectorCatalog &clg, const std::string &newcollname="")

Public Member Functions inherited from SensitiveDetector
void	EndOfEvent (G4HCofThisEvent *eventHC) override

const std::vector< std::string > &	getNames () const

void	Initialize (G4HCofThisEvent *eventHC) override

bool	isCaloSD () const

	SensitiveDetector (const std::string &iname, const SensitiveDetectorCatalog &, bool calo, const std::string &newcollname="")

	~SensitiveDetector () override

Public Member Functions inherited from Observer< const BeginOfJob *>
	Observer ()

void	slotForUpdate (const BeginOfJob * iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const BeginOfRun *>
	Observer ()

void	slotForUpdate (const BeginOfRun * iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const BeginOfEvent *>
	Observer ()

void	slotForUpdate (const BeginOfEvent * iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const EndOfEvent *>
	Observer ()

void	slotForUpdate (const EndOfEvent * iT)

virtual	~Observer ()

Protected Member Functions
void	update (const BeginOfJob *) override
	This routine will be called when the appropriate signal arrives. More...

void	update (const BeginOfRun *) override
	This routine will be called when the appropriate signal arrives. More...

void	update (const BeginOfEvent *) override
	This routine will be called when the appropriate signal arrives. More...

void	update (const ::EndOfEvent *) override

Protected Member Functions inherited from SensitiveDetector
TrackInformation *	cmsTrackInformation (const G4Track *aTrack)

Local3DPoint	ConvertToLocal3DPoint (const G4ThreeVector &point) const

Local3DPoint	FinalStepPosition (const G4Step *step, coordinates) const

Local3DPoint	InitialStepPosition (const G4Step *step, coordinates) const

Local3DPoint	LocalPostStepPosition (const G4Step *step) const

Local3DPoint	LocalPreStepPosition (const G4Step *step) const

void	NaNTrap (const G4Step *step) const

void	setNames (const std::vector< std::string > &)

Protected Member Functions inherited from Observer< const EndOfEvent *>
virtual void	update (const EndOfEvent *)=0
	This routine will be called when the appropriate signal arrives. More...

Private Attributes
FiberG4HitsCollection *	theHC

G4int	theHCID

HFShower *	theShower

Additional Inherited Members
Protected Types inherited from SensitiveDetector
enum	coordinates { WorldCoordinates, LocalCoordinates }

Detailed Description

Definition at line 29 of file FiberSD.h.

Constructor & Destructor Documentation

◆ FiberSD()

FiberSD::FiberSD	(	const std::string &	iname,
		const HcalSimulationConstants *	hsps,
		const HcalDDDSimConstants *	hdc,
		const SensitiveDetectorCatalog &	clg,
		edm::ParameterSet const &	p,
		const SimTrackManager *	manager
	)

explicit

Definition at line 22 of file FiberSD.cc.

References HcalSimulationConstants::hcalsimpar(), Calorimetry_cff::HFShower, AlCaHLTBitMon_ParallelJobs::p, and theShower.

     : 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::~FiberSD ( )

override

Definition at line 34 of file FiberSD.cc.

References theHC, and theShower.

                   {
   delete theShower;
   delete theHC;
 }

Member Function Documentation

◆ clear()

void FiberSD::clear ( void )

override

Definition at line 103 of file FiberSD.cc.

Referenced by EndOfEvent().

103 {}

◆ clearHits()

void FiberSD::clearHits ( )

overridevirtual

Implements SensitiveDetector.

Definition at line 117 of file FiberSD.cc.

117 {}

◆ DrawAll()

void FiberSD::DrawAll ( )

override

Definition at line 105 of file FiberSD.cc.

105 {}

◆ EndOfEvent()

void FiberSD::EndOfEvent ( G4HCofThisEvent * HCE )

override

Definition at line 97 of file FiberSD.cc.

References clear(), and theHC.

                                              {
   edm::LogVerbatim("FiberSim") << "FiberSD: finds " << theHC->entries() << " hits";
   clear();
   edm::LogVerbatim("FiberSim") << "theHC entries = " << theHC->entries();
 }

◆ fillHits()

void FiberSD::fillHits	(	edm::PCaloHitContainer &	,
		const std::string &
	)

overridevirtual

Implements SensitiveCaloDetector.

Definition at line 127 of file FiberSD.cc.

127 {}

◆ Initialize()

void FiberSD::Initialize ( G4HCofThisEvent * HCE )

override

Definition at line 39 of file FiberSD.cc.

References bysipixelclustmulteventfilter_cfi::collectionName, theHC, and theHCID.

                                              {
   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;
 }

◆ PrintAll()

void FiberSD::PrintAll ( )

override

Definition at line 107 of file FiberSD.cc.

107 {}

◆ ProcessHits()

G4bool FiberSD::ProcessHits	(	G4Step *	aStep,
		G4TouchableHistory *	ROhist
	)

overridevirtual

Implements SensitiveDetector.

Definition at line 49 of file FiberSD.cc.

References hcalRecHitTable_cff::depth, HFShower::getHits(), hfClusterShapes_cfi::hits, mps_fire::i, position, setDetUnitId(), theHC, and theShower.

                                                               {
   //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;
 }

◆ setDetUnitId()

uint32_t FiberSD::setDetUnitId ( const G4Step * aStep )

overridevirtual

Implements SensitiveDetector.

Definition at line 119 of file FiberSD.cc.

References l1tHGCalTowerProducer_cfi::tower.

Referenced by ProcessHits().

                                                   {
   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);
 }