#include <MuonG4Numbering.h>

Public Member Functions
	MuonG4Numbering (const MuonGeometryConstants &muonConstants, const MuonOffsetMap *offMap, bool dd4hep)

MuonBaseNumber	PhysicalVolumeToBaseNumber (const G4Step *aStep)

	~MuonG4Numbering ()

Private Member Functions
const bool	copyNoRelevant (const int)

const int	getCopyNoBaseNo (const int)

const int	getCopyNoLevel (const int)

const int	getCopyNoSuperNo (const int)

Private Attributes
const bool	dd4hep_

const MuonOffsetMap *	offMap_

int	theBasePart

int	theLevelPart

int	theStartCopyNo

int	theSuperPart

Detailed Description

class to handle the conversion to MuonBaseNumber from tree of G4 copy numbers extracted from G4VPhysicalVolume; needs DDD compatible format of G4 copy numbers; the format itself may in a later version also be defined by DDD constants

in the xml muon constant section one has to define level, super and base constants (eg. 1000,100,1) and the start value of the copy numbers (0 or 1)

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

Modification:

Definition at line 26 of file MuonG4Numbering.h.

Constructor & Destructor Documentation

◆ MuonG4Numbering()

MuonG4Numbering::MuonG4Numbering	(	const MuonGeometryConstants &	muonConstants,
		const MuonOffsetMap *	offMap,
		bool	dd4hep
	)

Definition at line 16 of file MuonG4Numbering.cc.

     : offMap_(offMap), dd4hep_(dd4hep) {
   theLevelPart = muonConstants.getValue("level");
   theSuperPart = muonConstants.getValue("super");
   theBasePart = muonConstants.getValue("base");
   theStartCopyNo = muonConstants.getValue("xml_starts_with_copyno");
  
   // some consistency checks
  
   if (theBasePart != 1) {
     edm::LogVerbatim("MuonSim") << "MuonGeometryNumbering finds unusual base constant:" << theBasePart;
   }
   if (theSuperPart < 100) {
     edm::LogVerbatim("MuonSim") << "MuonGeometryNumbering finds unusual super constant:" << theSuperPart;
   }
   if (theLevelPart < 10 * theSuperPart) {
     edm::LogVerbatim("MuonSim") << "MuonGeometryNumbering finds unusual level constant:" << theLevelPart;
   }
   if ((theStartCopyNo != 0) && (theStartCopyNo != 1)) {
     std::cout << "MuonGeometryNumbering finds unusual start value for copy numbers:" << theStartCopyNo << std::endl;
   }
  
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("MuonSim") << "StartCopyNo = " << theStartCopyNo;
   edm::LogVerbatim("MuonSim") << "MuonG4Numbering configured with"
                               << "Level = " << theLevelPart << " Super = " << theSuperPart << " Base = " << theBasePart
                               << " StartCopyNo = " << theStartCopyNo;
   edm::LogVerbatim("MuonSim") << "dd4hep flag set to " << dd4hep_ << " and offsetmap at " << offMap_;
 #endif
 }

References gather_cfg::cout, dd4hep_, MuonGeometryConstants::getValue(), offMap_, theBasePart, theLevelPart, theStartCopyNo, and theSuperPart.

◆ ~MuonG4Numbering()

MuonG4Numbering::~MuonG4Numbering ( )

inline

Definition at line 29 of file MuonG4Numbering.h.

29 {};

Member Function Documentation

◆ copyNoRelevant()

const bool MuonG4Numbering::copyNoRelevant ( const int copyno )

private

Definition at line 90 of file MuonG4Numbering.cc.

90 { return (copyno / theLevelPart) > 0; }

References theLevelPart.

Referenced by PhysicalVolumeToBaseNumber().

◆ getCopyNoBaseNo()

const int MuonG4Numbering::getCopyNoBaseNo ( const int copyno )

private

Definition at line 88 of file MuonG4Numbering.cc.

88 { return copyno % theSuperPart; }

References theSuperPart.

Referenced by PhysicalVolumeToBaseNumber().

◆ getCopyNoLevel()

const int MuonG4Numbering::getCopyNoLevel ( const int copyno )

private

Definition at line 84 of file MuonG4Numbering.cc.

84 { return copyno / theLevelPart; }

References theLevelPart.

Referenced by PhysicalVolumeToBaseNumber().

◆ getCopyNoSuperNo()

const int MuonG4Numbering::getCopyNoSuperNo ( const int copyno )

private

Definition at line 86 of file MuonG4Numbering.cc.

86 { return (copyno % theLevelPart) / theSuperPart; }

References theLevelPart, and theSuperPart.

Referenced by PhysicalVolumeToBaseNumber().

◆ PhysicalVolumeToBaseNumber()

MuonBaseNumber MuonG4Numbering::PhysicalVolumeToBaseNumber ( const G4Step * aStep )

Definition at line 47 of file MuonG4Numbering.cc.

                                                                               {
   MuonBaseNumber num;
   const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
  
   for (int ii = 0; ii < touch->GetHistoryDepth(); ii++) {
     G4VPhysicalVolume* vol = touch->GetVolume(ii);
     int copyno = vol->GetCopyNo();
     int extra(0);
     if (dd4hep_ && (offMap_ != nullptr)) {
       std::string namx = static_cast<std::string>(dd4hep::dd::noNamespace(vol->GetName()));
       std::size_t last = namx.rfind('_');
       std::string name = ((last == std::string::npos) ? namx : (namx.substr(0, last)));
       auto itr = offMap_->muonMap_.find(name);
       if (itr != offMap_->muonMap_.end())
         extra = (itr->second).first + (itr->second).second;
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("MuonSim") << "MuonG4Numbering: " << namx << ":" << name << " iterator "
                                   << (itr != offMap_->muonMap_.end()) << " Extra " << extra;
 #endif
     }
     copyno += extra;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("MuonSim") << "MuonG4Numbering: " << vol->GetName() << " " << copyno << " Split "
                                 << copyNoRelevant(copyno) << ":" << theLevelPart << ":" << theSuperPart << " ";
 #endif
     if (copyNoRelevant(copyno)) {
       num.addBase(getCopyNoLevel(copyno), getCopyNoSuperNo(copyno), getCopyNoBaseNo(copyno) - theStartCopyNo);
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("MuonSim") << " NoLevel " << getCopyNoLevel(copyno) << " Super " << getCopyNoSuperNo(copyno)
                                   << " Base " << getCopyNoBaseNo(copyno) << " Start " << theStartCopyNo;
 #endif
     }
   }
  
   return num;
 }