#include <MuonSensitiveDetector.h>

Inheritance diagram for MuonSensitiveDetector:

Public Member Functions
void	clearHits () override

void	EndOfEvent (G4HCofThisEvent *) override

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

const MuonSlaveSD *	GetSlaveMuon () const

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

G4bool	ProcessHits (G4Step , G4TouchableHistory ) override

uint32_t	setDetUnitId (const G4Step *) override

	~MuonSensitiveDetector () override

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

void	slotForUpdate (const BeginOfEvent * iT)

virtual	~Observer ()

Protected Member Functions
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
Local3DPoint	cmsUnits (const Local3DPoint &v)

void	createHit (const G4Step *)

Local3DPoint	FinalStepPositionVsParent (const G4Step *currentStep, G4int levelsUp)

Local3DPoint	InitialStepPositionVsParent (const G4Step *currentStep, G4int levelsUp)

bool	newHit (const G4Step *)

void	saveHit ()

void	updateHit (const G4Step *)

Private Attributes
bool	allMuonsPersistent

MuonSubDetector *	detector

float	ePersistentCutGeV

MuonG4Numbering *	g4numbering

uint32_t	newDetUnitId

MuonSimHitNumberingScheme *	numbering

bool	printHits

MuonSlaveSD *	slaveMuon

uint32_t	theDetUnitId

G4ProcessTypeEnumerator *	theG4ProcessTypeEnumerator

UpdatablePSimHit *	theHit

const SimTrackManager *	theManager

SimHitPrinter *	thePrinter

const G4VPhysicalVolume *	thePV

const MuonFrameRotation *	theRotation

int	theTrackID

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

Detailed Description

implementation of SensitiveDetector for the muon detector; a MuonSlaveSD handles the interfacing to the database; numbering scheme are booked according to the detector name

Author: Arno Straessner, CERN arno..nosp@m.stra.nosp@m.essne.nosp@m.r@ce.nosp@m.rn.ch

Modification: 19/05/03. P.Arce Add SimTracks selection

Definition at line 37 of file MuonSensitiveDetector.h.

Constructor & Destructor Documentation

◆ MuonSensitiveDetector()

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

explicit

Definition at line 39 of file MuonSensitiveDetector.cc.

     : SensitiveTkDetector(name, es, clg, p),
       thePV(nullptr),
       theHit(nullptr),
       theDetUnitId(0),
       newDetUnitId(0),
       theTrackID(0),
       theManager(manager) {
   edm::ParameterSet m_MuonSD = p.getParameter<edm::ParameterSet>("MuonSD");
   ePersistentCutGeV = m_MuonSD.getParameter<double>("EnergyThresholdForPersistency") / CLHEP::GeV;  //Default 1. GeV
   allMuonsPersistent = m_MuonSD.getParameter<bool>("AllMuonsPersistent");
   printHits = m_MuonSD.getParameter<bool>("PrintHits");
   bool dd4hep = p.getParameter<bool>("g4GeometryDD4hepSource");
   //
   // Here simply create 1 MuonSlaveSD for the moment
   //
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("MuonSim") << "create MuonSubDetector " << name << " with dd4hep flag " << dd4hep;
 #endif
   detector = new MuonSubDetector(name);
  
   //The constants take time to calculate and are needed by many helpers
   edm::ESHandle<MuonOffsetMap> mom;
   es.get<IdealGeometryRecord>().get(mom);
   const MuonOffsetMap* offmap = (mom.isValid()) ? mom.product() : nullptr;
   edm::LogVerbatim("MuonSim") << "Finds the offset map at " << offmap;
   edm::ESHandle<MuonGeometryConstants> mdc;
   es.get<IdealGeometryRecord>().get(mdc);
   if (!mdc.isValid())
     throw cms::Exception("MuonSensitiveDetector") << "Cannot find MuonGeometryConstants\n";
   MuonGeometryConstants constants = *(mdc.product());
   G4String sdet = "unknown";
   if (detector->isEndcap()) {
     theRotation = new MuonEndcapFrameRotation();
     sdet = "Endcap";
   } else if (detector->isRPC()) {
     theRotation = new MuonRPCFrameRotation(constants, offmap, dd4hep);
     sdet = "RPC";
   } else if (detector->isGEM()) {
     theRotation = new MuonGEMFrameRotation(constants);
     sdet = "GEM";
   } else if (detector->isME0()) {
     theRotation = new MuonME0FrameRotation(constants);
     sdet = "ME0";
   } else {
     theRotation = new MuonFrameRotation();
   }
   slaveMuon = new MuonSlaveSD(detector, theManager);
   numbering = new MuonSimHitNumberingScheme(detector, constants);
   g4numbering = new MuonG4Numbering(constants, offmap, dd4hep);
  
   if (printHits) {
     thePrinter = new SimHitPrinter("HitPositionOSCAR.dat");
   }
  
   edm::LogVerbatim("MuonSim") << " of type " << sdet << " <" << GetName() << "> EnergyThresholdForPersistency(GeV) "
                               << ePersistentCutGeV / CLHEP::GeV << " allMuonsPersistent: " << allMuonsPersistent;
  
   theG4ProcessTypeEnumerator = new G4ProcessTypeEnumerator;
 }

