#include <EmDQM.h>

Inheritance diagram for EmDQM:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &)

void	beginJob ()

void	beginRun (edm::Run const &, edm::EventSetup const &)

	EmDQM (const edm::ParameterSet &pset)
	Constructor. More...

void	endJob ()

void	endRun (edm::Run const &, edm::EventSetup const &)

	~EmDQM ()
	Destructor. More...

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

Private Member Functions
bool	checkGeneratedParticlesRequirement (const edm::Event &event)

bool	checkRecoParticlesRequirement (const edm::Event &event)

template<class T >
void	fillHistos (edm::Handle< trigger::TriggerEventWithRefs > &, const edm::Event &, unsigned int, std::vector< reco::Particle > &, bool &)

Private Attributes
DQMStore *	dbe

std::string	dirname_

MonitorElement *	etagen

std::vector< MonitorElement * >	etahist

std::vector< MonitorElement * >	etahistiso

std::vector< MonitorElement * >	etahistisomatch

std::vector< MonitorElement * >	etahistmatch

MonitorElement *	etgen

std::vector< MonitorElement * >	ethist

std::vector< MonitorElement * >	ethistiso

std::vector< MonitorElement * >	ethistisomatch

std::vector< MonitorElement * >	ethistmatch

unsigned int	gencut_

edm::InputTag	gencutCollection_

double	genEtaAcc

double	genEtAcc

std::vector< MonitorElement * >	histEtaIsoOfHltObjMatchToGen

std::vector< MonitorElement * >	histEtaOfHltObjMatchToGen

std::vector< MonitorElement * >	histEtIsoOfHltObjMatchToGen

std::vector< MonitorElement * >	histEtOfHltObjMatchToGen

std::vector< MonitorElement * >	histPhiIsoOfHltObjMatchToGen

std::vector< MonitorElement * >	histPhiOfHltObjMatchToGen

std::set< std::string >	hltCollectionLabelsFound

std::set< std::string >	hltCollectionLabelsMissed

HLTConfigProvider	hltConf_
	The instance of the HLTConfigProvider as a data member. More...

std::vector< std::vector < edm::InputTag > >	isoNames

bool	mcMatchedOnly

std::vector< unsigned int >	nCandCuts

bool	noIsolationPlots

bool	noPhiPlots

unsigned int	numOfHLTCollectionLabels

unsigned int	pathIndex

int	pdgGen

MonitorElement *	phigen

std::vector< MonitorElement * >	phihist

std::vector< MonitorElement * >	phihistiso

std::vector< MonitorElement * >	phihistisomatch

std::vector< MonitorElement * >	phihistmatch

unsigned int	plotBins

std::vector< std::pair< double, double > >	plotBounds

double	plotEtaMax

double	plotEtMin

std::vector< bool >	plotiso

unsigned int	plotMinEtForEtaEffPlot

double	plotPhiMax

double	plotPtMax

double	plotPtMin

GreaterByPt< reco::Particle >	pTComparator_

GreaterByPt< reco::GenParticle >	pTGenComparator_

unsigned int	reqNum

std::vector< std::string >	theHLTCollectionHumanNames

std::vector< edm::InputTag >	theHLTCollectionLabels

std::string	theHltName

std::vector< int >	theHLTOutputTypes

edm::InputTag	theL1Seed

MonitorElement *	total

MonitorElement *	totalmatch

edm::InputTag	triggerobjwithrefs

bool	useHumanReadableHistTitles

unsigned	verbosity

Static Private Attributes
static const unsigned	OUTPUT_ALL = 3

static const unsigned	OUTPUT_ERRORS = 1

static const unsigned	OUTPUT_SILENT = 0

static const unsigned	OUTPUT_WARNINGS = 2

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 24 of file EmDQM.h.

Constructor & Destructor Documentation

EmDQM::EmDQM ( const edm::ParameterSet & pset )

explicit

Constructor.

Definition at line 45 of file EmDQM.cc.

