#include <MtdSD.h>

Inheritance diagram for MtdSD:

Public Member Functions
	MtdSD (const std::string &, const SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager *)

uint32_t	setDetUnitId (const G4Step *) override

	~MtdSD () override

Public Member Functions inherited from TimingSD
void	clearHits () override

void	EndOfEvent (G4HCofThisEvent *eventHC) override

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

void	Initialize (G4HCofThisEvent *HCE) override

void	PrintAll () override

bool	ProcessHits (G4Step , G4TouchableHistory ) override

	TimingSD (const std::string &, const SensitiveDetectorCatalog &, const SimTrackManager *)

	~TimingSD () 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, const std::string &newcollname="")

	~SensitiveDetector () override

Public Member Functions inherited from Observer< const BeginOfEvent *>
	Observer ()

void	slotForUpdate (const BeginOfEvent * iT)

virtual	~Observer ()

Protected Member Functions
int	getTrackID (const G4Track *) override

Protected Member Functions inherited from TimingSD
virtual bool	checkHit (const G4Step , BscG4Hit )

const G4ThreeVector &	getGlobalEntryPoint () const

const G4ThreeVector &	getLocalEntryPoint () const

void	setCuts (double eCut, double historyCut)

void	setTimeFactor (double)

void	setToLocal (const G4StepPoint *stepPoint, const G4ThreeVector &globalPoint, G4ThreeVector &localPoint)

void	update (const BeginOfEvent *) 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 > &)

Private Member Functions
void	getBaseNumber (const G4Step *)

void	setNumberingScheme (MTDNumberingScheme *)

Private Attributes
double	energyCut

double	energyHistoryCut

bool	isBTL

bool	isETL

MTDNumberingScheme *	numberingScheme

MTDBaseNumber	theBaseNumber

Static Private Attributes
static constexpr unsigned int	k_idFromCaloOffset = 3

static constexpr unsigned int	k_idloopOffset = 2

static constexpr unsigned int	k_idsecOffset = 1

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

Detailed Description

Definition at line 19 of file MtdSD.h.

Constructor & Destructor Documentation

◆ MtdSD()

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

Definition at line 20 of file MtdSD.cc.

References energyCut, energyHistoryCut, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), isBTL, isETL, Skims_PA_cff::name, AlCaHLTBitMon_ParallelJobs::p, generator_cfi::scheme, TimingSD::setCuts(), setNumberingScheme(), and TimingSD::setTimeFactor().

     : TimingSD(name, clg, manager), numberingScheme(nullptr) {
   //Parameters
   edm::ParameterSet m_p = p.getParameter<edm::ParameterSet>("MtdSD");
   int verbn = m_p.getUntrackedParameter<int>("Verbosity");
 
   SetVerboseLevel(verbn);
 
   MTDNumberingScheme* scheme = nullptr;
   if (name == "FastTimerHitsBarrel") {
     scheme = dynamic_cast<MTDNumberingScheme*>(new BTLNumberingScheme());
     isBTL = true;
   } else if (name == "FastTimerHitsEndcap") {
     scheme = dynamic_cast<MTDNumberingScheme*>(new ETLNumberingScheme());
     isETL = true;
   } else {
     scheme = nullptr;
     edm::LogWarning("MtdSim") << "MtdSD: ReadoutName not supported";
   }
   if (scheme)
     setNumberingScheme(scheme);
 
   energyCut = m_p.getParameter<double>("EnergyThresholdForPersistencyInGeV") * CLHEP::GeV;     //default must be 0.5
   energyHistoryCut = m_p.getParameter<double>("EnergyThresholdForHistoryInGeV") * CLHEP::GeV;  //default must be 0.05
 
   setCuts(energyCut, energyHistoryCut);
 
   double newTimeFactor = 1. / m_p.getParameter<double>("TimeSliceUnit");
   edm::LogVerbatim("MtdSim") << "New time factor = " << newTimeFactor;
   setTimeFactor(newTimeFactor);
 
   edm::LogVerbatim("MtdSim") << "MtdSD: Instantiation completed for " << name;
 }

◆ ~MtdSD()

MtdSD::~MtdSD ( )

override

Definition at line 57 of file MtdSD.cc.

57 {}

Member Function Documentation

◆ getBaseNumber()

void MtdSD::getBaseNumber ( const G4Step * aStep )

private

Definition at line 80 of file MtdSD.cc.

References MTDBaseNumber::addLevel(), MTDBaseNumber::getCapacity(), cuy::ii, MTDBaseNumber::reset(), MTDBaseNumber::setSize(), and theBaseNumber.

Referenced by setDetUnitId().

                                              {
   theBaseNumber.reset();
   const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
   int theSize = touch->GetHistoryDepth() + 1;
   if (theBaseNumber.getCapacity() < theSize)
     theBaseNumber.setSize(theSize);
   //Get name and copy numbers
   if (theSize > 1) {
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("MtdSim") << "Building MTD basenumber:";
 #endif
     for (int ii = 0; ii < theSize; ii++) {
       theBaseNumber.addLevel(touch->GetVolume(ii)->GetName(), touch->GetReplicaNumber(ii));
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("MtdSim") << "MtdSD::getBaseNumber(): Adding level " << ii << ": "
                                  << touch->GetVolume(ii)->GetName() << "[" << touch->GetReplicaNumber(ii) << "]";
 #endif
     }
   }
 }

◆ getTrackID()

int MtdSD::getTrackID ( const G4Track * aTrack )

overrideprotectedvirtual

