QwMollerADC_Channel.cc File Reference

Implementation for Moller ADC channel decoding and management. More...

#include "QwMollerADC_Channel.h"
#include <stdexcept>
#include "QwLog.h"
#include "QwUnits.h"
#include "QwBlinder.h"
#include "QwHistogramHelper.h"

Include dependency graph for QwMollerADC_Channel.cc:

Go to the source code of this file.

Functions
	__attribute__ ((no_sanitize("signed-integer-overflow"))) void QwMollerADC_Channel
std::ostream &	operator<< (std::ostream &stream, const QwMollerADC_Channel &channel)

Detailed Description

Implementation for Moller ADC channel decoding and management.

Definition in file QwMollerADC_Channel.cc.

Function Documentation

__attribute__()

__attribute__

(

(no_sanitize("signed-integer-overflow"))

)

Definition at line 370 of file QwMollerADC_Channel.cc.

                                         {
  fNumberOfSamples = fNumberOfSamples_map;
  fHardwareBlockSum_raw = 0;
//  Double_t hwsum_test = 0.0;
//  std::cout <<  "*******In QwMollerADC_Channel::SetRawEventData for channel:\t" << this->GetElementName() << std::endl;
  for (Int_t i = 0; i < fBlocksPerEvent; i++) 
    {
     Double_t block_raw = (fBlock[i] / fCalibrationFactor + fPedestal) * fNumberOfSamples / (fBlocksPerEvent * 1.0);
     if (std::abs(block_raw) >= pow(2,29)) {
      block_raw = std::copysign(pow(2,29)-1, block_raw);
      QwWarning << "QwMollerADC_Channel::SetRawEventData: Overflow in conversion to raw data for channel "
                << this->GetElementName() << ": ("
                << "fBlock[i] = " << fBlock[i] << " / "
                << "fCalibrationFactor = " << fCalibrationFactor << " + "
                << "fPedestal = " << fPedestal << ") * "
                << "fNumberOfSamples = " << fNumberOfSamples << " / "
                << "fBlocksPerEvent = " << fBlocksPerEvent << ". "
                << "Capping value to " << block_raw << "."
                << QwLog::endl;
     }
     fBlock_raw[i] = Int_t(block_raw);
     fHardwareBlockSum_raw += fBlock_raw[i];
     
    double_t block = fBlock[i] / fCalibrationFactor;
    double_t sigma = fMockGaussianSigma / fCalibrationFactor;
    fBlockSumSq_raw[i] = (sigma*sigma + block*block)*fNumberOfSamples_map / (fBlocksPerEvent * 1.0);
    fBlock_min[i] = (block - 3.0 * sigma) * double_t(fNumberOfSamples_map) / (fBlocksPerEvent * 1.0);
    fBlock_max[i] = (block + 3.0 * sigma) * double_t(fNumberOfSamples_map) / (fBlocksPerEvent * 1.0);
    
    fBlockSumSq_raw[4] += fBlockSumSq_raw[i];
    fBlock_min[4] = TMath::Min(fBlock_min[i],fBlock_min[4]);
    fBlock_max[4] = TMath::Max(fBlock_max[i],fBlock_max[4]);
    }
 
 
 
  fSoftwareBlockSum_raw = fHardwareBlockSum_raw;
 
  return;
}

References QwLog::endl(), and QwWarning.

Here is the call graph for this function:

operator<<()

std::ostream & operator<<	(	std::ostream &	stream,
		const QwMollerADC_Channel &	channel )

Definition at line 1736 of file QwMollerADC_Channel.cc.

{
  stream << channel.GetHardwareSum();
  return stream;
}