#include <HcalTriggerPrimitiveAlgo.h>

Classes
struct	HFDetails

struct	HFUpgradeDetails

Public Member Functions
template<typename T , typename... Args>
void	addDigis (const T &collection, const Args &...digis)

template<typename T >
void	addDigis (const T &collection)

template<typename D >
void	addDigis (const HcalDataFrameContainer< D > &collection)

	HcalTriggerPrimitiveAlgo (bool pf, const std::vector< double > &w, int latency, uint32_t FG_threshold, uint32_t FG_HF_threshold, uint32_t ZS_threshold, int numberOfSamples, int numberOfPresamples, int numberOfSamplesHF, int numberOfPresamplesHF, uint32_t minSignalThreshold=0, uint32_t PMT_NoiseThreshold=0)

void	overrideParameters (const edm::ParameterSet &ps)

template<typename... Digis>
void	run (const HcalTPGCoder incoder, const HcalTPGCompressor outcoder, const HcalDbService conditions, HcalTrigPrimDigiCollection &result, const HcalTrigTowerGeometry trigTowerGeometry, float rctlsb, const HcalFeatureBit *LongvrsShortCut, const Digis &...digis)

void	runFEFormatError (const FEDRawDataCollection rawraw, const HcalElectronicsMap emap, HcalTrigPrimDigiCollection &result)

void	runZS (HcalTrigPrimDigiCollection &tp)

void	setNCTScaleShift (int)

void	setPeakFinderAlgorithm (int algo)

void	setRCTScaleShift (int)

void	setUpgradeFlags (bool hb, bool he, bool hf)

	~HcalTriggerPrimitiveAlgo ()

Private Types
typedef std::map< HcalTrigTowerDetId, std::vector< bool > >	FGbitMap

typedef std::vector< HcalFinegrainBit::Tower >	FGUpgradeContainer

typedef std::map< HcalTrigTowerDetId, FGUpgradeContainer >	FGUpgradeMap

typedef std::map< HcalTrigTowerDetId, std::map< uint32_t, HFDetails > >	HFDetailMap

typedef std::map< HcalTrigTowerDetId, std::map< uint32_t, std::array< HFUpgradeDetails, 4 > > >	HFUpgradeDetailMap

typedef std::vector< IntegerCaloSamples >	SumFGContainer

typedef std::map< HcalTrigTowerDetId, IntegerCaloSamples >	SumMap

typedef std::map< HcalTrigTowerDetId, SumFGContainer >	TowerMapFGSum

typedef std::map< uint32_t, std::vector< bool > >	TowerMapVeto

Private Member Functions
void	addFG (const HcalTrigTowerDetId &id, std::vector< bool > &msb)

void	addSignal (const HBHEDataFrame &frame)
	adds the signal to the map More...

void	addSignal (const HFDataFrame &frame)

void	addSignal (const QIE10DataFrame &frame)

void	addSignal (const QIE11DataFrame &frame)

void	addSignal (const IntegerCaloSamples &samples)

void	addUpgradeFG (const HcalTrigTowerDetId &id, int depth, const std::vector< std::bitset< 2 >> &bits)

void	analyze (IntegerCaloSamples &samples, HcalTriggerPrimitiveDigi &result)
	adds the actual RecHits More...

void	analyze2017 (IntegerCaloSamples &samples, HcalTriggerPrimitiveDigi &result, const HcalFinegrainBit &fg_algo)

void	analyzeHF (IntegerCaloSamples &samples, HcalTriggerPrimitiveDigi &result, const int hf_lumi_shift)

void	analyzeHF2016 (const IntegerCaloSamples &SAMPLES, HcalTriggerPrimitiveDigi &result, const int HF_LUMI_SHIFT, const HcalFeatureBit *HCALFEM)

void	analyzeHF2017 (const IntegerCaloSamples &SAMPLES, HcalTriggerPrimitiveDigi &result, const int HF_LUMI_SHIFT, const HcalFeatureBit *HCALFEM)

bool	validChannel (const QIE10DataFrame &digi, int ts) const

bool	validUpgradeFG (const HcalTrigTowerDetId &id, int depth) const

Private Attributes
const HcalDbService *	conditions_

uint32_t	FG_HF_threshold_

uint32_t	FG_threshold_

FGbitMap	fgMap_

FGUpgradeMap	fgUpgradeMap_

TowerMapVeto	HF_Veto

const HcaluLUTTPGCoder *	incoder_

int	latency_

HcalFeatureBit *	LongvrsShortCut

uint32_t	minSignalThreshold_

int	NCTScaleShift

int	numberOfPresamples_

int	numberOfPresamplesHF_

int	numberOfSamples_

int	numberOfSamplesHF_

const HcalTPGCompressor *	outcoder_

bool	override_adc_hf_ = false

uint32_t	override_adc_hf_value_

edm::ParameterSet	override_parameters_

bool	override_tdc_hf_ = false

unsigned long long	override_tdc_hf_value_

int	peak_finder_algorithm_

bool	peakfind_

uint32_t	PMT_NoiseThreshold_

int	RCTScaleShift

HFDetailMap	theHFDetailMap

HFUpgradeDetailMap	theHFUpgradeDetailMap

SumMap	theSumMap

double	theThreshold

TowerMapFGSum	theTowerMapFGSum

const HcalTrigTowerGeometry *	theTrigTowerGeometry

bool	upgrade_hb_ = false

bool	upgrade_he_ = false

bool	upgrade_hf_ = false

std::vector< double >	weights_

uint32_t	ZS_threshold_

int	ZS_threshold_I_

Static Private Attributes
static const int	FIRST_FINEGRAIN_TOWER = 30

static const int	HBHE_OVERLAP_TOWER = 16

static const int	LAST_FINEGRAIN_DEPTH = 6

static const int	LAST_FINEGRAIN_TOWER = 28

static const int	QIE10_LINEARIZATION_ET = HcaluLUTTPGCoder::QIE10_LUT_BITMASK

static const int	QIE10_MAX_LINEARIZATION_ET = 0x7FF

static const int	QIE11_LINEARIZATION_ET = HcaluLUTTPGCoder::QIE11_LUT_BITMASK

static const int	QIE11_MAX_LINEARIZATION_ET = 0x7FF

static const int	QIE8_LINEARIZATION_ET = HcaluLUTTPGCoder::QIE8_LUT_BITMASK

Detailed Description

Definition at line 25 of file HcalTriggerPrimitiveAlgo.h.

Member Typedef Documentation

typedef std::map<HcalTrigTowerDetId, std::vector<bool> > HcalTriggerPrimitiveAlgo::FGbitMap

private

Definition at line 180 of file HcalTriggerPrimitiveAlgo.h.

typedef std::vector<HcalFinegrainBit::Tower> HcalTriggerPrimitiveAlgo::FGUpgradeContainer

private

Definition at line 183 of file HcalTriggerPrimitiveAlgo.h.

typedef std::map<HcalTrigTowerDetId, FGUpgradeContainer> HcalTriggerPrimitiveAlgo::FGUpgradeMap

private

Definition at line 184 of file HcalTriggerPrimitiveAlgo.h.

typedef std::map<HcalTrigTowerDetId, std::map<uint32_t, HFDetails> > HcalTriggerPrimitiveAlgo::HFDetailMap

private

Definition at line 151 of file HcalTriggerPrimitiveAlgo.h.

typedef std::map<HcalTrigTowerDetId, std::map<uint32_t, std::array<HFUpgradeDetails, 4> > > HcalTriggerPrimitiveAlgo::HFUpgradeDetailMap

private

Definition at line 160 of file HcalTriggerPrimitiveAlgo.h.

typedef std::vector<IntegerCaloSamples> HcalTriggerPrimitiveAlgo::SumFGContainer

private

Definition at line 163 of file HcalTriggerPrimitiveAlgo.h.

typedef std::map<HcalTrigTowerDetId, IntegerCaloSamples> HcalTriggerPrimitiveAlgo::SumMap

private

Definition at line 142 of file HcalTriggerPrimitiveAlgo.h.

typedef std::map< HcalTrigTowerDetId, SumFGContainer > HcalTriggerPrimitiveAlgo::TowerMapFGSum

private

Definition at line 164 of file HcalTriggerPrimitiveAlgo.h.

