HistoFillerReco< T > Class Template Reference

#include <EmDQMReco.h>

Public Member Functions
void	fillHistos (edm::Handle< trigger::TriggerEventWithRefs > &triggerObj, const edm::Event &iEvent, unsigned int n, std::vector< reco::Particle > &sortedReco, bool plotReco, bool plotMonpath)
	HistoFillerReco (EmDQMReco *d)
	~HistoFillerReco ()

Public Attributes
std::vector< edm::EDGetTokenT< edm::AssociationMap< edm::OneToValue< T, float > > > >	isoNameTokens_

Private Attributes
EmDQMReco *	dqm

Detailed Description

template<class T>
class HistoFillerReco< T >

Definition at line 28 of file EmDQMReco.h.

Constructor & Destructor Documentation

HistoFillerReco()

template<class T>

HistoFillerReco< T >::HistoFillerReco ( EmDQMReco *  d )

inline

Definition at line 30 of file EmDQMReco.h.

: dqm(d) {}

~HistoFillerReco()

template<class T>

HistoFillerReco< T >::~HistoFillerReco ( )

inline

Definition at line 31 of file EmDQMReco.h.

{}

Member Function Documentation

fillHistos()

template<class T >

void HistoFillerReco< T >::fillHistos	(	edm::Handle< trigger::TriggerEventWithRefs > &	triggerObj,
		const edm::Event &	iEvent,
		unsigned int	n,
		std::vector< reco::Particle > &	sortedReco,
		bool	plotReco,
		bool	plotMonpath
	)

Definition at line 776 of file EmDQMReco.cc.

