EcalEBTrigPrimTestAlgo Class Reference

#include <EcalEBTrigPrimTestAlgo.h>

Public Member Functions
	EcalEBTrigPrimTestAlgo (const EcalTrigTowerConstituentsMap *eTTmap, const CaloGeometry *theGeometry, int nSamples, int binofmax, bool tcpFormat, bool debug, bool famos)
	EcalEBTrigPrimTestAlgo (int nSamples, int binofmax, bool tcpFormat, bool debug, bool famos)
void	run (const EBDigiCollection *col, EcalEBTrigPrimDigiCollection &result, EcalEBTrigPrimDigiCollection &resultTcp)
void	setPointers (const EcalTPGLinearizationConst *ecaltpLin, const EcalTPGPedestals *ecaltpPed, const EcalTPGCrystalStatus *ecaltpgBadX, const EcalTPGWeightIdMap *ecaltpgWeightMap, const EcalTPGWeightGroup *ecaltpgWeightGroup, const EcalTPGSlidingWindow *ecaltpgSlidW, const EcalTPGLutGroup *ecaltpgLutGroup, const EcalTPGLutIdMap *ecaltpgLut, const EcalTPGTowerStatus *ecaltpgBadTT, const EcalTPGSpike *ecaltpgSpike)
virtual	~EcalEBTrigPrimTestAlgo ()

Private Types
enum	{ nbMaxStrips_ = 5 }
enum	{ nbMaxXtals_ = 5 }

Private Member Functions
template<class T >
void	clean (std::vector< std::vector< std::pair< int, std::vector< T > > > > &towerMap)
template<class Coll >
void	fillMap (Coll const *col, std::vector< std::vector< std::pair< int, std::vector< typename Coll::Digi > > > > &towerMap)
int	findStripNr (const EBDetId &id)
int	findStripNr (const EEDetId &id)
EcalEBFenixAmplitudeFilter *	getFilter () const
EcalEBFenixTcpFormat *	getFormatter () const
EcalEBFenixStripFormatEB *	getFormatterEB () const
int	getIndex (const EBDigiCollection *, EcalTrigTowerDetId &id)
int	getIndex (const EEDigiCollection *, EcalTrigTowerDetId &id)
EcalEBFenixLinearizer *	getLinearizer (int i) const
EcalEBFenixPeakFinder *	getPeakFinder () const
void	init ()
template<class T >
void	initStructures (std::vector< std::vector< std::pair< int, std::vector< T > > > > &towMap)

Private Attributes
EcalEBFenixAmplitudeFilter *	amplitude_filter_
bool	barrelOnly_
int	binOfMaximum_
bool	debug_
const EcalTPGTowerStatus *	ecaltpgBadTT_
const EcalTPGCrystalStatus *	ecaltpgBadX_
const EcalTPGLutIdMap *	ecaltpgLut_
const EcalTPGLutGroup *	ecaltpgLutGroup_
const EcalTPGSlidingWindow *	ecaltpgSlidW_
const EcalTPGSpike *	ecaltpgSpike_
const EcalTPGWeightGroup *	ecaltpgWeightGroup_
const EcalTPGWeightIdMap *	ecaltpgWeightMap_
const EcalTPGLinearizationConst *	ecaltpLin_
const EcalTPGPedestals *	ecaltpPed_
const EcalTrigTowerConstituentsMap *	eTTmap_ = nullptr
bool	famos_
EcalEBFenixStripFormatEB *	fenixFormatterEB_
EcalEBFenixTcpFormat *	fenixTcpFormat_
std::vector< int >	fgvb_out_
std::vector< int >	fgvb_out_temp_
std::vector< int >	filt_out_
std::vector< int >	format_out_
std::vector< std::pair< int, EcalTrigTowerDetId > >	hitTowers_
std::vector< std::vector< int > >	lin_out_
std::vector< EcalEBFenixLinearizer * >	linearizer_
int	maxNrSamples_
int	nrTowers_
int	nSamples_
EcalEBFenixPeakFinder *	peak_finder_
std::vector< int >	peak_out_
std::vector< std::vector< int > >	striptp_
bool	tcpFormat_
std::vector< int >	tcpformat_out_
const CaloGeometry *	theGeometry_ = nullptr
const EcalElectronicsMapping *	theMapping_
float	threshold
std::vector< std::vector< std::pair< int, std::vector< EBDataFrame > > > >	towerMapEB_
std::vector< std::vector< std::pair< int, std::vector< EEDataFrame > > > >	towerMapEE_
std::vector< EcalEBTriggerPrimitiveSample >	towtp2_
std::vector< EcalEBTriggerPrimitiveSample >	towtp_

Static Private Attributes
static const unsigned int	maxNrSamplesOut_ = 10
static const unsigned int	maxNrTowers_ = 2448
static const unsigned int	nrSamples_ = 5

Detailed Description

Author

N. Marinelli - Univ. of Notre Dame forPhase II While the new digitization is not yet implemented, we use the old Digis to make TP per crystal