References filters, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), i, cppFunctionSkipper::operator, and PFRecoTauDiscriminationAgainstMuon2_cfi::verbosity.

 {
 
   dbe = edm::Service < DQMStore > ().operator->();
   dbe->setVerbose(0);
 
   //          Read from configuration file                  //
   dirname_="HLT/HLTEgammaValidation/"+pset.getParameter<std::string>("@module_label");
   dbe->setCurrentFolder(dirname_);
 
   triggerobjwithrefs = pset.getParameter<edm::InputTag>("triggerobject");
   pathIndex = pset.getUntrackedParameter<unsigned int>("pathIndex", 0);
   // parameters for generator study
   reqNum    = pset.getParameter<unsigned int>("reqNum");
   pdgGen    = pset.getParameter<int>("pdgGen");
   genEtaAcc = pset.getParameter<double>("genEtaAcc");
   genEtAcc  = pset.getParameter<double>("genEtAcc");
   // plotting parameters (untracked because they don't affect the physics)
   plotEtMin  = pset.getUntrackedParameter<double>("genEtMin",0.);
   plotPtMin  = pset.getUntrackedParameter<double>("PtMin",0.);
   plotPtMax  = pset.getUntrackedParameter<double>("PtMax",1000.);
   plotEtaMax = pset.getUntrackedParameter<double>("EtaMax", 2.7);
   plotPhiMax = pset.getUntrackedParameter<double>("PhiMax", 3.15);
   plotBins   = pset.getUntrackedParameter<unsigned int>("Nbins",40);
   plotMinEtForEtaEffPlot = pset.getUntrackedParameter<unsigned int>("minEtForEtaEffPlot", 15);
   useHumanReadableHistTitles = pset.getUntrackedParameter<bool>("useHumanReadableHistTitles", false);
   mcMatchedOnly = pset.getUntrackedParameter<bool>("mcMatchedOnly", true);
   noPhiPlots = pset.getUntrackedParameter<bool>("noPhiPlots", true);
   noIsolationPlots = pset.getUntrackedParameter<bool>("noIsolationPlots", true);
   verbosity = pset.getUntrackedParameter<unsigned int>("verbosity",0);
 
   //preselction cuts 
   gencutCollection_= pset.getParameter<edm::InputTag>("cutcollection");
   gencut_          = pset.getParameter<int>("cutnum");
 
   //         Read in the Vector of Parameter Sets.          //
   //           Information for each filter-step             //
   std::vector<edm::ParameterSet> filters = 
        pset.getParameter<std::vector<edm::ParameterSet> >("filters");
 
   int i = 0;
   for(std::vector<edm::ParameterSet>::iterator filterconf = filters.begin() ; filterconf != filters.end() ; filterconf++)
   {
 
     theHLTCollectionLabels.push_back(filterconf->getParameter<edm::InputTag>("HLTCollectionLabels"));
     theHLTOutputTypes.push_back(filterconf->getParameter<int>("theHLTOutputTypes"));
     // Grab the human-readable name, if it is not specified, use the Collection Label
     theHLTCollectionHumanNames.push_back(filterconf->getUntrackedParameter<std::string>("HLTCollectionHumanName",theHLTCollectionLabels[i].label()));
 
     std::vector<double> bounds = filterconf->getParameter<std::vector<double> >("PlotBounds");
     // If the size of plot "bounds" vector != 2, abort
     assert(bounds.size() == 2);
     plotBounds.push_back(std::pair<double,double>(bounds[0],bounds[1]));
     isoNames.push_back(filterconf->getParameter<std::vector<edm::InputTag> >("IsoCollections"));
     // If the size of the isoNames vector is not greater than zero, abort
     assert(isoNames.back().size()>0);
     if (isoNames.back().at(0).label()=="none") {
       plotiso.push_back(false);
     } else {
       if (!noIsolationPlots) plotiso.push_back(true);
       else plotiso.push_back(false);
     }
     nCandCuts.push_back(filterconf->getParameter<int>("ncandcut"));
     i++;
   } // END of loop over parameter sets
 
   // Record number of HLTCollectionLabels
   numOfHLTCollectionLabels = theHLTCollectionLabels.size();
   
 }

EmDQM::~EmDQM ( )

Destructor.

Definition at line 366 of file EmDQM.cc.

366 {

367 }

Member Function Documentation

void EmDQM::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

virtual

Implements edm::EDAnalyzer.

Definition at line 472 of file EmDQM.cc.

References eta(), event(), JetPDSkim_cfg::hltResults, i, edm::EventBase::isRealData(), edm::HandleBase::isValid(), n, phi, python.multivaluedict::sort(), pileupDistInMC::total, trigger::TriggerCluster, trigger::TriggerElectron, trigger::TriggerL1IsoEG, trigger::TriggerL1NoIsoEG, trigger::TriggerPhoton, and PFRecoTauDiscriminationAgainstMuon2_cfi::verbosity.

 {
   //           Check if there's enough gen particles        //
   //             of interest                                //
   edm::Handle< edm::View<reco::Candidate> > cutCounter;
   event.getByLabel(gencutCollection_,cutCounter);
   if (cutCounter->size() < (unsigned int)gencut_) {
     //edm::LogWarning("EmDQM") << "Less than "<< gencut_ <<" gen particles with pdgId=" << pdgGen;
     return;
   }
 
 
   // fill L1 and HLT info
   // get objects possed by each filter
   edm::Handle<trigger::TriggerEventWithRefs> triggerObj;
   event.getByLabel(triggerobjwithrefs,triggerObj);
   if(!triggerObj.isValid()) {
     if (verbosity >= OUTPUT_WARNINGS)
        edm::LogWarning("EmDQM") << "parameter triggerobject (" << triggerobjwithrefs << ") does not corresond to a valid TriggerEventWithRefs product. Please check especially the process name (e.g. when running over reprocessed datasets)";
     return;
   }
 
   // Were enough high energy gen particles found?
   if (event.isRealData())
     {
       // running validation on data.
       // TODO: we should check that the entire
       //       run is on the same type (all data or
       //       all MC). Otherwise one gets
       //       uninterpretable results...
       if (!checkRecoParticlesRequirement(event))
         return;
     }
   else
     {
       // MC
       if (!checkGeneratedParticlesRequirement(event))
         // if no, throw event away
         return;
     }
 
 
   // It was an event worth keeping. Continue.
 
   //  Fill the bin labeled "Total"                          //
   //   This will be the number of events looked at.         //
   if (!mcMatchedOnly) total->Fill(numOfHLTCollectionLabels+0.5);
   totalmatch->Fill(numOfHLTCollectionLabels+0.5);
 
 
   //               Fill generator info                      //
   // the gencut_ highest Et generator objects of the preselected type are our matches
 
   std::vector<reco::Particle> sortedGen;
   for(edm::View<reco::Candidate>::const_iterator genpart = cutCounter->begin(); genpart != cutCounter->end();genpart++){
     reco::Particle tmpcand(  genpart->charge(), genpart->p4(), genpart->vertex(),genpart->pdgId(),genpart->status() );
     if (tmpcand.et() >= plotEtMin) {
       sortedGen.push_back(tmpcand);
     }
   }
   std::sort(sortedGen.begin(),sortedGen.end(),pTComparator_ );
 
   // Now the collection of gen particles is sorted by pt.
   // So, remove all particles from the collection so that we 
   // only have the top "1 thru gencut_" particles in it
   if (sortedGen.size() < gencut_){
     return;
   }
   sortedGen.erase(sortedGen.begin()+gencut_,sortedGen.end());
 
   for (unsigned int i = 0 ; i < gencut_ ; i++ ) {
     etgen ->Fill( sortedGen[i].et()  ); //validity has been implicitily checked by the cut on gencut_ above
     if (sortedGen[i].et() > plotMinEtForEtaEffPlot) {
       etagen->Fill( sortedGen[i].eta() );
       if (!noPhiPlots) phigen->Fill( sortedGen[i].phi() );
     }
   } // END of loop over Generated particles
   if (gencut_ >= reqNum && !mcMatchedOnly) total->Fill(numOfHLTCollectionLabels+1.5); // this isn't really needed anymore keep for backward comp.
   if (gencut_ >= reqNum) totalmatch->Fill(numOfHLTCollectionLabels+1.5); // this isn't really needed anymore keep for backward comp.
       
   bool accepted = true;  // flags that the event has been accepted by all filters before
   edm::Handle<edm::TriggerResults> hltResults;
   event.getByLabel(edm::InputTag("TriggerResults","", triggerobjwithrefs.process()), hltResults);
   //            Loop over filter modules                    //
   for(unsigned int n=0; n < numOfHLTCollectionLabels ; n++) {
     // check that there are not less sortedGen particles than nCandCut requires for this filter
     if (sortedGen.size() < nCandCuts.at(n)) {
        if (verbosity >= OUTPUT_ERRORS)
           edm::LogError("EmDQM") << "There are less generated particles than the module '" << theHLTCollectionLabels[n].label() << "' requires.";
        continue;
     }
     std::vector<reco::Particle> sortedGenForFilter(sortedGen);
     sortedGenForFilter.erase(sortedGenForFilter.begin() + nCandCuts.at(n), sortedGenForFilter.end());
 
     // Fill only if this filter was run.
     if (pathIndex != 0 && hltConf_.moduleIndex(pathIndex, theHLTCollectionLabels[n].label()) > hltResults->index(pathIndex)) break;
     // These numbers are from the Parameter Set, such as:
     //   theHLTOutputTypes = cms.uint32(100)
     switch(theHLTOutputTypes[n]) 
     {
       case trigger::TriggerL1NoIsoEG: // Non-isolated Level 1
         fillHistos<l1extra::L1EmParticleCollection>(triggerObj,event,n,sortedGenForFilter,accepted);break;
       case trigger::TriggerL1IsoEG: // Isolated Level 1
         fillHistos<l1extra::L1EmParticleCollection>(triggerObj,event,n,sortedGenForFilter,accepted);break;
       case trigger::TriggerPhoton: // Photon 
         fillHistos<reco::RecoEcalCandidateCollection>(triggerObj,event,n,sortedGenForFilter,accepted);break;
       case trigger::TriggerElectron: // Electron 
         fillHistos<reco::ElectronCollection>(triggerObj,event,n,sortedGenForFilter,accepted);break;
       case trigger::TriggerCluster: // TriggerCluster
         fillHistos<reco::RecoEcalCandidateCollection>(triggerObj,event,n,sortedGenForFilter,accepted);break;
       default: 
         throw(cms::Exception("Release Validation Error") << "HLT output type not implemented: theHLTOutputTypes[n]" );
     }
   } // END of loop over filter modules
 }