Referenced by EmDQMReco::analyze().

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

  //      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
    return;
  }

  if (recoecalcands.size() >= dqm->reqNum)
    dqm->totalreco->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())) {
      edm::LogError("EmDQMReco") << "Event content inconsistent: TriggerEventWithRefs contains "
                                    "invalid Refs"
                                 << std::endl
                                 << "invalid refs for: " << dqm->theHLTCollectionLabels[n].label();
      return;
    }
  }

  //  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++) {
    dqm->standardHist[n]->fill(recoecalcands[i]->p4());

    //  Plot isolation variables (show the not-yet-cut        //
    //  isolation, i.e. associated to next filter)            //
    if (n + 1 < dqm->numOfHLTCollectionLabels) {  // can't plot beyond last
      if (dqm->plotiso[n + 1]) {
        for (unsigned int j = 0; j < isoNameTokens_.size(); j++) {
          edm::Handle<edm::AssociationMap<edm::OneToValue<T, float>>> depMap;
          iEvent.getByToken(isoNameTokens_.at(j), depMap);
          if (depMap.isValid()) {  // Map may not exist if only one candidate
                                   // passes a double filter
            typename edm::AssociationMap<edm::OneToValue<T, float>>::const_iterator mapi =
                depMap->find(recoecalcands[i]);
            if (mapi != depMap->end()) {  // found candidate in isolation map!
              dqm->etahistiso[n + 1]->Fill(recoecalcands[i]->eta(), mapi->val);
              dqm->ethistiso[n + 1]->Fill(recoecalcands[i]->et(), mapi->val);
              dqm->phiHistIso[n + 1]->Fill(recoecalcands[i]->phi(), mapi->val);
            }
          }
        }
      }
    }  // END of if n+1 < then the number of hlt collections
  }

  // Loop over the Reconstructed Particles, and find the        //
  // closest HLT object match.                              //
  if (plotReco == true) {
    for (unsigned int i = 0; i < dqm->recocut_; i++) {
      math::XYZVector currentRecoParticleMomentum = sortedReco[i].momentum();

      // float closestRecoDeltaR = 0.5;
      float closestRecoDeltaR = 1000.;
      int closestRecoEcalCandIndex = -1;
      for (unsigned int j = 0; j < recoecalcands.size(); j++) {
        float deltaR = DeltaR(recoecalcands[j]->momentum(), currentRecoParticleMomentum);

        if (deltaR < closestRecoDeltaR) {
          closestRecoDeltaR = deltaR;
          closestRecoEcalCandIndex = j;
        }
      }

      // If an HLT object was found within some delta-R
      // of this reco particle, store it in a histogram
      if (closestRecoEcalCandIndex >= 0) {
        //        histEtOfHltObjMatchToReco[n] ->Fill(
        //        recoecalcands[closestRecoEcalCandIndex]->et()  );
        //        histEtaOfHltObjMatchToReco[n]->Fill(
        //        recoecalcands[closestRecoEcalCandIndex]->eta() );
        //        histPhiOfHltObjMatchToReco[n]->Fill(
        //        recoecalcands[closestRecoEcalCandIndex]->phi() );

        dqm->histHltObjMatchToReco[n]->fill(recoecalcands[closestRecoEcalCandIndex]->p4());

        // Also store isolation info
        if (n + 1 < dqm->numOfHLTCollectionLabels) {  // can't plot beyond last
          if (dqm->plotiso[n + 1]) {                  // only plot if requested in config
            for (unsigned int j = 0; j < isoNameTokens_.size(); j++) {
              edm::Handle<edm::AssociationMap<edm::OneToValue<T, float>>> depMap;
              iEvent.getByToken(isoNameTokens_.at(j), depMap);
              if (depMap.isValid()) {  // Map may not exist if only one candidate
                                       // passes a double filter
                typename edm::AssociationMap<edm::OneToValue<T, float>>::const_iterator mapi =
                    depMap->find(recoecalcands[closestRecoEcalCandIndex]);
                if (mapi != depMap->end()) {  // found candidate in isolation map!
                  dqm->histEtaIsoOfHltObjMatchToReco[n + 1]->Fill(recoecalcands[closestRecoEcalCandIndex]->eta(),
                                                                  mapi->val);
                  dqm->histEtIsoOfHltObjMatchToReco[n + 1]->Fill(recoecalcands[closestRecoEcalCandIndex]->et(),
                                                                 mapi->val);
                  dqm->histPhiIsoOfHltObjMatchToReco[n + 1]->Fill(recoecalcands[closestRecoEcalCandIndex]->phi(),
                                                                  mapi->val);
                }
              }
            }
          }
        }
      }  // END of if closestEcalCandIndex >= 0
    }

    //        Fill reco matched objects into histograms         //
    unsigned int mtachedRecoParts = 0;
    float minrecodist = 0.3;
    if (n == 0)
      minrecodist = 0.5;  // low L1-resolution => allow wider matching
    for (unsigned int i = 0; i < dqm->recocut_; i++) {
      // match generator candidate
      bool matchThis = false;
      math::XYZVector candDir = sortedReco[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 > minrecodist) {  // it's not really a "match" if it's that far away
          closest = -1;
        } else {
          mtachedRecoParts++;
          matchThis = true;
        }
      }
      if (!matchThis)
        continue;  // only plot matched candidates
      // fill coordinates of mc particle matching trigger object
      //      ethistmatchreco[n] ->Fill( sortedReco[i].et()  );
      //      etahistmatchreco[n]->Fill( sortedReco[i].eta() );
      //      phiHistMatchReco[n]->Fill( sortedReco[i].phi() );
      dqm->histMatchReco[n]->fill(sortedReco[i].p4());

      if (plotMonpath) {
        //        ethistmatchrecomonpath[n]->Fill( sortedReco[i].et() );
        //        etahistmatchrecomonpath[n]->Fill( sortedReco[i].eta() );
        //        phiHistMatchRecoMonPath[n]->Fill( sortedReco[i].phi() );
        dqm->histMatchRecoMonPath[n]->fill(sortedReco[i].p4());
      }
      //  Plot isolation variables (show the not-yet-cut        //
      //  isolation, i.e. associated to next filter)            //
      if (n + 1 < dqm->numOfHLTCollectionLabels) {  // can't plot beyond last
        if (dqm->plotiso[n + 1]) {                  // only plot if requested in config
          for (unsigned int j = 0; j < isoNameTokens_.size(); j++) {
            edm::Handle<edm::AssociationMap<edm::OneToValue<T, float>>> depMapReco;
            iEvent.getByToken(isoNameTokens_.at(j), depMapReco);
            if (depMapReco.isValid()) {  // Map may not exist if only one
                                         // candidate passes a double filter
              typename edm::AssociationMap<edm::OneToValue<T, float>>::const_iterator mapi =
                  depMapReco->find(recoecalcands[closest]);
              if (mapi != depMapReco->end()) {  // found candidate in isolation map!
                dqm->etahistisomatchreco[n + 1]->Fill(sortedReco[i].eta(), mapi->val);
                dqm->ethistisomatchreco[n + 1]->Fill(sortedReco[i].et(), mapi->val);
                dqm->phiHistIsoMatchReco[n + 1]->Fill(sortedReco[i].eta(), mapi->val);
              }
            }
          }
        }
      }  // END of if n+1 < then the number of hlt collections
    }
    // fill total reco matched efficiency
    if (mtachedRecoParts >= dqm->reqNum)
      dqm->totalmatchreco->Fill(n + 0.5);
  }
}

Member Data Documentation

dqm

template<class T>

EmDQMReco* HistoFillerReco< T >::dqm

private

Definition at line 42 of file EmDQMReco.h.

isoNameTokens_

template<class T>

std::vector<edm::EDGetTokenT<edm::AssociationMap<edm::OneToValue<T, float> > > > HistoFillerReco< T >::isoNameTokens_

Definition at line 39 of file EmDQMReco.h.

Referenced by EmDQMReco::EmDQMReco(), and HistoFillerReco< l1extra::L1EmParticleCollection >::fillHistos().