Tool for constructing H62004 calib SDs for inactive material. More...

#include <H62004InactiveSDTool.h>

Inheritance diagram for LArG4::H62004InactiveSDTool:

Collaboration diagram for LArG4::H62004InactiveSDTool:

Public Member Functions
	H62004InactiveSDTool (const std::string &type, const std::string &name, const IInterface *parent)
	Constructor. More...

StatusCode	initialize () override final
	Initialize the tool. More...

StatusCode	Gather () override final
	Calls down to all the SDs to pack their hits into one collection. More...

virtual StatusCode	Gather (HitCollectionMap &) override

StatusCode	initializeSD () override
	Setup an SD in the current thread. More...

virtual StatusCode	SetupEvent () override
	Beginning of an athena event. More...

virtual StatusCode	SetupEvent (HitCollectionMap &) override

Protected Member Functions
std::unique_ptr< LArG4CalibSD >	makeOneSD (const std::string &name, ILArCalibCalculatorSvc *calc, const std::vector< std::string > &volumes) const
	Override helper method to create one SD. More...

StatusCode	assignSD (G4VSensitiveDetector *sd, const std::vector< std::string > &volumes) const
	Assign SD to a list of volumes. More...

G4VSensitiveDetector *	getSD ()
	Retrieve the current SD. More...

void	SetSensitiveDetector (G4LogicalVolume , G4VSensitiveDetector ) const
	Method stolen from G4VUserDetectorConstruction in G4 10.2. More...

Protected Attributes
Gaudi::Property< std::vector< std::string > >	m_volumeNames {this, "LogicalVolumeNames", {}}
	All the volumes to which this SD is assigned. More...

Gaudi::Property< std::vector< std::string > >	m_outputCollectionNames {this, "OutputCollectionNames", {}}
	Names of all output collections written out by this SD. More...

Gaudi::Property< bool >	m_noVolumes {this, "NoVolumes", false}
	This SensitiveDetector has no volumes associated with it. More...

Private Member Functions
StatusCode	initializeCalculators () override final
	Initialize Calculator Services. More...

G4VSensitiveDetector *	makeSD () const override final
	Create the SD wrapper for current worker thread. More...

void	setSD (G4VSensitiveDetector *)
	Set the current SD. More...

Static Private Member Functions
static bool	matchStrings (const char first, const char second)
	Match two strings with wildcard support. More...

Private Attributes
std::string	m_hitCollName {"LArCalibrationHitInactive"}
	Hit collection name. More...

ServiceHandle< ILArCalibCalculatorSvc >	m_emepiwcalc {this, "EMECPosIWCalibrationCalculator", "EMECPosInnerWheelCalibrationCalculator"}

ServiceHandle< ILArCalibCalculatorSvc >	m_heccalc {this, "HECWheelInactiveCalculator", "LocalCalibrationInactiveCalculator"}

ServiceHandle< ILArCalibCalculatorSvc >	m_fcal1calc {this, "FCAL1CalibCalculator", "LArFCAL1H62004CalibCalculator"}

ServiceHandle< ILArCalibCalculatorSvc >	m_fcal2calc {this, "FCAL2CalibCalculator", "LArFCAL2H62004CalibCalculator"}

G4VSensitiveDetector *	m_SD {}
	The sensitive detector to which this thing corresponds. More...

SD volumes
Gaudi::Property< std::vector< std::string > >	m_emecVolumes {this, "EMECVolumes"}

Gaudi::Property< std::vector< std::string > >	m_hecVolumes {this, "HECVolumes"}

Gaudi::Property< std::vector< std::string > >	m_fcal1Volumes {this, "FCAL1Volumes"}

Gaudi::Property< std::vector< std::string > >	m_fcal2Volumes {this, "FCAL2Volumes"}

Detailed Description

Tool for constructing H62004 calib SDs for inactive material.

Based on the previous LArG4H62004InactiveSDTool implementation.

This implementation uses the LAr SD wrapper design for managing multiple SDs when running multi-threaded. See ATLASSIM-2606 for discussions.

Definition at line 27 of file H62004InactiveSDTool.h.

Constructor & Destructor Documentation

◆ H62004InactiveSDTool()

