当前位置: 首页>>代码示例>>C++>>正文


C++ EventPtr::Electrons方法代码示例

本文整理汇总了C++中EventPtr::Electrons方法的典型用法代码示例。如果您正苦于以下问题:C++ EventPtr::Electrons方法的具体用法?C++ EventPtr::Electrons怎么用?C++ EventPtr::Electrons使用的例子?那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在EventPtr的用法示例。


在下文中一共展示了EventPtr::Electrons方法的9个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C++代码示例。

示例1: analyse

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_);
    }
}
开发者ID:kreczko,项目名称:AnalysisSoftware,代码行数:26,代码来源:ElectronAnalyser.cpp

示例2: goodElectrons

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;
}
开发者ID:,项目名称:,代码行数:13,代码来源:

示例3: nLooseElectrons

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;
}
开发者ID:senkin,项目名称:AnalysisSoftware,代码行数:14,代码来源:TopPairMuPlusJetsReferenceSelection2011.cpp

示例4: hasExactlyOneIsolatedLepton

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;
}
开发者ID:jjacob,项目名称:AnalysisSoftware,代码行数:16,代码来源:QCDPFRelIsoEPlusJetsSelection.cpp

示例5: analyse

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_);
    }
}
开发者ID:TopPairPlusGamma,项目名称:AnalysisTools,代码行数:23,代码来源:DiElectronAnalyser.cpp

示例6: signalLepton

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

}
开发者ID:jjacob,项目名称:AnalysisSoftware,代码行数:23,代码来源:QCDPFRelIsoEPlusJetsSelection.cpp

示例7: passesTriggerAnalysisSelection

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));
}
开发者ID:jjacob,项目名称:AnalysisSoftware,代码行数:36,代码来源:HLTriggerQCDAnalyser.cpp

示例8: analyse

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_++;
        }
//.........这里部分代码省略.........
开发者ID:jjacob,项目名称:AnalysisSoftware,代码行数:101,代码来源:HLTriggerQCDAnalyser.cpp

示例9: analyse

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_);
        }
    
//.........这里部分代码省略.........
开发者ID:TopPairPlusGamma,项目名称:AnalysisTools,代码行数:101,代码来源:PhotonAnalyser.cpp


注:本文中的EventPtr::Electrons方法示例由纯净天空整理自Github/MSDocs等开源代码及文档管理平台,相关代码片段筛选自各路编程大神贡献的开源项目,源码版权归原作者所有,传播和使用请参考对应项目的License;未经允许,请勿转载。