#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 48 of file EmDQM.h.

Constructor & Destructor Documentation

◆ HistoFiller()

template<class T>

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

inline

Definition at line 50 of file EmDQM.h.

50 : dqm(d){};

◆ ~HistoFiller()

template<class T>

HistoFiller< T >::~HistoFiller ( )

inline

Definition at line 51 of file EmDQM.h.

51 {};

Member Function Documentation

◆ fillHistos()

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.

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