DerivationFramework::EventInfoPixelDecorator Class Reference

#include <EventInfoPixelDecorator.h>

Inheritance diagram for DerivationFramework::EventInfoPixelDecorator:

Collaboration diagram for DerivationFramework::EventInfoPixelDecorator:

Public Member Functions
	EventInfoPixelDecorator (const std::string &type, const std::string &name, const IInterface *parent)
StatusCode	initialize ()
StatusCode	finalize ()
virtual StatusCode	addBranches () const

Private Types
enum	EIntDecor {   kperModuleMultiplicity, klayer, keta_module, kphi_module,   kNIntDecor }

Private Attributes
Gaudi::Property< std::string >	m_selectionString { this, "SelectionString", "" , "" }
std::atomic< unsigned int >	m_ntot {}
std::atomic< unsigned int >	m_npass {}
SG::ReadHandleKey< xAOD::EventInfo >	m_eventInfoKey { this, "EventInfoKey", "EventInfo", ""}
Gaudi::Property< std::string >	m_decorationPrefix { this, "DecorationPrefix", "", ""}
SG::ReadHandleKey< xAOD::TrackMeasurementValidationContainer >	m_pixelKey { this, "TrackMeasurementValidationKey", "PixelClusters", ""}
std::vector< SG::WriteDecorHandleKey< xAOD::EventInfo > >	m_intDecorKeys

Detailed Description

Definition at line 26 of file EventInfoPixelDecorator.h.

Member Enumeration Documentation

EIntDecor

enum DerivationFramework::EventInfoPixelDecorator::EIntDecor

private

Enumerator
kperModuleMultiplicity
klayer
keta_module
kphi_module
kNIntDecor

Definition at line 48 of file EventInfoPixelDecorator.h.

                     {kperModuleMultiplicity,
                      klayer,
                      keta_module,
                      kphi_module,
                      kNIntDecor};

Constructor & Destructor Documentation

EventInfoPixelDecorator()

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

Definition at line 23 of file EventInfoPixelDecorator.cxx.

                                : 
    base_class(type,name,parent)
  {
  }

Member Function Documentation

addBranches()

StatusCode DerivationFramework::EventInfoPixelDecorator::addBranches ( ) const

virtual

Definition at line 79 of file EventInfoPixelDecorator.cxx.

  {
    ATH_MSG_DEBUG("Adding per-module information to EventInfo");
    const EventContext& ctx = Gaudi::Hive::currentContext();
 
    SG::ReadHandle<xAOD::EventInfo> eventInfo(m_eventInfoKey, ctx);
    CHECK( eventInfo.isValid() ? StatusCode::SUCCESS : StatusCode::FAILURE);
 
    // Get the cluster container
    SG::ReadHandle<xAOD::TrackMeasurementValidationContainer> clusters(m_pixelKey, ctx);
    if (!clusters.isValid()) {
       ATH_MSG_ERROR ("Couldn't retrieve TrackMeasurementValidationContainer with key" << m_pixelKey.key());
      return StatusCode::FAILURE;
    }
 
    // Check the event contains clusters
    unsigned int nClusters = clusters->size();
    if (nClusters==0) return StatusCode::SUCCESS;
 
    // Set up a mask with the same entries as the full collection
    std::vector<bool> mask;
    mask.assign(nClusters,false); // default: don't keep anything
    m_ntot += nClusters;
 
    // Execute the text parser and update the mask
    if (m_parser) {
        std::vector<int> entries =  m_parser->evaluateAsVector();
        unsigned int nEntries = entries.size();
        // check the sizes are compatible
        if (nClusters != nEntries ) {
            ATH_MSG_ERROR("Sizes incompatible! Are you sure your selection string used the same collection as the input collection??");
            return StatusCode::FAILURE;
        } else {
            // set mask
            for (unsigned int i=0; i<nClusters; ++i) if (entries[i]==1) mask[i]=true;
        }
    }
    // Count the mask
    unsigned int n_pass=0;
    for (unsigned int i=0; i<nClusters; ++i) {
        if (mask[i]) ++n_pass;
    }
    m_npass += n_pass;
 
    // fill the per-module information
    std::array<std::vector<int>,kNIntDecor> vec;
 
    unsigned int i=0;
    std::vector<int> keys;
    for (const xAOD::TrackMeasurementValidation* clusIt : *clusters) {
 
      if( !mask[++i] ) continue;
 
      static const SG::AuxElement::Accessor<int> acc_layer("layer");
      int clus_layer = acc_layer(*clusIt);
      static const SG::AuxElement::Accessor<int> acc_eta_module("eta_module");
      int clus_eta_module = acc_eta_module(*clusIt);
      static const SG::AuxElement::Accessor<int> acc_phi_module("phi_module");
      int clus_phi_module = acc_phi_module(*clusIt);
 
      int key = clus_layer * 10000 + clus_eta_module * 100 + clus_phi_module;
      unsigned int index = 9999;
      for( unsigned int j = 0; j < keys.size() ; j++ )
        if( key == keys[j] ){
          index = j;
          break;
        }
      if( index == 9999 ){
        keys.push_back( key );
        vec[klayer].push_back( clus_layer );
        vec[keta_module].push_back( clus_eta_module );
        vec[kphi_module].push_back( clus_phi_module );
        vec[kperModuleMultiplicity].push_back( 1 );
      } else {
        assert(index < vec[kperModuleMultiplicity].size());
        vec[kperModuleMultiplicity][index]++;
      }
    }
 
    std::vector<SG::WriteDecorHandle<xAOD::EventInfo,std::vector<int> > >   int_decor_handles(createDecorators<xAOD::EventInfo,std::vector<int> >(m_intDecorKeys,ctx));
    assert(int_decor_handles.size() == kNIntDecor);
    for(unsigned int decorate_i=0; decorate_i<int_decor_handles.size(); ++decorate_i) {
       int_decor_handles[decorate_i](*eventInfo) = std::move(vec[decorate_i]);
    }
 
    return StatusCode::SUCCESS;
  }  

