#include <SimG4CMS/PPS/interface/PPSPixelSD.h>

Inheritance diagram for PPSPixelSD:

Public Member Functions
void	clear () override

void	DrawAll () override

void	EndOfEvent (G4HCofThisEvent *eventHC) override

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

void	Initialize (G4HCofThisEvent *HCE) override

	PPSPixelSD (const std::string &, const SensitiveDetectorCatalog &, edm::ParameterSet const &, SimTrackManager const *)

void	PrintAll () override

bool	ProcessHits (G4Step , G4TouchableHistory ) override

uint32_t	setDetUnitId (const G4Step *) override

	~PPSPixelSD () override

Public Member Functions inherited from SensitiveTkDetector
	SensitiveTkDetector (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 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 ()

Private Member Functions
void	clearHits () override

void	createNewHit ()

bool	hitExists ()

void	resetForNewPrimary ()

G4ThreeVector	setToLocal (const G4ThreeVector &globalPoint)

void	stepInfo (const G4Step *aStep)

void	storeHit (PPSPixelG4Hit *)

void	summarize ()

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

void	update (const ::EndOfEvent *) override

void	updateHit ()

Private Attributes
PPSPixelG4Hit *	currentHit_

G4VPhysicalVolume *	currentPV_

float	edeposit_

float	Eloss_

G4ThreeVector	entrancePoint_

int	eventno_

G4int	hcID_

G4ThreeVector	hitPoint_

float	incidentEnergy_

std::string	name_

std::unique_ptr < PPSVDetectorOrganization >	numberingScheme_

float	Pabs_

int	ParentId_

short	ParticleType_

float	PhiAtEntry_

G4ThreeVector	position_

G4StepPoint *	postStepPoint_

G4StepPoint *	preStepPoint_

uint32_t	previousUnitID_

int	primaryID_

G4int	primID_

std::unique_ptr< TrackingSlaveSD >	slave_

G4ThreeVector	theEntryPoint_

G4ThreeVector	theExitPoint_

PPSPixelG4HitCollection *	theHC_

const SimTrackManager *	theManager_

float	ThetaAtEntry_

G4Track *	theTrack_

float	Tof_

int	tsID_

double	tSlice_

int	tSliceID_

uint32_t	unitID_

float	Vx_

float	Vy_

float	Vz_

Static Private Attributes
static constexpr unsigned int	maxPixelHits_ = 15000

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

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

Detailed Description

Description: Stores hits of PPSPixel in appropriate container

Usage: Used in sensitive detector builder

Definition at line 46 of file PPSPixelSD.h.

Constructor & Destructor Documentation

PPSPixelSD::PPSPixelSD	(	const std::string &	name_,
		const SensitiveDetectorCatalog &	clg,
		edm::ParameterSet const &	p,
		SimTrackManager const *	manager
	)

Definition at line 37 of file PPSPixelSD.cc.

References ecaldqm::collectionName, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), name_, numberingScheme_, and slave_.

     : SensitiveTkDetector(name_, clg),
       numberingScheme_(nullptr),
       hcID_(-1),
       theHC_(nullptr),
       theManager_(manager),
       currentHit_(nullptr),
       theTrack_(nullptr),
       currentPV_(nullptr),
       unitID_(0),
       previousUnitID_(0),
       preStepPoint_(nullptr),
       postStepPoint_(nullptr),
       eventno_(0) {
   //Add PPS Sentitive Detector Names
   collectionName.insert(name_);
 
   //Parameters
   edm::ParameterSet m_p = p.getParameter<edm::ParameterSet>("PPSPixelSD");
   int verbn = m_p.getUntrackedParameter<int>("Verbosity");
   SetVerboseLevel(verbn);
   slave_ = std::make_unique<TrackingSlaveSD>(name_);
   if (name_ == "CTPPSPixelHits") {
     numberingScheme_ = std::make_unique<PPSPixelNumberingScheme>();
   } else {
     edm::LogWarning("PPSSim") << "PPSPixelSD: ReadoutName not supported\n";
   }
 
   edm::LogInfo("PPSSim") << "PPSPixelSD: Instantiation completed";
 }

PPSPixelSD::~PPSPixelSD ( )

override

Definition at line 71 of file PPSPixelSD.cc.

71 {}

Member Function Documentation

void PPSPixelSD::clear ( void )

override

Definition at line 129 of file PPSPixelSD.cc.

129 {}

void PPSPixelSD::clearHits ( )

overrideprivatevirtual

Implements SensitiveDetector.

Definition at line 152 of file PPSPixelSD.cc.

References slave_.

Referenced by update().

152 { slave_->Initialize(); }

PPSPixelSD::slave_

std::unique_ptr< TrackingSlaveSD > slave_

Definition: PPSPixelSD.h:81

void PPSPixelSD::createNewHit ( )