typedef std::map<uint32_t, std::vector<bool> > HcalTriggerPrimitiveAlgo::TowerMapVeto

private

Definition at line 177 of file HcalTriggerPrimitiveAlgo.h.

Constructor & Destructor Documentation

HcalTriggerPrimitiveAlgo::HcalTriggerPrimitiveAlgo	(	bool	pf,
		const std::vector< double > &	w,
		int	latency,
		uint32_t	FG_threshold,
		uint32_t	FG_HF_threshold,
		uint32_t	ZS_threshold,
		int	numberOfSamples,
		int	numberOfPresamples,
		int	numberOfSamplesHF,
		int	numberOfPresamplesHF,
		uint32_t	minSignalThreshold = `0`,
		uint32_t	PMT_NoiseThreshold = `0`
	)

Definition at line 23 of file HcalTriggerPrimitiveAlgo.cc.

References numberOfPresamples_, numberOfPresamplesHF_, numberOfSamples_, numberOfSamplesHF_, peakfind_, ZS_threshold_, and ZS_threshold_I_.

                                                    : incoder_(nullptr), outcoder_(nullptr),
                                                    theThreshold(0), peakfind_(pf), weights_(w), latency_(latency),
                                                    FG_threshold_(FG_threshold), FG_HF_threshold_(FG_HF_threshold), ZS_threshold_(ZS_threshold),
                                                    numberOfSamples_(numberOfSamples),
                                                    numberOfPresamples_(numberOfPresamples),
                                                    numberOfSamplesHF_(numberOfSamplesHF),
                                                    numberOfPresamplesHF_(numberOfPresamplesHF),
                                                    minSignalThreshold_(minSignalThreshold),
                                                    PMT_NoiseThreshold_(PMT_NoiseThreshold),
                                                    NCTScaleShift(0), RCTScaleShift(0),
                                                    peak_finder_algorithm_(2),
                                                    override_parameters_()
 {
    //No peak finding setting (for Fastsim)
    if (!peakfind_){
       numberOfSamples_ = 1; 
       numberOfPresamples_ = 0;
       numberOfSamplesHF_ = 1; 
       numberOfPresamplesHF_ = 0;
    }
    // Switch to integer for comparisons - remove compiler warning
    ZS_threshold_I_ = ZS_threshold_;
 }

HcalTriggerPrimitiveAlgo::~HcalTriggerPrimitiveAlgo ( )

Definition at line 54 of file HcalTriggerPrimitiveAlgo.cc.

54 {

55 }

Member Function Documentation

template<typename T , typename... Args>

void HcalTriggerPrimitiveAlgo::addDigis	(	const T &	collection,
		const Args &...	digis
	)

inline

Definition at line 44 of file HcalTriggerPrimitiveAlgo.h.

Referenced by run().

                                                            {
      addDigis(collection);
      addDigis(digis...);
   };

template<typename T >

void HcalTriggerPrimitiveAlgo::addDigis ( const T & collection )

inline

Definition at line 50 of file HcalTriggerPrimitiveAlgo.h.

References addSignal().

                                      {
      for (const auto& digi: collection) {
         addSignal(digi);
      }
   };

template<typename D >

void HcalTriggerPrimitiveAlgo::addDigis ( const HcalDataFrameContainer< D > & collection )

inline

Definition at line 57 of file HcalTriggerPrimitiveAlgo.h.

References addFG(), addSignal(), addUpgradeFG(), patPFMETCorrections_cff::algo, analyze(), analyze2017(), analyzeHF(), analyzeHF2016(), analyzeHF2017(), edm::DataFrameContainer::begin(), bits, particleFlowClusterECALTimeSelected_cfi::depth, edm::DataFrameContainer::end(), hcalSimParameters_cfi::hb, hcalSimParameters_cfi::he, photonIsolationHIProducer_cfi::hf, mps_fire::i, overrideParameters(), runFEFormatError(), runZS(), hcalTTPDigis_cfi::samples, setNCTScaleShift(), setPeakFinderAlgorithm(), setRCTScaleShift(), setUpgradeFlags(), validChannel(), and validUpgradeFG().

                                                              {
      for (auto i = collection.begin(); i != collection.end(); ++i) {
         D digi(*i);
         addSignal(digi);
      }
   };

void HcalTriggerPrimitiveAlgo::addFG	(	const HcalTrigTowerDetId &	id,
		std::vector< bool > &	msb
	)

private

Definition at line 727 of file HcalTriggerPrimitiveAlgo.cc.

References fgMap_, mps_fire::i, and triggerObjects_cff::id.

Referenced by addDigis(), and addSignal().

                                                                                       {
    FGbitMap::iterator itr = fgMap_.find(id);
    if (itr != fgMap_.end()){
       std::vector<bool>& _msb = itr->second;
       for (size_t i=0; i<msb.size(); ++i)
          _msb[i] = _msb[i] || msb[i];
    }
    else fgMap_[id] = msb;
 }

void HcalTriggerPrimitiveAlgo::addSignal ( const HBHEDataFrame & frame )

private

adds the signal to the map

Definition at line 83 of file HcalTriggerPrimitiveAlgo.cc.

References HcaluLUTTPGCoder::adc2Linear(), addFG(), HcalDetId::depth(), mps_fire::i, HBHEDataFrame::id(), electrons_cff::ids, incoder_, HcaluLUTTPGCoder::lookupMSB(), HBHEDataFrame::presamples(), HBHEDataFrame::size(), theTrigTowerGeometry, and HcalTrigTowerGeometry::towerIds().

Referenced by addDigis(), and addSignal().

                                                                     {
    // TODO: Need to add support for seperate 28, 29 in HE
    //Hack for 300_pre10, should be removed.
    if (frame.id().depth()==5) return;
 
    std::vector<HcalTrigTowerDetId> ids = theTrigTowerGeometry->towerIds(frame.id());
    assert(ids.size() == 1 || ids.size() == 2);
    IntegerCaloSamples samples1(ids[0], int(frame.size()));
 
    samples1.setPresamples(frame.presamples());
    incoder_->adc2Linear(frame, samples1);
 
    std::vector<bool> msb;
    incoder_->lookupMSB(frame, msb);
 
    if(ids.size() == 2) {
       // make a second trigprim for the other one, and split the energy
       IntegerCaloSamples samples2(ids[1], samples1.size());
       for(int i = 0; i < samples1.size(); ++i) {
          samples1[i] = uint32_t(samples1[i]*0.5);
          samples2[i] = samples1[i];
       }
       samples2.setPresamples(frame.presamples());
       addSignal(samples2);
       addFG(ids[1], msb);
    }
    addSignal(samples1);
    addFG(ids[0], msb);
 }

void HcalTriggerPrimitiveAlgo::addSignal ( const HFDataFrame & frame )

private

Definition at line 114 of file HcalTriggerPrimitiveAlgo.cc.