References allMuonsPersistent, detector, ePersistentCutGeV, g4numbering, edm::EventSetup::get(), get, edm::ParameterSet::getParameter(), GeV, MuonSubDetector::isEndcap(), MuonSubDetector::isGEM(), MuonSubDetector::isME0(), MuonSubDetector::isRPC(), edm::ESHandleBase::isValid(), Skims_PA_cff::name, numbering, AlCaHLTBitMon_ParallelJobs::p, printHits, edm::ESHandle< T >::product(), slaveMuon, theG4ProcessTypeEnumerator, theManager, thePrinter, and theRotation.

◆ ~MuonSensitiveDetector()

MuonSensitiveDetector::~MuonSensitiveDetector ( )

override

Definition at line 104 of file MuonSensitiveDetector.cc.

                                               {
   delete g4numbering;
   delete numbering;
   delete slaveMuon;
   delete theRotation;
   delete detector;
   delete theG4ProcessTypeEnumerator;
 }

References detector, g4numbering, numbering, slaveMuon, theG4ProcessTypeEnumerator, and theRotation.

Member Function Documentation

◆ clearHits()

void MuonSensitiveDetector::clearHits ( )

overridevirtual

Implements SensitiveDetector.

Definition at line 121 of file MuonSensitiveDetector.cc.

                                       {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("MuonSim") << "MuonSensitiveDetector::clearHits";
 #endif
   slaveMuon->Initialize();
 }

References TrackingSlaveSD::Initialize(), and slaveMuon.

Referenced by update().

◆ cmsUnits()

Local3DPoint MuonSensitiveDetector::cmsUnits ( const Local3DPoint & v )

inlineprivate

Definition at line 58 of file MuonSensitiveDetector.h.

58 { return Local3DPoint(v.x() * 0.1, v.y() * 0.1, v.z() * 0.1); }

References findQualityFiles::v.

Referenced by createHit(), and updateHit().

◆ createHit()

void MuonSensitiveDetector::createHit ( const G4Step * aStep )

private

