#include <SiStripRawProcessingAlgorithms.h>

Public Types
using	digivector_t = SiStripAPVRestorer::digivector_t

Public Member Functions
uint16_t	convertVirginRawToHybrid (const edm::DetSet< SiStripRawDigi > &rawDigis, edm::DetSet< SiStripDigi > &suppressedDigis)

uint16_t	convertVirginRawToHybrid (uint32_t detId, uint16_t firstAPV, digivector_t &inDigis, edm::DetSet< SiStripDigi > &rawDigis)

const std::vector< bool > &	getAPVFlags () const

const SiStripAPVRestorer::medians_t &	getAPVsCM () const

const SiStripAPVRestorer::baselinemap_t &	getBaselineMap () const

const std::map< uint16_t, SiStripAPVRestorer::digimap_t > &	getSmoothedPoints () const

void	initialize (const edm::EventSetup &)

void	initialize (const edm::EventSetup &, const edm::Event &)

uint16_t	suppressHybridData (const edm::DetSet< SiStripDigi > &inDigis, edm::DetSet< SiStripDigi > &suppressedDigis, uint16_t firstAPV=0)

uint16_t	suppressProcessedRawData (const edm::DetSet< SiStripRawDigi > &rawDigis, edm::DetSet< SiStripDigi > &output)

uint16_t	suppressProcessedRawData (uint32_t detId, uint16_t firstAPV, digivector_t &procRawDigis, edm::DetSet< SiStripDigi > &output)

uint16_t	suppressVirginRawData (const edm::DetSet< SiStripRawDigi > &rawDigis, edm::DetSet< SiStripDigi > &output)

uint16_t	suppressVirginRawData (uint32_t detId, uint16_t firstAPV, digivector_t &procRawDigis, edm::DetSet< SiStripDigi > &output)

Public Attributes
const std::unique_ptr< SiStripAPVRestorer >	restorer

const std::unique_ptr< SiStripCommonModeNoiseSubtractor >	subtractorCMN

const std::unique_ptr< SiStripPedestalsSubtractor >	subtractorPed

const std::unique_ptr< SiStripFedZeroSuppression >	suppressor

Private Member Functions
	SiStripRawProcessingAlgorithms (std::unique_ptr< SiStripPedestalsSubtractor > ped, std::unique_ptr< SiStripCommonModeNoiseSubtractor > cmn, std::unique_ptr< SiStripFedZeroSuppression > zs, std::unique_ptr< SiStripAPVRestorer > res, bool doAPVRest, bool useCMMap)

Private Attributes
const bool	doAPVRestore

const TrackerGeometry *	trGeo

const bool	useCMMeanMap

Friends
class	SiStripRawProcessingFactory

Detailed Description

Definition at line 15 of file SiStripRawProcessingAlgorithms.h.

Member Typedef Documentation

◆ digivector_t

using SiStripRawProcessingAlgorithms::digivector_t = SiStripAPVRestorer::digivector_t

Definition at line 19 of file SiStripRawProcessingAlgorithms.h.

Constructor & Destructor Documentation

◆ SiStripRawProcessingAlgorithms()

SiStripRawProcessingAlgorithms::SiStripRawProcessingAlgorithms	(	std::unique_ptr< SiStripPedestalsSubtractor >	ped,
		std::unique_ptr< SiStripCommonModeNoiseSubtractor >	cmn,
		std::unique_ptr< SiStripFedZeroSuppression >	zs,
		std::unique_ptr< SiStripAPVRestorer >	res,
		bool	doAPVRest,
		bool	useCMMap
	)

private

Definition at line 14 of file SiStripRawProcessingAlgorithms.cc.

     : subtractorPed(std::move(ped)),
       subtractorCMN(std::move(cmn)),
       suppressor(std::move(zs)),
       restorer(std::move(res)),
       doAPVRestore(doAPVRest),
       useCMMeanMap(useCMMap) {}

