SCT_CalibHvTool.cxx

Go to the documentation of this file.

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

 
#include "SCT_CalibHvTool.h"
#include "SCT_CalibUtilities.h"
#include "SCT_CalibNumbers.h"
 
using namespace SCT_CalibAlgs;
 
namespace {
void
initQueue(std::queue<int>& q, const int numvals, const int setval) {
   for (int i{0}; i<numvals; ++i) q.push(setval);
}
}
 
SCT_CalibHvTool::SCT_CalibHvTool(const std::string& type, const std::string& name, const IInterface* parent) :
   base_class(type, name, parent)
{
}
 
StatusCode
SCT_CalibHvTool::initialize() {
   ATH_CHECK(m_DCSConditionsTool.retrieve());
   return StatusCode::SUCCESS;
}
 
StatusCode
SCT_CalibHvTool::finalize() {
   return StatusCode::SUCCESS;
}
 
bool
SCT_CalibHvTool::book() {
   bool result{true};
   ATH_MSG_DEBUG("Book HVTrips");
   //now initialize the queues (fill 100 spots with 0):
   std::queue<int> qtemp;
   initQueue(qtemp, m_maxq,0);
   // need to keep previous bins lbn 4 deep.
   initQueue(m_prevLBN, 4, 0);
   std::vector<std::pair <int,int>> dummy;
   //first set num events processed to 0
   m_phvtripProcessedEventsInt.insert(m_phvtripProcessedEventsInt.end(), n_elements, 1);
   m_phvtripRunningTotalInt.insert(m_phvtripRunningTotalInt.end(), n_elements, 0);
   m_phvtripHasItTripped.insert(m_phvtripHasItTripped.end(),n_elements, 0);
   m_phvtripQueue.insert(m_phvtripQueue.end(), n_elements, qtemp);
   m_summarytrips.insert(m_summarytrips.end(),n_elements, dummy);
   m_summarytripslb.insert(m_summarytripslb.end(), n_elements, dummy);
   m_phvtripHasItTripped_prev = m_phvtripHasItTripped;
   return result;
}
 