void EmDQM::beginJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 125 of file EmDQM.cc.

 {
 
 }

void EmDQM::beginRun	(	edm::Run const &	iRun,
		edm::EventSetup const &	iSetup
	)

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 131 of file EmDQM.cc.

References i, diffTwoXMLs::label, NULL, pileupDistInMC::total, and PFRecoTauDiscriminationAgainstMuon2_cfi::verbosity.

 {
    bool changed(true);
    if (hltConf_.init(iRun, iSetup, triggerobjwithrefs.process(), changed)) {
 
       // if init returns TRUE, initialisation has succeeded!
    
       //edm::Service<TFileService> fs;
       dbe->setCurrentFolder(dirname_);
     
       //  Set up Histogram of Effiency vs Step.                 //
       //   theHLTCollectionLabels is a vector of InputTags      //
       //    from the configuration file.                        //
     
       std::string histName="total_eff";
       std::string histTitle = "total events passing";
       if (!mcMatchedOnly) {
          // This plot will have bins equal to 2+(number of
          //        HLTCollectionLabels in the config file)
          total = dbe->book1D(histName.c_str(),histTitle.c_str(),numOfHLTCollectionLabels+2,0,numOfHLTCollectionLabels+2);
          total->setBinLabel(numOfHLTCollectionLabels+1,"Total");
          total->setBinLabel(numOfHLTCollectionLabels+2,"Gen");
          for (unsigned int u=0; u<numOfHLTCollectionLabels; u++){total->setBinLabel(u+1,theHLTCollectionLabels[u].label().c_str());}
       }
     
       histName="total_eff_MC_matched";
       histTitle="total events passing (mc matched)";
       totalmatch = dbe->book1D(histName.c_str(),histTitle.c_str(),numOfHLTCollectionLabels+2,0,numOfHLTCollectionLabels+2);
       totalmatch->setBinLabel(numOfHLTCollectionLabels+1,"Total");
       totalmatch->setBinLabel(numOfHLTCollectionLabels+2,"Gen");
       for (unsigned int u=0; u<numOfHLTCollectionLabels; u++){totalmatch->setBinLabel(u+1,theHLTCollectionLabels[u].label().c_str());}
     
       MonitorElement* tmphisto;
       MonitorElement* tmpiso;
     
       // Set up generator-level histograms                      //
       std::string pdgIdString;
       switch(pdgGen) {
       case 11:
         pdgIdString="Electron";break;
       case 22:
         pdgIdString="Photon";break;
       default:
         pdgIdString="Particle";
       }
     
       histName = "gen_et";
       histTitle= "E_{T} of " + pdgIdString + "s" ;
       etgen =  dbe->book1D(histName.c_str(),histTitle.c_str(),plotBins,plotPtMin,plotPtMax);
       histName = "gen_eta";
       histTitle= "#eta of "+ pdgIdString +"s " ;
       etagen = dbe->book1D(histName.c_str(),histTitle.c_str(),plotBins,-plotEtaMax,plotEtaMax);
       histName = "gen_phi";
       histTitle= "#phi of "+ pdgIdString +"s " ;
       if (!noPhiPlots) phigen = dbe->book1D(histName.c_str(),histTitle.c_str(),plotBins,-plotPhiMax,plotPhiMax);
     
       
     
       //  Set up histograms of HLT objects                      //
     
       // Determine what strings to use for histogram titles
       std::vector<std::string> HltHistTitle;
       if ( theHLTCollectionHumanNames.size() == numOfHLTCollectionLabels && useHumanReadableHistTitles ) {
         HltHistTitle = theHLTCollectionHumanNames;
       } else {
         for (unsigned int i =0; i < numOfHLTCollectionLabels; i++) {
           HltHistTitle.push_back(theHLTCollectionLabels[i].label());
         }
       }
      
       for(unsigned int i = 0; i< numOfHLTCollectionLabels ; i++){
         if (!mcMatchedOnly) {
            // Et distribution of HLT objects passing filter i
            histName = theHLTCollectionLabels[i].label()+"et_all";
            histTitle = HltHistTitle[i]+" Et (ALL)";
            tmphisto =  dbe->book1D(histName.c_str(),histTitle.c_str(),plotBins,plotPtMin,plotPtMax);
            ethist.push_back(tmphisto);
            
            // Eta distribution of HLT objects passing filter i
            histName = theHLTCollectionLabels[i].label()+"eta_all";
            histTitle = HltHistTitle[i]+" #eta (ALL)";
            tmphisto =  dbe->book1D(histName.c_str(),histTitle.c_str(),plotBins,-plotEtaMax,plotEtaMax);
            etahist.push_back(tmphisto);          
 
            if (!noPhiPlots) {
              // Phi distribution of HLT objects passing filter i
              histName = theHLTCollectionLabels[i].label()+"phi_all";
              histTitle = HltHistTitle[i]+" #phi (ALL)";
              tmphisto =  dbe->book1D(histName.c_str(),histTitle.c_str(),plotBins,-plotPhiMax,plotPhiMax);
              phihist.push_back(tmphisto);
            }
     
      
            // Et distribution of HLT object that is closest delta-R match to sorted gen particle(s)
            histName  = theHLTCollectionLabels[i].label()+"et";
            histTitle = HltHistTitle[i]+" Et";
            tmphisto  = dbe->book1D(histName.c_str(),histTitle.c_str(),plotBins,plotPtMin,plotPtMax);
            histEtOfHltObjMatchToGen.push_back(tmphisto);
     
            // eta distribution of HLT object that is closest delta-R match to sorted gen particle(s)
            histName  = theHLTCollectionLabels[i].label()+"eta";
            histTitle = HltHistTitle[i]+" eta";
            tmphisto  = dbe->book1D(histName.c_str(),histTitle.c_str(),plotBins,-plotEtaMax,plotEtaMax);
            histEtaOfHltObjMatchToGen.push_back(tmphisto);
     
            if (!noPhiPlots) {
              // phi distribution of HLT object that is closest delta-R match to sorted gen particle(s)
              histName  = theHLTCollectionLabels[i].label()+"phi";
              histTitle = HltHistTitle[i]+" phi";
              tmphisto  = dbe->book1D(histName.c_str(),histTitle.c_str(),plotBins,-plotPhiMax,plotPhiMax);
              histPhiOfHltObjMatchToGen.push_back(tmphisto);
            }
        }
     
         // Et distribution of gen object matching HLT object passing filter i
         histName = theHLTCollectionLabels[i].label()+"et_MC_matched";
         histTitle = HltHistTitle[i]+" Et (MC matched)";
         tmphisto =  dbe->book1D(histName.c_str(),histTitle.c_str(),plotBins,plotPtMin,plotPtMax);
         ethistmatch.push_back(tmphisto);
         
         // Eta distribution of gen object matching HLT object passing filter i
         histName = theHLTCollectionLabels[i].label()+"eta_MC_matched";
         histTitle = HltHistTitle[i]+" #eta (MC matched)";
         tmphisto =  dbe->book1D(histName.c_str(),histTitle.c_str(),plotBins,-plotEtaMax,plotEtaMax);
         etahistmatch.push_back(tmphisto);
     
         if (!noPhiPlots) {
           // Phi distribution of gen object matching HLT object passing filter i
           histName = theHLTCollectionLabels[i].label()+"phi_MC_matched";
           histTitle = HltHistTitle[i]+" #phi (MC matched)";
           tmphisto =  dbe->book1D(histName.c_str(),histTitle.c_str(),plotBins,-plotPhiMax,plotPhiMax);
           phihistmatch.push_back(tmphisto);
         }
     
     
         if (!plotiso[i]) {
           tmpiso = NULL;
           if (!mcMatchedOnly) {
              etahistiso.push_back(tmpiso);
              phihistiso.push_back(tmpiso);
              ethistiso.push_back(tmpiso);
              histEtaIsoOfHltObjMatchToGen.push_back(tmpiso);
              histPhiIsoOfHltObjMatchToGen.push_back(tmpiso);
              histEtIsoOfHltObjMatchToGen.push_back(tmpiso);
           }
           etahistisomatch.push_back(tmpiso);
           phihistisomatch.push_back(tmpiso);
           ethistisomatch.push_back(tmpiso);
         } else {
           if (!mcMatchedOnly) {
              // 2D plot: Isolation values vs eta for all objects
              histName  = theHLTCollectionLabels[i].label()+"eta_isolation_all";
              histTitle = HltHistTitle[i]+" isolation vs #eta (all)";
              tmpiso    = dbe->book2D(histName.c_str(),histTitle.c_str(),plotBins,-plotEtaMax,plotEtaMax,plotBins,plotBounds[i].first,plotBounds[i].second);
              etahistiso.push_back(tmpiso);
     
              // 2D plot: Isolation values vs phi for all objects
              histName  = theHLTCollectionLabels[i].label()+"phi_isolation_all";
              histTitle = HltHistTitle[i]+" isolation vs #phi (all)";
              tmpiso    = dbe->book2D(histName.c_str(),histTitle.c_str(),plotBins,-plotPhiMax,plotPhiMax,plotBins,plotBounds[i].first,plotBounds[i].second);
              phihistiso.push_back(tmpiso);
     
              // 2D plot: Isolation values vs et for all objects
              histName  = theHLTCollectionLabels[i].label()+"et_isolation_all";
              histTitle = HltHistTitle[i]+" isolation vs Et (all)";
              tmpiso    = dbe->book2D(histName.c_str(),histTitle.c_str(),plotBins,plotPtMin,plotPtMax,plotBins,plotBounds[i].first,plotBounds[i].second);
              ethistiso.push_back(tmpiso);
      
              // 2D plot: Isolation values vs eta for HLT object that 
              // is closest delta-R match to sorted gen particle(s)
              histName  = theHLTCollectionLabels[i].label()+"eta_isolation";
              histTitle = HltHistTitle[i]+" isolation vs #eta";
              tmpiso    = dbe->book2D(histName.c_str(),histTitle.c_str(),plotBins,-plotEtaMax,plotEtaMax,plotBins,plotBounds[i].first,plotBounds[i].second);
              histEtaIsoOfHltObjMatchToGen.push_back(tmpiso);
     
              // 2D plot: Isolation values vs phi for HLT object that
              // is closest delta-R match to sorted gen particle(s)
              histName  = theHLTCollectionLabels[i].label()+"phi_isolation";
              histTitle = HltHistTitle[i]+" isolation vs #phi";
              tmpiso    = dbe->book2D(histName.c_str(),histTitle.c_str(),plotBins,-plotPhiMax,plotPhiMax,plotBins,plotBounds[i].first,plotBounds[i].second);
              histPhiIsoOfHltObjMatchToGen.push_back(tmpiso);
     
              // 2D plot: Isolation values vs et for HLT object that 
              // is closest delta-R match to sorted gen particle(s)
              histName  = theHLTCollectionLabels[i].label()+"et_isolation";
              histTitle = HltHistTitle[i]+" isolation vs Et";
              tmpiso    = dbe->book2D(histName.c_str(),histTitle.c_str(),plotBins,plotPtMin,plotPtMax,plotBins,plotBounds[i].first,plotBounds[i].second);
              histEtIsoOfHltObjMatchToGen.push_back(tmpiso);
           }
     
           // 2D plot: Isolation values vs eta for matched objects
           histName  = theHLTCollectionLabels[i].label()+"eta_isolation_MC_matched";
           histTitle = HltHistTitle[i]+" isolation vs #eta (mc matched)";
           tmpiso    = dbe->book2D(histName.c_str(),histTitle.c_str(),plotBins,-plotEtaMax,plotEtaMax,plotBins,plotBounds[i].first,plotBounds[i].second);
           etahistisomatch.push_back(tmpiso);
     
           // 2D plot: Isolation values vs phi for matched objects
           histName  = theHLTCollectionLabels[i].label()+"phi_isolation_MC_matched";
           histTitle = HltHistTitle[i]+" isolation vs #phi (mc matched)";
           tmpiso    = dbe->book2D(histName.c_str(),histTitle.c_str(),plotBins,-plotPhiMax,plotPhiMax,plotBins,plotBounds[i].first,plotBounds[i].second);
           phihistisomatch.push_back(tmpiso);
     
     
           // 2D plot: Isolation values vs et for matched objects
           histName  = theHLTCollectionLabels[i].label()+"et_isolation_MC_matched";
           histTitle = HltHistTitle[i]+" isolation vs Et (mc matched)";
           tmpiso    = dbe->book2D(histName.c_str(),histTitle.c_str(),plotBins,plotPtMin,plotPtMax,plotBins,plotBounds[i].first,plotBounds[i].second);
           ethistisomatch.push_back(tmpiso);
     
         } // END of HLT histograms
     
       }
 
       if (changed) {
          // The HLT config has actually changed wrt the previous Run, hence rebook your
          // histograms or do anything else dependent on the revised HLT config
       }
    } else {
       // if init returns FALSE, initialisation has NOT succeeded, which indicates a problem
       // with the file and/or code and needs to be investigated!
       if (verbosity >= OUTPUT_ERRORS)
          edm::LogError("EmDQM") << " HLT config extraction failure with process name '" << triggerobjwithrefs.process() << "'.";
       // In this case, all access methods will return empty values!
    }
 }

