#include <HBHEStatusBitSetter.h>

Public Member Functions
void	Clear ()

	HBHEStatusBitSetter ()

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

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 Member Functions
	HBHEStatusBitSetter (const HBHEStatusBitSetter &)

HBHEStatusBitSetter &	operator= (const 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 15 of file HBHEStatusBitSetter.h.

Constructor & Destructor Documentation

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::HBHEStatusBitSetter	(	double	nominalPedestal,
		double	hitEnergyMinimum,
		int	hitMultiplicityThreshold,
		const std::vector< edm::ParameterSet > &	pulseShapeParameterSets
	)

Definition at line 11 of file HBHEStatusBitSetter.cc.

References edm::ParameterSet::getParameter(), muonDTDigis_cfi::pset, and pulseShapeParameters_.

   : 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 ( )

Definition at line 28 of file HBHEStatusBitSetter.cc.

28 { }

HBHEStatusBitSetter::HBHEStatusBitSetter ( const HBHEStatusBitSetter & )

private

Member Function Documentation

void HBHEStatusBitSetter::Clear ( )

Definition at line 42 of file HBHEStatusBitSetter.cc.

References hpdMultiplicity_, and mps_fire::i.

Referenced by HcalHitReconstructor::produce().

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

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

private

void HBHEStatusBitSetter::rememberHit ( const HBHERecHit & hbhe )

Definition at line 48 of file HBHEStatusBitSetter.cc.

References CaloRecHit::detid(), CaloRecHit::energy(), frontEndMap_, hitEnergyMinimum_, hpdMultiplicity_, diffTreeTool::index, 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)++;
   }
 }

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_, HcalCaloFlagLabels::HBHEPulseShape, nominalPedestal_, pulseShapeParameters_, rememberHit(), CaloRecHit::setFlagField(), HBHEDataFrame::size(), and findQualityFiles::size.

Referenced by HcalHitReconstructor::produce().

                                                       {
   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;
   }
   
 }

void HBHEStatusBitSetter::SetFlagsFromRecHits ( HBHERecHitCollection & rec )

Definition at line 113 of file HBHEStatusBitSetter.cc.

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

Referenced by HcalHitReconstructor::produce().

                                                                        {
 
   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);
   }
 }

void HBHEStatusBitSetter::SetFrontEndMap ( const HcalFrontEndMap * m )

Definition at line 30 of file HBHEStatusBitSetter.cc.

References frontEndMap_, hpdMultiplicity_, funct::m, and HcalFrontEndMap::maxRMIndex().

Referenced by HcalHitReconstructor::beginRun().

                                                                  {
   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);
     }
   }
 }