#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 edm::EventSetup &, 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 edm::EventSetup &es, const SensitiveDetectorCatalog &clg, edm::ParameterSet const &p)

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 edm::EventSetup &es, const SensitiveDetectorCatalog &, edm::ParameterSet const &p, 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 ::EndOfEvent *) override

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

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...

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
const SimTrackManager *	m_trackManager

FiberG4HitsCollection *	theHC

G4int	theHCID

HFShower *	theShower

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

Detailed Description

Definition at line 26 of file FiberSD.h.

Constructor & Destructor Documentation

◆ FiberSD()

FiberSD::FiberSD	(	const std::string &	iname,
		const edm::EventSetup &	es,
		const SensitiveDetectorCatalog &	clg,
		edm::ParameterSet const &	p,
		const SimTrackManager *	manager
	)

explicit

Definition at line 22 of file FiberSD.cc.

     : SensitiveCaloDetector(iname, es, clg, p),
       m_trackManager(manager),
       theShower(nullptr),
       theHCID(-1),
       theHC(nullptr) {
   // Get pointer to HcalDDDConstant and HcalDDDSimulationConstants
   edm::ESHandle<HcalDDDSimulationConstants> hdsc;
   es.get<HcalSimNumberingRecord>().get(hdsc);
   if (!hdsc.isValid()) {
     edm::LogError("FiberSim") << "FiberSD : Cannot find HcalDDDSimulationConstant";
     throw cms::Exception("Unknown", "FiberSD") << "Cannot find HcalDDDSimulationConstant\n";
   }
   const HcalDDDSimulationConstants* hsps = hdsc.product();
   edm::ESHandle<HcalDDDSimConstants> hdc;
   es.get<HcalSimNumberingRecord>().get(hdc);
   if (hdc.isValid()) {
     const HcalDDDSimConstants* hcalConstants = hdc.product();
     theShower = new HFShower(iname, hcalConstants, hsps->hcalsimpar(), p, 1);
   } else {
     edm::LogError("FiberSim") << "FiberSD : Cannot find HcalDDDSimConstant";
     throw cms::Exception("Unknown", "FiberSD") << "Cannot find HcalDDDSimConstant\n";
   }
 }

References Exception, edm::EventSetup::get(), get, HcalDDDSimulationConstants::hcalsimpar(), Calorimetry_cff::HFShower, edm::ESHandleBase::isValid(), AlCaHLTBitMon_ParallelJobs::p, edm::ESHandle< T >::product(), and theShower.

◆ ~FiberSD()

FiberSD::~FiberSD ( )

override

Definition at line 51 of file FiberSD.cc.

                   {
   delete theShower;
   delete theHC;
 }

References theHC, and theShower.

Member Function Documentation

◆ clear()

void FiberSD::clear ( void )

override

Definition at line 114 of file FiberSD.cc.

114 {}

Referenced by EndOfEvent().

◆ clearHits()

void FiberSD::clearHits ( )

overridevirtual

Implements SensitiveDetector.

Definition at line 128 of file FiberSD.cc.

128 {}

◆ DrawAll()

void FiberSD::DrawAll ( )

override

Definition at line 116 of file FiberSD.cc.

116 {}

◆ EndOfEvent()

void FiberSD::EndOfEvent ( G4HCofThisEvent * HCE )

override

Definition at line 108 of file FiberSD.cc.

                                              {
   LogDebug("FiberSim") << "FiberSD: Sees" << theHC->entries() << " hits";
   clear();
   std::cout << "theHC entries = " << theHC->entries() << std::endl;
 }

References clear(), gather_cfg::cout, LogDebug, and theHC.

◆ fillHits()

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

overridevirtual

Implements SensitiveCaloDetector.

Definition at line 138 of file FiberSD.cc.

138 {}

◆ Initialize()

void FiberSD::Initialize ( G4HCofThisEvent * HCE )

override

Definition at line 56 of file FiberSD.cc.

                                              {
   LogDebug("FiberSim") << "FiberSD : Initialize called for " << GetName();
   theHC = new FiberG4HitsCollection(GetName(), collectionName[0]);
   if (theHCID < 0)
     theHCID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
   HCE->AddHitsCollection(theHCID, theHC);
 }

References bysipixelclustmulteventfilter_cfi::collectionName, LogDebug, theHC, and theHCID.

◆ PrintAll()

void FiberSD::PrintAll ( )

override

Definition at line 118 of file FiberSD.cc.

118 {}

◆ ProcessHits()

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

overridevirtual

Implements SensitiveDetector.

Definition at line 64 of file FiberSD.cc.

                                                               {
   //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++) {
       //std::cout<<"hit position z "<<hits[i].position.z()<<std::endl;
       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());
     //std::cout<<"presteppoint position z "<<preStepPoint->GetPosition().z()<<std::endl;
  
     FiberG4Hit* aHit = new FiberG4Hit(lv, detID, depth, trackID);
     std::cout << "hit size " << hits.size() << "  npe" << aHit->npe() << std::endl;
     std::cout << "pre hit position " << aHit->hitPos() << std::endl;
     aHit->setNpe(hits.size());
     aHit->setPos(theHitPos);
     aHit->setTime(preStepPoint->GetGlobalTime());
     aHit->setPhoton(thePE);
     std::cout << "ShowerPhoton position " << thePE[0].x() << " " << thePE[0].y() << " " << thePE[0].z() << std::endl;
  
     LogDebug("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++)
       LogDebug("FiberSim") << "FiberSD: PE[" << i << "] " << thePE[i];
  
     theHC->insert(aHit);
   }
   return true;
 }

References gather_cfg::cout, LEDCalibrationChannels::depth, HFShower::getHits(), hfClusterShapes_cfi::hits, mps_fire::i, LogDebug, position, setDetUnitId(), theHC, and theShower.

◆ setDetUnitId()

uint32_t FiberSD::setDetUnitId ( const G4Step * aStep )

overridevirtual

Implements SensitiveDetector.

Definition at line 130 of file FiberSD.cc.

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