bool EmDQM::checkGeneratedParticlesRequirement ( const edm::Event & event )

private

helper to check whether there were enough generator level electrons/photons (MC) or enough reco level electrons/photons to analyze this event.

Returns: if the event has enough of these candidates.

Definition at line 370 of file EmDQM.cc.

References abs, reco::LeafCandidate::et(), eta(), genParticleCandidates2GenParticles_cfi::genParticles, i, edm::HandleBase::isValid(), python.multivaluedict::sort(), and PFRecoTauDiscriminationAgainstMuon2_cfi::verbosity.

 {
    // Decide if this was an event of interest.               //
    //  Did the highest energy particles happen               //
    //  to have |eta| < 2.5 ?  Then continue.                 //
    edm::Handle< edm::View<reco::Candidate> > genParticles;
    event.getByLabel("genParticles", genParticles);
    if(!genParticles.isValid()) {
      if (verbosity >= OUTPUT_WARNINGS)
         edm::LogWarning("EmDQM") << "genParticles invalid.";
      return false;
    }
 
    std::vector<reco::LeafCandidate> allSortedGenParticles;
 
    for(edm::View<reco::Candidate>::const_iterator currentGenParticle = genParticles->begin(); currentGenParticle != genParticles->end(); currentGenParticle++){
 
      // TODO: do we need to check the states here again ?
      // in principle, there should collections produced with the python configuration
      // (other than 'genParticles') which fulfill these criteria
      if (  !( abs((*currentGenParticle).pdgId())==pdgGen  && (*currentGenParticle).status()==1 && (*currentGenParticle).et() > 2.0)  )  continue;
 
      reco::LeafCandidate tmpcand( *(currentGenParticle) );
 
      if (tmpcand.et() < plotEtMin) continue;
 
      allSortedGenParticles.push_back(tmpcand);
    }
 
    std::sort(allSortedGenParticles.begin(), allSortedGenParticles.end(),pTGenComparator_);
 
    // return false if not enough particles found
    if (allSortedGenParticles.size() < gencut_)
      return false;
 
    // additional check (this might be legacy code and we need to check
    // whether this should not be removed ?)
 
    // We now have a sorted collection of all generated particles
    // with pdgId = pdgGen.
    // Loop over them to see if the top gen particles have eta within acceptance
   // bool keepEvent = true;
    for (unsigned int i = 0 ; i < gencut_ ; i++ ) {
      bool inECALgap = fabs(allSortedGenParticles[i].eta()) > 1.4442 && fabs(allSortedGenParticles[i].eta()) < 1.556;
      if ( (fabs(allSortedGenParticles[i].eta()) > genEtaAcc) || inECALgap ) {
        //edm::LogWarning("EmDQM") << "Throwing event away. Gen particle with pdgId="<< allSortedGenParticles[i].pdgId() <<"; et="<< allSortedGenParticles[i].et() <<"; and eta="<< allSortedGenParticles[i].eta() <<" beyond acceptance.";
        return false;
      }
    }
 
    // all tests passed
    return true;
 }

