AFPSiLayerAlgorithm.cxx

Go to the documentation of this file.

/* 
  Copyright (C) 2002-2026 CERN for the benefit of the ATLAS collaboration
*
*
*   AFPSiLayerAlgorithm
*
*
*/
 
#include "Run3AFPMonitoring/AFPSiLayerAlgorithm.h"
#include "StoreGate/ReadHandleKey.h"
#include "xAODForward/AFPStationID.h"
 
 
namespace {
    constexpr int reorganizePlanes(const int station, const int layer)
    {
        bool reverse = station == 0 || station == 1;
        return station * 4 + (reverse ? 3 - layer : layer);
    }
}
 
 
AFPSiLayerAlgorithm::AFPSiLayerAlgorithm( const std::string& name, ISvcLocator* pSvcLocator )
:AthMonitorAlgorithm(name,pSvcLocator)
, m_afpHitContainerKey("AFPSiHitContainer")
{
    declareProperty("AFPSiHitContainer", m_afpHitContainerKey);
}
 
 
AFPSiLayerAlgorithm::~AFPSiLayerAlgorithm() {}
 
 
StatusCode AFPSiLayerAlgorithm::initialize() {
 
    using namespace Monitored;
 
    m_StationPlaneGroup = buildToolMap<std::map<std::string,int>>(m_tools,"AFPSiLayerTool", m_stationnames, m_pixlayers);
    m_StationGroup = buildToolMap<int>(m_tools, "AFPSiLayerTool", m_stationnames);
 
 
    // We must declare to the framework in initialize what SG objects we are going to use:
    SG::ReadHandleKey<xAOD::AFPSiHitContainer> afpHitContainerKey("AFPSiHits");
    ATH_CHECK(m_afpHitContainerKey.initialize());
    
    ATH_MSG_INFO( "BunchCrossingKey initialization (SiT)" );
    ATH_CHECK(m_bunchCrossingKey.initialize());
    ATH_MSG_INFO( "initialization completed (SiT)" );
    return AthMonitorAlgorithm::initialize();
}
 
StatusCode AFPSiLayerAlgorithm::fillHistograms( const EventContext& ctx ) const {
    using namespace Monitored;
    
    auto bcidAll    = Monitored::Scalar<int>("bcidAll", 0);
    auto bcidFront  = Monitored::Scalar<int>("bcidFront", 0);
    auto bcidMiddle = Monitored::Scalar<int>("bcidMiddle", 0);
    auto bcidEnd    = Monitored::Scalar<int>("bcidEnd", 0);
    
    auto numberOfEventsPerLumiblockFront  = Monitored::Scalar<int>("numberOfEventsPerLumiblockFront", 0);
    auto numberOfEventsPerLumiblockMiddle = Monitored::Scalar<int>("numberOfEventsPerLumiblockMiddle", 0);
    auto numberOfEventsPerLumiblockEnd    = Monitored::Scalar<int>("numberOfEventsPerLumiblockEnd", 0);
 
    SG::ReadHandle<xAOD::EventInfo> eventInfo = GetEventInfo(ctx);
    numberOfEventsPerLumiblockFront   = eventInfo->lumiBlock();
    numberOfEventsPerLumiblockMiddle  = eventInfo->lumiBlock();
    numberOfEventsPerLumiblockEnd     = eventInfo->lumiBlock();
 
    // BCX handler
    const unsigned int bcid = eventInfo->bcid();
    const unsigned int run_number = eventInfo->runNumber();
    SG::ReadCondHandle<BunchCrossingCondData> bcidHdl(m_bunchCrossingKey,ctx);
    if (!bcidHdl.isValid()) {
        ATH_MSG_ERROR( "Unable to retrieve BunchCrossing conditions object (SiT)" );
    }
    const BunchCrossingCondData* bcData{*bcidHdl};
 
    int bx_diff = int(run_number) >= 521800 ? 2 : 1;
    // Classifying bunches by position in train (Front, Middle, End)
    enum { FRONT, MIDDLE, END, NPOS } position = NPOS;
    if(bcData->isFilled(bcid))
    {
        bcidAll = bcid;
        fill("AFPSiLayerTool", bcidAll);
        if(!bcData->isFilled(bcid-bx_diff))
        {
            position = FRONT;
            bcidFront = bcid;
            fill("AFPSiLayerTool", bcidFront);
            fill("AFPSiLayerTool", numberOfEventsPerLumiblockFront);
        }
        else
        {
            if(bcData->isFilled(bcid+bx_diff))
            {
                position = MIDDLE;
                bcidMiddle = bcid;
                fill("AFPSiLayerTool", bcidMiddle);
                fill("AFPSiLayerTool", numberOfEventsPerLumiblockMiddle);
            }
            else
            {
                position = END;
                bcidEnd = bcid;
                fill("AFPSiLayerTool", bcidEnd);
                fill("AFPSiLayerTool", numberOfEventsPerLumiblockEnd);
            }
        }
    }
    
        
    // Declare the quantities which should be monitored:
    auto lb       = Monitored::Scalar<int>("lb", 0);
    auto muPerBX  = Monitored::Scalar<float>("muPerBX", 0.0);
    //auto run = Monitored::Scalar<int>("run",0);
 
    auto nSiHits = Monitored::Scalar<int>("nSiHits", 1);
 
    auto pixelRowIDChip = Monitored::Scalar<int>("pixelRowIDChip", 0); 
    auto pixelColIDChip = Monitored::Scalar<int>("pixelColIDChip", 0); 
 
    auto timeOverThreshold = Monitored::Scalar<int>("timeOverThreshold", 0);
 
    auto clusterX = Monitored::Scalar<float>("clusterX", 0.0);
    auto clusterY = Monitored::Scalar<float>("clusterY", 0.0); 
    auto clustersInPlanes = Monitored::Scalar<int>("clustersInPlanes", 0);
 
    auto trackX = Monitored::Scalar<float>("trackX", 0.0);
    auto trackY = Monitored::Scalar<float>("trackY", 0.0);
    
    auto planeHits        = Monitored::Scalar<int>("planeHits", 0);
    
    auto numberOfHitsPerStation = Monitored::Scalar<int>("numberOfHitsPerStation", 0);
    
    auto lbEvents             = Monitored::Scalar<int>("lbEvents", 0);
    auto lbHits               = Monitored::Scalar<int>("lbHits", 0);
    auto lbEventsStations     = Monitored::Scalar<int>("lbEventsStations", 0);
    auto lbEventsStationsAll  = Monitored::Scalar<int>("lbEventsStationsAll", 0);
    
    auto planes = Monitored::Scalar<int>("planes", 0);
    
    auto eventsPerStation = Monitored::Scalar<int>("eventsPerStation", 0);
    
    auto clusterToT = Monitored::Scalar<int>("clusterToT", 0);
    
    lb        = eventInfo->lumiBlock();
    lbEvents  = eventInfo->lumiBlock();
    //muPerBX   = lbAverageInteractionsPerCrossing(ctx);
    muPerBX   = lbInteractionsPerCrossing(ctx);
    if (muPerBX == 0.0) {
        ATH_MSG_DEBUG("AverageInteractionsPerCrossing is 0, forcing to 1.0");
        muPerBX=1.0;
    }
    fill("AFPSiLayerTool", lb, muPerBX);
    fill("AFPSiLayerTool", lbEvents);
    
    
    SG::ReadHandle<xAOD::AFPSiHitContainer> afpHitContainer(m_afpHitContainerKey, ctx);
    if(! afpHitContainer.isValid())
    {
        ATH_MSG_WARNING("evtStore() does not contain hits collection with name " << m_afpHitContainerKey);
        return StatusCode::SUCCESS;
    }
 
    ATH_CHECK( afpHitContainer.initialize() );
 
    nSiHits = afpHitContainer->size();
    fill("AFPSiLayerTool", lb, nSiHits);
 
    int eventsInStations[4] = {};
    int numberOfHitsPerPlane[4][4] = {};
 
    for(const xAOD::AFPSiHit *hitsItr: *afpHitContainer)
    {
        lb                  = eventInfo->lumiBlock();
        lbHits              = eventInfo->lumiBlock();
        lbEventsStations    = eventInfo->lumiBlock();
        lbEventsStationsAll = eventInfo->lumiBlock();
        pixelRowIDChip      = hitsItr->pixelRowIDChip();
        pixelColIDChip      = hitsItr->pixelColIDChip();
        timeOverThreshold   = hitsItr->timeOverThreshold();
        
        
        if (hitsItr->stationID()<4 && hitsItr->stationID()>=0 && hitsItr->pixelLayerID()<4 && hitsItr->pixelLayerID()>=0) 
        {
            ++eventsInStations[hitsItr->stationID()];
 
            fill(m_tools[m_StationPlaneGroup.at(m_stationnames.at(hitsItr->stationID())).at(m_pixlayers.at(hitsItr->pixelLayerID()))], pixelRowIDChip, pixelColIDChip);
            fill(m_tools[m_StationPlaneGroup.at(m_stationnames.at(hitsItr->stationID())).at(m_pixlayers.at(hitsItr->pixelLayerID()))], pixelRowIDChip);
            fill(m_tools[m_StationPlaneGroup.at(m_stationnames.at(hitsItr->stationID())).at(m_pixlayers.at(hitsItr->pixelLayerID()))], pixelColIDChip);
            fill(m_tools[m_StationPlaneGroup.at(m_stationnames.at(hitsItr->stationID())).at(m_pixlayers.at(hitsItr->pixelLayerID()))], timeOverThreshold);
 
            planeHits = hitsItr->pixelLayerID();
            fill(m_tools[m_StationGroup.at(m_stationnames.at(hitsItr->stationID()))], planeHits);
            
            ++numberOfHitsPerPlane[hitsItr->stationID()][hitsItr->pixelLayerID()];
            numberOfHitsPerStation = hitsItr->stationID();
            fill("AFPSiLayerTool", numberOfHitsPerStation);
 
            fill("AFPSiLayerTool", lbHits);
        }
        else ATH_MSG_WARNING("Unrecognised station index: " << hitsItr->stationID());
    }
 
    auto hitsPerPlaneProfile          = Monitored::Scalar<float>("hitsPerPlaneProfile", 0.0);
    auto lbhitsPerPlaneProfile       = Monitored::Scalar<int>("lbhitsPerPlaneProfile", 0);
    auto hitsPerPlaneEventsMu        = Monitored::Scalar<float>("hitsPerPlaneEventsMu", 0.0);
    auto hitPerPlaneEventMuIndex     = Monitored::Scalar<int>("hitPerPlaneEventMuIndex", 0);
 
    lbhitsPerPlaneProfile = eventInfo->lumiBlock();
    for(int i_station = 0; i_station < 4; i_station++)
        for(int j_layer = 0; j_layer < 4; j_layer++)
        {
            hitsPerPlaneProfile = numberOfHitsPerPlane[i_station][j_layer]/muPerBX;
            fill(m_tools[m_StationPlaneGroup.at(m_stationnames.at(i_station)).at(m_pixlayers.at(j_layer))], lbhitsPerPlaneProfile, hitsPerPlaneProfile);
            if (muPerBX != 0.0) {
                hitsPerPlaneEventsMu = numberOfHitsPerPlane[i_station][j_layer] / muPerBX;
            }
            hitPerPlaneEventMuIndex = reorganizePlanes(i_station, j_layer);
            fill("AFPSiLayerTool", hitPerPlaneEventMuIndex, hitsPerPlaneEventsMu);
        }
        
    bool noEventsInStations = true;
    for(int i=0; i<4; i++)
    {
        if(eventsInStations[i]>0) {
            fill(m_tools[m_StationGroup.at(m_stationnames.at(i))], lbEventsStations);
            
            eventsPerStation = i * 4;
            fill("AFPSiLayerTool", eventsPerStation);
            ++eventsPerStation;
            fill("AFPSiLayerTool", eventsPerStation);
            ++eventsPerStation;
            fill("AFPSiLayerTool", eventsPerStation);
            ++eventsPerStation;
            fill("AFPSiLayerTool", eventsPerStation);
 
            noEventsInStations = false;
        }
    }
    if(!noEventsInStations)
    {
        fill("AFPSiLayerTool", lbEventsStationsAll);
    }
    
    // Filling of cluster and track 2D histograms
    AFPMon::AFPFastReco fast(afpHitContainer.get());
    fast.reco();
 
    // Track histograms:
    unsigned int totalTracksAll[4] = {};
    unsigned int totalTracksFront[4] = {};
    unsigned int totalTracksMiddle[4] = {};
    unsigned int totalTracksEnd[4] = {};
        
    for (const auto& track : fast.tracks()) 
    {
        trackX = track.x * 1.0;
        trackY = track.y * 1.0;
        fill(m_tools[m_StationGroup.at(m_stationnames.at(track.station))], trackY, trackX);
        
        if (position == FRONT)
        {
            ++totalTracksFront[track.station];
            ++totalTracksAll[track.station];
        }
        else if (position == MIDDLE)
        {
            ++totalTracksMiddle[track.station];
            ++totalTracksAll[track.station];
        }
        else if (position == END)
        {
            ++totalTracksEnd[track.station];
            ++totalTracksAll[track.station];
        }
    }
    
    auto lbTracksAll        = Monitored::Scalar<int>("lbTracksAll", 0);
    auto lbTracksFront      = Monitored::Scalar<int>("lbTracksFront", 0);
    auto lbTracksMiddle     = Monitored::Scalar<int>("lbTracksMiddle", 0);
    auto lbTracksEnd        = Monitored::Scalar<int>("lbTracksEnd", 0);
 
    auto Total_tracks_All_profile     = Monitored::Scalar<float>("Total_tracks_All_profile", 0.0);
    auto Total_tracks_Front_profile   = Monitored::Scalar<float>("Total_tracks_Front_profile", 0.0);
    auto Total_tracks_Middle_profile  = Monitored::Scalar<float>("Total_tracks_Middle_profile", 0.0);
    auto Total_tracks_End_profile     = Monitored::Scalar<float>("Total_tracks_End_profile", 0.0);
    
    lbTracksAll     = eventInfo->lumiBlock();
    lbTracksFront   = eventInfo->lumiBlock();
    lbTracksMiddle  = eventInfo->lumiBlock();
    lbTracksEnd     = eventInfo->lumiBlock();
    
    for(int i = 0; i < 4; i++)
    {
        Total_tracks_All_profile = totalTracksAll[i] / muPerBX;
        fill(m_tools[m_StationGroup.at(m_stationnames.at(i))], lbTracksAll, Total_tracks_All_profile);
        totalTracksAll[i] = 0;
 
        Total_tracks_Front_profile = totalTracksFront[i] / muPerBX;
        if (position == FRONT)
            fill(m_tools[m_StationGroup.at(m_stationnames.at(i))], lbTracksFront, Total_tracks_Front_profile);
        totalTracksFront[i] = 0;
        
        Total_tracks_Middle_profile = totalTracksMiddle[i] / muPerBX;
        if (position == MIDDLE)
            fill(m_tools[m_StationGroup.at(m_stationnames.at(i))], lbTracksMiddle, Total_tracks_Middle_profile);
        totalTracksMiddle[i] = 0;
        
        Total_tracks_End_profile = totalTracksEnd[i] / muPerBX;
        if (position == END)
            fill(m_tools[m_StationGroup.at(m_stationnames.at(i))], lbTracksEnd, Total_tracks_End_profile);
        totalTracksEnd[i] = 0;
    }
    
    // Cluster histograms 
    unsigned int totalClustersAll[4][4] = {};
    unsigned int totalClustersFront[4][4] = {};
    unsigned int totalClustersMiddle[4][4] = {};
    unsigned int totalClustersEnd[4][4] = {};
 
    auto clustersPerPlaneAllPP        = Monitored::Scalar<float>("clustersPerPlaneAllPP", 0.0);
    auto clustersPerPlaneFrontPP        = Monitored::Scalar<float>("clustersPerPlaneFrontPP", 0.0);
    auto clustersPerPlaneMiddlePP       = Monitored::Scalar<float>("clustersPerPlaneMiddlePP", 0.0);
    auto clustersPerPlaneEndPP          = Monitored::Scalar<float>("clustersPerPlaneEndPP", 0.0);
 
    auto lbClustersPerPlanesAll        = Monitored::Scalar<int>("lbClustersPerPlanesAll", 0);
    auto lbClustersPerPlanesFront      = Monitored::Scalar<int>("lbClustersPerPlanesFront", 0);
    auto lbClustersPerPlanesMiddle     = Monitored::Scalar<int>("lbClustersPerPlanesMiddle", 0);
    auto lbClustersPerPlanesEnd        = Monitored::Scalar<int>("lbClustersPerPlanesEnd", 0);
    
    lbClustersPerPlanesAll     = eventInfo->lumiBlock();
    lbClustersPerPlanesFront   = eventInfo->lumiBlock();
    lbClustersPerPlanesMiddle  = eventInfo->lumiBlock();
    lbClustersPerPlanesEnd     = eventInfo->lumiBlock();
 
    for(const auto& cluster : fast.clusters()) 
    {
        clusterX = cluster.x * 1.0;
        clusterY = cluster.y * 1.0;
        fill(m_tools[m_StationPlaneGroup.at(m_stationnames.at(cluster.station)).at(m_pixlayers.at(cluster.layer))], clusterY, clusterX);
        if (cluster.station == 0 || cluster.station == 1)
        {
            clustersInPlanes = reorganizePlanes(cluster.station, cluster.layer);
        }
        else
        {
            clustersInPlanes = (cluster.station*4)+cluster.layer;
        }
        fill("AFPSiLayerTool", clustersInPlanes);
        
        clusterToT = cluster.sumToT;
        fill(m_tools[m_StationPlaneGroup.at(m_stationnames.at(cluster.station)).at(m_pixlayers.at(cluster.layer))], clusterToT);
 
        if (position == FRONT)
        {
            ++totalClustersFront[cluster.station][cluster.layer];
            ++totalClustersAll[cluster.station][cluster.layer];
        }
        else if (position == MIDDLE)
        {
            ++totalClustersMiddle[cluster.station][cluster.layer];
            ++totalClustersAll[cluster.station][cluster.layer];
        }
        else if (position == END)
        {
            ++totalClustersEnd[cluster.station][cluster.layer];
            ++totalClustersAll[cluster.station][cluster.layer];
        }
    }
 
    for(int i_station = 0; i_station < 4; i_station++)
        for(int j_layer = 0; j_layer < 4; j_layer++)
        {
            clustersPerPlaneAllPP = totalClustersAll[i_station][j_layer] / muPerBX;
            fill(m_tools[m_StationPlaneGroup.at(m_stationnames.at(i_station)).at(m_pixlayers.at(j_layer))], lbClustersPerPlanesAll, clustersPerPlaneAllPP);
            totalClustersAll[i_station][j_layer] = 0;
 
            clustersPerPlaneFrontPP = totalClustersFront[i_station][j_layer] / muPerBX;
            if (position == FRONT)
                fill(m_tools[m_StationPlaneGroup.at(m_stationnames.at(i_station)).at(m_pixlayers.at(j_layer))], lbClustersPerPlanesFront, clustersPerPlaneFrontPP);
            totalClustersFront[i_station][j_layer] = 0;
 
            clustersPerPlaneMiddlePP = totalClustersMiddle[i_station][j_layer] / muPerBX;
            if (position == MIDDLE)
                fill(m_tools[m_StationPlaneGroup.at(m_stationnames.at(i_station)).at(m_pixlayers.at(j_layer))], lbClustersPerPlanesMiddle, clustersPerPlaneMiddlePP);
            totalClustersMiddle[i_station][j_layer] = 0;
 
            clustersPerPlaneEndPP = totalClustersEnd[i_station][j_layer] / muPerBX;
            if (position == END)
                fill(m_tools[m_StationPlaneGroup.at(m_stationnames.at(i_station)).at(m_pixlayers.at(j_layer))], lbClustersPerPlanesEnd, clustersPerPlaneEndPP);
            totalClustersEnd[i_station][j_layer] = 0;
        }
 
    return(fillHistogramsPlaneEff(*afpHitContainer));
} // end of fillHistograms
 
StatusCode AFPSiLayerAlgorithm::fillHistogramsPlaneEff(const xAOD::AFPSiHitContainer& afpHitContainer) const {
    using namespace Monitored;
 
    // Define 2D histograms for tries and successes
    Monitored::Scalar<bool> sit_plane_eff_passed("sit_plane_eff_passed", false);
    Monitored::Scalar<float> sit_plane_eff_triedX("sit_plane_eff_triedX", 0.0);
    Monitored::Scalar<float> sit_plane_eff_triedY("sit_plane_eff_triedY", 0.0);
 
    auto triesX = Monitored::Scalar<float>("triesX", 0.0);
    auto triesY = Monitored::Scalar<float>("triesY", 0.0);
 
    auto successX = Monitored::Scalar<float>("successX", 0.0);
    auto successY = Monitored::Scalar<float>("successY", 0.0);
 
 
    auto clusterXLocal = Monitored::Scalar<float>("clusterXLocal", 0.0);
    auto clusterYLocal = Monitored::Scalar<float>("clusterYLocal", 0.0);
 
    auto clusterXTag = Monitored::Scalar<float>("clusterXTag", 0.0);
    auto clusterYTag = Monitored::Scalar<float>("clusterYTag", 0.0); 
 
    AFPMon::AFPFastReco fast(&afpHitContainer);
    fast.reco();
 
    int min_hits[4] = {2, 3, 3, 3};
    int numStations = 4;
    int numPlanes = 4;
    
    // Count clusters per station and plane
    std::vector<std::vector<int>> nclusters_planes(numStations, std::vector<int>(numPlanes, 0));
    for (const auto& cluster : fast.clusters()) {
        ++nclusters_planes[cluster.station][cluster.layer];
    }
    
    // Count planes with clusters at each station
    std::vector<std::vector<int>> nclusters_other_planes(numStations, std::vector<int>(numPlanes, 0));
    for (int iStation = 0; iStation < numStations; ++iStation) {
        for (int iPlane = 0; iPlane < numPlanes; ++iPlane) {
            for (int ip = 0; ip < numPlanes; ++ip) {
                if (iPlane != ip && nclusters_planes[iStation][ip] > 0) {
                    ++nclusters_other_planes[iStation][iPlane];
                }
            }
        }
    }
    
    // Precomputed tag planes
    std::array<std::set<int>, 4> precomputed_tag_planes = {
        std::set<int>{1, 2, 3},      
        std::set<int>{0, 1, 2, 3},   
        std::set<int>{0, 1, 2, 3},   
        std::set<int>{0, 1, 2, 3}       
    };
    
    // Precomputed clusters by Stations and Planes  
    using ClusterType = std::decay_t<decltype(*fast.clusters().begin())>;
    std::vector<std::vector<std::vector<ClusterType>>> clusters_by_station_layer(numStations, 
        std::vector<std::vector<ClusterType>>(numPlanes));
    
    for (const auto& cluster : fast.clusters()) {
        clusters_by_station_layer[cluster.station][cluster.layer].push_back(cluster);
    }
    
    for (int iStation = 0; iStation < numStations; ++iStation) {
        for (int iPlane = 0; iPlane < numPlanes; ++iPlane) {
            // Skip if no enough hits in other planes
            if (nclusters_other_planes[iStation][iPlane] < min_hits[iStation]) {continue;}
            
            std::set<int> tag_planes = precomputed_tag_planes[iStation];
            tag_planes.erase(iPlane);
 
            std::vector<int> v_tag_planes(tag_planes.begin(), tag_planes.end());
            int seed = v_tag_planes[0]; // Take first plane as a seed
            
            // Iterate only stations coresponding to the given Station and Layer
            for (const auto& cluster : clusters_by_station_layer[iStation][seed]) {
                int found_in_other_tag_planes = 0;
                
                // Check other tag planes using lookup
                for (size_t i = 1; i < v_tag_planes.size(); ++i) {
                    int tag_plane = v_tag_planes[i];
                    bool found = false;
                    
                    for (const auto& clusterTag : clusters_by_station_layer[iStation][tag_plane]) {
                        if (std::fabs(cluster.x - clusterTag.x) < 2.0 && std::fabs(cluster.y - clusterTag.y) < 2.0) {
                            found = true;
                            ++found_in_other_tag_planes;
                            break;
                        }
                    }
                    // if not found in current tag plane
                    if (!found) {break;}
                }
                
                // Check if valid tag
                if (found_in_other_tag_planes == static_cast<int>(v_tag_planes.size()) - 1) {
                    sit_plane_eff_triedX = cluster.x;
                    sit_plane_eff_triedY = cluster.y;
                    
                    // Check probe plane 
                    for (const auto& clusterTag : clusters_by_station_layer[iStation][iPlane]) {
                        if (std::fabs(cluster.x - clusterTag.x) < 2.0 && std::fabs(cluster.y - clusterTag.y) < 2.0) {
                            sit_plane_eff_passed = true;
                            break;
                        }
                    }
                    
                    fill(m_tools[m_StationPlaneGroup.at(m_stationnames.at(iStation)).at(m_pixlayers.at(iPlane))], 
                         sit_plane_eff_passed, sit_plane_eff_triedY, sit_plane_eff_triedX);
                    sit_plane_eff_passed = false;
                }
            }
        }
    }
    return StatusCode::SUCCESS;
}