private

Definition at line 246 of file PPSPixelSD.cc.

Referenced by ProcessHits().

                               {
 #ifdef debug
   LogDebug("PPSSim") << "PPSPixelSD CreateNewHit for"
                      << " PV " << currentPV_->GetName() << " PVid = " << currentPV_->GetCopyNo()
                      << " MVid = " << currentPV_->GetMother()->GetCopyNo() << " Unit " << unitID_ << "\n"
                      << " primary " << primaryID_ << " time slice " << tSliceID_ << " For Track  "
                      << theTrack_->GetTrackID() << " which is a " << theTrack_->GetDefinition()->GetParticleName();
 
   if (theTrack_->GetTrackID() == 1) {
     LogDebug("PPSSim") << " of energy " << theTrack_->GetTotalEnergy();
   } else {
     LogDebug("PPSSim") << " daughter of part. " << theTrack_->GetParentID();
   }
 
   if (theTrack_->GetCreatorProcess() != NULL)
     LogDebug("PPSSim") << theTrack_->GetCreatorProcess()->GetProcessName();
   else
     LogDebug("PPSSim") << "NO process";
 #endif
 
   currentHit_ = new PPSPixelG4Hit;
   currentHit_->setTrackID(primaryID_);
   currentHit_->setTimeSlice(tSlice_);
   currentHit_->setUnitID(unitID_);
   currentHit_->setIncidentEnergy(incidentEnergy_);
 
   currentHit_->setP(Pabs_);
   currentHit_->setTof(Tof_);
   currentHit_->setEnergyLoss(Eloss_);
   currentHit_->setParticleType(ParticleType_);
   currentHit_->setThetaAtEntry(ThetaAtEntry_);
   currentHit_->setPhiAtEntry(PhiAtEntry_);
 
   currentHit_->setMeanPosition(position_);
   currentHit_->setEntryPoint(theEntryPoint_);
   currentHit_->setExitPoint(theExitPoint_);
 
   currentHit_->setParentId(ParentId_);
   currentHit_->setVx(Vx_);
   currentHit_->setVy(Vy_);
   currentHit_->setVz(Vz_);
 
   updateHit();
 
   storeHit(currentHit_);
 }

void PPSPixelSD::DrawAll ( )

override

Definition at line 131 of file PPSPixelSD.cc.

131 {}

void PPSPixelSD::EndOfEvent ( G4HCofThisEvent * eventHC )

override

Definition at line 101 of file PPSPixelSD.cc.

References PPSPixelG4Hit::energyLoss(), PPSPixelG4Hit::entryPoint(), PPSPixelG4Hit::exitPoint(), dqmiolumiharvest::j, LogDebug, maxPixelHits_, PPSPixelG4Hit::p(), PPSPixelG4Hit::particleType(), PPSPixelG4Hit::phiAtEntry(), slave_, summarize(), theHC_, PPSPixelG4Hit::thetaAtEntry(), PPSPixelG4Hit::tof(), PPSPixelG4Hit::trackID(), PPSPixelG4Hit::unitID(), and PPSPixelG4Hit::z().

                                             {
   // here we loop over transient hits and make them persistent
   for (unsigned int j = 0; j < (unsigned int)theHC_->entries() && j < maxPixelHits_; j++) {
     PPSPixelG4Hit* aHit = (*theHC_)[j];
 #ifdef debug
     LogDebug("PPSSim") << "HIT NUMERO " << j << "unit ID = " << aHit->unitID() << "\n"
                        << "               "
                        << "enrty z " << aHit->entry().z() << "\n"
                        << "               "
                        << "theta   " << aHit->thetaAtEntry() << "\n";
 #endif
 
     Local3DPoint Enter(aHit->entryPoint().x(), aHit->entryPoint().y(), aHit->entryPoint().z());
     Local3DPoint Exit(aHit->exitPoint().x(), aHit->exitPoint().y(), aHit->exitPoint().z());
     slave_->processHits(PSimHit(Enter,
                                 Exit,
                                 aHit->p(),
                                 aHit->tof(),
                                 aHit->energyLoss(),
                                 aHit->particleType(),
                                 aHit->unitID(),
                                 aHit->trackID(),
                                 aHit->thetaAtEntry(),
                                 aHit->phiAtEntry()));
   }
   summarize();
 }

void PPSPixelSD::fillHits	(	edm::PSimHitContainer &	c,
		const std::string &	n
	)

overridevirtual

Implements SensitiveTkDetector.

Definition at line 138 of file PPSPixelSD.cc.

References dqmiodumpmetadata::n, and slave_.

                                                                     {
   if (slave_->name() == n) {
     c = slave_->hits();
   }
 }

bool PPSPixelSD::hitExists ( )

private

Definition at line 211 of file PPSPixelSD.cc.

