HcalNoiseHPD.cc

Go to the documentation of this file.

//
// HcalNoiseHPD.cc
//
//   description: container class of HPD information for analyzing HCAL Noise
//
//   author: J.P. Chou, Brown
//
//

#include "DataFormats/METReco/interface/HcalNoiseHPD.h"


using namespace reco;

// default constructor
HcalNoiseHPD::HcalNoiseHPD()
  : idnumber_(0), totalZeros_(0), maxZeros_(0),
    bigCharge_(HBHEDataFrame::MAXSAMPLES, 0.0),
    big5Charge_(HBHEDataFrame::MAXSAMPLES, 0.0)
{
  // reserve some space, so that there's no reallocation issues
  rechits_.reserve(19);
  calotowers_.reserve(19);
}
  
// destructor
HcalNoiseHPD::~HcalNoiseHPD()
{
}
  
// accessors
int HcalNoiseHPD::idnumber(void) const
{
  return idnumber_;
}
  
const std::vector<float> HcalNoiseHPD::bigCharge(void) const
{
  return bigCharge_;
}
  
float HcalNoiseHPD::bigChargeTotal(void) const
{
  float total=0;
  for(unsigned int i=0; i<bigCharge_.size(); i++) {
    total += bigCharge_[i];
  }
  return total;
}

float HcalNoiseHPD::bigChargeHighest2TS(unsigned int firstts) const
{
  float total=0;
  for(unsigned int i=firstts; i<firstts+2 && i<bigCharge_.size(); i++)
    total += bigCharge_[i];
  return total;
}


float HcalNoiseHPD::bigChargeHighest3TS(unsigned int firstts) const
{
  float total=0;
  for(unsigned int i=firstts; i<firstts+3 && i<bigCharge_.size(); i++)
    total += bigCharge_[i];
  return total;
}
  
const std::vector<float> HcalNoiseHPD::big5Charge(void) const
{
  return big5Charge_;
}
  
float HcalNoiseHPD::big5ChargeTotal(void) const
{
  float total=0;
  for(unsigned int i=0; i<big5Charge_.size(); i++) {
    total += big5Charge_[i];
  }
  return total;
}
  
float HcalNoiseHPD::big5ChargeHighest2TS(unsigned int firstts) const
{
  float total=0;
  for(unsigned int i=firstts; i<firstts+2 && i<big5Charge_.size(); i++)
    total += big5Charge_[i];
  return total;
}
  
float HcalNoiseHPD::big5ChargeHighest3TS(unsigned int firstts) const
{
  float total=0;
  for(unsigned int i=firstts; i<firstts+2 && i<big5Charge_.size(); i++)
    total += big5Charge_[i];
  return total;
}
  
int HcalNoiseHPD::totalZeros(void) const
{
  return totalZeros_;
}
  
int HcalNoiseHPD::maxZeros(void) const
{
  return maxZeros_;
}

const edm::RefVector<HBHERecHitCollection> HcalNoiseHPD::recHits(void) const
{
  return rechits_;
}
  
float HcalNoiseHPD::recHitEnergy(const float threshold) const
{
  double total=0.0;
  for(edm::RefVector<HBHERecHitCollection>::const_iterator it=rechits_.begin(); it!=rechits_.end(); ++it) {
    const float energy=(*it)->eraw();
    if(energy>=threshold) total+=energy;
  }
  return total;
}

float HcalNoiseHPD::recHitEnergyFailR45(const float threshold) const
{
  double total=0.0;
  for(edm::RefVector<HBHERecHitCollection>::const_iterator it=rechits_.begin(); it!=rechits_.end(); ++it) {
    const float energy=(*it)->eraw();
    if( (*it)->flagField(HcalCaloFlagLabels::HBHETS4TS5Noise) && 
       !(*it)->flagField(HcalCaloFlagLabels::HBHEOOTPU) ) 
      if(energy>=threshold) total+=energy;
  }
  return total;
}

float HcalNoiseHPD::minRecHitTime(const float threshold) const
{
  float mintime=9999999;
  for(edm::RefVector<HBHERecHitCollection>::const_iterator it=rechits_.begin(); it!=rechits_.end(); ++it) {
    if((*it)->energy()<threshold) continue;
    float time=(*it)->time();
    if(mintime>time) mintime=time;
  }
  return mintime;
}
  
float HcalNoiseHPD::maxRecHitTime(const float threshold) const
{
  float maxtime=-9999999;
  for(edm::RefVector<HBHERecHitCollection>::const_iterator it=rechits_.begin(); it!=rechits_.end(); ++it) {
    if((*it)->energy()<threshold) continue;
    float time=(*it)->time();
    if(maxtime<time) maxtime=time;
  }
  return maxtime;
}

int HcalNoiseHPD::numRecHits(const float threshold) const
{
  int count=0;
  for(edm::RefVector<HBHERecHitCollection>::const_iterator it=rechits_.begin(); it!=rechits_.end(); ++it) {
    // Exclude uncollapsed QIE11 channels
    if( CaloRecHitAuxSetter::getBit((*it)->auxPhase1(), HBHERecHitAuxSetter::OFF_TDC_TIME) &&
       !CaloRecHitAuxSetter::getBit((*it)->auxPhase1(), HBHERecHitAuxSetter::OFF_COMBINED) ) continue;
    if((*it)->eraw()>=threshold) ++count;
  }
  return count;
}

int HcalNoiseHPD::numRecHitsFailR45(const float threshold) const
{
  int count=0;
  for(edm::RefVector<HBHERecHitCollection>::const_iterator it=rechits_.begin(); it!=rechits_.end(); ++it)
    if( (*it)->flagField(HcalCaloFlagLabels::HBHETS4TS5Noise) && 
       !(*it)->flagField(HcalCaloFlagLabels::HBHEOOTPU)  )
      if((*it)->eraw()>=threshold) ++count;
  return count;
}

const edm::RefVector<CaloTowerCollection> HcalNoiseHPD::caloTowers(void) const
{
  return calotowers_;
}
  
double HcalNoiseHPD::caloTowerHadE(void) const
{
  double total=0;
  for(edm::RefVector<CaloTowerCollection>::const_iterator it=calotowers_.begin(); it!=calotowers_.end(); ++it)
    total += (*it)->hadEnergy();
  return total;
}
  
double HcalNoiseHPD::caloTowerEmE(void) const
{
  double total=0;
  for(edm::RefVector<CaloTowerCollection>::const_iterator it=calotowers_.begin(); it!=calotowers_.end(); ++it)
    total += (*it)->emEnergy();
  return total;
}

double HcalNoiseHPD::caloTowerTotalE(void) const
{
  double total=0;
  for(edm::RefVector<CaloTowerCollection>::const_iterator it=calotowers_.begin(); it!=calotowers_.end(); ++it)
    total += (*it)->emEnergy()+(*it)->hadEnergy();
  return total;
}
  
double HcalNoiseHPD::caloTowerEmFraction(void) const
{
  double h=0, e=0;
  for(edm::RefVector<CaloTowerCollection>::const_iterator it=calotowers_.begin(); it!=calotowers_.end(); ++it) {
    e += (*it)->emEnergy();
    h += (*it)->hadEnergy();
  }
  return (e+h)!=0 ? e/(e+h) : 999.;
}