References HcaluLUTTPGCoder::adc2Linear(), addSignal(), HcalDetId::depth(), HF_Veto, mps_fire::i, HFDataFrame::id(), electrons_cff::ids, incoder_, HcalTriggerPrimitiveAlgo::HFDetails::long_fiber, HcalTriggerPrimitiveAlgo::HFDetails::LongDigi, HcalDetId::maskDepth(), minSignalThreshold_, HFDataFrame::presamples(), hcalTTPDigis_cfi::samples, IntegerCaloSamples::setPresamples(), HcalTriggerPrimitiveAlgo::HFDetails::short_fiber, HcalTriggerPrimitiveAlgo::HFDetails::ShortDigi, HFDataFrame::size(), theHFDetailMap, theTowerMapFGSum, theTrigTowerGeometry, HcalTrigTowerGeometry::towerIds(), and HcalTrigTowerDetId::version().

                                                                   {
    if(frame.id().depth() == 1 || frame.id().depth() == 2) {
       std::vector<HcalTrigTowerDetId> ids = theTrigTowerGeometry->towerIds(frame.id());
       std::vector<HcalTrigTowerDetId>::const_iterator it;
       for (it = ids.begin(); it != ids.end(); ++it) {
          HcalTrigTowerDetId trig_tower_id = *it;
          IntegerCaloSamples samples(trig_tower_id, frame.size());
          samples.setPresamples(frame.presamples());
          incoder_->adc2Linear(frame, samples);
 
          // Don't add to final collection yet
          // HF PMT veto sum is calculated in analyzerHF()
          IntegerCaloSamples zero_samples(trig_tower_id, frame.size());
          zero_samples.setPresamples(frame.presamples());
          addSignal(zero_samples);
 
          // Pre-LS1 Configuration
          if (trig_tower_id.version() == 0) {
             // Mask off depths: fgid is the same for both depths
             uint32_t fgid = (frame.id().maskDepth());
 
             if ( theTowerMapFGSum.find(trig_tower_id) == theTowerMapFGSum.end() ) {
                SumFGContainer sumFG;
                theTowerMapFGSum.insert(std::pair<HcalTrigTowerDetId, SumFGContainer>(trig_tower_id, sumFG));
             }
 
             SumFGContainer& sumFG = theTowerMapFGSum[trig_tower_id];
             SumFGContainer::iterator sumFGItr;
             for ( sumFGItr = sumFG.begin(); sumFGItr != sumFG.end(); ++sumFGItr) {
                if (sumFGItr->id() == fgid) { break; }
             }
             // If find
             if (sumFGItr != sumFG.end()) {
                for (int i=0; i<samples.size(); ++i) {
                   (*sumFGItr)[i] += samples[i];
                }
             }
             else {
                //Copy samples (change to fgid)
                IntegerCaloSamples sumFGSamples(DetId(fgid), samples.size());
                sumFGSamples.setPresamples(samples.presamples());
                for (int i=0; i<samples.size(); ++i) {
                   sumFGSamples[i] = samples[i];
                }
                sumFG.push_back(sumFGSamples);
             }
 
             // set veto to true if Long or Short less than threshold
             if (HF_Veto.find(fgid) == HF_Veto.end()) {
                vector<bool> vetoBits(samples.size(), false);
                HF_Veto[fgid] = vetoBits;
             }
             for (int i=0; i<samples.size(); ++i) {
                if (samples[i] < minSignalThreshold_) {
                   HF_Veto[fgid][i] = true;
                }
             }
          }
          // HF 1x1
          else if (trig_tower_id.version() == 1) {
             uint32_t fgid = (frame.id().maskDepth());
             HFDetails& details = theHFDetailMap[trig_tower_id][fgid];
             // Check the frame type to determine long vs short
             if (frame.id().depth() == 1) { // Long
                details.long_fiber = samples;
                details.LongDigi = frame;
             } else if (frame.id().depth() == 2) { // Short
                details.short_fiber = samples;
                details.ShortDigi = frame;
             } else {
                 // Neither long nor short... So we have no idea what to do
                 edm::LogWarning("HcalTPAlgo") << "Unable to figure out what to do with data frame for " << frame.id();
                 return;
             }
          }
          // Uh oh, we are in a bad/unknown state! Things will start crashing.
          else {
              return;
          }
       }
    }
 }

void HcalTriggerPrimitiveAlgo::addSignal ( const QIE10DataFrame & frame )

private

Definition at line 198 of file HcalTriggerPrimitiveAlgo.cc.

References HcaluLUTTPGCoder::adc2Linear(), addSignal(), QIE10DataFrame::detid(), HcalForward, triggerObjects_cff::id, QIE10DataFrame::id(), electrons_cff::ids, training_settings::idx, incoder_, HcaluLUTTPGCoder::lookupMSB(), QIE10DataFrame::presamples(), hcalTTPDigis_cfi::samples, QIE10DataFrame::samples(), IntegerCaloSamples::setPresamples(), HcalDetId::subdet(), theHFUpgradeDetailMap, theTrigTowerGeometry, HcalTrigTowerGeometry::towerIds(), validChannel(), and jets_cff::version.

 {
    HcalDetId detId = frame.detid();
    // prevent QIE10 calibration channels from entering TP emulation
    if(detId.subdet() != HcalForward) return;
 
    auto ids = theTrigTowerGeometry->towerIds(frame.id());
    for (const auto& id: ids) {
       if (id.version() == 0) {
          edm::LogError("HcalTPAlgo") << "Encountered QIE10 data frame mapped to TP version 0:" << id;
          continue;
       }
 
       int nsamples=frame.samples();
 
       IntegerCaloSamples samples(id, nsamples);
       samples.setPresamples(frame.presamples());
       incoder_->adc2Linear(frame, samples);
 
       // Don't add to final collection yet
       // HF PMT veto sum is calculated in analyzerHF()
       IntegerCaloSamples zero_samples(id, nsamples);
       zero_samples.setPresamples(frame.presamples());
       addSignal(zero_samples);
 
       auto fid = HcalDetId(frame.id());
       auto& details = theHFUpgradeDetailMap[id][fid.maskDepth()];
       auto& detail = details[fid.depth()-1];
       detail.samples = samples;
       detail.digi = frame;
       detail.validity.resize(nsamples);
       incoder_->lookupMSB(frame, detail.fgbit);
       for (int idx = 0; idx < nsamples; ++idx){
          detail.validity[idx] = validChannel(frame, idx);
       }
    }
 }

void HcalTriggerPrimitiveAlgo::addSignal ( const QIE11DataFrame & frame )

private

Definition at line 237 of file HcalTriggerPrimitiveAlgo.cc.

References HcaluLUTTPGCoder::adc2Linear(), addSignal(), addUpgradeFG(), HcalBarrel, HcalEndcap, mps_fire::i, QIE11DataFrame::id(), electrons_cff::ids, incoder_, HcaluLUTTPGCoder::lookupMSB(), QIE11DataFrame::presamples(), QIE11DataFrame::samples(), theTrigTowerGeometry, and HcalTrigTowerGeometry::towerIds().

 {
    HcalDetId detId(frame.id());
    // prevent QIE11 calibration channels from entering TP emulation
    if(detId.subdet() != HcalEndcap && detId.subdet() != HcalBarrel) return;
 
    std::vector<HcalTrigTowerDetId> ids = theTrigTowerGeometry->towerIds(detId);
    assert(ids.size() == 1 || ids.size() == 2);
    IntegerCaloSamples samples1(ids[0], int(frame.samples()));
 
    samples1.setPresamples(frame.presamples());
    incoder_->adc2Linear(frame, samples1);
 
    std::vector<std::bitset<2>> msb(frame.samples(), 0);
    incoder_->lookupMSB(frame, msb);
 
    if(ids.size() == 2) {
       // make a second trigprim for the other one, and share the energy
       IntegerCaloSamples samples2(ids[1], samples1.size());
       for(int i = 0; i < samples1.size(); ++i) {
          samples1[i] = uint32_t(samples1[i]);
          samples2[i] = samples1[i];
       }
       samples2.setPresamples(frame.presamples());
       addSignal(samples2);
       addUpgradeFG(ids[1], detId.depth(), msb);
    }
    addSignal(samples1);
    addUpgradeFG(ids[0], detId.depth(), msb);
 }

void HcalTriggerPrimitiveAlgo::addSignal ( const IntegerCaloSamples & samples )

private

Definition at line 268 of file HcalTriggerPrimitiveAlgo.cc.

References mps_fire::i, triggerObjects_cff::id, IntegerCaloSamples::id(), IntegerCaloSamples::size(), and theSumMap.

                                                                            {
    HcalTrigTowerDetId id(samples.id());
    SumMap::iterator itr = theSumMap.find(id);
    if(itr == theSumMap.end()) {
       theSumMap.insert(std::make_pair(id, samples));
    }
    else {
       // wish CaloSamples had a +=
       for(int i = 0; i < samples.size(); ++i) {
          (itr->second)[i] += samples[i];
       }
    }
 }

void HcalTriggerPrimitiveAlgo::addUpgradeFG	(	const HcalTrigTowerDetId &	id,
		int	depth,
		const std::vector< std::bitset< 2 >> &	bits
	)

private

Definition at line 748 of file HcalTriggerPrimitiveAlgo.cc.

References bits, particleFlowClusterECALTimeSelected_cfi::depth, fgUpgradeMap_, mps_fire::i, and validUpgradeFG().

Referenced by addDigis(), and addSignal().

 {
    if (not validUpgradeFG(id, depth)) {
       return;
    }
 
    auto it = fgUpgradeMap_.find(id);
    if (it == fgUpgradeMap_.end()) {
       FGUpgradeContainer element;
       element.resize(bits.size());
       it = fgUpgradeMap_.insert(std::make_pair(id, element)).first;
    }
    for (unsigned int i = 0; i < bits.size(); ++i) {
       it->second[i][0][depth] = bits[i][0];
       it->second[i][1][depth] = bits[i][1];
    }
 }

