StripClusterAnalysisAlg.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#include "StripClusterAnalysisAlg.h"
#include "AthenaMonitoringKernel/Monitored.h"
 
namespace ActsTrk {
 
  StripClusterAnalysisAlg::StripClusterAnalysisAlg(const std::string& name, ISvcLocator *pSvcLocator)
    : AthMonitorAlgorithm(name, pSvcLocator) 
  {}
  
  StatusCode StripClusterAnalysisAlg::initialize() {
    ATH_MSG_DEBUG( "Initializing " << name() << " ... " );
    
    ATH_CHECK( m_stripClusterContainerKey.initialize() );
    ATH_CHECK(detStore()->retrieve(m_stripID,"SCT_ID"));
 
    ATH_MSG_DEBUG("Monitoring settings ...");
    ATH_MSG_DEBUG(m_monGroupName);
 
    return AthMonitorAlgorithm::initialize();
  }
 
  StatusCode StripClusterAnalysisAlg::fillHistograms(const EventContext& ctx) const {
    ATH_MSG_DEBUG(" In " << name() << "::fillHistograms()" );
    
    SG::ReadHandle< xAOD::StripClusterContainer > inputStripClusterContainer( m_stripClusterContainerKey, ctx );
    if (!inputStripClusterContainer.isValid()){
        ATH_MSG_FATAL("xAOD::StripClusterContainer with key " << m_stripClusterContainerKey.key() << " is not available...");
        return StatusCode::FAILURE;
    }
    
    auto monitor_barrelEndcap = Monitored::Collection("barrelEndcap", *inputStripClusterContainer,
                                                      [this] (const auto* cluster) -> int
                                                      {
                                                        const Identifier& id = m_stripID->wafer_id(cluster->identifierHash());
                                                        return m_stripID->barrel_ec(id);
                                                      });
    auto monitor_layerDisk = Monitored::Collection("layerDisk", *inputStripClusterContainer,
                           [this] (const auto* cluster) -> int
                                                   {
                                                     const Identifier& id = m_stripID->wafer_id(cluster->identifierHash());
                                                     return m_stripID->layer_disk(id);
                                                   });
    auto monitor_phiModule = Monitored::Collection("phiModule", *inputStripClusterContainer,
                           [this] (const auto* cluster) -> int
                           {
                             const Identifier& id = m_stripID->wafer_id(cluster->identifierHash());
                             return m_stripID->phi_module(id);
                           });
    auto monitor_etaModule = Monitored::Collection("etaModule", *inputStripClusterContainer,
                           [this] (const auto* cluster) -> int
                           {
                             const Identifier& id = m_stripID->wafer_id(cluster->identifierHash());
                             return m_stripID->eta_module(id);
                           });
    auto monitor_sideModule = Monitored::Collection("sideModule", *inputStripClusterContainer,
                            [this] (const auto* cluster) -> int
                            {
                              const Identifier& id = m_stripID->wafer_id(cluster->identifierHash());
                              return m_stripID->side(id);
                            });
 
 
    auto monitor_eta = Monitored::Collection("eta", *inputStripClusterContainer,
                                             [] (const auto* cluster) -> double
                                             {
                                               const auto& globalPos = cluster->globalPosition();
                           Amg::Vector3D globalPosition(globalPos(0, 0), globalPos(1, 0), globalPos(2, 0));
                                               return globalPosition.eta();
                                             });
    auto monitor_perp = Monitored::Collection("perp", *inputStripClusterContainer,
                          [] (const auto* cluster) -> float
                          {
                        const auto& globalPos = cluster->globalPosition();
                        Amg::Vector3D globalPosition(globalPos(0, 0), globalPos(1, 0), globalPos(2, 0));
                        return globalPosition.perp();
                          });
 
    auto monitor_globalX = Monitored::Collection("globalX", *inputStripClusterContainer,
                         [] (const auto* cluster) -> float
                         { return cluster->globalPosition()(0, 0); });
    auto monitor_globalY = Monitored::Collection("globalY", *inputStripClusterContainer,
                         [] (const auto* cluster) -> float
                         { return cluster->globalPosition()(1, 0); });
    auto monitor_globalZ = Monitored::Collection("globalZ", *inputStripClusterContainer,
                         [] (const auto* cluster) -> float
                         { return cluster->globalPosition()(2, 0); });
    
    auto monitor_localX = Monitored::Collection("localX", *inputStripClusterContainer,
                        [] (const auto cluster) -> float
                        {
                          const auto& localPos = cluster->template localPosition<1>();
                          return localPos(0,0);
                        });
    
    auto monitor_localCovXX = Monitored::Collection("localCovXX", *inputStripClusterContainer,
                            [] (const auto* cluster) -> float
                            { return cluster->template localCovariance<1>()(0, 0); });
 
    auto monitor_sizeX = Monitored::Collection("sizeX", *inputStripClusterContainer,
                           [] (const auto* cluster) -> int
                           { return cluster->channelsInPhi(); });
    
    fill(m_monGroupName.value(),
     monitor_barrelEndcap, monitor_layerDisk,
     monitor_phiModule, monitor_etaModule, monitor_sideModule,
     monitor_eta, monitor_perp,
     monitor_globalX, monitor_globalY, monitor_globalZ,
     monitor_localX, monitor_localCovXX,
     monitor_sizeX);
    
    return StatusCode::SUCCESS;
  }
  
}