Member Function Documentation

◆ convertVirginRawToHybrid() [1/2]

uint16_t SiStripRawProcessingAlgorithms::convertVirginRawToHybrid	(	const edm::DetSet< SiStripRawDigi > &	rawDigis,
		edm::DetSet< SiStripDigi > &	suppressedDigis
	)

Zero-suppress virgin raw data in "hybrid" mode

Subtracts pedestals (in 11bit mode, x->(x+1024-ped)/2) and common-mode noise, and inspects the digis then. If flagged by the hybrid APV inspector, the zero-suppression is performed as usual (evaluation and subtraction of the baseline, truncation). Otherwise, the pedestal-subtracted digis are saved in one 128-strip cluster.

Parameters

rawDigis	input (virgin) raw digis
output	zero-suppressed digis (or pedestal-subtracted digis, see above)

Returns: number of restored APVs

Definition at line 260 of file SiStripRawProcessingAlgorithms.cc.

                                                                                                            {
   digivector_t rawdigis;
   rawdigis.reserve(rawDigis.size());
   std::transform(std::begin(rawDigis), std::end(rawDigis), std::back_inserter(rawdigis), [](SiStripRawDigi digi) {
     return digi.adc();
   });
   return convertVirginRawToHybrid(rawDigis.id, 0, rawdigis, suppressedDigis);
 }

References SiStripRawDigi::adc(), begin, convertVirginRawToHybrid(), end, edm::DetSet< T >::id, edm::DetSet< T >::size(), and HcalDetIdTransform::transform().

◆ convertVirginRawToHybrid() [2/2]

uint16_t SiStripRawProcessingAlgorithms::convertVirginRawToHybrid	(	uint32_t	id,
		uint16_t	firstAPV,
		digivector_t &	procRawDigis,
		edm::DetSet< SiStripDigi > &	output
	)

Zero-suppress virgin raw data in "hybrid" mode

Subtracts pedestals (in 11bit mode, x->(x+1024-ped)/2) and common-mode noise, and inspects the digis then. If not flagged by the hybrid APV inspector, the zero-suppression is performed as usual (evaluation and subtraction of the baseline, truncation). Otherwise, the pedestal-subtracted digis (as above) are saved in one 128-strip cluster. Note: the APV restorer is used, it must be configured with APVInspectMode='HybridEmulation' if this method is called.

Parameters

id	module DetId
firstAPV	index of the first APV considered
procRawDigis	input (virgin raw) ADCs. Output: the ADCs after all subtractions, but before zero-suppression
output	zero-suppressed digis (or pedestal-subtracted digis, see above)

Returns: number of restored APVs

Definition at line 203 of file SiStripRawProcessingAlgorithms.cc.

                                                                                                   {
   digivector_t procRawDigisPedSubtracted;
  
   for (auto& digi : procRawDigis) {
     digi += 1024;
   }  // adding one MSB
  
   subtractorPed->subtract(id, firstAPV * 128, procRawDigis);  // all strips are pedestals subtracted
  
   for (auto& digi : procRawDigis) {
     digi /= 2;
   }
  
   procRawDigisPedSubtracted.assign(procRawDigis.begin(), procRawDigis.end());
  
   subtractorCMN->subtract(id, firstAPV, procRawDigis);
  
   const auto nAPVFlagged = restorer->inspect(id, firstAPV, procRawDigis, subtractorCMN->getAPVsCM());
  
   for (auto& digi : procRawDigis) {
     digi *= 2;
   }
  
   const std::vector<bool>& apvf = getAPVFlags();
   const std::size_t nAPVs = procRawDigis.size() / 128;
   for (uint16_t iAPV = firstAPV; iAPV < nAPVs + firstAPV; ++iAPV) {
     if (apvf[iAPV]) {
       //GB 23/6/08: truncation should be done at the very beginning
       for (uint16_t i = 0; i < 128; ++i) {
         const int16_t digi = procRawDigisPedSubtracted[128 * (iAPV - firstAPV) + i];
         output.push_back(SiStripDigi(128 * iAPV + i, (digi < 0 ? 0 : suppressor->truncate(digi))));
       }
     } else {
       const auto firstDigiIt = std::begin(procRawDigis) + 128 * (iAPV - firstAPV);
       std::vector<int16_t> singleAPVdigi(firstDigiIt, firstDigiIt + 128);
       suppressor->suppress(singleAPVdigi, iAPV, output);
     }
   }
  
   return nAPVFlagged;
 }

