CaloCellContainerSDTool Class Reference

#include <CaloCellContainerSDTool.h>

Inheritance diagram for CaloCellContainerSDTool:

Collaboration diagram for CaloCellContainerSDTool:

Public Member Functions
	CaloCellContainerSDTool (const std::string &type, const std::string &name, const IInterface *parent)
	~CaloCellContainerSDTool ()
StatusCode	initialize () override final
StatusCode	SetupEvent (HitCollectionMap &) override final
	Beginning of an athena event.
StatusCode	Gather (HitCollectionMap &) override final
	End of an athena event.
StatusCode	initializeSD () override
	Setup an SD in the current thread.

Protected Member Functions
G4VSensitiveDetector *	makeSD () const override final
StatusCode	assignSD (std::unique_ptr< G4VSensitiveDetector > sd, const std::vector< std::string > &volumes) const
	Assign SD to a list of volumes.
void	SetSensitiveDetector (G4LogicalVolume *, G4VSensitiveDetector *) const
	Method stolen from G4VUserDetectorConstruction in G4 10.2.

Protected Attributes
PublicToolHandle< ICaloCellMakerTool >	m_EmptyCellBuilderTool
PublicToolHandle< ICaloCellMakerTool >	m_FastHitConvertTool
Gaudi::Property< std::vector< std::string > >	m_volumeNames {this, "LogicalVolumeNames", {}}
	All the volumes to which this SD is assigned.
Gaudi::Property< std::vector< std::string > >	m_outputCollectionNames {this, "OutputCollectionNames", {}}
	Names of all output collections written out by this SD.
Gaudi::Property< bool >	m_noVolumes {this, "NoVolumes", false}
	This SensitiveDetector has no volumes associated with it.

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

Detailed Description

Definition at line 16 of file CaloCellContainerSDTool.h.

Constructor & Destructor Documentation

CaloCellContainerSDTool()

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

Definition at line 18 of file CaloCellContainerSDTool.cxx.

  : SensitiveDetectorBase( type , name , parent ),
    m_EmptyCellBuilderTool("EmptyCellBuilderTool/EmptyCellBuilderTool"),
    m_FastHitConvertTool ("FastHitConvertTool/FastHitConvertTool")
{
   declareProperty ("EmptyCellBuilderTool", m_EmptyCellBuilderTool, "EmptyCellBuilderTool");
   declareProperty ("FastHitConvertTool", m_FastHitConvertTool, "FastHitConvertTool");
}

~CaloCellContainerSDTool()

CaloCellContainerSDTool::~CaloCellContainerSDTool ( )

inline

Definition at line 22 of file CaloCellContainerSDTool.h.

{}

Member Function Documentation

assignSD()

StatusCode SensitiveDetectorBase::assignSD ( std::unique_ptr< 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 55 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();
  auto sdPtr = sd.get();
  // SDManager is now the SD owner
  //for later use
  auto sdName = sd->GetName();
  sdMgr->AddNewDetector(sd.release());
 
  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, sdPtr);
        }
        
      }
      // 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 " << sdName << " 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()

StatusCode CaloCellContainerSDTool::Gather ( HitCollectionMap & hitCollections )

finaloverridevirtual

End of an athena event.

Reimplemented from SensitiveDetectorBase.

Definition at line 56 of file CaloCellContainerSDTool.cxx.

{
  // Get current event context
  const EventContext& ctx = Gaudi::Hive::currentContext();
 
  const std::string& collectionName = m_outputCollectionNames[0];
  std::unique_ptr<CaloCellContainerBuilder> builder =
    hitCollections.Extract<CaloCellContainerBuilder>(collectionName);
  if (!builder || !builder->container) {
    ATH_MSG_ERROR("Gather: Failed to retrieve CaloCellContainerBuilder.");
    return StatusCode::FAILURE;
  }
 
  // Update the calo iterators of the calo cell container.
  builder->container->updateCaloIterators();
 
  // Convert FastCaloSim hits into HitCollections, taking into account sampling fractions.
  if(m_FastHitConvertTool->process(builder->container.get(), ctx).isFailure()){
    ATH_MSG_ERROR("Gather: Failed to process calo cell container with the fast hit convert tool.");
    return StatusCode::FAILURE;
  }
 
  SG::WriteHandle<CaloCellContainer> handle{collectionName, ctx};
  ATH_CHECK(handle.record(std::move(builder->container)));
  return StatusCode::SUCCESS;
}