void HcalTriggerPrimitiveAlgo::analyze	(	IntegerCaloSamples &	samples,
		HcalTriggerPrimitiveDigi &	result
	)

private

adds the actual RecHits

Definition at line 283 of file HcalTriggerPrimitiveAlgo.cc.

References HcalTPGCompressor::compress(), fgMap_, mps_fire::i, IntegerCaloSamples::id(), training_settings::idx, createfilelist::int, numberOfPresamples_, numberOfSamples_, outcoder_, convertSQLitetoXML_cfg::output, peak_finder_algorithm_, peakfind_, IntegerCaloSamples::presamples(), QIE8_LINEARIZATION_ET, edm::shift, IntegerCaloSamples::size(), theThreshold, and weights_.

Referenced by addDigis(), and run().

                                                                                                       {
    int shrink = weights_.size() - 1;
    std::vector<bool>& msb = fgMap_[samples.id()];
    IntegerCaloSamples sum(samples.id(), samples.size());
 
    //slide algo window
    for(int ibin = 0; ibin < int(samples.size())- shrink; ++ibin) {
       int algosumvalue = 0;
       for(unsigned int i = 0; i < weights_.size(); i++) {
          //add up value * scale factor
          algosumvalue += int(samples[ibin+i] * weights_[i]);
       }
       if (algosumvalue<0) sum[ibin]=0;            // low-side
                                                   //high-side
       //else if (algosumvalue>QIE8_LINEARIZATION_ET) sum[ibin]=QIE8_LINEARIZATION_ET;
       else sum[ibin] = algosumvalue;              //assign value to sum[]
    }
 
    // Align digis and TP
    int dgPresamples=samples.presamples(); 
    int tpPresamples=numberOfPresamples_;
    int shift = dgPresamples - tpPresamples;
    int dgSamples=samples.size();
    int tpSamples=numberOfSamples_;
    if(peakfind_){
        if((shift<shrink) || (shift + tpSamples + shrink > dgSamples - (peak_finder_algorithm_ - 1) )   ){
         edm::LogInfo("HcalTriggerPrimitiveAlgo::analyze") << 
         "TP presample or size from the configuration file is out of the accessible range. Using digi values from data instead...";
         shift=shrink;
         tpPresamples=dgPresamples-shrink;
         tpSamples=dgSamples-(peak_finder_algorithm_-1)-shrink-shift;
        }
    }
 
    std::vector<int> finegrain(tpSamples,false);
 
    IntegerCaloSamples output(samples.id(), tpSamples);
    output.setPresamples(tpPresamples);
 
    for (int ibin = 0; ibin < tpSamples; ++ibin) {
       // ibin - index for output TP
       // idx - index for samples + shift
       int idx = ibin + shift;
 
       //Peak finding
       if (peakfind_) {
          bool isPeak = false;
          switch (peak_finder_algorithm_) {
             case 1 :
                isPeak = (samples[idx] > samples[idx-1] && samples[idx] >= samples[idx+1] && samples[idx] > theThreshold);
                break;
             case 2:
                isPeak = (sum[idx] > sum[idx-1] && sum[idx] >= sum[idx+1] && sum[idx] > theThreshold);
                break;
             default:
                break;
          }
 
          if (isPeak){
             output[ibin] = std::min<unsigned int>(sum[idx],QIE8_LINEARIZATION_ET);
             finegrain[ibin] = msb[idx];
          }
          // Not a peak
          else output[ibin] = 0;
       }
       else { // No peak finding, just output running sum
          output[ibin] = std::min<unsigned int>(sum[idx],QIE8_LINEARIZATION_ET);
          finegrain[ibin] = msb[idx];
       }
 
       // Only Pegged for 1-TS algo.
       if (peak_finder_algorithm_ == 1) {
          if (samples[idx] >= QIE8_LINEARIZATION_ET)
             output[ibin] = QIE8_LINEARIZATION_ET;
       }
    }
    outcoder_->compress(output, finegrain, result);
 }

void HcalTriggerPrimitiveAlgo::analyze2017	(	IntegerCaloSamples &	samples,
		HcalTriggerPrimitiveDigi &	result,
		const HcalFinegrainBit &	fg_algo
	)

private

Definition at line 364 of file HcalTriggerPrimitiveAlgo.cc.

References HcalTPGCompressor::compress(), HcalFinegrainBit::compute(), fgUpgradeMap_, mps_fire::i, IntegerCaloSamples::id(), electrons_cff::ids, training_settings::idx, createfilelist::int, numberOfPresamples_, numberOfSamples_, outcoder_, convertSQLitetoXML_cfg::output, peak_finder_algorithm_, IntegerCaloSamples::presamples(), QIE11_MAX_LINEARIZATION_ET, edm::shift, IntegerCaloSamples::size(), theThreshold, theTrigTowerGeometry, HcalTrigTowerGeometry::towerIds(), and weights_.

Referenced by addDigis(), and run().

 {
    int shrink = weights_.size() - 1;
    auto& msb = fgUpgradeMap_[samples.id()];
    IntegerCaloSamples sum(samples.id(), samples.size());
 
    HcalDetId detId(samples.id());
    std::vector<HcalTrigTowerDetId> ids = theTrigTowerGeometry->towerIds(detId);
    //slide algo window
    for(int ibin = 0; ibin < int(samples.size())- shrink; ++ibin) {
       int algosumvalue = 0;
       for(unsigned int i = 0; i < weights_.size(); i++) {
     //add up value * scale factor
     // In addition, divide by two in the 10 degree phi segmentation region
     // to mimic 5 degree segmentation for the trigger
     if(ids.size()==2) algosumvalue += int(samples[ibin+i] * 0.5 * weights_[i]);
     else algosumvalue += int(samples[ibin+i] * weights_[i]);
       }
       if (algosumvalue<0) sum[ibin]=0;            // low-side
                                                   //high-side
       //else if (algosumvalue>QIE11_LINEARIZATION_ET) sum[ibin]=QIE11_LINEARIZATION_ET;
       else sum[ibin] = algosumvalue;              //assign value to sum[]
    }
 
    // Align digis and TP
    int dgPresamples=samples.presamples(); 
    int tpPresamples=numberOfPresamples_;
    int shift = dgPresamples - tpPresamples;
    int dgSamples=samples.size();
    int tpSamples=numberOfSamples_;
 
    if((shift<shrink) || (shift + tpSamples + shrink > dgSamples - (peak_finder_algorithm_ - 1) )   ){
       edm::LogInfo("HcalTriggerPrimitiveAlgo::analyze") << 
          "TP presample or size from the configuration file is out of the accessible range. Using digi values from data instead...";
       shift=shrink;
       tpPresamples=dgPresamples-shrink;
       tpSamples=dgSamples-(peak_finder_algorithm_-1)-shrink-shift;
    }
 
    std::vector<int> finegrain(tpSamples,false);
 
    IntegerCaloSamples output(samples.id(), tpSamples);
    output.setPresamples(tpPresamples);
 
    for (int ibin = 0; ibin < tpSamples; ++ibin) {
       // ibin - index for output TP
       // idx - index for samples + shift
       int idx = ibin + shift;
       bool isPeak = (sum[idx] > sum[idx-1] && sum[idx] >= sum[idx+1] && sum[idx] > theThreshold);
 
       if (isPeak){
          output[ibin] = std::min<unsigned int>(sum[idx],QIE11_MAX_LINEARIZATION_ET);
          finegrain[ibin] = fg_algo.compute(msb[idx]).to_ulong();
       } else {
          // Not a peak
          output[ibin] = 0;
          finegrain[ibin] = 0;
       }
    }
    outcoder_->compress(output, finegrain, result);
 }

void HcalTriggerPrimitiveAlgo::analyzeHF	(	IntegerCaloSamples &	samples,
		HcalTriggerPrimitiveDigi &	result,
		const int	hf_lumi_shift
	)

private

Definition at line 427 of file HcalTriggerPrimitiveAlgo.cc.