References begin, getAPVFlags(), mps_fire::i, convertSQLitetoXML_cfg::output, restorer, subtractorCMN, subtractorPed, and suppressor.

Referenced by convertVirginRawToHybrid().

◆ getAPVFlags()

const std::vector<bool>& SiStripRawProcessingAlgorithms::getAPVFlags ( ) const

inline

Definition at line 47 of file SiStripRawProcessingAlgorithms.h.

47 { return restorer->getAPVFlags(); }

References restorer.

Referenced by convertVirginRawToHybrid(), and suppressHybridData().

◆ getAPVsCM()

const SiStripAPVRestorer::medians_t& SiStripRawProcessingAlgorithms::getAPVsCM ( ) const

inline

Definition at line 52 of file SiStripRawProcessingAlgorithms.h.

52 { return subtractorCMN->getAPVsCM(); }

References subtractorCMN.

◆ getBaselineMap()

const SiStripAPVRestorer::baselinemap_t& SiStripRawProcessingAlgorithms::getBaselineMap ( ) const

inline

Definition at line 48 of file SiStripRawProcessingAlgorithms.h.

48 { return restorer->getBaselineMap(); }

References restorer.

◆ getSmoothedPoints()

const std::map<uint16_t, SiStripAPVRestorer::digimap_t>& SiStripRawProcessingAlgorithms::getSmoothedPoints ( ) const

inline

Definition at line 49 of file SiStripRawProcessingAlgorithms.h.

                                                                                         {
     return restorer->getSmoothedPoints();
   }

References restorer.

◆ initialize() [1/2]

void SiStripRawProcessingAlgorithms::initialize ( const edm::EventSetup & es )

Definition at line 27 of file SiStripRawProcessingAlgorithms.cc.

                                                                        {
   subtractorPed->init(es);
   subtractorCMN->init(es);
   suppressor->init(es);
   if (restorer.get())
     restorer->init(es);
  
   edm::ESHandle<TrackerGeometry> tracker;
   es.get<TrackerDigiGeometryRecord>().get(tracker);
   trGeo = tracker.product();
 }

References edm::EventSetup::get(), get, restorer, subtractorCMN, subtractorPed, suppressor, PbPb_ZMuSkimMuonDPG_cff::tracker, and trGeo.

Referenced by initialize().

◆ initialize() [2/2]

void SiStripRawProcessingAlgorithms::initialize	(	const edm::EventSetup &	es,
		const edm::Event &	e
	)

Definition at line 39 of file SiStripRawProcessingAlgorithms.cc.

                                                                                           {
   initialize(es);
   if (restorer.get() && doAPVRestore && useCMMeanMap)
     restorer->loadMeanCMMap(e);
 }

References doAPVRestore, MillePedeFileConverter_cfg::e, initialize(), restorer, and useCMMeanMap.

◆ suppressHybridData()

uint16_t SiStripRawProcessingAlgorithms::suppressHybridData	(	const edm::DetSet< SiStripDigi > &	hybridDigis,
		edm::DetSet< SiStripDigi > &	suppressedDigis,
		uint16_t	firstAPV = `0`
	)

Zero-suppress "hybrid" raw data

Subtracts common-mode noise, and inspects the digis then. If flagged by the APV inspector, the zero-suppression is performed as usual. Otherwise, the positive inputs are copied.

