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_);
}
}
}