Definition at line 173 of file MuonSensitiveDetector.cc.

                                                          {
   Local3DPoint theEntryPoint;
   Local3DPoint theExitPoint;
  
   if (detector->isBarrel()) {
     // 1 levels up
     theEntryPoint = cmsUnits(theRotation->transformPoint(InitialStepPositionVsParent(aStep, 1), aStep));
     theExitPoint = cmsUnits(theRotation->transformPoint(FinalStepPositionVsParent(aStep, 1), aStep));
   } else if (detector->isEndcap()) {
     // save local z at current level
     theEntryPoint = theRotation->transformPoint(InitialStepPosition(aStep, LocalCoordinates), aStep);
     theExitPoint = theRotation->transformPoint(FinalStepPosition(aStep, LocalCoordinates), aStep);
     float zentry = theEntryPoint.z();
     float zexit = theExitPoint.z();
     // 4 levels up
     Local3DPoint tempEntry = theRotation->transformPoint(InitialStepPositionVsParent(aStep, 4), aStep);
     Local3DPoint tempExit = theRotation->transformPoint(FinalStepPositionVsParent(aStep, 4), aStep);
     // reset local z from z wrt deep-parent volume to z wrt low-level volume
     theEntryPoint = cmsUnits(Local3DPoint(tempEntry.x(), tempEntry.y(), zentry));
     theExitPoint = cmsUnits(Local3DPoint(tempExit.x(), tempExit.y(), zexit));
   } else {
     theEntryPoint = cmsUnits(theRotation->transformPoint(InitialStepPosition(aStep, LocalCoordinates), aStep));
     theExitPoint = cmsUnits(theRotation->transformPoint(FinalStepPosition(aStep, LocalCoordinates), aStep));
   }
  
   const G4Track* theTrack = aStep->GetTrack();
   const G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
  
   float thePabs = preStepPoint->GetMomentum().mag() / CLHEP::GeV;
   float theTof = preStepPoint->GetGlobalTime() / CLHEP::nanosecond;
   float theEnergyLoss = aStep->GetTotalEnergyDeposit() / CLHEP::GeV;
   int theParticleType = G4TrackToParticleID::particleID(theTrack);
  
   theDetUnitId = newDetUnitId;
   thePV = preStepPoint->GetPhysicalVolume();
   theTrackID = theTrack->GetTrackID();
  
   // convert momentum direction it to local frame
   const G4ThreeVector& gmd = preStepPoint->GetMomentumDirection();
   G4ThreeVector lmd = static_cast<const G4TouchableHistory*>(preStepPoint->GetTouchable())
                           ->GetHistory()
                           ->GetTopTransform()
                           .TransformAxis(gmd);
   Local3DPoint lnmd = ConvertToLocal3DPoint(lmd);
   lnmd = theRotation->transformPoint(lnmd, aStep);
   float theThetaAtEntry = lnmd.theta();
   float thePhiAtEntry = lnmd.phi();
  
   theHit = new UpdatablePSimHit(theEntryPoint,
                                 theExitPoint,
                                 thePabs,
                                 theTof,
                                 theEnergyLoss,
                                 theParticleType,
                                 theDetUnitId,
                                 theTrackID,
                                 theThetaAtEntry,
                                 thePhiAtEntry,
                                 theG4ProcessTypeEnumerator->processId(theTrack->GetCreatorProcess()));
  
   // Make track persistent
   int thePID = std::abs(theTrack->GetDefinition()->GetPDGEncoding());
   //---VI - in parameters cut in energy is declared but applied to momentum
   if (thePabs > ePersistentCutGeV || (thePID == 13 && allMuonsPersistent)) {
     TrackInformation* info = cmsTrackInformation(theTrack);
     info->storeTrack(true);
   }
  
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("MuonSim") << "=== NEW Muon hit for " << GetName() << " Edep(GeV)= " << theEnergyLoss << " "
                               << thePV->GetLogicalVolume()->GetName();
   const G4VProcess* p = aStep->GetPostStepPoint()->GetProcessDefinedStep();
   const G4VProcess* p2 = aStep->GetPreStepPoint()->GetProcessDefinedStep();
   G4String sss = "";
   if (p)
     sss += " POST PROCESS: " + p->GetProcessName();
   if (p2)
     sss += ";  PRE  PROCESS: " + p2->GetProcessName();
   if (!sss.empty())
     edm::LogVerbatim("MuonSim") << sss;
   edm::LogVerbatim("MuonSim") << " theta= " << theThetaAtEntry << " phi= " << thePhiAtEntry << " Pabs(GeV/c)  "
                               << thePabs << " Eloss(GeV)= " << theEnergyLoss << " Tof(ns)=  " << theTof
                               << " trackID= " << theTrackID << " detID= " << theDetUnitId << "\n Local:  entry "
                               << theEntryPoint << " exit " << theExitPoint << " delta "
                               << (theExitPoint - theEntryPoint) << "\n Global: entry "
                               << aStep->GetPreStepPoint()->GetPosition() << " exit "
                               << aStep->GetPostStepPoint()->GetPosition();
 #endif
 }