bool EmDQM::checkRecoParticlesRequirement ( const edm::Event & event )

private

similar to checkGeneratedParticlesRequirement(..) but for reconstructed particles. For the moment, there are some additional requirements in the MC version so we can't use the same code for both cases.

Definition at line 427 of file EmDQM.cc.

References edm::HandleBase::isValid(), and PFRecoTauDiscriminationAgainstMuon2_cfi::verbosity.

 {
   // note that this code is very similar to the one in checkGeneratedParticlesRequirement(..)
   // and hopefully can be merged with it at some point in the future
 
   edm::Handle< edm::View<reco::Candidate> > referenceParticles;
   event.getByLabel(gencutCollection_,referenceParticles);
   if(!referenceParticles.isValid()) {
      if (verbosity >= OUTPUT_WARNINGS)
         edm::LogWarning("EmDQM") << "referenceParticles invalid.";
      return false;
   }
 
   std::vector<const reco::Candidate *> allSortedReferenceParticles;
 
   for(edm::View<reco::Candidate>::const_iterator currentReferenceParticle = referenceParticles->begin();
       currentReferenceParticle != referenceParticles->end();
       currentReferenceParticle++)
   {
      if ( currentReferenceParticle->et() <= 2.0)
        continue;
 
      // Note that for determining the overall efficiency,
      // we should only allow
      //
      // HOWEVER: for turn-on curves, we need to let
      //          more electrons pass
      if (currentReferenceParticle->et() < plotEtMin)
        continue;
 
      // TODO: instead of filling a new vector we could simply count here...
      allSortedReferenceParticles.push_back(&(*currentReferenceParticle));
   }
 
    // std::sort(allSortedReferenceParticles.begin(), allSortedReferenceParticles.end(),pTComparator_);
 
    // return false if not enough particles found
    return allSortedReferenceParticles.size() >= gencut_;
 }

