LArTBH6TriggerTimeTool.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#include "LArTBH6TriggerTimeTool.h"
#include "GaudiKernel/IIncidentSvc.h"
 
#include "LArSimEvent/LArHitContainer.h"
 
LArTBH6TriggerTimeTool::LArTBH6TriggerTimeTool(const std::string& type,
                                               const std::string& name,
                                               const IInterface* parent) :
  base_class(type,name,parent), m_time(0), m_newEvent(true),
  m_fixed(true)
{
  declareProperty("isFixed",m_fixed);
  declareProperty("FixedTime",m_time);
 
  m_hitcoll.push_back( SG::ReadHandle< LArHitContainer >("LArHitEMEC") );
  m_hitcoll.push_back( SG::ReadHandle< LArHitContainer >("LArHitHEC") );
  m_hitcoll.push_back( SG::ReadHandle< LArHitContainer >("LArHitFCAL") );
}
 
 
StatusCode LArTBH6TriggerTimeTool::initialize()
{
 
  SmartIF<IIncidentSvc> incsvc{Gaudi::svcLocator()->service("IncidentSvc")};
  if (!incsvc) {
    ATH_MSG_FATAL ( "Incident service not found" );
    return StatusCode::FAILURE ;
  }
 
  long int pri=100;
  incsvc->addListener(this,"BeginEvent",pri);
 
  ATH_MSG_DEBUG ( " LArTBH6TriggerTimeTool initialized " );
 
  return StatusCode::SUCCESS;
 
}
 
 
void LArTBH6TriggerTimeTool::handle(const Incident& incident){
 
  ATH_MSG_DEBUG ( " handle called  " << "for incident "<< incident. type() );
 
  m_newEvent = true ;
  return ;
}
 
double LArTBH6TriggerTimeTool::time(){
 
  if(m_fixed || !m_newEvent) return m_time;
 
  // new event, try to get it from FrontHit
  // double t1 = trackRecordTime();
  m_time = larTime() ;
 
  m_newEvent= false;
  return m_time;
}
 
 
double LArTBH6TriggerTimeTool::larTime(){
 
  double te = 0;
  double e = 0;
  int n=0;
 
  for(auto &it : m_hitcoll){
 
    for(const LArHit* hit : *it){
      e += hit->energy();
      te += hit->energy() * hit->time() ;
      ++n;
    }
  }
 
  if(n==0){
    ATH_MSG_INFO ( " no LArHit in this event" );
    return 0 ;
  }
  if(e==0){
    ATH_MSG_INFO ( " no LArHit energy in this event" );
    return 0;
  }
 
  double t = te/e;
  ATH_MSG_DEBUG ( " average time from LArHit =  " <<t );
 
  return t;
 
 
}