HiFJRhoProducer Class Reference

#include <HiJetBackground/HiFJRhoProducer/plugins/HiFJRhoProducer.cc>

Inheritance diagram for HiFJRhoProducer:

Public Member Functions
	HiFJRhoProducer (const edm::ParameterSet &)
	~HiFJRhoProducer () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
void	beginStream (edm::StreamID) override
double	calcMd (const reco::Jet *jet)
double	calcMedian (std::vector< double > &v)
void	endStream () override
bool	isPackedCandidate (const reco::Candidate *candidate)
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
bool	checkJetCand
double	etaMaxExcl2_
double	etaMaxExcl_
std::vector< double >	etaRanges
edm::EDGetTokenT< edm::View< reco::Jet > >	jetsToken_
unsigned int	nExcl2_
unsigned int	nExcl_
double	ptMinExcl2_
double	ptMinExcl_
edm::InputTag	src_
bool	usingPackedCand

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 26 of file HiFJRhoProducer.h.

Constructor & Destructor Documentation

HiFJRhoProducer()

HiFJRhoProducer::HiFJRhoProducer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 50 of file HiFJRhoProducer.cc.

    : src_(iConfig.getParameter<edm::InputTag>("jetSource")),
      nExcl_(iConfig.getParameter<int>("nExcl")),
      etaMaxExcl_(iConfig.getParameter<double>("etaMaxExcl")),
      ptMinExcl_(iConfig.getParameter<double>("ptMinExcl")),
      nExcl2_(iConfig.getParameter<int>("nExcl2")),
      etaMaxExcl2_(iConfig.getParameter<double>("etaMaxExcl2")),
      ptMinExcl2_(iConfig.getParameter<double>("ptMinExcl2")),
      checkJetCand(true),
      usingPackedCand(false) {
  jetsToken_ = consumes<edm::View<reco::Jet>>(src_);
 
  //register your products
  produces<std::vector<double>>("mapEtaEdges");
  produces<std::vector<double>>("mapToRho");
  produces<std::vector<double>>("mapToRhoM");
  produces<std::vector<double>>("ptJets");
  produces<std::vector<double>>("areaJets");
  produces<std::vector<double>>("etaJets");
  etaRanges = iConfig.getParameter<std::vector<double>>("etaRanges");
}

References etaRanges, edm::ParameterSet::getParameter(), jetsToken_, and src_.

~HiFJRhoProducer()

HiFJRhoProducer::~HiFJRhoProducer ( )

override

Definition at line 72 of file HiFJRhoProducer.cc.

                                  {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)
}

Member Function Documentation

beginStream()

void HiFJRhoProducer::beginStream ( edm::StreamID  )

overrideprivate

Definition at line 179 of file HiFJRhoProducer.cc.

{}

calcMd()

double HiFJRhoProducer::calcMd ( const reco::Jet *  jet )

private

Definition at line 184 of file HiFJRhoProducer.cc.

                                                 {
  //
  //get md as defined in http://arxiv.org/pdf/1211.2811.pdf
  //
 
  //Loop over the jet constituents
  double sum = 0.;
  for (auto daughter : jet->getJetConstituentsQuick()) {
    if (isPackedCandidate(daughter)) {  //packed candidate situation
      auto part = static_cast<const pat::PackedCandidate*>(daughter);
      sum += sqrt(part->mass() * part->mass() + part->pt() * part->pt()) - part->pt();
    } else {
      auto part = static_cast<const reco::PFCandidate*>(daughter);
      sum += sqrt(part->mass() * part->mass() + part->pt() * part->pt()) - part->pt();
    }
  }
 
  return sum;
}

References isPackedCandidate(), metsig::jet, and mathSSE::sqrt().

Referenced by produce().

calcMedian()

double HiFJRhoProducer::calcMedian ( std::vector< double > &  v )

private

Definition at line 231 of file HiFJRhoProducer.cc.

                                                       {
  //post-condition: After returning, the elements in v may be reordered and the resulting order is implementation defined.
  //works for even and odd collections
  if (v.empty()) {
    return 0.0;
  }
  auto n = v.size() / 2;
  std::nth_element(v.begin(), v.begin() + n, v.end());
  auto med = v[n];
  if (!(v.size() & 1)) {  //If the set size is even
    auto max_it = std::max_element(v.begin(), v.begin() + n);
    med = (*max_it + med) / 2.0;
  }
  return med;
}

