ZeroSuppressFP420 Class Reference

#include <ZeroSuppressFP420.h>

Inheritance diagram for ZeroSuppressFP420:

Public Member Functions
void	initParams (const edm::ParameterSet &conf_)
ZSuppressFP420::DigitalMapType	trkFEDclusterizer (const DigitalMapType &, int)
ZSuppressFP420::DigitalMapType	zeroSuppress (const DigitalMapType &, int) override
	ZeroSuppressFP420 (const edm::ParameterSet &conf, float noise)
	~ZeroSuppressFP420 () override
Public Member Functions inherited from ZSuppressFP420
virtual	~ZSuppressFP420 ()

Private Attributes
int	algoConf
edm::ParameterSet	conf_
double	highthreshConf
double	lowthreshConf
float	noiseInAdc
short	theFEDalgorithm
float	theFEDhighThresh
float	theFEDlowThresh
short	theNumFEDalgos
int	verbosity

Additional Inherited Members
Public Types inherited from ZSuppressFP420
typedef DConverterFP420::DigitalMapType	DigitalMapType

Detailed Description

Definition at line 7 of file ZeroSuppressFP420.h.

Constructor & Destructor Documentation

ZeroSuppressFP420::ZeroSuppressFP420	(	const edm::ParameterSet &	conf,
		float	noise
	)

Definition at line 10 of file ZeroSuppressFP420.cc.

References conf_, gather_cfg::cout, edm::ParameterSet::getUntrackedParameter(), initParams(), noiseInAdc, and verbosity.

                                                                             : conf_(conf), theNumFEDalgos(4) {
  noiseInAdc = noise;
  initParams(conf_);
  verbosity = conf_.getUntrackedParameter<int>("VerbosityLevel");
  // initParams();
  if (verbosity > 0) {
    std::cout << "ZeroSuppressFP420: constructor: noiseInAdc=  " << noiseInAdc << std::endl;
  }
}

ZeroSuppressFP420::~ZeroSuppressFP420 ( )

inlineoverride

Definition at line 11 of file ZeroSuppressFP420.h.

References conf_, initParams(), trkFEDclusterizer(), and zeroSuppress().

{}

Member Function Documentation

void ZeroSuppressFP420::initParams

(

const edm::ParameterSet &

conf_

)

Definition at line 26 of file ZeroSuppressFP420.cc.

References algoConf, gather_cfg::cout, cmsRelvalreport::exit, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), highthreshConf, lowthreshConf, noiseInAdc, theFEDalgorithm, theFEDhighThresh, theFEDlowThresh, theNumFEDalgos, and verbosity.

Referenced by ZeroSuppressFP420(), and ~ZeroSuppressFP420().

                                                               {
  verbosity = conf_.getUntrackedParameter<int>("VerbosityLevel");
  algoConf = conf_.getParameter<int>("FedFP420Algorithm");               // FedFP420Algorithm: =1  (,2,3,4)
  lowthreshConf = conf_.getParameter<double>("FedFP420LowThreshold");    // FedFP420LowThreshold  =3.
  highthreshConf = conf_.getParameter<double>("FedFP420HighThreshold");  // FedFP420HighThreshold  =4.

  /*
   * There are four possible algorithms, the default of which (4)
   * has different thresholds for isolated channels and ones in clusters.
   * It also merges clusters (single or multi channels) that are only separated
   * by one hole. This channel is selected as signal even though it is below
   * both thresholds.
   */

  theFEDalgorithm = algoConf;
  theFEDlowThresh = lowthreshConf * noiseInAdc;
  theFEDhighThresh = highthreshConf * noiseInAdc;

  if (verbosity > 0) {
    std::cout << "ZeroSuppressFP420: initParams: !!!  theFEDalgorithm=  " << theFEDalgorithm << std::endl;
    std::cout << " lowthreshConf=  " << lowthreshConf << " highthreshConf=  " << highthreshConf
              << " theFEDlowThresh=  " << theFEDlowThresh << " theFEDhighThresh=  " << theFEDhighThresh << std::endl;
  }

  // Check zero suppress algorithm
  if (theFEDalgorithm < 1 || theFEDalgorithm > theNumFEDalgos) {
    edm::LogError("FP420DigiInfo") << "ZeroSuppressFP420 FATAL ERROR: Unknown zero suppress algorithm "
                                   << theFEDalgorithm;
    exit(1);
  }

  // Check thresholds
  if (theFEDlowThresh > theFEDhighThresh) {
    edm::LogError("FP420DigiInfo") << "ZeroSuppressFP420 FATAL ERROR: Low threshold exceeds high "
                                      "threshold: "
                                   << theFEDlowThresh << " > " << theFEDhighThresh;
    exit(2);
  }
}