void EmDQM::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 852 of file EmDQM.cc.

 {
 
 }

void EmDQM::endRun	(	edm::Run const &	iRun,
		edm::EventSetup const &	iSetup
	)

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 813 of file EmDQM.cc.

References edm::InputTag::encode(), python.multivaluedict::sort(), GlobalPosition_Frontier_DevDB_cff::tag, and PFRecoTauDiscriminationAgainstMuon2_cfi::verbosity.

 {
   // print information about hltCollectionLabels which were not found
   // (but only those which were never found)
 
   // check which ones were never found
   std::vector<std::string> labelsNeverFound;
   
 
   // for (std::set<edm::InputTag>::const_iterator it = hltCollectionLabelsMissed.begin(); it != hltCollectionLabelsMissed.end(); ++it)
   BOOST_FOREACH(const edm::InputTag &tag, hltCollectionLabelsMissed)
   {
     if (hltCollectionLabelsFound.count(tag.encode()) == 0)
       // never found
       labelsNeverFound.push_back(tag.encode());
 
   } // loop over all tags which were missed at least once
 
   if (labelsNeverFound.empty())
     return;
 
   std::sort(labelsNeverFound.begin(), labelsNeverFound.end());
 
   // there was at least one label which was never found
   // (note that this could also be because the corresponding
   // trigger path slowly fades out to zero efficiency)
   if (verbosity >= OUTPUT_WARNINGS)
      edm::LogWarning("EmDQM") << "There were some HLTCollectionLabels which were never found:";
 
   BOOST_FOREACH(const edm::InputTag &tag, labelsNeverFound)
   {
     if (verbosity >= OUTPUT_ALL)
        edm::LogPrint("EmDQM") << "  " << tag;
   }
 }

template<class T >

void EmDQM::fillHistos	(	edm::Handle< trigger::TriggerEventWithRefs > &	triggerObj,
		const edm::Event &	iEvent,
		unsigned int	n,
		std::vector< reco::Particle > &	sortedGen,
		bool &	accepted
	)

private

Definition at line 601 of file EmDQM.cc.