References dqmiodumpmetadata::n, and findQualityFiles::v.

Referenced by produce().

endStream()

void HiFJRhoProducer::endStream ( )

overrideprivate

Definition at line 182 of file HiFJRhoProducer.cc.

{}

fillDescriptions()

void HiFJRhoProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 217 of file HiFJRhoProducer.cc.

                                                                                 {
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("jetSource", edm::InputTag("kt4PFJets"));
  desc.add<int>("nExcl", 2);
  desc.add<double>("etaMaxExcl", 2.);
  desc.add<double>("ptMinExcl", 20.);
  desc.add<int>("nExcl2", 2);
  desc.add<double>("etaMaxExcl2", 2.);
  desc.add<double>("ptMinExcl2", 20.);
  desc.add<std::vector<double>>("etaRanges", {});
  descriptions.add("hiFJRhoProducer", desc);
}

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and HLT_FULL_cff::InputTag.

isPackedCandidate()

bool HiFJRhoProducer::isPackedCandidate ( const reco::Candidate *  candidate )

private

Definition at line 204 of file HiFJRhoProducer.cc.

                                                                      {
  if (checkJetCand) {
    if (typeid(pat::PackedCandidate) == typeid(*candidate))
      usingPackedCand = true;
    else if (typeid(reco::PFCandidate) == typeid(*candidate))
      usingPackedCand = false;
    else
      throw cms::Exception("WrongJetCollection", "Jet constituents are not particle flow candidates");
    checkJetCand = false;
  }
  return usingPackedCand;
}

References checkJetCand, Exception, and usingPackedCand.

Referenced by calcMd().

produce()

void HiFJRhoProducer::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 78 of file HiFJRhoProducer.cc.

                                                                           {
  // Get the vector of jets
  edm::Handle<edm::View<reco::Jet>> jets;
  iEvent.getByToken(jetsToken_, jets);
 
  int neta = (int)etaRanges.size();
  auto mapEtaRangesOut = std::make_unique<std::vector<double>>(neta, -999.);
 
  for (int ieta = 0; ieta < neta; ieta++) {
    mapEtaRangesOut->at(ieta) = etaRanges[ieta];
  }
  auto mapToRhoOut = std::make_unique<std::vector<double>>(neta - 1, 1e-6);
  auto mapToRhoMOut = std::make_unique<std::vector<double>>(neta - 1, 1e-6);
 
  int njets = jets->size();
 
  auto ptJetsOut = std::make_unique<std::vector<double>>(njets, 1e-6);
  auto areaJetsOut = std::make_unique<std::vector<double>>(njets, 1e-6);
  auto etaJetsOut = std::make_unique<std::vector<double>>(njets, 1e-6);
 
  double rhoVec[999];
  double rhomVec[999];
  double etaVec[999];
 
  // int neta = (int)mapEtaRangesOut->size();
  int nacc = 0;
  unsigned int njetsEx = 0;
  unsigned int njetsEx2 = 0;
  for (auto jet = jets->begin(); jet != jets->end(); ++jet) {
    if (njetsEx < nExcl_ && fabs(jet->eta()) < etaMaxExcl_ && jet->pt() > ptMinExcl_) {
      njetsEx++;
      continue;
    }
    if (njetsEx2 < nExcl2_ && fabs(jet->eta()) < etaMaxExcl2_ && fabs(jet->eta()) > etaMaxExcl_ &&
        jet->pt() > ptMinExcl2_) {
      njetsEx2++;
      continue;
    }
    float pt = jet->pt();
    float area = jet->jetArea();
    float eta = jet->eta();
 
    if (eta < mapEtaRangesOut->at(0) || eta > mapEtaRangesOut->at(neta - 1))
      continue;
    if (area > 0.) {
      rhoVec[nacc] = pt / area;
      rhomVec[nacc] = calcMd(&*jet) / area;
      etaVec[nacc] = eta;
      ptJetsOut->at(nacc) = pt;
      areaJetsOut->at(nacc) = area;
      etaJetsOut->at(nacc) = eta;
      ++nacc;
    }
  }
 
  ptJetsOut->resize(nacc);
  areaJetsOut->resize(nacc);
  etaJetsOut->resize(nacc);
  //calculate rho and rhom in eta ranges
  double radius = 0.2;  //distance kt clusters needs to be from edge
  for (int ieta = 0; ieta < (neta - 1); ieta++) {
    std::vector<double> rhoVecCur;
    std::vector<double> rhomVecCur;
 
    double etaMin = mapEtaRangesOut->at(ieta) + radius;
    double etaMax = mapEtaRangesOut->at(ieta + 1) - radius;
 
    int naccCur = 0;
    double rhoCurSum = 0.;
    double rhomCurSum = 0.;
    for (int i = 0; i < nacc; i++) {
      if (etaVec[i] >= etaMin && etaVec[i] < etaMax) {
        rhoVecCur.push_back(rhoVec[i]);
        rhomVecCur.push_back(rhomVec[i]);
 
        rhoCurSum += rhoVec[i];
        rhomCurSum += rhomVec[i];
        ++naccCur;
      }  //eta selection
    }    //accepted jet loop
 
    if (naccCur > 0) {
      double rhoCur = calcMedian(rhoVecCur);
      double rhomCur = calcMedian(rhomVecCur);
      mapToRhoOut->at(ieta) = rhoCur;
      mapToRhoMOut->at(ieta) = rhomCur;
    }
  }  //eta ranges
 
  iEvent.put(std::move(mapEtaRangesOut), "mapEtaEdges");
  iEvent.put(std::move(mapToRhoOut), "mapToRho");
  iEvent.put(std::move(mapToRhoMOut), "mapToRhoM");
 
  iEvent.put(std::move(ptJetsOut), "ptJets");
  iEvent.put(std::move(areaJetsOut), "areaJets");
  iEvent.put(std::move(etaJetsOut), "etaJets");
 
  return;
}

