CommonHcalNoiseRBXData Class Reference

#include <HcalNoiseAlgo.h>

Public Member Functions
bool	CheckPassFilter (double Charge, double Discriminant, std::vector< std::pair< double, double > > const &Cuts, int Side)
	CommonHcalNoiseRBXData (const reco::HcalNoiseRBX &rbx, double minRecHitE, double minLowHitE, double minHighHitE, double TS4TS5EnergyThreshold, std::vector< std::pair< double, double > > const &TS4TS5UpperCut, std::vector< std::pair< double, double > > const &TS4TS5LowerCut, double MinRBXRechitR45E)
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
bool	PassTS4TS5 (void) const
int	r45Count (void) const
double	r45EnergyFraction (void) const
double	r45Fraction (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_
int	r45Count_
double	r45EnergyFraction_
double	r45Fraction_
double	RBXEMF_
edm::RefVector< CaloTowerCollection >	rbxtowers_
bool	TS4TS5Decision_

Detailed Description

Definition at line 11 of file HcalNoiseAlgo.h.

Constructor & Destructor Documentation

CommonHcalNoiseRBXData()

CommonHcalNoiseRBXData::CommonHcalNoiseRBXData	(	const reco::HcalNoiseRBX &	rbx,
		double	minRecHitE,
		double	minLowHitE,
		double	minHighHitE,
		double	TS4TS5EnergyThreshold,
		std::vector< std::pair< double, double > > const &	TS4TS5UpperCut,
		std::vector< std::pair< double, double > > const &	TS4TS5LowerCut,
		double	MinRBXRechitR45E
	)

Definition at line 3 of file HcalNoiseAlgo.cc.

References reco::HcalNoiseRBX::allCharge(), reco::HcalNoiseRBX::allChargeHighest2TS(), reco::HcalNoiseRBX::allChargeTotal(), reco::HcalNoiseRBX::caloTowerEmE(), CheckPassFilter(), e10ts_, e2ts_, custom_jme_cff::emf, energy(), energy_, highEHitTimeSqrd_, HPDEMF_, reco::HcalNoiseRBX::HPDsBegin(), reco::HcalNoiseRBX::HPDsEnd(), join(), lowEHitTimeSqrd_, maxHighEHitTime_, maxLowEHitTime_, minHighEHitTime_, minLowEHitTime_, numHighEHits_, numHPDHits_, numHPDNoOtherHits_, numLowEHits_, numRBXHits_, reco::HcalNoiseRBX::numRecHits(), reco::HcalNoiseRBX::numRecHitsFailR45(), numZeros_, r45Count_, r45EnergyFraction_, r45Fraction_, RBXEMF_, rbxtowers_, reco::HcalNoiseRBX::recHitEnergy(), reco::HcalNoiseRBX::recHitEnergyFailR45(), HI_PhotonSkim_cff::rechits, protons_cff::time, reco::HcalNoiseRBX::totalZeros(), TS4TS5Decision_, HLT_2022v15_cff::TS4TS5LowerCut, and HLT_2022v15_cff::TS4TS5UpperCut.

    : r45Count_(0), r45Fraction_(0), r45EnergyFraction_(0) {
  // energy
  energy_ = rbx.recHitEnergy(minRecHitE);

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

  // TS4TS5
  TS4TS5Decision_ = true;
  if (energy_ > TS4TS5EnergyThreshold)  // check filter
  {
    std::vector<float> AllCharge = rbx.allCharge();
    double BaseCharge = AllCharge[4] + AllCharge[5];
    if (BaseCharge < 1)
      BaseCharge = 1;
    double TS4TS5 = (AllCharge[4] - AllCharge[5]) / BaseCharge;

    if (CheckPassFilter(BaseCharge, TS4TS5, TS4TS5UpperCut, 1) == false)
      TS4TS5Decision_ = false;
    if (CheckPassFilter(BaseCharge, TS4TS5, TS4TS5LowerCut, -1) == false)
      TS4TS5Decision_ = false;
  } else
    TS4TS5Decision_ = true;

  // Rechit-wide R45
  int rbxHitCount = rbx.numRecHits(minRBXRechitR45E);
  r45Count_ = 0;
  r45Fraction_ = 0;
  r45EnergyFraction_ = 0;
  if (rbxHitCount > 0) {
    r45Count_ = rbx.numRecHitsFailR45(minRBXRechitR45E);
    r45Fraction_ = (double)(r45Count_) / (double)(rbxHitCount);
    r45EnergyFraction_ = rbx.recHitEnergyFailR45(minRBXRechitR45E) / rbx.recHitEnergy(minRBXRechitR45E);
  }

  // # 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 (get the one with minimum value)
  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)
      HPDEMF_ = emf;
  }
  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;
}

