|
|
|
#include <EmDQM.h>
| EmDQM::EmDQM | ( | const edm::ParameterSet & | pset | ) | [explicit] |
Constructor.
Definition at line 43 of file EmDQM.cc.
References filters, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), i, and cmsCodeRules::cppFunctionSkipper::operator.
{
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");
// paramters 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 paramters (untracked because they don't affect the physics)
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);
//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 {
plotiso.push_back(true);
}
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 333 of file EmDQM.cc.
References abs, eta(), event(), genParticleCandidates2GenParticles_cfi::genParticles, i, edm::HandleBase::isValid(), n, phi, python::multivaluedict::sort(), pileupDistInMC::total, trigger::TriggerCluster, trigger::TriggerElectron, trigger::TriggerL1IsoEG, trigger::TriggerL1NoIsoEG, and trigger::TriggerPhoton.
{
// 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()) {
edm::LogWarning("EmDQM") << "triggerobjwithrefs=" << triggerobjwithrefs;
return;
}
// 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::GenParticle> > genParticles;
event.getByLabel("genParticles", genParticles);
if(!genParticles.isValid()) {
edm::LogWarning("EmDQM") << "genParticles invalid.";
return;
}
std::vector<reco::GenParticle> allSortedGenParticles;
for(edm::View<reco::GenParticle>::const_iterator currentGenParticle = genParticles->begin(); currentGenParticle != genParticles->end(); currentGenParticle++){
if ( !( abs((*currentGenParticle).pdgId())==pdgGen && (*currentGenParticle).status()==1 && (*currentGenParticle).et() > 2.0) ) continue;
reco::GenParticle tmpcand( *(currentGenParticle) );
allSortedGenParticles.push_back(tmpcand);
}
std::sort(allSortedGenParticles.begin(), allSortedGenParticles.end(),pTGenComparator_);
// Were enough high energy gen particles found?
if (allSortedGenParticles.size() < gencut_) {
// if no, throw event away
return;
}
// 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.";
keepEvent=false;
break;
}
}
if (!keepEvent) {
return;
}
// It was an event worth keeping. Continue.
// Fill the bin labeled "Total" //
// This will be the number of events looked at. //
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() );
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
sortedGen.erase(sortedGen.begin()+gencut_,sortedGen.end());
if (gencut_ != sortedGen.size() ){
return;
}
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() );
phigen->Fill( sortedGen[i].phi() );
}
} // END of loop over Generated particles
if (gencut_ >= reqNum) 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.
// Loop over filter modules //
for(unsigned int n=0; n < numOfHLTCollectionLabels ; n++) {
// 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,sortedGen);break;
case trigger::TriggerL1IsoEG: // Isolated Level 1
fillHistos<l1extra::L1EmParticleCollection>(triggerObj,event,n,sortedGen);break;
case trigger::TriggerPhoton: // Photon
fillHistos<reco::RecoEcalCandidateCollection>(triggerObj,event,n,sortedGen);break;
case trigger::TriggerElectron: // Electron
fillHistos<reco::ElectronCollection>(triggerObj,event,n,sortedGen);break;
case trigger::TriggerCluster: // TriggerCluster
fillHistos<reco::RecoEcalCandidateCollection>(triggerObj,event,n,sortedGen);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 116 of file EmDQM.cc.
References i, label, NULL, and pileupDistInMC::total.
{
//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";
// 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 " ;
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++){
// 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);
// 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 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);
// 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);
// 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);
// 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);
if (!plotiso[i]) {
tmpiso = NULL;
etahistiso.push_back(tmpiso);
phihistiso.push_back(tmpiso);
ethistiso.push_back(tmpiso);
etahistisomatch.push_back(tmpiso);
phihistisomatch.push_back(tmpiso);
ethistisomatch.push_back(tmpiso);
histEtaIsoOfHltObjMatchToGen.push_back(tmpiso);
histPhiIsoOfHltObjMatchToGen.push_back(tmpiso);
histEtIsoOfHltObjMatchToGen.push_back(tmpiso);
} else {
// 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 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);
// 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);
} // END of HLT histograms
}
}
| void EmDQM::endJob | ( | void | ) | [virtual] |
| void EmDQM::fillHistos | ( | edm::Handle< trigger::TriggerEventWithRefs > & | triggerObj, |
| const edm::Event & | iEvent, | ||
| unsigned int | n, | ||
| std::vector< reco::Particle > & | sortedGen | ||
| ) | [private] |
Definition at line 476 of file EmDQM.cc.
References asciidump::at, deltaR(), eta(), edm::Event::getByLabel(), i, edm::HandleBase::isValid(), j, n, phi, pileupDistInMC::total, trigger::TriggerL1IsoEG, and trigger::TriggerL1NoIsoEG.
{
std::vector<edm::Ref<T> > recoecalcands;
if ( ( triggerObj->filterIndex(theHLTCollectionLabels[n])>=triggerObj->size() )){ // only process if available
return;
}
// 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
return;
}
if (recoecalcands.size() >= reqNum )
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())){
edm::LogError("EmDQM") << "Event content inconsistent: TriggerEventWithRefs contains invalid Refs" << std::endl << "invalid refs for: " << theHLTCollectionLabels[n].label();
return;
}
}
// Loop over the Generated Particles, and find the //
// closest HLT object match. //
for (unsigned int i=0; i < gencut_; 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() );
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++) {
/*
// See if this HLT object has a gen-level match
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 this HLT object did not have a gen particle match, go to next HLT object
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() );
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 < gencut_; 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 ) 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() );
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 >= reqNum )
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::vector<std::vector<edm::InputTag> > EmDQM::isoNames [private] |
unsigned int EmDQM::numOfHLTCollectionLabels [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] |
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] |
1.7.3