EcalEBTrigPrimTestAlgo starting point for Phase II: build TPs out of Phase I digis to start building the infrastructures

Definition at line 37 of file EcalEBTrigPrimTestAlgo.h.

Member Enumeration Documentation

anonymous enum

anonymous enum

private

Enumerator
nbMaxStrips_

Definition at line 124 of file EcalEBTrigPrimTestAlgo.h.

{ nbMaxStrips_ = 5 };

anonymous enum

anonymous enum

private

Enumerator
nbMaxXtals_

Definition at line 125 of file EcalEBTrigPrimTestAlgo.h.

{ nbMaxXtals_ = 5 };

Constructor & Destructor Documentation

EcalEBTrigPrimTestAlgo() [1/2]

EcalEBTrigPrimTestAlgo::EcalEBTrigPrimTestAlgo ( const EcalTrigTowerConstituentsMap *  eTTmap, const CaloGeometry *  theGeometry, int  nSamples, int  binofmax, bool  tcpFormat, bool  debug, bool  famos  )

explicit

Definition at line 43 of file EcalEBTrigPrimTestAlgo.cc.

References init(), and maxNrSamples_.

    : eTTmap_(eTTmap),
      theGeometry_(theGeometry),
      nSamples_(nSam),
      binOfMaximum_(binofmax),
      tcpFormat_(tcpFormat),
      barrelOnly_(false),
      debug_(debug),
      famos_(famos) {
  maxNrSamples_ = 10;
  this->init();
}

EcalEBTrigPrimTestAlgo() [2/2]

EcalEBTrigPrimTestAlgo::EcalEBTrigPrimTestAlgo ( int  nSamples, int  binofmax, bool  tcpFormat, bool  debug, bool  famos  )

explicit

Definition at line 63 of file EcalEBTrigPrimTestAlgo.cc.

References init(), and maxNrSamples_.

    : nSamples_(nSam),
      binOfMaximum_(binofmax),
      tcpFormat_(tcpFormat),
      barrelOnly_(true),
      debug_(debug),
      famos_(famos)

{
  maxNrSamples_ = 10;
  this->init();
}

~EcalEBTrigPrimTestAlgo()

EcalEBTrigPrimTestAlgo::~EcalEBTrigPrimTestAlgo ( )

virtual

Definition at line 110 of file EcalEBTrigPrimTestAlgo.cc.

References amplitude_filter_, fenixFormatterEB_, fenixTcpFormat_, mps_fire::i, linearizer_, nbMaxXtals_, peak_finder_, and theMapping_.

                                                {
  for (int i = 0; i < nbMaxXtals_; i++)
    delete linearizer_[i];
  delete amplitude_filter_;
  delete peak_finder_;
  delete fenixFormatterEB_;
  delete fenixTcpFormat_;
  delete theMapping_;
}

Member Function Documentation

clean()

template<class T >

void EcalEBTrigPrimTestAlgo::clean ( std::vector< std::vector< std::pair< int, std::vector< T > > > > &  towerMap )

private

Definition at line 167 of file EcalEBTrigPrimTestAlgo.h.

References mps_fire::i, dqmiolumiharvest::j, maxNrTowers_, and nbMaxStrips_.

Referenced by run().

                                                                                              {
  // clean internal data structures
  for (unsigned int i = 0; i < maxNrTowers_; ++i)
    for (int j = 0; j < nbMaxStrips_; ++j)
      (towMap[i])[j].first = 0;
  return;
}

fillMap()

template<class Coll >

void EcalEBTrigPrimTestAlgo::fillMap ( Coll const *  col, std::vector< std::vector< std::pair< int, std::vector< typename Coll::Digi > > > > &  towerMap  )

private

Definition at line 176 of file EcalEBTrigPrimTestAlgo.h.

References EcalTrigTowerConstituentsMap::barrelTowerOf(), cuy::col, gather_cfg::cout, debug_, TrendClient_cfi::Digi, eTTmap_, findStripNr(), dqmdumpme::first, getIndex(), hitTowers_, mps_fire::i, nbMaxXtals_, nrTowers_, EgammaValidation_cff::samples, edm::second(), and EcalTrigTowerConstituentsMap::towerOf().