References currentHit_, newFWLiteAna::found, dqmiolumiharvest::j, nhits, previousUnitID_, primaryID_, primID_, resetForNewPrimary(), theHC_, PPSPixelG4Hit::timeSliceID(), PPSPixelG4Hit::trackID(), tsID_, tSliceID_, PPSPixelG4Hit::unitID(), unitID_, and updateHit().

Referenced by ProcessHits().

                            {
   if (primaryID_ < 1) {
     edm::LogWarning("PPSSim") << "***** PPSPixelSD error: primaryID = " << primaryID_ << " maybe detector name changed";
   }
 
   // Update if in the same detector, time-slice and for same track
   if (tSliceID_ == tsID_ && unitID_ == previousUnitID_) {
     updateHit();
     return true;
   }
 
   // Reset entry point for new primary
   if (primaryID_ != primID_)
     resetForNewPrimary();
 
   //look in the HitContainer whether a hit with the same primID_, unitID_,
   //tSliceID_ already exists:
   bool found = false;
   int nhits = theHC_->entries();
   for (int j = 0; j < nhits && !found; j++) {
     PPSPixelG4Hit* aPreviousHit = (*theHC_)[j];
     if (aPreviousHit->trackID() == primaryID_ && aPreviousHit->timeSliceID() == tSliceID_ &&
         aPreviousHit->unitID() == unitID_) {
       currentHit_ = aPreviousHit;
       found = true;
     }
   }
 
   if (found) {
     updateHit();
     return true;
   }
   return false;
 }

void PPSPixelSD::Initialize ( G4HCofThisEvent * HCE )

override

Definition at line 87 of file PPSPixelSD.cc.

References ecaldqm::collectionName, hcID_, LogDebug, name_, primID_, theHC_, and tsID_.

                                                 {
   LogDebug("PPSSim") << "PPSPixelSD : Initialize called for " << name_;
 
   theHC_ = new PPSPixelG4HitCollection(GetName(), collectionName[0]);
   G4SDManager::GetSDMpointer()->AddNewCollection(name_, collectionName[0]);
 
   if (hcID_ < 0)
     hcID_ = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
   HCE->AddHitsCollection(hcID_, theHC_);
 
   tsID_ = -2;
   primID_ = -2;
 }

void PPSPixelSD::PrintAll ( )

override

Definition at line 133 of file PPSPixelSD.cc.

References LogDebug, and theHC_.

                           {
   LogDebug("PPSSim") << "PPSPixelSD: Collection " << theHC_->GetName();
   theHC_->PrintAllHits();
 }

bool PPSPixelSD::ProcessHits	(	G4Step *	aStep,
		G4TouchableHistory *
	)

overridevirtual

Implements SensitiveDetector.

Definition at line 73 of file PPSPixelSD.cc.

References createNewHit(), edeposit_, hitExists(), and stepInfo().

                                                                {
   if (!aStep)
     return true;
 
   stepInfo(aStep);
   if (!hitExists() && edeposit_ > 0.)
     createNewHit();
   return true;
 }

void PPSPixelSD::resetForNewPrimary ( )

private

Definition at line 318 of file PPSPixelSD.cc.

References entrancePoint_, hitPoint_, incidentEnergy_, preStepPoint_, and setToLocal().

Referenced by hitExists().

                                     {
   entrancePoint_ = setToLocal(hitPoint_);
 
   incidentEnergy_ = preStepPoint_->GetKineticEnergy();
 }

uint32_t PPSPixelSD::setDetUnitId ( const G4Step * aStep )

overridevirtual

Implements SensitiveDetector.

Definition at line 83 of file PPSPixelSD.cc.

References numberingScheme_.

Referenced by stepInfo().

                                                      {
   return (numberingScheme_ == nullptr ? 0 : numberingScheme_->unitID(aStep));
 }

G4ThreeVector PPSPixelSD::setToLocal ( const G4ThreeVector & globalPoint )

private

Definition at line 154 of file PPSPixelSD.cc.

References preStepPoint_.

Referenced by resetForNewPrimary().

                                                                 {
   G4ThreeVector localPoint;
   const G4VTouchable* touch = preStepPoint_->GetTouchable();
   localPoint = touch->GetHistory()->GetTopTransform().TransformPoint(global);
   return localPoint;
 }

void PPSPixelSD::stepInfo ( const G4Step * aStep )

private

Definition at line 161 of file PPSPixelSD.cc.

