HistoFiller< T > Class Template Reference

#include <EmDQM.h>

Public Member Functions
void	fillHistos (edm::Handle< trigger::TriggerEventWithRefs > &, const edm::Event &, unsigned int, unsigned int, std::vector< reco::Particle > &, bool &)
	HistoFiller (EmDQM *d)
	~HistoFiller ()

Private Attributes
EmDQM *	dqm

Detailed Description

template<class T>
class HistoFiller< T >

Definition at line 47 of file EmDQM.h.

Constructor & Destructor Documentation

template<class T>

HistoFiller< T >::HistoFiller ( EmDQM *  d )

inline

Definition at line 49 of file EmDQM.h.

: dqm(d){};

template<class T>

HistoFiller< T >::~HistoFiller ( )

inline

Definition at line 50 of file EmDQM.h.

{};

Member Function Documentation

template<class T >

void HistoFiller< T >::fillHistos	(	edm::Handle< trigger::TriggerEventWithRefs > &	triggerObj,
		const edm::Event &	iEvent,
		unsigned int	vPos,
		unsigned int	n,
		std::vector< reco::Particle > &	sortedGen,
		bool &	accepted
	)

If this HLT object did not have a gen particle match, go to next HLT object

Definition at line 911 of file EmDQM.cc.

References HLT_FULL_cff::DeltaR, HLT_FULL_cff::deltaR, HistoFillerReco< T >::dqm, runTauDisplay::dr, PVValHelper::eta, mps_fire::i, dqmiolumiharvest::j, dqmiodumpmetadata::n, PVValHelper::phi, EmDQMReco::theHLTCollectionLabels, EmDQMReco::theHLTOutputTypes, trigger::TriggerL1IsoEG, and trigger::TriggerL1NoIsoEG.

