#include <HcalTriggerPrimitiveAlgo.h>

Public Member Functions
	HcalTriggerPrimitiveAlgo (bool pf, const std::vector< double > &w, int latency, uint32_t FG_threshold1, uint32_t FG_threshold2, uint32_t ZS_threshold, int numberOfSamples, int numberOfPresamples, int numberOfSamplesHF, int numberOfPresamplesHF, uint32_t minSignalThreshold=0, uint32_t PMT_NoiseThreshold=0)

void	run (const HcalTPGCoder incoder, const HcalTPGCompressor outcoder, const HBHEDigiCollection &hbheDigis, const HFDigiCollection &hfDigis, HcalTrigPrimDigiCollection &result, const HcalTrigTowerGeometry *trigTowerGeometry, float rctlsb)

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

void	runZS (HcalTrigPrimDigiCollection &tp)

void	setPeakFinderAlgorithm (int algo)

	~HcalTriggerPrimitiveAlgo ()

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

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 IntegerCaloSamples &samples)

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

void	analyzeHF (int ieta, int iphi, int(&FGBit)[18][22], IntegerCaloSamples &samples, HcalTriggerPrimitiveDigi &result, float rctlsb)

Private Attributes
uint32_t	FG_threshold1_

uint32_t	FG_threshold2_

FGbitMap	fgMap_

TowerMapVeto	HF_Veto

const HcaluLUTTPGCoder *	incoder_

int	latency_

uint32_t	minSignalThreshold_

int	numberOfPresamples_

int	numberOfPresamplesHF_

int	numberOfSamples_

int	numberOfSamplesHF_

const HcalTPGCompressor *	outcoder_

int	peak_finder_algorithm_

bool	peakfind_

uint32_t	PMT_NoiseThreshold_

SumMap	theSumMap

double	theThreshold

TowerMapFGSum	theTowerMapFGSum

const HcalTrigTowerGeometry *	theTrigTowerGeometry

std::vector< double >	weights_

uint32_t	ZS_threshold_

int	ZS_threshold_I_

Detailed Description

Definition at line 20 of file HcalTriggerPrimitiveAlgo.h.

Member Typedef Documentation

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

private

Definition at line 103 of file HcalTriggerPrimitiveAlgo.h.

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

private

Definition at line 87 of file HcalTriggerPrimitiveAlgo.h.

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

private

Definition at line 84 of file HcalTriggerPrimitiveAlgo.h.

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

private

Definition at line 88 of file HcalTriggerPrimitiveAlgo.h.

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

private

Definition at line 100 of file HcalTriggerPrimitiveAlgo.h.

Constructor & Destructor Documentation

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

Definition at line 15 of file HcalTriggerPrimitiveAlgo.cc.

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

                                                    : incoder_(0), outcoder_(0),
                                                    theThreshold(0), peakfind_(pf), weights_(w), latency_(latency),
                                                    FG_threshold1_(FG_threshold1), FG_threshold2_(FG_threshold2), ZS_threshold_(ZS_threshold),
                                                    numberOfSamples_(numberOfSamples),
                                                    numberOfPresamples_(numberOfPresamples),
                                                    numberOfSamplesHF_(numberOfSamplesHF),
                                                    numberOfPresamplesHF_(numberOfPresamplesHF),
                                                    minSignalThreshold_(minSignalThreshold),
                                                    PMT_NoiseThreshold_(PMT_NoiseThreshold),
                                                    peak_finder_algorithm_(2)
 {
    //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 43 of file HcalTriggerPrimitiveAlgo.cc.

43 {

44 }

Member Function Documentation

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

private

Definition at line 817 of file HcalTriggerPrimitiveAlgo.cc.

References fgMap_, and i.

Referenced by 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 518 of file HcalTriggerPrimitiveAlgo.cc.

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

Referenced by addSignal(), and run().

                                                                     {
    //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 548 of file HcalTriggerPrimitiveAlgo.cc.

References HcaluLUTTPGCoder::adc2Linear(), addSignal(), assert(), HcalDetId::depth(), HF_Veto, i, HFDataFrame::id(), incoder_, HcalDetId::maskDepth(), minSignalThreshold_, HFDataFrame::presamples(), IntegerCaloSamples::setPresamples(), HFDataFrame::size(), theTowerMapFGSum, theTrigTowerGeometry, and HcalTrigTowerGeometry::towerIds().

                                                                   {
 
    if(frame.id().depth() == 1 || frame.id().depth() == 2) {
       std::vector<HcalTrigTowerDetId> ids = theTrigTowerGeometry->towerIds(frame.id());
       assert(ids.size() == 1);
       IntegerCaloSamples samples(ids[0], 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(ids[0], frame.size());
       zero_samples.setPresamples(frame.presamples());
       addSignal(zero_samples);
 
       // Mask off depths: fgid is the same for both depths
       uint32_t fgid = (frame.id().maskDepth());
 
       if ( theTowerMapFGSum.find(ids[0]) == theTowerMapFGSum.end() ) {
          SumFGContainer sumFG;
          theTowerMapFGSum.insert(std::pair<HcalTrigTowerDetId, SumFGContainer >(ids[0], sumFG));
       }
 
       SumFGContainer& sumFG = theTowerMapFGSum[ids[0]];
       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;
    }
 }

void HcalTriggerPrimitiveAlgo::addSignal ( const IntegerCaloSamples & samples )

private

Definition at line 600 of file HcalTriggerPrimitiveAlgo.cc.

References i, 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::analyze	(	IntegerCaloSamples &	samples,
		HcalTriggerPrimitiveDigi &	result
	)

private

adds the actual RecHits

Definition at line 615 of file HcalTriggerPrimitiveAlgo.cc.

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

Referenced by 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>0x3FF) sum[ibin]=0x3FF;
       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<bool> 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],0x3FF);
             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],0x3FF);
          finegrain[ibin] = msb[idx];
       }
 
       // Only Pegged for 1-TS algo.
       if (peak_finder_algorithm_ == 1) {
          if (samples[idx] >= 0x3FF)
             output[ibin] = 0x3FF;
       }
    }
    outcoder_->compress(output, finegrain, result);
 }