LArG4::H62004InactiveSDTool::H62004InactiveSDTool	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent
	)

Constructor.

Member Function Documentation

◆ assignSD()

StatusCode SensitiveDetectorBase::assignSD	(	G4VSensitiveDetector *	sd,
		const std::vector< std::string > &	volumes
	)		const

protectedinherited

Assign SD to a list of volumes.

This method supports wild card matching

Definition at line 64 of file SensitiveDetectorBase.cxx.

 {
   // Propagate verbosity setting to the SD
   if(msgLvl(MSG::VERBOSE)) sd->SetVerboseLevel(10);
   else if(msgLvl(MSG::DEBUG)) sd->SetVerboseLevel(5);
  
   // Add the sensitive detector to the SD manager in G4 for SDs,
   // even if it has no volumes associated to it.
   auto sdMgr = G4SDManager::GetSDMpointer();
   sdMgr->AddNewDetector(sd);
  
   if(!volumes.empty()) {
     bool gotOne = false;
     auto logicalVolumeStore = G4LogicalVolumeStore::GetInstance();
     for(const auto& volumeName : volumes) {
       // Keep track of how many volumes we find with this name string.
       // We allow for multiple matches.
       int numFound = 0;
  
       // Find volumes with this name
       for(auto* logVol : *logicalVolumeStore) {
  
         ATH_MSG_VERBOSE("Check whether "<<logVol->GetName()<<" belongs to the set of sensitive detectors "<<volumeName);
         if( matchStrings( volumeName.data(), logVol->GetName() ) ){
           ++numFound;
           SetSensitiveDetector(logVol, sd);
         }
         
       }
       // Warn if no volumes were found
       if(numFound == 0) {
         ATH_MSG_WARNING("Volume " << volumeName <<
                         " not found in G4LogicalVolumeStore.");
       }
       else {
         ATH_MSG_VERBOSE("Found " << numFound << " copies of LV " << volumeName <<
                         "; SD " << sd->GetName() << " assigned.");
         gotOne = true;
       }
  
     }
  
     // Abort if we have failed to assign any volume
     if(!gotOne) {
       ATH_MSG_ERROR( "Failed to assign *any* volume to SD " << name() <<
                      " and expected at least one. Size of the volume store "<<G4LogicalVolumeStore::GetInstance()->size() );
       return StatusCode::FAILURE;
     }
   }
  
   return StatusCode::SUCCESS;
 }

◆ Gather() [1/2]

StatusCode LArG4::CalibSDTool::Gather ( )

finaloverridevirtualinherited

Calls down to all the SDs to pack their hits into one collection.

Reimplemented from SensitiveDetectorBase.

◆ Gather() [2/2]

virtual StatusCode SensitiveDetectorBase::Gather ( HitCollectionMap & )

inlineoverridevirtualinherited

Reimplemented in TileGeoG4SDTool, TileGeoG4CalibSDTool, SctSensorSDTool, PixelSensorSDTool, SctSensor_CTBTool, BCMSensorSDTool, BLMSensorSDTool, and TRTSensitiveDetectorTool.

Definition at line 79 of file SensitiveDetectorBase.h.

                                                         {
     return Gather();
   };

◆ getSD()

G4VSensitiveDetector * SensitiveDetectorBase::getSD ( )

protectedinherited

Retrieve the current SD.

In AthenaMT, this means the thread-local SD. Otherwise, it is simply the single SD.

Definition at line 118 of file SensitiveDetectorBase.cxx.

 {
 #ifdef G4MULTITHREADED
   // Get current thread-ID
   const auto tid = std::this_thread::get_id();
   // Retrieve it from the SD map
   auto sdPair = m_sdThreadMap.find(tid);
   if(sdPair == m_sdThreadMap.end()) return nullptr;
   return sdPair->second;
 #else
   return m_SD;
 #endif
 }

◆ initialize()

StatusCode LArG4::CalibSDTool::initialize ( )

finaloverrideinherited

Initialize the tool.

◆ initializeCalculators()

StatusCode LArG4::H62004InactiveSDTool::initializeCalculators ( )

finaloverrideprivatevirtual

Initialize Calculator Services.

