MuonContainerMergingAlg.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/

//   MuonContainerMergingAlg   
//   author Maria Mironova, Simone Pagan Griso, Kehang Bai

// Generic algorithm for merging a list of muon containers
#include "MuonAnalysisAlgorithms/MuonContainerMergingAlg.h"
#include "xAODMuon/MuonContainer.h"
#include <xAODMuon/MuonAuxContainer.h>
#include <AthContainers/ConstDataVector.h>
#include <AsgDataHandles/WriteHandle.h>
 
namespace CP {
    MuonContainerMergingAlg::MuonContainerMergingAlg( const std::string& name, ISvcLocator* svcLoc )
    : EL::AnaReentrantAlgorithm( name, svcLoc ){
        //nothing to do here
    }
 
    StatusCode MuonContainerMergingAlg::initialize() {
 
        // initialize list of muon containters
        ATH_MSG_DEBUG("Initializing MuonContainerMergingAlg");
        ATH_CHECK( m_inputMuonContainers.initialize() );
       
        // Override m_outMuonLocationCopy key 
        if ( !m_createViewCollection.value() ) { m_outMuonLocationCopy = m_outMuonLocationView.key(); }
        
        // initialize output muon containters   
    ATH_CHECK( m_outMuonLocationCopy.initialize(!m_createViewCollection.value()) );
        ATH_CHECK( m_outMuonLocationView.initialize(m_createViewCollection.value()) );
    
        // Return gracefully:
        return StatusCode::SUCCESS;
    }
 
    StatusCode MuonContainerMergingAlg::execute(const EventContext &ctx) const  {
 
        // Setup containers for output, to avoid const conversions setup two different kind of containers
        auto outputViewCol = std::make_unique<ConstDataVector<xAOD::MuonContainer>>(SG::VIEW_ELEMENTS);
        auto outputCol = std::make_unique<xAOD::MuonContainer>();
 
        std::unique_ptr<xAOD::MuonAuxContainer> outputAuxCol;
        if(!m_createViewCollection) {
            outputAuxCol = std::make_unique<xAOD::MuonAuxContainer>();
            outputCol->setStore(outputAuxCol.get());
        }
 
        // retrieve muons from StoreGate
        std::vector<const xAOD::MuonContainer*> muonCollections;
        muonCollections.reserve(m_inputMuonContainers.size());
        size_t muonCountSum = 0;
        for (const auto& mcname : m_inputMuonContainers){
          // Retrieve tracks from StoreGate
          SG::ReadHandle<xAOD::MuonContainer> muCol (mcname, ctx);
          muonCollections.push_back(muCol.cptr());
          muonCountSum += muCol->size();
        }
        if (m_createViewCollection) {
          outputViewCol->reserve(muonCountSum);  
        } else {
          outputCol->reserve(muonCountSum);
        }

        for (auto& mciter : muonCollections) {
          if (mciter and !mciter->empty()) {
            ATH_MSG_DEBUG("Size of muon collection " << mciter->size());
            for (const auto* const muon : *mciter) {
              if (m_createViewCollection) {
                outputViewCol->push_back(muon);
              } else {
                xAOD::Muon* newMuon = new xAOD::Muon(*muon);
                outputCol->push_back(newMuon);
              }
            }
          }
        }
 
        // write output to SG
        if (m_createViewCollection) {
          SG::WriteHandle< ConstDataVector<xAOD::MuonContainer> > h_write(m_outMuonLocationView, ctx);
          ATH_CHECK(h_write.record(std::move(outputViewCol)));
        }
        else {
          SG::WriteHandle<xAOD::MuonContainer> h_write(m_outMuonLocationCopy, ctx);
          ATH_CHECK(h_write.record(std::move(outputCol), std::move(outputAuxCol)));
        }
 
        return StatusCode::SUCCESS;
 
    }
 
}