References asciidump::at, deltaR(), alcaDQMUpload::encode(), eta(), edm::Event::getByLabel(), i, edm::HandleBase::isValid(), j, n, phi, pileupDistInMC::total, trigger::TriggerL1IsoEG, trigger::TriggerL1NoIsoEG, and PFRecoTauDiscriminationAgainstMuon2_cfi::verbosity.

 {
   std::vector<edm::Ref<T> > recoecalcands;
   if ( ( triggerObj->filterIndex(theHLTCollectionLabels[n])>=triggerObj->size() )){ // only process if available
     hltCollectionLabelsMissed.insert(theHLTCollectionLabels[n].encode());
     accepted = false;
     return;
   }
 
   hltCollectionLabelsFound.insert(theHLTCollectionLabels[n].encode());
 
   //      Retrieve saved filter objects                     //
   triggerObj->getObjects(triggerObj->filterIndex(theHLTCollectionLabels[n]),theHLTOutputTypes[n],recoecalcands);
   //Danger: special case, L1 non-isolated
   // needs to be merged with L1 iso
   if (theHLTOutputTypes[n] == trigger::TriggerL1NoIsoEG){
     std::vector<edm::Ref<T> > isocands;
     triggerObj->getObjects(triggerObj->filterIndex(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() >= nCandCuts.at(n) && !mcMatchedOnly) 
     total->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 (verbosity >= OUTPUT_ERRORS)
          edm::LogError("EmDQMInvalidRefs") << "Event content inconsistent: TriggerEventWithRefs contains invalid Refs. Invalid refs for: " << theHLTCollectionLabels[n].label() << ". The collection that this module uses may has been dropped in the event.";
       return;
     }
   }
 
   if (!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 < 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 ) {
         histEtOfHltObjMatchToGen[n] ->Fill( recoecalcands[closestEcalCandIndex]->et()  );
         histEtaOfHltObjMatchToGen[n]->Fill( recoecalcands[closestEcalCandIndex]->eta() );
         if (!noPhiPlots) histPhiOfHltObjMatchToGen[n]->Fill( recoecalcands[closestEcalCandIndex]->phi() );
         
         // Also store isolation info
         if (n+1 < numOfHLTCollectionLabels){ // can't plot beyond last
           if (plotiso[n+1] ){  // only plot if requested in config
             for (unsigned int j =  0 ; j < isoNames[n+1].size() ;j++  ){
               edm::Handle<edm::AssociationMap<edm::OneToValue< T , float > > > depMap; 
               iEvent.getByLabel(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! 
                   histEtaIsoOfHltObjMatchToGen[n+1]->Fill( recoecalcands[closestEcalCandIndex]->eta(),mapi->val);
                   histPhiIsoOfHltObjMatchToGen[n+1]->Fill( recoecalcands[closestEcalCandIndex]->phi(),mapi->val);
                   histEtIsoOfHltObjMatchToGen[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
       ethist[n] ->Fill(recoecalcands[i]->et() );
       etahist[n]->Fill(recoecalcands[i]->eta() );
       if (!noPhiPlots) phihist[n]->Fill(recoecalcands[i]->phi() );
 
       //  Plot isolation variables (show the not-yet-cut        //
       //  isolation, i.e. associated to next filter)            //
       if ( n+1 < numOfHLTCollectionLabels ) { // can't plot beyond last
         if (plotiso[n+1]) {
           for (unsigned int j =  0 ; j < isoNames[n+1].size() ;j++  ){
             edm::Handle<edm::AssociationMap<edm::OneToValue< T , float > > > depMap; 
             iEvent.getByLabel(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! 
                 etahistiso[n+1]->Fill(recoecalcands[i]->eta(),mapi->val);
                 phihistiso[n+1]->Fill(recoecalcands[i]->phi(),mapi->val);
                 ethistiso[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 mtachedMcParts = 0;
   float mindist=0.3;
   if(n==0) mindist=0.5; //low L1-resolution => allow wider matching 
   for(unsigned int i =0; i < 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 {
     mtachedMcParts++;
     matchThis = true;
       }
     }
     if ( !matchThis ) {
       accepted = false;
       continue; // only plot matched candidates
     }
     // fill coordinates of mc particle matching trigger object
     ethistmatch[n] ->Fill( sortedGen[i].et()  );
     if (sortedGen[i].et() > plotMinEtForEtaEffPlot) {
       etahistmatch[n]->Fill( sortedGen[i].eta() );
       if (!noPhiPlots) phihistmatch[n]->Fill( sortedGen[i].phi() );
     }
     //  Plot isolation variables (show the not-yet-cut        //
     //  isolation, i.e. associated to next filter)            //
     if (n+1 < numOfHLTCollectionLabels){ // can't plot beyond last
       if (plotiso[n+1] ){  // only plot if requested in config
     for (unsigned int j =  0 ; j < isoNames[n+1].size() ;j++  ){
       edm::Handle<edm::AssociationMap<edm::OneToValue< T , float > > > depMap; 
       iEvent.getByLabel(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
           etahistisomatch[n+1]->Fill(sortedGen[i].eta(),mapi->val);
           phihistisomatch[n+1]->Fill(sortedGen[i].phi(),mapi->val);
           ethistisomatch[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 (mtachedMcParts >= nCandCuts.at(n) && accepted == true)
     totalmatch->Fill(n+0.5);
 }

Member Data Documentation

DQMStore* EmDQM::dbe

private

Definition at line 152 of file EmDQM.h.

std::string EmDQM::dirname_

private

Definition at line 153 of file EmDQM.h.

MonitorElement* EmDQM::etagen

private

Definition at line 148 of file EmDQM.h.

std::vector<MonitorElement*> EmDQM::etahist

private

Definition at line 124 of file EmDQM.h.

std::vector<MonitorElement*> EmDQM::etahistiso

private

Definition at line 134 of file EmDQM.h.

std::vector<MonitorElement*> EmDQM::etahistisomatch

private

Definition at line 137 of file EmDQM.h.

std::vector<MonitorElement*> EmDQM::etahistmatch

private

Definition at line 127 of file EmDQM.h.

MonitorElement* EmDQM::etgen

private

Definition at line 147 of file EmDQM.h.

std::vector<MonitorElement*> EmDQM::ethist

private

Definition at line 126 of file EmDQM.h.

std::vector<MonitorElement*> EmDQM::ethistiso

private

Definition at line 136 of file EmDQM.h.

std::vector<MonitorElement*> EmDQM::ethistisomatch

private

Definition at line 139 of file EmDQM.h.

std::vector<MonitorElement*> EmDQM::ethistmatch

private

Definition at line 129 of file EmDQM.h.

unsigned int EmDQM::gencut_

private

number of generator level particles (electrons/photons) required (for MC)

Definition at line 111 of file EmDQM.h.

edm::InputTag EmDQM::gencutCollection_

private

collection which should be used for generator particles (MC) or reconstructed particles (data).

This collection is used for matching the HLT objects against (e.g. match the HLT object to generated particles or reconstructed electrons/photons).

Definition at line 108 of file EmDQM.h.

double EmDQM::genEtaAcc

private

Definition at line 90 of file EmDQM.h.

double EmDQM::genEtAcc

private

Definition at line 91 of file EmDQM.h.

std::vector<MonitorElement*> EmDQM::histEtaIsoOfHltObjMatchToGen

private

Definition at line 141 of file EmDQM.h.

std::vector<MonitorElement*> EmDQM::histEtaOfHltObjMatchToGen

private

Definition at line 131 of file EmDQM.h.

std::vector<MonitorElement*> EmDQM::histEtIsoOfHltObjMatchToGen

private

Definition at line 140 of file EmDQM.h.

std::vector<MonitorElement*> EmDQM::histEtOfHltObjMatchToGen

private

Definition at line 130 of file EmDQM.h.

std::vector<MonitorElement*> EmDQM::histPhiIsoOfHltObjMatchToGen

private

Definition at line 142 of file EmDQM.h.

std::vector<MonitorElement*> EmDQM::histPhiOfHltObjMatchToGen

private

Definition at line 132 of file EmDQM.h.

std::set<std::string> EmDQM::hltCollectionLabelsFound

private

which hltCollectionLabels were SEEN at least once

Definition at line 114 of file EmDQM.h.

std::set<std::string> EmDQM::hltCollectionLabelsMissed

private

which hltCollectionLabels were MISSED at least once

Definition at line 117 of file EmDQM.h.

HLTConfigProvider EmDQM::hltConf_

private

The instance of the HLTConfigProvider as a data member.

Definition at line 58 of file EmDQM.h.

std::vector<std::vector<edm::InputTag> > EmDQM::isoNames

private

Definition at line 73 of file EmDQM.h.

bool EmDQM::mcMatchedOnly

private

Definition at line 66 of file EmDQM.h.

std::vector<unsigned int> EmDQM::nCandCuts

private

Definition at line 75 of file EmDQM.h.

bool EmDQM::noIsolationPlots

private

Definition at line 68 of file EmDQM.h.

bool EmDQM::noPhiPlots

private

Definition at line 67 of file EmDQM.h.

unsigned int EmDQM::numOfHLTCollectionLabels

private

Definition at line 64 of file EmDQM.h.

const unsigned EmDQM::OUTPUT_ALL = 3

staticprivate

Definition at line 82 of file EmDQM.h.

const unsigned EmDQM::OUTPUT_ERRORS = 1

staticprivate

Definition at line 80 of file EmDQM.h.

const unsigned EmDQM::OUTPUT_SILENT = 0

staticprivate

Definition at line 79 of file EmDQM.h.

const unsigned EmDQM::OUTPUT_WARNINGS = 2

staticprivate

Definition at line 81 of file EmDQM.h.

unsigned int EmDQM::pathIndex

private

Definition at line 62 of file EmDQM.h.

int EmDQM::pdgGen

private

Definition at line 89 of file EmDQM.h.

MonitorElement* EmDQM::phigen

private

Definition at line 149 of file EmDQM.h.

std::vector<MonitorElement*> EmDQM::phihist

private

Definition at line 125 of file EmDQM.h.

std::vector<MonitorElement*> EmDQM::phihistiso

private

Definition at line 135 of file EmDQM.h.

std::vector<MonitorElement*> EmDQM::phihistisomatch

private

Definition at line 138 of file EmDQM.h.

std::vector<MonitorElement*> EmDQM::phihistmatch

private

Definition at line 128 of file EmDQM.h.

unsigned int EmDQM::plotBins

private

Definition at line 98 of file EmDQM.h.

std::vector<std::pair<double,double> > EmDQM::plotBounds

private

Definition at line 74 of file EmDQM.h.

double EmDQM::plotEtaMax

private

Definition at line 94 of file EmDQM.h.

double EmDQM::plotEtMin

private

Definition at line 93 of file EmDQM.h.

std::vector<bool> EmDQM::plotiso

private

Definition at line 72 of file EmDQM.h.

unsigned int EmDQM::plotMinEtForEtaEffPlot

private

Definition at line 99 of file EmDQM.h.

double EmDQM::plotPhiMax

private

Definition at line 95 of file EmDQM.h.

double EmDQM::plotPtMax

private

Definition at line 97 of file EmDQM.h.

double EmDQM::plotPtMin

private

Definition at line 96 of file EmDQM.h.

GreaterByPt<reco::Particle> EmDQM::pTComparator_

private

Definition at line 156 of file EmDQM.h.

GreaterByPt<reco::GenParticle> EmDQM::pTGenComparator_

private

Definition at line 157 of file EmDQM.h.

unsigned int EmDQM::reqNum

private

Definition at line 88 of file EmDQM.h.

std::vector<std::string> EmDQM::theHLTCollectionHumanNames

private

Definition at line 69 of file EmDQM.h.

std::vector<edm::InputTag> EmDQM::theHLTCollectionLabels

private

Definition at line 63 of file EmDQM.h.

std::string EmDQM::theHltName

private

Definition at line 76 of file EmDQM.h.

std::vector<int> EmDQM::theHLTOutputTypes

private

Definition at line 71 of file EmDQM.h.

edm::InputTag EmDQM::theL1Seed

private

Definition at line 70 of file EmDQM.h.

MonitorElement* EmDQM::total

private

Definition at line 144 of file EmDQM.h.

MonitorElement* EmDQM::totalmatch

private

Definition at line 145 of file EmDQM.h.

edm::InputTag EmDQM::triggerobjwithrefs

private

Definition at line 61 of file EmDQM.h.

bool EmDQM::useHumanReadableHistTitles

private

Definition at line 65 of file EmDQM.h.

unsigned EmDQM::verbosity

private

Definition at line 77 of file EmDQM.h.

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation