HcalNoiseRBX.cc

Go to the documentation of this file.

//
// HcalNoiseRBX.cc
//
//   description: container class of RBX information for studying the HCAL Noise
//
//   author: J.P. Chou, Brown
//
//

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


using namespace reco;

// default constructor
HcalNoiseRBX::HcalNoiseRBX() :
  idnumber_(0), hpds_(4),
  allCharge_(HBHEDataFrame::MAXSAMPLES, 0.0)
{
}

// destructor
HcalNoiseRBX::~HcalNoiseRBX()
{
}

// accessors
int HcalNoiseRBX::idnumber(void) const
{
  return idnumber_;
}

const std::vector<HcalNoiseHPD> HcalNoiseRBX::HPDs(void) const
{
  return hpds_;
}

std::vector<HcalNoiseHPD>::const_iterator HcalNoiseRBX::maxHPD(double threshold) const
{
  std::vector<HcalNoiseHPD>::const_iterator maxit=hpds_.end();
  double maxenergy=-99999999.;
  for(std::vector<HcalNoiseHPD>::const_iterator it=hpds_.begin(); it!=hpds_.end(); ++it) {
    double tempenergy=it->recHitEnergy();
    if(tempenergy>maxenergy) {
      maxenergy=tempenergy;
      maxit=it;
    }
  }
  return maxit;
}

const std::vector<float> HcalNoiseRBX::allCharge(void) const
{
  return allCharge_;
}

float HcalNoiseRBX::allChargeTotal(void) const
{
  float total=0;
  for(unsigned int i=0; i<allCharge_.size(); i++)
    total += allCharge_[i];
  return total;
}
  
float HcalNoiseRBX::allChargeHighest2TS(unsigned int firstts) const
{
  float total=0;
  for(unsigned int i=firstts; i<firstts+2 && allCharge_.size(); i++)
    total += allCharge_[i];
  return total;
}
  
float HcalNoiseRBX::allChargeHighest3TS(unsigned int firstts) const
{
  float total=0;
  for(unsigned int i=firstts; i<firstts+3 && allCharge_.size(); i++)
    total += allCharge_[i];
  return total;
}
  

int HcalNoiseRBX::totalZeros(void) const
{
  int tot=0;
  for(unsigned int i=0; i<hpds_.size(); i++)
    tot += hpds_[i].totalZeros();
  return tot;
}

int HcalNoiseRBX::maxZeros(void) const
{
  int max=0;
  for(unsigned int i=0; i<hpds_.size(); i++)
    if(hpds_[i].maxZeros()>max)
      max=hpds_[i].maxZeros();
  return max;
}

double HcalNoiseRBX::recHitEnergy(double threshold) const
{
  double total=0;
  for(unsigned int i=0; i<hpds_.size(); i++)
    total += hpds_[i].recHitEnergy(threshold);
  return total;
}

double HcalNoiseRBX::minRecHitTime(double threshold) const
{
  double mintime=9999999.;
  for(unsigned int i=0; i<hpds_.size(); i++) {
    double temptime=hpds_[i].minRecHitTime(threshold);
    if(temptime<mintime) mintime=temptime;
  }
  return mintime;
}

double HcalNoiseRBX::maxRecHitTime(double threshold) const
{
  double maxtime=-9999999.;
  for(unsigned int i=0; i<hpds_.size(); i++) {
    double temptime=hpds_[i].maxRecHitTime(threshold);
    if(temptime>maxtime) maxtime=temptime;
  }
  return maxtime;
}

int HcalNoiseRBX::numRecHits(double threshold) const
{
  int total=0;
  for(unsigned int i=0; i<hpds_.size(); i++)
    total += hpds_[i].numRecHits(threshold);
  return total;
}

double HcalNoiseRBX::caloTowerHadE(void) const
{
  double h=0;
  towerset_t twrs;
  uniqueTowers(twrs);
  for(towerset_t::const_iterator it=twrs.begin(); it!=twrs.end(); ++it) {
    h += it->hadEnergy();
  }
  return h;
}

double HcalNoiseRBX::caloTowerEmE(void) const
{
  double e=0;
  towerset_t twrs;
  uniqueTowers(twrs);
  for(towerset_t::const_iterator it=twrs.begin(); it!=twrs.end(); ++it) {
    e += it->emEnergy();
  }
  return e;
}

double HcalNoiseRBX::caloTowerTotalE(void) const
{
  double e=0;
  towerset_t twrs;
  uniqueTowers(twrs);
  for(towerset_t::const_iterator it=twrs.begin(); it!=twrs.end(); ++it) {
    e += it->hadEnergy() + it->emEnergy();
  }
  return e;
}

double HcalNoiseRBX::caloTowerEmFraction(void) const
{
  double e=0, h=0;
  towerset_t twrs;
  uniqueTowers(twrs);
  for(towerset_t::const_iterator it=twrs.begin(); it!=twrs.end(); ++it) {
    h += it->hadEnergy();
    e += it->emEnergy();
  }
  return (e+h)==0 ? 999 : e/(e+h);
}

void HcalNoiseRBX::uniqueTowers(towerset_t& twrs_) const
{
  twrs_.clear();
  for(std::vector<HcalNoiseHPD>::const_iterator it1=hpds_.begin(); it1!=hpds_.end(); ++it1) {
    edm::RefVector<CaloTowerCollection> twrsref=it1->caloTowers();
    for(edm::RefVector<CaloTowerCollection>::const_iterator it2=twrsref.begin(); it2!=twrsref.end(); ++it2) {
      CaloTower twr=**it2;
      twrs_.insert(twr);
    }
  }
  return;
}