Referenced by run().

                                                                                                    {
  typedef typename Coll::Digi Digi;

  // implementation for Barrel
  if (col) {
    nrTowers_ = 0;
    if (debug_)
      std::cout << "Fill mapping, Collection size = " << col->size() << std::endl;
    ;
    for (unsigned int i = 0; i < col->size(); ++i) {
      Digi samples((*col)[i]);
      EcalTrigTowerDetId coarser =
          eTTmap_ ? (*eTTmap_).towerOf(samples.id()) : EcalTrigTowerConstituentsMap::barrelTowerOf(samples.id());
      int index = getIndex(col, coarser);
      int stripnr = findStripNr(samples.id());

      int filled = 0;
      for (unsigned int ij = 0; ij < towerMap[index].size(); ++ij)
        filled += towerMap[index][ij].first;
      if (!filled) {
        hitTowers_[nrTowers_++] = std::pair<int, EcalTrigTowerDetId>(index, coarser);
      }

      //FIXME: temporary protection
      int ncryst = towerMap[index][stripnr - 1].first;
      if (ncryst >= nbMaxXtals_) {
        edm::LogError("EcalTrigPrimFunctionAlgo")
            << "! Too many xtals for TT " << coarser << " stripnr " << stripnr << " xtalid " << samples.id();
        continue;
      }
      ((towerMap[index])[stripnr - 1].second)[ncryst] = samples;
      (towerMap[index])[stripnr - 1].first++;
    }

    if (debug_)
      std::cout << "fillMap"
                << "[EcalTrigPrimFunctionalAlgo] (found " << col->size() << " frames in " << towerMap.size()
                << " towers) " << std::endl;
  } else {
    if (debug_)
      std::cout << "FillMap - FillMap Collection size=0 !!!!" << std::endl;
    ;
  }
}

findStripNr() [1/2]

int EcalEBTrigPrimTestAlgo::findStripNr ( const EBDetId &  id )

private

Definition at line 319 of file EcalEBTrigPrimTestAlgo.cc.

References hcalRecHitTable_cff::ieta, dqmiodumpmetadata::n, and nbMaxStrips_.

Referenced by fillMap().

                                                         {
  int stripnr;
  int n = ((id.ic() - 1) % 100) / 20;  //20 corresponds to 4 * ecal_barrel_crystals_per_strip FIXME!!
  if (id.ieta() < 0)
    stripnr = n + 1;
  else
    stripnr = nbMaxStrips_ - n;
  return stripnr;
}

findStripNr() [2/2]

int EcalEBTrigPrimTestAlgo::findStripNr ( const EEDetId &  id )

private

getFilter()

EcalEBFenixAmplitudeFilter* EcalEBTrigPrimTestAlgo::getFilter ( ) const

inlineprivate

Definition at line 150 of file EcalEBTrigPrimTestAlgo.h.

References amplitude_filter_.

Referenced by run().

{ return amplitude_filter_; }

getFormatter()

EcalEBFenixTcpFormat* EcalEBTrigPrimTestAlgo::getFormatter ( ) const

inlineprivate

Definition at line 162 of file EcalEBTrigPrimTestAlgo.h.

References fenixTcpFormat_.

Referenced by run().

{ return fenixTcpFormat_; }

getFormatterEB()

EcalEBFenixStripFormatEB* EcalEBTrigPrimTestAlgo::getFormatterEB ( ) const

inlineprivate

Definition at line 160 of file EcalEBTrigPrimTestAlgo.h.

References fenixFormatterEB_.

Referenced by run().

{ return fenixFormatterEB_; }

getIndex() [1/2]

int EcalEBTrigPrimTestAlgo::getIndex ( const EBDigiCollection *  , EcalTrigTowerDetId &  id  )

inlineprivate

Definition at line 89 of file EcalEBTrigPrimTestAlgo.h.

Referenced by fillMap().

{ return id.hashedIndex(); }

getIndex() [2/2]

int EcalEBTrigPrimTestAlgo::getIndex ( const EEDigiCollection *  , EcalTrigTowerDetId &  id  )

inlineprivate

Definition at line 91 of file EcalEBTrigPrimTestAlgo.h.

References hcalRecHitTable_cff::iphi, and ecaldqm::zside().

                                                                 {
    int ind = (id.ietaAbs() - 18) * 72 + id.iphi();
    if (id.zside() < 0)
      ind += 792;
    return ind;
  }

getLinearizer()

EcalEBFenixLinearizer* EcalEBTrigPrimTestAlgo::getLinearizer ( int  i ) const

inlineprivate

Definition at line 147 of file EcalEBTrigPrimTestAlgo.h.

References mps_fire::i, and linearizer_.

Referenced by run().

{ return linearizer_[i]; }

getPeakFinder()

EcalEBFenixPeakFinder* EcalEBTrigPrimTestAlgo::getPeakFinder ( ) const

inlineprivate

Definition at line 159 of file EcalEBTrigPrimTestAlgo.h.

References peak_finder_.

Referenced by run().

{ return peak_finder_; }

init()

void EcalEBTrigPrimTestAlgo::init ( void  )

private

Definition at line 77 of file EcalEBTrigPrimTestAlgo.cc.

References amplitude_filter_, binOfMaximum_, debug_, famos_, fenixFormatterEB_, fenixTcpFormat_, fgvb_out_, fgvb_out_temp_, filt_out_, format_out_, hitTowers_, mps_fire::i, initStructures(), lin_out_, linearizer_, maxNrSamples_, maxNrTowers_, nbMaxXtals_, peak_finder_, peak_out_, tcpFormat_, tcpformat_out_, theMapping_, towerMapEB_, and findQualityFiles::v.