Referenced by EmDQM::analyze().

                                                {
  std::vector<edm::Ref<T>> recoecalcands;
  if ((triggerObj->filterIndex(dqm->theHLTCollectionLabels[n]) >= triggerObj->size())) {  // only process if available
    dqm->hltCollectionLabelsMissed.insert(dqm->theHLTCollectionLabels[n].encode());
    accepted = false;
    return;
  }

  dqm->hltCollectionLabelsFound.insert(dqm->theHLTCollectionLabels[n].encode());

  //      Retrieve saved filter objects                     //
  triggerObj->getObjects(
      triggerObj->filterIndex(dqm->theHLTCollectionLabels[n]), dqm->theHLTOutputTypes[n], recoecalcands);
  // Danger: special case, L1 non-isolated
  // needs to be merged with L1 iso
  if (dqm->theHLTOutputTypes[n] == trigger::TriggerL1NoIsoEG) {
    std::vector<edm::Ref<T>> isocands;
    triggerObj->getObjects(triggerObj->filterIndex(dqm->theHLTCollectionLabels[n]), trigger::TriggerL1IsoEG, isocands);
    if (!isocands.empty()) {
      for (unsigned int i = 0; i < isocands.size(); i++)
        recoecalcands.push_back(isocands[i]);
    }
  }  // END of if theHLTOutputTypes == 82

  if (recoecalcands.empty()) {  // stop if no object passed the previous filter
    accepted = false;
    return;
  }

  // if (recoecalcands.size() >= reqNum )
  if (recoecalcands.size() >= dqm->nCandCuts.at(n) && !dqm->mcMatchedOnly_)
    dqm->totals.at(vPos)->Fill(n + 0.5);

  // check for validity                            //
  // prevents crash in CMSSW_3_1_0_pre6            //
  for (unsigned int j = 0; j < recoecalcands.size(); j++) {
    if (!(recoecalcands.at(j).isAvailable())) {
      if (dqm->verbosity_ >= dqm->OUTPUT_ERRORS)
        edm::LogError("EmDQMInvalidRefs") << "Event content inconsistent: TriggerEventWithRefs contains "
                                             "invalid Refs. Invalid refs for: "
                                          << dqm->theHLTCollectionLabels[n].label()
                                          << ". The collection that this module uses may has been dropped in "
                                             "the event.";
      return;
    }
  }

  if (!dqm->mcMatchedOnly_) {
    // Loop over the Generated Particles, and find the        //
    // closest HLT object match.                              //
    // for (unsigned int i=0; i < gencut_; i++) {
    for (unsigned int i = 0; i < dqm->nCandCuts.at(n); i++) {
      math::XYZVector currentGenParticleMomentum = sortedGen[i].momentum();

      float closestDeltaR = 0.5;
      int closestEcalCandIndex = -1;
      for (unsigned int j = 0; j < recoecalcands.size(); j++) {
        float deltaR = DeltaR(recoecalcands[j]->momentum(), currentGenParticleMomentum);

        if (deltaR < closestDeltaR) {
          closestDeltaR = deltaR;
          closestEcalCandIndex = j;
        }
      }

      // If an HLT object was found within some delta-R
      // of this gen particle, store it in a histogram
      if (closestEcalCandIndex >= 0) {
        dqm->histEtOfHltObjMatchToGens.at(vPos).at(n)->Fill(recoecalcands[closestEcalCandIndex]->et());
        dqm->histEtaOfHltObjMatchToGens.at(vPos).at(n)->Fill(recoecalcands[closestEcalCandIndex]->eta());
        if (!dqm->noPhiPlots_)
          dqm->histPhiOfHltObjMatchToGens.at(vPos).at(n)->Fill(recoecalcands[closestEcalCandIndex]->phi());
        dqm->histEtaPhiOfHltObjMatchToGens.at(vPos).at(n)->Fill(recoecalcands[closestEcalCandIndex]->eta(),
                                                                recoecalcands[closestEcalCandIndex]->phi());
      }  // END of if closestEcalCandIndex >= 0
    }

    //  Loop over all HLT objects in this filter step, and    //
    //  fill histograms.                                      //
    //  bool foundAllMatches = false;
    //  unsigned int numOfHLTobjectsMatched = 0;
    for (unsigned int i = 0; i < recoecalcands.size(); i++) {
      // float closestGenParticleDr = 99.0;
      // for(unsigned int j =0; j < gencut_; j++) {
      //  math::XYZVector currentGenParticle = sortedGen[j].momentum();

      //  double currentDeltaR =
      //  DeltaR(recoecalcands[i]->momentum(),currentGenParticle); if (
      //  currentDeltaR < closestGenParticleDr ) {
      //    closestGenParticleDr = currentDeltaR;
      //  }
      //}
      // if ( !(fabs(closestGenParticleDr < 0.3)) ) continue;

      // numOfHLTobjectsMatched++;
      // if (numOfHLTobjectsMatched >= gencut_) foundAllMatches=true;

      // Fill HLT object histograms
      dqm->ethists.at(vPos).at(n)->Fill(recoecalcands[i]->et());
      dqm->etahists.at(vPos).at(n)->Fill(recoecalcands[i]->eta());
      if (!dqm->noPhiPlots_)
        dqm->phihists.at(vPos).at(n)->Fill(recoecalcands[i]->phi());
      dqm->etaphihists.at(vPos).at(n)->Fill(recoecalcands[i]->eta(), recoecalcands[i]->phi());
    }
  }

  //        Fill mc matched objects into histograms         //
  unsigned int matchedMcParts = 0;
  float mindist = 0.3;
  if (n == 0)
    mindist = 0.5;  // low L1-resolution => allow wider matching
  for (unsigned int i = 0; i < dqm->nCandCuts.at(n); ++i) {
    // match generator candidate
    bool matchThis = false;
    math::XYZVector candDir = sortedGen[i].momentum();
    // unsigned int closest = 0;
    double closestDr = 1000.;
    for (unsigned int trigOb = 0; trigOb < recoecalcands.size(); ++trigOb) {
      double dr = DeltaR(recoecalcands[trigOb]->momentum(), candDir);
      if (dr < closestDr) {
        closestDr = dr;
        // closest = trigOb;
      }
      if (closestDr > mindist) {  // it's not really a "match" if it's that far away
        // closest = -1;
      } else {
        matchedMcParts++;
        matchThis = true;
      }
    }
    if (!matchThis) {
      accepted = false;
      continue;  // only plot matched candidates
    }
    // fill coordinates of mc particle matching trigger object
    dqm->ethistmatchs.at(vPos).at(n)->Fill(sortedGen[i].et());
    if (sortedGen[i].et() > dqm->minEtForEtaEffPlot_) {
      dqm->etahistmatchs.at(vPos).at(n)->Fill(sortedGen[i].eta());
      if (!dqm->noPhiPlots_)
        dqm->phihistmatchs.at(vPos).at(n)->Fill(sortedGen[i].phi());
      dqm->etaphihistmatchs.at(vPos).at(n)->Fill(sortedGen[i].eta(), sortedGen[i].phi());
    }
  }
  // fill total mc matched efficiency

  if (matchedMcParts >= dqm->nCandCuts.at(n) && accepted == true)
    dqm->totalmatchs.at(vPos)->Fill(n + 0.5);
}

Member Data Documentation

template<class T>

EmDQM* HistoFiller< T >::dqm

private

Definition at line 62 of file EmDQM.h.