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

Constructor & Destructor Documentation

template<class T>

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

inline

Definition at line 33 of file EmDQMReco.h.

33 :dqm(d) {};

HistoFillerReco::dqm

EmDQMReco * dqm

Definition: EmDQMReco.h:40

template<class T>

HistoFillerReco< T >::~HistoFillerReco ( )

inline

Definition at line 34 of file EmDQMReco.h.

34 {};

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

References deltaR(), HistoFillerReco< T >::dqm, eta(), EmDQMReco::etahistiso, EmDQMReco::etahistisomatchreco, EmDQMReco::ethistiso, EmDQMReco::ethistisomatchreco, MonitorElement::Fill(), edm::Event::getByToken(), EmDQMReco::histEtaIsoOfHltObjMatchToReco, EmDQMReco::histEtIsoOfHltObjMatchToReco, EmDQMReco::histHltObjMatchToReco, EmDQMReco::histMatchReco, EmDQMReco::histMatchRecoMonPath, EmDQMReco::histPhiIsoOfHltObjMatchToReco, i, HistoFillerReco< T >::isoNameTokens_, edm::HandleBase::isValid(), j, n, EmDQMReco::numOfHLTCollectionLabels, p4, 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.size()>0)
       {
         for (unsigned int i=0; i < isocands.size(); i++)
           recoecalcands.push_back(isocands[i]);
       }
   } // END of if theHLTOutputTypes == 82
 
 
   if (recoecalcands.size() < 1){ // 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);
   }
 
 }