Referenced by EcalEBTrigPrimTestAlgo().

                                  {
  // initialise data structures
  theMapping_ = new EcalElectronicsMapping();
  initStructures(towerMapEB_);
  hitTowers_.resize(maxNrTowers_);

  linearizer_.resize(nbMaxXtals_);
  for (int i = 0; i < nbMaxXtals_; i++)
    linearizer_[i] = new EcalEBFenixLinearizer(famos_);

  //
  std::vector<int> v;
  v.resize(maxNrSamples_);
  lin_out_.resize(nbMaxXtals_);
  for (int i = 0; i < 5; i++)
    lin_out_[i] = v;
  //
  amplitude_filter_ = new EcalEBFenixAmplitudeFilter();
  filt_out_.resize(maxNrSamples_);
  peak_out_.resize(maxNrSamples_);
  // these two are dummy
  fgvb_out_.resize(maxNrSamples_);
  fgvb_out_temp_.resize(maxNrSamples_);
  //
  peak_finder_ = new EcalEBFenixPeakFinder();
  fenixFormatterEB_ = new EcalEBFenixStripFormatEB();
  format_out_.resize(maxNrSamples_);
  //
  fenixTcpFormat_ = new EcalEBFenixTcpFormat(tcpFormat_, debug_, famos_, binOfMaximum_);
  tcpformat_out_.resize(maxNrSamples_);
}

initStructures()

template<class T >

void EcalEBTrigPrimTestAlgo::initStructures ( std::vector< std::vector< std::pair< int, std::vector< T > > > > &  towMap )

private

Definition at line 223 of file EcalEBTrigPrimTestAlgo.h.

References mps_fire::i, maxNrSamples_, maxNrTowers_, nbMaxStrips_, nbMaxXtals_, and striptp_.

Referenced by init().

                                                                                                       {
  //initialise internal data structures

  std::vector<T> vec0(nbMaxXtals_);
  std::vector<std::pair<int, std::vector<T> > > vec1(nbMaxStrips_);
  for (int i = 0; i < nbMaxStrips_; ++i)
    vec1[i] = std::pair<int, std::vector<T> >(0, vec0);
  towMap.resize(maxNrTowers_);
  for (unsigned int i = 0; i < maxNrTowers_; ++i)
    towMap[i] = vec1;

  std::vector<int> vecint(maxNrSamples_);
  striptp_.resize(nbMaxStrips_);
  for (int i = 0; i < nbMaxStrips_; ++i)
    striptp_[i] = vecint;
}

run()

void EcalEBTrigPrimTestAlgo::run	(	const EBDigiCollection *	col,
		EcalEBTrigPrimDigiCollection &	result,
		EcalEBTrigPrimDigiCollection &	resultTcp
	)

Definition at line 120 of file EcalEBTrigPrimTestAlgo.cc.