initialize()

StatusCode CaloCellContainerSDTool::initialize ( )

finaloverride

Definition at line 27 of file CaloCellContainerSDTool.cxx.

{
  ATH_CHECK(m_EmptyCellBuilderTool.retrieve());
  ATH_CHECK(m_FastHitConvertTool.retrieve());
  return StatusCode::SUCCESS;
}

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 the SD stored by this tool
  auto sd = std::unique_ptr<G4VSensitiveDetector>(makeSD());
  if(!sd)
    {
      ATH_MSG_ERROR("Failed to create SD!");
      return StatusCode::FAILURE;
    }
  // Assign the SD to our list of volumes
  ATH_CHECK( assignSD( std::move(sd), m_volumeNames.value() ) );
 
  ATH_MSG_DEBUG( "Initialized and added SD " << name() );
  return StatusCode::SUCCESS;
}

makeSD()

G4VSensitiveDetector * CaloCellContainerSDTool::makeSD ( ) const

finaloverrideprotected

Definition at line 83 of file CaloCellContainerSDTool.cxx.

{
  ATH_MSG_DEBUG( "Initializing SD" );
 
  // Create a fresh SD
  return new CaloCellContainerSD(name(), m_outputCollectionNames[0]);
}

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 115 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;
}

SetSensitiveDetector()

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

protectedinherited

Method stolen from G4VUserDetectorConstruction in G4 10.2.

Definition at line 142 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();
          msd = new G4MultiSensitiveDetector(std::move(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()

StatusCode CaloCellContainerSDTool::SetupEvent ( HitCollectionMap & hitCollections )

finaloverridevirtual

Beginning of an athena event.

This is where collection initialization should happen.

Reimplemented from SensitiveDetectorBase.

Definition at line 34 of file CaloCellContainerSDTool.cxx.

{
  // Get current event context
  const EventContext& ctx = Gaudi::Hive::currentContext();
 
  const std::string& collectionName = m_outputCollectionNames[0];
  hitCollections.Emplace<CaloCellContainerBuilder>(collectionName);
  auto* builder = hitCollections.Find<CaloCellContainerBuilder>(collectionName);
  if (!builder) {
    ATH_MSG_ERROR("SetupEvent: Failed to create CaloCellContainerBuilder.");
    return StatusCode::FAILURE;
  }
 
  // Initialize cell container with empty cells.
  if(m_EmptyCellBuilderTool->process(builder->container.get(), ctx).isFailure()){
    ATH_MSG_ERROR("SetupEvent: Failed to process calo cell container with the empty cell builder tool.");
    return StatusCode::FAILURE;
  }
 
  return StatusCode::SUCCESS;
}

Member Data Documentation

m_EmptyCellBuilderTool

PublicToolHandle<ICaloCellMakerTool> CaloCellContainerSDTool::m_EmptyCellBuilderTool

protected

Definition at line 33 of file CaloCellContainerSDTool.h.

m_FastHitConvertTool

PublicToolHandle<ICaloCellMakerTool> CaloCellContainerSDTool::m_FastHitConvertTool

protected

Definition at line 34 of file CaloCellContainerSDTool.h.

m_noVolumes

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

protectedinherited

This SensitiveDetector has no volumes associated with it.

Definition at line 83 of file SensitiveDetectorBase.h.

{this, "NoVolumes", false};

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 80 of file SensitiveDetectorBase.h.

{this, "OutputCollectionNames", {}};

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 78 of file SensitiveDetectorBase.h.

{this, "LogicalVolumeNames", {}};

---

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

• CaloCellContainerSDTool.h
• CaloCellContainerSDTool.cxx