~CommonHcalNoiseRBXData()

CommonHcalNoiseRBXData::~CommonHcalNoiseRBXData ( )

inline

Definition at line 21 of file HcalNoiseAlgo.h.

{}

Member Function Documentation

CheckPassFilter()

bool CommonHcalNoiseRBXData::CheckPassFilter	(	double	Charge,
		double	Discriminant,
		std::vector< std::pair< double, double > > const &	Cuts,
		int	Side
	)

Definition at line 383 of file HcalNoiseAlgo.cc.

References PixelTestBeamValidation_cfi::Charge, first, mps_fire::i, createfilelist::int, PixelRegions::L1, PixelRegions::L2, and remoteMonitoring_LASER_era2018_cfg::limit.

Referenced by CommonHcalNoiseRBXData().

                                                       {
  //
  // Checks whether Discriminant value passes Cuts for the specified Charge.  True if pulse is good.
  //
  // The "Cuts" pairs are assumed to be sorted in terms of size from small to large,
  //    where each "pair" = (Charge, Discriminant)
  // "Side" is either positive or negative, which determines whether to discard the pulse if discriminant
  //    is greater or smaller than the cut value
  //

  if (Cuts.empty())  // safety check that there are some cuts defined
    return true;

  if (Charge <= Cuts[0].first)  // too small to cut on
    return true;

  int IndexLargerThanCharge = -1;  // find the range it is falling in
  for (int i = 1; i < (int)Cuts.size(); i++) {
    if (Cuts[i].first > Charge) {
      IndexLargerThanCharge = i;
      break;
    }
  }

  double limit = 1000000;

  if (IndexLargerThanCharge == -1)  // if charge is greater than the last entry, assume flat line
    limit = Cuts[Cuts.size() - 1].second;
  else  // otherwise, do a linear interpolation to find the cut position
  {
    double C1 = Cuts[IndexLargerThanCharge].first;
    double C2 = Cuts[IndexLargerThanCharge - 1].first;
    double L1 = Cuts[IndexLargerThanCharge].second;
    double L2 = Cuts[IndexLargerThanCharge - 1].second;

    limit = (Charge - C1) / (C2 - C1) * (L2 - L1) + L1;
  }

  if (Side > 0 && Discriminant > limit)
    return false;
  if (Side < 0 && Discriminant < limit)
    return false;

  return true;
}

e10ts()

double CommonHcalNoiseRBXData::e10ts ( void  ) const

inline

Definition at line 27 of file HcalNoiseAlgo.h.

References e10ts_.

{ return e10ts_; }

e2ts()

double CommonHcalNoiseRBXData::e2ts ( void  ) const

inline

Definition at line 26 of file HcalNoiseAlgo.h.

References e2ts_.

{ return e2ts_; }

energy()

double CommonHcalNoiseRBXData::energy ( void  ) const

inline

Definition at line 24 of file HcalNoiseAlgo.h.

References energy_.

Referenced by CommonHcalNoiseRBXData(), and Jet.Jet::rawEnergy().

{ return energy_; }

highEHitTimeSqrd()

double CommonHcalNoiseRBXData::highEHitTimeSqrd ( void  ) const

inline

Definition at line 39 of file HcalNoiseAlgo.h.

References highEHitTimeSqrd_.

{ return highEHitTimeSqrd_; }

HPDEMF()

double CommonHcalNoiseRBXData::HPDEMF ( void  ) const

inline

Definition at line 42 of file HcalNoiseAlgo.h.

References HPDEMF_.

{ return HPDEMF_; }

lowEHitTimeSqrd()

double CommonHcalNoiseRBXData::lowEHitTimeSqrd ( void  ) const

inline

Definition at line 35 of file HcalNoiseAlgo.h.

References lowEHitTimeSqrd_.

{ return lowEHitTimeSqrd_; }

maxHighEHitTime()

double CommonHcalNoiseRBXData::maxHighEHitTime ( void  ) const

inline

Definition at line 38 of file HcalNoiseAlgo.h.

References maxHighEHitTime_.

{ return maxHighEHitTime_; }

maxLowEHitTime()

double CommonHcalNoiseRBXData::maxLowEHitTime ( void  ) const

inline

Definition at line 34 of file HcalNoiseAlgo.h.