References binOfMaximum_, clean(), gather_cfg::cout, debug_, TauDecayModes::dec, ecaltpgBadTT_, ecaltpgBadX_, ecaltpgLut_, ecaltpgLutGroup_, ecaltpgSlidW_, ecaltpgSpike_, ecaltpgWeightGroup_, ecaltpgWeightMap_, ecaltpLin_, ecaltpPed_, fgvb_out_, fgvb_out_temp_, fillMap(), filt_out_, dqmdumpme::first, castor_dqm_sourceclient-live_cfg::firstSample, format_out_, getFilter(), getFormatter(), getFormatterEB(), getLinearizer(), getPeakFinder(), EcalElectronicsMapping::getTriggerElectronicsId(), hitTowers_, mps_fire::i, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), CastorRawToDigi_cfi::lastSample, lin_out_, nrSamples_, nrTowers_, peak_out_, EcalEBFenixPeakFinder::process(), EcalEBFenixAmplitudeFilter::process(), EcalEBFenixTcpFormat::process(), EcalEBFenixStripFormatEB::process(), EcalEBFenixLinearizer::process(), edm::SortedCollection< T, SORT >::push_back(), EcalTriggerElectronicsId::rawId(), nano_mu_digi_cff::rawId, DetId::rawId(), mps_fire::result, edm::second(), EcalEBFenixStripFormatEB::setParameters(), EcalEBFenixTcpFormat::setParameters(), EcalEBFenixAmplitudeFilter::setParameters(), EcalEBFenixLinearizer::setParameters(), edm::DataFrameContainer::size(), tcpFormat_, tcpformat_out_, theMapping_, ecalMatacq_cfi::timing, towerMapEB_, and cmsswSequenceInfo::tp.

                                                                          {
  //typedef typename Coll::Digi Digi;
  if (debug_) {
    std::cout << "  EcalEBTrigPrimTestAlgo: Testing that the algorythm with digis is well plugged " << std::endl;
    std::cout << "  EcalEBTrigPrimTestAlgo: digi size " << digi->size() << std::endl;
  }

  uint16_t etInADC;
  EcalEBTriggerPrimitiveDigi tp;
  int firstSample = binOfMaximum_ - 1 - nrSamples_ / 2;
  int lastSample = binOfMaximum_ - 1 + nrSamples_ / 2;

  if (debug_) {
    std::cout << "  binOfMaximum_ " << binOfMaximum_ << " nrSamples_" << nrSamples_ << std::endl;
    std::cout << " first sample " << firstSample << " last " << lastSample << std::endl;
  }

  clean(towerMapEB_);
  fillMap(digi, towerMapEB_);

  for (int itow = 0; itow < nrTowers_; ++itow) {
    int index = hitTowers_[itow].first;
    const EcalTrigTowerDetId &thisTower = hitTowers_[itow].second;
    if (debug_)
      std::cout << " Data for TOWER num " << itow << " index " << index << " TowerId " << thisTower << " size "
                << towerMapEB_[itow].size() << std::endl;
    // loop over all strips assigned to this trigger tower
    int nxstals = 0;
    for (unsigned int iStrip = 0; iStrip < towerMapEB_[itow].size(); ++iStrip) {
      if (debug_)
        std::cout << " Data for STRIP num " << iStrip << std::endl;
      std::vector<EBDataFrame> &dataFrames =
          (towerMapEB_[index])[iStrip].second;  //vector of dataframes for this strip, size; nr of crystals/strip

      nxstals = (towerMapEB_[index])[iStrip].first;
      if (nxstals <= 0)
        continue;
      if (debug_)
        std::cout << " Number of xTals " << nxstals << std::endl;

      const EcalTriggerElectronicsId elId = theMapping_->getTriggerElectronicsId(dataFrames[0].id());
      uint32_t stripid = elId.rawId() & 0xfffffff8;

      // loop over the xstals in a strip
      for (int iXstal = 0; iXstal < nxstals; iXstal++) {
        const EBDetId &myid = dataFrames[iXstal].id();
        tp = EcalEBTriggerPrimitiveDigi(myid);
        tp.setSize(nrSamples_);

        if (debug_) {
          std::cout << std::endl;
          std::cout << "iXstal= " << iXstal << " id " << dataFrames[iXstal].id() << " EBDataFrame is: " << std::endl;
          for (int i = 0; i < dataFrames[iXstal].size(); i++) {
            std::cout << " " << std::dec << dataFrames[iXstal][i].adc();
          }
          std::cout << std::endl;
        }
        //   Call the linearizer
        this->getLinearizer(iXstal)->setParameters(
            dataFrames[iXstal].id().rawId(), ecaltpPed_, ecaltpLin_, ecaltpgBadX_);
        this->getLinearizer(iXstal)->process(dataFrames[iXstal], lin_out_[iXstal]);

        for (unsigned int i = 0; i < lin_out_[iXstal].size(); i++) {
          if ((lin_out_[iXstal])[i] > 0X3FFFF)
            (lin_out_[iXstal])[i] = 0X3FFFF;
        }

        if (debug_) {
          std::cout << "output of linearizer for channel " << iXstal << std::endl;
          std::cout << " lin_out[iXstal].size()= " << std::dec << lin_out_[iXstal].size() << std::endl;
          for (unsigned int i = 0; i < lin_out_[iXstal].size(); i++) {
            std::cout << " " << std::dec << (lin_out_[iXstal])[i];
          }
          std::cout << std::endl;
        }

        // Call the amplitude filter
        this->getFilter()->setParameters(stripid, ecaltpgWeightMap_, ecaltpgWeightGroup_);
        this->getFilter()->process(lin_out_[iXstal], filt_out_, fgvb_out_temp_, fgvb_out_);

        if (debug_) {
          std::cout << "output of filter is a vector of size: " << std::dec << filt_out_.size() << std::endl;
          for (unsigned int ix = 0; ix < filt_out_.size(); ix++) {
            std::cout << std::dec << filt_out_[ix] << " ";
          }
          std::cout << std::endl;
        }

        // call peakfinder
        this->getPeakFinder()->process(filt_out_, peak_out_);

        if (debug_) {
          std::cout << "output of peakfinder is a vector of size: " << std::dec << peak_out_.size() << std::endl;
          for (unsigned int ix = 0; ix < peak_out_.size(); ix++) {
            std::cout << std::dec << peak_out_[ix] << " ";
          }
          std::cout << std::endl;
        }

        // call formatter
        this->getFormatterEB()->setParameters(stripid, ecaltpgSlidW_);
        this->getFormatterEB()->process(fgvb_out_, peak_out_, filt_out_, format_out_);

        if (debug_) {
          std::cout << "output of formatter is a vector of size: " << format_out_.size() << std::endl;
          for (unsigned int i = 0; i < format_out_.size(); i++) {
            std::cout << " " << std::dec << format_out_[i] << " ";
          }
          std::cout << std::endl;
        }

        // call final tcp formatter
        this->getFormatter()->setParameters(
            thisTower.rawId(), ecaltpgLutGroup_, ecaltpgLut_, ecaltpgBadTT_, ecaltpgSpike_);
        this->getFormatter()->process(format_out_, tcpformat_out_);
        // loop over the time samples and fill the TP
        int nSam = 0;
        for (int iSample = firstSample; iSample <= lastSample; ++iSample) {
          etInADC = tcpformat_out_[iSample];
          if (debug_)
            std::cout << " format_out " << tcpformat_out_[iSample] << " etInADC " << etInADC << std::endl;

          bool isASpike = false;  // no spikes for now
          int timing = 0;         //  set to 0  value for now
          tp.setSample(nSam, EcalEBTriggerPrimitiveSample(etInADC, isASpike, timing));

          nSam++;
          if (debug_)
            std::cout << "in TestAlgo"
                      << " tp size " << tp.size() << std::endl;
        }

        if (!tcpFormat_)
          result.push_back(tp);
        else
          resultTcp.push_back(tp);

      }  // Loop over the xStals

    }  //loop over strips in one tower
  }

  /*
  for (unsigned int i=0;i<digi->size();i++) {
    EBDataFrame myFrame((*digi)[i]);  
    const EBDetId & myid1 = myFrame.id();
    tp=  EcalTriggerPrimitiveDigi(  myid1);   
    tp.setSize( myFrame.size());
    int nSam=0;

    if (debug_) {
      std::cout << " data frame size " << myFrame.size() << " Id " <<  myFrame.id()  << std::endl;
      std::cout << " Sample data ADC: " << std::endl;
      for (int iSample=0; iSample<myFrame.size(); iSample++) {
    std::cout << " " << std::dec<< myFrame.sample(iSample).adc() ;
      }
      std::cout<<std::endl;
    }

    
    this->getLinearizer(i)->setParameters( myFrame.id().rawId(),ecaltpPed_,ecaltpLin_,ecaltpgBadX_) ; 
    //this->getLinearizer(i)->process( myFrame,lin_out_[i]);

    if (debug_) {
      std::cout<< "cryst: "<< i <<"  value : "<<std::dec<<std::endl;
      std::cout<<" lin_out[i].size()= "<<std::dec<<lin_out_[i].size()<<std::endl;
      for (unsigned int j =0; j<lin_out_[i].size();j++){
    std::cout <<" "<<std::dec<<(lin_out_[i])[j];
      }
      std::cout<<std::endl;
    }


    for (int iSample=0; iSample<myFrame.size(); iSample++) {
      etInADC= myFrame.sample(iSample).adc();
      EcalEBTriggerPrimitiveSample mysam(etInADC);
      tp.setSample(nSam, mysam );
      nSam++;
      if (debug_) std::cout << "in TestAlgo" <<" tp size "<<tp.size() << std::endl;
    }

    if (!tcpFormat_)
      result.push_back(tp);
    else 
      resultTcp.push_back(tp);
    
   
    if (debug_) std::cout << " result size " << result.size() << std::endl;
    
    
    
  }
  */
}

