#include #include "FWCore/Framework/interface/Frameworkfwd.h" #include "FWCore/Framework/interface/EDAnalyzer.h" #include "FWCore/Framework/interface/Event.h" #include "FWCore/Framework/interface/MakerMacros.h" #include "FWCore/ParameterSet/interface/ParameterSet.h" #include "FWCore/MessageLogger/interface/MessageLogger.h" #include "FWCore/ServiceRegistry/interface/Service.h" #include "CommonTools/UtilAlgos/interface/TFileService.h" #include "HighMassAnalysis/Analysis/interface/HiMassTauAnalysis.h" #include "EgammaAnalysis/ElectronTools/interface/EGammaCutBasedEleId.h" #include "DataFormats/Common/interface/ValueMap.h" #include "RecoEgamma/EgammaTools/interface/ConversionTools.h" #include "TFile.h" #include "TTree.h" class ElecTauAnalysis: public HiMassTauAnalysis{ public: //explicit ElecTauAnalysis(); explicit ElecTauAnalysis(const edm::ParameterSet&); ~ElecTauAnalysis(); private: void fillNtuple(); void setupBranches(); void initializeVectors(); void clearVectors(); bool passedLooseJetId(const pat::Jet&); //define tree variables TTree* _ElecTauTree; //event info vector _PDFWeights; double _ISRGluonWeight; double _ISRGammaWeight; double _FSRWeight; unsigned int _nVtx; unsigned int _nGoodVtx; float _MEt; float _elecRhoIsolation; bool _zEvent; float _zMassDiff; float _zPtAsymm; float _genWCandPt; float _genZCandMass; vector *_nBtagsHiEffTrkCnt ; vector *_nBtagsHiPurityTrkCnt ; vector *_nBTagsCombSecVtxLoose; vector *_nBTagsCombSecVtxMedium; vector *_nBTagsCombSecVtxTight; //unsigned int _nJets; //LeadTau Info vector *_hadTauMatched; vector *_hadTauMotherId; vector *_hadTauIsGenElec; vector *_hadTauGenPt; vector *_hadTauGenE; vector *_hadTauGenEta; vector *_hadTauGenPhi; vector *_hadTauPt; vector *_hadTauAltLVPt; vector *_hadTauAltLVEnergy; vector *_hadTauEnergy; vector *_hadTauEta; vector *_hadTauPhi; vector *_hadTauCharge; vector *_hadTauLeadChargedCandPt; vector *_hadTauLeadChargedCandCharge; vector *_hadTauLeadChargedCandRecHitsSize; vector *_hadTauLeadChargedCandDxyVtx; vector *_hadTauLeadChargedCandDxyBS; vector *_hadTauLeadChargedCandDxyError; vector *_hadTauLeadChargedCandDzVtx; vector *_hadTauLeadChargedCandDzBS; vector *_hadTauLeadChargedCandDzError; vector *_hadTauLeadChargedCandVz; vector *_hadTauLeadChargedCandVy; vector *_hadTauLeadChargedCandVx; vector *_hadTauDeltaZVtx; vector *_hadTauNProngs; vector *_hadTauEmFraction; vector *_hadTauHcalTotOverPLead; vector *_hadTauHcalMaxOverPLead; vector *_hadTauHcal3x3OverPLead; vector *_hadTauElectronPreId; vector *_hadTauModifiedEOverP; vector *_hadTauBremsRecoveryEOverPLead; vector *_hadTauLTSignedIp; vector *_hadTauTrkIsoSumPt; vector *_hadTauECALIsoSumEt; vector *_hadTauIdByDecayModeFinding; vector *_hadTauIdByDecayModeFindingNewDMs; vector *_hadTauIdByDecayModeFindingOldDMs; vector *_hadTauIdByLooseCombinedIsolationDBSumPtCorr; vector *_hadTauIdByMediumCombinedIsolationDBSumPtCorr; vector *_hadTauIdByTightCombinedIsolationDBSumPtCorr; vector *_hadTauIdByLooseCombinedIsolationDBSumPtCorr3Hits; vector *_hadTauIdByMediumCombinedIsolationDBSumPtCorr3Hits; vector *_hadTauIdByTightCombinedIsolationDBSumPtCorr3Hits; vector *_hadTauIdByVLooseIsolationMVA3newDMwLT; vector *_hadTauIdByLooseIsolationMVA3newDMwLT; vector *_hadTauIdByMediumIsolationMVA3newDMwLT; vector *_hadTauIdByTightIsolationMVA3newDMwLT; vector *_hadTauIdByVTightIsolationMVA3newDMwLT; vector *_hadTauIdByVVTightIsolationMVA3newDMwLT; vector *_hadTauIdByIsolationMVA3newDMwLTraw; vector *_hadTauIdByVLooseIsolationMVA3newDMwoLT; vector *_hadTauIdByLooseIsolationMVA3newDMwoLT; vector *_hadTauIdByMediumIsolationMVA3newDMwoLT; vector *_hadTauIdByTightIsolationMVA3newDMwoLT; vector *_hadTauIdByVTightIsolationMVA3newDMwoLT; vector *_hadTauIdByVVTightIsolationMVA3newDMwoLT; vector *_hadTauIdByIsolationMVA3newDMwoLTraw; vector *_hadTauIdByVLooseIsolationMVA3oldDMwLT; vector *_hadTauIdByLooseIsolationMVA3oldDMwLT; vector *_hadTauIdByMediumIsolationMVA3oldDMwLT; vector *_hadTauIdByTightIsolationMVA3oldDMwLT; vector *_hadTauIdByVTightIsolationMVA3oldDMwLT; vector *_hadTauIdByVVTightIsolationMVA3oldDMwLT; vector *_hadTauIdByVLooseIsolationMVA3oldDMwoLT; vector *_hadTauIdByLooseIsolationMVA3oldDMwoLT; vector *_hadTauIdByMediumIsolationMVA3oldDMwoLT; vector *_hadTauIdByTightIsolationMVA3oldDMwoLT; vector *_hadTauIdByVTightIsolationMVA3oldDMwoLT; vector *_hadTauIdByVVTightIsolationMVA3oldDMwoLT; vector *_hadTauIdByLooseElectronRejection; vector *_hadTauIdByMediumElectronRejection; vector *_hadTauIdByTightElectronRejection; vector *_hadTauIdByMVA5rawElectronRejection; vector *_hadTauIdByMVA5ElectronRejectionCat; vector *_hadTauIdByMVA5VLooseElectronRejection; vector *_hadTauIdByMVA5LooseElectronRejection; vector *_hadTauIdByMVA5MediumElectronRejection; vector *_hadTauIdByMVA5TightElectronRejection; vector *_hadTauIdByMVA5VTightElectronRejection; vector *_hadTauIdByElectronDeadECAL; vector *_hadTauIdByLooseMuonRejection; vector *_hadTauIdByMediumMuonRejection; vector *_hadTauIdByTightMuonRejection; vector *_hadTauIdByLooseMuonRejection2; vector *_hadTauIdByMediumMuonRejection2; vector *_hadTauIdByTightMuonRejection2; vector *_hadTauIdByLooseMuonRejection3; vector *_hadTauIdByTightMuonRejection3; vector *_hadTauIdByLooseMVAMuonRejection; vector *_hadTauIdByMediumMVAMuonRejection; vector *_hadTauIdByTightMVAMuonRejection; vector *_hadTauMEtMt; //Elec Info vector *_elecMatched; vector *_elecMotherId; vector *_elecGenPt; vector *_elecGenE; vector *_elecGenEta; vector *_elecGenPhi; vector *_elecPt; vector *_elecE; vector *_elecEta; vector *_elecPhi; vector *_elecCharge; vector *_elecDxyVtx; vector *_elecDxyBS; vector *_elecDxyError; vector *_elecDzVtx; vector *_elecDzBS; vector *_elecDzError; vector *_elecVx; vector *_elecVy; vector *_elecVz; vector *_elecDEtaIn; vector *_elecDPhiIn; vector *_elecSigmaIEtaIEta; vector *_elecHadOverE; vector *_elecFbrem; vector *_elecEOverPIn; vector *_elecOneOverEMinusOneOverPIn; vector *_elecConversionRejection; vector *_elecMissingHits; vector *_elecFiducialFlag; vector *_elecPFIsoDR03SumChargedHadronPt; vector *_elecPFIsoDR03SumNeutralHadronPt; vector *_elecPFIsoDR03SumPhotonPt; vector *_elecPFIsoDR03SumPUPt; vector *_elecPFIsoDR04SumChargedHadronPt; vector *_elecPFIsoDR04SumNeutralHadronPt; vector *_elecPFIsoDR04SumPhotonPt; vector *_elecPFIsoDR04SumPUPt; vector *_elecIsoDR03SumPtTrk; vector *_elecIsoDR03SumEtEcal; vector *_elecIsoDR03SumEtHcal; vector *_elecClass; vector *_elecMVAStatus; vector *_elecMVAOut; vector *_elecMEtMt; vector *_elecMetCosDPhi; //elecTau Info vector *_elecHadTauCosDPhi; vector *_elecHadTauDelatR; vector *_elecHadTauPZeta; vector *_elecHadTauPZetaVis; vector *_elecHadTauMass; vector *_elecHadTauMetMass; vector *_elecHadTauCollMass; //diJet Info vector *_nVBFJets; vector *_diJetMass; vector *_diJetDEta; vector *_jet1Pt; vector *_jet1Eta; vector *_jet2Pt; vector *_jet2Eta; vector *_elecDiJetMass; vector *_hadTauDiJetMass; vector *_elecTauDiJetMass; //vector *_jetPt; //vector *_jetEta; //vector *_diJetDEtaTest; //vector *_diJetMassTest; float _maxDiJetMass; float _maxDiJetDEta; vector *_jetEtSum; vector *_nJets; };