HBHEStatusBitSetter Class Reference

#include <HBHEStatusBitSetter.h>

Public Member Functions
void	Clear ()
	HBHEStatusBitSetter ()
	HBHEStatusBitSetter (double nominalPedestal, double hitEnergyMinimum, int hitMultiplicityThreshold, const std::vector< edm::ParameterSet > &pulseShapeParameterSets)
	HBHEStatusBitSetter (const HBHEStatusBitSetter &)=delete
HBHEStatusBitSetter &	operator= (const HBHEStatusBitSetter &)=delete
void	rememberHit (const HBHERecHit &hbhe)
void	SetFlagsFromDigi (HBHERecHit &hbhe, const HBHEDataFrame &digi, const HcalCoder &coder, const HcalCalibrations &calib)
void	SetFlagsFromRecHits (HBHERecHitCollection &rec)
void	SetFrontEndMap (const HcalFrontEndMap *m)
	~HBHEStatusBitSetter ()

Private Attributes
const HcalFrontEndMap *	frontEndMap_
double	hitEnergyMinimum_
int	hitMultiplicityThreshold_
std::vector< int >	hpdMultiplicity_
double	nominalPedestal_
std::vector< std::vector< double > >	pulseShapeParameters_

Detailed Description

Definition at line 14 of file HBHEStatusBitSetter.h.

Constructor & Destructor Documentation

HBHEStatusBitSetter() [1/3]

HBHEStatusBitSetter::HBHEStatusBitSetter

(

)

Definition at line 5 of file HBHEStatusBitSetter.cc.

References hitEnergyMinimum_, hitMultiplicityThreshold_, and nominalPedestal_.

                                         : frontEndMap_(nullptr) {
  nominalPedestal_ = 3.0;
  hitEnergyMinimum_ = 2.0;
  hitMultiplicityThreshold_ = 17;
}

HBHEStatusBitSetter() [2/3]

HBHEStatusBitSetter::HBHEStatusBitSetter	(	double	nominalPedestal,
		double	hitEnergyMinimum,
		int	hitMultiplicityThreshold,
		const std::vector< edm::ParameterSet > &	pulseShapeParameterSets
	)

Definition at line 11 of file HBHEStatusBitSetter.cc.

References submitPVValidationJobs::params, muonDTDigis_cfi::pset, pulseShapeParameters_, and hltHbhereco_cfi::pulseShapeParameterSets.

    : hitEnergyMinimum_(hitEnergyMinimum),
      hitMultiplicityThreshold_(hitMultiplicityThreshold),
      nominalPedestal_(nominalPedestal),
      frontEndMap_(nullptr) {
  const unsigned sz = pulseShapeParameterSets.size();
  pulseShapeParameters_.reserve(sz);
  for (unsigned iPSet = 0; iPSet < sz; iPSet++) {
    const edm::ParameterSet& pset = pulseShapeParameterSets[iPSet];
    const std::vector<double>& params = pset.getParameter<std::vector<double> >("pulseShapeParameters");
    pulseShapeParameters_.push_back(params);
  }
}

~HBHEStatusBitSetter()

HBHEStatusBitSetter::~HBHEStatusBitSetter

(

)

Definition at line 28 of file HBHEStatusBitSetter.cc.

{}

HBHEStatusBitSetter() [3/3]

HBHEStatusBitSetter::HBHEStatusBitSetter ( const HBHEStatusBitSetter &  )

delete

Member Function Documentation

Clear()

void HBHEStatusBitSetter::Clear

(

)

Definition at line 42 of file HBHEStatusBitSetter.cc.

References hpdMultiplicity_, and mps_fire::i.

                                {
  const unsigned sz = hpdMultiplicity_.size();
  for (unsigned i = 0; i < sz; i++)
    hpdMultiplicity_[i] = 0;
}

operator=()

HBHEStatusBitSetter& HBHEStatusBitSetter::operator= ( const HBHEStatusBitSetter &  )