setPointers()

void EcalEBTrigPrimTestAlgo::setPointers ( const EcalTPGLinearizationConst *  ecaltpLin, const EcalTPGPedestals *  ecaltpPed, const EcalTPGCrystalStatus *  ecaltpgBadX, const EcalTPGWeightIdMap *  ecaltpgWeightMap, const EcalTPGWeightGroup *  ecaltpgWeightGroup, const EcalTPGSlidingWindow *  ecaltpgSlidW, const EcalTPGLutGroup *  ecaltpgLutGroup, const EcalTPGLutIdMap *  ecaltpgLut, const EcalTPGTowerStatus *  ecaltpgBadTT, const EcalTPGSpike *  ecaltpgSpike  )

inline

Definition at line 54 of file EcalEBTrigPrimTestAlgo.h.

References ecaltpgBadTT_, ecaltpgBadX_, ecaltpgLut_, ecaltpgLutGroup_, ecaltpgSlidW_, ecaltpgSpike_, ecaltpgWeightGroup_, ecaltpgWeightMap_, ecaltpLin_, and ecaltpPed_.

                                                     {
    ecaltpPed_ = ecaltpPed;
    ecaltpLin_ = ecaltpLin;
    ecaltpgBadX_ = ecaltpgBadX;
    ecaltpgWeightMap_ = ecaltpgWeightMap;
    ecaltpgWeightGroup_ = ecaltpgWeightGroup;
    ecaltpgSlidW_ = ecaltpgSlidW;
    ecaltpgLutGroup_ = ecaltpgLutGroup;
    ecaltpgLut_ = ecaltpgLut;
    ecaltpgBadTT_ = ecaltpgBadTT;
    ecaltpgSpike_ = ecaltpgSpike;
  }

Member Data Documentation

amplitude_filter_

EcalEBFenixAmplitudeFilter* EcalEBTrigPrimTestAlgo::amplitude_filter_

private

Definition at line 130 of file EcalEBTrigPrimTestAlgo.h.

Referenced by getFilter(), init(), and ~EcalEBTrigPrimTestAlgo().

barrelOnly_

