LArCalibDigitsAccumulatorFreeGain.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
*/
 
#include "LArCalibDigitsAccumulatorFreeGain.h"
 
#include "CLHEP/Units/SystemOfUnits.h"
#include <cmath>
#include <cstdint>
 
using CLHEP::ns;
 
LArCalibDigitsAccumulatorFreeGain::LArCalibDigitsAccumulatorFreeGain (const std::string& name, ISvcLocator* pSvcLocator):
  AthAlgorithm(name, pSvcLocator),
  m_onlineHelper(nullptr),
  m_delay(0),
  m_isPulsed(false),
  m_calibAccuDigitContainerName("LArAccumulatedCalibDigits"),
  m_nStepTrigger(1),
  m_delayScale(1*ns),
  m_keepPulsed(false)
{
  declareProperty("LArAccuCalibDigitContainerName",m_calibAccuDigitContainerName);
  declareProperty("KeyList",m_keylist);
  declareProperty("StepOfTriggers",m_nStepTrigger);
  declareProperty("DelayScale",m_delayScale);
  declareProperty("KeepOnlyPulsed",m_keepPulsed);
 
  m_delay=-1;
  m_event_counter=0;
}
 
 
StatusCode LArCalibDigitsAccumulatorFreeGain::initialize(){
  
  StatusCode sc;
 
  // retrieve online ID helper
  sc = detStore()->retrieve(m_onlineHelper, "LArOnlineID");
  if (sc.isFailure()) {
    ATH_MSG_ERROR( "Could not get LArOnlineID helper !" );
    return StatusCode::FAILURE;
  }
 
  ATH_CHECK( m_calibMapKey.initialize() );
 
  m_AccumulatedHG.resize(m_onlineHelper->channelHashMax());
  m_AccumulatedMG.resize(m_onlineHelper->channelHashMax());
  m_AccumulatedLG.resize(m_onlineHelper->channelHashMax());
 
  return StatusCode::SUCCESS;
}
 
 
 
