#include <PPSDiamondSD.h>

Inheritance diagram for PPSDiamondSD:

Public Member Functions
void	clear () override

void	clearTrack (G4Track *Track)

void	DrawAll () override

void	EndOfEvent (G4HCofThisEvent *eventHC) override

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

void	Initialize (G4HCofThisEvent *HCE) override

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

void	PrintAll () override

void	printHitInfo ()

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

G4bool	hitExists ()

void	importInfoToHit ()

bool	isPrimary (const G4Track *track)

bool	ProcessHits (G4Step step, G4TouchableHistory tHistory) override

void	resetForNewPrimary ()

uint32_t	setDetUnitId (const G4Step *step) override

void	setNumberingScheme (PPSVDetectorOrganization *scheme)

G4ThreeVector	setToLocal (const G4ThreeVector &globalPoint)

void	stepInfo (const G4Step *aStep)

void	storeHit (PPSDiamondG4Hit *)

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
PPSDiamondG4Hit *	currentHit_

G4VPhysicalVolume *	currentPV_

double	Eloss_

G4ThreeVector	entrancePoint_

int	eventno_

G4ThreeVector	exitPoint_

double	Globaltimehit_

G4int	hcID_

G4ThreeVector	hitPoint_

double	incidentEnergy_

G4String	name_

std::unique_ptr < PPSVDetectorOrganization >	numberingScheme_

double	Pabs_

int	ParentId_

short	ParticleType_

double	PhiAtEntry_

G4StepPoint *	postStepPoint_

G4StepPoint *	preStepPoint_

G4int	primaryID_

std::unique_ptr< TrackingSlaveSD >	slave_

double	theglobaltimehit_

PPSDiamondG4HitCollection *	theHC_

G4ThreeVector	theLocalEntryPoint_

G4ThreeVector	theLocalExitPoint_

int	theMPDebug_

double	thePx_

double	thePy_

double	thePz_

double	ThetaAtEntry_

G4Track *	theTrack_

double	Tof_

G4double	tSlice_

G4int	tSliceID_

unsigned int	unitID_

int	verbosity_

double	Vx_

double	Vy_

double	Vz_

Static Private Attributes
static constexpr unsigned int	maxDiamondHits_ = 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

Definition at line 31 of file PPSDiamondSD.h.

Constructor & Destructor Documentation

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

Definition at line 28 of file PPSDiamondSD.cc.

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

     : SensitiveTkDetector(name_, clg),
       numberingScheme_(nullptr),
       hcID_(-1),
       theHC_(nullptr),
       currentHit_(nullptr),
       theTrack_(nullptr),
       currentPV_(nullptr),
       unitID_(0),
       preStepPoint_(nullptr),
       postStepPoint_(nullptr),
       eventno_(0) {
   collectionName.insert(name_);
   edm::ParameterSet m_Anal = p.getParameter<edm::ParameterSet>("PPSDiamondSD");
   verbosity_ = m_Anal.getParameter<int>("Verbosity");
 
   LogDebug("PPSSimDiamond") << "*******************************************************\n"
                             << "*                                                     *\n"
                             << "* Constructing a PPSDiamondSD  with name " << name_ << "\n"
                             << "*                                                     *\n"
                             << "*******************************************************"
                             << "\n";
 
   slave_ = std::make_unique<TrackingSlaveSD>(name_);
 
   if (name_ == "CTPPSTimingHits") {
     numberingScheme_ = std::make_unique<PPSDiamondNumberingScheme>();
     edm::LogInfo("PPSSimDiamond") << "Find CTPPSDiamondHits as name";
   } else {
     edm::LogWarning("PPSSimDiamond") << "PPSDiamondSD: ReadoutName not supported\n";
   }
 
   edm::LogInfo("PPSSimDiamond") << "PPSDiamondSD: Instantiation completed";
 }

PPSDiamondSD::~PPSDiamondSD ( )

override

Definition at line 66 of file PPSDiamondSD.cc.

66 {}

Member Function Documentation

void PPSDiamondSD::clear ( void )

override

Definition at line 262 of file PPSDiamondSD.cc.

262 {}

void PPSDiamondSD::clearHits ( )

overrideprivatevirtual

Implements SensitiveDetector.

Definition at line 292 of file PPSDiamondSD.cc.

References slave_.

Referenced by update().

292 { slave_->Initialize(); }

PPSDiamondSD::slave_

std::unique_ptr< TrackingSlaveSD > slave_

Definition: PPSDiamondSD.h:57

