void ElectronAnalyser::analyse(const EventPtr event) { histMan_->setCurrentHistogramFolder(histogramFolder_); weight_ = event->weight() * prescale_ * scale_; const ElectronCollection electrons = event->Electrons(); if (singleElectronOnly_) return; histMan_->H1D("Number_Of_Electrons")->Fill(electrons.size(), weight_); for (unsigned int index = 0; index < electrons.size(); ++index) { const ElectronPointer electron(electrons.at(index)); histMan_->H1D("All_Electron_Pt")->Fill(electron->pt(), weight_); histMan_->H1D("All_Electron_Eta")->Fill(electron->eta(), weight_); histMan_->H1D("All_Electron_AbsEta")->Fill(fabs(electron->eta()), weight_); histMan_->H1D("All_Electron_Phi")->Fill(electron->phi(), weight_); histMan_->H1D("All_Electron_pfIsolation_03_deltaBeta")->Fill(electron->PFRelIso03DeltaBeta(), weight_); histMan_->H1D("All_Electron_sigma_ietaieta")->Fill(electron->sigmaIEtaIEta(), weight_); histMan_->H1D("All_Electron_dPhi_in")->Fill(electron->dPhiIn(), weight_); histMan_->H1D("All_Electron_dEta_in")->Fill(electron->dEtaIn(), weight_); histMan_->H1D("All_Electron_HadOverEM")->Fill(electron->HadOverEm(), weight_); histMan_->H1D("All_Electron_mvaTrigV0")->Fill(electron->mvaTrigV0(), weight_); histMan_->H1D("All_Electron_mvaNonTrigV0")->Fill(electron->mvaNonTrigV0(), weight_); histMan_->H1D("All_Electron_dB")->Fill(electron->d0(), weight_); } }
const ElectronCollection TopPairEMuReferenceSelection::goodElectrons(const EventPtr event) const { const ElectronCollection allElectrons(event->Electrons()); ElectronCollection goodIsolatedElectrons; for (unsigned int index = 0; index < allElectrons.size(); ++index) { const ElectronPointer electron(allElectrons.at(index)); if (isGoodElectron(electron) && isIsolatedElectron(electron)) { goodIsolatedElectrons.push_back(electron); } } return goodIsolatedElectrons; }
bool TopPairMuPlusJetsReferenceSelection2011::passesLooseElectronVeto(const EventPtr event) const { unsigned int nLooseElectrons(0); const ElectronCollection allElectrons(event->Electrons()); for (unsigned int index = 0; index < allElectrons.size(); ++index) { const ElectronPointer electron(allElectrons.at(index)); if (isLooseElectron(electron)) ++nLooseElectrons; } //good isolated electron is always a loose electron as well return nLooseElectrons == 0; }
bool QCDPFRelIsoEPlusJetsSelection::hasExactlyOneIsolatedLepton(const EventPtr event) const { const ElectronCollection allElectrons(event->Electrons()); unsigned int nGoodElectrons(0), nGoodIsolatedElectrons(0); for (unsigned int index = 0; index < allElectrons.size(); ++index) { const ElectronPointer electron(allElectrons.at(index)); if (isGoodElectron(electron)) { ++nGoodElectrons; if (isIsolated(electron)) ++nGoodIsolatedElectrons; } } return nGoodElectrons > 0 && nGoodIsolatedElectrons < 2; }
void DiElectronAnalyser::analyse(const EventPtr event) { histMan_->setCurrentHistogramFolder(histogramFolder_); ElectronCollection electrons = event->Electrons(); weight_ = event->weight() * prescale_ * scale_; if (electrons.size() == 2) { ElectronPointer leadingElectron = electrons.front(); ElectronPointer secondElectron = electrons.at(1); histMan_->H1D_BJetBinned("diElectronMass")->Fill(leadingElectron->invariantMass(secondElectron), weight_); } ElectronCollection isolatedElectrons; for (unsigned int index = 0; index < electrons.size(); ++index) { const ElectronPointer electron(electrons.at(index)); if (electron->pfRelativeIsolation(0.3) < 0.1) isolatedElectrons.push_back(electron); } if (isolatedElectrons.size() == 2) { ElectronPointer leadingElectron = isolatedElectrons.front(); ElectronPointer secondElectron = isolatedElectrons.at(1); histMan_->H1D_BJetBinned("diElectronMass_iso")->Fill(leadingElectron->invariantMass(secondElectron), weight_); } }
const LeptonPointer QCDPFRelIsoEPlusJetsSelection::signalLepton(const EventPtr event) const { // if (!hasExactlyOneIsolatedLepton(event)) { // cerr << "An error occurred in QCD*Selection in event = " << event->eventnumber(); // cerr << ", run = " << event->runnumber() << ", lumi = " << event->lumiblock() << "!" << endl; // cerr << "File = " << event->file() << endl; // cerr // << "Access exception: No signal lepton could be found. Event doesn't pass 'hasExactlyOneIsolatedLepton' selection" // << endl; // throw "Access exception: No signal lepton could be found. Event doesn't pass 'hasExactlyOneIsolatedLepton' selection"; // } const ElectronCollection allElectrons(event->Electrons()); ElectronCollection goodElectrons; for (unsigned int index = 0; index < allElectrons.size(); ++index) { const ElectronPointer electron(allElectrons.at(index)); if (isGoodElectron(electron)) { goodElectrons.push_back(electron); } } return MostIsolatedElectron(goodElectrons); }
bool HLTriggerQCDAnalyser::passesTriggerAnalysisSelection(const EventPtr event) const { const ElectronCollection electrons(event->Electrons()); const JetCollection jets(event->Jets()); if (electrons.size() == 0 || jets.size() < 3) return false; unsigned int nElectrons(0); for (unsigned int index = 0; index < electrons.size(); ++index) { const ElectronPointer electron(electrons.at(index)); if (fabs(electron->eta()) < 2.5 && electron->pt() > 20) ++nElectrons; //if more than 2 electrons passing the selection of > 20GeV, reject event } const ElectronPointer mostEnergeticElectron(electrons.front()); //clean jets against electron JetCollection cleanedJets; for (unsigned int index = 0; index < jets.size(); ++index) { const JetPointer jet(jets.at(index)); if (!jet->isWithinDeltaR(0.3, mostEnergeticElectron)) cleanedJets.push_back(jet); } unsigned int nCleanedJetsAbove30GeV(0), nCleanedJetsAbove45GeV(0); for (unsigned int index = 0; index < cleanedJets.size(); ++index) { const JetPointer jet(cleanedJets.at(index)); if (jet->pt() > 45.) ++nCleanedJetsAbove45GeV; if (jet->pt() > 30.) ++nCleanedJetsAbove30GeV; } return nElectrons == 1 && (nCleanedJetsAbove45GeV >= 3 || (nCleanedJetsAbove45GeV >= 2 && nCleanedJetsAbove30GeV >= 3 && event->runnumber() >= 194270)); }
void HLTriggerQCDAnalyser::analyse(const EventPtr event) { unsigned long run(event->runnumber()); //only do this analysis for runs above 193834 as previous runs don't have all triggers //after run 194225 the trigger has been changed to HLT_Ele25_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_TriCentralPFNoPUJet30_30_20 if (!(run >= 193834 && run <= 194225 && event->isRealData())) return; // if (!(run >= 194270 && event->isRealData())) // return; if (event->Electrons().size() == 0) return; const ElectronPointer mostEnergeticElectron(event->Electrons().front()); bool passCaloIdVT_CaloIsoT_TrkIdT_TrkIsoT = event->HLT( HLTriggers::HLT_Ele25_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_TriCentralPFNoPUJet30) || event->HLT(HLTriggers::HLT_Ele25_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_TriCentralPFNoPUJet30_30_20); bool passCaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT = event->HLT( HLTriggers::HLT_Ele25_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_TriCentralPFNoPUJet30); bool passCaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT = event->HLT( HLTriggers::HLT_Ele25_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_TriCentralPFNoPUJet30); bool passCaloIdVT_TrkIdT = event->HLT(HLTriggers::HLT_Ele25_CaloIdVT_TrkIdT_TriCentralPFNoPUJet30) || event->HLT(HLTriggers::HLT_Ele25_CaloIdVT_TrkIdT_TriCentralPFNoPUJet30_30_20); if (passCaloIdVT_CaloIsoT_TrkIdT_TrkIsoT && (!studyExclusiveEvents_ || (!passCaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT && !passCaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT && !passCaloIdVT_TrkIdT))) { int prescale = event->HLTPrescale(HLTriggers::HLT_Ele25_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_TriCentralPFNoPUJet30); if (run >= 194270) prescale = event->HLTPrescale( HLTriggers::HLT_Ele25_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_TriCentralPFNoPUJet30_30_20); eleAnalyser_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_->setPrescale(prescale); metNonIsoRegionAnalyser_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_->setPrescale(prescale); metAntiIDRegionAnalyser_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_->setPrescale(prescale); if (passesTriggerAnalysisSelection(event)) { eleAnalyser_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_->analyse(event); eleAnalyser_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_->analyseElectron(mostEnergeticElectron, event->weight()); } if (passesNonIsoWithoutBtagAndHLT(event)) { QCDNonIsoRegionCount_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_++; metNonIsoRegionAnalyser_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_->analyse(event, qcdNonIsoSelection_->signalLepton(event)); // metNonIsoRegionAnalyser_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_->analyseTransverseMass(event, // qcdNonIsoSelection_->signalLepton(event)); } if (passesAntiIDWithoutBtagAndHLT(event)) { QCDAntiIDRegionCount_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_++; metAntiIDRegionAnalyser_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_->analyse(event, qcdAntiIDSelection_->signalLepton(event)); // metAntiIDRegionAnalyser_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_->analyseTransverseMass(event, // qcdAntiIDSelection_->signalLepton(event)); } if (passesSignalSelectionWithoutBtagAndHLT(event)) { TopSignalRegionCount_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_++; if (topSignalSelection_->hasAtLeastTwoGoodBJets(event)) TopSignal_TwoBtagsRegionCount_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_++; } } if (passCaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT && (!studyExclusiveEvents_ || (!passCaloIdVT_CaloIsoT_TrkIdT_TrkIsoT && !passCaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT && !passCaloIdVT_TrkIdT))) { int prescale = event->HLTPrescale( (HLTriggers::HLT_Ele25_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_TriCentralPFNoPUJet30)); eleAnalyser_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_->setPrescale(prescale); metNonIsoRegionAnalyser_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_->setPrescale(prescale); metAntiIDRegionAnalyser_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_->setPrescale(prescale); if (passesTriggerAnalysisSelection(event)) { eleAnalyser_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_->analyse(event); eleAnalyser_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_->analyseElectron(mostEnergeticElectron, event->weight()); } if (passesNonIsoWithoutBtagAndHLT(event)) { QCDNonIsoRegionCount_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_++; metNonIsoRegionAnalyser_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_->analyse(event, qcdNonIsoSelection_->signalLepton(event)); // metNonIsoRegionAnalyser_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_->analyseTransverseMass(event, // qcdNonIsoSelection_->signalLepton(event)); } if (passesAntiIDWithoutBtagAndHLT(event)) { QCDAntiIDRegionCount_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_++; metAntiIDRegionAnalyser_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_->analyse(event, qcdAntiIDSelection_->signalLepton(event)); // metAntiIDRegionAnalyser_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_->analyseTransverseMass(event, // qcdAntiIDSelection_->signalLepton(event)); } if (passesSignalSelectionWithoutBtagAndHLT(event)) { TopSignalRegionCount_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_++; if (topSignalSelection_->hasAtLeastTwoGoodBJets(event)) TopSignal_TwoBtagsRegionCount_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_++; } } if (passCaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT && (!studyExclusiveEvents_ || (passCaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT && !passCaloIdVT_CaloIsoT_TrkIdT_TrkIsoT && !passCaloIdVT_TrkIdT))) { int prescale = event->HLTPrescale( (HLTriggers::HLT_Ele25_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_TriCentralPFNoPUJet30)); eleAnalyser_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_->setPrescale(prescale); metNonIsoRegionAnalyser_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_->setPrescale(prescale); metAntiIDRegionAnalyser_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_->setPrescale(prescale); if (passesTriggerAnalysisSelection(event)) { eleAnalyser_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_->analyse(event); eleAnalyser_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_->analyseElectron(mostEnergeticElectron, event->weight()); } if (passesNonIsoWithoutBtagAndHLT(event)) { QCDNonIsoRegionCount_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_++; metNonIsoRegionAnalyser_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_->analyse(event, qcdNonIsoSelection_->signalLepton(event)); // metNonIsoRegionAnalyser_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_->analyseTransverseMass(event, // qcdNonIsoSelection_->signalLepton(event)); } if (passesAntiIDWithoutBtagAndHLT(event)) { QCDAntiIDRegionCount_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_++; metAntiIDRegionAnalyser_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_->analyse(event, qcdAntiIDSelection_->signalLepton(event)); // metAntiIDRegionAnalyser_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_->analyseTransverseMass(event, // qcdAntiIDSelection_->signalLepton(event)); } if (passesSignalSelectionWithoutBtagAndHLT(event)) { TopSignalRegionCount_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_++; if (topSignalSelection_->hasAtLeastTwoGoodBJets(event)) TopSignal_TwoBtagsRegionCount_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_++; } } if (passCaloIdVT_TrkIdT && (!studyExclusiveEvents_ || (!passCaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT && !passCaloIdVT_CaloIsoT_TrkIdT_TrkIsoT && !passCaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT))) { int prescale = event->HLTPrescale((HLTriggers::HLT_Ele25_CaloIdVT_TrkIdT_TriCentralPFNoPUJet30)); if (run >= 194270) prescale = event->HLTPrescale((HLTriggers::HLT_Ele25_CaloIdVT_TrkIdT_TriCentralPFNoPUJet30_30_20)); eleAnalyser_CaloIdVT_TrkIdT_->setPrescale(prescale); metNonIsoRegionAnalyser_CaloIdVT_TrkIdT_->setPrescale(prescale); metAntiIDRegionAnalyser_CaloIdVT_TrkIdT_->setPrescale(prescale); if (passesTriggerAnalysisSelection(event)) { eleAnalyser_CaloIdVT_TrkIdT_->analyse(event); eleAnalyser_CaloIdVT_TrkIdT_->analyseElectron(mostEnergeticElectron, event->weight()); } if (passesNonIsoWithoutBtagAndHLT(event)) { QCDNonIsoRegionCount_CaloIdVT_TrkIdT_++; metNonIsoRegionAnalyser_CaloIdVT_TrkIdT_->analyse(event, qcdNonIsoSelection_->signalLepton(event)); // metNonIsoRegionAnalyser_CaloIdVT_TrkIdT_->analyseTransverseMass(event, // qcdNonIsoSelection_->signalLepton(event)); } if (passesAntiIDWithoutBtagAndHLT(event)) { QCDAntiIDRegionCount_CaloIdVT_TrkIdT_++; metAntiIDRegionAnalyser_CaloIdVT_TrkIdT_->analyse(event, qcdAntiIDSelection_->signalLepton(event)); // metAntiIDRegionAnalyser_CaloIdVT_TrkIdT_->analyseTransverseMass(event, // qcdAntiIDSelection_->signalLepton(event)); } if (passesSignalSelectionWithoutBtagAndHLT(event)) { TopSignalRegionCount_CaloIdVT_TrkIdT_++; if (topSignalSelection_->hasAtLeastTwoGoodBJets(event)) TopSignal_TwoBtagsRegionCount_CaloIdVT_TrkIdT_++; } } }
void PhotonAnalyser::analyse(const EventPtr event){ histMan_->setCurrentHistogramFolder(histogramFolder_ + "/AllPhotons"); //const ElectronCollection electrons = topEERefSelection_->goodElectrons(event); //const MuonCollection muons = topMuMuRefSelection_->goodMuons(event); weight_ = event->weight() * prescale_ * scale_; const JetCollection jets = event->Jets(); const ElectronCollection electrons = event->Electrons(); const PhotonCollection photons = event->Photons(); const MuonCollection muons = event->Muons(); histMan_->H1D_BJetBinned("Number_Of_Photons")->Fill(photons.size(), weight_); for (unsigned int index = 0; index < photons.size(); ++index) { const PhotonPointer photon(photons.at(index)); histMan_->H1D_BJetBinned("Photon_Pt")->Fill(photon->pt(), weight_); histMan_->H1D_BJetBinned("Photon_Eta")->Fill(photon->eta(), weight_); histMan_->H1D_BJetBinned("Photon_AbsEta")->Fill(fabs(photon->eta()), weight_); histMan_->H1D_BJetBinned("Photon_Phi")->Fill(photon->phi(), weight_); histMan_->H1D_BJetBinned("Photon_ET")->Fill(photon->et(), weight_); if (photon->isInEndCapRegion()){ histMan_->H1D_BJetBinned("Photon_sigma_ietaieta_endcap")->Fill(photon->sigmaIEtaIEta(), weight_); histMan_->H1D_BJetBinned("Photon_PFChargedHadronIso_endcap")->Fill(photon->PFChargedHadronIso(), weight_); histMan_->H1D_BJetBinned("Photon_PFNeutralHadronIso_endcap")->Fill(photon->PFNeutralHadronIso(), weight_); histMan_->H1D_BJetBinned("Photon_PFPhotonIso_endcap")->Fill(photon->PFPhotonIso(), weight_); histMan_->H1D_BJetBinned("Photon_RhoCorrectedPFChargedHadronIso_endcap")->Fill(photon->RhoCorrectedPFChargedHadronIso(event->rho()), weight_); histMan_->H1D_BJetBinned("Photon_RhoCorrectedPFNeutralHadronIso_endcap")->Fill(photon->RhoCorrectedPFNeutralHadronIso(event->rho()), weight_); histMan_->H1D_BJetBinned("Photon_RhoCorrectedPFPhotonIso_endcap")->Fill(photon->RhoCorrectedPFPhotonIso(event->rho()), weight_); histMan_->H1D_BJetBinned("Photon_SCChIso_endcap")->Fill(photon->phoSCChIso(), weight_); histMan_->H1D_BJetBinned("Photon_SCNuIso_endcap")->Fill(photon->phoSCNuIso(), weight_); histMan_->H1D_BJetBinned("Photon_SCPhIso_endcap")->Fill(photon->phoSCPhIso(), weight_); histMan_->H1D_BJetBinned("Photon_RhoCorrectedSCChIso_endcap")->Fill(photon->RhoCorrectedSCChIso(event->rho()), weight_); histMan_->H1D_BJetBinned("Photon_RhoCorrectedSCNuIso_endcap")->Fill(photon->RhoCorrectedSCNuIso(event->rho()), weight_); histMan_->H1D_BJetBinned("Photon_RhoCorrectedSCPhIso_endcap")->Fill(photon->RhoCorrectedSCPhIso(event->rho()), weight_); histMan_->H2D_BJetBinned("RhoCorrectedSCFRChIso_v._Sigma_ietaieta_endcap")->Fill(photon->RhoCorrectedSCChIso(event->rho()), photon->sigmaIEtaIEta(), weight_); } else if (photon->isInBarrelRegion()) { histMan_->H1D_BJetBinned("Photon_sigma_ietaieta_barrel")->Fill(photon->sigmaIEtaIEta(), weight_); histMan_->H1D_BJetBinned("Photon_PFChargedHadronIso_barrel")->Fill(photon->PFChargedHadronIso(), weight_); histMan_->H1D_BJetBinned("Photon_PFNeutralHadronIso_barrel")->Fill(photon->PFNeutralHadronIso(), weight_); histMan_->H1D_BJetBinned("Photon_PFPhotonIso_barrel")->Fill(photon->PFPhotonIso(), weight_); histMan_->H1D_BJetBinned("Photon_RhoCorrectedPFChargedHadronIso_barrel")->Fill(photon->RhoCorrectedPFChargedHadronIso(event->rho()), weight_); histMan_->H1D_BJetBinned("Photon_RhoCorrectedPFNeutralHadronIso_barrel")->Fill(photon->RhoCorrectedPFNeutralHadronIso(event->rho()), weight_); histMan_->H1D_BJetBinned("Photon_RhoCorrectedPFPhotonIso_barrel")->Fill(photon->RhoCorrectedPFPhotonIso(event->rho()), weight_); histMan_->H1D_BJetBinned("Photon_SCChIso_barrel")->Fill(photon->phoSCChIso(), weight_); histMan_->H1D_BJetBinned("Photon_SCNuIso_barrel")->Fill(photon->phoSCNuIso(), weight_); histMan_->H1D_BJetBinned("Photon_SCPhIso_barrel")->Fill(photon->phoSCPhIso(), weight_); histMan_->H1D_BJetBinned("Photon_RhoCorrectedSCChIso_barrel")->Fill(photon->RhoCorrectedSCChIso(event->rho()), weight_); histMan_->H1D_BJetBinned("Photon_RhoCorrectedSCNuIso_barrel")->Fill(photon->RhoCorrectedSCNuIso(event->rho()), weight_); histMan_->H1D_BJetBinned("Photon_RhoCorrectedSCPhIso_barrel")->Fill(photon->RhoCorrectedSCPhIso(event->rho()), weight_); histMan_->H1D_BJetBinned("Photon_RandConeChIso_barrel")->Fill(photon->phoRandConeChIso(), weight_); histMan_->H1D_BJetBinned("Photon_RandConeNuIso_barrel")->Fill(photon->phoRandConeNuIso(), weight_); histMan_->H1D_BJetBinned("Photon_RandConePhIso_barrel")->Fill(photon->phoRandConePhIso(), weight_); histMan_->H2D_BJetBinned("RhoCorrectedSCFRChIso_v._Sigma_ietaieta_barrel")->Fill(photon->RhoCorrectedSCChIso(event->rho()), photon->sigmaIEtaIEta(), weight_); } histMan_->H1D_BJetBinned("Photon_HadOverEM")->Fill(photon->HadOverEm(), weight_); histMan_->H1D_BJetBinned("Photon_EcalIso")->Fill(photon->ecalIsolation(), weight_); histMan_->H1D_BJetBinned("Photon_HcalIso")->Fill(photon->hcalIsolation(), weight_); histMan_->H1D_BJetBinned("Photon_HcalIso2012")->Fill(photon->hcalIsolation2012(), weight_); histMan_->H1D_BJetBinned("Photon_TrckIso")->Fill(photon->trackerIsolation(), weight_); histMan_->H1D_BJetBinned("Photon_SCeta")->Fill(photon->superClusterEta(), weight_); histMan_->H1D_BJetBinned("Photon_SCphi")->Fill(photon->superClusterPhi(), weight_); histMan_->H1D_BJetBinned("Photon_SCenergy")->Fill(photon->superClusterEnergy(), weight_); histMan_->H1D_BJetBinned("Photon_SCSeedEnergy")->Fill(photon->superClusterSeedEnergy(), weight_); histMan_->H1D_BJetBinned("Photon_E3x3")->Fill(photon->Ecal3x3Cluster(), weight_); histMan_->H1D_BJetBinned("Photon_E5x5")->Fill(photon->Ecal5x5Cluster(), weight_); histMan_->H1D_BJetBinned("Photon_TrkVeto")->Fill(photon->TrackVeto(), weight_); histMan_->H1D_BJetBinned("Photon_ConvSEVeto")->Fill(photon->ConversionSafeElectronVeto(), weight_); histMan_->H1D_BJetBinned("Photon_HtowoE")->Fill(photon->SingleTowerHoE(), weight_); histMan_->H1D_BJetBinned("Photon_PFChargedHadronIso")->Fill(photon->PFChargedHadronIso(), weight_); histMan_->H1D_BJetBinned("Photon_PFNeutralHadronIso")->Fill(photon->PFNeutralHadronIso(), weight_); histMan_->H1D_BJetBinned("Photon_PFPhotonIso")->Fill(photon->PFPhotonIso(), weight_); histMan_->H1D_BJetBinned("Photon_SCChIso")->Fill(photon->phoSCChIso(), weight_); histMan_->H1D_BJetBinned("Photon_SCNuIso")->Fill(photon->phoSCNuIso(), weight_); histMan_->H1D_BJetBinned("Photon_SCPhIso")->Fill(photon->phoSCPhIso(), weight_); histMan_->H1D_BJetBinned("Photon_RandConeChIso")->Fill(photon->phoRandConeChIso(), weight_); histMan_->H1D_BJetBinned("Photon_RandConeNuIso")->Fill(photon->phoRandConeNuIso(), weight_); histMan_->H1D_BJetBinned("Photon_RandConePhIso")->Fill(photon->phoRandConePhIso(), weight_); histMan_->H1D_BJetBinned("Photon_RhoCorrectedSCChIso")->Fill(photon->RhoCorrectedSCChIso(event->rho()), weight_); histMan_->H1D_BJetBinned("Photon_RhoCorrectedSCNuIso")->Fill(photon->RhoCorrectedSCNuIso(event->rho()), weight_); histMan_->H1D_BJetBinned("Photon_RhoCorrectedSCPhIso")->Fill(photon->RhoCorrectedSCPhIso(event->rho()), weight_); histMan_->H2D_BJetBinned("RhoCorrectedSCFRChIso_v._Sigma_ietaieta")->Fill(photon->RhoCorrectedSCChIso(event->rho()), photon->sigmaIEtaIEta(), weight_); for (unsigned int index = 0; index < jets.size(); ++index) { const JetPointer jet(jets.at(index)); histMan_->H1D_BJetBinned("Photon_deltaR_jets")->Fill(photon->deltaR(jet), weight_); } for (unsigned int index = 0; index < electrons.size(); ++index) { const ElectronPointer electron(electrons.at(index)); histMan_->H1D_BJetBinned("Photon_deltaR_electrons")->Fill(photon->deltaR(electron), weight_); } for (unsigned int index = 0; index < muons.size(); ++index) { const MuonPointer muon(muons.at(index)); histMan_->H1D_BJetBinned("Photon_deltaR_muons")->Fill(photon->deltaR(muon), weight_); } } }