StatusCode LArCalibDigitsAccumulatorFreeGain::execute() 
{
  
  StatusCode sc;
 
  if ( m_event_counter < 100 || m_event_counter%100==0 )
    ATH_MSG_INFO( "Processing event " << m_event_counter );
  ++m_event_counter;
 
  SG::ReadCondHandle<LArCalibLineMapping> clHdl{m_calibMapKey};
  const LArCalibLineMapping *clcabling {*clHdl};
  if(!clcabling) {
     ATH_MSG_WARNING( "Do not have calib line mapping !!!" );
     return StatusCode::FAILURE;
  }
  
  // pointer to input container
  const LArCalibDigitContainer* calibDigitContainer=nullptr;
 
  // retrieve calibration settings
  const LArCalibParams* calibParams;
  sc=detStore()->retrieve(calibParams,"LArCalibParams");
  if (sc.isFailure())
    {ATH_MSG_ERROR( "Cannot load LArCalibParams from DetStore." );
      return StatusCode::FAILURE;
    }
 
  unsigned int sizeSteps = (m_nStepTrigger>1 ? (m_nStepTrigger+1):1);
 
  // retrieve input calibDigits
  
  //Loop over all containers that are to be processed (e.g. different gains)
  for (const std::string& key : m_keylist) {
    
    sc=evtStore()->retrieve(calibDigitContainer,key);
    if(sc.isFailure()) {
      ATH_MSG_ERROR( "Can't retrieve LArCalibDigitContainer with key " << key << "from StoreGate." );
      return StatusCode::SUCCESS;
    }else{
      ATH_MSG_DEBUG( "Retrieved LArCalibDigitContainer with key " << key << " from StoreGate." );
    }
 
    // Loop over CalibDigitContainer
 
    if(calibDigitContainer->empty()) {
      ATH_MSG_DEBUG( "LArCalibDigitContainer with key=" << key << " is empty " );
    }else{
      ATH_MSG_DEBUG( "LArCalibDigitContainer with key=" << key << " has size =  " << calibDigitContainer->size() );
    }
 
    // counter of triggers 
    std::vector<unsigned int> ntrigger, nTriggerPerStep, nStepTrigger, iStepTrigger;
    // effective number of triggers: per channel
    ntrigger.resize(calibDigitContainer->size(),0);
    // asked number of triggers (from calib settings): per FEB
    nTriggerPerStep.resize(m_onlineHelper->febHashMax(),0);
    nStepTrigger.resize(m_onlineHelper->febHashMax(),0);
    iStepTrigger.resize(m_onlineHelper->febHashMax(),0);
 
    // output container
    LArAccumulatedCalibDigitContainer* larAccuCalibDigitContainerHG = new LArAccumulatedCalibDigitContainer();
 
    LArAccumulatedCalibDigitContainer* larAccuCalibDigitContainerMG = new LArAccumulatedCalibDigitContainer();
 
    LArAccumulatedCalibDigitContainer* larAccuCalibDigitContainerLG = new LArAccumulatedCalibDigitContainer();
 
    //Loop over all cells
    for (const LArCalibDigit* digit : *calibDigitContainer) {
 
      if(m_keepPulsed && !digit->isPulsed()) continue;
      
      // identificators
      HWIdentifier chid=digit->hardwareID();      
      const HWIdentifier febid=m_onlineHelper->feb_Id(chid);
      const IdentifierHash febhash = m_onlineHelper->feb_Hash(febid);
      const IdentifierHash hashid = m_onlineHelper->channel_Hash(chid);
      
      // BELOW: DIRTY HACK BECAUSE THERE SEEMS TO BE A BUG IN THE CABLINGSVC CONCERNING THE CALIBLINES.
 
      // get calibration settings
      const std::vector<HWIdentifier>& calibLineID=clcabling->calibSlotLine(chid);
      HWIdentifier calibModuleID;
      if(!calibLineID.empty()){
        calibModuleID=m_onlineHelper->calib_module_Id(calibLineID[0]);
        nTriggerPerStep[febhash] = calibParams->NTrigger(calibModuleID);
        ATH_MSG_DEBUG( "Ntrigger per step = " << nTriggerPerStep[febhash] );
        if(nTriggerPerStep[febhash] > 1000) nTriggerPerStep[febhash]=100; // very dirty !!! 
      } else {
        nTriggerPerStep[febhash] = 100; // very dirty !! 
      }
 
      // cell is pulsed ? 
      m_isPulsed = digit->isPulsed();
      
      //First cell to be processed, set delay
      if (m_delay==-1) { 
        m_delay=digit->delay();
      } else {
        // next cells: should be the same delay
        if (m_delay!=digit->delay()) {
          ATH_MSG_DEBUG( "Delay is changing to " << digit->delay() << " from " << m_delay << ": book a new LArAccumulatedCalibDigitContainer" );
          m_delay=digit->delay();
        }
      }
      
      CaloGain::CaloGain gain=digit->gain();
      if (gain<0 || gain>CaloGain::LARNGAIN) {
    ATH_MSG_ERROR( "Found not-matching gain number ("<< (int)gain <<")" );
        delete larAccuCalibDigitContainerHG;
        delete larAccuCalibDigitContainerMG;
        delete larAccuCalibDigitContainerLG;
    return StatusCode::FAILURE;
      }
 
      /***************************** High Gain ***************************/
      // object to be filled for each cell
      LArAccumulated& cellAccumulatedHG = m_AccumulatedHG[hashid];
      if (gain==CaloGain::LARHIGHGAIN) {
    // trigger counter for each cell
    cellAccumulatedHG.m_ntrigger++;
    ATH_MSG_DEBUG( "HG chid = " << chid << ", trigger = " << cellAccumulatedHG.m_ntrigger << ", DAC = " << digit->DAC() );
    // at first trigger, initialize vectors
    unsigned int sizeSamples = digit->samples().size();
    ATH_MSG_DEBUG( "sizeSteps = " << sizeSteps << ", # of samples = " << sizeSamples );
    if(cellAccumulatedHG.m_ntrigger==1){
      cellAccumulatedHG.m_sum.clear();
      cellAccumulatedHG.m_sum2.clear();
      cellAccumulatedHG.m_sum.resize(sizeSamples,0);
      cellAccumulatedHG.m_sum2.resize(sizeSamples,0);
    }
    for(unsigned int j=0;j<sizeSamples;j++) {
      cellAccumulatedHG.m_sum[j] += digit->samples()[j];
      cellAccumulatedHG.m_sum2[j] += digit->samples()[j]*digit->samples()[j];
    }
    ATH_MSG_DEBUG( "Sum = " << cellAccumulatedHG.m_sum[2] );
    ATH_MSG_DEBUG( "Sum2 = " << cellAccumulatedHG.m_sum2[2] );
      } // end High Gain
 
      /***************************** Medium Gain ***************************/
      // object to be filled for each cell
      LArAccumulated& cellAccumulatedMG = m_AccumulatedMG[hashid];
      if (gain==CaloGain::LARMEDIUMGAIN) {
    // trigger counter for each cell
    cellAccumulatedMG.m_ntrigger++;
    ATH_MSG_DEBUG( "MG chid = " << chid << ", trigger = " << cellAccumulatedMG.m_ntrigger << ", DAC = " << digit->DAC() );
    // at first trigger, initialize vectors
    unsigned int sizeSamples = digit->samples().size();
    ATH_MSG_DEBUG( "sizeSteps = " << sizeSteps << ", # of samples = " << sizeSamples );
    if(cellAccumulatedMG.m_ntrigger==1){
      cellAccumulatedMG.m_sum.clear();
      cellAccumulatedMG.m_sum2.clear();
      cellAccumulatedMG.m_sum.resize(sizeSamples,0);
      cellAccumulatedMG.m_sum2.resize(sizeSamples,0);
    }
    for(unsigned int j=0;j<sizeSamples;j++){
      cellAccumulatedMG.m_sum[j] += digit->samples()[j];
      cellAccumulatedMG.m_sum2[j] += digit->samples()[j]*digit->samples()[j];
    }
    ATH_MSG_DEBUG( "Sum = " << cellAccumulatedMG.m_sum[2] );
    ATH_MSG_DEBUG( "Sum2 = " << cellAccumulatedMG.m_sum2[2] );
      } // end High Gain
 
      /***************************** Low Gain ***************************/
      // object to be filled for each cell
      LArAccumulated& cellAccumulatedLG = m_AccumulatedLG[hashid];
      if (gain==CaloGain::LARLOWGAIN) {
    if (chid==959628800) std::cout << "---> in loop LG - DAC = " << digit->DAC() << std::endl;
    cellAccumulatedLG.m_ntrigger++;
    unsigned int sizeSamples = digit->samples().size();
    if(cellAccumulatedLG.m_ntrigger==1){
      cellAccumulatedLG.m_sum.clear();
      cellAccumulatedLG.m_sum2.clear();
      cellAccumulatedLG.m_sum.resize(sizeSamples,0);
      cellAccumulatedLG.m_sum2.resize(sizeSamples,0);
    }
    for(unsigned int j=0;j<sizeSamples;j++){
      cellAccumulatedLG.m_sum[j] += digit->samples()[j];
      cellAccumulatedLG.m_sum2[j] += digit->samples()[j]*digit->samples()[j];
    }
    ATH_MSG_DEBUG( "Sum = " << cellAccumulatedLG.m_sum[2] );
    ATH_MSG_DEBUG( "Sum2 = " << cellAccumulatedLG.m_sum2[2] );
      } // end High Gain
 
 
      // when reached total number of triggers for this step, fill LArAccumulatedCalibDigit and reset number of triggers
      if( (cellAccumulatedHG.m_ntrigger+cellAccumulatedMG.m_ntrigger+cellAccumulatedLG.m_ntrigger)==nTriggerPerStep[febhash]){
    ATH_MSG_DEBUG( "filling LArAccumulatedCalibDigit " );
    ATH_MSG_DEBUG( "chid = " << chid << ", gain = " << gain << ", DAC = " << digit->DAC() << ", isPulsed = " << m_isPulsed << ", delay = " << m_delay << ", trigPerStep = " << nTriggerPerStep[febhash] << ", istep = " << iStepTrigger[febhash] );
    iStepTrigger[febhash]++;
    unsigned int sizeSamples = digit->samples().size();
 
    if (cellAccumulatedHG.m_ntrigger>0) {
      LArAccumulatedCalibDigit* accuCalibDigitHG;
      std::vector < uint64_t > sampleSumHG;
      std::vector < uint64_t > sample2SumHG;
      sampleSumHG.resize(sizeSamples,0);
      sample2SumHG.resize(sizeSamples,0);
      accuCalibDigitHG = new LArAccumulatedCalibDigit(chid, gain, sampleSumHG, sample2SumHG, 0, (uint16_t)digit->DAC(), (uint16_t)m_delay, m_isPulsed, 0, 0);
      accuCalibDigitHG->setAddSubStep(cellAccumulatedHG.m_sum,cellAccumulatedHG.m_sum2,cellAccumulatedHG.m_ntrigger);
      larAccuCalibDigitContainerHG->push_back(accuCalibDigitHG);
    }
 
    if (cellAccumulatedMG.m_ntrigger>0) {
      LArAccumulatedCalibDigit* accuCalibDigitMG;
      std::vector < uint64_t > sampleSumMG;
      std::vector < uint64_t > sample2SumMG;
      sampleSumMG.resize(sizeSamples,0);
      sample2SumMG.resize(sizeSamples,0);
      accuCalibDigitMG = new LArAccumulatedCalibDigit(chid, gain, sampleSumMG, sample2SumMG, 0, (uint16_t)digit->DAC(), (uint16_t)m_delay, m_isPulsed, 0, 0);
      accuCalibDigitMG->setAddSubStep(cellAccumulatedMG.m_sum,cellAccumulatedMG.m_sum2,cellAccumulatedMG.m_ntrigger);
      larAccuCalibDigitContainerMG->push_back(accuCalibDigitMG);
    }
 
    if (cellAccumulatedLG.m_ntrigger>0) {
      LArAccumulatedCalibDigit* accuCalibDigitLG;
      std::vector < uint64_t > sampleSumLG;
      std::vector < uint64_t > sample2SumLG;
      sampleSumLG.resize(sizeSamples,0);
      sample2SumLG.resize(sizeSamples,0);
      accuCalibDigitLG = new LArAccumulatedCalibDigit(chid, gain, sampleSumLG, sample2SumLG, 0, (uint16_t)digit->DAC(), (uint16_t)m_delay, m_isPulsed, 0, 0);
      accuCalibDigitLG->setAddSubStep(cellAccumulatedLG.m_sum,cellAccumulatedLG.m_sum2,cellAccumulatedLG.m_ntrigger);
      larAccuCalibDigitContainerLG->push_back(accuCalibDigitLG);
    }
 
    cellAccumulatedHG.m_ntrigger = 0;
    cellAccumulatedMG.m_ntrigger = 0;
    cellAccumulatedLG.m_ntrigger = 0;   
      }
 
 
      
    }// loop over cells in container
 
 
    larAccuCalibDigitContainerHG->setDelayScale(m_delayScale);
    larAccuCalibDigitContainerMG->setDelayScale(m_delayScale);
    larAccuCalibDigitContainerLG->setDelayScale(m_delayScale);
 
    sc = evtStore()->record(larAccuCalibDigitContainerHG,"HIGH");
    if (sc!=StatusCode::SUCCESS)
      ATH_MSG_WARNING( "Unable to record LArAccumulatedCalibDigitContainer with key HIGH from DetectorStore. " );
    else
      ATH_MSG_DEBUG( "Recorded succesfully LArAccumulatedCalibDigitContainer with key HIGH" );
    sc = evtStore()->setConst(larAccuCalibDigitContainerHG);
    if (sc.isFailure()) {
      ATH_MSG_ERROR( " Cannot lock LArAccumulatedCalibDigitContainerHG " );
      delete larAccuCalibDigitContainerMG;
      delete larAccuCalibDigitContainerLG;
      return(StatusCode::FAILURE);
    }
 
    sc = evtStore()->record(larAccuCalibDigitContainerMG,"MEDIUM");
    if (sc!=StatusCode::SUCCESS) 
      ATH_MSG_WARNING( "Unable to record LArAccumulatedCalibDigitContainer with key MEDIUM from DetectorStore. " );
    else
      ATH_MSG_DEBUG( "Recorded succesfully LArAccumulatedCalibDigitContainer with key MEDIUM" );
    sc = evtStore()->setConst(larAccuCalibDigitContainerMG);
    if (sc.isFailure()) {
      ATH_MSG_ERROR( " Cannot lock LArAccumulatedCalibDigitContainerMG " );
      delete larAccuCalibDigitContainerLG;
      return(StatusCode::FAILURE);
    }
 
    sc = evtStore()->record(larAccuCalibDigitContainerLG,"LOW");
    if (sc!=StatusCode::SUCCESS)
      ATH_MSG_WARNING( "Unable to record LArAccumulatedCalibDigitContainer with key LOW from DetectorStore. " );
    else
      ATH_MSG_DEBUG( "Recorded succesfully LArAccumulatedCalibDigitContainer with key LOW" );
    sc = evtStore()->setConst(larAccuCalibDigitContainerLG);
    if (sc.isFailure()) {
      ATH_MSG_ERROR( " Cannot lock LArAccumulatedCalibDigitContainerLG " );
      return(StatusCode::FAILURE);
    }
  } // loop over key container
  return StatusCode::SUCCESS;
}