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

Constructor & Destructor Documentation

template<class T>

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

inline

Definition at line 52 of file EmDQM.h.

:dqm(d) {};

template<class T>

HistoFiller< T >::~HistoFiller ( )

inline

Definition at line 53 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
	)

Definition at line 872 of file EmDQM.cc.

References deltaR(), HistoFillerReco< T >::dqm, eta(), i, j, n, 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.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
    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() );
        
        // commented out because uses data not included in HTLDEBUG and uses getByLabel
        //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 < dqm->isoNames[n+1].size() ;j++  ){
        //      edm::Handle<edm::AssociationMap<edm::OneToValue< T , float > > > depMap; 
        //      iEvent.getByLabel(dqm->isoNames[n+1].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[closestEcalCandIndex]);
        //        if (mapi!=depMap->end()) {  // found candidate in isolation map! 
        //          dqm->histEtaIsoOfHltObjMatchToGens.at(vPos).at(n+1)->Fill( recoecalcands[closestEcalCandIndex]->eta(),mapi->val);
        //          dqm->histPhiIsoOfHltObjMatchToGens.at(vPos).at(n+1)->Fill( recoecalcands[closestEcalCandIndex]->phi(),mapi->val);
        //          dqm->histEtIsoOfHltObjMatchToGens.at(vPos).at(n+1) ->Fill( recoecalcands[closestEcalCandIndex]->et(), mapi->val);
        //        }
        //      }
        //    }
        //  }
        //}
      } // 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() );

      // commented out because uses data not included in HTLDEBUG and uses getByLabel
      //  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 < dqm->isoNames[n+1].size() ;j++  ){
      //      edm::Handle<edm::AssociationMap<edm::OneToValue< T , float > > > depMap; 
      //      iEvent.getByLabel(dqm->isoNames[n+1].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->etahistisos.at(vPos).at(n+1)->Fill(recoecalcands[i]->eta(),mapi->val);
      //          dqm->phihistisos.at(vPos).at(n+1)->Fill(recoecalcands[i]->phi(),mapi->val);
      //          dqm->ethistisos.at(vPos).at(n+1)->Fill(recoecalcands[i]->et(),mapi->val);
      //        }
      //      }
      //    }
      //  }
      //} // END of if n+1 < then the number of hlt collections
    }
  }

  //        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() );
    }
    // commented out because uses data not included in HTLDEBUG and uses getByLabel
    //  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 < dqm->isoNames[n+1].size() ;j++  ){
    //      edm::Handle<edm::AssociationMap<edm::OneToValue< T , float > > > depMap; 
    //      iEvent.getByLabel(dqm->isoNames[n+1].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[closest]);
    //        if (mapi!=depMap->end()){  // found candidate in isolation map!
    //          // Only make efficiency plot using photons with some min Et
    //          dqm->etahistisomatchs.at(vPos).at(n+1)->Fill(sortedGen[i].eta(),mapi->val);
    //          dqm->phihistisomatchs.at(vPos).at(n+1)->Fill(sortedGen[i].phi(),mapi->val);
    //          dqm->ethistisomatchs.at(vPos).at(n+1)->Fill(sortedGen[i].et(),mapi->val);
    //        }
    //      }
    //    }
    //  }
    //} // END of if n+1 < then the number of hlt collections
  }
  // 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 59 of file EmDQM.h.