ZSuppressFP420::DigitalMapType ZeroSuppressFP420::trkFEDclusterizer	(	const DigitalMapType &	in,
		int	vrb
	)

Definition at line 74 of file ZeroSuppressFP420.cc.

References accept(), ecalMGPA::adc(), gather_cfg::cout, mps_fire::i, SiStripPI::max, indexGen::s2, AlCaHLTBitMon_QueryRunRegistry::string, digitizers_cfi::strip, theFEDalgorithm, theFEDhighThresh, and theFEDlowThresh.

Referenced by zeroSuppress(), and ~ZeroSuppressFP420().

                                                                                                   {
  const std::string s2("ZeroSuppressFP420::trkFEDclusterizer1");

  DigitalMapType selectedSignal;
  DigitalMapType::const_iterator i, iPrev, iNext, iPrev2, iNext2;

  if (vrb > 0) {
    std::cout << "Before For loop" << std::endl;
  }

  for (i = in.begin(); i != in.end(); i++) {
    // Find adc values for neighbouring strips
    int strip = i->first;
    int adc = i->second;
    iPrev = in.find(strip - 1);
    iNext = in.find(strip + 1);
    if (vrb > 0) {
      std::cout << "Inside For loop trkFEDclusterizer: strip= " << strip << " adc= " << adc << std::endl;
    }
    // Set values for channels just outside module to infinity.
    // This is to avoid losing channels at the edges,
    // which otherwise would pass cuts if strips were next to each other.
    int adcPrev = -99999;
    int adcNext = -99999;
    if (((strip) % 128) == 127)
      adcNext = 99999;
    if (((strip) % 128) == 0)
      adcPrev = 99999;
    // Otherwise if channel was found then find it's ADC count.
    if (iPrev != in.end())
      adcPrev = iPrev->second;
    if (iNext != in.end())
      adcNext = iNext->second;
    int adcMaxNeigh = std::max(adcPrev, adcNext);
    if (vrb > 0) {
      std::cout << "adcPrev= " << adcPrev << " adcNext= " << adcNext << " adcMaxNeigh= " << adcMaxNeigh << std::endl;
    }

    // Find adc values for next neighbouring channes
    iPrev2 = in.find(strip - 2);
    iNext2 = in.find(strip + 2);
    // See above
    int adcPrev2 = -99999;
    int adcNext2 = -99999;
    if (((strip) % 128) == 126)
      adcNext2 = 99999;
    if (((strip) % 128) == 1)
      adcPrev2 = 99999;
    if (iPrev2 != in.end())
      adcPrev2 = iPrev2->second;
    if (iNext2 != in.end())
      adcNext2 = iNext2->second;

    if (vrb > 0) {
      std::cout << "adcPrev2= " << adcPrev2 << " adcNext2= " << adcNext2 << std::endl;
      std::cout << "To be accepted or not?  adc= " << adc << " >= theFEDlowThresh=" << theFEDlowThresh << std::endl;
    }
    // Decide if this channel should be accepted.
    bool accept = false;
    switch (theFEDalgorithm) {
      case 1:
        accept = (adc >= theFEDlowThresh);
        break;

      case 2:
        accept = (adc >= theFEDhighThresh || (adc >= theFEDlowThresh && adcMaxNeigh >= theFEDlowThresh));
        break;

      case 3:
        accept = (adc >= theFEDhighThresh || (adc >= theFEDlowThresh && adcMaxNeigh >= theFEDhighThresh));
        break;

      case 4:
        accept = ((adc >= theFEDhighThresh) ||  // Test for adc>highThresh (same as
                                                // algorithm 2)
                  ((adc >= theFEDlowThresh) &&  // Test for adc>lowThresh, with neighbour
                                                // adc>lowThresh (same as algorithm 2)
                   (adcMaxNeigh >= theFEDlowThresh)) ||
                  ((adc < theFEDlowThresh) &&          // Test for adc<lowThresh
                   (((adcPrev >= theFEDhighThresh) &&  // with both neighbours>highThresh
                     (adcNext >= theFEDhighThresh)) ||
                    ((adcPrev >= theFEDhighThresh) &&  // OR with previous neighbour>highThresh and
                     (adcNext >= theFEDlowThresh) &&   // both the next neighbours>lowThresh
                     (adcNext2 >= theFEDlowThresh)) ||
                    ((adcNext >= theFEDhighThresh) &&  // OR with next neighbour>highThresh and
                     (adcPrev >= theFEDlowThresh) &&   // both the previous neighbours>lowThresh
                     (adcPrev2 >= theFEDlowThresh)) ||
                    ((adcNext >= theFEDlowThresh) &&   // OR with both next neighbours>lowThresh and
                     (adcNext2 >= theFEDlowThresh) &&  // both the previous neighbours>lowThresh
                     (adcPrev >= theFEDlowThresh) && (adcPrev2 >= theFEDlowThresh)))));
        break;
    }

    /*
     * When a channel satisfying only the lower threshold is at the edge of an
     * APV or module, the trkFEDclusterizer method assumes that every channel
     * just outside an APV or module has a hit on it. This is to avoid channel
     * inefficiencies at the edges of APVs and modules.
     */
    if (accept) {
      selectedSignal[strip] = adc;

      if (vrb > 0) {
        std::cout << "selected strips = " << strip << " adc= " << adc << std::endl;
      }
    }
  }

  if (vrb > 0) {
    std::cout << "last line of trkFEDclusterizer: return selectedSignal" << std::endl;
  }
  return selectedSignal;
}