delete

rememberHit()

void HBHEStatusBitSetter::rememberHit

(

const HBHERecHit &

hbhe

)

Definition at line 48 of file HBHEStatusBitSetter.cc.

References frontEndMap_, photonIsolationHIProducer_cfi::hbhe, hitEnergyMinimum_, hpdMultiplicity_, and HcalFrontEndMap::lookupRMIndex().

Referenced by SetFlagsFromDigi().

                                                            {
  if (frontEndMap_ == nullptr) {
    edm::LogError("HBHEStatusBitSetter") << "No HcalFrontEndMap in rememberHit";
    return;
  }
  //increment hit multiplicity
  if (hbhe.energy() > hitEnergyMinimum_) {
    const int index = frontEndMap_->lookupRMIndex(hbhe.detid());
    hpdMultiplicity_.at(index)++;
  }
}

SetFlagsFromDigi()

void HBHEStatusBitSetter::SetFlagsFromDigi	(	HBHERecHit &	hbhe,
		const HBHEDataFrame &	digi,
		const HcalCoder &	coder,
		const HcalCalibrations &	calib
	)

Definition at line 60 of file HBHEStatusBitSetter.cc.

References HcalCoder::adc2fC(), frontEndMap_, photonIsolationHIProducer_cfi::hbhe, HcalCaloFlagLabels::HBHEPulseShape, nominalPedestal_, pulseShapeParameters_, rememberHit(), HBHEDataFrame::size(), and findQualityFiles::size.

                                                                          {
  if (frontEndMap_ == nullptr) {
    edm::LogError("HBHEStatusBitSetter") << "No HcalFrontEndMap in SetFlagsFromDigi";
    return;
  }
  rememberHit(hbhe);

  //set pulse shape bits
  // Shuichi's algorithm uses the "correct" charge & pedestals, while Ted's uses "nominal" values.
  // Perhaps we should correct Ted's algorithm in the future, though that will mean re-tuning thresholds for his cuts. -- Jeff, 28 May 2010
  //double shuichi_charge_total=0.0;
  double nominal_charge_total = 0.0;
  double charge_max3 = -100.0;
  double charge_late3 = -100.0;
  unsigned int slice_max3 = 0;
  unsigned int size = digi.size();

  CaloSamples tool;
  coder.adc2fC(digi, tool);

  //  int capid=-1;
  for (unsigned int iSlice = 0; iSlice < size; iSlice++) {
    //      capid  = digi.sample(iSlice).capid();
    //shuichi_charge_total+=tool[iSlice]-calib.pedestal(capid);
    nominal_charge_total += digi[iSlice].nominal_fC() - nominalPedestal_;

    if (iSlice < 2)
      continue;
    // digi[i].nominal_fC() could be replaced by tool[iSlice], I think...  -- Jeff, 11 April 2011
    double qsum3 = digi[iSlice].nominal_fC() + digi[iSlice - 1].nominal_fC() + digi[iSlice - 2].nominal_fC() -
                   3 * nominalPedestal_;
    if (qsum3 > charge_max3) {
      charge_max3 = qsum3;
      slice_max3 = iSlice;
    }
  }

  if ((4 + slice_max3) > size)
    return;
  charge_late3 = digi[slice_max3 + 1].nominal_fC() + digi[slice_max3 + 2].nominal_fC() +
                 digi[slice_max3 + 3].nominal_fC() - 3 * nominalPedestal_;

  for (unsigned int iCut = 0; iCut < pulseShapeParameters_.size(); iCut++) {
    if (pulseShapeParameters_[iCut].size() != 6)
      continue;
    if (nominal_charge_total < pulseShapeParameters_[iCut].at(0) ||
        nominal_charge_total >= pulseShapeParameters_[iCut].at(1))
      continue;
    if (charge_late3 < (pulseShapeParameters_[iCut].at(2) + nominal_charge_total * pulseShapeParameters_[iCut].at(3)))
      continue;
    if (charge_late3 >= (pulseShapeParameters_[iCut].at(4) + nominal_charge_total * pulseShapeParameters_[iCut].at(5)))
      continue;
    hbhe.setFlagField(1, HcalCaloFlagLabels::HBHEPulseShape);
    return;
  }
}