void HcalTriggerPrimitiveAlgo::analyzeHF	(	int	ieta,
		int	iphi,
		int(&)	FGBit[18][22],
		IntegerCaloSamples &	samples,
		HcalTriggerPrimitiveDigi &	result,
		float	rctlsb
	)

private

Definition at line 695 of file HcalTriggerPrimitiveAlgo.cc.

References assert(), HcalTPGCompressor::compress(), converttoGCTetaFromTPG(), converttoGCTphiFromTPG(), HF_Veto, IntegerCaloSamples::id(), customizeTrackingMonitorSeedNumber::idx, cmsHarvester::index, numberOfPresamplesHF_, numberOfSamplesHF_, outcoder_, convertSQLitetoXML_cfg::output, PMT_NoiseThreshold_, IntegerCaloSamples::presamples(), edm::shift, IntegerCaloSamples::size(), theTowerMapFGSum, and TriggerAnalyzer::veto.

Referenced by 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<bool> finegrain(tpSamples, false);
 
    int ietaGCT = converttoGCTetaFromTPG(ieta);
    int iphiGCT = converttoGCTphiFromTPG(iphi); 
  
    for (int index = 0; index < tpSamples; ++index)
    {
      finegrain[index] = FGBit[iphiGCT][ietaGCT];
    }
 
    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];
      }
       }
    }
 
    IntegerCaloSamples output(samples.id(), tpSamples);
    output.setPresamples(tpPresamples);
 
    for (int ibin = 0; ibin < tpSamples; ++ibin) {
       int idx = ibin + shift;
       output[ibin] = samples[idx] / (rctlsb == 0.25 ? 4 : 8);
       if (output[ibin] > 0x3FF) output[ibin] = 0x3FF;
    }
    outcoder_->compress(output, finegrain, result);
 }

void HcalTriggerPrimitiveAlgo::run	(	const HcalTPGCoder *	incoder,
		const HcalTPGCompressor *	outcoder,
		const HBHEDigiCollection &	hbheDigis,
		const HFDigiCollection &	hfDigis,
		HcalTrigPrimDigiCollection &	result,
		const HcalTrigTowerGeometry *	trigTowerGeometry,
		float	rctlsb
	)

Definition at line 402 of file HcalTriggerPrimitiveAlgo.cc.

References addSignal(), analyze(), analyzeHF(), edm::SortedCollection< T, SORT >::back(), edm::SortedCollection< T, SORT >::begin(), converttoGCTeta(), converttoGCTphi(), edm::SortedCollection< T, SORT >::end(), FG_threshold1_, FG_threshold2_, fgMap_, HcalTrigTowerGeometry::firstHFTower(), HF_Veto, i, HcalTrigTowerDetId::ieta(), HcalDetId::ieta(), incoder_, HcalDetId::iphi(), j, relval_steps::k, outcoder_, edm::SortedCollection< T, SORT >::push_back(), theSumMap, theTowerMapFGSum, and theTrigTowerGeometry.