References maxLowEHitTime_.

{ return maxLowEHitTime_; }

minHighEHitTime()

double CommonHcalNoiseRBXData::minHighEHitTime ( void  ) const

inline

Definition at line 37 of file HcalNoiseAlgo.h.

References minHighEHitTime_.

{ return minHighEHitTime_; }

minLowEHitTime()

double CommonHcalNoiseRBXData::minLowEHitTime ( void  ) const

inline

Definition at line 33 of file HcalNoiseAlgo.h.

References minLowEHitTime_.

{ return minLowEHitTime_; }

numHighEHits()

int CommonHcalNoiseRBXData::numHighEHits ( void  ) const

inline

Definition at line 40 of file HcalNoiseAlgo.h.

References numHighEHits_.

{ return numHighEHits_; }

numHPDHits()

int CommonHcalNoiseRBXData::numHPDHits ( void  ) const

inline

Definition at line 29 of file HcalNoiseAlgo.h.

References numHPDHits_.

{ return numHPDHits_; }

numHPDNoOtherHits()

int CommonHcalNoiseRBXData::numHPDNoOtherHits ( void  ) const

inline

Definition at line 31 of file HcalNoiseAlgo.h.

References numHPDNoOtherHits_.

{ return numHPDNoOtherHits_; }

numLowEHits()

int CommonHcalNoiseRBXData::numLowEHits ( void  ) const

inline

Definition at line 36 of file HcalNoiseAlgo.h.

References numLowEHits_.

{ return numLowEHits_; }

numRBXHits()

int CommonHcalNoiseRBXData::numRBXHits ( void  ) const

inline

Definition at line 30 of file HcalNoiseAlgo.h.

References numRBXHits_.

{ return numRBXHits_; }

numZeros()

int CommonHcalNoiseRBXData::numZeros ( void  ) const

inline

Definition at line 32 of file HcalNoiseAlgo.h.

References numZeros_.

{ return numZeros_; }

PassTS4TS5()

bool CommonHcalNoiseRBXData::PassTS4TS5 ( void  ) const

inline

Definition at line 43 of file HcalNoiseAlgo.h.

References TS4TS5Decision_.

{ return TS4TS5Decision_; }

r45Count()

int CommonHcalNoiseRBXData::r45Count ( void  ) const

inline

Definition at line 45 of file HcalNoiseAlgo.h.

References r45Count_.

{ return r45Count_; }

r45EnergyFraction()

double CommonHcalNoiseRBXData::r45EnergyFraction ( void  ) const

inline

Definition at line 47 of file HcalNoiseAlgo.h.

References r45EnergyFraction_.

{ return r45EnergyFraction_; }

r45Fraction()

double CommonHcalNoiseRBXData::r45Fraction ( void  ) const

inline

Definition at line 46 of file HcalNoiseAlgo.h.

References r45Fraction_.

{ return r45Fraction_; }

ratio()

double CommonHcalNoiseRBXData::ratio ( void  ) const

inline

Definition at line 25 of file HcalNoiseAlgo.h.

References e10ts_, and e2ts_.

{ return e2ts_ / e10ts_; }

RBXEMF()

double CommonHcalNoiseRBXData::RBXEMF ( void  ) const

inline

Definition at line 41 of file HcalNoiseAlgo.h.

References RBXEMF_.

{ return RBXEMF_; }

rbxTowers()

edm::RefVector<CaloTowerCollection> CommonHcalNoiseRBXData::rbxTowers ( void  ) const

inline

Definition at line 44 of file HcalNoiseAlgo.h.

References rbxtowers_.

{ return rbxtowers_; }

validRatio()

bool CommonHcalNoiseRBXData::validRatio ( void  ) const

inline

Definition at line 28 of file HcalNoiseAlgo.h.

References e10ts_.

{ return e10ts_ != 0.0; }

Member Data Documentation

e10ts_

double CommonHcalNoiseRBXData::e10ts_

private

Definition at line 58 of file HcalNoiseAlgo.h.

Referenced by CommonHcalNoiseRBXData(), e10ts(), ratio(), and validRatio().

e2ts_

double CommonHcalNoiseRBXData::e2ts_

private

Definition at line 57 of file HcalNoiseAlgo.h.

Referenced by CommonHcalNoiseRBXData(), e2ts(), and ratio().

energy_

double CommonHcalNoiseRBXData::energy_

private

Definition at line 56 of file HcalNoiseAlgo.h.

Referenced by CommonHcalNoiseRBXData(), and energy().