SetFlagsFromRecHits()

void HBHEStatusBitSetter::SetFlagsFromRecHits

(

HBHERecHitCollection &

rec

)

Definition at line 120 of file HBHEStatusBitSetter.cc.

References edm::SortedCollection< T, SORT >::begin(), edm::SortedCollection< T, SORT >::end(), frontEndMap_, HcalCaloFlagLabels::HBHEHpdHitMultiplicity, hitMultiplicityThreshold_, hpdMultiplicity_, and HcalFrontEndMap::lookupRMIndex().

                                                                       {
  if (frontEndMap_ == nullptr) {
    edm::LogError("HBHEStatusBitSetter") << "No HcalFrontEndMap in SetFlagsFromRecHits";
    return;
  }

  for (HBHERecHitCollection::iterator iHBHE = rec.begin(); iHBHE != rec.end(); ++iHBHE) {
    const int index = frontEndMap_->lookupRMIndex(iHBHE->detid());
    if (hpdMultiplicity_.at(index) < hitMultiplicityThreshold_)
      continue;
    iHBHE->setFlagField(1, HcalCaloFlagLabels::HBHEHpdHitMultiplicity);
  }
}

SetFrontEndMap()

void HBHEStatusBitSetter::SetFrontEndMap

(

const HcalFrontEndMap *

m

)

Definition at line 30 of file HBHEStatusBitSetter.cc.

References frontEndMap_, hpdMultiplicity_, visualization-live-secondInstance_cfg::m, and HcalFrontEndMap::maxRMIndex().

                                                                 {
  frontEndMap_ = m;
  hpdMultiplicity_.clear();
  if (frontEndMap_) {
    const int sz = frontEndMap_->maxRMIndex();
    hpdMultiplicity_.reserve(sz);
    for (int iRm = 0; iRm < sz; iRm++) {
      hpdMultiplicity_.push_back(0);
    }
  }
}

Member Data Documentation

frontEndMap_

const HcalFrontEndMap* HBHEStatusBitSetter::frontEndMap_

private

Definition at line 39 of file HBHEStatusBitSetter.h.

Referenced by rememberHit(), SetFlagsFromDigi(), SetFlagsFromRecHits(), and SetFrontEndMap().

hitEnergyMinimum_

double HBHEStatusBitSetter::hitEnergyMinimum_

private

Definition at line 36 of file HBHEStatusBitSetter.h.

Referenced by HBHEStatusBitSetter(), and rememberHit().

hitMultiplicityThreshold_

int HBHEStatusBitSetter::hitMultiplicityThreshold_

private

Definition at line 37 of file HBHEStatusBitSetter.h.

Referenced by HBHEStatusBitSetter(), and SetFlagsFromRecHits().

hpdMultiplicity_

std::vector<int> HBHEStatusBitSetter::hpdMultiplicity_

private

Definition at line 40 of file HBHEStatusBitSetter.h.

Referenced by Clear(), rememberHit(), SetFlagsFromRecHits(), and SetFrontEndMap().

nominalPedestal_

double HBHEStatusBitSetter::nominalPedestal_

private

Definition at line 38 of file HBHEStatusBitSetter.h.

Referenced by HBHEStatusBitSetter(), and SetFlagsFromDigi().

pulseShapeParameters_

std::vector<std::vector<double> > HBHEStatusBitSetter::pulseShapeParameters_

private

Definition at line 41 of file HBHEStatusBitSetter.h.

Referenced by HBHEStatusBitSetter(), and SetFlagsFromDigi().