References custom_jme_cff::area, calcMd(), calcMedian(), MillePedeFileConverter_cfg::e, PVValHelper::eta, ALCARECOTkAlBeamHalo_cff::etaMax, etaMaxExcl2_, etaMaxExcl_, ALCARECOTkAlBeamHalo_cff::etaMin, etaRanges, mps_fire::i, LEDCalibrationChannels::ieta, iEvent, createfilelist::int, metsig::jet, singleTopDQM_cfi::jets, jetsToken_, eostools::move(), neta, nExcl2_, nExcl_, BTaggingMonitoring_cff::njets, DiDispStaMuonMonitor_cfi::pt, ptMinExcl2_, ptMinExcl_, and CosmicsPD_Skims::radius.

Member Data Documentation

checkJetCand

bool HiFJRhoProducer::checkJetCand

private

Definition at line 55 of file HiFJRhoProducer.h.

Referenced by isPackedCandidate().

etaMaxExcl2_

double HiFJRhoProducer::etaMaxExcl2_

private

Definition at line 52 of file HiFJRhoProducer.h.

Referenced by produce().

etaMaxExcl_

double HiFJRhoProducer::etaMaxExcl_

private

Definition at line 49 of file HiFJRhoProducer.h.

Referenced by produce().

etaRanges

std::vector<double> HiFJRhoProducer::etaRanges

private

Definition at line 54 of file HiFJRhoProducer.h.

Referenced by HiFJRhoProducer(), and produce().

jetsToken_

edm::EDGetTokenT<edm::View<reco::Jet> > HiFJRhoProducer::jetsToken_

private

Definition at line 44 of file HiFJRhoProducer.h.

Referenced by HiFJRhoProducer(), and produce().

nExcl2_

unsigned int HiFJRhoProducer::nExcl2_

private

Definition at line 51 of file HiFJRhoProducer.h.

Referenced by produce().

nExcl_

unsigned int HiFJRhoProducer::nExcl_

private

Definition at line 48 of file HiFJRhoProducer.h.

Referenced by produce().

ptMinExcl2_

double HiFJRhoProducer::ptMinExcl2_

private

Definition at line 53 of file HiFJRhoProducer.h.

Referenced by produce().

ptMinExcl_

double HiFJRhoProducer::ptMinExcl_

private

Definition at line 50 of file HiFJRhoProducer.h.

Referenced by produce().

src_

edm::InputTag HiFJRhoProducer::src_

private

Definition at line 47 of file HiFJRhoProducer.h.

Referenced by HiFJRhoProducer().

usingPackedCand

bool HiFJRhoProducer::usingPackedCand

private

Definition at line 55 of file HiFJRhoProducer.h.

Referenced by isPackedCandidate().