bool EcalEBTrigPrimTestAlgo::barrelOnly_

private

Definition at line 107 of file EcalEBTrigPrimTestAlgo.h.

binOfMaximum_

int EcalEBTrigPrimTestAlgo::binOfMaximum_

private

Definition at line 103 of file EcalEBTrigPrimTestAlgo.h.

Referenced by init(), and run().

debug_

bool EcalEBTrigPrimTestAlgo::debug_

private

Definition at line 108 of file EcalEBTrigPrimTestAlgo.h.

Referenced by fillMap(), init(), and run().

ecaltpgBadTT_

const EcalTPGTowerStatus* EcalEBTrigPrimTestAlgo::ecaltpgBadTT_

private

Definition at line 144 of file EcalEBTrigPrimTestAlgo.h.

Referenced by run(), and setPointers().

ecaltpgBadX_

const EcalTPGCrystalStatus* EcalEBTrigPrimTestAlgo::ecaltpgBadX_

private

Definition at line 138 of file EcalEBTrigPrimTestAlgo.h.

Referenced by run(), and setPointers().

ecaltpgLut_

const EcalTPGLutIdMap* EcalEBTrigPrimTestAlgo::ecaltpgLut_

private

Definition at line 143 of file EcalEBTrigPrimTestAlgo.h.

Referenced by run(), and setPointers().

ecaltpgLutGroup_

const EcalTPGLutGroup* EcalEBTrigPrimTestAlgo::ecaltpgLutGroup_

private

Definition at line 142 of file EcalEBTrigPrimTestAlgo.h.

Referenced by run(), and setPointers().

ecaltpgSlidW_

const EcalTPGSlidingWindow* EcalEBTrigPrimTestAlgo::ecaltpgSlidW_

private

Definition at line 141 of file EcalEBTrigPrimTestAlgo.h.

Referenced by run(), and setPointers().

ecaltpgSpike_

const EcalTPGSpike* EcalEBTrigPrimTestAlgo::ecaltpgSpike_

private

Definition at line 145 of file EcalEBTrigPrimTestAlgo.h.

Referenced by run(), and setPointers().

ecaltpgWeightGroup_

const EcalTPGWeightGroup* EcalEBTrigPrimTestAlgo::ecaltpgWeightGroup_

private

Definition at line 140 of file EcalEBTrigPrimTestAlgo.h.

Referenced by run(), and setPointers().

ecaltpgWeightMap_

const EcalTPGWeightIdMap* EcalEBTrigPrimTestAlgo::ecaltpgWeightMap_

private

Definition at line 139 of file EcalEBTrigPrimTestAlgo.h.

Referenced by run(), and setPointers().

ecaltpLin_

const EcalTPGLinearizationConst* EcalEBTrigPrimTestAlgo::ecaltpLin_

private

Definition at line 137 of file EcalEBTrigPrimTestAlgo.h.

Referenced by run(), and setPointers().

ecaltpPed_

const EcalTPGPedestals* EcalEBTrigPrimTestAlgo::ecaltpPed_

private

Definition at line 136 of file EcalEBTrigPrimTestAlgo.h.

Referenced by run(), and setPointers().

eTTmap_

const EcalTrigTowerConstituentsMap* EcalEBTrigPrimTestAlgo::eTTmap_ = nullptr

private

Definition at line 98 of file EcalEBTrigPrimTestAlgo.h.

Referenced by fillMap().

famos_

bool EcalEBTrigPrimTestAlgo::famos_

private

Definition at line 109 of file EcalEBTrigPrimTestAlgo.h.

Referenced by init().

fenixFormatterEB_

EcalEBFenixStripFormatEB* EcalEBTrigPrimTestAlgo::fenixFormatterEB_

private

Definition at line 132 of file EcalEBTrigPrimTestAlgo.h.

Referenced by getFormatterEB(), init(), and ~EcalEBTrigPrimTestAlgo().

fenixTcpFormat_

EcalEBFenixTcpFormat* EcalEBTrigPrimTestAlgo::fenixTcpFormat_

private

Definition at line 133 of file EcalEBTrigPrimTestAlgo.h.

Referenced by getFormatter(), init(), and ~EcalEBTrigPrimTestAlgo().

fgvb_out_

std::vector<int> EcalEBTrigPrimTestAlgo::fgvb_out_

private

Definition at line 155 of file EcalEBTrigPrimTestAlgo.h.

Referenced by init(), and run().

fgvb_out_temp_

std::vector<int> EcalEBTrigPrimTestAlgo::fgvb_out_temp_

private

Definition at line 156 of file EcalEBTrigPrimTestAlgo.h.

Referenced by init(), and run().

filt_out_

std::vector<int> EcalEBTrigPrimTestAlgo::filt_out_

private

Definition at line 151 of file EcalEBTrigPrimTestAlgo.h.

Referenced by init(), and run().

format_out_

std::vector<int> EcalEBTrigPrimTestAlgo::format_out_

private

Definition at line 153 of file EcalEBTrigPrimTestAlgo.h.