References HcalTPGCompressor::compress(), FG_threshold_, HF_Veto, IntegerCaloSamples::id(), training_settings::idx, createfilelist::int, numberOfPresamplesHF_, numberOfSamplesHF_, outcoder_, convertSQLitetoXML_cfg::output, PMT_NoiseThreshold_, IntegerCaloSamples::presamples(), QIE8_LINEARIZATION_ET, edm::shift, IntegerCaloSamples::size(), theTowerMapFGSum, and TriggerAnalyzer::veto.

Referenced by addDigis(), and run().

                                                                                                                                  {
    HcalTrigTowerDetId detId(samples.id());
 
    // Align digis and TP
    int dgPresamples=samples.presamples(); 
    int tpPresamples=numberOfPresamplesHF_;
    int shift = dgPresamples - tpPresamples;
    int dgSamples=samples.size();
    int tpSamples=numberOfSamplesHF_;
    if(shift<0 || shift+tpSamples>dgSamples){
     edm::LogInfo("HcalTriggerPrimitiveAlgo::analyzeHF") << 
         "TP presample or size from the configuration file is out of the accessible range. Using digi values from data instead...";
     tpPresamples=dgPresamples;
     shift=0;
     tpSamples=dgSamples;
    }
 
    std::vector<int> finegrain(tpSamples, false);
 
    TowerMapFGSum::const_iterator tower2fg = theTowerMapFGSum.find(detId);
    assert(tower2fg != theTowerMapFGSum.end());
 
    const SumFGContainer& sumFG = tower2fg->second;
    // Loop over all L+S pairs that mapped from samples.id()
    // Note: 1 samples.id() = 6 x (L+S) without noZS
    for (SumFGContainer::const_iterator sumFGItr = sumFG.begin(); sumFGItr != sumFG.end(); ++sumFGItr) {
       const std::vector<bool>& veto = HF_Veto[sumFGItr->id().rawId()];
       for (int ibin = 0; ibin < tpSamples; ++ibin) {
          int idx = ibin + shift;
          // if not vetod, add L+S to total sum and calculate FG
      bool vetoed = idx<int(veto.size()) && veto[idx];
          if (!(vetoed && (*sumFGItr)[idx] > PMT_NoiseThreshold_)) {
             samples[idx] += (*sumFGItr)[idx];
             finegrain[ibin] = (finegrain[ibin] || (*sumFGItr)[idx] >= FG_threshold_);
          }
       }
    }
 
    IntegerCaloSamples output(samples.id(), tpSamples);
    output.setPresamples(tpPresamples);
 
    for (int ibin = 0; ibin < tpSamples; ++ibin) {
       int idx = ibin + shift;
       output[ibin] = samples[idx] >> hf_lumi_shift;
       static const int MAX_OUTPUT = QIE8_LINEARIZATION_ET;  // QIE8_LINEARIZATION_ET = 1023
       if (output[ibin] > MAX_OUTPUT) output[ibin] = MAX_OUTPUT;
    }
    outcoder_->compress(output, finegrain, result);
 }

void HcalTriggerPrimitiveAlgo::analyzeHF2016	(	const IntegerCaloSamples &	SAMPLES,
		HcalTriggerPrimitiveDigi &	result,
		const int	HF_LUMI_SHIFT,
		const HcalFeatureBit *	HCALFEM
	)

private

Definition at line 477 of file HcalTriggerPrimitiveAlgo.cc.

References stringResolutionProvider_cfi::bin, HcalTPGCompressor::compress(), FG_HF_threshold_, HcalFeatureBit::fineGrainbit(), FIRST_FINEGRAIN_TOWER, IntegerCaloSamples::id(), min(), numberOfPresamplesHF_, numberOfSamplesHF_, outcoder_, convertSQLitetoXML_cfg::output, IntegerCaloSamples::presamples(), QIE8_LINEARIZATION_ET, IntegerCaloSamples::size(), and theHFDetailMap.

Referenced by addDigis(), and run().

           {
     // Align digis and TP
     const int SHIFT = samples.presamples() - numberOfPresamplesHF_;
     assert(SHIFT >= 0);
     assert((SHIFT + numberOfSamplesHF_) <= samples.size());
 
     // Try to find the HFDetails from the map corresponding to our samples
     const HcalTrigTowerDetId detId(samples.id());
     HFDetailMap::const_iterator it = theHFDetailMap.find(detId);
     // Missing values will give an empty digi
     if (it == theHFDetailMap.end()) {
         return;
     }
 
     std::vector<std::bitset<2>> finegrain(numberOfSamplesHF_, false);
 
     // Set up out output of IntergerCaloSamples
     IntegerCaloSamples output(samples.id(), numberOfSamplesHF_);
     output.setPresamples(numberOfPresamplesHF_);
 
     for (const auto& item: it->second) {
         auto& details = item.second;
         for (int ibin = 0; ibin < numberOfSamplesHF_; ++ibin) {
             const int IDX = ibin + SHIFT;
             int long_fiber_val = 0;
             if (IDX < details.long_fiber.size()) {
                 long_fiber_val = details.long_fiber[IDX];
             }
             int short_fiber_val = 0;
             if (IDX < details.short_fiber.size()) {
                 short_fiber_val = details.short_fiber[IDX];
             }
             output[ibin] += (long_fiber_val + short_fiber_val);
 
             uint32_t ADCLong = details.LongDigi[ibin].adc();
             uint32_t ADCShort = details.ShortDigi[ibin].adc();
 
             if (details.LongDigi.id().ietaAbs() >= FIRST_FINEGRAIN_TOWER) {
                finegrain[ibin][1] = (ADCLong > FG_HF_threshold_ || ADCShort > FG_HF_threshold_);
 
                if (embit != nullptr)
                   finegrain[ibin][0] = embit->fineGrainbit(details.ShortDigi, details.LongDigi, ibin);
             }
         }
     }
 
     for (int bin = 0; bin < numberOfSamplesHF_; ++bin) {
        static const unsigned int MAX_OUTPUT = QIE8_LINEARIZATION_ET;  // QIE8_LINEARIZATION_ET = 1023
        output[bin] = min({MAX_OUTPUT, output[bin] >> hf_lumi_shift});
     }
 
     std::vector<int> finegrain_converted;
     for (const auto& fg: finegrain)
        finegrain_converted.push_back(fg.to_ulong());
     outcoder_->compress(output, finegrain_converted, result);
 }

void HcalTriggerPrimitiveAlgo::analyzeHF2017	(	const IntegerCaloSamples &	SAMPLES,
		HcalTriggerPrimitiveDigi &	result,
		const int	HF_LUMI_SHIFT,
		const HcalFeatureBit *	HCALFEM
	)

private

Definition at line 564 of file HcalTriggerPrimitiveAlgo.cc.

References stringResolutionProvider_cfi::bin, HcalTPGCompressor::compress(), HcalFeatureBit::fineGrainbit(), FIRST_FINEGRAIN_TOWER, mps_fire::i, IntegerCaloSamples::id(), training_settings::idx, createfilelist::int, min(), numberOfPresamplesHF_, numberOfSamplesHF_, or, outcoder_, convertSQLitetoXML_cfg::output, IntegerCaloSamples::presamples(), QIE10_LINEARIZATION_ET, QIE10_MAX_LINEARIZATION_ET, edm::shift, IntegerCaloSamples::size(), and theHFUpgradeDetailMap.

