FixedGridEnergyDensity Class Reference

#include <FixedGridEnergyDensity.h>

Public Types
enum	EtaRegion { Central, Forward, All }

Public Member Functions
	FixedGridEnergyDensity (const reco::PFCandidateCollection *input)
float	fixedGridRho (EtaRegion etaRegion=Central)
float	fixedGridRho (std::vector< float > &etabins, std::vector< float > &phibins)
	~FixedGridEnergyDensity ()

Private Attributes
const reco::PFCandidateCollection *	pfCandidates

Detailed Description

Definition at line 6 of file FixedGridEnergyDensity.h.

Member Enumeration Documentation

EtaRegion

enum FixedGridEnergyDensity::EtaRegion

Enumerator
Central
Forward
All

Definition at line 10 of file FixedGridEnergyDensity.h.

{ Central, Forward, All };

Constructor & Destructor Documentation

FixedGridEnergyDensity()

FixedGridEnergyDensity::FixedGridEnergyDensity ( const reco::PFCandidateCollection *  input )

inline

Definition at line 8 of file FixedGridEnergyDensity.h.

References input, and pfCandidates.

{ pfCandidates = input; }

~FixedGridEnergyDensity()

FixedGridEnergyDensity::~FixedGridEnergyDensity ( )

inline

Definition at line 9 of file FixedGridEnergyDensity.h.

{}

Member Function Documentation

fixedGridRho() [1/2]

float FixedGridEnergyDensity::fixedGridRho

(

EtaRegion

etaRegion = Central

)

Definition at line 11 of file FixedGridEnergyDensity.cc.

References muon::All, defaults_cfi::etabins, Forward, mps_fire::i, Pi, and TwoPi.

                                                              {
  //define the phi bins
  vector<float> phibins;
  phibins.reserve(10);
  for (int i = 0; i < 10; i++)
    phibins.push_back(-TMath::Pi() + (2 * i + 1) * TMath::TwoPi() / 20.);
  //define the eta bins
  vector<float> etabins;
  if (etaRegion == Central) {
    for (int i = 0; i < 8; ++i)
      etabins.push_back(-2.1 + 0.6 * i);
  } else if (etaRegion == Forward) {
    for (int i = 0; i < 10; ++i) {
      if (i < 5)
        etabins.push_back(-5.1 + 0.6 * i);
      else
        etabins.push_back(2.7 + 0.6 * (i - 5));
    }
  } else if (etaRegion == All) {
    for (int i = 0; i < 18; ++i)
      etabins.push_back(-5.1 + 0.6 * i);
  }
  return fixedGridRho(etabins, phibins);
}

fixedGridRho() [2/2]

float FixedGridEnergyDensity::fixedGridRho	(	std::vector< float > &	etabins,
		std::vector< float > &	phibins
	)

Definition at line 36 of file FixedGridEnergyDensity.cc.

References reco::deltaPhi(), defaults_cfi::etabins, hcalRecHitTable_cff::ieta, hcalRecHitTable_cff::iphi, zmumugammaAnalyzer_cfi::pfCandidates, and jetUpdater_cfi::sort.

                                                                                               {
  float etadist = etabins[1] - etabins[0];
  float phidist = phibins[1] - phibins[0];
  float etahalfdist = (etabins[1] - etabins[0]) / 2.;
  float phihalfdist = (phibins[1] - phibins[0]) / 2.;
  vector<float> sumPFNallSMDQ;
  sumPFNallSMDQ.reserve(etabins.size() * phibins.size());
  for (unsigned int ieta = 0; ieta < etabins.size(); ++ieta) {
    for (unsigned int iphi = 0; iphi < phibins.size(); ++iphi) {
      float pfniso_ieta_iphi = 0;
      for (PFCandidateCollection::const_iterator pf_it = pfCandidates->begin(); pf_it != pfCandidates->end(); pf_it++) {
        if (fabs(etabins[ieta] - pf_it->eta()) > etahalfdist)
          continue;
        if (fabs(reco::deltaPhi(phibins[iphi], pf_it->phi())) > phihalfdist)
          continue;
        pfniso_ieta_iphi += pf_it->pt();
      }
      sumPFNallSMDQ.push_back(pfniso_ieta_iphi);
    }
  }
  float evt_smdq = 0;
  sort(sumPFNallSMDQ.begin(), sumPFNallSMDQ.end());
  if (sumPFNallSMDQ.size() % 2)
    evt_smdq = sumPFNallSMDQ[(sumPFNallSMDQ.size() - 1) / 2];
  else
    evt_smdq = (sumPFNallSMDQ[sumPFNallSMDQ.size() / 2] + sumPFNallSMDQ[(sumPFNallSMDQ.size() - 2) / 2]) / 2.;
  return evt_smdq / (etadist * phidist);
}

Member Data Documentation

pfCandidates

const reco::PFCandidateCollection* FixedGridEnergyDensity::pfCandidates

private

Definition at line 15 of file FixedGridEnergyDensity.h.

Referenced by FixedGridEnergyDensity().