Internally used data structure. More...

Collaboration diagram for LArNoiseCorrelationMonAlg::perFeb_t:

Public Member Functions
	perFeb_t (const std::string &n)
void	sumSamples (const LArOnlineID *lArOnlineIDHelper)

Public Attributes
std::string	m_febName
std::vector< std::pair< const LArDigit *, double > >	m_digitsAndPed
std::vector< std::pair< int, double > >	m_meanSum
std::vector< std::pair< std::pair< int, int >, double > >	m_partSum

Detailed Description

Internally used data structure.

Definition at line 93 of file LArNoiseCorrelationMonAlg.h.

Constructor & Destructor Documentation

◆ perFeb_t()

LArNoiseCorrelationMonAlg::perFeb_t::perFeb_t ( const std::string & n )

inline

Definition at line 94 of file LArNoiseCorrelationMonAlg.h.

94: m_febName(n) {};

LArNoiseCorrelationMonAlg::perFeb_t::m_febName

std::string m_febName

Definition LArNoiseCorrelationMonAlg.h:95

Member Function Documentation

◆ sumSamples()

void LArNoiseCorrelationMonAlg::perFeb_t::sumSamples ( const LArOnlineID * lArOnlineIDHelper )

Definition at line 207 of file LArNoiseCorrelationMonAlg.cxx.

                                                                                       {
  const size_t S=m_digitsAndPed.size();
  for (size_t i1=0;i1<S;++i1) {
    const auto& [pDig1,pedestal1]=m_digitsAndPed[i1];
    const int ch1 = lArOnlineIDHelper->channel(pDig1->channelID());
    //Sum mean:
    for (const short adc : pDig1->samples()) {
      m_meanSum.emplace_back(ch1,adc-pedestal1);
    }  
    //sum of squares
    for (size_t i2=i1;i2<S;++i2) {
      const auto& [pDig2,pedestal2]=m_digitsAndPed[i2];
      const int ch2 = lArOnlineIDHelper->channel(pDig2->channelID());
      double sumSquare=0;
      const unsigned nADC=std::min(pDig1->nsamples(),pDig2->nsamples());    
      for (unsigned i=0;i<nADC;++i) {
        sumSquare+=((pDig1->samples().at(i)-pedestal1)*(pDig2->samples().at(i)-pedestal2));
      }//end loop over samples   
     m_partSum.emplace_back(std::make_pair(std::make_pair(ch1,ch2),sumSquare));
    }//end loop pDig2
  }//end loop pDig1 
}