References currentPV_, edeposit_, Eloss_, SensitiveDetector::FinalStepPosition(), hitPoint_, SensitiveDetector::InitialStepPosition(), SensitiveDetector::LocalCoordinates, LogDebug, name_, Pabs_, ParentId_, HLT_FULL_cff::particleType, ParticleType_, PhiAtEntry_, position_, postStepPoint_, preStepPoint_, primaryID_, setDetUnitId(), theEntryPoint_, theExitPoint_, ThetaAtEntry_, theTrack_, Tof_, tSlice_, tSliceID_, unitID_, Vx_, Vy_, Vz_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by ProcessHits().

                                              {
   preStepPoint_ = aStep->GetPreStepPoint();
   postStepPoint_ = aStep->GetPostStepPoint();
   theTrack_ = aStep->GetTrack();
   Local3DPoint TheEntryPoint = SensitiveDetector::InitialStepPosition(aStep, LocalCoordinates);
   Local3DPoint TheExitPoint = SensitiveDetector::FinalStepPosition(aStep, LocalCoordinates);
 
 #ifdef _PRINT_HITS_
   LogDebug("PPSSim") << "theEntryPoint_ " << TheEntryPoint << "\n";
   LogDebug("PPSSim") << "position " << preStepPoint_->GetPosition() << "\n";
 #endif
   hitPoint_ = preStepPoint_->GetPosition();
   currentPV_ = preStepPoint_->GetPhysicalVolume();
 
   G4String name_ = currentPV_->GetName();
   name_.assign(name_, 0, 4);
   G4String particleType = theTrack_->GetDefinition()->GetParticleName();
   edeposit_ = aStep->GetTotalEnergyDeposit();
 
   tSlice_ = (postStepPoint_->GetGlobalTime()) / nanosecond;
   tSliceID_ = (int)tSlice_;
   unitID_ = setDetUnitId(aStep);
 #ifdef debug
   LogDebug("PPSSim") << "UNIT ID " << unitID_;
 #endif
   primaryID_ = theTrack_->GetTrackID();
 
   theEntryPoint_.setX(TheEntryPoint.x());
   theEntryPoint_.setY(TheEntryPoint.y());
   theEntryPoint_.setZ(TheEntryPoint.z());
   theExitPoint_.setX(TheExitPoint.x());
   theExitPoint_.setY(TheExitPoint.y());
   theExitPoint_.setZ(TheExitPoint.z());
 
   position_ = hitPoint_;
   Pabs_ = aStep->GetPreStepPoint()->GetMomentum().mag() / GeV;
   Tof_ = aStep->GetPostStepPoint()->GetGlobalTime() / nanosecond;
 
   Eloss_ = aStep->GetTotalEnergyDeposit() / GeV;
   ParticleType_ = theTrack_->GetDefinition()->GetPDGEncoding();
 
   ThetaAtEntry_ = aStep->GetPreStepPoint()->GetPosition().theta();
   PhiAtEntry_ = aStep->GetPreStepPoint()->GetPosition().phi();
 
   ParentId_ = theTrack_->GetParentID();
   Vx_ = theTrack_->GetVertexPosition().x();
   Vy_ = theTrack_->GetVertexPosition().y();
   Vz_ = theTrack_->GetVertexPosition().z();
 }

void PPSPixelSD::storeHit ( PPSPixelG4Hit * hit )

private

Definition at line 307 of file PPSPixelSD.cc.

References primID_, and theHC_.

Referenced by createNewHit().

                                             {
   if (primID_ < 0)
     return;
   if (hit == nullptr) {
     edm::LogWarning("PPSSim") << "PPSPixelSD: hit to be stored is NULL !!";
     return;
   }
 
   theHC_->insert(hit);
 }

void PPSPixelSD::summarize ( )

private

Definition at line 324 of file PPSPixelSD.cc.

Referenced by EndOfEvent().

324 {}

void PPSPixelSD::update ( const BeginOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 144 of file PPSPixelSD.cc.

References clearHits(), eventno_, and LogDebug.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                              {
   LogDebug("PPSSim") << " Dispatched BeginOfEvent for " << GetName() << " !";
   clearHits();
   eventno_ = (*i)()->GetEventID();
 }

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

overrideprivate

Definition at line 150 of file PPSPixelSD.cc.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

150 {}

void PPSPixelSD::updateHit ( )

private

Definition at line 293 of file PPSPixelSD.cc.

References currentHit_, Eloss_, LogDebug, postStepPoint_, previousUnitID_, primaryID_, primID_, PPSPixelG4Hit::setEnergyLoss(), tsID_, tSliceID_, and unitID_.

Referenced by createNewHit(), and hitExists().

                            {
   if (Eloss_ > 0.) {
 #ifdef debug
     LogDebug("PPSSim") << "G4PPSPixelSD updateHit: add eloss " << Eloss_ << "\nCurrentHit=" << currentHit_
                        << ", PostStepPoint=" << postStepPoint_->GetPosition();
 #endif
     currentHit_->setEnergyLoss(Eloss_);
   }
   // buffer for next steps:
   tsID_ = tSliceID_;
   primID_ = primaryID_;
   previousUnitID_ = unitID_;
 }