void PPSDiamondSD::clearTrack ( G4Track * Track )

Definition at line 290 of file PPSDiamondSD.cc.

290 { track->SetTrackStatus(fStopAndKill); }

HLT_FULL_cff.track

tuple track

Definition: HLT_FULL_cff.py:11953

void PPSDiamondSD::DrawAll ( )

override

Definition at line 264 of file PPSDiamondSD.cc.

264 {}

void PPSDiamondSD::EndOfEvent ( G4HCofThisEvent * eventHC )

override

Definition at line 239 of file PPSDiamondSD.cc.

References PPSDiamondG4Hit::energyLoss(), mps_splice::entry, beamvalidation::exit(), dqmiolumiharvest::j, PPSDiamondG4Hit::localEntry(), PPSDiamondG4Hit::localExit(), maxDiamondHits_, PPSDiamondG4Hit::p(), PPSDiamondG4Hit::particleType(), PPSDiamondG4Hit::phiAtEntry(), slave_, summarize(), theHC_, PPSDiamondG4Hit::thetaAtEntry(), PPSDiamondG4Hit::tof(), PPSDiamondG4Hit::trackID(), and PPSDiamondG4Hit::unitID().

                                               {
   // here we loop over transient hits and make them persistent
   for (unsigned int j = 0; j < (unsigned int)theHC_->entries() && j < maxDiamondHits_; j++) {
     PPSDiamondG4Hit* aHit = (*theHC_)[j];
 
     Local3DPoint entry(aHit->localEntry().x(), aHit->localEntry().y(), aHit->localEntry().z());
     Local3DPoint exit(aHit->localExit().x(), aHit->localExit().y(), aHit->localExit().z());
     slave_->processHits(PSimHit(entry,
                                 exit,
                                 aHit->p(),
                                 aHit->tof(),
                                 aHit->energyLoss(),
                                 aHit->particleType(),
                                 aHit->unitID(),
                                 aHit->trackID(),
                                 aHit->thetaAtEntry(),
                                 aHit->phiAtEntry()));
   }
   summarize();
 }

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

overridevirtual

Implements SensitiveTkDetector.

Definition at line 271 of file PPSDiamondSD.cc.

References dqmiodumpmetadata::n, and slave_.

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

G4bool PPSDiamondSD::hitExists ( )

private

void PPSDiamondSD::importInfoToHit ( )

private

Definition at line 202 of file PPSDiamondSD.cc.

Referenced by ProcessHits().

                                    {
   currentHit_ = new PPSDiamondG4Hit;
   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_->setEntry(hitPoint_);
   currentHit_->setExit(exitPoint_);
   currentHit_->setLocalEntry(theLocalEntryPoint_);
   currentHit_->setLocalExit(theLocalExitPoint_);
   currentHit_->setParentId(ParentId_);
   currentHit_->setVx(Vx_);
   currentHit_->setVy(Vy_);
   currentHit_->setVz(Vz_);
   currentHit_->setPx(thePx_);
   currentHit_->setPy(thePy_);
   currentHit_->setPz(thePz_);
   currentHit_->setGlobalTimehit(Globaltimehit_);
 
   storeHit(currentHit_);
   LogDebug("PPSSimDiamond") << "STORED HIT IN: " << unitID_ << "\n";
 }

void PPSDiamondSD::Initialize ( G4HCofThisEvent * HCE )

override

Definition at line 68 of file PPSDiamondSD.cc.

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

                                                   {
   LogDebug("PPSSimDiamond") << "PPSDiamondSD : Initialize called for " << name_;
 
   theHC_ = new PPSDiamondG4HitCollection(name_, collectionName[0]);
   G4SDManager::GetSDMpointer()->AddNewCollection(name_, collectionName[0]);
   if (hcID_ < 0)
     hcID_ = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
   HCE->AddHitsCollection(hcID_, theHC_);
 }

bool PPSDiamondSD::isPrimary ( const G4Track * track )

private

Definition at line 294 of file PPSDiamondSD.cc.

References info(), and TrackInformation::isPrimary().

Referenced by stepInfo().

                                                  {
   TrackInformation* info = dynamic_cast<TrackInformation*>(track->GetUserInformation());
   return info && info->isPrimary();
 }

void PPSDiamondSD::PrintAll ( )

override

Definition at line 266 of file PPSDiamondSD.cc.

References LogDebug, and theHC_.

                             {
   LogDebug("PPSSimDiamond") << "PPSDiamond: Collection " << theHC_->GetName() << "\n";
   theHC_->PrintAllHits();
 }

void PPSDiamondSD::printHitInfo ( )

Definition at line 78 of file PPSDiamondSD.cc.

References currentPV_, Eloss_, LogDebug, postStepPoint_, preStepPoint_, primaryID_, gpuVertexFinder::printf(), theTrack_, tSliceID_, and unitID_.

                                 {
   LogDebug("PPSSimDiamond") << theTrack_->GetDefinition()->GetParticleName() << " PPS_Timing_SD CreateNewHit for"
                             << " PV " << currentPV_->GetName() << " PVid = " << currentPV_->GetCopyNo() << " Unit "
                             << unitID_ << "\n";
   LogDebug("PPSSimDiamond") << " primary " << primaryID_ << " time slice " << tSliceID_ << " of energy "
                             << theTrack_->GetTotalEnergy() << " Eloss_ " << Eloss_ << " positions "
                             << "\n";
   printf(" PreStepPoint(%10f,%10f,%10f)",
          preStepPoint_->GetPosition().x(),
          preStepPoint_->GetPosition().y(),
          preStepPoint_->GetPosition().z());
   printf(" PosStepPoint(%10f,%10f,%10f)\n",
          postStepPoint_->GetPosition().x(),
          postStepPoint_->GetPosition().y(),
          postStepPoint_->GetPosition().z());
   LogDebug("PPSSimDiamond") << " positions "
                             << "(" << postStepPoint_->GetPosition().x() << "," << postStepPoint_->GetPosition().y()
                             << "," << postStepPoint_->GetPosition().z() << ")"
                             << " For Track  " << theTrack_->GetTrackID() << " which is a "
                             << theTrack_->GetDefinition()->GetParticleName() << " ParentID is "
                             << theTrack_->GetParentID() << "\n";
 
   if (theTrack_->GetTrackID() == 1) {
     LogDebug("PPSSimDiamond") << " primary particle ";
   } else {
     LogDebug("PPSSimDiamond") << " daughter of part. " << theTrack_->GetParentID();
   }
 
   LogDebug("PPSSimDiamond") << " and created by ";
 
   if (theTrack_->GetCreatorProcess() != nullptr)
     LogDebug("PPSSimDiamond") << theTrack_->GetCreatorProcess()->GetProcessName();
   else
     LogDebug("PPSSimDiamond") << "NO process";
 
   LogDebug("PPSSimDiamond") << "\n";
 }

bool PPSDiamondSD::ProcessHits	(	G4Step *	step,
		G4TouchableHistory *	tHistory
	)

overrideprivatevirtual

Implements SensitiveDetector.

Definition at line 116 of file PPSDiamondSD.cc.

References importInfoToHit(), LogDebug, name_, stepInfo(), and theTrack_.

                                                                  {
   if (aStep == nullptr) {
     LogDebug("PPSSimDiamond") << " There is no hit to process ";
     return true;
   } else {
     LogDebug("PPSSimDiamond") << "*******************************************************\n"
                               << "*                                                     *\n"
                               << "* PPS Diamond Hit initialized  with name " << name_ << "\n"
                               << "*                                                     *\n"
                               << "*******************************************************"
                               << "\n";
 
     stepInfo(aStep);
 
     if (theTrack_->GetDefinition()->GetPDGEncoding() == 2212) {
       importInfoToHit();  //in addtion to import info to hit it STORE hit as well
       LogDebug("PPSSimDiamond") << " information imported to the hit ";
     }
 
     return true;
   }
 }

void PPSDiamondSD::resetForNewPrimary ( )

private

uint32_t PPSDiamondSD::setDetUnitId ( const G4Step * step )

overrideprivatevirtual

Implements SensitiveDetector.

Definition at line 187 of file PPSDiamondSD.cc.

References numberingScheme_.

Referenced by stepInfo().

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

void PPSDiamondSD::setNumberingScheme ( PPSVDetectorOrganization * scheme )

private

Definition at line 275 of file PPSDiamondSD.cc.

References LogDebug, and numberingScheme_.

                                                                       {
   if (scheme) {
     LogDebug("PPSDiamond") << "PPSDiamondSD: updates numbering scheme for " << GetName();
     numberingScheme_.reset(scheme);
   }
 }

G4ThreeVector PPSDiamondSD::setToLocal ( const G4ThreeVector & globalPoint )

private

Definition at line 231 of file PPSDiamondSD.cc.

References preStepPoint_.

Referenced by stepInfo().

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

void PPSDiamondSD::stepInfo ( const G4Step * aStep )