ZSuppressFP420::DigitalMapType ZeroSuppressFP420::zeroSuppress ( const DigitalMapType &  notZeroSuppressedMap, int  vrb  )

overridevirtual

Implements ZSuppressFP420.

Definition at line 66 of file ZeroSuppressFP420.cc.

References gather_cfg::cout, and trkFEDclusterizer().

Referenced by FP420DigiMain::run(), and ~ZeroSuppressFP420().

                                                                                                                {
  return trkFEDclusterizer(notZeroSuppressedMap, vrb);
  if (vrb > 0) {
    std::cout << "zeroSuppress: return trkFEDclusterizer(notZeroSuppressedMap)" << std::endl;
  }
}

Member Data Documentation

int ZeroSuppressFP420::algoConf

private

Definition at line 28 of file ZeroSuppressFP420.h.

Referenced by initParams().

edm::ParameterSet ZeroSuppressFP420::conf_

private

Definition at line 25 of file ZeroSuppressFP420.h.

Referenced by ZeroSuppressFP420(), and ~ZeroSuppressFP420().

double ZeroSuppressFP420::highthreshConf

private

Definition at line 30 of file ZeroSuppressFP420.h.

Referenced by initParams().

double ZeroSuppressFP420::lowthreshConf

private

Definition at line 29 of file ZeroSuppressFP420.h.

Referenced by initParams().

float ZeroSuppressFP420::noiseInAdc

private

Definition at line 21 of file ZeroSuppressFP420.h.

Referenced by initParams(), and ZeroSuppressFP420().

short ZeroSuppressFP420::theFEDalgorithm

private

Definition at line 22 of file ZeroSuppressFP420.h.

Referenced by initParams(), and trkFEDclusterizer().

float ZeroSuppressFP420::theFEDhighThresh

private

Definition at line 24 of file ZeroSuppressFP420.h.

Referenced by initParams(), and trkFEDclusterizer().

float ZeroSuppressFP420::theFEDlowThresh

private

Definition at line 23 of file ZeroSuppressFP420.h.

Referenced by initParams(), and trkFEDclusterizer().

short ZeroSuppressFP420::theNumFEDalgos

private

Definition at line 26 of file ZeroSuppressFP420.h.

Referenced by initParams().

int ZeroSuppressFP420::verbosity

private

Definition at line 28 of file ZeroSuppressFP420.h.

Referenced by initParams(), and ZeroSuppressFP420().