#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 DDCompactView &, 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
	SensitiveCaloDetector (const std::string &iname, const DDCompactView &cpv, 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

	SensitiveDetector (const std::string &iname, const DDCompactView &cpv, const SensitiveDetectorCatalog &, edm::ParameterSet const &p)

	~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
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	(	const std::string &	iname,
		const DDCompactView &	cpv,
		const SensitiveDetectorCatalog &	clg,
		edm::ParameterSet const &	p,
		const SimTrackManager *	manager
	)

explicit

Definition at line 21 of file FiberSD.cc.

References g4SimHits_cfi::HFShower, and theShower.

                                          :
   SensitiveCaloDetector(iname, cpv, clg, p),
   m_trackManager(manager), theHCID(-1), theHC(nullptr) {
 
   theShower = new HFShower(iname, cpv, p, 1);
 
 }

FiberSD::~FiberSD ( )

override

Definition at line 31 of file FiberSD.cc.

References theHC, and theShower.

                   {
  
   delete theShower;
   delete theHC;
 }

Member Function Documentation

void FiberSD::clear ( void )

override

Definition at line 105 of file FiberSD.cc.

Referenced by EndOfEvent().

105 {}

void FiberSD::clearHits ( )

overridevirtual

Implements SensitiveDetector.

Definition at line 132 of file FiberSD.cc.

132 {}

void FiberSD::DrawAll ( )

override

Definition at line 107 of file FiberSD.cc.

107 {}

void FiberSD::EndOfEvent ( G4HCofThisEvent * HCE )

override

Definition at line 98 of file FiberSD.cc.

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

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

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

overridevirtual

Implements SensitiveCaloDetector.

Definition at line 142 of file FiberSD.cc.

142 {}

void FiberSD::Initialize ( G4HCofThisEvent * HCE )

override

Definition at line 37 of file FiberSD.cc.

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

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

void FiberSD::PrintAll ( )

override

Definition at line 109 of file FiberSD.cc.

109 {}

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

overridevirtual

Implements SensitiveDetector.

Definition at line 47 of file FiberSD.cc.

References gather_cfg::cout, particleFlowClusterECALTimeSelected_cfi::depth, HFShower::getHits(), hfClusterShapes_cfi::hits, mps_fire::i, LogDebug, 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++) {
       //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;
 }

uint32_t FiberSD::setDetUnitId ( const G4Step * aStep )

overridevirtual

Implements SensitiveDetector.

Definition at line 134 of file FiberSD.cc.

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

void FiberSD::update ( const BeginOfJob * )

overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 111 of file FiberSD.cc.

References Exception, edm::EventSetup::get(), HFShower::initRun(), and theShower.

                                            {
 
   const edm::EventSetup* es = (*job)();
   edm::ESHandle<HcalDDDSimConstants>    hdc;
   es->get<HcalSimNumberingRecord>().get(hdc);
   if (hdc.isValid()) {
     HcalDDDSimConstants *hcalConstants = (HcalDDDSimConstants*)(&(*hdc));
     theShower->initRun(nullptr, hcalConstants);
   } else {
     edm::LogError("HcalSim") << "HCalSD : Cannot find HcalDDDSimConstant";
     throw cms::Exception("Unknown", "HCalSD") << "Cannot find HcalDDDSimConstant" << "\n";
   }
 
 }

void FiberSD::update ( const BeginOfRun * )

overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 126 of file FiberSD.cc.

126 {}

void FiberSD::update ( const BeginOfEvent * )

overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 128 of file FiberSD.cc.

128 {}

void FiberSD::update ( const ::EndOfEvent * )

overrideprotected

Definition at line 130 of file FiberSD.cc.

130 {}

Member Data Documentation

const SimTrackManager* FiberSD::m_trackManager

private

Definition at line 59 of file FiberSD.h.

FiberG4HitsCollection* FiberSD::theHC

private

Definition at line 63 of file FiberSD.h.

Referenced by EndOfEvent(), Initialize(), ProcessHits(), and ~FiberSD().

G4int FiberSD::theHCID

private

Definition at line 62 of file FiberSD.h.

Referenced by Initialize().

HFShower* FiberSD::theShower

private

Definition at line 60 of file FiberSD.h.

Referenced by FiberSD(), ProcessHits(), update(), and ~FiberSD().

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation