#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)

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)

	~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	(	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(), 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 ( )

override

Definition at line 34 of file FiberSD.cc.

References theHC, and theShower.

                   {
   delete theShower;
   delete theHC;
 }

Member Function Documentation

void FiberSD::clear ( void )

override

Definition at line 104 of file FiberSD.cc.

Referenced by EndOfEvent().

104 {}

void FiberSD::clearHits ( )

overridevirtual

Implements SensitiveDetector.

Definition at line 118 of file FiberSD.cc.

118 {}

void FiberSD::DrawAll ( )

override

Definition at line 106 of file FiberSD.cc.

106 {}

void FiberSD::EndOfEvent ( G4HCofThisEvent * HCE )

override

Definition at line 98 of file FiberSD.cc.

References clear(), and theHC.

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

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

overridevirtual

Implements SensitiveCaloDetector.

Definition at line 128 of file FiberSD.cc.

128 {}

void FiberSD::Initialize ( G4HCofThisEvent * HCE )

override

Definition at line 39 of file FiberSD.cc.

References ecaldqm::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;
 }

void FiberSD::PrintAll ( )

override

Definition at line 108 of file FiberSD.cc.

108 {}

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

overridevirtual

Implements SensitiveDetector.

Definition at line 50 of file FiberSD.cc.

References HLT_FULL_cff::depth, HFShower::getHits(), 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;
 }

uint32_t FiberSD::setDetUnitId ( const G4Step * aStep )

overridevirtual

Implements SensitiveDetector.

Definition at line 120 of file FiberSD.cc.

References hgcalTowerProducer_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);
 }