#include <HcalNoiseAlgo.h>

Public Member Functions
	CommonHcalNoiseRBXData (const reco::HcalNoiseRBX &rbx, double minRecHitE, double minLowHitE, double minHighHitE)
double	e10ts (void) const
double	e2ts (void) const
double	energy (void) const
double	highEHitTimeSqrd (void) const
double	HPDEMF (void) const
double	lowEHitTimeSqrd (void) const
double	maxHighEHitTime (void) const
double	maxLowEHitTime (void) const
double	minHighEHitTime (void) const
double	minLowEHitTime (void) const
int	numHighEHits (void) const
int	numHPDHits (void) const
int	numHPDNoOtherHits (void) const
int	numLowEHits (void) const
int	numRBXHits (void) const
int	numZeros (void) const
double	ratio (void) const
double	RBXEMF (void) const
edm::RefVector < CaloTowerCollection >	rbxTowers (void) const
bool	validRatio (void) const
	~CommonHcalNoiseRBXData ()
Private Attributes
double	e10ts_
double	e2ts_
double	energy_
double	highEHitTimeSqrd_
double	HPDEMF_
double	lowEHitTimeSqrd_
double	maxHighEHitTime_
double	maxLowEHitTime_
double	minHighEHitTime_
double	minLowEHitTime_
int	numHighEHits_
int	numHPDHits_
int	numHPDNoOtherHits_
int	numLowEHits_
int	numRBXHits_
int	numZeros_
double	RBXEMF_
edm::RefVector < CaloTowerCollection >	rbxtowers_

Detailed Description

Definition at line 11 of file HcalNoiseAlgo.h.

Constructor & Destructor Documentation

CommonHcalNoiseRBXData::CommonHcalNoiseRBXData	(	const reco::HcalNoiseRBX &	rbx,
		double	minRecHitE,
		double	minLowHitE,
		double	minHighHitE
	)

Definition at line 3 of file HcalNoiseAlgo.cc.

{
  // energy
  energy_ = rbx.recHitEnergy(minRecHitE); 

  // ratio
  e2ts_ = rbx.allChargeHighest2TS();
  e10ts_ = rbx.allChargeTotal();

  // # of hits
  numHPDHits_ = 0;
  for(std::vector<reco::HcalNoiseHPD>::const_iterator it1=rbx.HPDsBegin(); it1!=rbx.HPDsEnd(); ++it1) {
    int nhpdhits=it1->numRecHits(minRecHitE);
    if(numHPDHits_ < nhpdhits) numHPDHits_ = nhpdhits;
  }
  numRBXHits_ = rbx.numRecHits(minRecHitE);
  numHPDNoOtherHits_ = (numHPDHits_ == numRBXHits_) ? numHPDHits_ : 0;    

  // # of ADC zeros
  numZeros_ = rbx.totalZeros();

  // timing
  minLowEHitTime_ = minHighEHitTime_ = 99999.;
  maxLowEHitTime_ = maxHighEHitTime_ = -99999.;
  lowEHitTimeSqrd_ = highEHitTimeSqrd_ = 0;
  numLowEHits_ = numHighEHits_ = 0;
  for(std::vector<reco::HcalNoiseHPD>::const_iterator it1=rbx.HPDsBegin(); it1!=rbx.HPDsEnd(); ++it1) {
    edm::RefVector<HBHERecHitCollection> rechits=it1->recHits();
    for(edm::RefVector<HBHERecHitCollection>::const_iterator it2=rechits.begin(); it2!=rechits.end(); ++it2) {
      float energy=(*it2)->energy();
      float time=(*it2)->time();
      if(energy>=minLowHitE) {
        if(minLowEHitTime_ > time) minLowEHitTime_ = time;
        if(maxLowEHitTime_ < time) maxLowEHitTime_ = time;
        lowEHitTimeSqrd_ += time*time;
        ++numLowEHits_;
      }
      if(energy>=minHighHitE) {
        if(minHighEHitTime_ > time) minHighEHitTime_ = time;
        if(maxHighEHitTime_ < time) maxHighEHitTime_ = time;
        highEHitTimeSqrd_ += time*time;
        ++numHighEHits_;
      }
    }
  }

  // emf
  HPDEMF_ = 999.;
  for(std::vector<reco::HcalNoiseHPD>::const_iterator it1=rbx.HPDsBegin(); it1!=rbx.HPDsEnd(); ++it1) {
    double eme=it1->caloTowerEmE();
    double hade=it1->recHitEnergy(minRecHitE);
    double emf=(eme+hade)==0 ? 999 : eme/(eme+hade);
    if(HPDEMF_ > emf) emf = HPDEMF_;
  }
  double eme=rbx.caloTowerEmE();
  RBXEMF_ = (eme+energy_)==0 ? 999 : eme/(eme+energy_);

  // calotowers
  rbxtowers_.clear();
  JoinCaloTowerRefVectorsWithoutDuplicates join;
  for(std::vector<reco::HcalNoiseHPD>::const_iterator it1=rbx.HPDsBegin(); it1!=rbx.HPDsEnd(); ++it1) {
    join(rbxtowers_, it1->caloTowers());
  }

  return;
}