Referenced by addDigis(), and run().

 {
     // Align digis and TP
     const int shift = samples.presamples() - numberOfPresamplesHF_;
     assert(shift >= 0);
     assert((shift + numberOfSamplesHF_) <= samples.size());
     assert(hf_lumi_shift>=2);
 
     // Try to find the HFDetails from the map corresponding to our samples
     const HcalTrigTowerDetId detId(samples.id());
     auto it = theHFUpgradeDetailMap.find(detId);
     // Missing values will give an empty digi
     if (it == theHFUpgradeDetailMap.end()) {
         return;
     }
 
     std::vector<std::bitset<2>> finegrain(numberOfSamplesHF_, false);
 
     // Set up out output of IntergerCaloSamples
     IntegerCaloSamples output(samples.id(), numberOfSamplesHF_);
     output.setPresamples(numberOfPresamplesHF_);
 
     for (const auto& item: it->second) {
         auto& details = item.second;
         for (int ibin = 0; ibin < numberOfSamplesHF_; ++ibin) {
             const int idx = ibin + shift;
 
             int long_fiber_val = 0;
             int long_fiber_count = 0;
             int short_fiber_val = 0;
             int short_fiber_count = 0;
 
             bool saturated = false;
 
             for (auto i: {0, 2}) {
                if (idx < details[i].samples.size() and details[i].validity[idx]) {
                   long_fiber_val += details[i].samples[idx];
                   saturated = saturated || (details[i].samples[idx] == QIE10_LINEARIZATION_ET);
                   ++long_fiber_count;
                }
             }
             for (auto i: {1, 3}) {
                if (idx < details[i].samples.size() and details[i].validity[idx]) {
                   short_fiber_val += details[i].samples[idx];
                   saturated = saturated || (details[i].samples[idx] == QIE10_LINEARIZATION_ET);
                   ++short_fiber_count;
                }
             }
 
             if (saturated) {
                output[ibin] = QIE10_MAX_LINEARIZATION_ET;
             } else {
         // For details of the energy handling, see:
         // https://cms-docdb.cern.ch/cgi-bin/DocDB/ShowDocument?docid=12306
         // If both readouts are valid, average of the two energies is taken
         // division by 2 is compensated by adjusting the total scale shift in the end
         if (long_fiber_count == 2) long_fiber_val >>=1;
         if (short_fiber_count == 2) short_fiber_val >>=1;
         
         auto sum = long_fiber_val + short_fiber_val;
         // Similar to above, if both channels are valid, 
         // average of the two energies is calculated
         // division by 2 here is also compensated by adjusting the total scale shift in the end
         if (long_fiber_count > 0 and short_fiber_count > 0) sum >>=1;
 
         output[ibin] += sum;
             }
 
             for (const auto& detail: details) {
                if (idx < int(detail.digi.size()) and detail.validity[idx] and HcalDetId(detail.digi.id()).ietaAbs() >= FIRST_FINEGRAIN_TOWER) {
                   finegrain[ibin][1] = finegrain[ibin][1] or detail.fgbit[idx];
                }
             }
 
             if (embit != nullptr) {
                finegrain[ibin][0] = embit->fineGrainbit(
                      details[1].digi, details[3].digi,
                      details[0].digi, details[2].digi,
                      details[1].validity[idx], details[3].validity[idx],
                      details[0].validity[idx], details[2].validity[idx],
                      idx
                );
             }
         }
     }
 
     for (int bin = 0; bin < numberOfSamplesHF_; ++bin) {
        output[bin] = min({(unsigned int) QIE10_MAX_LINEARIZATION_ET, output[bin] >> (hf_lumi_shift-2)});
     }
     std::vector<int> finegrain_converted;
     for (const auto& fg: finegrain)
        finegrain_converted.push_back(fg.to_ulong());
     outcoder_->compress(output, finegrain_converted, result);
 }

void HcalTriggerPrimitiveAlgo::overrideParameters ( const edm::ParameterSet & ps )

Definition at line 68 of file HcalTriggerPrimitiveAlgo.cc.

References edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), override_adc_hf_, override_adc_hf_value_, override_parameters_, override_tdc_hf_, and override_tdc_hf_value_.

Referenced by addDigis(), and HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

 {
    override_parameters_ = ps;
 
    if (override_parameters_.exists("ADCThresholdHF")) {
       override_adc_hf_ = true;
       override_adc_hf_value_ = override_parameters_.getParameter<uint32_t>("ADCThresholdHF");
    }
    if (override_parameters_.exists("TDCMaskHF")) {
       override_tdc_hf_ = true;
       override_tdc_hf_value_ = override_parameters_.getParameter<unsigned long long>("TDCMaskHF");
    }
 }

template<typename... Digis>

void HcalTriggerPrimitiveAlgo::run	(	const HcalTPGCoder *	incoder,
		const HcalTPGCompressor *	outcoder,
		const HcalDbService *	conditions,
		HcalTrigPrimDigiCollection &	result,
		const HcalTrigTowerGeometry *	trigTowerGeometry,
		float	rctlsb,
		const HcalFeatureBit *	LongvrsShortCut,
		const Digis &...	digis
	)

Definition at line 215 of file HcalTriggerPrimitiveAlgo.h.

References addDigis(), analyze(), analyze2017(), analyzeHF(), analyzeHF2016(), analyzeHF2017(), edm::SortedCollection< T, SORT >::back(), conditions_, edm::ParameterSet::exists(), fgMap_, fgUpgradeMap_, HcalTrigTowerGeometry::firstHFTower(), HcalTPParameters::getFGVersionHBHE(), HcalDbService::getHcalTPParameters(), edm::ParameterSet::getParameter(), HF_Veto, incoder_, LongvrsShortCut, NCTScaleShift, or, outcoder_, override_parameters_, edm::SortedCollection< T, SORT >::push_back(), RCTScaleShift, theHFDetailMap, theHFUpgradeDetailMap, theSumMap, theTowerMapFGSum, theTrigTowerGeometry, upgrade_hb_, upgrade_he_, upgrade_hf_, and jets_cff::version.

Referenced by HcalTrigPrimDigiProducer::produce().

                                                           {
    theTrigTowerGeometry = trigTowerGeometry;
 
    incoder_ = dynamic_cast<const HcaluLUTTPGCoder*>(incoder);
    outcoder_ = outcoder;
    conditions_ = conditions;
 
    theSumMap.clear();
    theTowerMapFGSum.clear();
    HF_Veto.clear();
    fgMap_.clear();
    fgUpgradeMap_.clear();
    theHFDetailMap.clear();
    theHFUpgradeDetailMap.clear();
 
    // Add all digi collections
    addDigis(digis...);
 
    // Prepare the fine-grain calculation algorithm for HB/HE
    int version = 0;
    if (upgrade_he_ or upgrade_hb_)
       version = conditions_->getHcalTPParameters()->getFGVersionHBHE();
    if (override_parameters_.exists("FGVersionHBHE"))
       version = override_parameters_.getParameter<uint32_t>("FGVersionHBHE");
    HcalFinegrainBit fg_algo(version);
 
    // VME produces additional bits on the front used by lumi but not the
    // trigger, this shift corrects those out by right shifting over them.
    for (auto& item: theSumMap) {
       result.push_back(HcalTriggerPrimitiveDigi(item.first));
       HcalTrigTowerDetId detId(item.second.id());
       if(detId.ietaAbs() >= theTrigTowerGeometry->firstHFTower(detId.version())) { 
          if (detId.version() == 0) {
             analyzeHF(item.second, result.back(), RCTScaleShift);
          } else if (detId.version() == 1) {
             if (upgrade_hf_)
                analyzeHF2017(item.second, result.back(), NCTScaleShift, LongvrsShortCut);
             else
                analyzeHF2016(item.second, result.back(), NCTScaleShift, LongvrsShortCut);
          } else {
             // Things are going to go poorly
          }
       }
       else {
          // Determine which energy reconstruction path to take based on the
          // fine-grain availability:
          //  * QIE8 TP add entries into fgMap_
          //  * QIE11 TP add entries into fgUpgradeMap_
          //    (not for tower 16 unless HB is upgraded, too)
          if (fgMap_.find(item.first) != fgMap_.end()) {
             analyze(item.second, result.back());
          } else if (fgUpgradeMap_.find(item.first) != fgUpgradeMap_.end()) {
             analyze2017(item.second, result.back(), fg_algo);
          }
       }
    }
 
    // Free up some memory
    theSumMap.clear();
    theTowerMapFGSum.clear();
    HF_Veto.clear();
    fgMap_.clear();
    fgUpgradeMap_.clear();
    theHFDetailMap.clear();
    theHFUpgradeDetailMap.clear();
 
    return;
 }

void HcalTriggerPrimitiveAlgo::runFEFormatError	(	const FEDRawDataCollection *	rawraw,
		const HcalElectronicsMap *	emap,
		HcalTrigPrimDigiCollection &	result
	)

Definition at line 675 of file HcalTriggerPrimitiveAlgo.cc.