Reimplemented from TimingSD.

Definition at line 101 of file MtdSD.cc.

References PSimHit::addTrackIdOffset(), SensitiveDetector::cmsTrackInformation(), Exception, TrackInformation::getIDonCaloSurface(), TrackInformation::isBTLlooper(), TrackInformation::isExtSecondary(), TrackInformation::isInTrkFromBackscattering(), k_idFromCaloOffset, k_idloopOffset, k_idsecOffset, TrackInformation::mcTruthID(), TrackInformation::Print(), and rname.

                                            {
   int theID = aTrack->GetTrackID();
   TrackInformation* trkInfo = cmsTrackInformation(aTrack);
   const G4String& rname = aTrack->GetVolume()->GetLogicalVolume()->GetRegion()->GetName();
   if (trkInfo != nullptr) {
 #ifdef EDM_ML_DEBUG
     trkInfo->Print();
 #endif
     if (rname == "FastTimerRegionSensBTL") {
       theID = trkInfo->mcTruthID();
       if (trkInfo->isExtSecondary() && !trkInfo->isInTrkFromBackscattering()) {
         theID = PSimHit::addTrackIdOffset(theID, k_idsecOffset);
       } else if (trkInfo->isInTrkFromBackscattering()) {
         theID = PSimHit::addTrackIdOffset(theID, k_idFromCaloOffset);
       } else if (trkInfo->isBTLlooper()) {
         theID = PSimHit::addTrackIdOffset(theID, k_idloopOffset);
       }
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("MtdSim") << "MtdSD: Track ID: " << aTrack->GetTrackID()
                                  << " BTL Track ID: " << trkInfo->mcTruthID() << ":" << theID;
 #endif
     } else if (rname == "FastTimerRegionSensETL") {
       theID = trkInfo->getIDonCaloSurface();
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("MtdSim") << "MtdSD: Track ID: " << aTrack->GetTrackID()
                                  << " ETL Track ID: " << trkInfo->mcTruthID() << ":" << theID;
 #endif
     } else {
       throw cms::Exception("MtdSDError") << "MtdSD called in incorrect region " << rname;
     }
   } else {
 #ifdef EDM_ML_DEBUG
     edm::LogWarning("MtdSim") << "MtdSD: Problem with primaryID **** set by force to TkID **** " << theID;
 #endif
   }
   return theID;
 }

◆ setDetUnitId()

uint32_t MtdSD::setDetUnitId ( const G4Step * aStep )

overridevirtual

Implements SensitiveDetector.

Definition at line 59 of file MtdSD.cc.

References getBaseNumber(), MTDNumberingScheme::getUnitID(), numberingScheme, and theBaseNumber.

                                                 {
   if (numberingScheme == nullptr) {
     return MTDDetId();
   } else {
     getBaseNumber(aStep);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("MtdSim") << "DetId = " << numberingScheme->getUnitID(theBaseNumber);
 #endif
     return numberingScheme->getUnitID(theBaseNumber);
   }
 }

◆ setNumberingScheme()

void MtdSD::setNumberingScheme ( MTDNumberingScheme * scheme )

private

Definition at line 71 of file MtdSD.cc.

References numberingScheme, and generator_cfi::scheme.

Referenced by MtdSD().

                                                          {
   if (scheme != nullptr) {
     edm::LogVerbatim("MtdSim") << "MtdSD: updates numbering scheme for " << GetName();
     if (numberingScheme)
       delete numberingScheme;
     numberingScheme = scheme;
   }
 }

Member Data Documentation

◆ energyCut

double MtdSD::energyCut

private

Definition at line 35 of file MtdSD.h.

Referenced by MtdSD().

◆ energyHistoryCut

double MtdSD::energyHistoryCut

private

Definition at line 36 of file MtdSD.h.

Referenced by MtdSD().

◆ isBTL

bool MtdSD::isBTL

private

Definition at line 43 of file MtdSD.h.

Referenced by MtdSD().

◆ isETL

bool MtdSD::isETL

private

Definition at line 44 of file MtdSD.h.

Referenced by MtdSD().

◆ k_idFromCaloOffset

constexpr unsigned int MtdSD::k_idFromCaloOffset = 3

staticprivate

Definition at line 33 of file MtdSD.h.

Referenced by getTrackID().

◆ k_idloopOffset

constexpr unsigned int MtdSD::k_idloopOffset = 2

staticprivate

Definition at line 32 of file MtdSD.h.

Referenced by getTrackID().

◆ k_idsecOffset

constexpr unsigned int MtdSD::k_idsecOffset = 1

staticprivate

Definition at line 31 of file MtdSD.h.

Referenced by getTrackID().

◆ numberingScheme

MTDNumberingScheme* MtdSD::numberingScheme

private

Definition at line 41 of file MtdSD.h.

Referenced by setDetUnitId(), and setNumberingScheme().

◆ theBaseNumber

MTDBaseNumber MtdSD::theBaseNumber

private

Definition at line 42 of file MtdSD.h.

Referenced by getBaseNumber(), and setDetUnitId().

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MtdSD()

◆ ~MtdSD()

Member Function Documentation

◆ getBaseNumber()

◆ getTrackID()

◆ setDetUnitId()

◆ setNumberingScheme()

Member Data Documentation

◆ energyCut

◆ energyHistoryCut

◆ isBTL

◆ isETL

◆ k_idFromCaloOffset

◆ k_idloopOffset

◆ k_idsecOffset

◆ numberingScheme

◆ theBaseNumber