#include <SiStripFEDEmulator.h>

Public Member Functions
	FEDEmulator ()

uint32_t	fedIndex (const uint16_t aFedChannel)

void	fillMedians (const std::map< uint32_t, std::vector< uint32_t > >::const_iterator &meds)

void	fillNoises (const edm::DetSetVector< SiStripProcessedRawDigi >::const_iterator &noise)

void	fillPeds (const edm::DetSetVector< SiStripRawDigi >::const_iterator &peds)

void	initialise (const bool byModule)

void	initialiseModule (const uint32_t aDetId, const uint32_t aNPairs, const uint32_t aPair)

void	print (std::ostream &aOs)

void	printMeds (std::ostream &aOs)

void	printNoises (std::ostream &aOs)

void	printPeds (std::ostream &aOs)

void	retrieveNoises (const edm::ESHandle< SiStripNoises > &aHandle)

void	retrievePedestals (const edm::ESHandle< SiStripPedestals > &aHandle)

void	subtractCM (const std::vector< SiStripRawDigi > &pedSubtrDetSetData, std::vector< SiStripRawDigi > &cmSubtrDetSetData)

void	subtractPedestals (const edm::DetSetVector< SiStripRawDigi >::const_iterator &inputChannel, std::vector< SiStripRawDigi > &pedsDetSetData, std::vector< SiStripProcessedRawDigi > &noiseDetSetData, std::vector< SiStripRawDigi > &pedSubtrDetSetData, std::vector< uint32_t > &medsDetSetData, const bool fillApvsForCM)

void	zeroSuppress (const std::vector< SiStripRawDigi > &cmSubtrDetSetData, edm::DetSet< SiStripDigi > &zsDetSetData, const std::unique_ptr< SiStripRawProcessingAlgorithms > &algorithms)

	~FEDEmulator ()

Private Attributes
bool	byModule_

uint32_t	detId_

uint32_t	maxStrip_

std::vector< uint32_t >	medians_

uint32_t	minStrip_

std::vector< float >	noises_

uint32_t	nPairs_

uint32_t	pair_

std::vector< int >	pedestals_

Static Private Attributes
static const char *const	messageLabel_ = "SiStripFEDEmulator"

Detailed Description

Definition at line 35 of file SiStripFEDEmulator.h.

Constructor & Destructor Documentation

sistrip::FEDEmulator::FEDEmulator ( )

Definition at line 15 of file SiStripFEDEmulator.cc.

References byModule_, detId_, maxStrip_, medians_, minStrip_, noises_, nPairs_, pair_, and pedestals_.

                            {
     byModule_ = false;
     minStrip_ = 0;
     maxStrip_ = 0;
     pedestals_.clear();
     noises_.clear();
     medians_.clear();
 
     detId_ = 0;
     nPairs_ = 0;
     pair_ = 0;
   }

sistrip::FEDEmulator::~FEDEmulator ( )

Definition at line 28 of file SiStripFEDEmulator.cc.

28 {}

Member Function Documentation

uint32_t sistrip::FEDEmulator::fedIndex ( const uint16_t aFedChannel )

void sistrip::FEDEmulator::fillMedians ( const std::map< uint32_t, std::vector< uint32_t > >::const_iterator & meds )

Definition at line 215 of file SiStripFEDEmulator.cc.

References SplitLinear::begin, dataset::end, and medians_.

                                                                                                 {
     std::vector<uint32_t>::const_iterator iApv = (meds->second).begin();
     for (; iApv != (meds->second).end(); ++iApv) {
       medians_.push_back(*iApv);
     }
   }

void sistrip::FEDEmulator::fillNoises ( const edm::DetSetVector< SiStripProcessedRawDigi >::const_iterator & noise )

Definition at line 209 of file SiStripFEDEmulator.cc.

References edm::DetSetVector< T >::begin(), edm::DetSetVector< T >::end(), and noises_.

                                                                                                    {
     for (edm::DetSet<SiStripProcessedRawDigi>::const_iterator iApv = noises->begin(); iApv != noises->end(); ++iApv) {
       noises_.push_back(iApv->adc());
     }
   }

void sistrip::FEDEmulator::fillPeds ( const edm::DetSetVector< SiStripRawDigi >::const_iterator & peds )