highEHitTimeSqrd_

double CommonHcalNoiseRBXData::highEHitTimeSqrd_

private

Definition at line 69 of file HcalNoiseAlgo.h.

Referenced by CommonHcalNoiseRBXData(), and highEHitTimeSqrd().

HPDEMF_

double CommonHcalNoiseRBXData::HPDEMF_

private

Definition at line 71 of file HcalNoiseAlgo.h.

Referenced by CommonHcalNoiseRBXData(), and HPDEMF().

lowEHitTimeSqrd_

double CommonHcalNoiseRBXData::lowEHitTimeSqrd_

private

Definition at line 65 of file HcalNoiseAlgo.h.

Referenced by CommonHcalNoiseRBXData(), and lowEHitTimeSqrd().

maxHighEHitTime_

double CommonHcalNoiseRBXData::maxHighEHitTime_

private

Definition at line 68 of file HcalNoiseAlgo.h.

Referenced by CommonHcalNoiseRBXData(), and maxHighEHitTime().

maxLowEHitTime_

double CommonHcalNoiseRBXData::maxLowEHitTime_

private

Definition at line 64 of file HcalNoiseAlgo.h.

Referenced by CommonHcalNoiseRBXData(), and maxLowEHitTime().

minHighEHitTime_

double CommonHcalNoiseRBXData::minHighEHitTime_

private

Definition at line 67 of file HcalNoiseAlgo.h.

Referenced by CommonHcalNoiseRBXData(), and minHighEHitTime().

minLowEHitTime_

double CommonHcalNoiseRBXData::minLowEHitTime_

private

Definition at line 63 of file HcalNoiseAlgo.h.

Referenced by CommonHcalNoiseRBXData(), and minLowEHitTime().

numHighEHits_

int CommonHcalNoiseRBXData::numHighEHits_

private

Definition at line 70 of file HcalNoiseAlgo.h.

Referenced by CommonHcalNoiseRBXData(), and numHighEHits().

numHPDHits_

int CommonHcalNoiseRBXData::numHPDHits_

private

Definition at line 59 of file HcalNoiseAlgo.h.

Referenced by CommonHcalNoiseRBXData(), and numHPDHits().

numHPDNoOtherHits_

int CommonHcalNoiseRBXData::numHPDNoOtherHits_

private

Definition at line 61 of file HcalNoiseAlgo.h.

Referenced by CommonHcalNoiseRBXData(), and numHPDNoOtherHits().

numLowEHits_

int CommonHcalNoiseRBXData::numLowEHits_

private

Definition at line 66 of file HcalNoiseAlgo.h.

Referenced by CommonHcalNoiseRBXData(), and numLowEHits().

numRBXHits_

int CommonHcalNoiseRBXData::numRBXHits_

private

Definition at line 60 of file HcalNoiseAlgo.h.

Referenced by CommonHcalNoiseRBXData(), and numRBXHits().

numZeros_

int CommonHcalNoiseRBXData::numZeros_

private

Definition at line 62 of file HcalNoiseAlgo.h.

Referenced by CommonHcalNoiseRBXData(), and numZeros().

r45Count_

int CommonHcalNoiseRBXData::r45Count_

private

Definition at line 75 of file HcalNoiseAlgo.h.

Referenced by CommonHcalNoiseRBXData(), and r45Count().

r45EnergyFraction_

double CommonHcalNoiseRBXData::r45EnergyFraction_

private

Definition at line 77 of file HcalNoiseAlgo.h.

Referenced by CommonHcalNoiseRBXData(), and r45EnergyFraction().

r45Fraction_

double CommonHcalNoiseRBXData::r45Fraction_

private

Definition at line 76 of file HcalNoiseAlgo.h.

Referenced by CommonHcalNoiseRBXData(), and r45Fraction().

RBXEMF_

double CommonHcalNoiseRBXData::RBXEMF_

private

Definition at line 72 of file HcalNoiseAlgo.h.

Referenced by CommonHcalNoiseRBXData(), and RBXEMF().

rbxtowers_

edm::RefVector<CaloTowerCollection> CommonHcalNoiseRBXData::rbxtowers_

private

Definition at line 74 of file HcalNoiseAlgo.h.

Referenced by CommonHcalNoiseRBXData(), and rbxTowers().

TS4TS5Decision_

bool CommonHcalNoiseRBXData::TS4TS5Decision_

private

Definition at line 73 of file HcalNoiseAlgo.h.

Referenced by CommonHcalNoiseRBXData(), and PassTS4TS5().