finalize()

StatusCode DerivationFramework::EventInfoPixelDecorator::finalize

(

)

Definition at line 71 of file EventInfoPixelDecorator.cxx.

  {
    ATH_MSG_VERBOSE("finalize() ...");
    ATH_MSG_INFO("Processed "<< m_ntot <<" measurements, "<< m_npass<< " were used for the per-module summary");
    ATH_CHECK( finalizeParser() );
    return StatusCode::SUCCESS;
  }

initialize()

StatusCode DerivationFramework::EventInfoPixelDecorator::initialize

(

)

Definition at line 30 of file EventInfoPixelDecorator.cxx.

  {
    ATH_MSG_VERBOSE("initialize() ...");
 
    if (m_selectionString.empty()) {
        ATH_MSG_FATAL("No inner detector track selection string provided!");
        return StatusCode::FAILURE;
    } else {ATH_MSG_INFO("Selection string for the per-module pixel cluster decoration: " << m_selectionString);}
 
    if (m_decorationPrefix.empty()) {
      ATH_MSG_WARNING("No decoration prefix name provided for the output of EventInfoPixelDecorator!");
    }
    
    if (m_eventInfoKey.key().empty()) {
      ATH_MSG_ERROR("No collection provided for EventInfoPixelDecorator!");
      return StatusCode::FAILURE;
    }
    ATH_CHECK(m_eventInfoKey.initialize());
 
    //check xAOD::TrackMeasurementValidation collection
    ATH_CHECK(m_pixelKey.initialize());
    ATH_MSG_INFO("Using " << m_pixelKey.key() << " as the source collection for EventInfoPixelDecorator.");
 
    // Set up the text-parsing machinery
    if (!m_selectionString.empty()) {
       ATH_CHECK( initializeParser(m_selectionString));
    }
 
    {
       std::vector<std::string> names;
       names.resize(kNIntDecor);
       names[kperModuleMultiplicity] ="perModuleMultiplicity";
       names[klayer]                 ="layer";
       names[keta_module]            ="eta_module";
       names[kphi_module]            ="phi_module";
       createDecoratorKeys(*this,m_eventInfoKey,m_decorationPrefix.value()+"_",names, m_intDecorKeys);
    }
 
    return StatusCode::SUCCESS;
  }

Member Data Documentation

m_decorationPrefix

Gaudi::Property<std::string> DerivationFramework::EventInfoPixelDecorator::m_decorationPrefix { this, "DecorationPrefix", "", ""}

private

Definition at line 43 of file EventInfoPixelDecorator.h.

m_eventInfoKey

SG::ReadHandleKey<xAOD::EventInfo> DerivationFramework::EventInfoPixelDecorator::m_eventInfoKey { this, "EventInfoKey", "EventInfo", ""}

private

Definition at line 41 of file EventInfoPixelDecorator.h.

m_intDecorKeys

std::vector<SG::WriteDecorHandleKey<xAOD::EventInfo> > DerivationFramework::EventInfoPixelDecorator::m_intDecorKeys

private

Definition at line 53 of file EventInfoPixelDecorator.h.

m_npass

std::atomic<unsigned int> DerivationFramework::EventInfoPixelDecorator::m_npass {}

mutableprivate

Definition at line 39 of file EventInfoPixelDecorator.h.

m_ntot

std::atomic<unsigned int> DerivationFramework::EventInfoPixelDecorator::m_ntot {}

mutableprivate

Definition at line 38 of file EventInfoPixelDecorator.h.

m_pixelKey

SG::ReadHandleKey<xAOD::TrackMeasurementValidationContainer> DerivationFramework::EventInfoPixelDecorator::m_pixelKey { this, "TrackMeasurementValidationKey", "PixelClusters", ""}

private

Definition at line 45 of file EventInfoPixelDecorator.h.

m_selectionString

Gaudi::Property<std::string> DerivationFramework::EventInfoPixelDecorator::m_selectionString { this, "SelectionString", "" , "" }

private

Definition at line 35 of file EventInfoPixelDecorator.h.

---

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

• EventInfoPixelDecorator.h
• EventInfoPixelDecorator.cxx