Referenced by HcalTrigPrimDigiProducer::produce().

                                                  {
    theTrigTowerGeometry = trigTowerGeometry;
     
    incoder_=dynamic_cast<const HcaluLUTTPGCoder*>(incoder);
    outcoder_=outcoder;
 
    theSumMap.clear();
    theTowerMapFGSum.clear();
    HF_Veto.clear();
    fgMap_.clear();
 
    // do the HB/HE digis
    for(HBHEDigiCollection::const_iterator hbheItr = hbheDigis.begin();
    hbheItr != hbheDigis.end(); ++hbheItr) {
       addSignal(*hbheItr);
    }
 
    // and the HF digis
    for(HFDigiCollection::const_iterator hfItr = hfDigis.begin();
    hfItr != hfDigis.end(); ++hfItr) {
       addSignal(*hfItr);
 
    }
 
    int phiBin = 18;
    int etaBin = 22;
    int thresholds = 2;
 
    int fiberQIEThresh[phiBin][etaBin][thresholds];
    int FGBit[18][22];
 
    for (int i = 0; i < phiBin; i++)
    {
      for (int j = 0; j < etaBin; j++)
      {
        FGBit[i][j] = 0;
        for (int k = 0; k < thresholds; k++)
        {
          fiberQIEThresh[i][j][k] = 0;
        }
      }
    }
 
    for(HFDigiCollection::const_iterator hfdigi = hfDigis.begin(); hfdigi != hfDigis.end(); hfdigi++)
    {
      HcalDetId hcalid = HcalDetId(hfdigi->id()) ;
 
      int ieta = hcalid.ieta();
      int iphi = hcalid.iphi();
      int presample = hfdigi->presamples();
  
      double qieval = hfdigi->sample(presample).adc();
  
      int GCTeta = converttoGCTeta(ieta);
      int GCTphi = converttoGCTphi(iphi);
 
      if (qieval > FG_threshold1_ )
      {
        fiberQIEThresh[GCTphi][GCTeta][0]++;
      }
      if (qieval > FG_threshold2_)
      {
        fiberQIEThresh[GCTphi][GCTeta][1]++; 
      }
    }
 
    // Calculate the final FG bit
    for (int i = 0; i < phiBin; i++)
    {
      for (int j = 0; j < etaBin; j++)
      {
        if (fiberQIEThresh[i][j][0] > 0 && j < 4)
        {
          FGBit[i][0] = 1;
          FGBit[i][2] = 1;
        }
        if (fiberQIEThresh[i][j][0] > 0 && j > 17)
        {
          FGBit[i][19] = 1;
          FGBit[i][21] = 1;
        }
        if (fiberQIEThresh[i][j][1] > 0 && j < 4)
        {
          FGBit[i][1] = 1;
          FGBit[i][3] = 1;
        }
        if (fiberQIEThresh[i][j][1] > 0 && j > 17)
        {
          FGBit[i][18] = 1;
          FGBit[i][20] = 1;
        }
      }
    }
 
  
   for(SumMap::iterator mapItr = theSumMap.begin(); mapItr != theSumMap.end(); ++mapItr) {
       result.push_back(HcalTriggerPrimitiveDigi(mapItr->first));
       HcalTrigTowerDetId detId(mapItr->second.id());
       if(detId.ietaAbs() >= theTrigTowerGeometry->firstHFTower())
       { 
         int ietaTow = detId.ieta();
     int iphiTow = detId.iphi();
     analyzeHF(ietaTow, iphiTow, FGBit, mapItr->second, result.back(), rctlsb);
       }
       else{ analyze(mapItr->second, result.back()); }
    }        
    return;
 }

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

Definition at line 765 of file HcalTriggerPrimitiveAlgo.cc.

References edm::SortedCollection< T, SORT >::begin(), HcalHTRData::check(), FEDRawData::data(), DetId::det(), edm::SortedCollection< T, SORT >::end(), FEDRawDataCollection::FEDData(), HcalHTRData::getErrorsWord(), HcalDCCHeader::getSourceId(), HcalDCCHeader::getSpigotData(), HcalDCCHeader::getSpigotPresent(), HcalBarrel, HcalEndcap, HcalForward, HcalHTRData::htrSlot(), HcalHTRData::htrTopBottom(), i, 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 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 751 of file HcalTriggerPrimitiveAlgo.cc.

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

Referenced by 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::setPeakFinderAlgorithm ( int algo )

Definition at line 827 of file HcalTriggerPrimitiveAlgo.cc.

References ecalcalib_dqm_sourceclient-live_cfg::algo, Exception, and peak_finder_algorithm_.

Referenced by HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

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