Parameters

hybridDigis	input ADCs in ZS format (regular ZS or "hybrid", i.e. processed as x->(x+1024-ped)/2)
suppressedDigis	zero-suppressed digis
firstAPV	(optional) number of the first APV for which digis are should be handled (otherwise all present)

Returns: number of restored APVs

Definition at line 58 of file SiStripRawProcessingAlgorithms.cc.

                                                                                {
   uint16_t nAPVFlagged = 0;
   auto beginAPV = hybridDigis.begin();
   const auto indigis_end = hybridDigis.end();
   auto iAPV = firstAPV;
   while (beginAPV != indigis_end) {
     const auto endAPV = std::lower_bound(beginAPV, indigis_end, SiStripDigi((iAPV + 1) * 128, 0));
     const auto nDigisInAPV = std::distance(beginAPV, endAPV);
     if (nDigisInAPV > 64) {
       digivector_t workDigis(128, -1024);
       for (auto it = beginAPV; it != endAPV; ++it) {
         workDigis[it->strip() - 128 * iAPV] = it->adc() * 2 - 1024;
       }
       digivector_t workDigisPedSubtracted(workDigis);
       subtractorCMN->subtract(hybridDigis.id, iAPV, workDigis);
       const auto apvFlagged = restorer->inspectAndRestore(
           hybridDigis.id, iAPV, workDigisPedSubtracted, workDigis, subtractorCMN->getAPVsCM());
       nAPVFlagged += apvFlagged;
       if (getAPVFlags()[iAPV]) {
         suppressor->suppress(workDigis, iAPV, suppressedDigis);
       } else {  // bad APV: more than 64 but not flagged
         for (uint16_t i = 0; i != 128; ++i) {
           const auto digi = workDigisPedSubtracted[i];
           if (digi > 0) {
             suppressedDigis.push_back(SiStripDigi(iAPV * 128 + i, suppressor->truncate(digi)));
           }
         }
       }
     } else {  // already zero-suppressed, copy and truncate
       std::transform(beginAPV, endAPV, std::back_inserter(suppressedDigis), [this](const SiStripDigi inDigi) {
         return SiStripDigi(inDigi.strip(), suppressor->truncate(inDigi.adc()));
       });
     }
     beginAPV = endAPV;
     ++iAPV;
   }
   return nAPVFlagged;
 }

References SiStripDigi::adc(), edm::DetSet< T >::begin(), HLT_2018_cff::distance, edm::DetSet< T >::end(), getAPVFlags(), mps_fire::i, edm::DetSet< T >::id, cuda_std::lower_bound(), edm::DetSet< T >::push_back(), restorer, SiStripDigi::strip(), subtractorCMN, suppressor, and HcalDetIdTransform::transform().

◆ suppressProcessedRawData() [1/2]

uint16_t SiStripRawProcessingAlgorithms::suppressProcessedRawData	(	const edm::DetSet< SiStripRawDigi > &	rawDigis,
		edm::DetSet< SiStripDigi > &	output
	)

Zero-suppress processed (pedestals-subtracted) raw data.

Subtracts common-mode noise and (optionally, if doAPVRestore) re-evaluates and subtracts the baseline.

Parameters

rawDigis	input (processed) raw digis
output	zero-suppressed digis

Returns: number of restored APVs

Definition at line 178 of file SiStripRawProcessingAlgorithms.cc.

                                                                                                   {
   digivector_t rawdigis;
   rawdigis.reserve(rawDigis.size());
   std::transform(std::begin(rawDigis), std::end(rawDigis), std::back_inserter(rawdigis), [](SiStripRawDigi digi) {
     return digi.adc();
   });
   return suppressProcessedRawData(rawDigis.id, 0, rawdigis, output);
 }

References SiStripRawDigi::adc(), begin, end, edm::DetSet< T >::id, convertSQLitetoXML_cfg::output, edm::DetSet< T >::size(), suppressProcessedRawData(), and HcalDetIdTransform::transform().