Definition at line 203 of file SiStripFEDEmulator.cc.

References edm::DetSetVector< T >::begin(), edm::DetSetVector< T >::end(), and pedestals_.

                                                                                       {
     for (edm::DetSet<SiStripRawDigi>::const_iterator iApv = peds->begin(); iApv != peds->end(); ++iApv) {
       pedestals_.push_back(iApv->adc());
     }
   }

void sistrip::FEDEmulator::initialise ( const bool byModule )

Definition at line 30 of file SiStripFEDEmulator.cc.

References MeasureLA_cff::byModule, byModule_, medians_, noises_, and pedestals_.

Referenced by sistrip::FEDEmulatorModule::FEDEmulatorModule().

                                                   {
     byModule_ = byModule;
     pedestals_.clear();
     noises_.clear();
     medians_.clear();
   }

void sistrip::FEDEmulator::initialiseModule	(	const uint32_t	aDetId,
		const uint32_t	aNPairs,
		const uint32_t	aPair
	)

Definition at line 37 of file SiStripFEDEmulator.cc.

References byModule_, detId_, maxStrip_, medians_, minStrip_, noises_, nPairs_, pair_, pedestals_, and sistrip::STRIPS_PER_FEDCH.

Referenced by sistrip::FEDEmulatorModule::produce().

                                                                                                         {
     detId_ = aDetId;
     nPairs_ = aNPairs;
     pair_ = aPair;
 
     minStrip_ = 0;
     maxStrip_ = nPairs_ * sistrip::STRIPS_PER_FEDCH;
     if (!byModule_) {
       minStrip_ = sistrip::STRIPS_PER_FEDCH * pair_;
       maxStrip_ = sistrip::STRIPS_PER_FEDCH * (pair_ + 1);
     }
 
     //need resize because it will not be filled by push_back...
     pedestals_.clear();
     pedestals_.resize(nPairs_ * sistrip::STRIPS_PER_FEDCH);
 
     noises_.clear();
     noises_.resize(nPairs_ * sistrip::STRIPS_PER_FEDCH);
 
     medians_.clear();
     medians_.reserve(nPairs_ * 2);
   }

void sistrip::FEDEmulator::print ( std::ostream & aOs )

Definition at line 222 of file SiStripFEDEmulator.cc.

References byModule_, detId_, maxStrip_, medians_, minStrip_, noises_, nPairs_, pair_, and pedestals_.

                                          {
     aOs << "===============================" << std::endl
         << " ===== FEDEmulator::print =====" << std::endl
         << " === byModule = " << byModule_ << std::endl
         << " === minStrip = " << minStrip_ << std::endl
         << " === maxStrip = " << maxStrip_ << std::endl
         << " === size of pedestals = " << pedestals_.size() << std::endl
         << " === size of noises = " << noises_.size() << std::endl
         << " === size of medians = " << medians_.size() << std::endl
         << " === detId = " << detId_ << std::endl
         << " === nPairs = " << nPairs_ << std::endl
         << " === pair = " << pair_ << std::endl
         << "===============================" << std::endl;
   }

void sistrip::FEDEmulator::printMeds ( std::ostream & aOs )

Definition at line 259 of file SiStripFEDEmulator.cc.

References mps_fire::i, and medians_.

                                              {
     aOs << "=========================================" << std::endl
         << " ===== FEDEmulator::printMedians =====" << std::endl;
     for (unsigned int i(0); i < medians_.size(); i++) {
       if (i % 32 == 0)
         aOs << std::endl;
       aOs << medians_[i] << " ";
     }
     aOs << "=========================================" << std::endl;
   }

void sistrip::FEDEmulator::printNoises ( std::ostream & aOs )

Definition at line 248 of file SiStripFEDEmulator.cc.

References mps_fire::i, and noises_.

                                                {
     aOs << "=========================================" << std::endl
         << " ===== FEDEmulator::printNoises =====" << std::endl;
     for (unsigned int i(0); i < noises_.size(); i++) {
       if (i % 32 == 0)
         aOs << std::endl;
       aOs << noises_[i] << " ";
     }
     aOs << "=========================================" << std::endl;
   }

void sistrip::FEDEmulator::printPeds ( std::ostream & aOs )

Definition at line 237 of file SiStripFEDEmulator.cc.

References mps_fire::i, and pedestals_.

                                              {
     aOs << "=========================================" << std::endl
         << " ===== FEDEmulator::printPedestals =====" << std::endl;
     for (unsigned int i(0); i < pedestals_.size(); i++) {
       if (i % 32 == 0)
         aOs << std::endl;
       aOs << pedestals_[i] << " ";
     }
     aOs << "=========================================" << std::endl;
   }

void sistrip::FEDEmulator::retrieveNoises ( const edm::ESHandle< SiStripNoises > & aHandle )

Definition at line 72 of file SiStripFEDEmulator.cc.

References detId_, alignCSCRings::e, messageLabel_, noises_, and cms::Exception::what().

Referenced by sistrip::FEDEmulatorModule::produce().

                                                                             {
     try {
       SiStripNoises::Range noisesRange = aHandle->getRange(detId_);
       aHandle->allNoises(noises_, noisesRange);
     } catch (const cms::Exception& e) {
       LogError(messageLabel_) << " Something wrong whilst getting noises: size = " << noises_.size()
                               << " detId = " << detId_ << std::endl
                               << e.what();
       //throw cms::Exception("ProductNotFound") << "Problem with noises..." ;
     }
   }  //end of retrieveNoises method.

void sistrip::FEDEmulator::retrievePedestals ( const edm::ESHandle< SiStripPedestals > & aHandle )

Definition at line 60 of file SiStripFEDEmulator.cc.

References detId_, alignCSCRings::e, messageLabel_, pedestals_, and cms::Exception::what().

Referenced by sistrip::FEDEmulatorModule::produce().

                                                                                   {
     try {
       SiStripPedestals::Range pedestalsRange = aHandle->getRange(detId_);
       aHandle->allPeds(pedestals_, pedestalsRange);
     } catch (const cms::Exception& e) {
       LogError(messageLabel_) << " Something wrong whilst getting pedestals: size = " << pedestals_.size()
                               << " detId = " << detId_ << std::endl
                               << e.what();
       //throw cms::Exception("ProductNotFound") << "Problem with pedestals..." ;
     }
   }  //end of retrievePedestals method.

void sistrip::FEDEmulator::subtractCM	(	const std::vector< SiStripRawDigi > &	pedSubtrDetSetData,
		std::vector< SiStripRawDigi > &	cmSubtrDetSetData
	)

Definition at line 169 of file SiStripFEDEmulator.cc.

References maxStrip_, medians_, minStrip_, sistrip::STRIPS_PER_APV, and relativeConstraints::value.

Referenced by sistrip::FEDEmulatorModule::produce().

                                                                              {
     //subtract Medians
     std::vector<SiStripRawDigi>::const_iterator lDigi = pedSubtrDetSetData.begin();
     uint32_t lCount = 0;
 
     for (uint32_t st = minStrip_; st < maxStrip_; st++) {
       uint32_t lApvNum = static_cast<uint32_t>(lCount * 1. / sistrip::STRIPS_PER_APV);
 
       int value = lDigi->adc() - medians_.at(lApvNum);
       if (value < 0) {
         //FED doesn't handle negative values
         value = 0;
       }
       cmSubtrDetSetData.push_back(SiStripRawDigi(static_cast<uint16_t>(value)));
 
       ++lDigi;
       lCount++;
     }
   }

void sistrip::FEDEmulator::subtractPedestals	(	const edm::DetSetVector< SiStripRawDigi >::const_iterator &	inputChannel,
		std::vector< SiStripRawDigi > &	pedsDetSetData,
		std::vector< SiStripProcessedRawDigi > &	noiseDetSetData,
		std::vector< SiStripRawDigi > &	pedSubtrDetSetData,
		std::vector< uint32_t > &	medsDetSetData,
		const bool	fillApvsForCM
	)

Definition at line 84 of file SiStripFEDEmulator.cc.

References edm::DetSetVector< T >::begin(), detId_, edm::DetSetVector< T >::end(), LogDebug, maxStrip_, medians_, messageLabel_, minStrip_, gpuVertexFinder::noise, noises_, nPairs_, pedestals_, edm::DetSet< T >::push_back(), and sistrip::STRIPS_PER_APV.

Referenced by sistrip::FEDEmulatorModule::produce().

                                                                 {
     edm::DetSet<SiStripRawDigi>::const_iterator digi_it = inputChannel->begin();
     uint32_t lCount = 0;
     std::vector<std::vector<uint16_t> > apvs;
     apvs.reserve(nPairs_ * 2);
 
     //subtract pedestals and fill apv vector
     //for median calculation
     //in a module, 256*nPairs, but here we are dealing
     //channel by channel, so 256 values at a time in the pedestal object.
     //reordered = 0-127 APV0 128-255 APV1
 
     std::vector<uint16_t> apvVec;
     ;
     apvVec.reserve(sistrip::STRIPS_PER_APV);
 
     for (uint32_t st = minStrip_; st < maxStrip_; st++) {
       uint16_t ped = static_cast<uint16_t>(pedestals_[st]);
       pedsDetSetData.push_back(SiStripRawDigi(ped));
 
       float noise = noises_[st];
       noiseDetSetData.push_back(SiStripProcessedRawDigi(noise));
 
       if (digi_it == inputChannel->end()) {
         LogError(messageLabel_) << " Error, end of inputchannel reached for detid " << detId_ << "! Processed "
                                 << lCount << " strips." << std::endl;
         break;
       }
 
       int lVal = digi_it->adc() - ped;
       if (lVal < 0) {
         if (digi_it->adc() > 0)
           LogDebug(messageLabel_) << " Negative value after pedestal subtraction ! Detid = " << detId_ << ", strip "
                                   << st << ", digi = " << digi_it->adc() << ", ped = " << ped << std::endl;
         lVal = 0;
       }
       pedSubtrDetSetData.push_back(SiStripRawDigi(static_cast<uint16_t>(lVal)));
 
       if (fillApvsForCM) {
         apvVec.push_back(static_cast<uint16_t>(lVal));
       }
 
       lCount++;
       ++digi_it;
 
       if (lCount % sistrip::STRIPS_PER_APV == 0 && fillApvsForCM) {
         if (!apvVec.empty())
           apvs.push_back(apvVec);
         apvVec.clear();
         apvVec.reserve(sistrip::STRIPS_PER_APV);
       }
 
     }  //end of loop over strips
 
     if (fillApvsForCM) {
       if (apvs.size() != nPairs_ * 2) {
         LogError(messageLabel_) << " Error, wrong output size for median calculation for detid " << detId_
                                 << "! Processed " << apvs.size() << " APVs." << std::endl;
       }
 
       //calculate common mode values
       for (uint32_t iapv(0); iapv < apvs.size(); iapv++) {
         std::vector<uint16_t> lVec = apvs[iapv];
 
         std::vector<uint16_t>::iterator mid = lVec.begin() + (lVec.size() / 2 - 1);
         std::nth_element(lVec.begin(), mid, lVec.end());
 
         //std::sort(lVec.begin(),lVec.end());
         //get the 64th element
         //if (*mid != lVec.at(63)) exit(0);
         //medians_.push_back(lVec.at(63));
         //medsDetSetData.push_back(lVec.at(63));
         medians_.push_back(*mid);
         medsDetSetData.push_back(*mid);
       }
     }
 
   }  //end of FEDEmulator::subtractPedestals method.

void sistrip::FEDEmulator::zeroSuppress	(	const std::vector< SiStripRawDigi > &	cmSubtrDetSetData,
		edm::DetSet< SiStripDigi > &	zsDetSetData,
		const std::unique_ptr< SiStripRawProcessingAlgorithms > &	algorithms
	)

Definition at line 190 of file SiStripFEDEmulator.cc.

References SiStripRawDigi::adc(), and HcalDetIdTransform::transform().

Referenced by sistrip::FEDEmulatorModule::produce().

                                                                                                   {
     //transform the input digis to a vector of integers
     std::vector<int16_t> cmSubtrRawDigis;
     transform(cmSubtrDetSetData.begin(),
               cmSubtrDetSetData.end(),
               back_inserter(cmSubtrRawDigis),
               std::bind(&SiStripRawDigi::adc, std::placeholders::_1));
     algorithms->suppressor->suppress(cmSubtrRawDigis, 0, zsDetSetData);
 
   }  //end of FEDEmulator::zeroSuppress method.