LArCaliWaveSubsetCnv_p3.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
#include "LArRawConditions/LArConditionsSubset.h"
#include "LArCondTPCnv/LArCaliWaveSubsetCnv_p3.h"
 
 
void
LArCaliWaveSubsetCnv_p3::persToTrans(const LArCaliWaveSubset_p3* persObj,  LArCWTransType* transObj, MsgStream & log) const
{
  log<<MSG::DEBUG<<"LArCaliWaveSubsetCNV_p3  begin persToTrans"<<endmsg;
 
  // Copy basic metadata
  transObj->setChannel       (persObj->m_subset.m_channel);
  transObj->setGroupingType  (persObj->m_subset.m_groupingType);
  transObj->initialize (persObj->m_subset.m_febIds, persObj->m_subset.m_gain);
  unsigned int nfebids = persObj->m_subset.m_febIds.size();         
  log<<MSG::DEBUG<<"Total febs:"<<nfebids;
  unsigned int ncorrs    = persObj->m_subset.m_corrChannels.size();     
  log<<MSG::DEBUG<<"\tCorrections:"<<ncorrs<<endmsg;
  const unsigned int nChannelsPerFeb  = persObj->m_subset.subsetSize();
  unsigned int cwvIndex=0, waveIndex=0, chIndex = 0;
   
  unsigned int ifebWithData = 0; // counter for febs with data
 
  auto subsetIt = transObj->subsetBegin();
  for (unsigned int i = 0; i < nfebids; ++i, ++subsetIt){
    unsigned int febid = subsetIt->first;
    log<<MSG::DEBUG<<"feb index: "<<i<<"\tfebID: "<<febid;
        
    bool hasSparseData       = false;
    unsigned int chansSet    = 0;
    unsigned int chansOffset = 0;
        
    if (ifebWithData < persObj->m_subset.m_febsWithSparseData.size() && febid == persObj->m_subset.m_febsWithSparseData[ifebWithData]) { 
      hasSparseData = true;// Found feb with sparse data
      ifebWithData++;
      chansSet  = persObj->m_subset.m_febsWithSparseData[ifebWithData];
      chansOffset = 0;
      ifebWithData++;
    }
      
    log<<MSG::DEBUG<<"\tFeb sparse? "<< hasSparseData <<endmsg;
  
    for (unsigned int j = 0; j < nChannelsPerFeb; ++j){
      bool copyChannel = true;
      if (hasSparseData) {            
        if (!(chansSet & (1 << (j - chansOffset)))) {   // Channel is missing data - skip
          copyChannel = false;
        }
        if (j%32 == 31 && j < nChannelsPerFeb-2) {
          chansSet     = persObj->m_subset.m_febsWithSparseData[ifebWithData];
          chansOffset += 32;
          ifebWithData++;
        }
      }
                      
      if (copyChannel) { // This channel has wave, resize vectors     
        LArCaliWaveVec& CWV = subsetIt->second[j]; 
                   
        unsigned int nDAC=persObj->m_vDAC[cwvIndex];
        cwvIndex++;
        for (unsigned int nD=0;nD<nDAC;nD++){
          double time       =persObj->m_dt[chIndex];
          unsigned int f  =persObj->m_flag[chIndex];
          int dac           =persObj->m_DAC[chIndex];
          int pulsed            =persObj->m_isPulsed[chIndex];
          chIndex++;
                
          std::vector<double> val, err;
          std::vector<int>  tri;
        
          for (unsigned int k = 0; k < persObj->m_samples; ++k){
            val.push_back(persObj->m_vAmplitudes[waveIndex]);
            err.push_back(persObj->m_vErrors[waveIndex]);
            tri.push_back(persObj->m_vTriggers[waveIndex]);
            waveIndex++;
          }   
                
          LArCaliWave cv(val,err,tri, time,dac, pulsed, f);
          val.clear();
          err.clear();
          tri.clear();
          
          CWV.push_back(cv);
        }
      }
    } // loop on channels
  } // loop on FEBs
    
 
  log<<MSG::DEBUG << "Correction reading, ncorr:" << ncorrs << "\twaveIndex=" <<  waveIndex << "\tsize:" << persObj->m_vAmplitudes.size()<<endmsg;
  LArCWTransType::CorrectionVec corrs;
  corrs.resize(ncorrs);
   
  // Loop over corrections
  for (unsigned int i = 0; i < ncorrs; ++i){
    if (cwvIndex    >= persObj->m_vDAC.size()) {// check indexes
      log << MSG::ERROR 
          << "LArCaliWaveSubsetCnv_p3::persToTrans - CaliWave index too large: cwvIndex/sizeInFile " 
          << cwvIndex << " " << persObj->m_vDAC.size() << " " << endmsg;
      return;
    }
        
        
    corrs[i].first = persObj->m_subset.m_corrChannels[i];
    
    LArCaliWaveVec& CWV = corrs[i].second;
        
    unsigned int nDAC=persObj->m_vDAC[cwvIndex];
    cwvIndex++;
    for (unsigned int nD=0;nD<nDAC;nD++){   
      if (chIndex    >= persObj->m_dt.size()) {// check indexes
        log << MSG::ERROR 
            << "LArCaliWaveSubsetCnv_p3::persToTrans - CaliWave index too large: WaveIndex/sizeInFile " 
            << chIndex << " " << persObj->m_dt.size() << " " << endmsg;
        return;
      }
        
      double time   = persObj->m_dt[chIndex];
      unsigned int f  = persObj->m_flag[chIndex];
      int dac           = persObj->m_DAC[chIndex];
      int pulsed        = persObj->m_isPulsed[chIndex];
      chIndex++;
      std::vector<double> val, err;
      std::vector<int>  tri;
      for (unsigned int k = 0; k < persObj->m_samples; ++k){
            
        if (waveIndex>=persObj->m_vAmplitudes.size()) 
          log << MSG::ERROR << "Persistent LArCaliWave object is inconsistent. i=" << i << " WaveIndes=" 
              << waveIndex << " size=" << persObj->m_vAmplitudes.size() << "samples=" << persObj->m_samples <<endmsg;
        
        val.push_back(persObj->m_vAmplitudes[waveIndex]);
        err.push_back(persObj->m_vErrors[waveIndex]);
        tri.push_back(persObj->m_vTriggers[waveIndex]);
        waveIndex++;
      }     
    
      LArCaliWave cv(val,err,tri, time,dac, pulsed, f);
      val.clear();
      err.clear();
      tri.clear();
                    
      CWV.push_back(cv);
            
    }// over DAC values
  }//over corrections
  transObj->insertCorrections (std::move (corrs));
    
 
  log<< MSG::DEBUG <<"CaliWave ver p2 successfully read in."<<endmsg;
}
 
 
 
 
 
 
void
LArCaliWaveSubsetCnv_p3::transToPers(const LArCWTransType* transObj,  LArCaliWaveSubset_p3* persObj, MsgStream & log) const
{
  log<<MSG::DEBUG<<" LArCaliWaveSubsetCNV_p3  begin Writing"<<endmsg;
    
  unsigned int nfebs        = transObj->subsetSize();
  log<<MSG::DEBUG<<"total febs:"<<nfebs;
  unsigned int ncorrs       = transObj->correctionVecSize();
  log<<MSG::DEBUG<<"\tCorrections: "<<ncorrs<<endmsg;
  const unsigned int nChannelsPerFeb  = transObj->channelVectorSize();
 
  unsigned int nsubsetsNotEmpty = 0;
  unsigned int nchans         = 0;
 
  persObj->m_samples=0;
 
  std::vector<unsigned int> febsWithSparseData; // collects IDs of febs with sparse data
    
  const auto subsetEnd = transObj->subsetEnd();
  for (auto subsetIt = transObj->subsetBegin();
       subsetIt != subsetEnd;
       ++subsetIt)
  {
    unsigned int nfebChans = subsetIt->second.size();
 
    if (nfebChans != 0 && nfebChans != nChannelsPerFeb) {
      log << MSG::ERROR << "LArCaliWaveSubsetCnv_p3::transToPers - found incorrect number of channels per feb: " << nfebChans<< endmsg;
      return;
    }
        
    if (nfebChans) ++nsubsetsNotEmpty; // count number of non-empty FEBS
                
    bool isSparse = false;
    for (unsigned int j = 0; j < nfebChans; ++j) { // Loop over channels and check if this subset has sparse data
            
      const LArCaliWaveVec& CWV = subsetIt->second[j];
            
      if (!isSparse && CWV.size() == 0) {
        isSparse = true;
        febsWithSparseData.push_back(subsetIt->first);// save febids for sparse subsets
      }
            
      if (CWV.size()){  // should be like this no?
        //else { 
        nchans++; // count number of channels
        persObj->m_samples=CWV[0].getSize(); // completely ineficient but fast to fix :-(
      }
            
    }
            
  } // loop over febs
    
  persObj->m_subset.m_febIds.reserve(nsubsetsNotEmpty);
    
  persObj->m_subset.m_corrChannels.reserve(ncorrs);
    
  if (febsWithSparseData.size())
    persObj->m_subset.m_febsWithSparseData.reserve(febsWithSparseData.size()*5); // this stores FEBID + 128 bitmask of existing channels
  
  persObj->m_dt.reserve(nchans);
  persObj->m_vAmplitudes.reserve(nchans*persObj->m_samples);
  persObj->m_vErrors.reserve(nchans*persObj->m_samples);
  persObj->m_vTriggers.reserve(nchans*persObj->m_samples);
  persObj->m_flag.reserve(nchans);
    
  int counterSparse=0;
  for (auto subsetIt = transObj->subsetBegin();
       subsetIt != subsetEnd;
       ++subsetIt)
  {
    unsigned int nfebChans = subsetIt->second.size();
        
    if (nfebChans == 0) continue; // skip subsets without any channels
        
    unsigned int febid = subsetIt->first;
    persObj->m_subset.m_febIds.push_back(febid);   // stores FEBID's of nonempty FEBs
        
    bool isSparse=false; // this is just to avoid search loop
    if (counterSparse < (int)febsWithSparseData.size() && febid == febsWithSparseData[counterSparse]){
      counterSparse++;
      isSparse=true;
      persObj->m_subset.m_febsWithSparseData.push_back(febid); 
    }
            
    unsigned int chansSet    = 0;
    unsigned int chansOffset = 0;
    for (unsigned int j = 0; j < nfebChans; ++j){
 
      bool saveAmplitudes=true;
      if (isSparse) { // subset with sparse data
                
        if (subsetIt->second[j].size() > 0) { // channel exists
                    
          assert (j >= chansOffset && (j - chansOffset) <= 31);
          chansSet |= (1 << (j - chansOffset)); //store the channel number in the bit map
                    
        }
        else { // channel does not exist
          saveAmplitudes = false;
        }
                
        // Save chansSet
        if  (j == (chansOffset + 31) || j == nfebChans-1 ) {
          persObj->m_subset.m_febsWithSparseData.push_back(chansSet);
          chansSet    =   0;
          chansOffset += 32;
        }
      } // if sparse
            
      if (saveAmplitudes) {             // save amplitudes, errors and triggers
        const LArCaliWaveVec& CWV = subsetIt->second[j];
        int dacValues=CWV.size();
        persObj->m_vDAC.push_back(dacValues);   
                
        for(int dv=0;dv<dacValues;dv++){
          persObj->m_dt.push_back(CWV[dv].getDt());
          persObj->m_flag.push_back(CWV[dv].getFlag());
          persObj->m_DAC.push_back(CWV[dv].getDAC());
          persObj->m_isPulsed.push_back(CWV[dv].getIsPulsedInt());
          const std::vector<double>& w=CWV[dv].getWave();
          const std::vector<double>& e=CWV[dv].getErrors();
          const std::vector<int>& t=CWV[dv].getTriggers();
          for (unsigned int k = 0; k< w.size(); ++k){
            persObj->m_vAmplitudes.push_back(w[k]);
            persObj->m_vErrors.push_back(e[k]);
            persObj->m_vTriggers.push_back(t[k]);
          } // over samples
        }// over DAC values
      }
            
    } // over channels
  }// over febs
 
  // Copy corrections
  const auto corrEnd = transObj->correctionVecEnd();
  for (auto corrIt = transObj->correctionVecBegin();
       corrIt != corrEnd;
       ++corrIt)
  {
    persObj->m_subset.m_corrChannels.push_back(corrIt->first);// Save channel id in febid vector
    
    const LArCaliWaveVec& CWV = corrIt->second; 
        
    int dacValues=CWV.size();
    persObj->m_vDAC.push_back(dacValues);    
    log<<MSG::DEBUG<<"\t CWV size: "<<dacValues;    
        
    for(int dv=0;dv<dacValues;dv++){
      persObj->m_dt.push_back(CWV[dv].getDt());
      persObj->m_flag.push_back(CWV[dv].getFlag());
      persObj->m_DAC.push_back(CWV[dv].getDAC());
      persObj->m_isPulsed.push_back(CWV[dv].getIsPulsedInt());
      const std::vector<double>& w=CWV[dv].getWave();
      const std::vector<double>& e=CWV[dv].getErrors();
      const std::vector<int>&   t=CWV[dv].getTriggers();
      persObj->m_samples=w.size();
      for (unsigned int k = 0; k< w.size(); ++k){
        persObj->m_vAmplitudes.push_back(w[k]);
        persObj->m_vErrors.push_back(e[k]);
        persObj->m_vTriggers.push_back(t[k]);
      }// over samples
    }// over dac values
  }// over corrections
 
  // Copy the rest
  persObj->m_subset.m_gain          = transObj->gain(); 
  persObj->m_subset.m_channel       = transObj->channel();
  persObj->m_subset.m_groupingType  = transObj->groupingType();
  log<< MSG::DEBUG <<" Cali Wave p2 successfully written down."<<endmsg;
}