References funct::abs(), allMuonsPersistent, SensitiveDetector::cmsTrackInformation(), cmsUnits(), SensitiveDetector::ConvertToLocal3DPoint(), detector, ePersistentCutGeV, SensitiveDetector::FinalStepPosition(), FinalStepPositionVsParent(), GeV, info(), SensitiveDetector::InitialStepPosition(), InitialStepPositionVsParent(), MuonSubDetector::isBarrel(), MuonSubDetector::isEndcap(), SensitiveDetector::LocalCoordinates, PV3DBase< T, PVType, FrameType >::mag(), newDetUnitId, AlCaHLTBitMon_ParallelJobs::p, p2, G4TrackToParticleID::particleID(), PV3DBase< T, PVType, FrameType >::phi(), G4ProcessTypeEnumerator::processId(), theDetUnitId, theG4ProcessTypeEnumerator, theHit, thePV, theRotation, PV3DBase< T, PVType, FrameType >::theta(), theTrackID, MuonFrameRotation::transformPoint(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by ProcessHits().

◆ EndOfEvent()

void MuonSensitiveDetector::EndOfEvent ( G4HCofThisEvent * )

override

Definition at line 315 of file MuonSensitiveDetector.cc.

315 { saveHit(); }

References saveHit().

◆ fillHits()

void MuonSensitiveDetector::fillHits	(	edm::PSimHitContainer &	cc,
		const std::string &	hname
	)

overridevirtual

Implements SensitiveTkDetector.

Definition at line 317 of file MuonSensitiveDetector.cc.

                                                                                     {
   if (slaveMuon->name() == hname) {
     cc = slaveMuon->hits();
   }
 }

References TrackingSlaveSD::hits(), TrackingSlaveSD::name(), and slaveMuon.

◆ FinalStepPositionVsParent()

Local3DPoint MuonSensitiveDetector::FinalStepPositionVsParent	(	const G4Step *	currentStep,
		G4int	levelsUp
	)

private

Definition at line 336 of file MuonSensitiveDetector.cc.

                                                                                                        {
   const G4StepPoint* postStepPoint = currentStep->GetPostStepPoint();
   const G4StepPoint* preStepPoint = currentStep->GetPreStepPoint();
   const G4ThreeVector& globalCoordinates = postStepPoint->GetPosition();
  
   const G4TouchableHistory* theTouchable = (const G4TouchableHistory*)(preStepPoint->GetTouchable());
  
   G4int depth = theTouchable->GetHistory()->GetDepth();
   G4ThreeVector localCoordinates =
       theTouchable->GetHistory()->GetTransform(depth - levelsUp).TransformPoint(globalCoordinates);
  
   return ConvertToLocal3DPoint(localCoordinates);
 }

References SensitiveDetector::ConvertToLocal3DPoint(), and LEDCalibrationChannels::depth.

Referenced by createHit(), and updateHit().

◆ GetSlaveMuon()

const MuonSlaveSD* MuonSensitiveDetector::GetSlaveMuon ( ) const

inline

Definition at line 52 of file MuonSensitiveDetector.h.

52 { return slaveMuon; }

References slaveMuon.

◆ InitialStepPositionVsParent()

Local3DPoint MuonSensitiveDetector::InitialStepPositionVsParent	(	const G4Step *	currentStep,
		G4int	levelsUp
	)

private

Transform from local coordinates of a volume to local coordinates of a parent volume one or more levels up the volume hierarchy: e.g. levelsUp = 1 for immediate parent. This is done by moving from local_1 -> global -> local_2.

Definition at line 323 of file MuonSensitiveDetector.cc.

                                                                                                          {
   const G4StepPoint* preStepPoint = currentStep->GetPreStepPoint();
   const G4ThreeVector& globalCoordinates = preStepPoint->GetPosition();
  
   const G4TouchableHistory* theTouchable = (const G4TouchableHistory*)(preStepPoint->GetTouchable());
  
   G4int depth = theTouchable->GetHistory()->GetDepth();
   G4ThreeVector localCoordinates =
       theTouchable->GetHistory()->GetTransform(depth - levelsUp).TransformPoint(globalCoordinates);
  
   return ConvertToLocal3DPoint(localCoordinates);
 }

References SensitiveDetector::ConvertToLocal3DPoint(), and LEDCalibrationChannels::depth.

Referenced by createHit().

◆ newHit()

bool MuonSensitiveDetector::newHit ( const G4Step * aStep )

private

Definition at line 168 of file MuonSensitiveDetector.cc.

                                                       {
   return (!theHit || (aStep->GetTrack()->GetTrackID() != theTrackID) ||
           (aStep->GetPreStepPoint()->GetPhysicalVolume() != thePV) || newDetUnitId != theDetUnitId);
 }

References newDetUnitId, theDetUnitId, theHit, thePV, and theTrackID.

Referenced by ProcessHits().

◆ ProcessHits()

bool MuonSensitiveDetector::ProcessHits	(	G4Step *	aStep,
		G4TouchableHistory *	ROhist
	)

overridevirtual

Implements SensitiveDetector.

Definition at line 128 of file MuonSensitiveDetector.cc.

                                                                                  {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("MuonSim") << " MuonSensitiveDetector::ProcessHits " << InitialStepPosition(aStep, WorldCoordinates);
 #endif
  
   if (aStep->GetTotalEnergyDeposit() > 0.) {
     newDetUnitId = setDetUnitId(aStep);
  
     if (newHit(aStep)) {
       saveHit();
       createHit(aStep);
     } else {
       updateHit(aStep);
     }
     thePV = aStep->GetPreStepPoint()->GetPhysicalVolume();
     theTrackID = aStep->GetTrack()->GetTrackID();
     theDetUnitId = newDetUnitId;
   }
   return true;
 }

References createHit(), SensitiveDetector::InitialStepPosition(), newDetUnitId, newHit(), saveHit(), setDetUnitId(), theDetUnitId, thePV, theTrackID, updateHit(), and SensitiveDetector::WorldCoordinates.

◆ saveHit()

void MuonSensitiveDetector::saveHit ( )

private

Definition at line 301 of file MuonSensitiveDetector.cc.

                                     {
   if (theHit) {
     if (printHits) {
       thePrinter->startNewSimHit(detector->name());
       thePrinter->printId(theHit->detUnitId());
       thePrinter->printLocal(theHit->entryPoint(), theHit->exitPoint());
     }
     // hit is included into hit collection
     slaveMuon->processHits(*theHit);
     delete theHit;
     theHit = nullptr;
   }
 }

References detector, PSimHit::detUnitId(), PSimHit::entryPoint(), PSimHit::exitPoint(), MuonSubDetector::name(), printHits, SimHitPrinter::printId(), SimHitPrinter::printLocal(), TrackingSlaveSD::processHits(), slaveMuon, SimHitPrinter::startNewSimHit(), theHit, and thePrinter.

Referenced by EndOfEvent(), and ProcessHits().

◆ setDetUnitId()

uint32_t MuonSensitiveDetector::setDetUnitId ( const G4Step * aStep )

overridevirtual

Implements SensitiveDetector.

Definition at line 149 of file MuonSensitiveDetector.cc.

                                                                 {
   MuonBaseNumber num = g4numbering->PhysicalVolumeToBaseNumber(aStep);
  
 #ifdef EDM_ML_DEBUG
   std::stringstream MuonBaseNumber;
   MuonBaseNumber << "MuonNumbering :: number of levels = " << num.getLevels() << std::endl;
   MuonBaseNumber << "Level \t SuperNo \t BaseNo" << std::endl;
   for (int level = 1; level <= num.getLevels(); level++) {
     MuonBaseNumber << level << " \t " << num.getSuperNo(level) << " \t " << num.getBaseNo(level) << std::endl;
   }
   std::string MuonBaseNumbr = MuonBaseNumber.str();
  
   edm::LogVerbatim("MuonSim") << "MuonSensitiveDetector::setDetUnitId :: " << MuonBaseNumbr;
   edm::LogVerbatim("MuonSim") << "MuonSensitiveDetector::setDetUnitId :: MuonDetUnitId = "
                               << (numbering->baseNumberToUnitNumber(num));
 #endif
   return numbering->baseNumberToUnitNumber(num);
 }

References MuonSimHitNumberingScheme::baseNumberToUnitNumber(), g4numbering, personalPlayback::level, EgammaValidation_cff::num, numbering, MuonG4Numbering::PhysicalVolumeToBaseNumber(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by ProcessHits().

◆ update()

void MuonSensitiveDetector::update ( const BeginOfEvent * )

overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 113 of file MuonSensitiveDetector.cc.

                                                         {
   clearHits();
   //----- Initialize variables to check if two steps belong to same hit
   thePV = nullptr;
   theDetUnitId = 0;
   theTrackID = 0;
 }

References clearHits(), theDetUnitId, thePV, and theTrackID.

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

◆ updateHit()

void MuonSensitiveDetector::updateHit ( const G4Step * aStep )

private

Definition at line 263 of file MuonSensitiveDetector.cc.

                                                          {
   Local3DPoint theExitPoint;
  
   if (detector->isBarrel()) {
     theExitPoint = cmsUnits(theRotation->transformPoint(FinalStepPositionVsParent(aStep, 1), aStep));
   } else if (detector->isEndcap()) {
     // save local z at current level
     theExitPoint = theRotation->transformPoint(FinalStepPosition(aStep, LocalCoordinates), aStep);
     float zexit = theExitPoint.z();
     Local3DPoint tempExitPoint = theRotation->transformPoint(FinalStepPositionVsParent(aStep, 4), aStep);
     theExitPoint = cmsUnits(Local3DPoint(tempExitPoint.x(), tempExitPoint.y(), zexit));
   } else {
     theExitPoint = cmsUnits(theRotation->transformPoint(FinalStepPosition(aStep, LocalCoordinates), aStep));
   }
  
   float theEnergyLoss = aStep->GetTotalEnergyDeposit() / CLHEP::GeV;
  
   theHit->updateExitPoint(theExitPoint);
   theHit->addEnergyLoss(theEnergyLoss);
  
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("MuonSim") << "=== NEW Update muon hit for " << GetName() << " Edep(GeV)= " << theEnergyLoss << " "
                               << thePV->GetLogicalVolume()->GetName();
   const G4VProcess* p = aStep->GetPostStepPoint()->GetProcessDefinedStep();
   const G4VProcess* p2 = aStep->GetPreStepPoint()->GetProcessDefinedStep();
   G4String sss = "";
   if (p)
     sss += " POST PROCESS: " + p->GetProcessName();
   if (p2)
     sss += ";  PRE  PROCESS: " + p2->GetProcessName();
   if (!sss.empty())
     edm::LogVerbatim("MuonSim") << sss;
   edm::LogVerbatim("MuonSim") << " delEloss(GeV)= " << theEnergyLoss
                               << " Tof(ns)=  " << aStep->GetPreStepPoint()->GetGlobalTime() / CLHEP::nanosecond
                               << " trackID= " << theTrackID << " detID= " << theDetUnitId << " exit " << theExitPoint;
 #endif
 }

References UpdatablePSimHit::addEnergyLoss(), cmsUnits(), detector, SensitiveDetector::FinalStepPosition(), FinalStepPositionVsParent(), GeV, MuonSubDetector::isBarrel(), MuonSubDetector::isEndcap(), SensitiveDetector::LocalCoordinates, AlCaHLTBitMon_ParallelJobs::p, p2, theDetUnitId, theHit, thePV, theRotation, theTrackID, MuonFrameRotation::transformPoint(), UpdatablePSimHit::updateExitPoint(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by ProcessHits().

Member Data Documentation

◆ allMuonsPersistent

bool MuonSensitiveDetector::allMuonsPersistent

private

Definition at line 90 of file MuonSensitiveDetector.h.

Referenced by createHit(), and MuonSensitiveDetector().

◆ detector

MuonSubDetector* MuonSensitiveDetector::detector

private

Definition at line 62 of file MuonSensitiveDetector.h.

Referenced by createHit(), MuonSensitiveDetector(), saveHit(), updateHit(), and ~MuonSensitiveDetector().

◆ ePersistentCutGeV

float MuonSensitiveDetector::ePersistentCutGeV

private

Definition at line 89 of file MuonSensitiveDetector.h.

Referenced by createHit(), and MuonSensitiveDetector().

◆ g4numbering

MuonG4Numbering* MuonSensitiveDetector::g4numbering

private

Definition at line 64 of file MuonSensitiveDetector.h.

Referenced by MuonSensitiveDetector(), setDetUnitId(), and ~MuonSensitiveDetector().

◆ newDetUnitId

uint32_t MuonSensitiveDetector::newDetUnitId

private

Definition at line 82 of file MuonSensitiveDetector.h.

Referenced by createHit(), newHit(), and ProcessHits().

◆ numbering

MuonSimHitNumberingScheme* MuonSensitiveDetector::numbering

private

Definition at line 61 of file MuonSensitiveDetector.h.

Referenced by MuonSensitiveDetector(), setDetUnitId(), and ~MuonSensitiveDetector().

◆ printHits

bool MuonSensitiveDetector::printHits

private

Definition at line 85 of file MuonSensitiveDetector.h.

Referenced by ntuplePrintersDiff.TrackingParticlePrinter::__call__(), ntuplePrintersDiff.SeedPrinter::diff(), ntuplePrintersDiff.TrackPrinter::diff(), ntuplePrintersDiff.TrackingParticlePrinter::diff(), MuonSensitiveDetector(), ntuplePrintersDiff.SeedPrinter::printSeed(), ntuplePrintersDiff.TrackPrinter::printTrack(), and saveHit().

◆ slaveMuon

MuonSlaveSD* MuonSensitiveDetector::slaveMuon

private

Definition at line 60 of file MuonSensitiveDetector.h.

Referenced by clearHits(), fillHits(), GetSlaveMuon(), MuonSensitiveDetector(), saveHit(), and ~MuonSensitiveDetector().

◆ theDetUnitId

uint32_t MuonSensitiveDetector::theDetUnitId

private

Definition at line 81 of file MuonSensitiveDetector.h.

Referenced by createHit(), newHit(), ProcessHits(), update(), and updateHit().

◆ theG4ProcessTypeEnumerator

G4ProcessTypeEnumerator* MuonSensitiveDetector::theG4ProcessTypeEnumerator

private

Definition at line 92 of file MuonSensitiveDetector.h.

Referenced by createHit(), MuonSensitiveDetector(), and ~MuonSensitiveDetector().

◆ theHit

UpdatablePSimHit* MuonSensitiveDetector::theHit

private

Definition at line 80 of file MuonSensitiveDetector.h.

Referenced by createHit(), newHit(), saveHit(), and updateHit().

◆ theManager

const SimTrackManager* MuonSensitiveDetector::theManager

private

Definition at line 94 of file MuonSensitiveDetector.h.

Referenced by MuonSensitiveDetector().

◆ thePrinter

SimHitPrinter* MuonSensitiveDetector::thePrinter

private

Definition at line 86 of file MuonSensitiveDetector.h.

Referenced by MuonSensitiveDetector(), and saveHit().

◆ thePV

const G4VPhysicalVolume* MuonSensitiveDetector::thePV

private

Definition at line 79 of file MuonSensitiveDetector.h.

Referenced by createHit(), newHit(), ProcessHits(), update(), and updateHit().

◆ theRotation

const MuonFrameRotation* MuonSensitiveDetector::theRotation

private

Definition at line 63 of file MuonSensitiveDetector.h.

Referenced by createHit(), MuonSensitiveDetector(), updateHit(), and ~MuonSensitiveDetector().

◆ theTrackID

int MuonSensitiveDetector::theTrackID

private

Definition at line 83 of file MuonSensitiveDetector.h.

Referenced by createHit(), newHit(), ProcessHits(), update(), and updateHit().

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MuonSensitiveDetector()

◆ ~MuonSensitiveDetector()

Member Function Documentation

◆ clearHits()

◆ cmsUnits()

◆ createHit()

◆ EndOfEvent()

◆ fillHits()

◆ FinalStepPositionVsParent()

◆ GetSlaveMuon()

◆ InitialStepPositionVsParent()

◆ newHit()

◆ ProcessHits()

◆ saveHit()

◆ setDetUnitId()

◆ update()

◆ updateHit()

Member Data Documentation

◆ allMuonsPersistent

◆ detector

◆ ePersistentCutGeV

◆ g4numbering

◆ newDetUnitId

◆ numbering

◆ printHits

◆ slaveMuon

◆ theDetUnitId

◆ theG4ProcessTypeEnumerator

◆ theHit

◆ theManager

◆ thePrinter

◆ thePV

◆ theRotation

◆ theTrackID