◆ suppressProcessedRawData() [2/2]

uint16_t SiStripRawProcessingAlgorithms::suppressProcessedRawData	(	uint32_t	id,
		uint16_t	firstAPV,
		digivector_t &	procRawDigis,
		edm::DetSet< SiStripDigi > &	output
	)

Zero-suppress processed (pedestals-subtracted) raw data.

Subtracts common-mode noise and (optionally, if doAPVRestore) re-evaluates and subtracts the baseline.

Parameters

id	module DetId
firstAPV	index of the first APV to consider
procRawDigis	input (processed raw) ADCs. Output: the ADCs after all subtractions, but before zero-suppression
output	zero-suppressed digis

Returns: number of restored APVs

Definition at line 151 of file SiStripRawProcessingAlgorithms.cc.

                                                                                                   {
   digivector_t procRawDigisPedSubtracted;
  
   int16_t nAPVFlagged = 0;
   if (doAPVRestore)
     procRawDigisPedSubtracted.assign(procRawDigis.begin(), procRawDigis.end());
   subtractorCMN->subtract(id, firstAPV, procRawDigis);
   if (doAPVRestore)
     nAPVFlagged =
         restorer->inspectAndRestore(id, firstAPV, procRawDigisPedSubtracted, procRawDigis, subtractorCMN->getAPVsCM());
   suppressor->suppress(procRawDigis, firstAPV, output);
   return nAPVFlagged;
 }

References doAPVRestore, convertSQLitetoXML_cfg::output, restorer, subtractorCMN, and suppressor.

Referenced by suppressProcessedRawData(), and suppressVirginRawData().

◆ suppressVirginRawData() [1/2]

uint16_t SiStripRawProcessingAlgorithms::suppressVirginRawData	(	const edm::DetSet< SiStripRawDigi > &	rawDigis,
		edm::DetSet< SiStripDigi > &	output
	)

Zero-suppress virgin raw data.

Subtracts pedestals and common-mode noise, and (optionally, if doAPVRestore) re-evaluates and subtracts the baseline.

Parameters

rawDigis	input (virgin) raw digis
output	zero-suppressed digis

Returns: number of restored APVs

Definition at line 129 of file SiStripRawProcessingAlgorithms.cc.

                                                                                                {
   digivector_t rawdigis;
   rawdigis.reserve(rawDigis.size());
   std::transform(std::begin(rawDigis), std::end(rawDigis), std::back_inserter(rawdigis), [](SiStripRawDigi digi) {
     return digi.adc();
   });
   return suppressVirginRawData(rawDigis.id, 0, rawdigis, output);
 }

References SiStripRawDigi::adc(), begin, end, edm::DetSet< T >::id, convertSQLitetoXML_cfg::output, edm::DetSet< T >::size(), suppressVirginRawData(), and HcalDetIdTransform::transform().

◆ suppressVirginRawData() [2/2]

uint16_t SiStripRawProcessingAlgorithms::suppressVirginRawData	(	uint32_t	id,
		uint16_t	firstAPV,
		digivector_t &	procRawDigis,
		edm::DetSet< SiStripDigi > &	output
	)

Zero-suppress virgin raw data.

Subtracts pedestals and common-mode noise, and (optionally, if doAPVRestore) re-evaluates and subtracts the baseline.

Parameters

id	module DetId
firstAPV	index of the first APV to consider
procRawDigis	input (virgin raw) ADCs. Output: the ADCs after all subtractions, but before zero-suppression
output	zero-suppressed digis

Returns: number of restored APVs

Definition at line 111 of file SiStripRawProcessingAlgorithms.cc.

                                                                                                {
   subtractorPed->subtract(id, firstAPV * 128, procRawDigis);
   return suppressProcessedRawData(id, firstAPV, procRawDigis, output);
 }