Referenced by init(), and run().

hitTowers_

std::vector<std::pair<int, EcalTrigTowerDetId> > EcalEBTrigPrimTestAlgo::hitTowers_

private

Definition at line 120 of file EcalEBTrigPrimTestAlgo.h.

Referenced by fillMap(), init(), and run().

lin_out_

std::vector<std::vector<int> > EcalEBTrigPrimTestAlgo::lin_out_

private

Definition at line 148 of file EcalEBTrigPrimTestAlgo.h.

Referenced by init(), and run().

linearizer_

std::vector<EcalEBFenixLinearizer *> EcalEBTrigPrimTestAlgo::linearizer_

private

Definition at line 129 of file EcalEBTrigPrimTestAlgo.h.

Referenced by getLinearizer(), init(), and ~EcalEBTrigPrimTestAlgo().

maxNrSamples_

int EcalEBTrigPrimTestAlgo::maxNrSamples_

private

Definition at line 104 of file EcalEBTrigPrimTestAlgo.h.

Referenced by EcalEBTrigPrimTestAlgo(), init(), and initStructures().

maxNrSamplesOut_

const unsigned int EcalEBTrigPrimTestAlgo::maxNrSamplesOut_ = 10

staticprivate

Definition at line 113 of file EcalEBTrigPrimTestAlgo.h.

maxNrTowers_

const unsigned int EcalEBTrigPrimTestAlgo::maxNrTowers_ = 2448

staticprivate

Definition at line 112 of file EcalEBTrigPrimTestAlgo.h.

Referenced by clean(), init(), and initStructures().

nrSamples_

const unsigned int EcalEBTrigPrimTestAlgo::nrSamples_ = 5

staticprivate

Definition at line 114 of file EcalEBTrigPrimTestAlgo.h.

Referenced by run().

nrTowers_

int EcalEBTrigPrimTestAlgo::nrTowers_

private

Definition at line 111 of file EcalEBTrigPrimTestAlgo.h.

Referenced by fillMap(), and run().

nSamples_

int EcalEBTrigPrimTestAlgo::nSamples_

private

Definition at line 102 of file EcalEBTrigPrimTestAlgo.h.

peak_finder_

EcalEBFenixPeakFinder* EcalEBTrigPrimTestAlgo::peak_finder_

private

Definition at line 131 of file EcalEBTrigPrimTestAlgo.h.

Referenced by getPeakFinder(), init(), and ~EcalEBTrigPrimTestAlgo().

peak_out_

std::vector<int> EcalEBTrigPrimTestAlgo::peak_out_

private

Definition at line 152 of file EcalEBTrigPrimTestAlgo.h.

Referenced by init(), and run().

striptp_

std::vector<std::vector<int> > EcalEBTrigPrimTestAlgo::striptp_

private

Definition at line 117 of file EcalEBTrigPrimTestAlgo.h.

Referenced by initStructures().

tcpFormat_

bool EcalEBTrigPrimTestAlgo::tcpFormat_

private

Definition at line 106 of file EcalEBTrigPrimTestAlgo.h.

Referenced by init(), and run().

tcpformat_out_

std::vector<int> EcalEBTrigPrimTestAlgo::tcpformat_out_

private

Definition at line 163 of file EcalEBTrigPrimTestAlgo.h.

Referenced by init(), and run().

theGeometry_

const CaloGeometry* EcalEBTrigPrimTestAlgo::theGeometry_ = nullptr

private

Definition at line 99 of file EcalEBTrigPrimTestAlgo.h.

theMapping_

const EcalElectronicsMapping* EcalEBTrigPrimTestAlgo::theMapping_

private

Definition at line 127 of file EcalEBTrigPrimTestAlgo.h.

Referenced by init(), run(), and ~EcalEBTrigPrimTestAlgo().

threshold

float EcalEBTrigPrimTestAlgo::threshold

private

Definition at line 101 of file EcalEBTrigPrimTestAlgo.h.

Referenced by utils.StatisticalTest::get_status().

towerMapEB_

std::vector<std::vector<std::pair<int, std::vector<EBDataFrame> > > > EcalEBTrigPrimTestAlgo::towerMapEB_

private

Definition at line 118 of file EcalEBTrigPrimTestAlgo.h.

Referenced by init(), and run().

towerMapEE_

std::vector<std::vector<std::pair<int, std::vector<EEDataFrame> > > > EcalEBTrigPrimTestAlgo::towerMapEE_

private

Definition at line 119 of file EcalEBTrigPrimTestAlgo.h.

towtp2_

std::vector<EcalEBTriggerPrimitiveSample> EcalEBTrigPrimTestAlgo::towtp2_

private

Definition at line 122 of file EcalEBTrigPrimTestAlgo.h.

towtp_

std::vector<EcalEBTriggerPrimitiveSample> EcalEBTrigPrimTestAlgo::towtp_

private

Definition at line 121 of file EcalEBTrigPrimTestAlgo.h.