Reimplemented from LArG4::CalibSDTool.

◆ initializeSD()

StatusCode SensitiveDetectorBase::initializeSD ( )

overrideinherited

Setup an SD in the current thread.

Separate from the AthAlgTool initialize() method because it needs to be called once per worker thread in AthenaMT. Don't confuse this with the G4 SD method Initialize which is called at the beginning of every G4 event.

Definition at line 25 of file SensitiveDetectorBase.cxx.

 {
   ATH_MSG_VERBOSE( name() << "::initializeSD()" );
  
   // Sanity check for volume configuration problems.
   // It would be better to have a more robust solution for this.
   if(m_volumeNames.empty() != m_noVolumes) {
     ATH_MSG_ERROR("Initializing SD from " << name() << ", NoVolumes = "
                   << (m_noVolumes? "true" : "false") << ", but LogicalVolumeNames = "
                   << m_volumeNames.value());
     return StatusCode::FAILURE;
   }
  
   // Make sure SD isn't already registered
   if(getSD())
     {
       ATH_MSG_ERROR("Trying to create a SD which already exists!");
       return StatusCode::FAILURE;
     }
  
   // Make the SD stored by this tool
   auto* sd = makeSD();
   if(!sd)
     {
       ATH_MSG_ERROR("Failed to create SD!");
       return StatusCode::FAILURE;
     }
   setSD(sd);
  
   // Assign the SD to our list of volumes
   ATH_CHECK( assignSD( getSD(), m_volumeNames.value() ) );
  
   ATH_MSG_DEBUG( "Initialized and added SD " << name() );
   return StatusCode::SUCCESS;
 }

◆ makeOneSD()

std::unique_ptr<LArG4CalibSD> LArG4::H62004CalibSDTool::makeOneSD	(	const std::string &	name,
		ILArCalibCalculatorSvc *	calc,
		const std::vector< std::string > &	volumes
	)		const

protectedinherited

Override helper method to create one SD.

◆ makeSD()

G4VSensitiveDetector* LArG4::H62004InactiveSDTool::makeSD ( ) const

finaloverrideprivate

Create the SD wrapper for current worker thread.

◆ matchStrings()

bool SensitiveDetectorBase::matchStrings	(	const char *	first,
		const char *	second
	)

staticprivateinherited

Match two strings with wildcard support.

Compares two strings character by character with optional * wildcard in the first argument

Definition at line 145 of file SensitiveDetectorBase.cxx.

 {
   // If we reach at the end of both strings, we are done
   if (*first == '\0' && *second == '\0')
     return true;
  
   // If there are consecutive '*' present in the first string
   // advance to the next character
   if(*first == '*' && *(first + 1) == '*')
     return matchStrings(first + 1, second);
  
   // Make sure that the characters after '*' are present in second string.
   if (*first == '*' && *(first + 1) != '\0' && *second == '\0')
     return false;
  
   // If the current characters of both strings match
   if (*first == *second)
     return matchStrings(first + 1, second + 1);
  
   // If there is *, then there are two possibilities
   // a) We consider current character of second string
   // b) We ignore current character of second string.
   if (*first == '*')
     return matchStrings(first + 1, second) || matchStrings(first, second + 1);
   return false;
 }

◆ setSD()

void SensitiveDetectorBase::setSD ( G4VSensitiveDetector * sd )

privateinherited

Set the current SD.

In hive, this gets assigned as the thread-local SD.

Definition at line 132 of file SensitiveDetectorBase.cxx.

 {
 #ifdef G4MULTITHREADED
   const auto tid = std::this_thread::get_id();
   ATH_MSG_DEBUG("Creating and registering SD " << sd << " in thread " << tid);
   m_sdThreadMap.insert( std::make_pair(tid, sd) );
 #else
   m_SD = sd;
 #endif
 }

◆ SetSensitiveDetector()

void SensitiveDetectorBase::SetSensitiveDetector	(	G4LogicalVolume *	logVol,
		G4VSensitiveDetector *	aSD
	)		const

protectedinherited

Method stolen from G4VUserDetectorConstruction in G4 10.2.

Definition at line 172 of file SensitiveDetectorBase.cxx.

 {
   // New Logic: allow for "multiple" SDs being attached to a single LV.
   // To do that we use a special proxy SD called G4MultiSensitiveDetector
  
   // Get existing SD if already set and check if it is of the special type
   G4VSensitiveDetector* originalSD = logVol->GetSensitiveDetector();
   if ( originalSD == nullptr )
     {
       logVol->SetSensitiveDetector(aSD);
     }
   else
     {
       G4MultiSensitiveDetector* msd = dynamic_cast<G4MultiSensitiveDetector*>(originalSD);
       if ( msd != nullptr )
         {
           msd->AddSD(aSD);
         }
       else
         {
           // Construct a unique name using the volume address
           std::stringstream ss;
           ss << static_cast<const void*>(logVol);
           const G4String msdname = "/MultiSD_" + logVol->GetName() + ss.str();
           //ATH_MSG_INFO("MultiSD name: " << msdname);
           msd = new G4MultiSensitiveDetector(msdname);
           // We need to register the proxy to have correct handling of IDs
           G4SDManager::GetSDMpointer()->AddNewDetector(msd);
           msd->AddSD(originalSD);
           msd->AddSD(aSD);
           logVol->SetSensitiveDetector(msd);
         }
     }
 }

◆ SetupEvent() [1/2]

virtual StatusCode SensitiveDetectorBase::SetupEvent ( )

inlineoverridevirtualinherited

Beginning of an athena event.

This is where collection initialization should happen. If we are using a WriteHandle, then this could be empty.

Reimplemented in MuonWallSDTool, AFP_SensitiveDetectorTool, AFP_SiDSensitiveDetectorTool, AFP_TDSensitiveDetectorTool, CaloCellContainerSDTool, and ALFA_SensitiveDetectorTool.

Definition at line 70 of file SensitiveDetectorBase.h.

70 { return StatusCode::SUCCESS; }

◆ SetupEvent() [2/2]

virtual StatusCode SensitiveDetectorBase::SetupEvent ( HitCollectionMap & )

inlineoverridevirtualinherited

Reimplemented in TileGeoG4SDTool, TileGeoG4CalibSDTool, SctSensorSDTool, PixelSensorSDTool, SctSensor_CTBTool, BCMSensorSDTool, BLMSensorSDTool, and TRTSensitiveDetectorTool.

Definition at line 71 of file SensitiveDetectorBase.h.

                                                             {
     return SetupEvent();
   }

Member Data Documentation

◆ m_caloDmID

const CaloDM_ID* LArG4::CalibSDTool::m_caloDmID {nullptr}

protectedinherited

Definition at line 83 of file CalibSDTool.h.

◆ m_doPID

Gaudi::Property<G4bool> LArG4::CalibSDTool::m_doPID {this, "ParticleID", false}

protectedinherited

Are we set up to run with PID hits?

Definition at line 74 of file CalibSDTool.h.

◆ m_emecVolumes

Gaudi::Property<std::vector<std::string> > LArG4::H62004InactiveSDTool::m_emecVolumes {this, "EMECVolumes"}

private

Definition at line 53 of file H62004InactiveSDTool.h.

◆ m_emepiwcalc

ServiceHandle<ILArCalibCalculatorSvc> LArG4::H62004InactiveSDTool::m_emepiwcalc {this, "EMECPosIWCalibrationCalculator", "EMECPosInnerWheelCalibrationCalculator"}

private

Definition at line 46 of file H62004InactiveSDTool.h.

◆ m_fcal1calc

ServiceHandle<ILArCalibCalculatorSvc> LArG4::H62004InactiveSDTool::m_fcal1calc {this, "FCAL1CalibCalculator", "LArFCAL1H62004CalibCalculator"}

private

Definition at line 48 of file H62004InactiveSDTool.h.

◆ m_fcal1Volumes

Gaudi::Property<std::vector<std::string> > LArG4::H62004InactiveSDTool::m_fcal1Volumes {this, "FCAL1Volumes"}

private

Definition at line 55 of file H62004InactiveSDTool.h.

◆ m_fcal2calc

ServiceHandle<ILArCalibCalculatorSvc> LArG4::H62004InactiveSDTool::m_fcal2calc {this, "FCAL2CalibCalculator", "LArFCAL2H62004CalibCalculator"}

private

Definition at line 49 of file H62004InactiveSDTool.h.

◆ m_fcal2Volumes

Gaudi::Property<std::vector<std::string> > LArG4::H62004InactiveSDTool::m_fcal2Volumes {this, "FCAL2Volumes"}

private

Definition at line 56 of file H62004InactiveSDTool.h.

◆ m_heccalc

ServiceHandle<ILArCalibCalculatorSvc> LArG4::H62004InactiveSDTool::m_heccalc {this, "HECWheelInactiveCalculator", "LocalCalibrationInactiveCalculator"}

private

Definition at line 47 of file H62004InactiveSDTool.h.

◆ m_hecVolumes

Gaudi::Property<std::vector<std::string> > LArG4::H62004InactiveSDTool::m_hecVolumes {this, "HECVolumes"}

private

Definition at line 54 of file H62004InactiveSDTool.h.

◆ m_hitCollName

std::string LArG4::H62004InactiveSDTool::m_hitCollName {"LArCalibrationHitInactive"}

private

Hit collection name.

Definition at line 44 of file H62004InactiveSDTool.h.

◆ m_id_helper

const AtlasDetectorID* LArG4::CalibSDTool::m_id_helper {nullptr}

protectedinherited

Definition at line 84 of file CalibSDTool.h.

◆ m_larEmID

const LArEM_ID* LArG4::CalibSDTool::m_larEmID {nullptr}

protectedinherited

Definition at line 80 of file CalibSDTool.h.

◆ m_larFcalID

const LArFCAL_ID* LArG4::CalibSDTool::m_larFcalID {nullptr}

protectedinherited

Definition at line 81 of file CalibSDTool.h.

◆ m_larHecID

const LArHEC_ID* LArG4::CalibSDTool::m_larHecID {nullptr}

protectedinherited

Definition at line 82 of file CalibSDTool.h.

◆ m_noVolumes

Gaudi::Property<bool> SensitiveDetectorBase::m_noVolumes {this, "NoVolumes", false}

protectedinherited

This SensitiveDetector has no volumes associated with it.

Definition at line 101 of file SensitiveDetectorBase.h.

◆ m_outputCollectionNames

Gaudi::Property<std::vector<std::string> > SensitiveDetectorBase::m_outputCollectionNames {this, "OutputCollectionNames", {}}

protectedinherited

Names of all output collections written out by this SD.

Definition at line 98 of file SensitiveDetectorBase.h.

◆ m_SD

G4VSensitiveDetector* SensitiveDetectorBase::m_SD {}

privateinherited

The sensitive detector to which this thing corresponds.

Definition at line 125 of file SensitiveDetectorBase.h.

◆ m_volumeNames

Gaudi::Property<std::vector<std::string> > SensitiveDetectorBase::m_volumeNames {this, "LogicalVolumeNames", {}}

protectedinherited

All the volumes to which this SD is assigned.

Definition at line 96 of file SensitiveDetectorBase.h.

The documentation for this class was generated from the following file:

H62004InactiveSDTool.h

Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Static Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ H62004InactiveSDTool()

Member Function Documentation

◆ assignSD()

◆ Gather() [1/2]

◆ Gather() [2/2]

◆ getSD()

◆ initialize()

◆ initializeCalculators()

◆ initializeSD()

◆ makeOneSD()

◆ makeSD()

◆ matchStrings()

◆ setSD()

◆ SetSensitiveDetector()

◆ SetupEvent() [1/2]

◆ SetupEvent() [2/2]

Member Data Documentation

◆ m_caloDmID

◆ m_doPID

◆ m_emecVolumes

◆ m_emepiwcalc

◆ m_fcal1calc

◆ m_fcal1Volumes

◆ m_fcal2calc

◆ m_fcal2Volumes

◆ m_heccalc

◆ m_hecVolumes

◆ m_hitCollName

◆ m_id_helper

◆ m_larEmID

◆ m_larFcalID

◆ m_larHecID

◆ m_noVolumes

◆ m_outputCollectionNames

◆ m_SD

◆ m_volumeNames