private

Definition at line 139 of file PPSDiamondSD.cc.

References SensitiveDetector::ConvertToLocal3DPoint(), currentPV_, Eloss_, exitPoint_, GeV, hitPoint_, incidentEnergy_, isPrimary(), LogDebug, Pabs_, ParentId_, ParticleType_, PV3DBase< T, PVType, FrameType >::phi(), PhiAtEntry_, postStepPoint_, preStepPoint_, primaryID_, setDetUnitId(), setToLocal(), theglobaltimehit_, theLocalEntryPoint_, theLocalExitPoint_, thePx_, thePy_, thePz_, PV3DBase< T, PVType, FrameType >::theta(), ThetaAtEntry_, theTrack_, Tof_, tSlice_, tSliceID_, unitID_, verbosity_, Vx_, Vy_, and Vz_.

Referenced by ProcessHits().

                                                {
   theTrack_ = aStep->GetTrack();
   preStepPoint_ = aStep->GetPreStepPoint();
   postStepPoint_ = aStep->GetPostStepPoint();
   hitPoint_ = preStepPoint_->GetPosition();
   exitPoint_ = postStepPoint_->GetPosition();
   currentPV_ = preStepPoint_->GetPhysicalVolume();
   theLocalEntryPoint_ = setToLocal(hitPoint_);
   theLocalExitPoint_ = setToLocal(exitPoint_);
   theglobaltimehit_ = preStepPoint_->GetGlobalTime() / nanosecond;
   incidentEnergy_ = (aStep->GetPreStepPoint()->GetTotalEnergy() / eV);
   tSlice_ = (postStepPoint_->GetGlobalTime()) / nanosecond;
   tSliceID_ = (int)tSlice_;
   unitID_ = setDetUnitId(aStep);
 
   if (verbosity_)
     LogDebug("PPSSimDiamond") << "UNIT " << unitID_ << "\n";
 
   primaryID_ = theTrack_->GetTrackID();
   Pabs_ = (aStep->GetPreStepPoint()->GetMomentum().mag()) / GeV;
   thePx_ = (aStep->GetPreStepPoint()->GetMomentum().x()) / GeV;
   thePy_ = (aStep->GetPreStepPoint()->GetMomentum().y()) / GeV;
   thePz_ = (aStep->GetPreStepPoint()->GetMomentum().z()) / GeV;
   Tof_ = aStep->GetPreStepPoint()->GetGlobalTime() / nanosecond;
   Eloss_ = (aStep->GetPreStepPoint()->GetTotalEnergy() / eV);  //pps added
   ParticleType_ = theTrack_->GetDefinition()->GetPDGEncoding();
 
   //corrected phi and theta treatment
   G4ThreeVector gmd = aStep->GetPreStepPoint()->GetMomentumDirection();
   // convert it to local frame
   G4ThreeVector lmd = ((G4TouchableHistory*)(aStep->GetPreStepPoint()->GetTouchable()))
                           ->GetHistory()
                           ->GetTopTransform()
                           .TransformAxis(gmd);
   Local3DPoint lnmd = ConvertToLocal3DPoint(lmd);
   ThetaAtEntry_ = lnmd.theta();
   PhiAtEntry_ = lnmd.phi();
 
   if (isPrimary(theTrack_))
     ParentId_ = 0;
   else
     ParentId_ = theTrack_->GetParentID();
 
   Vx_ = theTrack_->GetVertexPosition().x() / mm;
   Vy_ = theTrack_->GetVertexPosition().y() / mm;
   Vz_ = theTrack_->GetVertexPosition().z() / mm;
 }

void PPSDiamondSD::storeHit ( PPSDiamondG4Hit * hit )

private

Definition at line 191 of file PPSDiamondSD.cc.

References LogDebug, theHC_, and verbosity_.

Referenced by importInfoToHit().

                                                 {
   if (hit == nullptr) {
     if (verbosity_)
       LogDebug("PPSSimDiamond") << "PPSDiamond: hit to be stored is NULL !!"
                                 << "\n";
     return;
   }
 
   theHC_->insert(hit);
 }

void PPSDiamondSD::summarize ( )

private

Definition at line 260 of file PPSDiamondSD.cc.

Referenced by EndOfEvent().

260 {}

void PPSDiamondSD::update ( const BeginOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 281 of file PPSDiamondSD.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("PPSSimDiamond") << " Dispatched BeginOfEvent !"
                             << "\n";
   clearHits();
   eventno_ = (*i)()->GetEventID();
 }