|
|
|
#include <EmDQM.h>
| 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, cmsCodeRules::cppFunctionSkipper::operator, and 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();
}
| 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 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] |
| 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, label, NULL, pileupDistInMC::total, and 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.
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 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 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] |
| 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 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;
}
}
| 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 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);
}
DQMStore* EmDQM::dbe [private] |
std::string EmDQM::dirname_ [private] |
MonitorElement* EmDQM::etagen [private] |
std::vector<MonitorElement*> EmDQM::etahist [private] |
std::vector<MonitorElement*> EmDQM::etahistiso [private] |
std::vector<MonitorElement*> EmDQM::etahistisomatch [private] |
std::vector<MonitorElement*> EmDQM::etahistmatch [private] |
MonitorElement* EmDQM::etgen [private] |
std::vector<MonitorElement*> EmDQM::ethist [private] |
std::vector<MonitorElement*> EmDQM::ethistiso [private] |
std::vector<MonitorElement*> EmDQM::ethistisomatch [private] |
std::vector<MonitorElement*> EmDQM::ethistmatch [private] |
unsigned int EmDQM::gencut_ [private] |
edm::InputTag EmDQM::gencutCollection_ [private] |
double EmDQM::genEtaAcc [private] |
double EmDQM::genEtAcc [private] |
std::vector<MonitorElement*> EmDQM::histEtaIsoOfHltObjMatchToGen [private] |
std::vector<MonitorElement*> EmDQM::histEtaOfHltObjMatchToGen [private] |
std::vector<MonitorElement*> EmDQM::histEtIsoOfHltObjMatchToGen [private] |
std::vector<MonitorElement*> EmDQM::histEtOfHltObjMatchToGen [private] |
std::vector<MonitorElement*> EmDQM::histPhiIsoOfHltObjMatchToGen [private] |
std::vector<MonitorElement*> EmDQM::histPhiOfHltObjMatchToGen [private] |
std::set<std::string> EmDQM::hltCollectionLabelsFound [private] |
std::set<std::string> EmDQM::hltCollectionLabelsMissed [private] |
HLTConfigProvider EmDQM::hltConf_ [private] |
The instance of the HLTConfigProvider as a data member.
std::vector<std::vector<edm::InputTag> > EmDQM::isoNames [private] |
bool EmDQM::mcMatchedOnly [private] |
std::vector<unsigned int> EmDQM::nCandCuts [private] |
bool EmDQM::noIsolationPlots [private] |
bool EmDQM::noPhiPlots [private] |
unsigned int EmDQM::numOfHLTCollectionLabels [private] |
const unsigned EmDQM::OUTPUT_ALL = 3 [static, private] |
const unsigned EmDQM::OUTPUT_ERRORS = 1 [static, private] |
const unsigned EmDQM::OUTPUT_SILENT = 0 [static, private] |
const unsigned EmDQM::OUTPUT_WARNINGS = 2 [static, private] |
unsigned int EmDQM::pathIndex [private] |
int EmDQM::pdgGen [private] |
MonitorElement* EmDQM::phigen [private] |
std::vector<MonitorElement*> EmDQM::phihist [private] |
std::vector<MonitorElement*> EmDQM::phihistiso [private] |
std::vector<MonitorElement*> EmDQM::phihistisomatch [private] |
std::vector<MonitorElement*> EmDQM::phihistmatch [private] |
unsigned int EmDQM::plotBins [private] |
std::vector<std::pair<double,double> > EmDQM::plotBounds [private] |
double EmDQM::plotEtaMax [private] |
double EmDQM::plotEtMin [private] |
std::vector<bool> EmDQM::plotiso [private] |
unsigned int EmDQM::plotMinEtForEtaEffPlot [private] |
double EmDQM::plotPhiMax [private] |
double EmDQM::plotPtMax [private] |
double EmDQM::plotPtMin [private] |
GreaterByPt<reco::Particle> EmDQM::pTComparator_ [private] |
GreaterByPt<reco::GenParticle> EmDQM::pTGenComparator_ [private] |
unsigned int EmDQM::reqNum [private] |
std::vector<std::string> EmDQM::theHLTCollectionHumanNames [private] |
std::vector<edm::InputTag> EmDQM::theHLTCollectionLabels [private] |
std::string EmDQM::theHltName [private] |
std::vector<int> EmDQM::theHLTOutputTypes [private] |
edm::InputTag EmDQM::theL1Seed [private] |
MonitorElement* EmDQM::total [private] |
MonitorElement* EmDQM::totalmatch [private] |
edm::InputTag EmDQM::triggerobjwithrefs [private] |
bool EmDQM::useHumanReadableHistTitles [private] |
unsigned EmDQM::verbosity [private] |
1.7.3