bool
SCT_CalibHvTool::fill(const bool fromData) {
   const EventContext& ctx = Gaudi::Hive::currentContext();
 
   if (fromData) {
      return fillFromData();
   }
   bool result{true};
   int lumi_block{0}; //fix me!
   const int wafersize{static_cast<int>(m_sct_waferHash->size())};
   int time_stamp{static_cast<int>(ctx.eventID().time_stamp())};
   int curr_time{time_stamp};
   int dtime{curr_time - m_phvtripPrevTime};
   int totalHits{0};
   if (curr_time<m_phvtripPrevTime) {
      ATH_MSG_ERROR("Events not sorted properly (time is going backwards!) ");
      return false;
   }
   // for first event
   if (m_phvtripFirstTime > curr_time) {
      m_phvtripFirstTime = curr_time;
   }
   bool newbin{false};
   // check if we have a new time bin, if we do pop off the back (oldest) part
   // of the queue and enqueue 0 for every wafer
   int maxtbins{5};
   if (dtime > maxtbins) {
      newbin = true;
      ATH_MSG_DEBUG(" new t bin " << dtime <<" since start "<< (curr_time-(m_phvtripFirstTime-1.01)));
      m_phvtripHasItTripped_prev = m_phvtripHasItTripped;
      for (int iwaf{0}; iwaf!=n_elements; ++iwaf) {
         m_phvtripQueue[iwaf].pop();
         m_phvtripQueue[iwaf].push(0);
         m_phvtripHasItTripped[iwaf]=0;
      }
   }
   // Loop over all hits in the event
   bool isgoodnow;
   for( int itrk{0}; itrk!=wafersize; ++itrk) {
      int waferhash{(*m_sct_waferHash)[itrk]};
      Identifier waferId{m_pSCTHelper->wafer_id(waferhash)};
      Identifier moduleId{m_pSCTHelper->module_id(waferId)};
      //step one is to make sure this one isn't already know to be messed up:
      isgoodnow = m_DCSConditionsTool->isGood(moduleId,InDetConditions::SCT_MODULE);
      ATH_MSG_DEBUG("checked is good "<< isgoodnow);
      if (isgoodnow) {
         int numhits{(*m_sct_numHitsInWafer)[itrk]};
         totalHits += numhits;
         double limit{m_relativetriplimit * static_cast<double>(m_phvtripRunningTotalInt[waferhash]) / (static_cast<double>(curr_time)-(static_cast<double>(m_phvtripFirstTime-0.01)))};
 
         ATH_MSG_DEBUG("waferhash "<<waferhash<<" itrk "<<itrk <<" run tot "<< m_phvtripRunningTotalInt[waferhash]
                       <<" num hits "<<numhits<<" evnts processed "<<  m_phvtripProcessedEventsInt[waferhash] <<" limit "<<limit
                       <<" abs limit, rel limit "<<m_absolutetriplimit<<","<< m_relativetriplimit<<" max bins "<<maxtbins);
 
         if ( ((numhits + m_phvtripQueue[waferhash].back() )/ static_cast<double>(maxtbins) ) > m_absolutetriplimit
               and ((numhits + m_phvtripQueue[waferhash].back() )/ static_cast<double>(maxtbins) ) > limit) {
            //read back queue and see if the high hit rate is persistant for 3 previous bins (a bin is 5 seconds by default)
            for (int iq{0}; iq<m_maxq; ++iq) {
               m_tq[iq]=m_phvtripQueue[waferhash].front();
               m_phvtripQueue[waferhash].pop();
            }
            if (m_tq[m_maxq-2]/static_cast<double>(maxtbins) > limit && m_tq[m_maxq-3]/static_cast<double>(maxtbins) > limit && m_tq[m_maxq-4]/static_cast<double>(maxtbins) > limit) {
               // found a persistent trip/thing
               ATH_MSG_DEBUG("Found a potential trip in SCT wafer " << waferhash);
               ATH_MSG_DEBUG("TRIP more info: limit = " << limit << " abslimit = " << m_absolutetriplimit
                             << " hits/secs now = " << ((numhits + m_tq[m_maxq-1])/ maxtbins )
                             << " running total = " << m_phvtripRunningTotalInt[waferhash]
                             << " running average =  " <<  m_phvtripRunningTotalInt[waferhash] / (curr_time-(m_phvtripFirstTime-1.01)));
               m_phvtripHasItTripped[waferhash] = (numhits + m_tq[m_maxq-1] );
            } // end trip found
 
            //put queue back the way we found it:
            for (int iq{0}; iq < m_maxq; ++iq) {
               m_phvtripQueue[waferhash].push(m_tq[iq]);
            }
         } // end potential trip id
         //incriment event counter
         m_phvtripProcessedEventsInt[waferhash] += 1;
         //recalcuate running total:
         m_phvtripRunningTotalInt[waferhash] += numhits;
         m_phvtripQueue[waferhash].back() += numhits;
      } //end is good check
      else {
         ATH_MSG_DEBUG("Was a bad module already: " << waferhash);
      }
   } //end loop over tracks
 
   // text output, summarizing the event
   //if new bin summarzie the previous event.
   if (newbin) {
      SCT_ID::const_id_iterator waferItr{m_waferItrBegin};
      SCT_ID::const_id_iterator waferItrE{m_waferItrEnd};
      for (; waferItr not_eq waferItrE; ++waferItr) {
         Identifier waferId{*waferItr};
         IdentifierHash waferHash{m_pSCTHelper->wafer_hash(waferId)};
         if (m_phvtripHasItTripped_prev[waferHash]>0) {
            std::pair<int, int> wp;
            wp.first = (m_phvtripPrevTime - 3*maxtbins);
            wp.second = (m_phvtripPrevTime + maxtbins);
            m_summarytrips[waferHash].push_back(wp);
            std::pair<int,int> lbn;
            lbn.first = m_prevLBN.front();
            lbn.second = m_prevLBN.back();
            m_summarytripslb[waferHash].push_back(lbn);
         }
      }
      m_phvtripPrevTime = curr_time; //reset prev and current time
      m_prevLBN.pop();
      m_prevLBN.push(lumi_block);
   }
 
   for (int iwaf{0}; iwaf!=n_elements; ++iwaf) {
      ATH_MSG_DEBUG("mod "<< iwaf <<" events with that waf " <<m_phvtripProcessedEventsInt[iwaf] << " numhits "<<m_phvtripRunningTotalInt[iwaf]);
      ATH_MSG_DEBUG("did it trip "<<m_phvtripHasItTripped[iwaf]);
      if (m_phvtripHasItTripped[iwaf]) {
         ATH_MSG_DEBUG("TRIP more info: abslimit = "<<  m_absolutetriplimit
                       <<" hits/secs now = "<< ( m_phvtripQueue[iwaf].back() / maxtbins )
                       <<" running total = "<< m_phvtripRunningTotalInt[iwaf]
                       <<" running average =  "<<  m_phvtripRunningTotalInt[iwaf] / (curr_time-(m_phvtripFirstTime-1.01)));
      }
   }
   if (m_outputLowHits && (totalHits < m_lowHitCut) ) {
   }
   return result;
}
 
bool
SCT_CalibHvTool::fillFromData() {
   bool result{true};
   return result;
}