References edm::SortedCollection< T, SORT >::begin(), HcalHTRData::check(), FEDRawData::data(), DetId::det(), runTauDisplay::eid, edm::SortedCollection< T, SORT >::end(), FEDRawDataCollection::FEDData(), HcalHTRData::getErrorsWord(), HcalDCCHeader::getSourceId(), HcalDCCHeader::getSpigotData(), HcalDCCHeader::getSpigotPresent(), HcalBarrel, HcalEndcap, HcalForward, HcalHTRData::htrSlot(), HcalHTRData::htrTopBottom(), mps_fire::i, electrons_cff::ids, HcalHTRData::isHistogramEvent(), HcalElectronicsMap::lookup(), FEDNumbering::MAXHCALFEDID, FEDNumbering::MINHCALFEDID, DetId::null(), HcalHTRData::readoutVMECrateId(), HcalElectronicsId::setHTR(), FEDRawData::size(), HcalDCCHeader::SPIGOT_COUNT, DetId::subdetId(), theTrigTowerGeometry, and HcalTrigTowerGeometry::towerIds().

Referenced by addDigis(), and HcalTrigPrimDigiProducer::produce().

                                                  {
   std::set<uint32_t> FrontEndErrors;
 
   for(int i=FEDNumbering::MINHCALFEDID; i<=FEDNumbering::MAXHCALFEDID; ++i) {
     const FEDRawData& raw = rawraw->FEDData(i);
     if (raw.size()<12) continue;
     const HcalDCCHeader* dccHeader=(const HcalDCCHeader*)(raw.data());
     if(!dccHeader) continue;
     HcalHTRData htr;
     for (int spigot=0; spigot<HcalDCCHeader::SPIGOT_COUNT; spigot++) {
       if (!dccHeader->getSpigotPresent(spigot)) continue;
       dccHeader->getSpigotData(spigot,htr,raw.size());
       int dccid = dccHeader->getSourceId();
       int errWord = htr.getErrorsWord() & 0x1FFFF;
       bool HTRError = (!htr.check() || htr.isHistogramEvent() || (errWord & 0x800)!=0);
 
       if(HTRError) {
         bool valid =false;
         for(int fchan=0; fchan<3 && !valid; fchan++) {
           for(int fib=0; fib<9 && !valid; fib++) {
             HcalElectronicsId eid(fchan,fib,spigot,dccid-FEDNumbering::MINHCALFEDID);
             eid.setHTR(htr.readoutVMECrateId(),htr.htrSlot(),htr.htrTopBottom());
             DetId detId = emap->lookup(eid);
             if(detId.null()) continue;
             HcalSubdetector subdet=(HcalSubdetector(detId.subdetId()));
             if (detId.det()!=4||
               (subdet!=HcalBarrel && subdet!=HcalEndcap &&
               subdet!=HcalForward )) continue;
             std::vector<HcalTrigTowerDetId> ids = theTrigTowerGeometry->towerIds(detId);
             for (std::vector<HcalTrigTowerDetId>::const_iterator triggerId=ids.begin(); triggerId != ids.end(); ++triggerId) {
               FrontEndErrors.insert(triggerId->rawId());
             }
             //valid = true;
           }
         }
       }
     }
   }
 
   // Loop over TP collection
   // Set TP to zero if there is FE Format Error
   HcalTriggerPrimitiveSample zeroSample(0);
   for (HcalTrigPrimDigiCollection::iterator tp = result.begin(); tp != result.end(); ++tp){
     if (FrontEndErrors.find(tp->id().rawId()) != FrontEndErrors.end()) {
       for (int i=0; i<tp->size(); ++i) tp->setSample(i, zeroSample);
     }
   }
 }

void HcalTriggerPrimitiveAlgo::runZS ( HcalTrigPrimDigiCollection & tp )

Definition at line 661 of file HcalTriggerPrimitiveAlgo.cc.

References edm::SortedCollection< T, SORT >::begin(), edm::SortedCollection< T, SORT >::end(), mps_fire::i, CaloTPGTranscoder_cfi::ZS, and ZS_threshold_I_.

Referenced by addDigis(), and HcalTrigPrimDigiProducer::produce().

                                                                        {
    for (HcalTrigPrimDigiCollection::iterator tp = result.begin(); tp != result.end(); ++tp){
       bool ZS = true;
       for (int i=0; i<tp->size(); ++i) {
          if (tp->sample(i).compressedEt()  > ZS_threshold_I_) {
             ZS=false;
             break;
          }
       }
       if (ZS) tp->setZSInfo(false,true);
       else tp->setZSInfo(true,false);
    }
 }

void HcalTriggerPrimitiveAlgo::setNCTScaleShift ( int shift )

Definition at line 772 of file HcalTriggerPrimitiveAlgo.cc.

References NCTScaleShift, and edm::shift.

Referenced by addDigis(), and HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

                                                         {
    NCTScaleShift = shift;
 }

void HcalTriggerPrimitiveAlgo::setPeakFinderAlgorithm ( int algo )

Definition at line 766 of file HcalTriggerPrimitiveAlgo.cc.

References patPFMETCorrections_cff::algo, Exception, and peak_finder_algorithm_.

Referenced by addDigis(), and HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

                                                              {
    if (algo <=0 && algo>2)
       throw cms::Exception("ERROR: Only algo 1 & 2 are supported.") << std::endl;
    peak_finder_algorithm_ = algo;
 }

void HcalTriggerPrimitiveAlgo::setRCTScaleShift ( int shift )

Definition at line 776 of file HcalTriggerPrimitiveAlgo.cc.

References RCTScaleShift, and edm::shift.

Referenced by addDigis(), and HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

                                                         {
    RCTScaleShift = shift;
 }

void HcalTriggerPrimitiveAlgo::setUpgradeFlags	(	bool	hb,
		bool	he,
		bool	hf
	)

Definition at line 59 of file HcalTriggerPrimitiveAlgo.cc.

References hcalSimParameters_cfi::hb, hcalSimParameters_cfi::he, photonIsolationHIProducer_cfi::hf, upgrade_hb_, upgrade_he_, and upgrade_hf_.

Referenced by addDigis(), and HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

 {
    upgrade_hb_ = hb;
    upgrade_he_ = he;
    upgrade_hf_ = hf;
 }

bool HcalTriggerPrimitiveAlgo::validChannel	(	const QIE10DataFrame &	digi,
		int	ts
	)		const

private

Definition at line 540 of file HcalTriggerPrimitiveAlgo.cc.

References ecalMGPA::adc(), conditions_, HcalDbService::getHcalTPChannelParameter(), HcalDbService::getHcalTPParameters(), HcalTPChannelParameter::getMask(), QIE10DataFrame::id(), QIE10DataFrame::linkError(), RecoTauDiscriminantConfiguration::mask, override_adc_hf_, override_adc_hf_value_, override_tdc_hf_, override_tdc_hf_value_, metProducer_cfi::parameters, and QIE10DataFrame::samples().

Referenced by addDigis(), and addSignal().

 {
    // channels with invalid data should not contribute to the sum
    if(digi.linkError() || ts>=digi.samples() || !digi[ts].ok()) return false;
 
    auto mask = conditions_->getHcalTPChannelParameter(HcalDetId(digi.id()))->getMask();
    if (mask)
       return false;
 
    auto parameters = conditions_->getHcalTPParameters();
    auto adc_threshold = parameters->getADCThresholdHF();
    auto tdc_mask = parameters->getTDCMaskHF();
 
    if (override_adc_hf_)
       adc_threshold = override_adc_hf_value_;
    if (override_tdc_hf_)
       tdc_mask = override_tdc_hf_value_;
 
    if (digi[ts].adc() < adc_threshold)
       return true;
 
    return (1ul << digi[ts].le_tdc()) & tdc_mask;
 }

bool HcalTriggerPrimitiveAlgo::validUpgradeFG	(	const HcalTrigTowerDetId &	id,
		int	depth
	)		const

private

Definition at line 738 of file HcalTriggerPrimitiveAlgo.cc.

References LAST_FINEGRAIN_DEPTH, and LAST_FINEGRAIN_TOWER.

Referenced by addDigis(), and addUpgradeFG().

 {
    if (depth > LAST_FINEGRAIN_DEPTH)
       return false;
    if (id.ietaAbs() > LAST_FINEGRAIN_TOWER)
       return false;
    return true;
 }

Member Data Documentation

const HcalDbService* HcalTriggerPrimitiveAlgo::conditions_

private

Definition at line 113 of file HcalTriggerPrimitiveAlgo.h.

Referenced by run(), and validChannel().

uint32_t HcalTriggerPrimitiveAlgo::FG_HF_threshold_

private

Definition at line 119 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHF2016().

uint32_t HcalTriggerPrimitiveAlgo::FG_threshold_

private

Definition at line 118 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHF().

FGbitMap HcalTriggerPrimitiveAlgo::fgMap_

private

Definition at line 181 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addFG(), analyze(), and run().

FGUpgradeMap HcalTriggerPrimitiveAlgo::fgUpgradeMap_

private

Definition at line 185 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addUpgradeFG(), analyze2017(), and run().

const int HcalTriggerPrimitiveAlgo::FIRST_FINEGRAIN_TOWER = 30

staticprivate

Definition at line 204 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHF2016(), and analyzeHF2017().

const int HcalTriggerPrimitiveAlgo::HBHE_OVERLAP_TOWER = 16

staticprivate

Definition at line 199 of file HcalTriggerPrimitiveAlgo.h.

TowerMapVeto HcalTriggerPrimitiveAlgo::HF_Veto

private

Definition at line 178 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal(), analyzeHF(), and run().

const HcaluLUTTPGCoder* HcalTriggerPrimitiveAlgo::incoder_

private

Definition at line 111 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal(), and run().

const int HcalTriggerPrimitiveAlgo::LAST_FINEGRAIN_DEPTH = 6

staticprivate

Definition at line 200 of file HcalTriggerPrimitiveAlgo.h.

Referenced by validUpgradeFG().

const int HcalTriggerPrimitiveAlgo::LAST_FINEGRAIN_TOWER = 28

staticprivate

Definition at line 201 of file HcalTriggerPrimitiveAlgo.h.

Referenced by validUpgradeFG().

int HcalTriggerPrimitiveAlgo::latency_

private

Definition at line 117 of file HcalTriggerPrimitiveAlgo.h.

HcalFeatureBit* HcalTriggerPrimitiveAlgo::LongvrsShortCut

private

Definition at line 176 of file HcalTriggerPrimitiveAlgo.h.

Referenced by run().

uint32_t HcalTriggerPrimitiveAlgo::minSignalThreshold_

private

Definition at line 126 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal().

int HcalTriggerPrimitiveAlgo::NCTScaleShift

private

Definition at line 128 of file HcalTriggerPrimitiveAlgo.h.

Referenced by run(), and setNCTScaleShift().

int HcalTriggerPrimitiveAlgo::numberOfPresamples_

private

Definition at line 123 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyze(), analyze2017(), and HcalTriggerPrimitiveAlgo().

int HcalTriggerPrimitiveAlgo::numberOfPresamplesHF_

private

Definition at line 125 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHF(), analyzeHF2016(), analyzeHF2017(), and HcalTriggerPrimitiveAlgo().

int HcalTriggerPrimitiveAlgo::numberOfSamples_

private

Definition at line 122 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyze(), analyze2017(), and HcalTriggerPrimitiveAlgo().

int HcalTriggerPrimitiveAlgo::numberOfSamplesHF_

private

Definition at line 124 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHF(), analyzeHF2016(), analyzeHF2017(), and HcalTriggerPrimitiveAlgo().

const HcalTPGCompressor* HcalTriggerPrimitiveAlgo::outcoder_

private

Definition at line 112 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyze(), analyze2017(), analyzeHF(), analyzeHF2016(), analyzeHF2017(), and run().

bool HcalTriggerPrimitiveAlgo::override_adc_hf_ = false

private

Definition at line 193 of file HcalTriggerPrimitiveAlgo.h.

Referenced by overrideParameters(), and validChannel().

uint32_t HcalTriggerPrimitiveAlgo::override_adc_hf_value_

private

Definition at line 194 of file HcalTriggerPrimitiveAlgo.h.

Referenced by overrideParameters(), and validChannel().

edm::ParameterSet HcalTriggerPrimitiveAlgo::override_parameters_

private

Definition at line 191 of file HcalTriggerPrimitiveAlgo.h.

Referenced by overrideParameters(), and run().

bool HcalTriggerPrimitiveAlgo::override_tdc_hf_ = false

private

Definition at line 195 of file HcalTriggerPrimitiveAlgo.h.

Referenced by overrideParameters(), and validChannel().

unsigned long long HcalTriggerPrimitiveAlgo::override_tdc_hf_value_

private

Definition at line 196 of file HcalTriggerPrimitiveAlgo.h.

Referenced by overrideParameters(), and validChannel().

int HcalTriggerPrimitiveAlgo::peak_finder_algorithm_

private

Definition at line 135 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyze(), analyze2017(), and setPeakFinderAlgorithm().

bool HcalTriggerPrimitiveAlgo::peakfind_

private

Definition at line 115 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyze(), and HcalTriggerPrimitiveAlgo().

uint32_t HcalTriggerPrimitiveAlgo::PMT_NoiseThreshold_

private

Definition at line 127 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHF().

const int HcalTriggerPrimitiveAlgo::QIE10_LINEARIZATION_ET = HcaluLUTTPGCoder::QIE10_LUT_BITMASK

staticprivate

Definition at line 207 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHF2017().

const int HcalTriggerPrimitiveAlgo::QIE10_MAX_LINEARIZATION_ET = 0x7FF

staticprivate

Definition at line 210 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHF2017().

const int HcalTriggerPrimitiveAlgo::QIE11_LINEARIZATION_ET = HcaluLUTTPGCoder::QIE11_LUT_BITMASK

staticprivate

Definition at line 208 of file HcalTriggerPrimitiveAlgo.h.

const int HcalTriggerPrimitiveAlgo::QIE11_MAX_LINEARIZATION_ET = 0x7FF

staticprivate

Definition at line 211 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyze2017().

const int HcalTriggerPrimitiveAlgo::QIE8_LINEARIZATION_ET = HcaluLUTTPGCoder::QIE8_LUT_BITMASK

staticprivate

Definition at line 206 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyze(), analyzeHF(), and analyzeHF2016().

int HcalTriggerPrimitiveAlgo::RCTScaleShift

private

Definition at line 129 of file HcalTriggerPrimitiveAlgo.h.

Referenced by run(), and setRCTScaleShift().

HFDetailMap HcalTriggerPrimitiveAlgo::theHFDetailMap

private

Definition at line 152 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal(), analyzeHF2016(), and run().

HFUpgradeDetailMap HcalTriggerPrimitiveAlgo::theHFUpgradeDetailMap

private

Definition at line 161 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal(), analyzeHF2017(), and run().

SumMap HcalTriggerPrimitiveAlgo::theSumMap

private

Definition at line 143 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal(), and run().

double HcalTriggerPrimitiveAlgo::theThreshold

private

Definition at line 114 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyze(), and analyze2017().

TowerMapFGSum HcalTriggerPrimitiveAlgo::theTowerMapFGSum

private

Definition at line 165 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal(), analyzeHF(), and run().

const HcalTrigTowerGeometry* HcalTriggerPrimitiveAlgo::theTrigTowerGeometry

private

Definition at line 140 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal(), analyze2017(), run(), and runFEFormatError().

bool HcalTriggerPrimitiveAlgo::upgrade_hb_ = false

private

Definition at line 187 of file HcalTriggerPrimitiveAlgo.h.

Referenced by run(), and setUpgradeFlags().

bool HcalTriggerPrimitiveAlgo::upgrade_he_ = false

private

Definition at line 188 of file HcalTriggerPrimitiveAlgo.h.

Referenced by run(), and setUpgradeFlags().

bool HcalTriggerPrimitiveAlgo::upgrade_hf_ = false

private

Definition at line 189 of file HcalTriggerPrimitiveAlgo.h.

Referenced by run(), and setUpgradeFlags().

std::vector<double> HcalTriggerPrimitiveAlgo::weights_

private

Definition at line 116 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyze(), and analyze2017().

uint32_t HcalTriggerPrimitiveAlgo::ZS_threshold_

private

Definition at line 120 of file HcalTriggerPrimitiveAlgo.h.

Referenced by HcalTriggerPrimitiveAlgo().

int HcalTriggerPrimitiveAlgo::ZS_threshold_I_

private

Definition at line 121 of file HcalTriggerPrimitiveAlgo.h.

Referenced by HcalTriggerPrimitiveAlgo(), and runZS().

Classes

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation