#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 29 of file EmDQMReco.h.

Constructor & Destructor Documentation

template<class T>

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

inline

Definition at line 31 of file EmDQMReco.h.

31 : dqm(d){};

HistoFillerReco::dqm

EmDQMReco * dqm

Definition: EmDQMReco.h:43

template<class T>

HistoFillerReco< T >::~HistoFillerReco ( )

inline

Definition at line 32 of file EmDQMReco.h.

32 {};

Member Function Documentation

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 778 of file EmDQMReco.cc.

References HLT_FULL_cff::DeltaR, HLT_FULL_cff::deltaR, HistoFillerReco< T >::dqm, runTauDisplay::dr, PVValHelper::eta, EmDQMReco::etahistiso, EmDQMReco::etahistisomatchreco, EmDQMReco::ethistiso, EmDQMReco::ethistisomatchreco, dqm::impl::MonitorElement::Fill(), edm::Event::getByToken(), EmDQMReco::histEtaIsoOfHltObjMatchToReco, EmDQMReco::histEtIsoOfHltObjMatchToReco, EmDQMReco::histHltObjMatchToReco, EmDQMReco::histMatchReco, EmDQMReco::histMatchRecoMonPath, EmDQMReco::histPhiIsoOfHltObjMatchToReco, mps_fire::i, HistoFillerReco< T >::isoNameTokens_, edm::HandleBase::isValid(), dqmiolumiharvest::j, dqmiodumpmetadata::n, EmDQMReco::numOfHLTCollectionLabels, phi, EmDQMReco::phiHistIso, EmDQMReco::phiHistIsoMatchReco, EmDQMReco::plotiso, EmDQMReco::recocut_, EmDQMReco::reqNum, EmDQMReco::standardHist, EmDQMReco::theHLTCollectionLabels, EmDQMReco::theHLTOutputTypes, EmDQMReco::totalmatchreco, EmDQMReco::totalreco, trigger::TriggerL1IsoEG, and trigger::TriggerL1NoIsoEG.

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);
   }
 }